Difference between revisions of "Reserves estimation"

From AAPG Wiki
Jump to navigation Jump to search
 
(19 intermediate revisions by 3 users not shown)
Line 55: Line 55:
 
* ø = [[porosity]] (decimal) from log and/or core data
 
* ø = [[porosity]] (decimal) from log and/or core data
 
* ''S''<sub>w</sub> = connate water saturation (decimal) from log and/or core data
 
* ''S''<sub>w</sub> = connate water saturation (decimal) from log and/or core data
* B<sub>oi</sub> = formation volume factor for oil at initial conditions (reservoir bbl/STB) from lab data; a quick estimate is ''B''<sub>oi</sub> = 1.05 + (''N'' × 0.05), where ''N'' is the number of hundreds of ft<sup>3</sup> of gas produced per bbl of oil [for example, in a well with a GOR of 1000, ''B''<sub>oi</sub> = 1.05 + (10 × 0.05)]
+
* B<sub>oi</sub> = formation volume factor for oil at initial conditions (reservoir bbl/STB) from lab data; a quick estimate is <math>B_{oi} = 1.05 + (N \times 0.05)</math>, where ''N'' is the number of hundreds of ft<sup>3</sup> of gas produced per bbl of oil [for example, in a well with a GOR of 1000, ''B''<sub>oi</sub> = 1.05 + (10 × 0.05)]
  
 
Another basic volumetric equation is
 
Another basic volumetric equation is
Line 97: Line 97:
 
|     Edge || 35–60
 
|     Edge || 35–60
 
|-
 
|-
| Gravity || 50–70
+
| [[Gravity]] || 50–70
 
|}
 
|}
  
Line 120: Line 120:
 
One general equation is
 
One general equation is
  
:Change in pore volume = Change in oil volume + change in free gas volume + change in water volume
+
:<math>\text{Change in pore volume} = \text{ Change in oil volume } + \text{ change in free gas volume } + \text{ change in water volume}</math>
  
 
where
 
where
  
* Change in pore volume = ''NB''<sub>oi</sub>/(1 – ''S''<sub>wi</sub>)''c''<sub>f</sub>''P''
+
* <math>\text{Change in pore volume} = \frac{NB_{oi}}{(1 - S_{wi})}c_fP</math>
* Change in oil volume = ''NB''<sub>oi</sub> –(''N'' – ''N''<sub>p</sub>)''B''<sub>oi</sub>
+
* <math>\text{Change in oil volume} = NB_{oi} - (N - N_p)B_{oi}</math>
* Change in gas volume = (''GB''<sub>gi</sub> – ''GB''<sub>g</sub>) + [''N''<sub>p</sub>''R''<sub>p</sub> ''(N'' – ''N''<sub>p</sub>) – NR<sub>si</sub>]''B''<sub>g</sub> due to gas produced, evolved, and encroached from a gas cap
+
* <math>\text{Change in gas volume} = (GB_{gi} - GB_g) + [N_p R_p (N - N_p) - NR_{si}] B_g</math> due to gas produced, evolved, and encroached from a gas cap
* Change in water volume = –''NB''<sub>oi</sub>''S''<sub>wi</sub>/(1 – S<sub>wi</sub>)''c''<sub>w</sub>''P'' – ''W''<sub>e</sub> + ''W''<sub>p</sub>''B''<sub>w</sub>, due to connate water volume change, encroached water, and produced water
+
* <math>\text{Change in water volume} = \frac{-NB_{oi}S_{wi}}{(1 - S_{wi})}c_wP - W_e + W_pB_w</math>, due to connate water volume change, encroached water, and produced water
  
 
where
 
where
Line 175: Line 175:
 
The material balance technique for calculating gas reserves, like material balance for oil, attempts to mathematically equilibrate changes in reservoir volume as a result of production. The basic equation is
 
The material balance technique for calculating gas reserves, like material balance for oil, attempts to mathematically equilibrate changes in reservoir volume as a result of production. The basic equation is
  
:Weight (or SCF) of gas produced = weight (or SCF) of gas initially in the reservoir weight (or SCF) of gas remaining in the reservoir
+
:<math>\text{Weight (or SCF) of gas produced} = \text{ weight (or SCF) of gas initially in the reservoir } - \text{ weight (or SCF) of gas remaining in the reservoir}</math>
  
 
The equations used to calculate OGIP are
 
The equations used to calculate OGIP are
Line 194: Line 194:
 
===Reservoir simulation===
 
===Reservoir simulation===
  
The material balance method is actually a subset of the mathematical techniques that are available to modern petroleum engineers. Reservoir simulators use material balance as well as fluid flow equations to model the reservoir as a group of interconnected tanks. The advent of powerful computers has made the use of numerical simulation quite common for estimating reserves and recovery as well as initial volume in place. Since reservoir simulation can account for performance history through history matching, this method incorporates facets of all the techniques discussed. With sufficient data and prudent use of simulators, this method provides the best recovery estimates for complex reservoirs.
+
The material balance method is actually a subset of the mathematical techniques that are available to modern [[petroleum]] engineers. Reservoir simulators use material balance as well as fluid flow equations to model the reservoir as a group of interconnected tanks. The advent of powerful computers has made the use of numerical simulation quite common for estimating reserves and recovery as well as initial volume in place. Since reservoir simulation can account for performance history through history matching, this method incorporates facets of all the techniques discussed. With sufficient data and prudent use of simulators, this method provides the best recovery estimates for complex reservoirs.
  
 
==Production history analysis==
 
==Production history analysis==
Line 225: Line 225:
 
| Rate of production || <math>q_t = q_{\text{i}}e^{-Dt}</math> || <math> q_{t} = q_{\rm i}(1 + nD_{\rm i}t)^{-1/n}</math>
 
| Rate of production || <math>q_t = q_{\text{i}}e^{-Dt}</math> || <math> q_{t} = q_{\rm i}(1 + nD_{\rm i}t)^{-1/n}</math>
 
|-
 
|-
| Cumulative production || <math>N_p = \frac{(q_{\text{i}} - q_t)}{D}</math> || <math> N_{\rm p} = q_{\rm i}^{n}/[(1 - n)D_{\rm i}](q_{\rm i}^{1 - n} - q_{\rm t}^{1-n})</math>
+
| Cumulative production || <math>N_p = \frac{(q_{\text{i}} - q_t)}{D}</math> || <math> N_{\rm p} = \frac{q_{\rm i}^{n}}{[(1 - n)D_{\rm i}]}(q_{\rm i}^{1 - n} - q_{\rm t}^{1-n})</math>
 
|-
 
|-
| Life of reservoir || <math>t = \left( \frac{1}{D} \right) \ln \left( \frac{q_{\text{i}}}{q_{ec}} \right)</math> || <math>t = [(q_{\rm i}/q_{\rm ec})^{n} - 1/nD_{\rm i}</math>
+
| Life of reservoir || <math>t = \left( \frac{1}{D} \right) \ln \left( \frac{q_{\text{i}}}{q_{ec}} \right)</math> || <math>t = \left( \frac{q_{\rm i}}{q_{\rm ec}} \right)^n -\frac{1}{nD_{\rm i}}</math>
 
|}
 
|}
 +
<sup>where</sup><br>
 +
:<sup>q<sub>t</sub> = Rate of production at time t (BOPD).</sup><br>
 +
:<sup>q<sub>i</sub> = Rate of initial production (BOPD).</sup><br>
 +
<sup>q<sub>ec</sub> = Economic limit rate of production (BOPD).</sup><br>
 +
:<sup>D = Decine rate (decimal).</sup><br>
 +
<sup>D<sub>i</sub> = Initial decline rate (decimal).</sup><br>
 +
:<sup>t= Time (years).</sup><br>
 +
:<sup>n = Exponent usually between 0 and 0.7.</sup><br>
 +
<sup>N<sub>p</sub> = Cumulative production (STBO).
  
 
==Analogy method==
 
==Analogy method==
Line 251: Line 260:
 
* [http://archives.datapages.com/data/alt-browse/aapg-special-volumes/me10.htm Original content in Datapages]
 
* [http://archives.datapages.com/data/alt-browse/aapg-special-volumes/me10.htm Original content in Datapages]
 
* [http://store.aapg.org/detail.aspx?id=612 Find the book in the AAPG Store]
 
* [http://store.aapg.org/detail.aspx?id=612 Find the book in the AAPG Store]
 +
* [http://www.aapg.org/science/discipline/business-and-economics/reserve-estimation#3547338-new-publications Reserves Estimation Articles]
  
 
[[Category:Reservoir engineering methods]]
 
[[Category:Reservoir engineering methods]]
 +
[[Category:Methods in Exploration 10]]

Latest revision as of 17:47, 21 January 2022

Development Geology Reference Manual
Series Methods in Exploration
Part Reservoir engineering methods
Chapter Reserves estimation
Author Rick J. Sustakoski, Diana Morton-Thompson
Link Web page
Store AAPG Store

To better understand reserves estimation, a few important terms require definition. Original oil in place (OOIP) and original gas in place (OGIP) refer to the total volume of hydrocarbon stored in a reservoir prior to production. Reserves or recoverable reserves are the volume of hydrocarbons that can be profitably extracted from a reservoir using existing technology. Resources are reserves plus all other hydrocarbons that may eventually become producible; this includes known oil and gas deposits present that cannot be technologically or economically recovered (OOIP and OGIP) as well as other undiscovered potential reserves.

Estimating hydrocarbon reserves is a complex process that involves integrating geological and engineering data. Depending on the amount and quality of data available, one or more of the following methods may be used to estimate reserves:

  • Volumetric
  • Material balance
  • Production history
  • Analogy

These methods are summarized in Table 1.

Table 1 Summary of methods used to derive hydrocarbon reserves
Method Application Accuracy
Volumetric OOIP, OGIP, recoverable reserves. Use early in life of field. Dependent on quality of reservoir description. Reserves estimates often high because this method does not consider problems of reservoir heterogeneity.
Material balance OOIP, OGIP (assumes adequate production history available), recoverable reserves (assumes OOIP and OGIP known). Use in a mature field with abundant geological, petrophysical, and engineering data. Highly dependent on quality of reservoir description and amount of production data available. Reserve estimates variable.
Production history Recoverable reserves. Use after a moderate amount of production data is available. Dependent on amount of production history available. Reserve estimates tend to be realistic.
Analogy OOIP, OGIP, recoverable reserves. Use early in exploration and initial field development. Highly dependent on similarity of reservoir characteristics. Reserve estimates are often very general.

Volumetric estimation[edit]

Volumetric estimates of OOIP and OGIP are based on a geological model that geometrically describes the volume of hydrocarbons in the reservoir. However, due mainly to gas evolving from the oil as pressure and temperature are decreased, oil at the surface occupies less space than it does in the subsurface. Conversely, gas at the surface occupies more space than it does in the subsurface because of expansion. This necessitates correcting subsurface volumes to standard units of volume measured at surface conditions.

One basic volumetric equation is

where

  • N = OOIP (STB)
  • 7758 = conversion factor from acre-ft to bbl
  • A = area of reservoir (acres) from map data
  • h = height or thickness of pay zone (ft) from log and/or core data
  • ø = porosity (decimal) from log and/or core data
  • Sw = connate water saturation (decimal) from log and/or core data
  • Boi = formation volume factor for oil at initial conditions (reservoir bbl/STB) from lab data; a quick estimate is , where N is the number of hundreds of ft3 of gas produced per bbl of oil [for example, in a well with a GOR of 1000, Boi = 1.05 + (10 × 0.05)]

Another basic volumetric equation is

where

  • G = OGIP(SCF)
  • 43560 = conversion factor from acre-ft to ft3
  • Bgi = formation volume factor for gas at initial conditions (RES ft3/SCF)

Recoverable reserves are a fraction of the OOIP or OGIP and are dependent on the efficiency of the reservoir drive mechanism. The basic equation used to calculate recoverable oil reserves is

where RF = recovery factor, which equals RFP + RFS

The primary recovery factor, RFP, is estimated from the type of drive mechanism (Table 2). The secondary recovery factor, RFS, equals

Table 2 Estimation of primary recovery factor
Drive Mechanidm Primary Recovery Factor Drive Mechanism (%)
Depletion
    Solution gas 18–25
    Expansion 2–5
Gas cap drive 20–40
Water drive
    Bottom 20–40
    Edge 35–60
Gravity 50–70

where

  • ED = displacement efficiency
  • EA = areal sweep efficiency
  • EV = vertical sweep efficiency

These efficiency terms are influenced by such factors as residual oil saturation, relative permeability, reservoir heterogeneity, and operational limitations that govern reservoir production and management. Thus, it is difficult to calculate the recovery factor directly using these terms, and other methods, such as decline curves, are often applied.

The basic equation to calculate recoverable gas reserves is

In this case, the recovery factor (RF) is typically higher than for oil reservoirs; it is often near unity for dry gas reservoirs.

Material balance estimation for oil[edit]

The material balance technique mathematically models the reservoir as a tank. This method uses limiting assumptions and attempts to equilibrate changes in reservoir volume as a result of production. Aquifer support and gas cap expansion can be accounted for by using this method.

One general equation is

where

  • due to gas produced, evolved, and encroached from a gas cap
  • , due to connate water volume change, encroached water, and produced water

where

  • Bg = formation volume factor of free gas
  • Bgi = formation volume factor of free gas at initial conditions
  • cf = formation (rock) compressibility (psi–1)
  • cw = water compressibility (psi–1)
  • N = OOIP (STB)
  • Np = cumulative oil produced (STB); from production history data
  • P = Change in reservoir pressure due to production, that is, initial pressure minus current pressure; taken from field pressure surveys
  • Rp = cumulative gas-oil ratio, or total produced gas (in SCF)/ total produced oil (in STB); from production history data
  • Rsi = inital solution gas-oil ratio (SCF/STB)
  • Swi = initial connate water saturation (decimal)
  • We = cumulative amount of water encroachment; from map and field data
  • Wp = cumulative water produced; from production history data

Another general equation is

where

  • Bt = total (two-phase) formation volume factor
  • Bti = total formation volume factor at initial conditions
  • M = gas cap size expressed as a fraction of initial reservoir oil volume; from map data

This equation assumes thermodynamic equilibrium between oil and gas, a uniform pressure distribution, and a uniform saturation distribution in the reservoir. Additional equations can be derived from the general material balance equation for specific reservoir types.

A simplified equation can be used for a quick estimate of initial oil in place. This equation assumes a closed reservoir system (no active water drive), no initial gas cap, and initial reservoir pressure close to the bubblepoint:

where

  • 5.61 = conversion factor from volume/volume to ft3/bbl
  • Bob = formation volume factor for oil at the bubblepoint; determined for specific separator conditions
  • R = gas-oil ratio, or GOR, equal to produced gas (in SCF)/produced oil (in STB); from production history data
  • Rs = solution gas-oil ratio (SCF/STB) or gas solubility in oil
  • Vi = initial volume of oil plus liberated gas as a function of pressure measured at reservoir temperature
  • Vt = volume of oil plus liberated gas as a function of pressure measured at reservoir temperature; determined under flash liberation conditions

This equation can also be used to predict Np (how much a reservoir can produce, or recoverable reserves) assuming N is determined by an independent method and R, the gas-oil ratio, can be controlled throughout the life of the field.

Material balance estimation for gas[edit]

The material balance technique for calculating gas reserves, like material balance for oil, attempts to mathematically equilibrate changes in reservoir volume as a result of production. The basic equation is

The equations used to calculate OGIP are

Gas reservoir with active water drive:

Gas reservoir with no water drive (We = 0):

where

  • G = OGIP (SCF)
  • Gp = cumulative gas produced (SCF)

These equations can also be used to predict Gp (recoverable reserves) assuming G is determined by an independent method and the production conditions remain constant.

Reservoir simulation[edit]

The material balance method is actually a subset of the mathematical techniques that are available to modern petroleum engineers. Reservoir simulators use material balance as well as fluid flow equations to model the reservoir as a group of interconnected tanks. The advent of powerful computers has made the use of numerical simulation quite common for estimating reserves and recovery as well as initial volume in place. Since reservoir simulation can account for performance history through history matching, this method incorporates facets of all the techniques discussed. With sufficient data and prudent use of simulators, this method provides the best recovery estimates for complex reservoirs.

Production history analysis[edit]

Production history analysis is used to estimate economic ultimate recovery (or recoverable reserves) and the expected economic life of a reservoir. The rate of production and cumulative production at any point in time can also be estimated. This method relies on historical production data to extrapolate future production performance. A variety of curves can be used (Figure 1), the most common being a semilog plot of rate of production versus time (Figure 2). These data are easily obtained through operator records or state regulatory agencies.

Three mathematical models can be used to describe decline curve (usually rate versus time) behavior. They are

  • Exponential decline
  • Hyperbolic decline
  • Harmonic decline

Exponential and hyperbolic decline are commonly used to describe reservoirs. Harmonic decline is an infrequently applied special case of hyperbolic decline.

The different types of decline behavior are not necessarily mutually exclusive. Often different decline curve characteristics are related to different stages of reservoir development, and the overall trends can be significantly affected by workovers or stimulation, infill drilling, a change in lift mechanics, or secondary or tertiary flood initiation (Figure 3).

Formulas used to calculate the rate of production, cumulative production, and economic life of a reservoir are given in Table 3.

Table 3 Decline equations
Solving for Exponential Hyperbolic
Rate of production
Cumulative production
Life of reservoir

where

qt = Rate of production at time t (BOPD).
qi = Rate of initial production (BOPD).

qec = Economic limit rate of production (BOPD).

D = Decine rate (decimal).

Di = Initial decline rate (decimal).

t= Time (years).
n = Exponent usually between 0 and 0.7.

Np = Cumulative production (STBO).

Analogy method[edit]

The analogy method for estimating reserves directly compares a newly discovered or poorly defined reservoir to a known reservoir thought to have similar geological or petrophysical properties (depth, lithology, porosity, and so on). While analogy is the least accurate of the methods presented, it is often used early in the life of a reservoir to establish an order-of-magnitude recovery estimate. As the field matures and data become available to make volumetric OOIP or OGIP estimates, analogy is often used to establish a range of recovery factors to apply to the in-place volumes. Evaluating recovery in this fashion is particularly useful when some performance history is available but a decline rate has yet to be established. Analogy should always be used in conjunction with other techniques to ensure that the results of the more computationally intensive methods make sense within the geological framework.

See also[edit]

References[edit]

  1. 1.0 1.1 IHRDC, 1982, Production rate decline curves: PE107, Boston, MA, IHRDC.
  2. Carr, B. S., Viret, J., 1986, PSD860030: Houston, TX, Chevron Exploration and Production Services, March.

External links[edit]

find literature about
Reserves estimation
Datapages button.png GeoScienceWorld button.png OnePetro button.png Google button.png