Capillary Pressure and Saturation-Height Functions - Fluid Contacts, Initial Water Saturation, and Reservoir Characterization

Capillary Pressure and Saturation - Height Functions - Fluid Contacts, Initial Water Saturation, and Reservoir Characterization

Capillary pressure is the mechanism that determines the initial distribution of fluids in a reservoir before production begins. It answers the questions that govern every volumetric calculation: where exactly is the free water level (FWL), how does water saturation vary from the FWL up to the top of the reservoir, and what is the initial water saturation at any depth within the pay zone? These questions cannot be answered from resistivity logs alone without a capillary pressure model - the Archie equation applied to a resistivity log gives Sw at a point, but it cannot explain why Sw at 100 ft above the FWL is different from Sw at 20 ft above the FWL in the same rock type. Capillary pressure provides the physical model that connects these observations. The saturation-height function derived from capillary pressure measurements allows the engineer to predict Sw at any height above the FWL for any rock type in the reservoir and therefore to calculate OOIP with a spatial Sw distribution rather than a single average assumption.

1. Capillary Pressure Physics

1.1 The Young-Laplace Equation - What Creates Capillary Pressure

Capillary pressure arises at any curved interface between two immiscible fluids. In a reservoir pore throat, the curvature of the oil-water interface generates a pressure difference between the two phases - the non-wetting phase (oil in a water-wet rock) is at higher pressure than the wetting phase (water):

Young-Laplace equation for capillary pressure:
Pc = 2 x sigma x cos(theta) / r

Where:
Pc = capillary pressure (psi or Pa)
sigma = interfacial tension between oil and water (dynes/cm or mN/m)
theta = contact angle (degrees) - measures wettability. theta = 0° fully water-wet, theta = 180° fully oil-wet
r = pore throat radius (cm or m)

Key relationships:
Small pore throat (small r) → high Pc → more energy required for oil to enter → oil stays out of smallest pores
Large pore throat (large r) → low Pc → oil enters easily
Oil-wet rock (theta > 90°) → cos(theta) is negative → Pc sign reverses → oil is the wetting phase

Example: Water-wet sandstone, sigma = 30 dynes/cm, theta = 35°, r = 5 micrometers = 5 x 10^-4 cm:
Pc = 2 x 30 x cos(35°) / (5 x 10^-4)
= 2 x 30 x 0.819 / 5 x 10^-4
= 49.14 / 5 x 10^-4 = 98,280 dynes/cm2
= 98,280 / 68,948 psi/dyne/cm2 = 1.43 psi capillary pressure at this pore throat

Converting Pc to height above free water level:
Pc (psi) = (rho_water - rho_oil) x 0.052 x H
H (ft above FWL) = Pc / ((rho_water - rho_oil) x 0.052)

For rho_water = 8.7 ppg, rho_oil = 7.2 ppg:
delta_rho = 8.7 - 7.2 = 1.5 ppg
H = 1.43 / (1.5 x 0.052) = 1.43 / 0.078 = 18.3 ft above FWL where this pore throat fills with oil

1.2 The Capillary Pressure Curve : Laboratory Measurement

The Pc-Sw curve is measured in the laboratory by incrementally increasing the pressure of a non-wetting phase (mercury or oil) and measuring the saturation of the non-wetting phase that enters the core plug at each pressure step. The curve shape characterizes the pore size distribution of the rock:

Pc Curve Feature	Physical Meaning	Reservoir Quality Indicator
Entry pressure (Pd)	Minimum pressure required for non-wetting phase to enter the rock. Determined by the largest connected pore throat.	Low Pd: Good reservoir quality (large pore throats). High Pd: Tight rock (small pore throats).
Plateau (flat section)	Large pore throats of similar size fill over a small pressure range. Indicates good sorting.	Long flat plateau: Well-sorted, high permeability rock. Short plateau: Poorly sorted rock.
Transition zone length	Vertical extent of the zone where Sw decreases from ~1.0 to Swi. Reflects range of pore throat sizes.	Short transition zone: Uniform pore size, sharp fluid contact. Long transition zone: Wide pore size distribution, gradual Sw decline with height.
Irreducible water saturation (Swi)	Minimum water saturation achievable - water in the smallest pores that oil cannot displace regardless of pressure.	Low Swi: Clean, well-sorted rock with fewer small pores. High Swi: Clayey or poorly sorted rock with many small pores retaining water.

2. Converting Lab Pc to Reservoir Conditions

2.1 The Leverett J-Function: Normalizing Pc for Reservoir Use

Laboratory Pc measurements are made with mercury-air (for MICP) or oil-brine systems at room temperature and pressure. Before using these curves to predict initial Sw in the reservoir, they must be converted to reservoir fluid-rock conditions using the Leverett J-function:

Leverett J-function (dimensionless capillary pressure):
J(Sw) = Pc x sqrt(k/phi) / (sigma x cos(theta))

Where Pc in consistent units (psi or Pa), k in md, phi as fraction, sigma in dynes/cm or mN/m

Converting from laboratory to reservoir conditions:
Pc_reservoir = Pc_lab x (sigma_res x cos(theta_res)) / (sigma_lab x cos(theta_lab))

For mercury-air to oil-brine conversion:
sigma_mercury-air = 480 dynes/cm, theta_mercury-air = 140°
sigma_oil-brine_reservoir = 25 dynes/cm (typical), theta_oil-brine = 35° (water-wet)

Conversion factor = (25 x cos(35°)) / (480 x cos(140°))
= (25 x 0.819) / (480 x (-0.766))
= 20.475 / (-367.68)
= -0.0557

Since mercury is non-wetting and oil is non-wetting in water-wet rock, both work against the wetting phase - use absolute values:
Conversion factor = 20.475 / 367.68 = 0.0557

Pc_reservoir = Pc_mercury x 0.0557

Example: MICP mercury pressure = 150 psi at Sw = 0.35:
Pc_reservoir = 150 x 0.0557 = 8.36 psi reservoir capillary pressure at Sw = 0.35

Height above FWL at Sw = 0.35:
H = 8.36 / (1.5 ppg x 0.052) = 8.36 / 0.078 = 107 ft above FWL where Sw = 0.35 in this rock type

3. Saturation-Height Function: Predicting Initial Sw Throughout the Reservoir

3.1 Building the Saturation-Height Function

The saturation-height function (SHF) combines the J-function normalization with the reservoir fluid-density contrast to predict water saturation at any height above the free water level for any rock quality in the reservoir:

Saturation-height function (SHF) using J-function:
Sw(H, k, phi) = f(J) where J = Pc x sqrt(k/phi) / (sigma_res x cos(theta_res))
and Pc = H x (rho_water - rho_oil) x 0.052

So: J = H x delta_rho x 0.052 x sqrt(k/phi) / (sigma_res x cos(theta))

Example SHF calculation for three rock types at 80 ft above FWL:
Reservoir: delta_rho = 1.5 ppg, sigma_res = 25 dynes/cm, theta = 35°
At H = 80 ft: Pc_reservoir = 80 x 1.5 x 0.052 = 6.24 psi

Rock Type A (good quality): k = 150 md, phi = 0.22
J_A = 6.24 x sqrt(150/0.22) / (25 x 0.819) = 6.24 x sqrt(682) / 20.475 = 6.24 x 26.12 / 20.475 = J = 7.96
From J-function correlation: Sw_A at J = 7.96 ≈ 0.22 (22% Sw in good rock at 80 ft above FWL)

Rock Type B (moderate quality): k = 25 md, phi = 0.16
J_B = 6.24 x sqrt(25/0.16) / 20.475 = 6.24 x sqrt(156.25) / 20.475 = 6.24 x 12.5 / 20.475 = J = 3.81
From J-function: Sw_B at J = 3.81 ≈ 0.38 (38% Sw in moderate rock at same height)

Rock Type C (poor quality): k = 3 md, phi = 0.12
J_C = 6.24 x sqrt(3/0.12) / 20.475 = 6.24 x sqrt(25) / 20.475 = 6.24 x 5.0 / 20.475 = J = 1.52
From J-function: Sw_C at J = 1.52 ≈ 0.65 (65% Sw in poor rock at same height)

At the same height (80 ft above FWL), the same reservoir produces:
- 78% oil saturation in good rock (Rock A)
- 62% oil saturation in moderate rock (Rock B)
- 35% oil saturation in poor rock (Rock C)

This dramatic difference in Sw at the same height explains why OOIP calculations that use a single average Sw for the entire reservoir significantly overestimate recoverable hydrocarbons in heterogeneous reservoirs.

3.2 SHF Application to OOIP Calculation

Sw Estimation Method	Approach	Typical Error Range	When Appropriate
Single average Sw from logs	Average Archie Sw over all perforated intervals. Apply uniformly to entire reservoir volume.	±15-25% OOIP	Early exploration with few wells and no capillary pressure data. Acceptable for order-of-magnitude estimates only.
Log-derived Sw by zone	Archie Sw calculated at each well for each flow unit. Applied to respective reservoir zones.	±10-15% OOIP	Development planning with multiple wells. Adequate for most field development decisions.
Saturation-height function	Sw(H, k, phi) applied throughout reservoir volume. Sw varies continuously with height above FWL and with rock quality at each point.	±3-7% OOIP	Field development, reserves certification, facility sizing, reservoir simulation initialization. Required for heterogeneous reservoirs with significant transition zones.

4. Free Water Level vs Oil-Water Contact

4.1 The Critical Distinction Between FWL and OWC

Free Water Level (FWL): The depth at which Pc = 0. At this depth, oil and water are in equilibrium - no capillary pressure difference. This is the theoretical contact.

Oil-Water Contact (OWC): The shallowest depth at which 100% water saturation is observed on logs. Due to the entry pressure (Pd), there is a zone below the FWL where water saturation is still 100% even though oil exists above. The OWC is always AT or ABOVE the FWL.

OWC depth (TVD) = FWL depth (TVD) - H_entry
H_entry = Pd / (delta_rho x 0.052)

Example: Pd = 1.8 psi (entry pressure of reservoir rock), delta_rho = 1.5 ppg:
H_entry = 1.8 / (1.5 x 0.052) = 1.8 / 0.078 = 23.1 ft

If OWC on logs = 9,850 ft TVD:
FWL = 9,850 + 23.1 = 9,873 ft TVD (23 ft deeper than the logged OWC)

Why this matters for OOIP:
If OOIP is calculated using OWC = 9,850 ft instead of FWL = 9,873 ft as the base of the oil column:
Area of reservoir = 2,000 acres, phi = 0.20, (1-Sw) ≈ 0.75, Bo = 1.30
OOIP correction for 23 ft depth error = 7,758 x 2,000 x 23 x 0.20 x 0.75 / 1.30 = 41,200 STB = 41,200 bbls underestimate

On a large field: this error scales directly with reservoir area. At 50,000 acres: 1,030,000 STB ≈ 1 MMstb underestimate from the FWL/OWC confusion alone.

5. Capillary Pressure and EOR: IFT Reduction

5.1 Capillary Number: Quantifying the Mobilization Criterion

Capillary number (Nca) - ratio of viscous to capillary forces:
Nca = (mu x v) / (sigma x cos(theta))

Where mu = displacing fluid viscosity (cp), v = Darcy velocity (cm/sec), sigma = IFT (dynes/cm)

Residual oil is mobilized when Nca exceeds a critical threshold:
Nca_critical ≈ 10^-5 to 10^-4 (depends on rock type)

In a conventional waterflood:
mu_water = 1 cp, v = 1 ft/day = 3.53 x 10^-4 cm/sec, sigma = 30 dynes/cm, theta = 35°
Nca = (1 x 3.53e-4) / (30 x 0.819) = 3.53e-4 / 24.57 = 1.44 x 10^-5

Conventional waterflood barely reaches the lower end of the mobilization threshold → leaves Sorw in place.

With surfactant injection (ultra-low IFT):
sigma_surfactant = 0.001 dynes/cm (ultra-low IFT achieved)
Nca_surfactant = (1 x 3.53e-4) / (0.001 x 0.819) = 3.53e-4 / 8.19e-4 = 0.431

Nca_surfactant / Nca_waterflood = 0.431 / 1.44e-5 = 30,000x higher capillary number

At Nca = 0.431: Residual oil saturation approaches zero → essentially all capillary-trapped oil is mobilized. This is the engineering basis for surfactant EOR: reduce IFT by 30,000x, increase Nca by 30,000x, eliminate residual oil saturation.

Conclusion

The saturation-height function calculation in this article - Sw = 22% in good rock versus Sw = 65% in poor rock at the same height of 80 ft above the free water level - quantifies the single most important reason why heterogeneous reservoir OOIP calculations require the J-function approach rather than a simple average Sw. The 43-percentage-point difference in Sw at the same height is not an artifact of the calculation - it is the physical consequence of the Young-Laplace equation applied to rocks with different pore throat size distributions. Poor quality rock (3 md) has much smaller pore throats than good quality rock (150 md), so capillary forces hold much more water in the pore space at any given height above the free water level.

The capillary number calculation - Nca increasing 30,000x from conventional waterflood to surfactant EOR - provides the engineering explanation for why surfactant flooding works and why it requires ultra-low IFT. At Nca = 1.44 x 10^-5 in a conventional waterflood, viscous forces are barely strong enough to approach the mobilization threshold. At Nca = 0.431 with 0.001 dynes/cm surfactant IFT, viscous forces overwhelm capillary forces and residual oil saturation approaches zero. The chemical cost of achieving that 0.001 dynes/cm IFT - typically $0.50-2.00 per barrel of produced water injected - must be justified against the incremental oil recovery from Sorw mobilization. On reservoirs where Sorw = 0.25 and OOIP = 100 MMstb, the Sorw represents 25 MMstb of oil that conventional waterflood leaves in the reservoir. At $60/bbl recovery value, that is $1.5B of oil that surfactant flooding can potentially recover.

Want to access our capillary pressure toolkit with MICP-to-reservoir conversion, J-function fitting, saturation-height function generator, and capillary number calculator, or discuss initial water saturation characterization for a specific reservoir? Join our Telegram group for reservoir engineering discussions, or visit our YouTube channel for step-by-step tutorials on capillary pressure analysis and saturation-height function development.