Pore Pressure Prediction: Compaction Trends, Seismic Velocity Methods, and Mud Weight Window Design

Pore Pressure Prediction: Compaction Trends, Seismic Velocity Methods, and Mud Weight Window Design

Pore pressure prediction is the pre-drill engineering process that determines the mud weight required to maintain wellbore stability at each depth in the planned well. An accurate pore pressure prediction is the difference between drilling a well safely to total depth and encountering a kick, a blowout, or a catastrophic wellbore collapse that results in the loss of the well and potentially the drilling unit. The challenge is that pore pressure cannot be directly measured before the well is drilled - it must be predicted from geophysical data (seismic velocity, gravity) and geological analogs, then refined in real time as the well is drilled using drilling parameters and wireline log data. This guide covers the quantitative methods for pore pressure prediction: the compaction trend analysis that identifies overpressured zones from sonic or resistivity log departures, the Eaton method that converts these departures to pore pressure gradients, and the drilling engineering application that converts pore pressures to the mud weight window that governs casing program design.

1. Normal Compaction Trend: The Baseline for Overpressure Detection

1.1 What Creates Normal and Abnormal Pore Pressure

In normally pressured sediments, pore pressure at any depth equals the hydrostatic pressure of a column of formation water from the surface to that depth. As sediments compact during burial, pore fluids are expelled upward. Normal pore pressure gradient in marine environments is approximately 0.45-0.465 psi/ft (8.65-8.95 ppg equivalent).

Overpressure develops when compaction is faster than fluid expulsion - fluids become trapped in the pore space and the rock compacts less than normal for that burial depth. The most common cause is rapid burial of fine-grained sediments (shales) where low permeability prevents timely fluid expulsion. The undercompacted shale retains more porosity and higher pore fluid pressure than a normally compacted shale at the same depth:

Terzaghi effective stress relationship:
Sigma_eff = Sigma_v - Pp

Where:
Sigma_eff = effective (grain-to-grain) stress (psi)
Sigma_v = vertical (overburden) stress (psi)
Pp = pore pressure (psi)

In a normally pressured zone at 10,000 ft (OBG = 18.5 ppg, PP = 8.9 ppg):
Sigma_v = 18.5 x 0.052 x 10,000 = 9,620 psi
Pp = 8.9 x 0.052 x 10,000 = 4,628 psi
Sigma_eff = 9,620 - 4,628 = 4,992 psi effective stress (normal)

In an overpressured zone at same depth (PP = 14.5 ppg):
Pp_over = 14.5 x 0.052 x 10,000 = 7,540 psi
Sigma_eff_over = 9,620 - 7,540 = 2,080 psi effective stress (reduced)

The overpressured shale at 10,000 ft has effective stress equivalent to a normally pressured shale at approximately 4,500 ft burial depth - it is "undercompacted" relative to its actual depth.

Key diagnostic: The overpressured shale will have the same sonic velocity, resistivity, and porosity as a normally compacted shale at 4,500 ft - not at 10,000 ft. This departure from the normal compaction trend is the detection method.

1.2 Establishing the Normal Compaction Trend (NCT)

Normal compaction trend from sonic (DtC) log in shale intervals:
DtC (normally compacted shale) decreases with depth as shale compacts and pore space reduces.
Typical empirical relationship: DtC_normal = DtC_surface x exp(-c x TVD)

Where DtC_surface = transit time at surface (typically 200 microsec/ft for marine shales), c = compaction coefficient (typically 0.00015 - 0.00025 per ft), TVD in feet

Example: DtC_surface = 200 microsec/ft, c = 0.0002/ft

DtC_normal at 5,000 ft = 200 x exp(-0.0002 x 5,000) = 200 x exp(-1.0) = 200 x 0.368 = 73.6 microsec/ft
DtC_normal at 10,000 ft = 200 x exp(-0.0002 x 10,000) = 200 x exp(-2.0) = 200 x 0.135 = 27.1 microsec/ft
DtC_normal at 8,000 ft = 200 x exp(-1.6) = 200 x 0.202 = 40.4 microsec/ft

Actual measured DtC in shale at 8,000 ft = 68 microsec/ft
68 >> 40.4 → shale is slower than normal → undercompacted → overpressure detected at 8,000 ft

2. Eaton Method - Converting Log Departure to Pore Pressure

2.1 Eaton Pore Pressure Equation

The Eaton method is the industry standard for converting sonic or resistivity log departures from the normal compaction trend into pore pressure estimates. It was developed empirically from Gulf of Mexico overpressure data and is widely applied globally:

Eaton equation for pore pressure from sonic transit time:
PP_gradient (psi/ft) = OBG - (OBG - PPn) x (DtC_normal / DtC_observed)^3

Where:
OBG = overburden gradient (psi/ft) at the depth of interest
PPn = normal pore pressure gradient (psi/ft) ≈ 0.45 psi/ft (8.65 ppg)
DtC_normal = normal compaction trend sonic transit time at that depth (microsec/ft)
DtC_observed = actual measured sonic transit time in shale (microsec/ft)
Exponent = 3 for sonic (Eaton's original; some regions use 2 or 4)

Eaton equation from resistivity:
PP_gradient = OBG - (OBG - PPn) x (Rt_observed / Rt_normal)^1.2

Where Rt = shale resistivity (ohm-m), exponent = 1.2 for resistivity

Example calculation at 8,000 ft TVD:
OBG at 8,000 ft = 1.0 psi/ft (19.23 ppg, typical offshore)
Wait - typical OBG: OBG = rho_bulk x 0.052 ≈ 18.5-19.5 ppg → OBG = 0.962 psi/ft using 18.5 ppg
PPn = 0.468 psi/ft (9.0 ppg)
DtC_normal at 8,000 ft = 40.4 microsec/ft (from NCT above)
DtC_observed = 68 microsec/ft (measured in overpressured shale)

PP_gradient = 0.962 - (0.962 - 0.468) x (40.4/68)^3
= 0.962 - 0.494 x (0.594)^3
= 0.962 - 0.494 x 0.2097
= 0.962 - 0.1036
= 0.858 psi/ft = 16.5 ppg equivalent pore pressure at 8,000 ft

This is severely overpressured (16.5 ppg vs normal 9.0 ppg). The mud weight at 8,000 ft must be at least 16.5 ppg equivalent, with a safety margin of +0.3-0.5 ppg.

2.2 Pore Pressure Prediction from Seismic Velocity (Pre-Drill)

Before the well is drilled, seismic interval velocity from velocity analysis (PSDM - pre-stack depth migration) provides the only direct subsurface measurement available for pore pressure prediction. The seismic velocity is converted to sonic transit time and then the Eaton method is applied:

Seismic velocity to DtC conversion:
DtC (microsec/ft) = 1,000,000 / V_interval (ft/sec)

Example: Seismic interval velocity in target shale = 6,200 ft/sec
DtC_seismic = 1,000,000 / 6,200 = 161.3 microsec/ft

Normal DtC at this depth from NCT = 40.4 microsec/ft (same depth 8,000 ft)
161.3 >> 40.4 → very large departure → severe overpressure predicted

PP_gradient = 0.962 - 0.494 x (40.4/161.3)^3 = 0.962 - 0.494 x (0.250)^3
= 0.962 - 0.494 x 0.01563 = 0.962 - 0.00772 = 0.954 psi/ft ≈ 18.3 ppg pre-drill pore pressure prediction

Seismic velocity uncertainty in pore pressure prediction:
Seismic velocity analysis has ±5-10% accuracy in interval velocity estimation.
At V = 6,200 ft/sec with ±10% error: V_range = 5,580 to 6,820 ft/sec
DtC range = 146.6 to 179.2 microsec/ft → PP range = 0.950 to 0.957 psi/ft (small difference due to flat Eaton curve at high DtC)

In this severely overpressured case, the seismic uncertainty does not change the conclusion (well-above-normal pressure). But in the transition zone between normal and overpressured, ±10% velocity uncertainty can mean ±1-2 ppg mud weight uncertainty.

3. Mud Weight Window Design from Pore Pressure and Fracture Gradient

3.1 Building the Full Pressure Profile

Depth (ft TVD)	PP (ppg)	FG (ppg)	Window (ppg)	MW Required (ppg)	Casing Required?
0-2,000	8.7	14.5	5.8	9.5	Conductor at 500 ft. Surface casing at 2,000 ft.
2,000-6,500	9.0	15.2	6.2	9.7	Drill to 6,500 ft with 9.7 ppg MW. No casing needed within this interval.
6,500-8,000	9.0-16.5	17.2	0.7	16.8	NARROW WINDOW. Intermediate casing REQUIRED before penetrating 7,000 ft where PP > 14.5 ppg. Cannot drill open hole at 9.7 ppg through this transition.
8,000-10,500 (TD)	16.5	18.0	1.5	17.0	Drill with 17.0 ppg MW. Production casing at TD. ECD management critical (window = 1.5 ppg).

3.2 Real-Time Pore Pressure Monitoring While Drilling

The pre-drill pore pressure prediction is updated in real time as the well is drilled using drilling parameters and MWD/LWD measurements. The most sensitive real-time indicators:

Real-Time Indicator	Normal Behavior	Overpressure Warning Sign
dc-exponent (d-exponent corrected for MW)	Increases with depth as rock compacts (harder drilling)	dc-exponent decreases → undercompacted rock → overpressure onset
LWD sonic velocity	Velocity increases (DtC decreases) with depth	DtC increases while depth increases → shale slower than normal → overpressure
LWD resistivity	Resistivity increases with depth as porosity decreases	Resistivity decreases or plateaus → abnormally high porosity → overpressure
Pit gain / flow check	No pit gain when pumps running or during flow check	Pit gain → formation fluid influx → mud weight insufficient → SHUT IN WELL
Gas shows (C1, C2, C3)	Background gas level stable	Gas shows increasing toward top of reservoir → approaching pay → increase MW

4. Wellbore Stability - Collapse and Breakout Prevention

4.1 Minimum Mud Weight for Shale Stability

Mohr-Coulomb shale stability criterion (minimum MW to prevent breakout):
In a vertical well:
MW_min (ppg) = (3 x Sh_min - Sh_max - UCS - Pp) / (2 x 0.052 x TVD) + Pp/(0.052 x TVD)

Simplified: MW_min = Pp_gradient + [(3 x Sh_min - Sh_max - UCS - Pp)] / (0.052 x TVD x 2)

Example at 8,000 ft: Sh_min = 15.8 ppg, Sh_max = 17.2 ppg, UCS = 3,500 psi, Pp = 16.5 ppg:
Sh_min_psi = 15.8 x 0.052 x 8,000 = 6,573 psi
Sh_max_psi = 17.2 x 0.052 x 8,000 = 7,155 psi
Pp_psi = 16.5 x 0.052 x 8,000 = 6,864 psi

MW_min = (3 x 6,573 - 7,155 - 3,500 - 6,864) / (2 x 0.052 x 8,000) + 16.5
= (19,719 - 7,155 - 3,500 - 6,864) / 832 + 16.5
= 2,200 / 832 + 16.5 = 2.64 + 16.5 = 19.14 ppg minimum MW to prevent shale breakout

Problem: Fracture gradient = 18.0 ppg and minimum stability MW = 19.14 ppg → impossible to drill this interval safely with any conventional mud weight!

Engineering response options:
1. Set casing above this interval and reduce mud weight before entering
2. Use OBM/SBM to reduce effective stress on wellbore (reduces required MW by 0.5-2 ppg)
3. Reduce Sh_max by chemical treatment or wellbore orientation change
4. Drill quickly through the unstable zone to minimize exposure time
5. Consider alternate well path (horizontal wells in some geomechanical settings can avoid worst stability conditions)

Conclusion

The Eaton calculation in this article - converting a measured DtC of 68 microsec/ft versus the normal compaction trend value of 40.4 microsec/ft at 8,000 ft into a pore pressure of 16.5 ppg - demonstrates the complete workflow from log observation to mud weight design. The 68 microsec/ft measurement is not anomalous-looking by itself: it is a reasonable sonic transit time for a shale. What makes it diagnostic is its position relative to the normal compaction trend. Without the NCT established from shallower normally-pressured shales, the 68 microsec/ft value has no context and cannot be interpreted as overpressure. This is why pore pressure prediction requires the NCT to be established first, and why wells drilled without adequate shale data to establish the NCT produce unreliable pre-drill pore pressure predictions.

The wellbore stability calculation that produces MW_min = 19.14 ppg in a formation with FG = 18.0 ppg identifies the most dangerous drilling scenario in deep overpressured wells: a stability window that is physically impossible to satisfy. No mud weight simultaneously above PP (16.5 ppg), below FG (18.0 ppg), and above stability requirement (19.14 ppg) exists. The engineering responses - set casing first, use OBM, drill quickly - are not optional preferences. They are the only physical options when the geomechanical constraints eliminate any viable mud weight. This analysis must be performed before spud, not discovered at 8,000 ft when stuck pipe or lost circulation makes the point for the first time.

Want to access our pore pressure prediction toolkit with NCT establishment, Eaton method calculator, mud weight window design, and stability analysis, or discuss pore pressure prediction for a specific well or basin? Join our Telegram group for drilling engineering discussions, or visit our YouTube channel for step-by-step tutorials on pore pressure prediction and wellbore stability analysis.