Types of Casing: Choosing the Right Casing for Every Drilling Scenario

Reservoir Pressure Management - Decline Curve Analysis, Material Balance, and Enhanced Recovery Screening

Every reservoir has a finite pressure drive, and every production decision is simultaneously a reservoir management decision. The engineer who increases the choke size to maximize short-term rate without calculating the impact on reservoir pressure and ultimate recovery is making an irreversible decision with long-term consequences. The engineer who understands the material balance of the reservoir - how much pressure decline corresponds to how much produced volume - can optimize the production rate to balance rate of return against ultimate recovery. This guide covers the quantitative tools that connect production decisions to reservoir outcomes: decline curve analysis for forecasting, material balance for understanding drive mechanisms, and enhanced recovery screening for extending field life beyond primary production.

1. Decline Curve Analysis - Production Forecasting

1.1 Arps Decline Curve Models

Arps (1945) established the empirical decline curve framework that remains the standard production forecasting tool in the industry. Three models cover the observed range of decline behavior:

Exponential decline (b = 0):
q(t) = qi x exp(-Di x t)
EUR = qi / Di

Hyperbolic decline (0 < b < 1):
q(t) = qi / (1 + b x Di x t)^(1/b)
EUR = qi^b / ((1-b) x Di) x (qi^(1-b) - q_abandonment^(1-b))

Harmonic decline (b = 1):
q(t) = qi / (1 + Di x t)
EUR = (qi / Di) x ln(qi / q_abandonment)

Where:
qi = initial rate (bbl/day or Mscf/day)
Di = initial decline rate (fraction/year or fraction/month)
b = hyperbolic exponent (dimensionless)
t = time

Worked example - hyperbolic decline:
Well data: qi = 1,200 bbl/day, Di = 0.85/year (85% annual decline), b = 0.6, q_abandonment = 20 bbl/day

Rate after 3 years: q(3) = 1,200 / (1 + 0.6 x 0.85 x 3)^(1/0.6)
= 1,200 / (1 + 1.53)^1.667 = 1,200 / (2.53)^1.667 = 1,200 / 4.76 = 252 bbl/day after 3 years

Time to abandonment (q = 20): 20 = 1,200/(1 + 0.51t)^1.667
(1 + 0.51t)^1.667 = 60 → 1 + 0.51t = 60^0.6 = 12.09 → t = 11.09/0.51 = 21.7 years well life

EUR ≈ 1,200^0.6 / ((1-0.6) x 0.85) x (1,200^0.4 - 20^0.4) = calculated as approximately 1.85 MMbbl EUR

1.2 Decline Exponent (b) Interpretation

b Value	Physical Meaning	Reservoir Type	EUR Ranking vs Exponential
b = 0 (exponential)	Constant fractional decline rate. Decline rate does not change with time.	Solution gas drive, single-layer depletion	Baseline
0 < b < 0.5	Decline rate decreases over time - slowing decline. Production lasts longer than exponential prediction.	Gravity drainage, solution gas with some compaction drive	Higher
b = 0.5-1.0	Significant decline rate decrease over time. Production tail is much longer.	Unconventional shale/tight gas (b often 1.2-2.0 early, transitioning to <1 long term)	Significantly higher
b > 1.0	Theoretically possible in transient flow (fracture-dominated early flow). Not physical at long term - must switch to exponential or b<1 as boundary-dominated flow begins.	Early transient flow in fractured/unconventional wells	Much higher (often over-estimated if not corrected)

2. Material Balance Equation - Understanding the Reservoir Drive

2.1 The Havlena-Odeh Material Balance Plot

The material balance equation (MBE) is the volumetric accounting equation for a reservoir: the expansion of fluids and rock due to pressure depletion must equal the cumulative production. The Havlena-Odeh linearization transforms the general MBE into a straight-line plot that identifies the drive mechanism and quantifies the original fluids in place:

Simplified material balance (undersaturated oil reservoir):
Np x Bo + Wp x Bw = N x Boi x ct x dP

Where:
Np = cumulative oil produced (STB)
Bo = oil formation volume factor (RB/STB)
Wp = cumulative water produced (STB)
N = OOIP (STB) - the unknown to solve for
Boi = initial Bo
ct = total compressibility (1/psi)
dP = pressure drop (Pi - current P)

Havlena-Odeh form - plot F vs Eo:
F = Np x Bo + Wp x Bw (production voidage)
Eo = Boi x ct x dP (oil + rock expansion term)

Plot F (Y-axis) vs Eo (X-axis): straight line through origin with slope = N (OOIP)

Example calculation:
At a given time: Np = 2,500,000 STB, Wp = 180,000 STB, Bo = 1.285 RB/STB, Bw = 1.02 RB/STB
F = 2,500,000 x 1.285 + 180,000 x 1.02 = 3,212,500 + 183,600 = 3,396,100 RB

If dP = 485 psi, Boi = 1.275, ct = 18 x 10^-6 psi^-1:
Eo = 1.275 x 18e-6 x 485 = 1.275 x 0.00873 = 0.01113 RB/STB

N = F/Eo = 3,396,100 / 0.01113 = 305,000,000 STB = 305 MMbbl OOIP

If volumetric estimate was 280 MMbbl: material balance suggests 9% higher OOIP → re-evaluate well data or model.

2.2 Drive Mechanism Identification from Material Balance

Drive Mechanism	Havlena-Odeh Plot Behavior	Recovery Factor Range	Implication for EOR
Rock/liquid expansion (undersaturated)	Straight line through origin with F vs Eo	5-15%	Low recovery - pressure maintenance by water injection is cost-effective at early stage
Solution gas drive	Curve - F/Eo changes with time as gas evolves	5-30%	Gas injection to maintain pressure above bubble point dramatically improves recovery
Gas cap drive	F vs (Eo + m x Eg) gives straight line - must include gas cap expansion term m x Eg	20-40%	Maintain gas cap integrity. Avoid producing from gas cap. Horizontal wells below GOC.
Water drive (aquifer)	F vs Eo deviates upward from N-slope line → aquifer influx We supplements expansion	30-60%	Strong water drive can provide good recovery without injection. Manage water coning.

3. Enhanced Oil Recovery - Screening and Selection

3.1 EOR Screening Criteria

Enhanced oil recovery extends production beyond what primary and secondary recovery can achieve by introducing energy or materials into the reservoir that improve displacement efficiency or reduce residual oil saturation. Screening criteria match the EOR method to the reservoir properties:

EOR Method	Oil API Gravity	Viscosity (cp)	Permeability (md)	Depth (ft)	Mechanism
Miscible CO2 injection	>26°	<10	>1	>2,500	CO2 dissolves in oil, reduces viscosity, swells oil, achieves miscibility → eliminates capillary pressure. Recovery: 10-20% of OOIP additional.
Polymer flooding	>15°	10-150	>10	<9,000	Polymer thickens injected water → improves mobility ratio → reduces viscous fingering → better sweep efficiency. Recovery: 5-15% additional vs waterflood.
Surfactant flooding (CEOR)	>20°	<30	>5	<9,000	Ultra-low IFT reduces capillary number → mobilizes residual oil. Very high chemical cost. Recovery: 10-20% additional.
Steam injection (thermal)	<20°	>100	>200	<4,500	Steam heats reservoir → reduces oil viscosity 10-1,000x → enables flow. Required for heavy oil. Recovery: 40-60% of OOIP.
In-situ combustion (THAI)	<25°	>50	>50	<11,500	Inject air → burn small fraction of oil in reservoir → generates heat in-situ → upgrades remaining oil. Applicable to deep heavy oil where steam injection is uneconomic.

3.2 Mobility Ratio - The Key Waterflood Parameter

Mobility ratio (M):
M = (krw / mu_water) / (kro / mu_oil) = (krw x mu_oil) / (kro x mu_water)

M < 1: Favorable - water is less mobile than oil. Piston-like displacement. High sweep efficiency.
M = 1: Neutral
M > 1: Unfavorable - water more mobile than oil. Viscous fingering. Poor sweep efficiency.

Example: krw = 0.40, kro = 0.75, mu_oil = 15 cp, mu_water = 1.0 cp:
M = (0.40 x 15) / (0.75 x 1.0) = 6.0 / 0.75 = M = 8.0 → Very unfavorable

At M = 8, waterflood areal sweep efficiency ≈ 50-60% (significant unswept areas)
vs M = 1: sweep efficiency ≈ 85-90%

Polymer flooding correction:
Add polymer to water: mu_polymer_solution = 5 cp (vs 1 cp water)
M_polymer = (0.40 x 15) / (0.75 x 5.0) = 6.0 / 3.75 = M = 1.6 → Near-favorable

Polymer flooding reduces M from 8.0 to 1.6 → Sweep efficiency from 55% to approximately 80%
Recovery improvement from polymer: (0.80 - 0.55) x N x Sor_reduction = significant incremental recovery

4. Pressure Maintenance - The Economic Case for Early Injection

4.1 Timing of Water Injection - Early vs Late

Injection Timing	Reservoir Pressure Maintained?	Ultimate Recovery vs No Injection	Additional Economic Consideration
Early injection (before bubble point)	Yes - maintains above bubble point	+15-25% additional recovery	Avoids gas-oil ratio increase. Oil produced above bubble point has maximum Bo (most volume). Infrastructure cost offset by earlier plateau production.
Late injection (after pressure drops to bubble point)	Partially - arrests further decline	+8-15% additional recovery	Gas has already evolved. Residual gas saturation reduces oil relative permeability. Some of the best reservoir drive has already been wasted.
Very late injection (well into solution gas depletion)	No - too late to maintain pressure	+3-8% additional recovery	High GOR. Increased free gas saturation locks oil in reservoir. Water injection primarily provides pressure support for artificial lift rather than enhanced recovery.

Conclusion

The material balance calculation in this article - OOIP = 305 MMbbl from the Havlena-Odeh plot versus 280 MMbbl from volumetrics - illustrates the most important feature of the material balance equation: it is calibrated by actual production history and actual pressure measurements. Unlike the volumetric estimate, which depends on the accuracy of the geological model and log interpretation, the material balance result is constrained by what the reservoir has actually produced and at what pressure. When the two estimates agree, confidence in the OOIP is high. When they disagree by 9% as in this example, one of the underlying assumptions requires investigation - either the volumetric model overestimates the bulk volume or net-to-gross, or the material balance misidentifies the drive mechanism (aquifer influx not accounted for, for example).

The mobility ratio calculation - M = 8.0 for 15 cp oil in a standard waterflood, reducing to M = 1.6 with polymer - shows exactly why polymer flooding improves sweep efficiency and why screening criteria specify viscosity limits for waterflood application. At M = 8, sweep efficiency is approximately 55% regardless of injection volume, because the injected water channels through the high-permeability paths. Increasing the injection rate does not fix M = 8 - it only accelerates the channeling. Polymer reduces M to near-1 and changes the physics from channeling-dominated displacement to piston-like displacement. The incremental recovery from polymer is not a chemical curiosity - it is the mathematical consequence of reducing the mobility ratio from 8.0 to 1.6 in a reservoir where sweep efficiency is the primary recovery limit.

Want to access our reservoir management calculator with Arps decline curve fitting, Havlena-Odeh material balance, mobility ratio, and EOR screening matrix, or discuss production optimization for a specific reservoir type? Join our Telegram group for reservoir engineering discussions, or visit our YouTube channel for step-by-step tutorials on decline curve analysis, material balance, and EOR selection.