Module IV - Dca
Module IV - Dca
Analysis
An Overview
• One of the most extensively used technique of
data analysis for evaluating oil and gas reserves
and for prediction of future production
• Based on assumption that past production trend
with it’s controlling factors will continue in the
Decline future and therefore can be extrapolated
• The method of extrapolation of a ‘trend’ for
Curve predicting future production must satisfy the
condition that the factors which caused changes
• The instantaneous decline rate ‘D’ is defined as the rate of change of natural logarithm of
the production rate i.e. ln(q) with respect to time
• Parameters calculated using classical fit of the historical data, namely instantaneous
decline rate D and the exponent b, can be used to predict future production
• Based on the decline behaviour, the value of b ranges from 0 to 1, and accordingly Arps’
equation can be conveniently expressed in following three forms
Decline Curve Analysis
Decline Curve Analysis
• Should be noted that above set equations is only applicable when the
reservoir is flowing under pseudo steady state conditions
• The area under the decline curve q vs time between times t1 and t2 is a measure of
cumulative gas (or oil) production during this period expressed mathematically as:
• Replacing qt in above expression with the three individual expressions that describe the
decline curves and integrating gives:
Decline Curve Analysis
Decline Curve Analysis
• Care must be taken while using the expressions for cumulative gas (or oil) and flow rate
such that units for rates qi, qt must match with the time unit of decline rate
• For eg: if the flow rate is measured in scf/month (or stb/month), the decline rate D must have unit of
1/month
• Note that traditional Arps’ equation gives reasonable estimate of reserves but completely
ignores the flowing pressure data
• As a result, it can overestimate or underestimate the reserves
Exponential Decline (b = 0)
Or linearly as
Similarly,
Or linearly as
Exponential Decline
• Graphical representation of such plot indicates that the plot of qt vs t on a semilog scale or
a plot of qt vs Gp(t) on a Cartesian scale will result in a straight line
Exponential Decline
• Steps to be followed:
1. Plot qt vs Gp on a cartesian plot and qt vs time on semi log scale
2. For both plots, draw best straight line possible
3. Extrapolate st line qt vs Gp to Gp = 0, which intercepts the Y axis with a flow rate value
that is identified as qi
4. Calculate the initial decline rate Di by selecting a point on the Cartesian st line with
coordinates (qt, Gpt) or on a semilog st line with coordinates (qt, t) and solve for Di
OR
5. Calculate time to reach economic flow rate qa and the corresponding cum gas (or oil)
production
Harmonic Decline (b = 1)
• The above two expressions can be rearranged and expressed respectively as:
Harmonic Decline (b = 1)
• The above two expressions suggest that the plot of 1/qt vs t on a cartesian scale will yield a
st line with slope (Di/qi) and intercept of 1/qt
• Plot of qt vs Gp(t) on a semilog scale will yield a st line with negative slope of (Di/qi) and
intercept of qi
Harmonic Decline (b = 1)
• Time to reach economic flow rate and the corresponding cumulative production is calculated as:
Hyperbolic Decline (0<b<1)
• Two governing relations for a reservoir or well when its production follows hyperbolic
decline behavior are given by equations:
Hyperbolic Decline (0<b<1)
• Following simplified iterative method to be used to determine Di and b from historical
production data:
1. Plot qt vs t on a semilog scale and draw a smooth curve through the points
2. Extend the curve to intercept Y axis at t = 0 and read qi
3. Select the other end point as (t2, q2)
4. Determine the coordinates of middle point on the smooth curve that corresponds to (t1,
q1) with q1 obtained as follows: . The corresponding value of t1 is read from the curve.
5. Solve the following expression iteratively for b as:
The Newton-Raphson iterative method can be employed to solve the above non-linear
expression by using following recursion technique:
Hyperbolic Decline (0<b<1)
Where:
Starting with initial value of b = 0.5, the method will usually converge after 4-5 iterations
6. Solve for Di using calculated value of b from step 5 and coordinates ( t2, q2) to give:
Hyperbolic Decline (0<b<1)
Effect of Restricted Rate
3. Regardless the type of decline, calculate the total time of the restricted flow rate from
4. Generate well performance profile as a function of time with the help of appropriate equations