1. The Differential Foundation of Production Decline
Arps established that during boundary-dominated pseudo-steady-state flow, the fractional loss of production rate per unit time is related to current flow rate by an empirical exponent b:
Where:
- q_i: Initial stabilized production flow rate (BOPD or MCFD).
- D_i: Initial nominal decline rate (fraction per year or day).
- b: Arps empirical decline exponent (governs rate of change of decline).
2. The Three Classical Decline Models
- Exponential Decline (b = 0): The decline rate remains constant over time.
Rate: q(t) = q_i × e^(-D_i × t)
Cumulative: N_p(t) = (q_i - q(t)) / D_i - Harmonic Decline (b = 1.0): The decline rate is directly proportional to production rate.
Rate: q(t) = q_i / (1 + D_i × t)
Cumulative: N_p(t) = (q_i / D_i) × ln(q_i / q(t)) - Hyperbolic Decline (0 < b < 1): The decline rate slows down continuously as reservoir pressure drops.
Rate: q(t) = q_i / (1 + b × D_i × t)^(1 / b)
Cumulative: N_p(t) = [ q_i^b / (D_i × (1 - b)) ] × [ q_i^(1 - b) - q(t)^(1 - b) ]
Modified Hyperbolic Decline in Unconventional Shale Plays
In multi-fractured horizontal shale wells (Permian, Eagle Ford, Bakken), early transient linear fracture flow produces apparent b-factors greater than 1.0 (often 1.2 to 1.8). If left unconstrained, pure hyperbolic decline projects infinite, mathematically unphysical reserves. Engineers apply Modified Hyperbolic Decline, which switches to a terminal exponential decline (D_lim, typically 6% - 10%/yr) when the instantaneous decline rate drops to that threshold.
A conventional vertical oil well in a water-drive reservoir begins production at q_i = 850 BOPD with an annual nominal exponential decline D_i = 0.22 per year. The economic limit for the well (operating cost equals revenue) is q_econ = 25 BOPD:
- Step 1: Calculate Economic Well Life (t_econ)
q(t) = q_i × e^(-D_i × t)
25 = 850 × e^(-0.22 × t)
ln(25 / 850) = -0.22 × t → -3.5264 = -0.22 × t
t_econ = 16.03 years - Step 2: Calculate Estimated Ultimate Recovery (EUR)
EUR = (q_i - q_econ) / D_i = (850 - 25) / 0.22
EUR = 825 / 0.22 = 3,750 daily-rate units × 365.25 days = 1,369,687 STB (1.37 Million Barrels).
Frequently Asked Questions
What is the physical cause of hyperbolic decline?
In depletion drive (solution gas) reservoirs, as reservoir pressure falls below the bubble point, gas comes out of solution. Gas expands and provides partial energy support to the liquid phase, causing the percentage decline rate to diminish over time.
What is the difference between nominal and effective decline rates?
Nominal decline (D) is continuous instantaneous fractional loss (used in calculus integration). Effective decline (d) is the stepped percentage loss over a discrete calendar period (e.g. 20% annual drop): d = 1 - e^(-D).