Depleted Reservoir Underground Gas Storage (UGS): Cushion-to-Working Ratio, p/Z Material Balance & API RP 1171 Compliance Guide

Interactive Engineering Calculation Tool: Model reservoir pressure cycles, calculate cushion-to-working gas inventory splits, deliverability curves, and compressor horsepower using the companion Depleted Reservoir Underground Gas Storage (UGS) Modeler.

1. Executive Summary & Storage Asset Physics

Underground Natural Gas Storage (UGS) serves as the primary strategic and commercial shock absorber for regional and international natural gas grids. While salt caverns provide high deliverability for high-frequency peak shaving (10–30 day inventory cycling), depleted gas and condensate reservoirs account for over 80% of global underground storage working capacity, providing the seasonal multi-month inventory required to bridge summer injection baseloads with winter space-heating demand spikes.

Converting a depleted hydrocarbon reservoir into a cyclic storage facility transforms a one-way depletion asset into a dynamic, bi-directional mechanical spring. The reservoir pore volume ($\phi$), caprock capillary threshold, wellbores, and gathering surface network undergo cyclical stress regimes between Maximum Storage Pressure ($p_{max}$)—traditionally bounded at or slightly below original discovery pressure ($p_i$)—and Abandonment / Base Pressure ($p_{min}$).

       SEASONAL INVENTORY CYCLE (DEPLETED RESERVOIR UGS)
Pressure / Gas Inventory
  ^
  |                  Working Gas (Vw)
p_max |-------------/----------------\------------- (Oct/Nov: Peak Inventory)
      |            /                  \
      |           / Injection          \ Withdrawal
      |          /   Period             \   Period
      |         /  (Summer)              \ (Winter)
p_min |--------/                          \-------- (Mar/Apr: Turnover Point)
      |============================================
      |          Cushion / Base Gas (Vc)
      |          (Permanent Reservoir Spring)
  0   +---------------------------------------------> Time (Months)

The reservoir gas inventory is divided into two distinct components: 1. Cushion Gas (Base Gas): The immobile volume of gas required to maintain sufficient reservoir pressure to sustain contract deliverability rates until the final day of the withdrawal cycle. Cushion gas represents 40% to 60% of total reservoir capacity in volumetric depleted fields and constitutes the single largest capital expenditure of a storage project. 2. Working Gas (Top Gas): The cyclic volume injected during off-peak periods (typically April through October in the Northern Hemisphere) and withdrawn during peak market demand periods (November through March).

Managing a depleted reservoir storage asset requires adherence to API RP 1171 (Functional Integrity of Natural Gas Storage in Depleted Hydrocarbon Reservoirs and Aquifer Reservoirs), rigorous monitoring of pressure hysteresis curves, assessment of active bottom-water encroachment, and mitigation of wellbore tubular thermal-mechanical fatigue.

2. Thermodynamics & Reservoir Material Balance Formulations

2.1 The Generalized Gas Storage Material Balance ($p/Z$)

In a closed, volumetric depleted gas reservoir (no water influx, no rock compaction), the relationship between average reservoir pressure ($p_R$) and cumulative inventory is governed by the classical real gas law:

$$\frac{p_R}{Z(p_R, T_R)} = \frac{p_{max}}{Z_{max}} \left[ 1 - \frac{G_p}{G} \right] = \frac{p_{sc} T_R}{Z_{sc} T_{sc} V_p (1 - S_{wi})} \cdot G_R$$

Where: * $p_R$ = Average volumetric reservoir shut-in pressure ($\text{psia}$) * $Z(p_R, T_R)$ = Gas compressibility factor at reservoir conditions (dimensionless, via Hall-Yarborough, Dranchuk-Abu-Kassem, or AGA-8 EOS) * $V_p$ = Total reservoir pore volume ($\text{ft}^3$ or $\text{m}^3$) * $S_{wi}$ = Irreducible connate water saturation (fraction) * $G_R$ = Total standard volume of gas in the reservoir ($\text{MMscf}$) * $G_p$ = Cumulative gas withdrawn ($\text{MMscf}$)

2.2 Working Gas and Cushion Gas Determination

The total gas capacity at maximum inventory ($G_{max}$) and the cushion gas requirement ($G_c$) are defined by:

$$G_{max} = \frac{V_p (1 - S_{wi}) T_{sc}}{p_{sc} T_R} \left( \frac{p_{max}}{Z_{max}} \right)$$

$$G_c = \frac{V_p (1 - S_{wi}) T_{sc}}{p_{sc} T_R} \left( \frac{p_{min}}{Z_{min}} \right)$$

$$V_w = G_{max} - G_c = \frac{V_p (1 - S_{wi}) T_{sc}}{p_{sc} T_R} \left( \frac{p_{max}}{Z_{max}} - \frac{p_{min}}{Z_{min}} \right)$$

The Cushion-to-Working Ratio ($R_{cw}$) is given by:

$$R_{cw} = \frac{G_c}{V_w} = \frac{\frac{p_{min}}{Z_{min}}}{\frac{p_{max}}{Z_{max}} - \frac{p_{min}}{Z_{min}}}$$

2.3 Water Drive Dynamics & Pressure Hysteresis

In water-drive reservoirs, aquifer encroachment occurs during the withdrawal season as $p_R$ declines, invading pore volume ($W_e$). During summer injection, the injected gas must push back the water front. Because gas-water relative permeability exhibits severe capillary hysteresis ($k_{rg}$ during drainage $\ll k_{rg}$ during imbibition), water trapping creates an elliptical pressure-volume hysteresis loop:

$$\frac{p_R}{Z_R} \left( 1 - \frac{W_e - W_p B_w}{V_p} \right) = \frac{p_{max}}{Z_{max}} \left( 1 - \frac{G_{withdrawn}}{G_{max}} \right)$$

Trapped gas behind the imbibing water front can temporarily lock 10% to 25% of the working gas inventory into immobile residual gas saturation ($S_{gr} \approx 20\text{--}35\%$), delaying seasonal inventory recovery.

                  P/Z HYSTERESIS IN WATER-DRIVE STORAGE
       p/Z
        ^             /  Injection Leg (Gas Displacing Water)
        |            / /
        |           / /  Hysteresis
        |          / /   Envelope
        |         / /
        |        / /   Withdrawal Leg (Water Encroachment)
        |       /_/
        +----------------------------------------> Cumulative Gas Inventory (G)

3. Reservoir Geomechanical Integrity & API RP 1171 Compliance

3.1 Maximum Allowable Operating Pressure (MAOP) & Delta-Pressuring

A central engineering question in depleted field storage is Delta-Pressuring—operating the reservoir above its initial discovery pressure ($p_i$) to expand working gas capacity without drilling new wells. Under API RP 1171 Section 6, delta-pressuring requires proving: 1. Caprock Capillary Breakthrough Pressure ($P_{cb}$): Mercury injection capillary pressure (MICP) core testing must confirm: $$p_{max} - p_w < P_{cb} = \frac{2 \gamma_{gw} \cos \theta}{r_{pore, max}}$$ 2. Minimum In-Situ Stress & Hydraulic Fracture Gradient ($g_f$): Step-rate injection tests must establish the fracture breakdown pressure ($p_{frac}$). Regulators enforce a minimum safety margin: $$p_{max} \le 0.80 \times p_{frac} \approx 0.80 \times \left( \nu / (1 - \nu) (\sigma_v - p_p) + p_p \right)$$ Where $\sigma_v$ is total vertical overburden stress ($\approx 1.0\text{--}1.1\ \text{psi/ft}$) and $\nu$ is Poisson's ratio ($0.22\text{--}0.28$ for sandstones).

3.2 Subsurface Barrier Verification (API RP 1171 Workflow)

+--------------------------------------------------------------------------+
|                  API RP 1171 STORAGE INTEGRITY FRAMEWORK                 |
+--------------------------------------------------------------------------+
  |-- Geological Characterization: Caprock thickness > 50 ft, permeability < 10^-6 mD
  |-- Well Integrity Review: Ultrasonic cement logging (USIT/CBL) across storage interval
  |-- Annular Pressure Monitoring: Continuous surface casing & intermediate casing pressure
  |-- Inventory Verification: Biannual shut-in static pressure surveys for p/Z tracking
  |-- Microseismic Monitoring: Real-time array detection during peak injection cycles

4. Well Deliverability & Surface Facilities Architecture

4.1 Inflow Performance Relationship (IPR)

Deliverability during withdrawal is modeled using the Forchheimer equation (accounting for non-Darcy turbulence factor $\beta$) or the empirical Houpeurt / Rawlins-Schellhardt equation:

$$q_{sc} = C \left( p_R^2 - p_{wf}^2 \right)^n$$

Where: * $q_{sc}$ = Gas withdrawal rate ($\text{MMscf/d}$) * $C$ = Performance deliverability coefficient ($\text{MMscf/d}/(\text{psi}^2)^n$) * $n$ = Turbulence exponent ($0.5 \le n \le 1.0$; $n \approx 0.75\text{--}0.85$ typical) * $p_{wf}$ = Flowing bottom-hole pressure ($\text{psia}$)

At the end of the withdrawal season, $p_R$ declines to $p_{min}$, severely compressing the driving differential $(p_R^2 - p_{wf}^2)$. To maintain sales pipeline delivery pressure ($p_{line} \approx 800\text{--}1,200\ \text{psig}$), field booster compression must be staged.

       TYPICAL DEPLETED FIELD UGS SURFACE INFRASTRUCTURE
                     +----------------------------+
                     | Sales Gas Pipeline Network |
                     +--------------+-------------+
                                    |
          +-------------------------+-------------------------+
          | (Withdrawal Mode)                                 | (Injection Mode)
          v                                                   v
+-------------------+                               +-------------------+
| Free-Water Knock  |                               | Bi-Directional    |
| Out & Separators  |                               | Fiscal Metering   |
+---------+---------+                               +---------+---------+
          |                                                   |
          v                                                   v
+-------------------+                               +-------------------+
| TEG Dehydration   |                               | Multi-Stage Gas   |
| Skid (or MolSieve)|                               | Injection Compr.  |
+---------+---------+                               +---------+---------+
          |                                                   |
          v                                                   v
+-------------------+                               +-------------------+
| Withdrawal Gas    |                               | High-Pressure     |
| Compression (Late)|                               | Distribution Man. |
+---------+---------+                               +---------+---------+
          |                                                   |
          +-------------------------+-------------------------+
                                    |
                            [ Storage Wells ]
                                    |
                        +-----------v-----------+
                        | Storage Reservoir     |
                        | (Pore Space Spring)   |
                        +-----------------------+

5. Technology Selection & Field Type Comparison Matrix

6. Step-by-Step Engineering Sizing & Operational Workflow

   STEP 1: RESERVOIR SELECTION & VOLUMETRIC SCREENING
   * Determine discovery pressure (p_i), pore volume (V_p), S_wi, and permeability (k).
   * Confirm no active thief zones, open faults, or un-isolated legacy abandoned wells.
                         |
                         v
   STEP 2: MATERIAL BALANCE & PRESSURE ENVELOPE DEFINITION
   * Set p_max <= 0.80 * p_frac (or <= p_i for non-delta projects).
   * Set p_min based on minimum contract deliverability and compressor suction limits.
   * Compute total gas capacity G_max, cushion gas G_c, and working gas V_w.
                         |
                         v
   STEP 3: WATER INFLUX & HYSTERESIS CHARACTERIZATION
   * Fit Carter-Tracy or van Everdingen-Hurst aquifer model to production history.
   * Calculate trapped gas inventory behind the advancing water front.
                         |
                         v
   STEP 4: WELL DELIVERABILITY & NETWORK NODAL ANALYSIS
   * Run multi-rate deliverability tests on candidate storage wells.
   * Determine required well count N_wells = q_peak_demand / q_well_min.
                         |
                         v
   STEP 5: INJECTION COMPRESSOR & SURFACE SIZING
   * Calculate compression ratio CR = p_max / p_suction_grid.
   * Size multi-stage reciprocating compressors and TEG dehydration contactors.
                         |
                         v
   STEP 6: API RP 1171 INTEGRITY & SURVEILLANCE PLAN
   * Establish baseline logging, annular pressure alarms, and static fall-off schedules.

7. Fully Worked Industrial Example

Problem Statement:

An operator evaluates converting a depleted volumetric sandstone reservoir into a regional peak gas storage facility.

• Reservoir Depth ($D$): $6,500\ \text{ft}$
• Reservoir Temperature ($T_R$): $140^\circ\text{F} = 600\ \text{R}$
• Initial Discovery Pressure ($p_i$): $3,100\ \text{psia}$
• Caprock Fracture Gradient ($g_f$): $0.78\ \text{psi/ft} \implies p_{frac} = 5,070\ \text{psia}$
• Effective Pore Volume ($V_p \cdot (1 - S_{wi})$): $4.85 \times 10^8\ \text{ft}^3$
• Target Maximum Pressure ($p_{max}$): $3,000\ \text{psia}$ (below discovery)
• Target Minimum Pressure ($p_{min}$): $1,200\ \text{psia}$
• Sales Grid Interconnect Pressure: $900\ \text{psig}$ ($915\ \text{psia}$)
• Compressor Pipeline Suction Pressure: $600\ \text{psig}$ ($615\ \text{psia}$)
• Injection Season: 180 days; Peak Withdrawal Period: 90 days

Step 1: Compressibility Factors

Using standard Standing-Katz correlations for $0.60$ specific gravity gas ($T_c = 358\ \text{R}$, $p_c = 668\ \text{psia}$): * At $p_{max} = 3,000\ \text{psia}$: $p_{pr} = 4.49$, $T_{pr} = 1.68 \implies Z_{max} = 0.842$ * At $p_{min} = 1,200\ \text{psia}$: $p_{pr} = 1.80$, $T_{pr} = 1.68 \implies Z_{min} = 0.895$

Step 2: Storage Gas Inventory Calculation

Using standard surface reference conditions ($T_{sc} = 520\ \text{R}$, $p_{sc} = 14.7\ \text{psia}$):

$$G_{max} = \frac{4.85 \times 10^8 \times 520}{14.7 \times 600} \cdot \left( \frac{3,000}{0.842} \right) = 28,594 \times 3,563 = 101.88 \times 10^6\ \text{Mscf} = \mathbf{101.88\ \text{Bscf}}$$

$$G_c = 28,594 \times \left( \frac{1,200}{0.895} \right) = 28,594 \times 1,340.8 = 38.34 \times 10^6\ \text{Mscf} = \mathbf{38.34\ \text{Bscf}}$$

$$V_w = G_{max} - G_c = 101.88 - 38.34 = \mathbf{63.54\ \text{Bscf}}$$

$$\text{Cushion Gas Fraction} = \frac{38.34}{101.88} = \mathbf{37.6\%}; \quad R_{cw} = \frac{38.34}{63.54} = \mathbf{0.603}$$

Step 3: Injection & Withdrawal Rates

• Average Daily Injection Rate ($180\ \text{days}$): $$q_{inj, avg} = \frac{63,540\ \text{MMscf}}{180\ \text{days}} = \mathbf{353\ \text{MMscf/d}}$$
• Average Daily Withdrawal Rate ($90\ \text{days}$): $$q_{wdr, avg} = \frac{63,540\ \text{MMscf}}{90\ \text{days}} = \mathbf{706\ \text{MMscf/d}}$$

Step 4: Compressor Power Sizing (Injection Peak)

• Compressor Suction: $615\ \text{psia}$; Maximum Discharge: $3,200\ \text{psia}$ (accounting for surface friction).
• Overall Pressure Ratio: $r_p = 3,200 / 615 = 5.20$.
• Sizing for 2 stages of reciprocating compression: $r_{p, stg} = \sqrt{5.20} = 2.28$.
• Theoretical Adiabatic Horsepower ($k = 1.30$, $\eta_{mech} = 0.85$): $$\text{BHP} \approx 22 \times \frac{k}{k-1} \times q_{inj} \times (r_{p, stg}^{(k-1)/k} - 1) \times N_{stg} \approx \mathbf{28,400\ \text{HP}}$$

8. Diagnostic Troubleshooting, Operational Hazards & Failure Modes

Engineering Design Verification: To perform detailed real gas $p/Z$ iterations, deliverability nodal matching, and compressor horsepower curves for your field assets, run the Depleted Reservoir Underground Gas Storage Modeler.