OilAuthority.com

Bounded Reservoir Depletion: Pseudosteady-State Pressure Decline, Drainage Area Calculation, Material Balance, and Pool Development Implications in WCSB Conventional Pools

A bounded reservoir is a reservoir volume enclosed by no-flow boundaries on all sides — sealing faults, stratigraphic pinch-outs, low-permeability barriers, or the outer limit of the aquifer influence — such that produced fluids cannot be replaced by influx from outside the bounded volume and the reservoir pressure declines continuously as production proceeds, the rate of pressure decline being directly proportional to the production rate and inversely proportional to the pore volume, compressibility, and initial pressure of the trapped fluid system. The pressure behavior during production from a bounded reservoir transitions through three distinct regimes as time progresses: the early-time infinite-acting radial flow regime (where the pressure disturbance has not yet reached any boundary and the well produces as if the reservoir were infinite), the transition regime (where the transient front reaches successive boundaries and the derivative shows doubling or unit-slope increases on the log-log plot), and the pseudosteady-state (PSS) regime (where the transient front has reached all boundaries and the entire drainage volume is depleting uniformly at a constant rate), described by the PSS equation dP/dt = -q × Bo / (Vp × ct) where q is the production rate (m³/day), Bo is the oil formation volume factor, Vp is the pore volume of the bounded drainage area, and ct is the total compressibility. The PSS regime is diagnostically useful: the slope of pressure vs. time (or pressure vs. cumulative production) at constant rate gives a direct measurement of Vp × ct, from which — with known ct from PVT analysis and porosity from logs — the drained pore volume and consequently the drainage area and hydrocarbon pore volume (HCPV) of the bounded system can be computed. For gas reservoirs, the bounded reservoir behavior is revealed on the P/z plot (cumulative gas production Gp vs. P/z where z is the gas compressibility factor): a perfectly bounded, volumetric gas reservoir produces a straight-line P/z decline from initial conditions to abandonment P/z, and the x-intercept of this straight line (where P/z = 0) gives the gas initially in place (GIIP) directly — deviations from linearity indicate aquifer support (concave upward, apparent GIIP overstated by volumetrics) or reservoir compartmentalization (concave downward, early boundary detection). In WCSB conventional pool development, the identification of bounded reservoir behavior is crucial for reserves booking under NI 51-101: producing reserves (1P, 2P, 3P) for a bounded pool depend critically on whether the pool is drained as a single connected volume or as isolated fault-bounded compartments, since a compartmentalized pool with barriers between wells requires infill drilling to drain each compartment individually, while a fully connected bounded pool can be drained with fewer wells at wider spacing. WCSB regulatory practice under AER Directive 040 requires that pool boundary determination — from a combination of pressure transient analysis, fluid contact mapping, and production correlation — be documented in the pool scheme approval process, with boundary conditions explicitly stated as supporting or modifying the drainage area assigned to each well.

Key Takeaways

  • Pseudosteady-state pressure decline rate and its use for pore volume calculation: During PSS, the pressure decline rate at constant production rate is dP/dt = -q × Bo / (Vp × ct). Rearranging: Vp = -q × Bo / (ct × dP/dt). For a WCSB Cardium pool producing at q = 50 m³/day Bo = 1.12 with measured pressure decline of 0.025 MPa/day at constant rate, and ct = 1.5×10^-5 /kPa = 1.5×10^-2 /MPa: Vp = -(50 × 1.12) / (1.5×10^-2 × (-0.025)) = 56 / 3.75×10^-4 = 149,000 m³. If phi = 0.20 and Swi = 0.30 from logs: HCPV = Vp × (1-Swi) = 149,000 × 0.70 = 104,000 m³. At Bo = 1.12: OOIP = 104,000 / 1.12 = 92,860 m³ (stock tank) = 584,000 bbl. This PSS-derived OOIP from production data independently confirms the volumetric calculation and can be used to calibrate the net pay thickness used in the volumetric model when the well drainage area is mapped from production correlation with offset wells.
  • P/z plot linearity as diagnostic for bounded volumetric gas reservoirs: For a bounded, water-free gas reservoir (volumetric depletion), the material balance equation simplifies to: P/z = (Pi/zi) × (1 - Gp/G), where G is GIIP, Pi is initial pressure, zi is initial z-factor, and Gp is cumulative production. This is a straight line on the P/z vs Gp plot with slope -Pi/(zi×G) and x-intercept at G (GIIP). WCSB Foothills tight gas pools in the Cadomin and Nikanassin formations often show nearly perfectly linear P/z trends, confirming bounded volumetric behavior with no aquifer support — this means 100% of the gas in place must be produced by depletion, with recovery factor limited by abandonment pressure (typically 20-30% of initial pressure in Foothills tight gas), giving recovery factors of 70-85% GIIP for well-connected bounded tight gas pools. A concave-up P/z trend (apparent GIIP from x-intercept higher than volumetric GIIP) indicates aquifer influx supporting pressure and reducing net gas production per unit of pressure decline — which in carbonates can appear beneficial early in life but leads to water encroachment and loss of producible gas reserves later.
  • Reservoir compartmentalization within a bounded system: impact on infill spacing and well count: A nominally bounded pool may contain internal no-flow barriers (smaller faults, diagenetic cementation zones, shale baffles) that subdivide the drainage volume into pressure-isolated compartments communicating only through limited cross-fault flow. Compartmentalization is identified by pressure inequality between wells in the same mapped pool: two wells at the same structural elevation with different shut-in pressures (more than 0.5 MPa difference adjusted for GWC depth) are separated by a barrier with transmissibility lower than required for pressure equilibration over the producing life. In WCSB Devonian Leduc reefs, internal fault compartments separated by calcite-sealed fracture zones create reservoir blocks of 0.5-5 million m³ pore volume — each requiring a dedicated producer to drain, versus the unfaulted reef that might be drained by a single central well. Compartmentalization can increase the required well count by 2-5x compared to a simply-connected bounded pool of the same total OOIP, materially affecting the economics of pool development and the AER-approved well spacing for the pool scheme.
  • Aquifer boundary versus no-flow boundary: distinguishing bounded from supported reservoirs through production history: The production decline character of a bounded reservoir (exponential, Arps b = 0 for single-phase liquid production under constant BHP or PSS) differs from an aquifer-supported reservoir (hyperbolic or harmonic decline as aquifer influx partially replaces voidage). A bounded Cardium pool producing under PSS shows a straight line on the reciprocal productivity index (1/PI) vs cumulative oil plot (Craft-Hawkins method), with slope = (1 / (Vp × ct × PI_ref)) confirming bounded volumetric depletion. A Cardium pool with partial aquifer support shows a convex shape on the same plot (slower 1/PI increase per unit cumulative production than pure depletion predicts). Distinguishing the two is critical for WCSB secondary recovery decisions: an aquifer-supported pool may not need waterflood if natural water influx achieves adequate voidage replacement, while a fully bounded pool requires active waterflood injection at voidage replacement ratio (VRR) = 1.0 from early in pool life to maintain reservoir pressure and prevent solution gas exsolution below bubble point.
  • MBH correction for BHSIP in bounded reservoirs: converting shut-in pressure to average reservoir pressure: After a production shut-in, the wellbore pressure recovers toward the average reservoir pressure but does not reach it in finite time because of the bounded geometry — pressure buildup overshoots the drainage-area average pressure as pressure equalizes from the wellbore outward but is bounded by the outer no-flow walls. The Matthews-Brons-Hazebroek (MBH) correction converts the extrapolated Horner straight-line pressure (p*) to true average reservoir pressure (P_bar) using dimensionless correction charts that depend on the shape and size of the drainage area. For a circular bounded reservoir: P_bar = p* - m × MBH(tDAmax), where tDAmax is the dimensionless time at shut-in, m is the Horner slope, and MBH is a tabulated correction. Without the MBH correction, using p* as a proxy for P_bar overestimates average pressure by 0.2-1.0 MPa in WCSB Cardium and Viking bounded pools — leading to optimistic material balance OOIP estimates and underestimation of the degree of reservoir depletion, which in turn affects injection pressure requirements for any planned waterflood.

Material Balance Confirms Bounded Behavior in a WCSB Cardium Pool

A seven-well WCSB Pembina Cardium pool with a mapped pool area of 1,250 ha and average net pay of 6.5 m produces for 8 years from initial pressure of 13.6 MPa. Cumulative oil production: 285,000 m³. At year 8, all-wells average shut-in BHP: 9.8 MPa (after MBH correction applied). P/z not applicable (oil reservoir). Volumetric OOIP: phi × (1-Swi) × net_pay × area / Bo = 0.21 × 0.65 × 6.5 × 1.25×10^7 / 1.13 = 9.84 × 10^6 m³. Pressure-depletion material balance: Np × Bo = Vp × ct × (Pi - P): 285,000 × 1.13 = Vp × 1.4×10^-5 × 1,000,000/kPa × 3,800 kPa; 322,050 = Vp × 0.0532; Vp = 6.05 × 10^6 m³. HCPV from material balance = 6.05 × 10^6 × 0.65 = 3.93 × 10^6 m³. OOIP from MB = 3.93 × 10^6 / 1.13 = 3.48 × 10^6 m³. The material balance OOIP (3.48 × 10^6 m³) is 35% lower than the volumetric OOIP (9.84 × 10^6 m³) — indicating the pool is compartmentalized: only 35% of the mapped pool volume is in pressure communication with the seven producing wells after 8 years of production, consistent with the internal fault pattern mapped on 3D seismic. The undrilled compartments (65% of pool area) represent proved undeveloped reserves requiring 4-5 infill wells to drain, confirmed by seismic attribute analysis identifying two isolated sub-pools separated from the main pool by a calcite-sealed fault corridor.

Fast Facts

The P/z plot method for estimating gas initially in place from the linear pressure-depletion trend of a bounded volumetric gas reservoir was formalized by Cole (1969) and is so reliably diagnostic that WCSB gas producers use it as the primary reserves tool for tight gas pools in the Foothills and Deep Basin where conventional volumetric uncertainty from sparse well control and complex reservoir geometry makes static GIIP estimation unreliable. The P/z straight line can identify aquifer support as early as 5-10% of GIIP produced — when the observed P/z is measurably above the extrapolated depletion line — giving the operator early warning that natural water encroachment will reduce ultimate recovery below the volumetric GIIP and that compression or pressure maintenance may be required to meet contractual gas delivery obligations.

The boundary conditions that define a bounded reservoir — and the pressure transient analysis methods for detecting no-flow boundaries and estimating their distance from a producing well — are described under boundary conditions, where the derivative doubling signature of a no-flow boundary, the constant-pressure boundary derivative downturn, the method of images, and WCSB pool edge identification from well tests are covered. The pseudosteady-state behavior of the bounded reservoir and its role in well test interpretation — including the MBH pressure correction and Horner analysis for average reservoir pressure determination — are described under pressure transient analysis. The material balance methods used to quantify OOIP or GIIP from the pressure-production history of bounded reservoirs — including the P/z method for gas, the Havlena-Odeh plot for oil, and the role of aquifer model selection in distinguishing bounded from supported reservoir behavior — are covered under material balance.

Browse the Oil and Gas Glossary