Modified Hyperbolic Decline Model

Overview

The Modified Hyperbolic Decline Model addresses a fundamental limitation of classical hyperbolic decline: when b>0b > 0, the hyperbolic model predicts a continuously decreasing decline rate that can lead to physically unrealistic EUR estimates. The modified hyperbolic model solves this by transitioning from hyperbolic to exponential decline at a specified terminal (limiting) decline rate.

Key Characteristics

Feature Hyperbolic Modified Hyperbolic
Early behavior Hyperbolic Hyperbolic
Late behavior Decreasing D Constant D (exponential)
EUR for b > 1 Unbounded Bounded
Physical realism Limited for long-term Improved
Parameters qiq_i, DiD_i, bb qiq_i, DiD_i, bb, DlimD_{lim}

When to Use

  • Conventional reservoirs with established decline but b-factor approaching or exceeding 1
  • Long-term forecasts where hyperbolic alone would overestimate EUR
  • Tight gas and unconventional wells during transition to boundary-dominated flow
  • Economic limit studies requiring bounded EUR estimates

Historical Context

The modified hyperbolic model emerged from practical engineering experience showing that the decline rate for mature wells typically stabilizes at 5-10% per year, regardless of the early-time b-factor. This observation led to the "terminal decline" or "limiting decline" concept, which has become standard practice in many regulatory filings and reserve estimates.


Theory

The Problem with Standard Hyperbolic Decline

From the classical hyperbolic decline, the instantaneous decline rate is:

D(t)=Di1+bDitD(t) = \frac{D_i}{1 + b D_i t}

For b>0b > 0, as tt \to \infty, the decline rate D(t)0D(t) \to 0. This means:

  • Production rate never quite reaches zero
  • Cumulative production continues to infinity
  • EUR depends entirely on the economic limit rate

When b1b \geq 1, the cumulative production integral does not converge, leading to theoretically infinite EUR. This is physically unrealistic.

The Modified Hyperbolic Solution

The modified hyperbolic model is a piecewise function that switches from hyperbolic to exponential decline when the instantaneous decline rate reaches a specified terminal value DlimD_{lim}:

$$q(t) = \begin \frac{(1 + b D_i t){1/b}} & t < t* \[1em] q_i^ \cdot \exp[-D_(t - t^)] & t \geq t^ \end$$

where tt^* is the transition time and qiexpq_i^{exp} is the rate at the transition point.

Transition Time

The transition from hyperbolic to exponential occurs when the instantaneous decline rate equals the terminal decline rate:

D(t)=Di1+bDit=DlimD(t^*) = \frac{D_i}{1 + b D_i t^*} = D_{lim}

Solving for tt^*:

t=1bDi(DiDlim1)t^* = \frac{1}{b D_i}\left(\frac{D_i}{D_{lim}} - 1\right)

This can be simplified to:

t=DiDlimbDiDlimt^* = \frac{D_i - D_{lim}}{b D_i D_{lim}}

Rate at Transition

The rate at the transition point is:

q(t)=qi(1+bDit)1/bq(t^*) = \frac{q_i}{(1 + b D_i t^*)^{1/b}}

For continuity, this becomes the initial rate for the exponential portion:

qiexp=q(t)=qi(DiDlim)1/bq_i^{exp} = q(t^*) = \frac{q_i}{\left(\frac{D_i}{D_{lim}}\right)^{1/b}}

Terminal Decline Rate from Nominal Annual Decline

In practice, the terminal decline rate is often specified as a percentage per year. For a nominal annual decline of pp (as a decimal):

Dlim=ln(1p)365D_{lim} = -\frac{\ln(1 - p)}{365}

For example:

  • 5% per year: Dlim=0.000141D_{lim} = 0.000141 per day
  • 10% per year: Dlim=0.000289D_{lim} = 0.000289 per day
  • 15% per year: Dlim=0.000445D_{lim} = 0.000445 per day

Equations

Rate-Time Relation

Hyperbolic period (t<tt < t^*):

q(t)=qi(1+bDit)1/bq(t) = \frac{q_i}{(1 + b D_i t)^{1/b}}

Exponential period (ttt \geq t^*):

q(t)=qiexpexp[Dlim(tt)]q(t) = q_i^{exp} \cdot \exp[-D_{lim}(t - t^*)]

where:

qiexp=qi(1+bDit)1/b1exp[Dlimt]q_i^{exp} = \frac{q_i}{(1 + b D_i t^*)^{1/b}} \cdot \frac{1}{\exp[-D_{lim} t^*]}

Note: The second factor accounts for the exponential function being evaluated from time zero.

Cumulative Production

Hyperbolic period (t<tt < t^*):

Np(t)=qi(1b)Di[1(1+bDit)1(1/b)]N_p(t) = \frac{q_i}{(1-b)D_i}\left[1 - (1 + b D_i t)^{1-(1/b)}\right]

Exponential period (ttt \geq t^*):

Np(t)=Np(t)+qiexpDlim[exp(Dlimt)exp(Dlimt)]N_p(t) = N_p(t^*) + \frac{q_i^{exp}}{D_{lim}}\left[\exp(-D_{lim} t^*) - \exp(-D_{lim} t)\right]

where Np(t)N_p(t^*) is the cumulative production at the transition time.

EUR Calculation

To a specified economic limit rate qeconq_{econ}:

If qeconq(t)q_{econ} \geq q(t^*) (limit reached during hyperbolic phase):

EUR=qi(1b)Di[1(qeconqi)1b]EUR = \frac{q_i}{(1-b)D_i}\left[1 - \left(\frac{q_{econ}}{q_i}\right)^{1-b}\right]

If qecon<q(t)q_{econ} < q(t^*) (limit reached during exponential phase):

EUR=Np(t)+q(t)qeconDlimEUR = N_p(t^*) + \frac{q(t^*) - q_{econ}}{D_{lim}}

Time to Economic Limit

If limit reached during hyperbolic phase:

tecon=1bDi[(qiqecon)b1]t_{econ} = \frac{1}{b D_i}\left[\left(\frac{q_i}{q_{econ}}\right)^b - 1\right]

If limit reached during exponential phase:

tecon=t+1Dlimln(q(t)qecon)t_{econ} = t^* + \frac{1}{D_{lim}}\ln\left(\frac{q(t^*)}{q_{econ}}\right)


Applicability & Limitations

Typical Parameter Ranges

Parameter Symbol Typical Range Units
Initial rate qiq_i 10 - 10,000 STB/d or Mscf/d
Initial decline rate DiD_i 0.001 - 0.1 1/day
b-factor bb 0 - 2 dimensionless
Terminal decline rate DlimD_{lim} 0.05 - 0.15 1/year

Terminal Decline Rate Guidelines

Well/Reservoir Type Typical DlimD_{lim} (annual) Notes
Conventional oil 5-8% Well-established fields
Conventional gas 5-10% Depletion drive
Tight gas 8-12% May still be in transient flow
CBM 5-8% Long stabilization period
Shale oil 10-15% Higher due to completion degradation
Shale gas 8-12% Conservative for long-term

When to Apply Terminal Decline

  1. b > 0.3: Always consider using modified hyperbolic
  2. b > 1.0: Must use modified hyperbolic or alternative model
  3. Long forecast periods: > 10 years
  4. Regulatory filings: Many agencies require bounded EUR

Selection of Terminal Decline Rate

The terminal decline rate should be based on:

  1. Analogous mature fields in the same basin
  2. Company/operator experience with similar reservoirs
  3. Regulatory requirements (some jurisdictions specify limits)
  4. Physical reasoning about ultimate drainage

Limitations

  1. Arbitrary transition: The switch point is user-specified, not derived from physics
  2. Discontinuous derivative: The rate-time derivative has a discontinuity at tt^*
  3. Not suitable for transient flow: Assumes boundary-dominated behavior
  4. Single-phase assumption: Does not account for changing fluid properties
  5. Stationary parameters: DiD_i and bb assumed constant before transition

Selection Guidance


References

  1. Arps, J.J. (1945). "Analysis of Decline Curves." Transactions of the AIME, 160: 228-247.

  2. Johnson, R.H. and Bollens, A.L. (1927). "The Loss Ratio Method of Extrapolating Oil Well Decline Curves." Transactions of the AIME, 77: 771. SPE-927771-G.

  3. Currie, S.M. (2010). User Guide for the Rate Time Relations Spreadsheet. Texas A&M University.

  4. Ilk, D., Perego, A.D., Rushing, J.A., and Blasingame, T.A. (2008). "Exponential vs. Hyperbolic Decline in Tight Gas Sands — Understanding the Origin and Implications for Reserve Estimates Using Arps' Decline Curves." SPE 116731, SPE Annual Technical Conference and Exhibition, Denver, Colorado.

  5. Robertson, S. (1988). "Generalized Hyperbolic Equation." SPE 18731, unpublished.

  6. Society of Petroleum Evaluation Engineers (SPEE). (2011). Monograph 3: Guidelines for the Practical Evaluation of Undeveloped Reserves in Resource Plays.

  7. Harmony Enterprise. Online Help: Decline Analysis - Modified Hyperbolic.

An unhandled error has occurred. Reload X