Step-by-Step PVT Analysis for Reservoir Engineering and Excel Application

1. Introduction

Pressure–Volume–Temperature (PVT) analysis is one of the most fundamental components of petroleum reservoir engineering. It provides a quantitative description of the phase behavior of reservoir fluids and serves as a bridge between what exists in the reservoir and what is measured at the surface. To accurately perform material-balance calculations, reservoir simulations, well testing, and production forecasting, engineers rely heavily on precise PVT data.

This article introduces the essential concepts of PVT analysis, discusses key laboratory experiments, and explains how field PVT parameters are obtained by combining flash expansion, differential liberation, and separator tests. The article concludes with a full, step-by-step walkthrough of Exercise 2.2 from L. P. Dake’s Fundamentals of Reservoir Engineering textbook.

2. Importance of PVT Analysis in Reservoir Engineering

PVT analysis is essential for:

• Understanding phase changes during reservoir depletion

  As reservoir pressure drops, oil may release dissolved gas, gas may expand, and liquid–gas ratios change. PVT data help engineers predict exactly how fluids behave at each pressure stage.

• Estimating stock-tank oil and produced gas volumes

  Formation volume factors (Bo and Bg) and solution GOR (Rs) allow conversion of reservoir fluid volumes into surface (stock-tank) measurements, enabling accurate production forecasting.

• Evaluating recovery mechanisms

  PVT data reveal whether a reservoir is dominated by solution-gas drive, gas-cap drive, water drive, or a combination, helping engineers understand expected pressure decline and recovery behavior.

• Designing surface separation systems

  Separator pressure and temperature influence shrinkage, GOR, and stock-tank oil yield. PVT flash tests help determine optimum separator conditions for maximizing oil recovery.

• Providing accurate input to reservoir simulators

  Reservoir simulations require pressure-dependent functions such as Bo, Rs, Bg, viscosity, and compressibility. PVT data ensure models accurately reflect real fluid behavior during forecasting.

• Calculating hydrocarbons initially in place (HCIIP)

  Oil and gas in place calculations depend on formation volume factors and fluid densities. Reliable PVT data allow correct estimation of the total hydrocarbons stored in the reservoir.

3. Overview of Key PVT Parameters

• Oil Formation Volume Factor (Bo)

Bo (rb/stb) quantifies the volumetric change that occurs when oil moves from reservoir conditions to stock tank (surface) conditions. Defined as the ratio of the volume of oil at reservoir pressure and temperature to the volume of that same oil after it has been brought to the surface and stabilized at standard conditions (typically 14.7 psia and 60°F), Bo​ is always greater than 1.0 because oil shrinks as dissolved gas is liberated during pressure reduction. For example, a Bo​ value of 1.25 reservoir barrels per stock tank barrel means that for every barrel of oil measured in the stock tank, 1.25 barrels existed in the reservoir, highlighting the significant impact of solution gas on subsurface volume. This parameter is essential for converting surface production data into reservoir voidage volumes for material balance calculations and reserves estimation. Bo typically increases as pressure decreases toward the bubble point due to liquid expansion, then decreases below the bubble point as gas liberation dominates.

• Solution Gas–Oil Ratio (Rs)

Rs represents the amount of gas dissolved in a unit volume of oil at specific reservoir pressure and temperature conditions, expressed in standard cubic feet per stock tank barrel of oil (scf/stb). Rs remains constant above the bubble point since no gas is liberated. As reservoir pressure decreases below the bubble point, gas comes out of solution, causing Rs​ to decline from its maximum value at the bubble point (Rsi) to zero at atmospheric pressure (14.7 psia). This parameter essentially describes the gas content of the oil at any given pressure; for instance, an Rs of 500 scf/stb indicates that each stock tank barrel of oil originally contained 500 standard cubic feet of dissolved gas while underground. Rs is critical for predicting surface gas production rates, designing separation facilities, and modeling the compositional changes in the reservoir during depletion.

• Gas Formation Volume Factor (Bg)

Bg (rb/scf) expresses the dramatic expansion that occurs when gas moves from high-pressure, high-temperature reservoir conditions to standard surface conditions. Defined as the ratio of the volume occupied by a given mass of gas at reservoir conditions to its volume at standard conditions (14.7 psia, 60°F), Bg values are typically very small, often between 0.0001 and 0.01 reservoir barrels per standard cubic foot, reflecting the substantial compressibility of natural gas. This parameter mathematically captures the effects of pressure, temperature, and gas compressibility (Z-factor) through the real gas law. Bg​ enables engineers to relate surface gas production measurements back to the reservoir voidage created by gas production, which is vital for gas material balance, reservoir simulation, and understanding gas cap expansion or gas injection projects. Bg increases continuously as pressure declines due to gas expansion. Bg is not recorded above the bubble point pressure (Pb) because above Pb, the reservoir fluid is considered a single-phase oil (or liquid) system, and gas is not liberated.

4. Laboratory PVT Experiments

Three PVT experiments define how we understand hydrocarbon behavior: the flash expansion, the differential liberation, and the surface separator tests. These experiments capture profoundly different physical processes that occur in both the reservoir and surface facilities. Understanding their differences and their respective limitations is essential for correctly interpreting PVT data and applying it to field operations.

4.1 Flash Expansion Test

The flash expansion test, also known as constant-composition expansion, is a PVT experiment designed primarily to determine the bubble point pressure of a reservoir fluid and measure its total compressibility. In this test, a recombined sample of reservoir oil and gas is placed in a PV cell at reservoir temperature, and the pressure is gradually reduced in stages from an initial value well above the expected bubble point down to lower pressures. Crucially, throughout the entire process, any gas that evolves from the oil is not removed; it remains in continuous contact with the liquid phase inside the closed cell. This means the overall composition of the hydrocarbon system remains unchanged from start to finish. The test measures the total fluid volume relative to the volume at bubble point, producing a characteristic pressure-volume curve that shows a distinct increase in compressibility once gas begins to form. Flash expansion simulates scenarios where gas and oil remain together in equilibrium, such as in a closed reservoir volume before gas mobility begins or within a single separator vessel where phases are temporarily in contact.

Key Characteristics of Flash Expansion:

• No removal of gas: At each pressure step, any liberated gas remains in contact with the oil in the PV cell.

• Constant overall composition: The total mass and composition of hydrocarbons in the cell remain unchanged throughout the experiment.

• Primary purpose: To determine the bubble point pressure (the pressure at which the first gas bubble appears).

• What it measures: Total fluid volume relative to the volume at bubble point.

4.2 Differential Liberation Test

The differential liberation test is a more complex PVT experiment that simulates the depletion process in a producing reservoir where liberated gas migrates away from the oil phase. Starting at bubble point pressure, the pressure is reduced in steps, but after each pressure reduction, all the gas that has just come out of solution is immediately and completely removed from the PV cell before proceeding to the next lower pressure stage. This continuous removal of gas causes the composition of the remaining oil to change progressively, becoming heavier as light components are stripped away. The test measures not only the shrinking oil volume at each stage but also the cumulative volume of gas liberated and its properties (Z-factor, expansion factor). Differential liberation better represents actual reservoir behavior once gas saturation exceeds the critical value and gas begins to flow independently, since in a producing reservoir, gas does not remain in equilibrium with the oil but instead moves toward the wellbore due to mobility differences.

Key Characteristics of Differential Liberation:

• Gas removal after each step: Mercury injection maintains constant cell volume while removing liberated gas.

• Changing composition: As lighter components (methane, ethane) are progressively removed, the remaining oil becomes heavier and more viscous.

• Primary purpose: To simulate reservoir depletion where gas, once liberated, migrates separately from oil due to mobility differences.

• What it measures:

  • Relative oil volume (\(v_o\)): This is the volume of oil measured at the current reservoir pressure and temperature (P and T) after the liberated gas has been removed at each stage of the differential liberation experiment. It is expressed relative to the unit volume of oil at the bubble point pressure ( \(rb/rb_b\)).\(v_o\) is the differential liberation parameter used to derive the Oil Formation Volume Factor (Bo).

  • Relative gas volume at reservoir conditions (\(v_g \) ): This is the incremental volume of gas liberated from the oil during a specific pressure drop (P <Pb​) in the differential liberation experiment, measured at the lower pressure and the fixed reservoir temperature. These gas volumes are expressed relative to the unit volume of bubble point oil.

  • Relative gas volume at standard conditions (\(V_g\)): This is the volume of the incremental liberated gas (which was measured as \(v_g\)) after it has been expanded and measured at standard conditions. Like \(v_g\), this volume is expressed in relative terms.

  • Cumulative relative gas volume (F): This is the total amount of gas liberated from the oil below the bubble point pressure, expressed in relative volumes, but measured at standard conditions. F

    (\(stb/rb_b\)) is defined as the sum of all Vg values (i.e., F=ΣVg​). F is the differential liberation parameter used to derive the Solution Gas-Oil Ratio (Rs)

  • Gas expansion factor (E): This factor is determined by dividing the relative gas volume at standard conditions (\(V_g \) ​) by the relative gas volume at reservoir conditions (\(v_g\)​), calculated as E=Vg​/vg​. E is expressed in units of (scf/rcf). This parameter is used to determine the Gas Formation Volume Factor (Bg​).

  • Gas compressibility factor (Z): This factor is determined by explicitly solving the gas volume equation using the pressure (P), temperature (T), and the Gas Expansion Factor (E) for the liberated gas sample. The relationship used is

$$Z = 35.37\frac{P}{ET}$$

where the values for P and T are the reservoir conditions at which E was measured.

4.3 Surface Separator Tests

Laboratory differential liberation data provides a scientifically rigorous understanding of how oil and gas separate under controlled conditions. However, the journey from reservoir to stock tank involves a crucial intermediate step: surface separation. This process fundamentally alters the final volumes of oil and gas measured at standard conditions, creating what engineers often call "the separator gap", the difference between what differential liberation predicts and what actually reaches the stock tank. Understanding and quantifying this gap is essential for converting laboratory PVT data into parameters that accurately reflect field operations.

When reservoir fluids reach the surface, they don't immediately flash to atmospheric conditions. Instead, they pass through one or more separators operating at specific pressures and temperatures. This staged separation process affects both the quantity and quality of the final products:

1. Volumetric effects: Different separator conditions yield different stock tank oil volumes from the same reservoir fluid.

2. Compositional effects: The distribution of intermediate hydrocarbons (C₃-C₅) between oil and gas phases changes with separator conditions.

3. Economic implications: More stock tank oil means higher revenue, making separator optimization a critical operational decision.

To quantify separator effects, PVT laboratories conduct a series of separator flash tests as part of the overall PVT analysis. These tests provide the essential correction factors needed to adapt differential liberation data to field conditions. The experimental setup involves:

1. Sample preparation: Bubble point oil is placed in the PV cell at reservoir temperature and pressure.

2. System configuration: The cell is connected to a model separator system, either single-stage or multi-stage, with each separator at fixed pressure and temperature.

3. Flash process: The oil is flashed through the separator system to stock tank conditions.

4. Volume measurement: The resulting stock tank oil volume and separator gas volumes are carefully measured.

The essential correction factors needed to adapt differential liberation data to field conditions are:

1. The Shrinkage Factor (Cbf): This factor represents the volume of oil collected in the stock tank relative to the unit volume of oil at the bubble point pressure (\(stb/rb_b\)). This factor is critical because the final volume of oil collected in the stock tank (and thus the field PVT parameters) is dependent on the manner of surface separation.

2. The Initial Solution Gas-Oil Ratio (Rsif): This is the solution gas-oil ratio corresponding to the specific separators used, and it is measured in the experiments in scf/stb.

The results of these separator flash tests (Cbf and Rsif) are then used in conjunction with the fixed differential liberation data to convert the absolute, laboratory PVT parameters (like \(v_o\)​ and F) into the required field parameters (Bo​ and Rs​).

N/B: The ‘f’ subscript indicates that the value comes from the separator flash test.

5. Conversion Equations for Field Parameters

• Oil formation volume factor (Bo)

$$B_o = \frac{v_o}{c_{bf}}$$

• Solution gas-oil ratio (Rs)

$$R_s = R_{sif} - \frac{5.615F}{c_{bf}}$$

• Gas formation volume factor (Bg)

$$B_g = \frac{1}{5.615E}$$

Where:

• \(v_o\) = differential oil volume (\(rb/rb_b\))

• F = cumulative relative gas volume (\(stb/rb_b\))

• E = gas expansion factor (scf/rcf)

• Cbf = flash shrinkage factor (\(stb/rb_b\))

• Rsif = initial solution GOR from separator test (scf/stb)

6. Exercise: Conversion of Differential Data

This exercise is gotten from Exercise 2.2 from L. P. Dake’s Fundamentals of Reservoir Engineering textbook.

Problem Statement

Convert Table 2 (differential liberation data) to field PVT parameters using separator conditions from Table 3.

Table 1: Results of isothermal flash expansion at 200°F

Table 2: Results of isothermal differential liberation at 200º F (All volumes are measured relative to the unit volume of oil at the bubble point pressure of 3330 psi)

Table 3: Separator flash expansion experiments performed on the oil sample whose properties are listed in tables 1 and 2

Solution

Before we delve into the solution, let’s point out some important observations from the differential liberation data in Table 2:

1. The oil volume decreases as pressure drops - oil shrinks as gas comes out of solution.

2. The shrinkage is relatively small at first (only 2.3% by 3000 psia) but accelerates at lower pressures.

3. The oil relative volume at atmospheric pressure (14.7 psia) and reservoir temperature (200°F) is 0.8296 \(rb/rb_b\), showing that oil shrinks by about 17% from bubble point to atmospheric pressure at reservoir temperature.

4. After cooling to 60°F, the oil shrinks further to 0.7794 \(rb/rb_b\), an additional 6% shrinkage due to thermal contraction.

5. This final oil relative volume at atmospheric pressure and standard temperature,\(0.7794 \text{ stb-residual}/rb_b\), represents what we call \(c_{bd}\) - the differential shrinkage factor. This is the residual oil volume after all possible gas has been liberated and the oil has cooled to standard temperature (60°F), relative to the original bubble point oil.

Step 1: Identify Optimum Separator Conditions

From Table 3, among the tested conditions, the separator pressure at 150 psia yields the maximum shrinkage factor (0.7993 \(stb/rb_b\)) and the minimum GOR (510 scf/stb), meaning we recover the most stock tank oil with the least associated gas, generally the economically favorable configuration. Hence, we get the optimum separator conditions as:

• Optimum separator pressure: 150 psia

• Shrinkage factor: \(c_{bf}\) = 0.7993 \(stb/rb_b\)

• Initial solution GOR: \(R_{sif}\) = 510 scf/stb

Step 2: Step-By-Step Conversion

Using the conversion equations in Section 5 of this article, we get the results as shown in Table 4.

Table 4: Field PVT parameters adjusted for single stage, surface separation at 150 psia and 80°F

Case 1: At 4000 psia (Pi)

From Table 1 (Flash Expansion):

• Vt = 0.9925 \(rb/rb_b\)

• Note that above bubble point, \(v_o = V_t\) since all fluid is single-phase oil.

Calculations:

1. Bo:

$$B_o = \frac{V_o}{c_{bf}} = \frac{0.9925}{0.7993} = 1.2417\text{ rb/stb}$$

2. Rs: Since no gas is liberated above bubble point,

$$Rs = Rsi =510 scf/stb$$

2. Bg​: Not applicable above bubble point for oil systems.

Case 2: At Bubble Point (3330 psia)

From Table 2:

• \(v_o\) = 1.0000 (no gas has been liberated)

• F = 0 (no gas liberated yet)

• No E value (no gas liberated)

Calculations:

1. Bo:

$$B_o = \frac{V_o}{c_{bf}} = \frac{1.0000}{0.7993} = 1.2511\text{ rb/stb}$$

2. Rs​: Since no gas is liberated at bubble point,

$$Rs = Rsi =510 scf/stb$$

3. Bg​: Since no E exists at bubble point, we must estimate Bg​ from the gas compressibility equation:

$$B_g = 0.02827\frac{ZT}{5.615P}$$

Using Z = 0.868 (from the first differential step at 3000 psia, Table 2):

$$B_g = 0.02827*\frac{0.868*(200+460)}{5.615*3330} = 0.00087 rb/scf$$

Case 3: At 3000 psia (below bubble point)

From Table 2:

• \(v_o\) = 0.9769 \(rb/rb_b\)

• F = 8.5211 \(stb/rb_b\)

• E = 185.24 rcf/scf

Calculations:

1. Bo:

$$B_o = \frac{V_o}{c_{bf}} = \frac{0.9769}{0.7993} = 1.2222\ \text{rb/stb}$$

2. Rs:

$$R_s = R_{sif} - \frac{5.615 F}{c_{bf}} = 510 - \frac{5.615 \times 8.5211}{0.7993} = 450\ \text{scf/stb}$$

3. Bg:

$$B_g = \frac{1}{5.615 E} = \frac{1}{5.615 \times 185.24} = 0.00096\ \text{rb/scf}$$

7. Differential Data as Reported by Laboratories

Most laboratory reports present Bo and Rs relative to residual oil at 60°F, not bubble-point oil. These are:

• Bod: differential oil FVF relative to residual stock-tank oil

$$B_{od} = \frac{v_o}{c_{bd}}$$

• Rsd: differential solution GOR relative to residual oil

$$R_{sd} = R_{sid}-\frac{5.615F}{c_{bd}}$$

where Rsid is the initial GOR relative to residual oil.

$$R_{sid} = \frac{5.615*\text{Maximum Value of F}}{c_{bd}}$$

Using the differential data in Table 2, the values of Bod and Rsd are calculated as shown in Table 5.

Table 5: Differential PVT parameters as conventionally presented by laboratories, in which Bo and Rs are measured relative to the residual oil volume at 60°F

Case 1: At Initial Reservoir Pressure (4000 psia)

• \(v_o\) = Vt = 0.9925 \(rb/rb_b\) (From the flash expansion data, Table 1)

• Maximum F = 74.955 \(stb/rb_b\) (From differential liberation data, Table 2)

• \(c_{bd} = 0.7794 \text{ stb-residual/}rb_b\)

Calculations:

1. Bod:

$$B_{od} = \frac{v_o}{c_{bd}} = \frac{0.9925}{0.7794} = 1.2734 \text{ rb/stb-residual}$$

2. Rsd: At or above the bubble point,

$$R_{sd} = R_{sid} = \frac{5.615*\text{Maximum Value of F}}{c_{bd}} = \frac{5.615*74.955}{0.7794} = 540 \text{ scf/stb-residual}$$

Case 2: At Bubble Point (3330 psia)

From Table 2:

• \(v_o\) = 1.0000 \(rb/rb_b\)

• F=0 (no gas liberated yet)

Calculations:

1. Bod:

$$B_{od} = \frac{v_o}{c_{bd}} = \frac{1.0000}{0.7794} = 1.2830 \text{ rb/stb-residual}$$

2. Rsd:

$$R_{sd} = R_{sid} = 540 \text{ scf/stb-residual}$$

Case 3: At 3000 psia (below bubble point)

From Table 2:

• \(v_o\) = 0.9769 \(rb/rb_b\)

• F = 8.5211 \(stb/rb_b\)

Calculations:

1. Bod:

$$B_{od} = \frac{v_o}{c_{bd}} = \frac{0.9769}{0.7794} = 1.2534 \text{ rb/stb-residual}$$

2. Rsd:

$$R_{sd} = R_{sid}-\frac{5.615F}{c_{bd}} = 540-\frac{5.615*8.5211}{0.7794} = 479 \text{ scf/stb-residual}$$

However, these values(Bod and Rsd) are not corrected for surface separator conditions and, therefore, cannot be used directly in reservoir calculations unless converted using the following conversion expressions:

Convert Bobd to field Bo:

$$B_o = B_{od} \times \frac{B_{obf}}{B_{obd}}$$

Convert Rsd to field Rs:

$$R_s = R_{sif} - (R_{sid} - R_{sd})\times \frac{B_{obf}}{B_{obd}}$$

Where:

• Rsif = FVF of bubble-point oil from flash separator (stb basis)

• Bobd = FVF of bubble-point oil from differential test (residual-oil basis)

• Rsid = initial dissolved GOR in differential test

As an example, we will correct the Bod and the Rsd values at 2400 psia for surface separator conditions. Given:

• From Table 5: Bod=1.2123 rb/stb-residual, Rsd=378 scf/stb-residual, Bobd=1.2830 rb/stb-residual, Rsid=540 scf/stb-residual

• From Table 3 (at the optimum separator pressure of 150 psia): Bobf=1.2511 \(rb_b/stb\), Rsif=510 scf/stb

Calculate Bo​:

$$B_o = B_{od} \times \frac{B_{obf}}{B_{obd}} = 1.2123 \times \frac{1.2511}{1.2830} = 1.822 \text{ rb/stb}$$

Calculate Rs​:

$$R_s = R_{sif} - (R_{sid} - R_{sd})\times \frac{B_{obf}}{B_{obd}} = 510 - (540-378)\times \frac{1.2511}{1.2830} = 352 \text{ scf/stb}$$

The corrected values, Bo = 1.822 rb/stb and Rs = 352 scf/stb, can be confirmed in Table 4.

8. Application in Microsoft Excel

I will be sharing screenshots of my spreadsheet to provide an overview of how I performed the PVT analysis in Excel.

8.1 Input the Flash Expansion Data

8.2 Input the Differential Liberation Data

8.3 Input the Separator Flash Data

8.4 Calculate the Required PVT Parameters

8.5 Conversion of the Differential PVT Parameters to their Conventional Form and their Correction to Surface Separator Conditions

8.6 Plots of PVT Parameters (Bo, Rs, and Bg) Against Pressure

From the plots, as pressure decreases from the initial reservoir pressure (4000 psia) to atmospheric conditions, the oil formation volume factor Bo (rb/stb) initially increases slightly above the bubble point due to oil expansion, reaching a maximum at the bubble point pressure (3330 psia), then declines as gas comes out of solution. The solution gas-oil ratio Rs (scf/stb) remains constant above the bubble point since no gas is liberated, then decreases linearly with pressure below the bubble point as gas is liberated. Meanwhile, the gas formation volume factor Bg (rb/scf) increases sharply at lower pressures, reflecting gas expansion as pressure drops, consistent with real gas behavior and confirming the expected inverse relationship between pressure and gas volume.

9. Conclusion

PVT analysis remains a foundational element of reservoir engineering, providing the essential link between reservoir fluid behavior and measurable surface volumes. By combining differential liberation data, which accurately represents phase changes within the reservoir, with separator flash data that reflects surface processing conditions, engineers can derive reliable field parameters, such as Bo, Rs, and Bg. Exercise 2.2 from L. P. Dake’s Fundamentals of Reservoir Engineering textbook clearly demonstrates how these conversions transform laboratory measurements into practical engineering values, ensuring accurate material balance calculations, production forecasting, and reservoir simulation. Ultimately, a proper understanding of PVT behavior enables better decision-making, improved reservoir management, and more realistic predictions of hydrocarbon recovery.