Chapter 3

Laboratory Analysis
of Reservoir Fluids
Accurate laboratory studies of PVT and phase-equilibria behavior of reservoir
fluids are necessary for characterizing these fluids and evaluating their
volumet- ric performance at various pressure levels. There are many
laboratory analyses that can be made on a reservoir fluid sample. The amount
of data desired deter- mines the number of tests performed in the laboratory.
In general, there are three types of laboratory tests used to measure
hydrocarbon reservoir samples:
1. Primary tests
These are simple, routine field (on-site) tests involving the
measurements of the specific gravity and the gas-oil ratio of the produced
hydrocarbon fluids.
2. Routine laboratory tests
These are several laboratory tests that are routinely conducted to char-
acterize the reservoir hydrocarbon fluid. They include:
○ Compositional analysis of the system
○ Constant-composition expansion
○ Differential liberation
○ Separator tests
○ Constant-volume depletion
3. Special laboratory PVT tests
These types of tests are performed for very specific applications. If a
res- ervoir is to be depleted under miscible gas injection or a gas cycling
scheme, the following tests may be performed:
○ Slim-tube test
○ Swelling test
The objective of this chapter is to review the PVT laboratory tests and to illus-
trate the proper use of the information contained in PVT reports.

COMPOSITION OF THE RESERVOIR FLUID

It is desirable to obtain a fluid sample as early in the life of a field as possible
so that the sample will closely approximate the original reservoir fluid.
Collection of a fluid sample early in the life of a field reduces the chances of
free gas exist- ing in the oil zone of the reservoir.
Reservoir Engineering Handbook. https://doi.org/10.1016/B978-0-12-813649-2.00003-7
© 2019 Elsevier Inc. All rights reserved. 123
TABLE 3-1 Hydrocarbon Analysis of Reservoir Fluid Sample

Composition of Reservoir Fluid Sample

(by Flash, Extended-Capillary Chromatography)
Liquid
Component Name Mol % Wt % Density MW
(gm/cc)
Hydrogen Sulfide 0.00 0.00 0.800 34.08
6
Carbon Dioxide 0.25 0.11 0.817 44.01
2
Nitrogen 0.88 0.25 0.808 28.013
6
Methane 23.94 3.82 0.299 16.043
7
Ethane 11.67 3.49 0.356 30.07
2
Propane 9.36 4.11 0.507 44.097 Total Sample Properties
0
iso-Butane 1.39 0.81 0.562 58.123 Molecular Weight……100.55
9
n-Butane 4.61 2.66 0.584 58.123 Equivalent Liquid Density,
0 gm/scc ..0.7204
iso-Pentane 1.50 1.07 0.624 72.15
4
n-Pentane 2.48 1.78 0.631 72.15
1
Hexanes 3.26 2.73 0.685 84
0
Heptanes 5.83 5.57 0.722 96
0
Octanes 5.52 5.88 0.745 107
0
Nonanes 3.74 4.50 0.764 121
 0
Decanes 3.38 4.50 0.778 134
0
Undecanes 2.57 3.76 0.789 147
0
Dodecanes 2.02 3.23 0.800 161
0
Tridecanes 2.02 3.52 0.811 175
0
Tetradecanes 1.65 3.12 0.8220 190
Pentadecanes 1.48 3.03 0.8320 206
Hexadecanes 1.16 2.57 0.8390 222
Heptadecanes 1.06 2.50 0.8470 237
Octadecanes 0.93 2.31 0.8520 251
Nonadecanes 0.88 2.31 0.8570 263
Eicosanes 0.77 2.11 0.8620 275
Heneicosanes 0.68 1.96 0.8670 291
Docosanes 0.60 1.83 0.8720 305

Tricosanes 0.55 1.74 0.8770 318 Plus Fractions Mol% Wt% Density MW
Tetracosanes 0.48 1.57 0.8810 331 Heptanes plus 40.66 79.17 0.8494 196
Pentacosanes 0.47 1.60 0.8850 345 Undecanes plus 22.19 58.72 0.8907 266
Hexacosanes 0.41 1.46 0.8890 359 Pentadecanes plus 13.93 45.09 0.9204 326
Heptacosanes 0.36 1.33 0.8930 374 Eicosanes plus 8.42 32.37 0.9540 387
Octacosanes 0.37 1.41 0.8960 388 Pentacosanes plus 5.34 23.16 0.9916 437
Nonacosanes 0.34 1.34 0.8990 402 Triacontanes plus 3.39 16.02 1.0440 474
Triacontanes plus 3.39 16.02 1.0440 474

Totals 100.00 100.00

126 Reservoir Engineering Handbook

Most of the parameters measured in a reservoir fluid study can be

calculated with some degree of accuracy from the composition. It is the most
complete description of reservoir fluid that can be made. In the past, reservoir
fluid com- positions were usually measured to include separation of the
component meth- ane through hexane, with the heptanes and heavier
components grouped as a single component reported with the average
molecular weight and density.
With the development of more sophisticated equations-of-state to calculate
fluid properties, it was learned that a more complete description of the heavy
components was necessary. It is recommended that compositional analyses
of the reservoir fluid should include a separation of components through C 10
as a minimum. The more sophisticated research laboratories now use
equations-of-state that require compositions through C30 or higher.
Table 3-1 shows a chromatographic “fingerprint” compositional analysis
of the Big Butte crude oil system. The table includes the mole fraction, weight
fraction, density, and molecular weight of the individual component.

CONSTANT-COMPOSITION EXPANSION TESTS

Constant-composition expansion experiments are performed on gas conden-
sates or crude oil to simulate the pressure-volume relations of these hydrocar-
bon systems. The test is conducted for the purposes of determining:
○ Saturation pressure (bubble-point or dew-point pressure)
○ Isothermal compressibility coefficients of the single-phase fluid in excess
of saturation pressure
○ Compressibility factors of the gas phase
○ Total hydrocarbon volume as a function of pressure
The experimental procedure, as shown schematically in Figure 3-1 involves
placing a hydrocarbon fluid sample (oil or gas) in a visual PVT cell at
reservoir temperature and at a pressure in excess of the initial reservoir
pressure (Figure 3-1, Section A). The pressure is reduced in steps at constant
temperature by removing mercury from the cell, and the change in the total
hydrocarbon volume Vt is measured for each pressure increment.
The saturation pressure (bubble-point or dew-point pressure) and the
corre- sponding volume are observed and recorded and used as a reference
volume Vsat (Figure 3-1, Section C). The volume of the hydrocarbon system as
a function of the cell pressure is reported as the ratio of the reference volume.
This volume is termed the relative volume and is expressed mathematically by
the following equation:
Vt
Vrel ¼ V (3-1)
sat

where
Vrel ¼ relative volume
Vt ¼ total hydrocarbon volume
Vsat ¼ volume at the saturation pressure
 Laboratory Analysis of Reservoir Fluids Chapter 3
127

FIGURE 3-1 Constant-composition expansion test.

The relative volume is equal to one at the saturation pressure. This test is
com- monly called pressure-volume relations, flash liberation, flash
vaporization, or flash expansion.
It should be noted that no hydrocarbon material is removed from the cell,
thus, the composition of the total hydrocarbon mixture in the cell remains
fixed at the original composition.
Table 3-2 shows the results of the flash liberation test (the constant compo-
sition expansion test) for the Big Butte crude oil system. The bubble-point
pres- sure of the hydrocarbon system is 1930 psi at 247°F. In addition to the
reported values of the relative volume, the table includes the measured values
of the oil density at and above the saturation pressure.
The density of the oil at the saturation pressure is 0.6484 gm/cc and is
deter- mined from direct weight-volume measurements on the sample in the
PVT cell. Above the bubble-point pressure, the density of the oil can be
calculated by using the recorded relative volume:
ρ ¼ ρsat=Vrel (3-2)
where
ρ ¼ density at any pressure above the saturation pressure
ρsat ¼ density at the saturation pressure
Vrel ¼ relative volume at the pressure of interest

Example 3-1
Given the experimental data in Table 3-2, verify the oil density values at 4,000
and 6,500 psi.

Solution
Using Equation 3-2 gives:
TABLE 3-2 Constant Composition Expansion Data

Pressure-Volume Relations (at 247°F)

Relative Y- Densit
Pressure psig Volume (A) Function y
(B) gm/cc
6500 0.9371 0.6919
6000 0.9422 0.6882
5500 0.9475 0.6843
5000 0.9532 0.6803
4500 0.9592 0.6760
4000 0.9657 0.6714
3500 0.9728 0.6665
3000 0.9805 0.6613
2500 0.9890 0.6556
2400 0.9909 0.6544
2300 0.9927 0.6531
2200 0.9947 0.6519
2100 0.9966 0.6506
2000 0.9987 0.6493
b≫ 1936 1.0000 0.6484
1930 1.0014
1928 1.0018
1923 1.0030
1918 1.0042
1911 1.0058
1878 1.0139
1808 1.0324
1709 1.0625 2.108
1600 1.1018 2.044
1467 1.1611 1.965
1313 1.2504 1.874
1161 1.3694 1.784
1035 1.5020 1.710
782 1.9283 1.560
600 2.4960 1.453
437 3.4464 1.356

(A) Relative volume: V/Vsat or volume at indicated pressure per volume at saturation pressure.
ðpsat — pÞ
(B) Where Y-function
ðpabs Þ · ðV =Vsat — 1Þ
 Laboratory Analysis of Reservoir Fluids Chapter 3
129
○ At 4,000 psi
0:6484
ρo ¼ ¼ 0:6714 gm=cc
0:9657
○ At 6500 psi
0:6484
ρo ¼ ¼ 0:6919
0:9371
The relative volume data frequently require smoothing to correct for
laboratory inaccuracies in measuring the total hydrocarbon volume just below
the satura- tion pressure and also at lower pressures. A dimensionless
compressibility func- tion, commonly called the Y-function, is used to
smooth the values of the relative volume. The function in its mathematical
form is only defined below the saturation pressure and is given by the
following expression:

psat — p (3-3)
Y ¼ pðVrel —
1Þ
where
psat ¼ saturation pressure, psia
p ¼ pressure, psia
Vrel ¼ relative volume at pressure p
Column 3 in Table 3-2 lists the computed values of the Y-function as
calculated by using Equation 3-3. To smooth the relative volume data below
the saturation pressure, the Y-function is plotted as a function of pressure on a
Cartesian scale. When plotted, the Y-function forms a straight line or has only
a small curvature. Figure 3-2 shows the Y-function versus pressure for the Big
Butte crude oil sys- tem. The figure illustrates the erratic behavior of the data
near the bubble-point pressure.
The following steps summarize the simple procedure of smoothing and
cor- recting the relative volume data:
Step 1. Calculate the Y-function for all pressures below the saturation
pressure by using Equation 3-3.
Step 2. Plot the Y-function versus pressure on a Cartesian scale.
Step 3. Determine the coefficients of the best straight fit of the data, or:

Y ¼ a + bp (3-4)

where a and b are the intercept and slope of the lines, respectively.
Step 4. Recalculate the relative volume at all pressure below the saturation
pressure from the following expression:

Vrel ¼ 1+ psat — p (3-5)

pða + bpÞ
130 Reservoir Engineering
Handbook

FIGURE 3-2 Y-function versus pressure.

Example 3-2
The best straight fit of the Y-function as a function of pressure for the Big
Butte oil system is given by:
where
Y ¼ a + bp
a ¼ 1.0981
b ¼ 0.000591
Smooth the recorded relative volume data of Table 3-2.
 Laboratory Analysis of Reservoir Fluids Chapter 3
131
Solution
Smoothed Vrel
Pressure Measured Vrel Equation 3-5
193 — —
6
193 — 1.001
0 4
192 — 1.001
8 8
192 — 1.003
3 0
191 — 1.004
8 2
191 — 1.005
1 8
187 — 1.013
8 9
180 — 1.032
8 4
170 1.0625 1.063
9 0
160 1.1018 1.102
0 8
146 1.1611 1.162
7 6
131 1.2504 1.253
3 2
116 1.3696 1.374
1 1
103 1.5020 1.509
5 1
782 1.9283 1.945
8
600 2.4960 2.532
8
437 3.4464 3.529
0

The oil compressibility coefficient co above the bubble-point pressure is

also obtained from the relative volume data as listed in Table 3-3 for the Big
Butte oil system.
The oil compressibility is defined by Equations 2-94 through 2-96 and
equivalently can be written in terms of the relative volume, as:
—1 ∂Vrel
c¼ (3-6)
o Vrel ∂p
Commonly, the relative volume data above the bubble-point pressure is
plotted as a function of pressure as shown in Figure 3-3. To evaluate co at
any pressure p, it is only necessary to graphically differentiate the curve by
drawing a tangent line and determining the slope of the line, i.e., ∂Vrel/∂p.

Example 3-3
Using Figure 3-3, evaluate co at 3,000 psi.
132 Reservoir Engineering
Handbook
Solution
○ Draw a tangent line to the curve and determine the slope.
∂Vrel=∂p ¼ —14:92 × 10—6
○ Apply Equation 3-6 to give
. Σ
—1 . Σ
c ¼o —14:92 × 10—6 ¼ 15:23 × 10—6 psi—1
0:98
 Laboratory Analysis of Reservoir Fluids Chapter 3
133

TABLE 3-3 Undersaturated Compressibility Data

Volumetric Data
(at 247°F)
Saturation Pressure (Psat)......................1936 psig
Density at Psat...................................... 0.6484 gm/cc
Thermal Exp @ 6500 psig....................1.10401 V at 247°F/V at 60°F
Average Single-Phase Compressibilities
Single-Phase
Pressure Compressibili
Range ty v/v/psi
psig
6500 to 6000 10.73 E-6
6000 to 5500 11.31 E-6
5500 to 5000 11.96 E-6
5000 to 4500 12.70 E-6
4500 to 4000 13.57 E-6
4000 to 3500 14.61 E-6
3500 to 3000 15.86 E-6
3000 to 2500 17.43 E-6
2500 to 2000 19.47 E-6
2000 to 1936 20.79 E-6

It should be noted that Table 3-3 lists the compressibility coefficient at several
ranges of pressure, e.g. 6,500–6,000. These values are determined by calculat-
ing the changes in the relative volume at the indicated pressure interval and
evaluating the relative volume at the lower pressure, or
1
c ðVrelÞ1 —
o ¼ ½Vrel ]2 (3-7)
ðVrelÞ2
p1 — p2
where the subscripts 1 and 2 represent the corresponding values at the
higher and lower pressure range, respectively.

Example 3-4
Using the measured relative volume data in Table 3-2 for the Big Butte crude
oil system, calculate the average oil compressibility in the pressure range of
2,500 to 2,000 psi.
134 Reservoir Engineering
Handbook

FIGURE 3-3 Relative volume data above the bubble-point pressure.

Solution
Apply Equation 3-7 to give
—1 0:9890 — 0:9987 —6 —1
co¼ 0:9987 2500 — 2000 ¼ 19:43 × 10 psi

DIFFERENTIAL LIBERATION (VAPORIZATION) TEST

In the differential liberation process, the solution gas that is liberated from an
oil sample during a decline in pressure is continuously removed from contact
with
 Laboratory Analysis of Reservoir Fluids Chapter 3
135
the oil, and before establishing equilibrium with the liquid phase. This type of
liberation is characterized by a varying composition of the total hydrocarbon
system.
The experimental data obtained from the test include:
○ Amount of gas in solution as a function of pressure
○ The shrinkage in the oil volume as a function of pressure
○ Properties of the evolved gas including the composition of the liberated gas,
the gas compressibility factor, and the gas specific gravity
○ Density of the remaining oil as a function of pressure
The differential liberation test is considered to better describe the separation
process taking place in the reservoir and is also considered to simulate the
flow- ing behavior of hydrocarbon systems at conditions above the critical gas
satu- ration. As the saturation of the liberated gas reaches the critical gas
saturation, the liberated gas begins to flow, leaving behind the oil that
originally contained it. This is attributed to the fact that gases have, in general,
higher mobility than oils. Consequently, this behavior follows the differential
liberation sequence. The test is carried out on reservoir oil samples and
involves charging a visual PVT cell with a liquid sample at the bubble-point
pressure and at reservoir tem- perature. As shown schematically in Figure 3-4,
the pressure is reduced in steps, usually 10 to 15 pressure levels, and all the
liberated gas is removed and its vol- ume is measured at standard conditions.
The volume of oil remaining VL is also measured at each pressure level. It
should be noted that the remaining oil is sub- jected to continual compositional
changes as it becomes progressively richer in
the heavier components.
The above procedure is continued to atmospheric pressure where the
volume of the residual (remaining) oil is measured and converted to a volume
at 60°F, Vsc. The differential oil formation volume factors Bod (commonly
called the

FIGURE 3-4 Differential vaporization test at constant temperature “T.”

136 Reservoir Engineering
Handbook

relative oil volume factors) at all the various pressure levels are calculated by
dividing the recorded oil volumes VL by the volume of residual oil Vsc, or:
VL
Bod ¼ (3-8)

Vsc
The differential solution gas-oil ratio R sd is also calculated by dividing the
volume of gas in solution by the residual oil volume.
Table 3-4 shows the results of a differential liberation test for the Big
Butte crude. The test indicates that the differential gas-oil ratio and differential
relative oil volume at the bubble-point pressure are 933 scf/STB and 1.730
bbl/STB, respectively. The symbols Rsdb and Bodb are used to represent these
two values, i.e.:
Rsdb ¼ 933 scf=STB and Bodb ¼ 1:730 bbl=STB

TABLE 3-4 Differential Liberation Data

Differential Vaporization (at 247°F)
Gas
Solutio Relativ Relativ Formatio Increment
n e Oil e Total Oil Deviatio n al Gas
Pressur Gas/Oil Volume Volum Densit n Factor Volume Gravity
e Ratio Bod(B) e y Z Factor (Air 5
psig Rsd(A) Btd(C) gm/cc (D) 1.000)
b≫ 1936 933 1.730 1.730 0.6484

1700 841 1.679 1.846 0.6577 0.864 0.01009 0.885

1500 766 1.639 1.982 0.6650 0.869 0.01149 0.894

1300 693 1.600 2.171 0.6720 0.876 0.01334 0.901

1100 622 1.563 2.444 0.6790 0.885 0.01591 0.909

900 551 1.525 2.862 0.6863 0.898 0.01965 0.927

700 479 1.486 3.557 0.6944 0.913 0.02559 0.966

500 400 1.440 4.881 0.7039 0.932 0.03626 1.051

300 309 1.382 8.138 0.7161 0.955 0.06075 1.230

185 242 1.335 13.302 0.7256 0.970 0.09727 1.423

120 195 1.298 20.439 0.7328 0.979 0.14562 1.593

0 0 1.099 0.7745 2.375

@ 60°F ¼ 1.000

Gravity of residual oil ¼ 34.6°API at 60°F

Density of residual oil ¼ 0.8511 gm/cc at 60°F
(A) Cubic feet of gas at 14.73 psia and 60°F per barrel of residual oil at 60°F.
(B) Barrels of oil at indicated pressure and temperature per barrel of residual oil at 60°F.
(C) Barrels of oil plus liberated gas at indicated pressure and temperature per barrel of residual
oil at 60°F.
(D) Cubic feet of gas at indicated pressure and temperature per cubic feet at 14.73 psia and 60°F.
 Laboratory Analysis of Reservoir Fluids Chapter 3
137

Column C of Table 3-4 shows the relative total volume Btd from differential
liberation as calculated from the following expression:

Btd ¼ Bod + ðRsdb — RsdÞ Bg (3-9)

where
Btd ¼ relative total volume, bbl/STB
Bg ¼ gas formation volume factor, bbl/scf
The gas deviation z-factor listed in column 6 of Table 3-4 represents the z-
factor of the liberated (removed) solution gas at the specific pressure and
these values are calculated from the recorded gas volume measurements as
follows:
. Σ. Σ
Vp Tsc
z¼ (3-10)
T Vscpsc

where
V ¼ volume of the liberated gas in the PVT cell at p and T
Vsc ¼ volume of the removed gas at standard column 7 of Table 3-4 con-
tains the gas formation volume factor B g as expressed by the following
equation:
. Σ
psc zT
B ¼ (3-11)
g
Tsc p
where
Bg ¼ gas formation volume factor, ft3/scf
T ¼ temperature, °R
p ¼ cell pressure, psia
Tsc ¼ standard temperature, °R
psc ¼ standard pressure, psia
Moses (1986)) pointed out that reporting the experimental data in relation to
the residual oil volume at 60°F (as shown graphically in Figures 3-5 and 3-6)
gives the relative oil volume Bod and that the differential gas-oil ratio Rsd
curves the appearance of the oil formation volume factor B o and the solution
gas solubility Rs curves, leading to their misuse in reservoir calculations.
It should be pointed out that the differential liberation test represents the
behavior of the oil in the reservoir as the pressure declines. We must find a
way of bringing this oil to the surface through separators and into the stock
tank. This process is a flash or separator process.
138 Reservoir Engineering
Handbook

FIGURE 3-5 Relative volume versus pressure.

SEPARATOR TESTS
Separator tests are conducted to determine the changes in the volumetric
behav- ior of the reservoir fluid as the fluid passes through the separator (or
separators) and then into the stock tank. The resulting volumetric behavior is
influenced to a large extent by the operating conditions, i.e., pressures and
temperatures, of the surface separation facilities. The primary objective of
conducting separator tests, therefore, is to provide the essential laboratory
information necessary for determining the optimum surface separation
conditions, which in turn will maximize the stock-tank oil production. In
addition, the results of the test, when appropriately combined with the
differential liberation test data, provide a means of obtaining the PVT
parameters (Bo, Rs, and Bt) required for petroleum engineering calculations.
These separator tests are performed only on the orig- inal oil at the bubble
point.
The test, as shown schematically in Figure 3-6A, involves placing a hydro-
carbon sample at its saturation pressure and reservoir temperature in a PVT
cell. The volume of the sample is measured as V sat. The hydrocarbon sample is
then displaced and flashed through a laboratory multistage separator system—
commonly one to three stages. The pressure and temperature of these stages
are set to represent the desired or actual surface separation facilities. The gas
lib- erated from each stage is removed and its specific gravity and volume at
standard conditions are measured. The volume of the remaining oil in the last
stage (repre- senting the stock-tank condition) is measured and recorded as
(Vo)st. These
 Laboratory Analysis of Reservoir Fluids Chapter 3
139

GOR1 GOR2

Pb

(V)Pb,T
VST

T, Pb Tsc

Tsep Psc
(B) Psep

FIGURE 3-6 (A) Solution gas-oil ratio versus pressure. (B) Flash liberation (separator) tests.

experimental, measured data can then be used to determine the oil formation
volume factor and gas solubility at the bubble-point pressure as follows:
Vsat
Bofb ¼ (3-12)
ðVoÞst

. VΣ g sc
Rsfb (3-13)
ðV¼oÞst

where
Bofb bubble-point oil formation volume factor, as measured by flash lib-
¼
eration, bbl of the bubble-point oil/STB
Rsfb ¼ bubble-point solution gas-oil ratio as measured by flash liberation,
scf/STB
(Vg)sc ¼ total volume of gas removed from separators, scf (Figure 3-6B)
140 Reservoir Engineering
Handbook

The above laboratory procedure is repeated at a series of different separator

pressures and at a fixed temperature. It is usually recommended that four of
these tests be used to determine the optimum separator pressure, which is usu-
ally considered the separator pressure that results in minimum oil formation
vol- ume factor. At the same pressure, the stocktank oil gravity will be a
maximum and the total evolved gas, i.e., the separator gas and the stock-tank
gas will be at a minimum.
A typical example of a set of separator tests for a two-stage separation sys-
tem, as reported by Moses (1986)), is shown in Table 3-5. By examining the
laboratory results reported in Table 3-5, it should be noted that the optimum
separator pressure is 100 psia, considered to be the separator pressure that
results in the minimum oil formation volume factor. It is important to notice
that the oil formation volume factor varies from 1.474 bbl/STB to 1.495 bbl/
STB while the gas solubility ranges from 768 scf/STB to 795 scf/STB.
Table 3-5 indicates that the values of the crude oil PVT data are dependent
on the method of surface separation. Table 3-6 presents the results of
performing a separator test on the Big Butte crude oil. The differential
liberation data, as expressed in Table 3-4, show that the solution gas-oil ratio
at the bubble point is 933 scf/STB as compared with the measured value of
646 scf/STB from the

TABLE 3-5 Separator Tests (Permission to publish by the

Society
of Petroleum Engineers of AIME. Copyright SPE-AIME.)

Separator Stock-Tank
Pressure Temperature Oil Gravity
(psig) (°F) GOR Rstb* (°API at FVF, Botb**
60°F)
50 75 737
to 0 75 41 40.5 1.481
778
100 75 676
to 0 75 92 40.7 1.474
768
200 75 602
to 0 75 178 40.4 1.483
780
300 75 549
to 0 75 246 40.1 1.495
795
*GOR in cubic feet of gas at 14.65 psia and 60°Fper barrel of stock-tank oil at
60°F.
**
FVF is barrels of saturated oil at 2.620 psig and 220°F per barrel of stock-
 Laboratory Analysis of Reservoir Fluids Chapter 3
141

TABLE 3-6 Separator Tests Data

Separator Flash Analysis
Stoc
Gas/Oil Gas/Oil k Formatio Specific
Flash Ratio Ratio Tank n Separato Gravity of Oil
Conditio (scf/bbl (scf/STbbl Oil Volume r Flashed Phase
ns psig ) (A) ) (B) at Factor Volume Gas (Air 5 Densit
°F 60°F Bofb(C) Factor(D 1.000) y
(°API )
)
1936 247 0.6484

28 130 593 632 1.066 1.132* 0.7823

0 80 13 13 38.8 1.527 1.010 ** 0.8220

Rsfb ¼
646
*Collected and analyzed in the laboratory by gas chromatography.
**
Insufficient quantity for measurement.
(A) Cubic feet of gas at 14.73 psia and 60°F per barrel of oil at indicated pressure and
temperature.
(B) Cubic feet of gas at 14.73 psia and 60°F per barrel of stock-tank oil at 60°F.
(C) Barrels of saturated oil at 1936 psig and 247°F per barrel of stock-tank oil at 60°F.
(D) Barrels of oil at indicated pressure and temperature per barrel of stock-tank oil at 60°F.

separator test. This significant difference is attributed to the fact that the
proces- ses of obtaining residual oil and stock-tank oil from bubble-point oil
are different. The differential liberation is considered as a multiple series of
flashes at the elevated reservoir temperatures. The separator test is generally a
one-or two- stage flash at low pressure and low temperature. The quantity of
gas released will be different and the quantity of final liquid will be different.
Again, it should be pointed out that oil formation volume factor, as expressed
by Equation 3-12, is defined as “the volume of oil at reservoir pressure and
temperature divided by the
resulting stock-tank oil volume after it passes through the surface separators.”

Adjustment of Differential Liberation Data to Separator

Conditions
To perform material balance calculations, the oil formation volume factor Bo
and gas solubility Rs as a function of the reservoir pressure must be available.
The ideal method of obtaining these data is to place a large crude oil sample in
a PVT cell at its bubble-point pressure and reservoir temperature. At some
pres- sure a few hundred psi below the bubble-point pressure, a small portion
of the oil is removed and flashed at temperatures and pressures equal to those
in the surface separators and stock tank. The liberated gas volume and stock-
tank oil volume are measured to obtain Bo and Rs. This process is repeated at
several progressively lower reservoir pressures until complete curves of B o
and Rs ver- sus pressure have been obtained. This procedure is occasionally
conducted in the laboratory. This experimental methodology was originally
proposed by Dodson (1953) and is called the Dodson Method.
142 Reservoir Engineering
Handbook
Amyx et al. (1960) and Dake (1978) proposed a procedure for constructing
the oil formation volume factor and gas solubility curves by using the
 Laboratory Analysis of Reservoir Fluids Chapter 3
143
differential liberation data (as shown in Table 3-4) in conjunction with the
experimental separator flash data (as shown in Table 3-6) for a given set of
sep- arator conditions. The method is summarized in the following steps:
Step 1. Calculate the differential shrinkage factors at various pressures by
dividing each relative oil volume factors Bod by the relative oil
volume factor at the bubble-point Bodb, or:
Bod
Sod ¼ B (3-14)
odb

where
Bod ¼ differential relative oil volume factor at pressure p, bbl/STB
Bodb ¼ differential relative oil volume factor at the bubble-point pres-
sure pb, psia, bbl/STB
Sod ¼ differential oil shrinkage factor, bbl/bbl of bubble-point oil
The differential oil shrinkage factor has a value of one at the
bubble-point and a value less than one at subsequent pressures below pb.
Step 2. Adjust the relative volume data by multiplying the separator (flash)
formation volume factor at the bubble-point Bofb (as defined by
Equation 3-12) by the differential oil shrinkage factor S od (as defined
by Equation 3-14) at various reservoir pressures. Mathematically,
this relationship is expressed as follows:

Bo ¼ BofbSod (3-15)
where
Bo ¼ oil formation volume factor, bbl/STB
Bofb ¼ bubble-point oil formation volume factor, bbl of the bubble-
point oil/STB (as obtained from the separator test)
Sod ¼ differential oil shrinkage factor, bbl/bbl of bubblepoint oil
Step 3. Calculate the oil formation volume factor at pressures above the
bubble- point pressure by multiplying the relative oil volume data
Vrel, as generated from the constant-composition expansion test, by
Bofb, or:

Bo ¼ ðVrelÞ ðBofbÞ (3-16)

where
Bo ¼ oil formation volume factor above the bubble-point pressure,
bbl/STB
Vrel ¼ relative oil volume, bbl/bbl
Step 4. Adjust the differential gas solubility data Rsd to give the required gas
solubility factor Rs
R¼R
B
—
ðR — R ofb
Þ (3-17)
144 Reservoir Engineering
Handbook s sfb sdb sd Bodb
 Laboratory Analysis of Reservoir Fluids Chapter 3
145

where
Rs ¼ gas solubility, scf/STB
Rsfb ¼ bubble-point solution gas-oil ratio from the separator test,
scf/STB
¼
Rsdb solution gas-oil at the bubble-point pressure as measured by the
differential liberation test, scf/STB
¼
Rsd solution gas-oil ratio at various pressure levels as measured by
the differential liberation test, scf/STB
These adjustments will typically produce lower formation volume
factors and gas solubilities than the differential liberation data.
Step 5. Obtain the two-phase (total) formation volume factor Bt by
multiplying values of the relative oil volume V rel below the bubble-
point pressure by Bofb, or:
Bt ¼ ðBofbÞ ðVrelÞ (3-18)
where
B ¼ two-phase formation volume factor, bbl/STB
Vrel ¼ relative oil volume below the pb, bbl/bbl
Similar values for Bt can be obtained from the differential liberation test by
mul- tiplying the relative total volume Btd (see Table 3-4, Column C) by Bofb,
or
Bt ¼ ðBtdÞðBofbÞ=Bodb (3-19)
It should be pointed out that Equations 3-16 and 3-17 usually produce
values less than one for B o and negative values for R s at low pressures. The
calculated curves of Bo and Rs versus pressures must be manually drawn to B o
¼ 1.0 and Rs ¼ 0 at atmospheric pressure.

Example 3-5
The constant-composition expansion test, differential liberation test, and sepa-
rator test for the Big Butte crude oil system are given in Tables 3-2, 3-4, and
3-6, respectively. Calculate:
○ Oil formation volume factor at 4,000 and 1,100 psi
○ Gas solubility at 1,100 psi
○ The two-phase formation volume factor at 1,300 psi

Solution
Step 1. Determine Bodb, Rsdb, Bofb, and Rsfb from Tables 3-4 and 3-6
Bodb ¼ 1:730 bbl=STB Rsdb ¼ 933 scf=STB

Bofb ¼ 1:527 bbl=STB Rsfb ¼ 646 scf=STB

Step 2. Calculate Bo at 4,000 by applying Equation 3-16

Bo ¼ ð0:9657Þ ð1:57Þ ¼ 1:4746 bbl=STB

146 Reservoir Engineering
Handbook

Step 3. Calculate Bo at 1,100 psi by applying Equations 3-14 and 3-15.

1:563
Sod ¼ ¼ 0:9035
1:730

Bo ¼ ð0:9035Þ ð1:527Þ ¼ 1:379 bbl=STB

Step 4. Calculate the gas solubility at 1,100 psi by using Equation 3-17.
. Σ
1:527
Rs ¼ 646 — ð933 — 622Þ ¼ 371 scf=STB
1:730

Step 5. From the pressure-volume relations (i.e., constant-composition data)

of Table 3-2 the relative volume at 1,300 PSI in 1.2579 bbl/bbl.
Using Equation 3-18, calculate Bt to give:

Bt ¼ ð1:527Þ ð1:2579Þ ¼ 1:921 bbl=STB

Applying Equation 3-19 gives:

Bt ¼ ð2:171Þ ð1:527Þ=1:73 ¼ 1:916 bbl=STB

Table 3-7 presents a complete documentation of the adjusted differential

vapor- ization data for the Big Butte crude oil system. Figures 3-7 and 3-8
compare graphically the adjusted values of R sand Bo with those of the
unadjusted PVT data. It should be noted that no adjustments are needed for
the gas forma- tion volume factor, oil density, or viscosity data.

EXTRAPOLATION OF RESERVOIR FLUID DATA

In partially depleted reservoirs or in fields that originally existed at the
bubble- point pressure, it is difficult to obtain a fluid sample, which usually
represents the original oil in the reservoir at the time of discovery. Also, in
collecting fluid samples from oil wells, the possibility exists of obtaining
samples with a satu- ration pressure that might be lower than or higher than
the actual saturation pressure of the reservoir. In these cases, it is necessary to
correct or adjust the laboratory PVT measured data to reflect the actual
saturation pressure. The proposed correction procedure for adjusting the
following laboratory test data is described in the subsequent sections:
○ Constant-composition expansion (CCE) test
○ Differential expansion (DE) test
○ Oil viscosity test
○ Separator tests
TABLE 3-7 Adjusted Differential Liberation Data
Laboratory Analysis of Reservoir Fluids Chapter 3
147
Differential Vaporization
Adjusted to Separator Conditions*
Solutio Gas
n Formation Formatio Oil Oil/Gas
Pressur Gas/Oil Volume n Densit Viscosit
e Ratio Factor Volume y y Ratio
Psig Rs(A) Bo(B) Factor gm/cc
(C)
6500 646 1.431 0.6919
6000 646 1.439 0.6882
5500 646 1.447 0.6843
5000 646 1.456 0.6803
4500 646 1.465 0.6760
4000 646 1.475 0.6714
3500 646 1.486 0.6665
3000 646 1.497 0.6613
2500 646 1.510 0.6556
2400 646 1.513 0.6544
2300 646 1.516 0.6531
2200 646 1.519 0.6519
2100 646 1.522 0.6506
2000 646 1.525 0.6493
b≫ 1936 646 1.527 0.6484
1700 564 1.482 0.01009 0.6577 19.0
1500 498 1.446 0.01149 0.6650 21.3
1300 434 1.412 0.01334 0.6720 23.8
1100 371 1.379 0.01591 0.6790 26.6
900 309 1.346 0.01965 0.6863 29.8
700 244 1.311 0.02559 0.6944 33.7
500 175 1.271 0.03626 0.7039 38.6
300 95 1.220 0.06075 0.7161 46.0
185 36 1.178 0.09727 0.7256 52.8
120 1.146 0.14562 0.7328 58.4
0 0.7745

*Separator Conditions
Fist Stage 28 psig at 130°F
Stock 0 psig at
Tank 80°F
(A) Cubic feet of gas at 14.73 psia and 60°F per barrel of stock-tank oil at 60°F.
(B) Barrel of oil at indicated pressure and temperature per barrel of stock-tank oil at 60°F.
(C) Cubic feet of gas at indicated pressure and temperature per cubic feet at 14.73 psia and 60°F.
 Laboratory Analysis of Reservoir Fluids Chapter 3
145

FIGURE 3-7 Adjusted gas solubility versus pressure.

FIGURE 3-8 Adjusted oil formation volume factor versus pressure.

Correcting Constant-Composition Expansion Data

The correction procedure, summarized in the following steps, is based on cal-
culating the Y-function value for each point below the “old” saturation
pressure.
Step 1. Calculate the Y-function, as expressed by Equation 3-3, for each point
by using the old saturation pressure.
146 Reservoir Engineering
Handbook

Step 2. Plot the values of the Y-function versus pressure on a Cartesian scale
and draw the best straight line. Points in the neighborhood of the sat-
uration pressure may be erratic and need not be used.
Step 3. Calculate the coefficients a and b of the straight-line equation, i.e.:
Y ¼ a+ bp
Step 4. Recalculate the relative volume V rel values by applying Equation 3-5
and using the “new” saturation pressure, or:
Vrel ¼ 1+ pnew
sat — p (3-20)
pða +
bpÞ

To determine points above the “new” saturation pressure, apply the following
steps:
Step 1. Plot the “old” relative volume values above the “old” saturation
pressure versus pressure on a regular scale and draw the best straight
line through these points.
Step 2. Calculate the slope of the Line S. It should be noted that the slope is
negative, i.e., S < 0.
Step 3. Draw a straight line that passes through the point (Vrel¼ 1, psat new) and
parallel to the line of Step 1.
Step 4. Relative volume data above the new saturation pressure are read from
the straight line or determined from the following expression at any
pressure p:
. Σ
Vrel ¼ 1 — S psatnew — p (3-21)

where
S ¼ slope of the line
p ¼ pressure

Example 3-6
The pressure-volume relations of the Big Butte crude oil system is given in
Table 3-2. The test indicates that the oil has a bubble-point pressure of 1,930
psig at 247°F. The Y-function for the oil system is expressed by the following
linear equation:
Y ¼ 1:0981 + 0:000591p
Above the bubble-point pressure, the relative volume data versus pressure
exhibit a straight-line relationship with a slope of –0.0000138.
The surface production data of the field suggest that the actual bubble-
point pressure is approximately 2,500 psig. Reconstruct the pressure-volume
data using the new reported saturation pressure.
 Laboratory Analysis of Reservoir Fluids Chapter 3
147
Solution
Using Equations 3-30 and 3-31, gives:
Pressure Old New
Psig Vrel Vrel Comments
6500 0.9371 0.944 Equation 3-21
8
6000 0.9422 0.951
7
5000 0.9532 0.965
5
4000 0.9657 0.979
3
3000 0.9805 0.993
1
pb new ¼ 2500 0.9890 1.0000
2000 0.9987 1.1096 Equation 3-20
pb new ¼ 1936 1.0000 1.1299
191 1.0058 1.1384
1
180 1.0324 1.1767
8
160 1.1018 1.1018
0
600 2.4960 2.4960
437 3.4404 3.4404

Correcting Differential Liberation Data

Relative oil volume Bod versus pressure:
The laboratory measured Bod data must be corrected to account for the new
bubble-point pressure pbnew. The proposed procedure is summarized in the fol-
lowing steps:
Step 1. Plot the Bod data versus gauge pressure on a regular scale.
Step 2. Draw the best straight line through the middle pressure range of
30%–90% pb.
Step 3. Extend the straight line to the new bubble-point pressure, as shown
schematically in Figure 3-9.
Step 4. Transfer any curvature at the end of the original curve, i.e., ΔBol at
pbold, to the new bubble-point pressure by placing ΔBol above or
new
below the straight line
b at p .
Step 5. Select any differential pressure Δp below the pbold and transfer the
new
cor- responding curvature to the pressure (pb— Δp).
Step 6. Repeat the above process and draw a curve that connects the
generated Bod points with original curve at the point of intersection
with the straight line. Below this point, no change is needed.

Solution gas-oil ratio:

The correction procedure for the isolation gas-oil ratio R sd data is identical to
that of the relative oil volume data.
148 Reservoir Engineering
Handbook
Correcting Oil Viscosity Data
The oil viscosity data can be extrapolated to a new higher bubble-point
pressure by applying the following steps:
 Laboratory Analysis of Reservoir Fluids Chapter 3
149

Original (Bod)@ old Pb +

New (Bod)@ new Pb ΔB1
ΔB1
Original (Bod)@ old Pb ΔB2
Original Bo curve ΔB1

ΔB2

Bod

Δp
Δp

0.3 Pb 0.9 Pb

(Pb)old (Pb)new

FIGURE 3-9 Adjusting Bod curve to reflect new bubble-point pressure.

Step 1. Defining the fluidity as the reciprocal of the oil viscosity, i.e., 1/μo, cal-
culate the fluidity for each point below the original saturation pressure.
Step 2. Plot fluidity versus pressure on a Cartesian scale (see Figure 3-10).
Step 3. Draw the best straight line through the points and extend it to the new
saturation pressure pbold.
Step 4. New oil viscosity values above pbold are read from the straight line.
To obtain the oil viscosity for pressures above the new bubble-point pressure
pbnew, follow these steps:
Step 1. Plot the viscosity values for all points above the old saturation
pressure on a Cartesian coordinate as shown schematically in Figure
3-11, and draw the best straight line through them, as Line A.
Step 2. Through the point on the extended viscosity curve at p bnew, draw a
straight line (Line B) parallel to A.
Step 3. Viscosities above the new saturation pressure are then read from
Line A.

Correcting the Separator Tests Data

Stock-tank gas-oil ratio and gravity:
No corrections are needed for the stock-tank gas-oil ratio and the stock-tank
API gravity.
150 Reservoir Engineering
Handbook

Extrapolation of crude oil viscosity

.
.
. .
1/o o
.

(ob)old

(ob)New

0 Pressure
(pb)old (pb)new
FIGURE 3-10 Extrapolation of oil viscosity to new Pb.

Extrapolation of crude oil viscosity

µo
A

(µob)old
B
new
(µob)

0 Pressure
(pb)old (pb)new
FIGURE 3-11 Extrapolation of oil viscosity above the new Pb.
 Laboratory Analysis of Reservoir Fluids Chapter 3
151
Separator gas-oil ratio:
The total gas-oil ratio Rsfb is changed in the same proportion as the differential
ratio was changed, or
. Σ
Rnew
sfb
old new
¼ Rsfb R sdb =Rold
sdb
(3-22)

The separator gas-oil ratio is then the difference between the new (cor-
rected) gas solubility Rsfbnew and the unchanged stock-tank gas-oil ratio.

Formation volume factor:

The separator oil formation volume factor B ofb is adjusted in the same propor-
tion as the differential liberation values:
. new old Σ
Bnew
ofb
¼ B old
ofb
Bodb =Bodb (3-23)

Example 3-7
Results of the differential liberation and the separator tests on the Big Butte
crude oil system are given in Tables 3-4 and 3-6, respectively. New field and
production data indicate that the bubble-point pressure is better described by
a value of 2,500 psi as compared with the laboratory reported value of 1,936
psi. The correction procedure for B od and Rsd as described previously was
applied, to give the following values at the new bubble point:

Bnew new
odb ¼ 2:013 bbl=STB R sbd ¼ 1; 134 scf=STB

Using the separator test data as given in Table 3-6, calculate the gas solu-
bility and the oil formation volume factor at the new bubble-point pressure.

Solution
○ Gas solubility: from Equation 3-22
. Σ
1134
Rsb ¼ 646 ¼ 785 scf=STB
933

Separator GOR ¼ 785 — 13 ¼ 772 scf=STB

○ Oil formation volume factor

Applying Equation 3-23, gives
. Σ 2:013
Bob 1:527 1:777 bbl=STB
1:730
¼¼
152 Reservoir Engineering
Handbook

LABORATORY ANALYSIS OF GAS CONDENSATE SYSTEMS

In the laboratory, a standard analysis of a gas-condensate sample consists of:
○ Recombination and analysis of separator samples
○ Measuring the pressure-volume relationship, i.e., constant-composition
expansion test
○ Constant-volume depletion test (CVD)

Recombination of Separator Samples

Obtaining a representative sample of the reservoir fluid is considerably more
dif- ficult for a gas-condensate fluid than for a conventional black-oil
reservoir. The principal reason for this difficulty is that liquid may condense
from the reservoir fluid during the sampling process, and if representative
proportions of both liquid and gas are not recovered then an erroneous
composition will be calculated.
Because of the possibility of erroneous compositions and also because of
the limited volumes obtainable, subsurfacesampling is seldom used in gas-
condensate reservoirs. Instead, surface samplingtechniques are used, and
samples are obtained only after long stabilized flow periods. During this
stabilized flow period, volumes of liquid and gas produced in the surface
separation facilities are accurately mea- sured, and the fluid samples are then
recombined in these proportions.
The hydrocarbon composition of separator samples is also determined by
chromatography or low-temperature fractional distillation or a combination of
both. Table 3-7 shows the hydrocarbon analyses of the separator liquid and
gas samples taken from the Nameless Field. The gas and liquid samples are
recom- bined in the proper ratio to obtain the well stream composition as
given in Table 3-8. The laboratory data indicates that the overall well-stream
system con- tains 63.71 mol% Methane and 10.75 mol% Heptanes-plus.
Frequently, the surface gas is processed to remove and liquefy all
hydrocar- bon components that are heavier than methane, i.e., ethane,
propanes, etc. These liquids are called plant products. These quantities of
liquid products are expressed in gallons of liquid per thousand standard cubic
feet of gas processed, i.e., gal/Mscf, or GPM. McCain (1990) derived the
following expression for cal- culating the anticipated GPM for each
component in the gas phase: . Σ. Σ
psc yi Mi
GPM ¼ 11:173 (3-24)
i Tsc γoi
where
psc ¼ standard pressure, psia
Tsc ¼ standard temperature,
°R
yi ¼ mole fraction of component i in the gas phase
Mi ¼ molecular weight of component i
 Laboratory Analysis of Reservoir Fluids Chapter 3
γoi ¼ specific gravity153
of component i as a liquid at standard conditions
(Chapter 1, Table 1-1, Column E)
154 Reservoir Engineering
Handbook

TABLE 3-8 Hydrocarbon Analyses of Separator Products and

Calculated Wellstream

Separator Separator Gas Well Stream

Component mol % mol % GPM mol % GPM
Hydrogen Sulfide Nil Nil Nil
Carbon Dioxide 0.29 1.17 0.92
Nitrogen 0.13 0.38 0.31
Methane 18.02 81.46 63.71
Ethane 12.08 11.46 11.63
Propane 11.40 3.86 1.083 5.97 1.675
iso-Butane 3.05 0.49 0.163 1.21 0.404
n-Butane 5.83 0.71 0.228 2.14 0.688
iso-Pentane 3.07 0.18 0.067 0.99 0.369
Pentane 2.44 0.12 0.044 0.77 0.284
Hexanes 5.50 0.09 0.037 1.60 0.666
Heptanes-plus 38.19 0.08 0.037 10.75 7.944
100.00 100.00 1.659 100.00 12.03
0
Properties of Heptanes-plus
API gravity @ 60°F 43.4
Specific gravity @
60/60°F 0.8091 0.809
Molecular weight 185 103 185
Calculated separator gas gravity (air ¼ 1.000)
¼ 0.687
Calculated gross heating value for
separator gas ¼ 1209 BTU per cubic foot
of dry gas @ 15.025 psia and 60°F.

Primary separator gas collected@ 745

psigand 74°F.
Primary separator liquid collected @ 745
psig and 74°F.

Primary separator gas/separator liquid 2413 scf/bbl @ 60°F

ratio
Primary separator liquid/stock-tank liquid ratio 1.360 bbl @ 60°F
Primary separator gas/wellstream ratio 720.13 Mscf/MMscf
Stock-tank liquid/wellstream ratio 219.4 bbl/MMscf
 Laboratory Analysis of Reservoir Fluids Chapter 3
155
McCain pointed out that the complete recovery of these products is not
feasible. He proposed that, as a rule of thumb, 5 to 25% of ethane, 80 to 90%
of the pro- pane, 95% or more of the butanes, and 100% of the heavier
components can be recovered from a simple surface facility.
Example 3-8
Table 3-8 shows the wellstream compositional analysis of the Nameless Field.
Using Equation 3-24, calculate the maximum available liquid products assum-
ing 100% plant efficiency.
Solution
○ Using the standard conditions as given in Table 3-8, gives:
. Σ. Σ . Σ
15:025 yi Mi yi Mi
GPM ¼ 11:173 γoi ¼ 0:3228 γoi
520

○ Construct the following working table:

Component yi Mi γoi GPMi
CO2 0.009
2
N2 0.003
1
C1 0.637
1
C2 0.116 30.07 0.3561 1.069
3 0 9
C3 0.059 44.09 0.5069 1.676
7 7 9
i-C4 0.012 58.12 0.5628 0.403
1 3 7
n-C4 0.021 58.12 0.5840 0.688
4 3 1
i-C5 0.009 72.15 0.6311 0.284
9 0 2
n-C5 0.007 72.15 0.6311 0.284
7 0 2
C6 0.016 86.17 0.6638 0.670
0 7 3
C7 +
0.1075 185.00 0.809 7.936
15.20 GPM

Constant-Composition Test
This test involves measuring the pressure-volume relations of the reservoir
fluid at reservoir temperature with a visual cell. This usual PVT cell allows
the visual observation of the condensation process that results from changing
the pres- sures. The experimental test procedure is similar to that conducted
on crude oil systems. The CCE test is designed to provide the dew-point
pressure pd at reservoir temperature and the total relative volume V rel of the
reservoir fluid (relative to the dew-point volume) as a function of pressure.
The relative vol- ume is equal to one at pd. The gas compressibility factor at
pressures greater than or equal to the saturation pressure is also reported. It is
156 Reservoir Engineering
Handbook
only necessary to experimentally measure the z-factor at one pressure p 1 and
determine the gas deviation factor at the other pressure p from:
. Σ
p Vrel (3-25)
z ¼ z1
p1 ðVrelÞ 1
 Laboratory Analysis of Reservoir Fluids Chapter 3
157
where
z ¼ gas deviation factor at p
Vrel ¼ relative volume at pressure p
(Vrel)1 ¼ relative volume at pressure p1
If the gas compressibility factor is measured
. Σ at the dew-point pressure, then:
p
z¼zd ðV Þ (3-26)
p
d rel
where
zd ¼ gas compressibility factor at the dew-point pressure pd
pd ¼ dew-point pressure, psia
p ¼ pressure, psia
Table 3-9 shows the dew-point determination and the pressure-volume
relations of the Nameless Field. The dew-point pressure of the system is
reported as 4,968 psi at 262°F. The measured gas compressibility factor at the
dew point is 1.043.

TABLE 3-9 Pressure-Volume Relations of Reservoir Fluid at

262°F (Constant-Composition Expansion)

Pressure psig Relative Volume Deviation Factor Z

8100 0.8733 1.484
7800 0.8806 1.441
7500 0.8880 1.397
7000 0.9036 1.327
6500 0.9195 1.254
6000 0.9397 1.184
5511 0.9641 1.116
5309 0.9764 1.089
5100 0.9909 1.061
5000 0.9979 1.048
4968 Dew-point 1.0000 1.043
Pressure
4905 1.0057
4800 1.0155
4600 1.0369
4309 1.0725
4000 1.1177
158 Reservoir Engineering
Handbook

TABLE 3-9 Pressure-Volume Relations of Reservoir Fluid at

262°F (Constant-Composition Expansion)—cont’d

Pressure psig Relative Volume Deviation Factor Z

3600 1.1938
3200 1.2970
2830 1.4268
2400 1.6423
2010 1.9312
1600 2.4041
1230 3.1377
1000 3.8780
861 4.5249
770 5.0719

Example 3-9
Using Equation 3-26 and the data in Table 3-9, calculate the gas deviation
factor at 6,000 and 8,100 psi.

Solution
○ At 6,000 psi
. Σ
6000 + 15:025
z 1:043 ¼ 0:9397
ð Þ 1:184
4968 + 15:025
○ At 8,000 psi ¼
. Σ
8; 100 + 15:025
z 1:043 ð Þ 0:8733
¼¼ 1:483
4; 968 + 15:025

Constant-Volume Depletion (CVD) Test

Constant-volume depletion (CVD) experiments are performed on gas conden-
sates and volatile oils to simulate reservoir depletion performance and compo-
sitional variation. The test provides a variety of useful and important
information that is used in reservoir engineering calculations.
The laboratory procedure of the test is shown schematically in Figure 3-12
and is summarized in the following steps:
 Laboratory Analysis of Reservoir Fluids Chapter 3
159

FIGURE 3-12 A schematic illustration of the constant-volume depletion (CVD) test.

Step 1. A measured amount of a representative sample of the original

reservoir fluid with a known overall composition of zi is charged to a
visual PVT cell at the dew-point pressure pd (“a” in Figure 3-12). The
temperature of the PVT cell is maintained at the reservoir
temperature T throughout the experiment. The initial volume Vi of
the saturated fluid is used as a reference volume.
Step 2. The initial gas compressibility factor is calculated from the real gas
equation
pdV i (3-27)
zd ni RT
¼
where
pd ¼ dew-point pressure, psia
Vi ¼ initial gas volume, ft3
ni ¼ initial number of moles of the gas ¼ m/Ma
R ¼ gas constant, 10.73
T ¼ temperature, °R
zd ¼ compressibility factor at dew-point pressure
Step 3. The cell pressure is reduced from the saturation pressure to a
predeter- mined level P. This can be achieved by withdrawing
mercury from the cell, as illustrated in column b of Figure 3-12.
During the process, a second phase (retrograde liquid) is formed. The
fluid in the cell is brought to equilibrium and the gas volume V g and
volume of the ret- rograde liquid VL are visually measured. This
retrograde volume is reported as a percent of the initial volume V i,
which basically repre- sents the
. retrograde
Σ liquid saturation SL:
VL
SL ¼ 100
Vi
160 Reservoir Engineering
Handbook

Step 4. Mercury is reinjected into the PVT cell at constant pressure P while
an equivalent volume of gas is simultaneously removed. When the
initial volume Vi is reached, mercury injection is ceased, as
illustrated in col- umn c of Figure 3-12. This step simulates a
reservoir producing only gas, with retrograde liquid remaining
immobile in the reservoir.
Step 5. The removed gas is charged to analytical equipment where its com-
position yi is determined, and its volume is measured at standard
conditions and recorded as (Vgp)sc. The corresponding moles of gas
produced can be calculated from the expression
psc
. VΣgp sc
np R Tsc (3-28)
¼
where
np ¼ moles of gas produced
(Vgp)sc ¼ volume of gas produced measured at standard conditions,
scf Tsc ¼ standard temperature, °R
psc ¼ standard pressure, psia
R ¼ 10.73
Step 6. The gas compressibility factor at cell pressure and temperature is cal-
culated from the real gas equation-of-state as follows:

p Vg
z . npΣ RT (3-29)
¼
Another property, the two-phase compressibility factor, is also
cal- culated. The two-phase compressibility factor represents the total
compressibility of all the remaining fluid (gas and retrograde liquid)
in the cell and is computed from the real gas law as
pVi
ztwo—phase ¼ . Σ (3-30)
n—i n p RT
where
(ni – np) ¼ the remaining moles of fluid in the cell
ni ¼ initial moles in the cell
np ¼ cumulative moles of gas removed
The two-phase z-factor is a significant property because it is used
when the p/z versus cumulative-gas produced plot is constructed for
evaluating gas-condensate production.
Equation 3-30 can be expressed in a more convenient form by
repla- cing moles of gas, i.e., n and np, with their equivalent gas
volumes, or:
. Σ" #
zd p
ztwo—phase ¼ (3-31)
Pd 1 —
. Σ
Gp =GIIP
 Laboratory Analysis of Reservoir Fluids Chapter 3
161
where
zd ¼ gas deviation factor at the dew-point pressure
Pd ¼ dew-point pressure, psia
P ¼ reservoir pressure, psia
GIIP ¼ initial gas in place, scf
Gp ¼ cumulative gas produced at pressure p, scf
Step 7. The volume of gas produced as a percentage of gas initially in place
is calculated by dividing the cumulative volume of the produced gas
2 both
by the gas initially in place, 3
. atΣ standard conditions
%G ¼
∑ Vgp sc
5100 (3-32)
p
4
GIIP
or

%Gp
" ∑ np#
100
ðniÞoriginal
¼
The above experimental procedure is repeated several times until a minimum
test pressure is reached, after which the quantity and composition of the gas
and retrograde liquid remaining in the cell are determined.
The test procedure can also be conducted on a volatile oil sample. In this
case, the PVT cell initially contains liquid, instead of gas, at its bubble-point
pressure. The results of the pressure-depletion study for the Nameless Field
are illus- trated in Tables 3-10 and 3-11. Note that the composition listed in
the 4,968 psi pressure column in Table 3-10 is the composition of the
reservoir fluid at the dew point and exists in the reservoir in the gaseous
state. Table 3-10 and Figure 3-13 show the changing composition of the
well-stream during deple- tion. Notice the progressive reduction of C7+ below
the dew point and increase
in the Methane fraction, i.e., C1.
The concentrations of intermediates, i.e., C2–C6, are also seen to decrease
(they condense) as pressure drops down to about 2,000 psi, then increase as
they revaporize at the lower pressures. The final column shows the
composition of the liquid remaining in the cell (or reservoir) at the
abandonment pressure of 700 psi; the predominance of C7+ components in
the liquid is apparent.
The z-factor of the equilibrium gas and the two-phase z are presented.
(Note: if a (p/z) versus Gp analysis is to be done, the two-phase
compressibility factors are the appropriate values to use.)
The row in the table, “Wellstream Produced, % of initial GPM from
smooth compositions,” gives the fraction of the total moles (of scf) in the cell
(or res- ervoir) that has been produced. This is total recovery of wellstream
and has not been separated here into surface gas and oil recoveries.
In addition to the composition of the produced wellstream at the final
depletion pressure, the composition of the retrograde liquid is also measured.
The composition of the liquid is reported in the last column of Table 3-10 at
162 Reservoir Engineering
Handbook

TABLE 3-10 Depletion Study at 262°F

Hydrocarbon Analyses of Produced Wellstream-Mol Percent

Reservoir Pressure–
psig
Component 4968 4300 3500 2800 2000 1300 700 700*
Carbon dioxide 0.92 0.97 0.99 1.01 1.02 1.03 1.03 0.30
Nitrogen 0.31 0.34 0.37 0.39 0.39 0.37 0.31 0.02

159
Laboratory Analysis of Reservoir Fluids Chapter 3
Methane 63.71 69.14 71.96 73.24 73.44 72.48 69.74 12.09
Ethane 11.63 11.82 11.87 11.92 12.25 12.67 13.37 5.86
Propane 5.97 5.77 5.59 5.54 5.65 5.98 6.80 5.61
iso-Butane 1.21 1.14 1.07 1.04 1.04 1.13 1.32 1.61
n-Butane 2.14 1.99 1.86 1.79 1.76 1.88 2.24 3.34
iso-Pentane 0.99 0.88 0.79 0.73 0.72 0.77 0.92 2.17
n-Pentane 0.77 0.68 0.59 0.54 0.53 0.56 0.68 1.88
Hexanes 1.60 1.34 1.12 0.98 0.90 0.91 1.07 5.34
Heptanes plus 10.75 5.93 3.79 2.82 2.30 2.22 2.52 61.78
100.00 100.00 100.00 100.00 100.00 100.00 100.00 100.00
Molecular weight of 185 143 133 125 118 114 112 203
heptanes- plus
Specific gravity of heptanes- 0.809 0.777 0.768 0.760 0.753 0.749 0.747 0.819
plus
Continued
 Handbook
160 Reservoir Engineering
TABLE 3-10 Depletion Study at 262°F—cont’d

Hydrocarbon Analyses of Produced Wellstream-Mol Percent

Reservoir Pressure–
psig
Component 4968 4300 3500 2800 2000 1300 700 700
Deviation Factor-Z 1.043 0.927 0.874 0.862 0.879 0.908 0.946
Equilibrium gas
Two-phase 1.043 0.972 0.897 0.845 0.788 0.720 0.603
Wellstream produced— 0.000 7.021 17.957 30.268 46.422 61.745 75.172
Cumulative percent of
initial
GPM from smooth 12.030 7.303 5.623 4.855 4.502 4.624 5.329
compositions Propane-plus
Butanes-plus 10.354 5.683 4.054 3.301 2.916 2.946 3.421
Pentanes-plus 9.263 4.664 3.100 2.378 2.004 1.965 2.261

*Equilibrium liquid phase, representing 13.323 percent of original well stream.

 Laboratory Analysis of Reservoir Fluids Chapter 3
161

TABLE 3-11 Retrograde Condensation During Gas Depletion at 262°F

Pressure Retrograde Liquid Volume

psig Percent of Hydrocarbon Pore Space
4968 Dew-point Pressure 0.0
4905 19.3
4800 25.0
4600 29.9
4300 First Depletion Level 33.1
3500 34.4
2800 34.1
2000 32.5
1300 30.2
700 27.3
0 21.8

700 psi. These data are included as a control composition in the event the
study is used for compositional material-balance purposes.
The volume of the retrograde liquid, i.e., liquid dropout, measured during
the course of the depletion study is shown in Table 3-11. The data are
reshown as a percent of hydrocarbon pore space. The measurements indicate
that the maximum liquid dropout of 34.4% occurs at 3,500 psi. The liquid
dropout can be expressed as a percent of the pore volume, i.e., saturation, by
adjusting the reported values to account for the presence of the initial water
saturation, or
ðSo ¼ ðLDOÞ ð1 — (3-33)
SwiÞ
where
So ¼ retrograde liquid (oil) saturation, %
LDO ¼ liquid dropout, %
Swi ¼ initial water saturation, fraction

Example 3-10
Using the experimental data of the Nameless gas-condensate field given in
Table 3-10, calculate the two-phase compressibility factor at 2,000 psi by
applying Equation 3-31.
Solution
The laboratory report indicates that the base (standard) pressure is 15.025 psia.
Applying Equation 3-31 gives:
162 Reservoir Engineering Handbook

1.0
9
8
7
6 C1

Dew-point pressure
4

C2
1.0
9 C7+
8
7
6
C3
5

4
Mol %

nC4
2

C6

iC4
1.0 iC5
9
CO2
8
nC5
7
6
5

3 N2

0 1000 2000 3000 4000 5000

Pressure, psi

FIGURE 3-13 Hydrocarbon analysis during depletion.

Σ ΣΣ Σ
1:043 2000 + 15:025
Z 2—phase¼ ¼ 0:787
4968 + 1 — 0:46422
15:025

PROBLEMS
Table 3-12 shows the experimental results performed on a crude oil sample
taken from the Mtech field. The results include the CCE, DE, and separator
tests.
 Laboratory Analysis of Reservoir Fluids Chapter 3
163
○ Select the optimum separator conditions and generate Bo, Rs, and Bt values
for the crude oil system. Plot your results and compare with the unadjusted
○ values. Assume that new field indicates that the bubble-point pressure is
better described by a value of 2,500 psi. Adjust the PVT to reflect for the
new bubble-point pressure.

TABLE 3-12 Pressure-Volume Relations of Reservoir Fluid at

260°F (Constant-Composition Expansion)

Pressure psig Relative Volume

5000 0.9460
4500 0.9530
4000 0.9607
3500 0.9691
3000 0.9785
2500 0.9890
2300 0.9938
2200 0.9962
2100 0.9987
2051 1.0000
2047 1.0010
2041 1.0025
2024 1.0069
2002 1.0127
1933 1.0320
1843 1.0602
1742 1.0966
1612 1.1524
1467 1.2299
1297 1.3431
1102 1.5325
862 1.8992
653 2.4711
482 3.4050
Continued
 Handbook
164 Reservoir Engineering
TABLE 3-12 Pressure-Volume Relations of Reservoir Fluid at 260°F (Constant-Composition Expansion)—
cont’d Differential Vaporization at 260°F

Solution Relative Relative Deviation Gas Formation Incremental

Pressure Gas/Oil Oil Total Oil Density Factor Volume Gas
psig Ratio (1) Volume (2) Volume (3) gm/cc Z Factor(4) Gravity
205 1004 1.808 1.808 0.5989
1
190 930 1.764 1.887 0.6063 0.880 0.00937 0.843
0
170 838 1.708 2.017 0.6165 0.884 0.01052 0.840
0
150 757 1.660 2.185 0.6253 0.887 0.01194 0.844
0
130 678 1.612 2.413 0.6348 0.892 0.01384 0.857
0
110 601 1.566 2.743 0.6440 0.899 0.01644 0.876
0
900 529 1.521 3.229 0.6536 0.906 0.02019 0.901
700 456 1.476 4.029 0.6635 0.917 0.02616 0.948
500 379 1.424 5.537 0.6755 0.933 0.03695 0.018
300 291 1.362 9.214 0.6896 0.955 0.06183 1.188
170 223 1.309 16.246 0.7020 0.974 0.10738 1.373
0 0 1.110 0.7298 2.230
at 60°F ¼ 1.000
Gravity of Residual Oil ¼Laboratory Analysis
43.1 °API at 60°F of Reservoir Fluids Chapter 3
165
(1) Cubic feet of gas at 14.73 psia and 60°F per barrel of residual oil at 60°F.
(2) Barrels of oil at indicated pressure and temperature per barrel of residual oil at 60°F.
(3) Barrels of oil plus liberated gas at indicated pressure and temperature per barrel of residual oil at 60°F.
(4) Cubic feet of gas at indicated pressure and temperature per cubic foot at 14.73 psia and 60°F.
TABLE 3-12 Pressure-Volume Relations of Reservoir Fluid at 260°F (Constant-Composition Expansion)—
cont’d Separator Tests of Reservoir Fluid Sample

Separator Separator Stock Tank Formation Separator Specific

Gas/Oil Gas/Oil
Pressure PSI Temperature Gravity °API Volume Volume Gravity of
Ratio (1) Ratio (2)
Gauge °F @ 60°F Factor (3) Factor (4) Flashed Gas
200 to 0 71 431 490 1.138 0.739
*

165
Laboratory Analysis of Reservoir Fluids Chapter 3
71 222 223 48.2 1.549 1.006 1.367
100 to 0 72 522 566 1.083 0.801
*
72 126 127 48.6 1.529 1.006 1.402
50 to 0 71 607 632 1.041 0.869
*
71 54 54 48.6 1.532 1.006 1.398
25 to 0 70 669 682 1.020 0.923
*
70 25 25 48.4 1.558 1.006 1.340

*Collected and analyzed in the laboratory

(1) Gas/oil ratio in cubic feet of gas @ 60°F and 14.75 psi absolute per barrel of oil @ indicated pressure and temperature.
(2) Gas/oil ratio in cubic feet of gas @ 60°F and 14.73 psi absolute per barrel of stock-tank oil @ 60°F.
(3) Formation volume factor in barrels of saturated oil @ 2051 psi gauge and 260°F per barrel of stock-tank oil @ 60°F.
(4) Separator volume factor in barrels of oil @ indicated pressure and temperature per barrel of stock-tank oil @ 60°F.
166 Reservoir Engineering Handbook

REFERENCES
Amyx, J.M., Bass, D.M., Whiting, R., 1960a. Petroleum Reservoir Engineering-Physical
Properties.
McGraw-Hill Book Company, New York.
Dake, L.P., 1978a. Fundamentals of Reservoir Engineering. Elsevier Scientific Publishing Com-
pany, Amsterdam.
Dodson, L.P., 1953. Application of laboratory PVT data to reservoir engineering problems. J. Pet.
Technol. 287–298.
McCain, W., 1990. The Properties of Petroleum Fluids. PennWell Publishing Company, Tulsa,
OK. Moses, P., 1986. Engineering application of phase behavior of crude oil and condensate
systems. J.
Pet. Technol. 38, 715–723.