OG
AshToolbox Oil & Gas Petroleum Engineering Suite
SPE • PRMS Reserve Estimation Compliant

Arps Decline Curve Analysis (DCA) & EUR Forecast

Forecast oil and natural gas production trajectories, determine Estimated Ultimate Recovery (EUR), and calculate remaining recoverable reserves. Support for pure exponential, hyperbolic, harmonic, and modern modified hyperbolic decline with terminal decline rate (D_lim) switching for unconventional shale plays.

Engineering Disclaimer:Calculations are estimates for screening, planning, and academic evaluation. Not a substitute for certified commercial simulators, certified well-control software, API/ASME standards, or qualified engineering judgement.

Arps Decline Curve Analysis & EUR Suite

100% In-Browser • Field Confidential • Offline-Ready

Exponential, Hyperbolic, Harmonic & Modified Hyperbolic (terminal Dlim) forecasting and reserve estimation

Decline Model Parameters

Arps Model ClassificationUnconventional Transient Hyperbolic (b = 1.2 > 1.0, Dlim = 6%)
Cross-Tool Workflow: Pipe Results To...

Seamlessly transfer calculated fluid properties, EUR forecasts, or hydraulics into the next engineering phase (zero server upload).

Estimated Ultimate Recovery (EUR)
1555.9 MSTB
1,555,933 STB
Remaining Reserves
1430.9 MSTB
EUR minus 125.0 MSTB produced
Time to Abandonment
46.8 years
@ 15 BOPD economic cutoff

Production Forecast (Rate & Cumulative Curve)

Year 0 (850 BOPD)• Production Rate q(t)• Econ Limit (15 BOPD)Year 20

EUR Sensitivity Impact (% ΔEUR relative to base)

Impact on Target Output when inputs swing ±20% while holding other parameters constant.

Initial Rate (qi)637.5 → 1,062.5 (Δ 28)
Initial Decline (Di)52 → 78 (Δ 22)
b-Factor0.9 → 1.5 (Δ 34)
Terminal Decline (Dlim)4 → 8 (Δ 16)

Principles of Empirical Production Decline Analysis in Petroleum Reserves

Decline Curve Analysis (DCA) is the most widely utilized empirical method in reservoir engineering and petroleum economics for forecasting well performance, evaluating remaining commercial life, and classifying proved (1P), probable (2P), and possible (3P) reserves under the SPE Petroleum Resources Management System (PRMS).

Introduced by J.J. Arps in 1945, empirical decline analysis models the loss of reservoir drive energy and drainage area depressurization under pseudo-steady-state boundary-dominated flow conditions.

1. The Arps Family of Decline Equations

The fundamental differential definition of fractional loss of production rate per unit time is:

D(t) = - (1 / q) × (dq / dt) = D_i × (q / q_i)^b

Depending on the value of the empirical b-factor, three classical solutions arise:

Exponential Decline (b = 0)
q(t) = q_i × e^(-D_i × t)
N_p(t) = (q_i - q(t)) / D_i

Constant nominal decline. Characteristic of incompressible liquid systems above bubble point, constant bottomhole pressure depletion, and strong waterfloods.

Harmonic Decline (b = 1.0)
q(t) = q_i / (1 + D_i × t)
N_p(t) = (q_i / D_i) × ln(q_i / q(t))

Decline rate is directly proportional to production rate. Occurs in gravity drainage reservoirs and high-permeability gas depletion wells.

Hyperbolic Decline (0 < b < 1)
q(t) = q_i / (1 + b × D_i × t)^(1/b)
N_p(t) = [q_i^b / (D_i(1-b))] × [q_i^(1-b) - q(t)^(1-b)]

Decline rate diminishes continuously over time. Standard for solution-gas drive reservoirs with increasing free gas saturation.

2. The Unconventional Challenge & Modified Hyperbolic (Terminal D_lim)

In ultra-low permeability unconventional tight-oil and shale-gas reservoirs (Permian Basin, Eagle Ford, Bakken, Marcellus), wells experience multi-year transient linear and bilinear flow before boundary-dominated flow ever occurs. When fitted with pure Arps equations, the empirical b-factor frequently exceeds 1.0 (often 1.1 to 1.6).

If an unconstrained hyperbolic curve with b > 1 is extrapolated over a 20-to-30-year well life, the mathematical integral diverges or yields physically impossible infinite reserves.

To rectify this fatal flaw, reservoir engineers utilize the Modified Hyperbolic Model (Robertson 1988). The production trajectory follows the hyperbolic curve until the instantaneous decline rate reaches a specified terminal floor, D_lim (typically 5% to 8% per year). At that exact inflection point, the curve switches to an exponential decline:

Instantaneous Decline: D(t) = D_i / (1 + b × D_i × t)
Switching Condition: D(t*) = D_lim  ⇒  t* = (D_i / D_lim - 1) / (b × D_i)
Post-Switch Rate: q(t) = q(t*) × e^(-D_lim × (t - t*))   for t ≥ t*

Frequently Asked Questions (Answer-First)

Why does the Arps b-factor exceed 1.0 in unconventional shale wells?

Arps derived b ≤ 1.0 assuming boundary-dominated flow. In tight-shale nanodarcy rock, multi-stage hydraulic fractures induce extended transient linear flow that can persist for 5 to 10 years, mathematically resulting in apparent b-values between 1.1 and 1.8.

What is terminal decline (D_lim) and why is it necessary?

Terminal decline is an imposed lower boundary (typically 5% to 8%/yr) where hyperbolic decline switches to exponential decline. Without D_lim, an unconstrained hyperbolic curve with b > 1 will vastly overpredict reserves, predicting unrealistically high production rates decades into the future.

How is the economic limit (q_ab) determined?

The economic limit is the minimum production rate where monthly well revenue equals monthly direct lease operating expenses (LOE). When net oil/gas revenues after royalties and severance taxes fail to cover lifting costs, the well becomes uneconomic and is plugged.