Gas Reservoir Material Balance & Reserves Estimator (p/Z)

The material balance equation (MBE) formulated by Schilthuis and popularized by Craft, Hawkins, and Terry treats the hydrocarbon reservoir as a giant, isothermal pressure vessel (tank model). For a dry gas reservoir under pressure depletion, the gas formation volume factor (\(B_g\)) converts surface standard volumes (SCF or Sm\(^3\)) into downhole reservoir volumes (RCF or Rm\(^3\)):

\[ B_g = \frac{p_{sc}}{T_{sc}} \cdot \frac{Z \cdot T}{p} = 0.02827 \cdot \frac{Z \cdot T_{res} (\text{R})}{p (\text{psia})} \quad \left[\frac{\text{rcf}}{\text{scf}}\right] \]

In metric field units, \(B_g\) is expressed as \(\text{m}^3/\text{Sm}^3\):

\[ B_g = \frac{T (\text{K}) \cdot Z}{288.15 \cdot p (\text{bar})} \times \frac{1.01325}{1} \]

Derivation of the Linear p/Z Equation

The initial reservoir pore volume occupied by gas is \(V_i = G \cdot B_{gi}\). After cumulative production of \(G_p\) standard volumes of gas, the remaining gas occupying the same pore volume \(V_i\) (assuming no water encroachment and rigid rock) is \((G - G_p)\). Equating pore volumes:

\[ G \cdot B_{gi} = (G - G_p) \cdot B_g \]

Substituting the definition of \(B_g\) and canceling constants yields the famous linear relation:

\[ \frac{p}{Z} = \frac{p_i}{Z_i} \left( 1 - \frac{G_p}{G} \right) \]

Extrapolating this straight line to \(\frac{p}{Z} = 0\) reveals the dynamic volume of gas originally in contact with the wellbores, known as the Connected Original Gas In Place (OGIP). When historical pressure survey points lie precisely along a straight downward trajectory, the reservoir engineer has positive verification that the reservoir is volumetric (closed) with zero significant water encroachment.

A dangerous pitfall in reservoir engineering is applying straight-line \(p/Z\) extrapolation to a reservoir supported by an active water drive. In water-drive reservoirs, encroaching water from an underlying or edge aquifer sustains reservoir pressure. This causes observed \(p/Z\) points to plot above the true depletion line. Extrapolating this shallow trend to zero results in a massive, fictitious overestimation of OGIP (frequently 200% to 300% of actual reserves).

Havlena-Odeh Linearization Framework

Havlena and Odeh (1963) rearranged the generalized material balance equation into straight-line forms depending on the underlying drive mechanism. Defining underground gas expansion \(E_g = B_g - B_{gi}\) and total reservoir fluid withdrawal \(F = G_p B_g + W_p B_w\):

• Closed Volumetric Reservoir: \(F = G \cdot E_g\). A plot of \(F\) versus \(E_g\) yields a straight line passing through the origin \((0,0)\) with slope equal to \(G\).
• Water-Drive Reservoir: \(F = G \cdot E_g + W_e\). Plotting \(\frac{F}{E_g}\) versus \(G_p\) results in an upward-turning curve due to the term \(\frac{W_e}{E_g}\) growing over time as aquifer water encroaches into the gas cap.

The Cole Diagnostic Plot

The Cole plot graphs \(\frac{F}{E_g}\) on the vertical axis against cumulative gas production \(G_p\) on the horizontal axis:

• Horizontal Line: Confirms a closed volumetric gas reservoir with 100% gas expansion drive. The constant y-value represents true OGIP \(G\).
• Upward Convex Curve: Indicates active, strong water drive. As production continues, water encroachment replaces voidage, causing \(\frac{F}{E_g}\) to trend upward.
• Downward Trajectory: Diagnostic of geopressured depletion with high initial formation rock compaction or low-permeability compartmentalization where early wells drain only the immediate completion vicinity before pressure waves reach sealing faults.

While dry gas volumetric material balance provides an elegant linear relationship, real-world field developments encounter geological and fluid complexities that require specialized corrections:

Overpressured / Geopressured Gas Reservoirs

In ultra-deep or overpressured formations (initial pore pressure gradient \(> 0.65\text{ psi/ft}\)), rock matrix pore volume compressibility (\(c_f\)) and connate water expansion (\(c_w\)) provide substantial energy during the early depletion stage. An uncorrected \(p/Z\) plot displays two distinct slopes:

1. An early apparent straight line with a shallow slope reflecting combined gas expansion, rock compaction, and connate water expansion.
2. A distinct downward "knee" inflection once reservoir pressure drops below critical pore collapse thresholds, transitioning to true gas-only expansion.

To linearize geopressured data, the effective compressibility term \(S_{wi} c_w + c_f\) is integrated into the cumulative expansion term per the Ramagost-Farshad or Roach methods.

Retrograde Gas Condensate Reservoirs

If initial reservoir fluid is a gas condensate (rich in \(C_5+\) pentanes-plus fractions), dropping below the dew point pressure triggers liquid condensation within the reservoir pore throats. This "retrograde condensation" immobilizes valuable heavy hydrocarbons as residual oil saturation, permanently altering the flowing gas composition and gas gravity. For retrograde fluids, material balance requires two-phase \(Z\)-factor modeling using compositional EOS flash calculations to account for lean gas shrinkage.

Economic Abandonment Pressure Optimization

Field abandonment pressure (\(p_{ab}\)) is not a physical property of the reservoir; it is an economic threshold determined by surface compression capabilities, pipeline delivery pressures, and operating expenses. Typical abandonment pressures range from 15 to 35 barg (200 to 500 psia). Implementing multi-stage wellhead booster compressors can lower \(p_{ab}\) from 50 barg down to 10 barg, unlocking an additional 10% to 15% recovery factor over the field's productive life.