Gas Reservoir Material Balance: The Complete Guide to p/Z Plots, Drive Mechanism Diagnostics & Reserves Estimation

Executive Summary & Petroleum Reservoir Fundamentals

In petroleum reservoir engineering, accurate quantification of the Original Gas In Place ($G$ or OGIP) and Estimated Ultimate Recovery ($EUR$) is the fundamental technical basis for field development planning, facility sizing, project debt financing, and commercial gas sales contracts.

While 3D geological seismic geomodeling and petrophysical log-core analyses provide volumetric static estimates ($G_{vol} = A \cdot h \cdot \phi \cdot (1 - S_{wi}) / B_{gi}$), static methods suffer from substantial uncertainties in reservoir compartmentalization, fault sealing integrity, and drainage area boundaries.

The Material Balance Equation (MBE)—first formulated by Schilthuis (1936) and perfected by Havlena and Odeh (1963) and Craft and Hawkins (2014)—acts as a dynamic reservoir scale. By tracking reservoir fluid expansion and voidage as static shut-in bottomhole pressure declines with cumulative gas production, material balance determines the connected, drainable dynamic reservoir volume.

For a volumetric gas reservoir producing solely by pressure depletion, plotting $(p/Z)$ versus cumulative gas production ($G_p$) yields a fundamental straight line whose extrapolation directly reveals the Original Gas In Place.

1. First-Principles Derivation of the Gas Material Balance Equation

The reservoir is modeled as a zero-dimensional tank of initial pore volume $V_p$. Applying the law of conservation of mass:

$$\text{Mass of Gas Initially in Place} - \text{Mass of Gas Produced} = \text{Mass of Gas Remaining in Reservoir}$$

At standard surface conditions ($P_{sc} = 14.696\text{ psia} = 1.01325\text{ bar}$, $T_{sc} = 60.0^\circ\text{F} = 15.56^\circ\text{C}$): * Let $G$ = Original Gas In Place ($\text{SCF}$ or $\text{Sm}^3$). * Let $G_p$ = Cumulative gas produced at standard conditions. * Let $(G - G_p)$ = Remaining gas volume in the reservoir at standard conditions.

Gas Formation Volume Factor ($B_g$)

The gas formation volume factor $B_g$ relates the reservoir volume of gas at downhole pressure $p$ and temperature $T$ to its equivalent volume at standard surface conditions:

$$B_g = \frac{V_{res}}{V_{sc}} = \frac{p_{sc}}{T_{sc}} \cdot \frac{Z(p, T) \cdot T}{p}$$

In oilfield units ($p$ in $\text{psia}$, $T$ in ${}^\circ\text{R}$, $B_g$ in $\text{rcf/scf}$): $$B_g = \frac{14.696}{519.67} \cdot \frac{Z \cdot T}{p} = 0.02827 \cdot \frac{Z \cdot T}{p} \quad \left[\frac{\text{rcf}}{\text{scf}}\right]$$

In metric units ($p$ in $\text{kPa}$, $T$ in $\text{K}$, $B_g$ in $\text{m}^3/\text{Sm}^3$): $$B_g = \frac{101.325}{288.15} \cdot \frac{Z \cdot T}{p} = 0.3516 \cdot \frac{Z \cdot T}{p} \quad \left[\frac{\text{m}^3}{\text{Sm}^3}\right]$$

The General Gas Material Balance Formulation

In the general case including water influx ($W_e$), water production ($W_p$), and formation/connate water compressibility ($c_f, c_w$):

$$G (B_g - B_{gi}) + G B_{gi} \left[ \frac{c_w S_{wi} + c_f}{1 - S_{wi}} \right] \Delta p + W_e - B_w W_p = G_p B_g$$

Where: * $B_{gi}$ is the initial gas formation volume factor at initial reservoir pressure $p_i$. * $\Delta p = (p_i - p)$ is the cumulative reservoir pressure drop. * $S_{wi}$ is the initial connate water saturation. * $c_f$ and $c_w$ are the rock matrix and water compressibilities ($\text{psi}^{-1}$). * $W_e$ is cumulative water influx from an adjoining aquifer ($\text{bbl}$ or $\text{m}^3$). * $W_p$ is cumulative produced water.

2. The Classic Volumetric Depletion Straight Line ($p/Z$ vs $G_p$)

In normally pressured, volumetric gas reservoirs where: 1. The reservoir is completely sealed (no active aquifer encroachment: $W_e = 0$, $W_p = 0$). 2. Rock and connate water expansion are negligible compared to gas expansion ($c_f, c_w \ll c_g$).

The general equation simplifies to: $$G (B_g - B_{gi}) = G_p B_g$$

Dividing through by $G \cdot B_g$: $$1 - \frac{B_{gi}}{B_g} = \frac{G_p}{G}$$

Substituting the definition of $B_g = 0.02827 \frac{Z T}{p}$ and recognizing that reservoir temperature $T$ is isothermal: $$\frac{B_{gi}}{B_g} = \frac{(Z_i T / p_i)}{(Z T / p)} = \frac{p / Z}{p_i / Z_i}$$

Substituting into the balance equation yields the classical Craft & Hawkins linear material balance equation:

$$\frac{p}{Z} = \frac{p_i}{Z_i} \left( 1 - \frac{G_p}{G} \right) = \frac{p_i}{Z_i} - m \cdot G_p$$

Where the regression slope $m$ is defined as: $$m = \frac{p_i / Z_i}{G}$$

   (p/Z)
     ^
p_i/Z_i |*
        |  *
        |    *
        |      *  Linear Volumetric Depletion
        |        *
(p/Z)_ab|..........* (Abandonment Point)
        |            *
        |              *
      0 +----------------*----------------> G_p
        0              G_p,ab      G (OGIP)

Extrapolation to Original Gas In Place ($G$)

By measuring static bottomhole pressures over time and calculating $Z(p, T)$ at each point, a linear regression of $(p/Z)$ against $G_p$ allows direct determination of $G$:

$$G = \frac{p_i / Z_i}{|m|}$$

Estimated Ultimate Recovery ($EUR$) & Recovery Factor ($RF$)

At the field's economic abandonment pressure $p_{ab}$ (dictated by separator operating pressure, pipeline delivery pressure, and wellbore liquid loading limits):

$$(p/Z)_{ab} = \frac{p_{ab}}{Z(p_{ab}, T_{res})}$$

$$EUR = G_{p,ab} = G \left[ 1 - \frac{(p/Z)_{ab}}{(p_i / Z_i)} \right]$$

$$RF\% = \frac{EUR}{G} \times 100\% = \left[ 1 - \frac{(p/Z)_{ab}}{(p_i / Z_i)} \right] \times 100\%$$

In typical volumetric gas reservoirs with $p_i = 3,000-5,000\text{ psia}$ and $p_{ab} = 250-400\text{ psia}$, the theoretical recovery factor is exceptionally high: $85\%$ to $92\%$.

3. Reservoir Drive Mechanism Diagnostics & Curve Shapes

Real-world reservoirs frequently deviate from ideal volumetric depletion. Analyzing the curvature and shape of the $p/Z$ plot provides critical early detection of the reservoir drive mechanism:

  (p/Z)
    ^
    |  /-- Water Drive (Concave Upward) -> Traps Gas!
    | /
    |/
    |* 
    |  * \  Geopressured / Overpressured (Apparent 2-Slope)
    |    *  \
    |      *  \
    |        *  \
    |          *  \ Volumetric Straight Line
    +-----------------------------------------> G_p

A. Water Drive / Aquifer Encroachment

When a gas reservoir is in hydraulic communication with an active edge-water or bottom-water aquifer, water encroaches into the reservoir pore space as gas is produced ($W_e > 0$). * Mechanism: The encroaching water replaces voidage and supports reservoir pressure. As a result, the measured pressure $p$ and $(p/Z)$ decline much more slowly than in a closed volumetric tank. * Plot Signature: The $(p/Z)$ vs $G_p$ plot exhibits a characteristic concave upward curvature. * Critical Danger: If an engineer forces a straight-line extrapolation through early production data points of a water-drive reservoir, the extrapolated OGIP ($G$) can be $50\%$ to $200\%$ higher than the true gas in place! * Ultimate Recovery Penalty: Despite higher operating pressures, water drive reservoirs achieve substantially lower recovery factors ($50\%$ to $65\%$ vs $90\%$ for volumetric). Encroaching water invades permeable thief zones, drowning out production wells prematurely, and traps significant residual gas ($S_{gr} = 25-35\%$) at high reservoir pressure under capillary water-block.

B. Geopressured / Overpressured Gas Reservoirs ($p_i > 0.70\text{ psi/ft}$)

In deep, overpressured reservoirs (often exceeding $10,000\text{ psia}$ where pore pressure gradient exceeds normal hydrostatic $0.465\text{ psi/ft}$): * Mechanism: At extreme initial pressure, the rock matrix supports a huge portion of the overburden. As gas is produced and pressure drops, the rock pores compress significantly ($c_f$ is high), and connate water expands ($c_w$). * Plot Signature: The $p/Z$ curve displays an apparent two-slope behavior: an initial steep slope supported by rock and water expansion, followed by a flattening as the rock reaches its terminal compaction. * Correction: The Ramagost-Farshad or Roach modified material balance equation must be applied, plotting $p/Z [1 - (c_w S_{wi} + c_f)/(1 - S_{wi}) \Delta p]$ to restore true linearity.

C. Retrograde Gas Condensate Reservoirs

In rich gas condensate reservoirs, when bottomhole pressure declines below the dew point pressure ($p < p_{dew}$): * Mechanism: Heavy hydrocarbon fractions ($C_5+$) condense out of the vapor phase into liquid within the reservoir pore matrix. Because this retrograde liquid is below its critical saturation, it remains immobile and cannot flow to the wellbore. * Plot Signature: The produced gas becomes progressively leaner, and the effective single-phase $Z$-factor fails. * Correction: A two-phase deviation factor $Z_{2\phi}$ must be computed using composition tracking or constant volume depletion (CVD) laboratory PVT data.

4. Advanced Linearization Methods: Havlena-Odeh & Cole Diagnostics

To mathematically isolate drive mechanisms without relying on visual curve estimation, petroleum engineers employ specialized linear transformation groupings:

A. Havlena-Odeh Linear Formulation

Havlena and Odeh (1963) rearranged the general material balance equation into dynamic reservoir voidage ($F$) and expansion terms:

$$F = G \cdot E_g + W_e$$

Where: * Underground Voidage ($F$): $$F = G_p B_g + W_p B_w \quad [\text{res bbl or rcf}]$$ * Gas Expansion Term ($E_g$): $$E_g = B_g - B_{gi} \quad [\text{rcf/scf}]$$

If the reservoir is purely volumetric ($W_e = 0$), plotting $F$ versus $E_g$ must yield a straight line passing directly through the coordinate origin $(0,0)$ with slope equal to OGIP ($G$).

If the plot curves upward away from the origin, water influx $W_e$ is definitively present.

B. Cole Diagnostic Plot ($F / E_g$ vs $G_p$)

Dividing the Havlena-Odeh equation by $E_g$:

$$\frac{F}{E_g} = G + \frac{W_e}{E_g}$$

• For Volumetric Depletion: $W_e = 0$, so $F / E_g = G = \text{Constant}$. The Cole plot forms a horizontal straight line at value $G$.
• For Water Drive: Water influx causes $W_e / E_g$ to increase over time, resulting in a distinct upward curvature on the Cole plot. This gives reservoir engineers an unmistakable diagnostic fingerprint long before water actually breaks through into production wells.

5. Economic Abandonment Pressure ($p_{ab}$) & Recovery Optimization

The recovery factor ($RF\%$) of a volumetric gas field is ultimately governed by the economic abandonment pressure $p_{ab}$:

$$RF = 1 - \frac{(p/Z)_{ab}}{(p_i / Z_i)}$$

To maximize ultimate recovery, production engineers employ three primary field strategies to lower $p_{ab}$:

1. Wellhead Compression / Depletion Compression Skids:
2. As reservoir pressure depletes, wellhead flowing pressure drops below the minimum required pipeline transmission pressure (typically $60-80\text{ barg} = 870-1,160\text{ psia}$).
3. Installing multi-stage reciprocating or centrifugal field compressors boosts low-pressure wellhead gas ($50-150\text{ psia}$) up to pipeline sales pressure, allowing bottomhole abandonment pressure to drop from $1,000\text{ psia}$ down to $250\text{ psia}$, frequently unlocking an additional $10-15\%$ of total field reserves.
4. Velocity String Retrofits & Artificial Lift:
5. In low-pressure wells, declining gas velocity triggers liquid loading (condensation of water and hydrocarbon liquids dropping back down the tubing).
6. Installing smaller diameter velocity strings or plunger lift restores critical lift velocity ($v_{crit}$ per Turner/Coleman), keeping wells unloaded down to true reservoir abandonment.
7. Surface Flowline De-bottlenecking:
8. Minimizing pressure drop across surface manifolds, separators, and piping lowers the backpressure on the reservoir face.

6. Comparison of Major Gas Reservoir Drive Mechanisms