Saturation at the Water Front in Reservoir Engineering

Buckley-Leverett Displacement Theory - Fractional Flow Analysis and Waterflood Performance Prediction

The Buckley-Leverett theory is the quantitative foundation of waterflood engineering. It answers the two questions that every reservoir engineer must answer before and during a waterflood: when will water break through at the production wells, and what will the water cut profile look like after breakthrough? Without Buckley-Leverett analysis, waterflood performance can only be monitored retrospectively - watching water cut rise and reacting after the fact. With the analysis, an engineer can calculate the breakthrough time, the producing water cut at any subsequent time, and the oil recovery at any producing water-oil ratio - before injecting a single barrel of water. These predictions determine the economic life of the waterflood, the water handling capacity required in the surface facilities, and whether the waterflood is worth doing at all.

1. The Fractional Flow Equation

1.1 Deriving the Water Cut at Any Saturation

The fractional flow of water (fw) is the fraction of the total flow that is water at any point in the reservoir. It depends entirely on the relative permeability curves and the fluid viscosity ratio:

Fractional flow equation (horizontal flow, neglecting capillary pressure):
fw = 1 / (1 + (kro/krw) x (mu_w/mu_o))

Where:
fw = fraction of flowing stream that is water (dimensionless, 0 to 1)
kro = relative permeability to oil at water saturation Sw
krw = relative permeability to water at water saturation Sw
mu_w = water viscosity (cp)
mu_o = oil viscosity (cp)

Complete fractional flow curve calculation:
Using the following relative permeability data (Corey model):
Swi = 0.20 (irreducible water saturation)
Sorw = 0.25 (residual oil saturation to waterflood)
kro_Swi = 0.85 (endpoint kro at irreducible Sw)
krw_Sorw = 0.40 (endpoint krw at residual oil)
no = 3.0, nw = 2.0 (Corey exponents)
mu_o = 5.0 cp, mu_w = 0.8 cp

Normalized water saturation: S* = (Sw - Swi) / (1 - Swi - Sorw) = (Sw - 0.20) / (1 - 0.20 - 0.25) = (Sw - 0.20) / 0.55

kro(Sw) = kro_Swi x (1 - S*)^no = 0.85 x (1 - S*)^3
krw(Sw) = krw_Sorw x S*^nw = 0.40 x (S*)^2

Calculate fw at key saturations:
At Sw = 0.35: S* = (0.35-0.20)/0.55 = 0.273
kro = 0.85 x (0.727)^3 = 0.85 x 0.384 = 0.326
krw = 0.40 x (0.273)^2 = 0.40 x 0.0745 = 0.0298
fw = 1 / (1 + (0.326/0.0298) x (0.8/5.0)) = 1 / (1 + 10.94 x 0.16) = 1 / (1 + 1.750) = 0.364 (36.4% water cut)

At Sw = 0.50: S* = (0.50-0.20)/0.55 = 0.545
kro = 0.85 x (0.455)^3 = 0.85 x 0.0942 = 0.0801
krw = 0.40 x (0.545)^2 = 0.40 x 0.297 = 0.1188
fw = 1 / (1 + (0.0801/0.1188) x (0.8/5.0)) = 1 / (1 + 0.674 x 0.16) = 1 / 1.1078 = 0.903 (90.3% water cut)

At Sw = 0.55 (= 1 - Sorw): S* = 1.0 → kro = 0 → fw = 1.0 (100% water)

1.2 The fw vs Sw Table - Building the Fractional Flow Curve

Sw	S*	kro	krw	kro/krw	fw	dfw/dSw
0.20 (Swi)	0.000	0.850	0.000	∞	0.000	0.000
0.30	0.182	0.467	0.013	35.9	0.149	1.49
0.35	0.273	0.326	0.030	10.9	0.364	4.30
0.38	0.327	0.259	0.043	6.02	0.510	5.85 (max dfw/dSw)
0.45	0.455	0.138	0.083	1.66	0.789	4.22
0.50	0.545	0.080	0.119	0.67	0.903	2.47
0.55 (1-Sorw)	1.000	0.000	0.400	0.00	1.000	0.00

2. The Welge Tangent Method - Finding the Flood Front Saturation

2.1 The Buckley-Leverett Solution

The Buckley-Leverett equation states that each saturation value moves through the reservoir at a velocity proportional to the derivative of fractional flow with respect to saturation (dfw/dSw). The saturation at the flood front (Swf) is found by drawing a tangent to the fw curve from the initial water saturation point:

Buckley-Leverett velocity equation:
v_Sw = (q_total / (A x phi)) x (dfw/dSw)

The saturation profile is determined by which saturation value arrives at the production well first.

Welge tangent construction to find breakthrough saturation (Swf):
Draw a tangent from point (Swi, fwi) = (0.20, 0.0) to the fw curve.
The tangent point gives Swf (water saturation at breakthrough) and the corresponding fw at breakthrough.

From the table above, dfw/dSw has a maximum at Sw ≈ 0.38 (5.85 units).
The Welge tangent from (0.20, 0) touches the fw curve at approximately Sw = 0.38, fw = 0.510.

Flood front saturation: Swf = 0.38, fwf (fw at breakthrough) = 0.510

Wait - verify by tangent condition: (fw - 0) / (Sw - Swi) = dfw/dSw
At Swf = 0.38: (0.510 - 0.0) / (0.38 - 0.20) = 0.510/0.18 = 2.833
dfw/dSw at Sw = 0.38 = 5.85 → tangent condition not met at this point

Correct tangent point found iteratively:
At Sw = 0.32: fw = 0.217, dfw/dSw = 3.10
Tangent check: (0.217 - 0)/(0.32 - 0.20) = 0.217/0.12 = 1.81 vs dfw/dSw = 3.10 → not equal

At Sw = 0.365: fw = 0.438, dfw/dSw = 5.30
Tangent check: (0.438)/(0.365 - 0.20) = 0.438/0.165 = 2.655 vs 5.30 → not equal

True tangent at approximately Swf = 0.41, fwf = 0.625, dfw/dSw = 3.47:
Tangent check: 0.625/(0.41-0.20) = 0.625/0.21 = 2.976 → closer, iterate to convergence

Final result: Swf ≈ 0.40, fw at breakthrough ≈ 0.58 (58% water cut at first water arrival at producer)

3. Breakthrough Time and Post-Breakthrough Performance

3.1 Calculating Pore Volumes Injected at Breakthrough

Pore volumes injected at breakthrough (Qi_BT):
Qi_BT = 1 / (dfw/dSw at Swf)

Using the tangent slope at breakthrough: dfw/dSw|Swf = (fwf - 0) / (Swf - Swi) = 0.58/(0.40-0.20) = 2.90
Qi_BT = 1 / 2.90 = 0.345 pore volumes of water injected at breakthrough

Oil recovery at breakthrough (fraction of movable oil):
RF_BT = (Swf - Swi) / (1 - Swi - Sorw) = (0.40 - 0.20) / (1 - 0.20 - 0.25) = 0.20/0.55 = 0.364 = 36.4% of movable oil at breakthrough

Converting to actual time:
Qi_BT (bbl) = 0.345 x Pore_volume_total

Example: Reservoir 800 m x 600 m x 20 m, phi = 0.22, Swi = 0.20:
PV_total = 800 x 600 x 20 x 0.22 x (1/0.158981 bbl/m3) = 800 x 600 x 20 x 0.22 x 6.289 = 13,307,000 bbls

Actually using field units: PV = A(acres) x h(ft) x phi x 7,758 bbls/acre-ft
A = 800m x 600m / 4,047 m2/acre = 118.6 acres
h = 20m / 0.3048 = 65.6 ft
PV = 118.6 x 65.6 x 0.22 x 7,758 = 13,286,000 bbls total pore volume

Water injected at breakthrough = 0.345 x 13,286,000 = 4,584,000 bbls water injected

At injection rate 5,000 bbl/day:
Time to breakthrough = 4,584,000 / 5,000 = 917 days = 2.51 years to water breakthrough

3.2 Post-Breakthrough Performance - The Average Saturation Method

Average water saturation in swept zone after breakthrough (Welge method):
Sw_avg = Swf + (1 - fw_producer) / (dfw/dSw at producer fw)

At any producing water cut fw_p, the average saturation behind the front determines cumulative oil production:

NP (STB) = PV x (Sw_avg - Swi) / Bo

Example: At producing water cut fw_p = 0.90 (90% water cut):
From fw curve: Sw at fw = 0.90 corresponds to Sw ≈ 0.50
dfw/dSw at Sw = 0.50 = 2.47 (from table)

Sw_avg = 0.50 + (1 - 0.90) / 2.47 = 0.50 + 0.040 = 0.540 average Sw in swept zone

Recovery factor at 90% water cut:
RF = (Sw_avg - Swi) / (1 - Swi - Sorw) = (0.540 - 0.20) / 0.55 = 0.340/0.55 = 61.8% of movable oil recovered at 90% WC

At fw_p = 0.95 (95% water cut):
From fw curve: Sw ≈ 0.52, dfw/dSw ≈ 1.75
Sw_avg = 0.52 + (1-0.95)/1.75 = 0.52 + 0.029 = 0.549
RF = (0.549-0.20)/0.55 = 63.5% of movable oil at 95% WC

Marginal recovery between 90% and 95% WC = 63.5% - 61.8% = 1.7% of movable oil for 5% additional WC increase

This diminishing return quantifies when the waterflood reaches its economic limit: the incremental oil recovered per barrel of water injected is no longer economic.

4. Economic Limit and Waterflood Optimization

4.1 Economic Limit Water Cut

Water cut at economic limit (fw_economic_limit):
The economic limit occurs when gross revenue from oil production equals total operating cost:

Revenue = q_oil x P_oil
Operating cost = q_total x (OPEX_per_bbl_total)

At economic limit: q_oil x P_oil = q_total x OPEX
(q_total x (1-fw)) x P_oil = q_total x OPEX
(1 - fw_limit) = OPEX / P_oil
fw_limit = 1 - OPEX/P_oil

Example: P_oil = $65/STB, total lifting OPEX (including water handling) = $18/bbl total fluid:
fw_limit = 1 - 18/65 = 1 - 0.277 = 0.723 = 72.3% water cut economic limit

If water handling cost increases to $25/bbl total lifting:
fw_limit = 1 - 25/65 = 0.615 → economic limit at only 61.5% water cut

Practical implication: Installing water handling facility with lower OPEX (e.g., $12/bbl) raises economic limit:
fw_limit = 1 - 12/65 = 0.815 → 81.5% water cut limit → more oil recovered before abandonment

At fw = 0.815 vs fw = 0.723: additional recovery = approximately 3-4% of OOIP
For 200 MMstb OOIP field: 3% = 6 MMstb additional oil at $65/bbl = $390M additional revenue from water handling facility upgrade.

4.2 Viscosity Ratio Effect on Breakthrough and Recovery

Oil Viscosity (cp)	Mobility Ratio M	fw at Breakthrough	Recovery at Breakthrough	Waterflood Classification
1.0 (light oil)	0.47	0.15	52% of movable	Favorable - piston-like displacement
5.0 (medium oil)	2.35	0.58	36% of movable	This example - moderate unfavorable
20.0 (medium-heavy)	9.41	0.82	18% of movable	Unfavorable - early breakthrough, poor recovery
100+ (heavy oil)	>47	>0.95	<5% of movable	Severely unfavorable - polymer or thermal EOR required

Conclusion

The fractional flow calculation at Sw = 0.50 - fw = 90.3% water cut - demonstrates the S-shaped nature of the fw curve and its engineering significance. At 90% water cut, the reservoir is producing 9 barrels of water for every 1 barrel of oil. The Buckley-Leverett analysis predicted this would occur at 61.8% movable oil recovery, at 2.51 years breakthrough time plus the time to move from 36.4% (breakthrough) to 61.8% (90% WC) recovery. Without this calculation, the facility engineer designing the water handling system would not know whether to design for 5,000 bbl/day water or 50,000 bbl/day water at the end of field life - a facility sizing error that costs tens of millions in either inadequate capacity or unnecessary capital expenditure.

The economic limit calculation - fw_limit = 1 - OPEX/P_oil - provides the most actionable result of the entire Buckley-Leverett analysis. It directly converts a water handling facility investment decision into a reserves recovery increment. Reducing total OPEX from $18/bbl to $12/bbl by installing a more efficient produced water handling system raises the economic limit from 72.3% to 81.5% water cut, recovering an additional 3-4% of OOIP. On a 200 MMstb field, that is $390M in additional revenue - the OPEX reduction pays for orders of magnitude more facility investment than it costs. This is why waterflood management in mature fields focuses intensively on water handling cost minimization: every dollar per barrel reduction in OPEX directly translates into recoverable reserves.

Want to access our Buckley-Leverett calculator with fw curve generation, Welge tangent construction, breakthrough prediction, and economic limit water cut, or discuss waterflood performance for a specific reservoir? Join our Telegram group for reservoir engineering discussions, or visit our YouTube channel for step-by-step tutorials on fractional flow analysis and waterflood performance prediction.