PVT Analysis Techniques for Reservoir Fluids
PVT Analysis Techniques for Reservoir Fluids
88 PHASE BEHAVIOR
TABLE 6.1-LABORATORY ANALYS~S PERFORMED ON TABLE 6.2-EXAMPLE GENERAL INFORMATION SHEET
RESERVOIR-OIL ANO GAS-CONDE~SATE SYSTEMS FOR GOOD OIL CO. WELL 4 OIL SAMPLE
Laboratory Analysis Formation Characteristics
Bottomhole sample Name Cretaceous
Recombined composition First well completed /19 (m/d/y)
C7 + TBP distillation o o Original reservoir pressure at 8,692 ft, psig 4,100
C7 + simulated distillation o o Original produced GOR, scf/bbl 600
Constant-composition expansion • • Production rate, B/D 300
Multistage surface separation • o
Separator temperature, ºF 75
Differential liberation • N Separator pressure, psig 200
CVD o • Oil gravity at 60ºF, ºAPI
Multicontact gas injection o o
Datum 8,000
•=standard, o = can be performed, and N = not performed.
Original gas cap No
where Rs = GOR Vg/V0 in scf/STB from the single-stage flash; Last reservoir pressure at 8,500 ft, psig 3,954*
Date I /19 (m/d/y)
w'KJM; .
. . ....... (6.3) Reservoir temperature at 8,500 ft, ºF 217*
Y; = -I-(------=-)-(- - ) '
w-./M. + W-c /M-c
El 1 g 7+ g 7+ Well status Shut in 72 hours
}"C1+
Pressure gauge Amerada
Normal production rate, B/D 300
X¡= . . . . . . . . . . (6.4)
GOR, scf/bbl 600
""
~ (w-/M
Oj 1·) + (w-c
O 7+
/M-c
O 7+
)'
}"C1+ Separator pressure, psig 200
Separator temperature, ºF 75
and M-c Base pressure, psia 14.65
o 7+
Well making water, % water cut o
Sampling Conditions
Surface gas usually contains less than 1 mol% C 7 + material con-
sisting mainly of heptanes and octanes; M-c = 105 is usually a Sample depth, ft 8,500
good assumption. Surface oil contains less lha1+¡ mol% of the light Well status Shut in 72 hours
constituents C 1, C 2 , and nonhydrocarbons. Low-temperature dis-
GOR
tillation can be used to improve the accuracy of reported weight
fractions for intermediate components in the surface oil ( C 3 through Separator pressure, psig
C 6); however, gas chromatography is more widely used. Separator temperature, ºF
The most probable source of error in wellstream composition of a
Tubing pressure, psig 1,400
bottornhole sample is the surface-oil molecular weight, M 0 , which
appears in Eq. 6.2 for Fg and Eq. 6.4 for X¡. M 0 is usually accurate Casing pressure, psig
within ± 4 to 10%. In Chap. 5, we showed that the Watson character- Sampled by
ization factor, Kw, of surface oil (Eq. 5.35) should be constant (to
Sampler type Wofford
within ± 0.03 of the deterrnined value) for a given reservoir. Oncean
"Pressure and temperatura extrapolated to the midpoint of the producing
average has been established for a reservoir (usually requiring three
interval = 4,01 O psig and 220ºF.
separate measurements), potential errors in M0 can be checked. A
calculated Kw that deviates from the field-average Kw by more than
± 0.03 may indicate an erroneous molecular-weight measurement. and compositions of the resulting surface oil and gas, and recombin-
Eqs. 6.1 through 6.4 show that ali component compositions are ing these compositions to give the separator-oil composition; Tables
affected by M 0 c , which is backcalculated from M 0 with Eq. 6.5 and 6.6 report the results.
6.5. Fortunately, th1+amount oflighter components (particularly C 1) Separator gas is introduced directly into a gas chromatograph,
in the surface oil are small, so the real effect on conversion from which yields weight fractions, w 8. These weight fractions are con-
weight to mole fractions of the surface oil usually is not significant. verted to mole fractions, y¡, by use of appropriate molecular
weights; Tables 6.5 and 6.6 show the results. C 7 + molecular weight
6.2.2 Recombined Samples. Tables 6.5 and 6.6 present the separa- is backcalculated with measured separator-gas specific gravity, Y¡·
tor-oil and -gas compositional analyses of a gas-condensate fluid
and recombined wellstream composition. The separator-oil com-
position is obtained by use of the same procedure as that used for
bottornhole oil samples (Eqs. 6.1 through 6.5). This involves bring- ....... (6.6)
ing the separator oil to standard conditions, measuring properties
Formation Characteristics
Formation name Pay sand
Date first well completed /19 (m/d/y)
Original reservoir pressure at 10, 148 ft, psig 5,713
Original produced-gas/liquid ratio, scf/bbl
Production rate, B/D
Separator pressure, psig
Separator temperature, ºF
Liquid gravity at 60ºF, ºAPI
Datum, ft subsea 8,000
Well Characteristics
Elevation, ft KB 2,214
Total depth, ft 10,348
Producing interval, ft 10,124to 10,176
Tubing size, in. 2
Tubing depth, ft 10, 100
Open-flow potential, MMscf/D
Last reservoir pressure at 1O,148 ft, psig 5,713
Date / /19 (m/d/y)
Reservoir temperature at 1O,148 ft, ºF 186
Status of well status Shut in
Pressure gauge Amerada
Sampling Conditions
Flowing tubing pressure, psig 3,375
FBHP, psig 5,500
Primary-separator pressure, psig 300
Primary-separator temperature, ºF 62
Secondary-separator pressure, psig 20
Secondary-separator temperature, ºF 60
Field stock-tank-liquid gravity at 60ºF, ºAPI 58.5
Primary-separator-gas production rate, Mscf/D 762.14
Pressure base, psia 14.696
Temperature base, ºF 60
Compressibility factor, Fpv 1.043
Gas gravity (laboratory) 0.737
Gas-gravity factor, Fg 0.902
Stock-tank-liquid production rate at 60ºF, B/D 127.3
Primary-separator-gas/stock-tank-liquid ratio
In scf/bbl 5,987
In bbl/MMscf 167.0
Sampled by
For the example PVT report (Tables 6.5 and 6.6), the separator separator gas/wellstream ratio" (801.66 Mscf/MMscf), which is
gas/oil ratio, R,P, during sampling is reported as standard gas vol- equivalen! to mole per mole (Fgsp =0.80166 mol/mol). The re-
ume per separator-oil volume (4,428 scf/bbl). In this report, the units ported value of Fgsp can be checked with
are incorrectly labeled scf/bbl at 60ºF, where in fact the separator-oil -1
volume is measured at separator pressure (300 psig) and tempera-
ture (62ºF). The separator-oil forrnation volume factor (FVF), Bosp• Fgsp = (] + 2~30~osp)
osp sp
' . , ....... . , ......... , (6.8)
is 1.352 bbl/STB and represents the volume of separator oil required
to yield 1 STB of oil (i.e., condensate). N
The equation used to calculate wellstream composition, Z;, is where Mosp = L X ¡M;. . • . . . . . . . . . . . . . . . . . . . . . . . . . (6.9)
i=l
Z; = FgspY; + (1 - Fgsp)X;, . . . . . . . . . . . . . . . . . . . . . (6.7)
p osp in lbm/ft3 is calculated with a correlation (e.g., Standing-Katz9)
where Fgsp = mole fraction of wellstream mixture that becomes or with the relation (62.4y 0 + 0.0136y 8R,)/B 0 , where R, and
separator gas. In the laboratory report, Fgsp is reported as "primary- B0 = separator-oil values in scf/STB and bbl/STB, respectively;
90 PHASE BEHAVIOR
TABLE 6.4-WELLSTREAM (RESERVOIR-FLUID) removed gas, n 8; and specific gravity of removed gas, y-. If re-
COMPOSITION FOR GOOD OIL CO. WELL 4 quested, the gas samples can be analyzed chromatographi~ally to
BOTT9MHGl;E-0J!:. SAMPLE give molar composition, y.
Lmol% wt%
o~·
(g/cm~ ºAPI*
Molecular
Weight
The oil remaining after gas removal is brought to the conditions
of the next separator stage. The gas is removed again and quantified
by moles and specific gravity. Oil volume is noted, and the process
Nil Nil
is repeated until stock-tank conditions are reached. Final oil volume,
0.91 0.43 V0 , and specific gravity, y 0 , are measured at 60ºF.
N2 0.16 0.05 Table 6.7 gives results from four separator tests, each consisting
Methane 36.47 6.24 of two stages of separation. The first-stage-separator pressure is var-
Ethane 9.67 3.10 ied from 50 to 300 psig, and stock-tank conditions are held constant
at Opsig and 75ºF. GOR's are reported as standard gas volume per
Propane 6.95 3.27
separator-oil volume, Rsp• and as standard gas volume per stock-
í-butane 1.44 0.89 tank-oil volume, R,, respectively.
n-butane 3.93 2.44
!!..V-
i-pentane 1.44 1.11 !!..Rsp = V 8 ............ .... .................. (6.11)
n-pentane 1.41 1.09 osp
Hexanes 4.33 3.97 !!..V-8
and !!..Rs = V- . . . . . • . . . . . . . • • . • • . • • • • . • • • • • • . (6.12)
Heptanes plus 33.29 77.41 0.8515 34.5 218 o
Total 100.00 100.00 Total GOR is calculated by adding the stock-tank-oil-based GOR's
"At60ºF. from each separator stage.
Nsp
y 0 = stock-tank-oil density; and Yg- = gravity of gas released from Rs = L(!!..Rs)k · · ·........... . .....•......... (6.13)
the separator oil. k=l
Finally, the value of stock-tank-liquid/wellstream ratio in bbl/MMscf
Separator-oil FVF's, Bosp• are reported as the ratio of separator-oil
represents the separator barreis produced per 1 MMscf of wellstream.
volume to stock-tank-oil volume.
ln terms of F8sp and separator properties, this value equals
- vosp
bbJ _ 470(1-Fgsp)(Mosp/Posp) Bosp - V-' ...... .. .. ...... .......... ....... (6.14)
- B , .. .. .. .. .. . .. . (6.10) o
MM se f osp
Accordingly, the relation between separator gas/oil ratio and stock-
where 470= (1 million scf/MMscf)/[(379 scf/lbm mol)(5.615 ft3/bbl)]. tank gas/oil ratio at a given stage is
The separator-oil and -gas compositions can be checked for con-
sistency with the Hoffman et a/. 10 K-value method and Standing' si 1 !!..Rs
!!..Rsp = -¡¡-· .. . .. . ................ . ...... ... (6.15)
Jow-pressure K-value equations. osp
Because Bosp > 1, it follows that Rsp < R,.
6.3 Multistage-Separator Test Bubblepoint-oil FVF, B 0 b, is the ratio of bubblepoint-oil volume
The multistage-separator test is performed on an oil sample primari- to stock-tank-oil volume.
ly to provide a basis for converting differential-liberation data from
V
a residual-oil to a stock-tank=oil basis. Occasionally, severa) separa- B0 b = y~· .......... ....... ................. (6.16)
tor tests are run to help choose separator conditions that maximize o
stock-tank-oil production. Usually, two or three stages of separation The average gas gravity, y8 , is used in oil PVT correlations and
are used, with the Jast stage at atmospheric pressure and near-ambi- to calculate reservoir densities on the basis of black-oil properties.
ent temperature (60 to 80ºF). The multistage-separator test can also The average gas gravity is calculated from
be conducted for high-liquid-yield gas-condensate fluids.
Fig. 6.2 illustrates schematically how the separator test is per- Nsp
formed. Initially, the reservoir sample is at saturation conditions and L(yg)J!!..Rs)k
the volume is measured ( V0 b or V8d). The sample is then brought to y = _k=_l- - -- .......................... (6.17)
the pressure and temperature of the first-stage separator. Ali the gas g Nsp
is removed, and the oil volume at the separator stage, Vosp• is noted L(!!..Rs)k
together with the volume of removed gas, !!.. V8 ; number of moles of k= J
!
1
/ y;
!
/ Wellstream
Composition
,/ X¡ Z;
(/)
I
ID M¡¡
Reservoir Standard
Conditions Conditions
92 PHASE BEHAVIOR
where M; =molecular weight and p; = componen! liquid density
in lbm/ft 3 at standard conditions (Table A-1). The C7+ material in
separator gases is usually treated as normal heptane.
..... .... . . ment. The sample initially is brought to a condition somewhat above
initial reservo ir pressure, ensuring that the fluid is single phase. As the
pressure is lowered, oil volume expands and is recorded.
The fluid is agitated at each pressure by rotating the ce!!. This
avoids the phenomenon of supersaturation, or metastable equilibri-
_ . . . . . ....1111···-::::::::: : ::-
Ps1=14.7 psia
100 psi below actual bubblepoint pressure. By agitating the mixture
at each new pressure, the condition of supersaturation is avoided, al-
lowing more accurate determination of the bubblepoint.
T5 1=60ºF
Just below the bubblepoint, the measured volume will increase
Fig. 6.2-Schematic of a multistage-separator test.
more rapidly because gas evolves from theoil, yielding a higher sys-
tem compressibility. The total vol u me, V,, is recorded after the two-
where (y8)k =separator-gas gravity at Stage k. This relation is based
phase mixture is brought to equilibrium. Pressure is lowered in steps
on the ideal gas law at standard conditions, where moles of gas are di-
of 5 to 200 psi, where equilibrium is obtained at each pressure.
rectly proportional with standard gas volume ( v8 = 379 scf/lbm mol).
Table 6.8 gives the composition of the first-stage-separator gas When the lowest pressure is reached, total volume is three to five
at 50 psig and 75ºF. The gross heating value, H 8, ofthis gas is calcu- times larger than the original bubblepoint volume.
lated by Kay's 12 mixing rule and componen! heating values, H;, The recorded cell volumes are plotted vs. pressure, and the result-
given in Table A-1. ing curve should be similar to one of the curves in Fig. 6.4.16 For a
black oil (far from its critica! temperature), the discontinuity in vol-
N
ume at the bubblepoint is sharp and the bubblepoint pressure and
Hg = LY;H¡ . . ....................... .. . .. (6.18)
i=I
volume are easily read from the intersection of the p-V trends in the
single- and two-phase regions.
Component liquid yields, L;, represent the liquid volumes of a
Volatile oils do not exhibit the same clear discontinuity in volu-
componen! or group of components that can theoretically be pro-
metric behavior at the bubblepoint pressure. Instead, the p-V curve
cessed from 1 Mscf of separator gas (gallons per million standard
is practically continuous in the region of the bubblepoint because
cubic feet). L¡ can be calculated from
the undersaturated-oil compressibility is similar to the effective
= 19.73y;(~J
two-phase compressibility. This makes determining the bubble-
L; ........................... (6.19)
point of volatile oils in a blind cell difficult. Instead, a windowed cell
Separator
Separator Separator Stock-Tank Vol ume Flashed-Gas
Pressure Temperature GORb GQRC Gravity FVFd Facto re Specific
(psia) (ºF) (ft3/bbl) (ft3/bbl) (ºAPI) (bbl/bbl) (bbl/bbl) Gravity
e= _virrP
(ªav,.,) =_viru (ªavro)
7 P
; P > Pb ·7
......... (6.21)
" Saturation pressure.
1 Relativa volume = VIV a in barreis at indicated pressure per barre! at saturation
pressure.
5 1
VOLATILE/[Link] OIL 1
The V,., function at undersaturated conditions may be fit with a se-
cond-degree polynomial , resulting in an explicit relation for under-
GAS saturated-oil compressibility (see Chap. 3) .
OL
... Total volumes below the bubblepoint can be correlated by the Y
function, 16· 17 defined as
y= Pb - P Pb - P
. . . . . . . . . . . . . . (6.22)
6000
paia
BUBBLEPOINT
5000 paio
41100
pala
4500
p•ia
p(Vrr-1)
94 PHASE BEHAVIOR
IH
•• '
.1
P)
E
u r..,
.
cu t
E ¡;¡_
-
:::1
o
>
- u._, + ·- ~
~
t:
H
..L .,
~~ t-
,.,. -
;
-
70 -t:t .I
j h.., '-.-8:_- r.--
._._.... 'T
!:h: ~~
1 ~c.¡:_ .ttt:t :I -:-
' ~ -~Ú
•• H+ 1
-++-• .,.,... .... ·~!-
:t:+i .+,.. H
; ~, -+
., --:·_ ~t:+i H
f!;t:ff.;:
i •.
¡¡~~
1
,+ , ;.
'¡
~r
·H-
~t .... ~- ~+
-·-·--
:H =:1t:±:ti -1 L f-F 1.t~+ rt+d 3:~
IOOO .. 6000 #00
Pressure, psia
Fig. 6.4-Volume vs. pressure for an oil during a OLE test (after Standing16).
tors' are reported at pressures greater than and equal to the dewpoint calculate reservoir performance.16,l9- 21Fig.6.6 illustrates the labo-
pressure. Table 6.10 gives these data for a gas-condensate example. ratory procedure of a OLE experiment. Figs. 6.7A through 6.7C
Reciproca! wet-gas FVF, b8 w, is reported at dewpoint and initial and Table 6.11 give DLE data for an oil sample.
reservoir pressures, where these values represent the gas equivalent A blind cell is filled with an oil sample, which is brought to a
or wet-gas volume at standard conditions produced from 1 bbl of single phase at reservoir temperature. Pressure is decreased until the
reservoir gas volume. fluid reaches its bubblepoint, where the oil volume, V 0 b , is recorded.
Because the initial mass of the sample is known, bubblepoint densi-
bgw = (5.615 X 10- 3) ~;: IT = 0.198 IT' ........ (6.23) ty, p ob• can be calculated.
The pressure is decreased below the bubblepoint, and the cell is
with b 8 w in Mscf/bbl, pin psia, and T in ºR. agitated until equilibrium is reached. Ali gas is removed at constant
Most CCE experiments are conducted in a visual cell for gas con- pressure. Then, the volume, L1 V8 ; moles, Lln 8 ; and specific gravity,
densates, and relative oil (condensate) volumes, Vro, are reported at y 8 , of the removed gas are measured. The remaining oil volume, V 0,
pressures below the dewpoint. V,0 normally is defined as the oil vol- is also recorded. This procedure is repeated 10 to 15 times at de-
ume divided by the total volume of gas and oil, although sorne re- creasing pressures and finally at atmospheric pressure. Residual-oil
ports define itas the oil volume divided by the dewpoint volume. volume, V0 ,, and specific gravity, y 0,, are measured at 60ºF.
Other properties are calculated on the basis of measured data
6.5 Ditferential Liberation Expansion (LlV8 , V0 , Lln , y , V0,, and y 0,), including differential solution
8 8
The DLE experiment is designed to approximate the depletion pro- gas/oil ratio, R,d; differential oil FVF, B0 d; oil density, Po; and gas
cess of an oil reservoir 18 and thereby provide suitable PVT data to Z factor, Z. For Stage k, these properties can be determined from
~
185 2,520 87.92
205 2,615 89.05
"'
~ (
(11)
E-
c.> --
CD
. .,
E
:J
o•
>
·~~ .,.,
.o
~
~,,
~
• ,,
,-,¡._~
., · ~
(6.25)
. ..... .. .. . . .. ...... . . ........... (6.29)
96 PHASE BEHAVIOR
TABLE 6.10-CCE DATA FOR GOOD OIL CO. DIFFERENTIAL LIBERATlON EXPERIMENT
WELL 7 GAS-CONDENSATE SAMPLE
Pressure
(psig) Relative volume
Deviation Factor
z
""r----1
•' L_____J ... ······· .:o
6,000 0.8808 1.144
OL
5,713* 0.8948 1.107**
5,300 0.9158 1.051 BUB!l..EPOINT 4000
5090 pala
5,000 0.9317 1.009 pma
4,200
4,100
0.9832
0.9914
0.895
0.881
= Gas is bled down
to Psc at TA
~flashed)
~:fJ 14.7 60 F
plio
a multistage separator to give the volume ratios, R, and B 0 • This lab- (';; )k = k¿(Llng)j. .. .. .. . . .. . . .. . . .. . . . .. .. (6.35)
oratory procedure is costly and time-consuming and is seldom used. j=I
However, the method is readily incorporated into an equation-of- where j = l corresponds to saturation pressure and (Lln 8 ) 1 = O. The
state (EOS) -based PVT program. initial amount (in moles) of the saturated fluid is known when the cell
is charged. The quantity np/ n is usually reported as cumulative wet
6.6 Constant Volume Depletion gas produced in MMscf/MMscf, which is equivalent to mol/mol.
The CVD experiment is designed to provide volumetric and com- Surface compositions Y¡¡ and x 0 of the removed reservoir gas and
positional data for gas-condensate and volatile-oil reservoirs pro- properties of the removed gas are not reported directly in the labora-
ducing by pressure depletion. Fig. 6.8 shows the stepwise procedure tory report but are recombined to yield the equilibrium gas (well-
of a CVD experiment schematically, and Figs. 6.9A through 6.9D stream) composition, y¡, which also represents the equilibrium gas
and Table 6.12 give CVD data for an example gas-condensate fluid. remaining in the cell. The C 7 + molecular weight of the wellstream,
The CVD experiment provides data that can be used directly by M 8c7+, is backcalculated from measured specific gravity
the reservoir engineer, including (1) a reservoir material balance ( Yw = y8 ) and reservoir-gas composition, y. C7+ specific gravity of
that gives average reservoir pressure vs. recovery of total well- the produced gas, y 8c 7 +,is also reported, but this value is calculated
strearn (wet-gas recovery), sales gas, condensate, and natural gas from a correlation.
liquids; (2) produced-wellstream composition and surface products Knowing the cumulative moles removed and its volume occupied
vs. reservoir pressure; and (3) average oíl saturation in the reservoir as a single-phase gas at the removal pressure, we can calculate the
(liquid dropout and revaporization) that occurs during pressure equilibrium gas Z factor from
depletion. For many gas-condensate reservoirs, the recoveries and pilV
oil saturation vs. pressure data from the CVD analysis closely z 8
= -- . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . (6.36)
iln 8 RT
approximate actual field performance for reservoirs producing by
pressure depletion. When other recovery mechanisms, such as wa- A "two-phase" Z factor is also reported that is calculated assum-
terdrive and gas cycling, are considered, the basic data required for ing that the gas-condensate reservoir depletes according to the ma-
reservoir engineering are still taken mainly from a CVD report. This terial balance for a dry gas and that the initial condition of the reser-
section provides a description of the data provided in a standard voir is at dewpoint pressure.
:e
..Q
¡¡¡ \ZOO
:::1
"C
·¡;;
- ...CD
CJ
11)
.t,
1000
l:C'° \
ri
o 800
o
...eo
:::1
o
(1) óOO '
...
¡¡¡
e
...CDCD
-
i5
400
zoo
Pressure, psig
Fig. 6.7A-DLE data for an oil sample from Good Oil Co. Well 4; differential solution gas/oil
ratio, Rsd·
98 PHASE BEHAVIOR
l. 7
:E
.J:l
m
:::J
1. ó
"O
..
·¡;;
Q)
.....
:E
.J:l l. 5
-l:io
a:i
o
.
·~
ca 1. 4
a:
CD
E
:::J
o
> l. 3
o
m
·~
.
e
-
CD
CI)
i5
l. z.
.,.
l. 1
Pressure. psig
Fig. 6.78-DLE data for an oil sample from Good Oil Co. Well 4; differential oil FVF (relative
volume), Bod·
illustrates the process schematically. The calculated recoveries are B = Bgw = (Psc!Tsc)(ZT/p)
based on multistage-separator calculations that use low-pressure K gd Fgg Fgg
values anda set of separator conditions chosen arbitrarily or speci-
fied when the PVT study is requested. (14.7 /520)[[0.867(186 + 460)]/40151
0.8788
6.6.1 Recoveries: "Normal Temperature Separation." Column
1: lnitial in Place. In Column 1, Row 2a the stock-tank oil in solu- = 4.487 x 10- 3 ft 3 /scf. .. . ... ..... ..... . .. (6.44)
tion in the initial dewpoint fluid (N = 135.7 STB) is calculated by The producing GOR of the dewpoint mixture for the specified
flashing 1 MMscf of the original dewpoint fluid, Gw, through a separator conditions can be calculated as
multistage separator.
Rows 2b through 2d give the volumes of separator gas at each RP = ~ = [(757.87 + 96.68 + 24.23 Mscf/lbm mol)
stage of a three-stage flash of the initial dewpoint fluid: 757.87,
96.68, and 24.23 Mscf, respectively. The mole fraction of well- x (1 x 10 3 scf/Mscf)J/135.7 STB/lbm mol
stream resulting as a surface gas Fgg is given by
= 6, 476 scf/STB. (6.45)
G
F-gg = Gd = (757 .87 + 96.68 + 24.23 Mscf/lbm mol) The dewpoint solution oil/gas ratio, rsd• is simply the inverse of Rp.
w
l. 4
l. 2.
1.0
Q.
.
CJ
~
\
\
·¡¡; 0.8
oCJ
U)
>
o 0.6 1
Saturation
Pressure
0.4
o.z
o.o
o 1000 zoo o 3000 4000 5000 6000
Pressure. psig
Fig. 6.7C-DLE data for an oil sample from Good Oil Co. Well 4; oil viscosity, µ 0 •
Column 2 and Higher. On the basis of 1 MMscf of initial dew- At a given pressure, the mole fraction of the removed CVD gas
point fluid, Rows 2a through 2d give cumulative volumes of separa- wellstream that becomes dry separator gas is given by
tor products at each depletion pressure (Np, GP 1' GP 2' and GP 3).
+ GP1 + + GP 1 +
(F88 t
The producing GOR of the wellstream produced during a depletion (GP 1 GP3) - (GP 1 GP3)
= k . k-1
stage is given by Gw [(np/n)k - (np/n)k_ 1]
(GP 1 + GP2 + GP 3) - (GP 1 + GP2 + GP 3) (6.50)
(RP)k = k k-1
(NP)k - (NP)k-1 For p = 2, 100 psig, this gives
i=J
600 0.654
0.783
0.0126
0.0121
'1J
350 J'n"3tlfi
159
o
0.855
1.286
0.0114
0.0093
Vi
Gravity of residual oil = 35.1 ºAPI at 60ºF.
"Barreis o! residual oi l.
00
At60º F.
Mscf dry gas originally in place, or 83.4%. These recoveries can be CONSTANT VOLUME DEPLE'TlON
compared with the reported wet-gas (or molar) recovery of76. 787%
at 605 psig. In addition to recoveries, the calculated results in this
section can be used to calculate solution oil/gas ratio, r,, and dry-gas .......·· /c=::J . ··· ··o ···:c=J
FVF, B8d, for modified black-oil applications.
are for the calculated separator gases from Stages 1 and 2 of the
.. . . .~·· · LJ.::::::::::::c=J
three-stage surface separation. Recoveries in Row 5 are for the pro-
.. .······· · ··· o
~. ... •
duced wellstreams from the CVD experiment and represent the ab-
solute maximum liquid recoveries that can be expected if the reser-
voir is produced by pressure depletion. Fig. 6.10 illustrates the
recovery calculations schematically. Liquid volumes (in gal/MMscf .......,,, , , ,
of initial dewpoint fluid) at CVD Stage k are calculated from
*.. ~
?:."
o
> 38
>e
.;
:::i 25
o
a.
o
Q
"C
·:; 2B
a"
:::¡
o
o>
15
..:;
U)
U)
l!!
Q.
e:
·c;
111
e
..
Q.
3::
11
B 500 11100 1500
Pressure, psig
Fig. 6.9A-CVD data for gas-condensate sample from Good Oil Co. Well 7; liquid-dropout curve, V,.0 .
Fig. 6.9B-CVD data for gas-condensate sample from Good Oil Co. Well 7; equilibrium gas
compositions, y¡.
1 at dewpoint pressure, and Q as actual cumulative recoveries based should be increased by 117 .3 Mscf for ali depletion pressures in the
on hydrocarbons in place at the initial pressure, CVD table.
For stock-tank-oil recovery, Qd = 135.7 STB, so liQd = 15.9
STB. Stock-tank-oil recovery at 4,000 psig is 15.9 + O = 15.9
STB; at 3,500 psig the recovery should be 15.9 + 6.4 = 22.3 STB,
and so on.
On the basis of 1 MMscf wet gas at the dewpoint or 1.1173 MMscf
Q = QTable + f1Qd; P < Pd• · · · · · · · · • · · • · · · · · · · (6.57) at initial reservoir pressure, the laboratory hydrocarbon pore vol-
ume (HCPV), VpHCJab• is the same.
(p/Z); ]
and liQd = Qd (p/Z)d - 1 , . . . . . . . . . . . . . . . . . . . . (6.58)
[ VpHClab = ([Link])d
l. 409
e
:s
l.308
.."'
a:"'
~
·5
..."'
e:
..
N l. 2mJ
s
a:
...u
u.. e
e :s
o "'"'e
·~ 1.1ee
·s:
CD
a.
'E
o ·5
o
"'... .
Q,
Q
~
l. 082
il.809
il. 71!0
QI .
Pressure. psig
Fig. 6.9C-CVD data for gas-condensate sample from Good Oil Co. Well 7; equilibrium gas Z
factor,~ ·
VpHCactual reservoir. Correcting for water saturation, Sw, the reservoir-oil satu-
Qacrual = Qlab V • · · · · · · · · · · · · · · · · · · · . . . . (6.61) ration can be calculated from V,0 with
pHClab
So = (1 - Sw)Vro . . . . . . . . . . . . . . . . . . . . . . . . . . . . (6.64)
where Q 1ab = laboratory value given by Eqs. 6.55 and6.57. As anex-
ample, suppose geological data indicate a HCPV of 625,000 bbl For most gas condensates, V,0 shows a maximum near 2,000 to
(82.45 acre-ft), or 3.509 x 106 ft 3 . Then, original wet gas in place is 2,500 psia. Cho et al. 27 give a correlation for maximum liquid drop-
out as a function of temperature and C 7 + mole percent in the dew-
106
G = l 1173 X 10 6 3 ·509 X point mixture.
w • 3, 943
(V, )max = 93.404 + 4.799zc7+ - 19.73lnT, ..... . (6.65)
0
98
'#-
j
<.:I
~a. 7iJ
<.:I
.¿
CD
u
:::1
"C
o 88
d:
~
.."'
<.:I
QI Sil
--
! U)
E
..
"'...
GI
.!!
-¡
3:
Cll
48
::>
>
º+:: [Link]
Cll
:;
38
e(
E
:::1 e"'"'
a.
CJ
28
~
Ul
a .J
a 1aaa !Sila 2aila
Pressure, psig
25ilZ •• 00
Fig. 6.90-CVD data for gas-condensate sample from Good Oil Co. Well 7; wet-gas material
balance.
tude of liquid dropout reflects the change in producing oil/gas ratio, should be emphasized. In fact, the effect of shape and magnitude of
rp = r,. A tail on a liquid-dropout curve implies that the producing liquid dropout on fluid flow in the reservoir is negligible, and any
wellstream is becoming only slightly leaner (i.e., r, is decreasing EOS match will probably have the same effect on fluid flow from the
only slightly). The cumulative condensate recovery is given by reservoir into the wellbore (i.e., inflow performance).
Gp
6.6.S Consistency Check of CVD Data. Reudelhuber and Hinds24
NP = f r,dGP, ... .. . . ....... ... .. ..... .. .... (6.66) give a detailed procedure for checking CVD data consistency that
involves a material-balance check on components and phases and
o
yields oil compositions, density, molecular weight, and Me?+. To-
where Gp = cumulative dry gas produced. Cumulative condensate gether with reported data, these calculated properties allow K values
production is readily evaluated from a plot of r, vs. Gp. to be calculated and checked for consistency with the Hoffman et
One of the most important checks of an EOS characterization for aL.10 method.11,28 Whitson and Torp's23 material-balance equations
any gas condensate, particularly one with a tail, is NP calculated
are summarized later. Similar equations can also be derived for a
from CVD data vs. Np calculated from the EOS characterization. It
DLE experiment when equilibrium gas compositions and oil rela-
is alarming how much the surface condensate recovery can be un-
tive volumes are reported. Reported CVD data include temperature,
derestimated if the tail is not matched properly. We do not recom-
mend matching the dewpoint exactly with a liquid-dropout curve T; dewpoint pressure, p d• or bubblepoint pressure, p b; dewpoint Z
that is severely overpredicted in the region where measured results factor, Zd, or bubblepoint-oil density, p ob. Additional data at each
indicate little dropout. If the EOS characterization cannot be modi- Depletion Stage k include oil relative volume, V,0 ; initial fraction
fied to honor the tail of liquid-dropout curve, it is preferable to of cumulative moles produced, np/n; gas Z factor (not the two-
underpredict the measured dewpoint pressure and match only the phase Z factor), Z; equilibrium gas composition, y;; and equilibrium
higher liquid-dropout volumes. gas (wellstream) C7+ molecular weight, M,c,,.
In summary, oil relative volume, V,0 , is not importantper se; how- The equilibrium gas density, p 8 ; molecular weight, M8 ; and well-
ever, the effect of liquid dropout on surface condensate production stream gravity, Yw = M8 /Mair, are readily calculated at each
(6.68) and (11 0 )k = (11 ,)k - (n 8)k, ......... . . . .. .. ... . ... . (6.73)
and moles and mass of the individual phases remaining in the cell at
(Mgt
and (yg)k = (y .. )k = 28 .97. . ............. . .... . (6 .69) Stage k are given by
M, - ~( Ll:g)}Mg)1.
1
On a basis of 1 mol initial dewpoint fluid (11 = 1), the cell vol-
ume is (m,)k =
ZdRT
vcell = ---¡¡;¡-- . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . (6.70)
(mg)k = (1~g)JMg)k.
for a gas condensate and
M,,b
Vce11 = Pub .......... . . - .. . ... . .. .. ... ........ (6.71)
.. .. . .......... . ...... /
In Eqs. 6.73 and 6.74,
for a volatile oil. Oil and gas volumes, respectively, at Stage k are
Moles and mass of the total material remaining in the cell at Stage k M, = saturated-fluid molecular weight, and (np/11) 1 = O.
are given by Densities and molecular weights of the oil phase are calculated from
K values can be calculated from K; = yjx;. and z; = overall com- Table 6.6 summarizes these calculations for the sample gas-conden-
position of the mixture remaining in the cell at Stage k. sate mixture.
OR
L¡1 ¡,; = liquid yields
(Separator Gas 1)
Separator 1 c,.(gas)
o-o
¡···~~;;,o~;d···¡
Separator 2
(tank gas) } no liquid yields
Fig. 6.10-Schematic of method of calculating plant recoveries in a evo report for a gas
condensate.
-20+---"'-......,.::_-+-~~~-+-~~~--r-~~~-+~~~-+~~~-+~~~-t~~~---t
o 2 4 6 8 10 12 14 16
Heptanes Plus, mol%
Fig. 6.11-Calculated maximum retrograde oil relative vol u mes trom the Cho et a/. 27 correlation.
The oil composition at the last depletion state (605 psig forthe ex- 14. Silvey, F.C., Reamer, H.H., and Sage, B.H.: "Supersaturation in Hydrocar-
ample condensate) can be measured, but it must be requested specif- bon Systems: Methane-n-Decane," Ind. Eng. Chem (1958) 3, No. 2, 181.
15. Tindy, R. and Raynal, M.: "Are Test-Cell Saturation Pressures Accurate
ically. Also, the residual-oil molecular weight, Mor, and specific
Enough?," Oil & Gas J. (December 1966) 126.
gravity, yor, remaining after depletion at atmospheric pressure are 16. Standing, M .B .: Volumetric and Phase Behavior of Oil Field Hydrocar-
typically measured and reported as shown in Table 6.12. These val- bon Systems, eighth edition, SPE, Richardson, Texas (1977).
ues can be compared with calculated values by use ofthe material- 17. Clark, N.J.: "Adjusting Oil Sample Data for Reservoir Studies," JPT
balance equations shown earlier. (February 1962) 143.
The material-balance calculations are more accurate for rich gas 18. Moses, P.L.: "Engineering Applications of Phase Behavior of Crude-Oil
and Condensate Systems," JPT (July 1986) 715.
condensates and volatile oils. In fact, obtaining reasonable material- 19. Amyx, J.W., Bass, D .M. Jr., and Whiting, R.L.: PetroleumReservoir En-
balance oil properties for lean gas condensates is difficult. Some- gineering , McGraw-Hill Book Co. Inc., New York City (1960).
times it is useful to modify the reported oil relative volumes (partic- 20. Craft, B .C . and Hawkins, M.: Applied PetroleumReservoir Engineering,
ularly those close to the dewpoint) to monitor the effect on first edition, Prentice-Hall Inc., Englewood Cliffs, New Jersey (1959) .
calculated oíl properties. 21. Dake, L.P.: Fundamentals of Reservoir Engineering, Elsevier Scientific
Publishing Co., Amsterdam (1978).
An altemative material-balance check that may be even more 22. Dodson, C .R., Goodwill, D. , andMayer, E.H.: "Application ofLabora-
useful for determining data consistency (particularly for leaner gas tory PVT Data to Reservoir Engineering Problems," Trans., AIME
condensates) involves starting with reported final-stage condensate (1953) 198, 287 .
composition, (x¡)k=N• and adding back the removed gases, (y¡h, for 23. Whitson, C.H. and Torp, S.B .: "Evaluating Constant-Volume-Depletion
each stage from k = N to k = l. This results in the original gas Data," JPT (March 1983) 610; Trans., AIME, 275.
24. Drohm, J.K., Goldthorpe, W.H. , and Trengove, R.: "Enhancing the Eval-
composition, (z¡)k=t• which can be compared quantitatively with uation of PVT Data," paper OSEA 88174 presented at the 1988 Offshore
the laboratory-reported composition. Southeast Asia Conference, Singapore, 2-5 February.
25. Drohm, J.K., Trengove, R., and Goldthorpe, W.H.: "On the Quality of
References Data From Standard Gas-Condensate PVT Experiments," paper SPE
17768 presented at the 1988 Gas Technology Symposium, Dallas,
1. "Core Laboratories Good Oil Company Oil Well No . 4 PVT Study," Core
13-15 June.
Laboratories, Houston.
26. Reudelhuber, F.O. and Hinds, R.F.: "Compositional Material Balance
2. "Core Laboratories Good Oil Company Condensate Well No. 7 PVT
Method for Prediction of Recovery From Volatile-Oil Depletion-Drive
Study," Core Laboratories, Houston.
Reservoirs," JPT (January 1957) 19; Trans., AIME, 210.
3. Flaitz, J.M. and Parks, A .S.: "Sampling Gas-Condensate Wells," Trans. ,
27. Cho, S.J., Civan, F., and Starling, K.E.: "A Correlation To PredictMaxi-
AIME (1942) 146, 13. mum Condensation for Retrograde Condensation Fluids and Its Use in
4. Katz, D.L., Brown, G.G ., and Parks, A.S .: "NGAA Report on Sampling Pressure-Depletion Calculations," paper SPE 14268 presented at the
Two-Phase Gas Streams from High Pressure Condensate Wells," (Sep- 1985 SPE AnnualTechnical Conference and Exhibition, Las Vegas, Ne-
tember 1945). vada, 22-25 September.
5. Reudelhuber, F.O.: "Sampling Procedures forOil Reservoir Fluids," JPT 28. Clark, N.J.: "Theoretical Aspects of Oil and Gas Equilibrium Calcula-
(December 1957) 15. tions," JPT (April 1962) 373.
6. Clark, N.J.: "Sampling and Testing Oil Reservoir Samples," JPT (Jan.
1962) 12.
SI M etric Conversion Factors
7. Clark, N.J .: "Sampling and Testing Gas Reservoir Samples," JPT
(March 1962) 266. º API 141.5/(131.5 + º API) =g/cm3
8. Recommended Practice for Sampling Petroleum Reservoir Fluids, API, bbl X 1.589 873 E- 01 = m3
Dallas (1966) 44. Btu x 1.055 056 E+OO =kJ
9. Standing, M .B. and Katz, D.L. : "Density of Natural Gases," Trans. ,
cp X 1.0* E - 03 = Pa · s
AIME, (1942) 146, 140.
10. Hoffmann, A.E., Crump, J.S., and Hocott, C .R. : "Equilibrium Constants ft X3.048* E-01 =m
for a Gas-Condensate System," Trans., AIME (1953) 198, l. ft3 X2.831685 E-02 =m3
11. Standing, M.B.: "A Set of Equations for Computing Equilibrium Ratios ºF ('F-32)/1.8 = ºC
of a Crude Oil/Natural Gas System at Pressures Below 1,000 psia," JPT gal x3.785 412 E-03 =m3 •
(September 1979) 1193. in. x 2.54* E+OO =cm
12. Kay, W.B.: "The Ethane-Heptane System," Ind. & Eng. Chem. (1938) lbmmol x4.535 924 E-01 =kmol
30, 459.
13. Kennedy, H.T. and Olson, C.R .: "Bubble Formation in Supersaturated
psi x 6.894 757 E+ 00 = kPa
Hydrocarbon Mixtures," Oil & Gas J. (October 1952) 271. 'Conversion factor is exact.