Wellbore Stability Analysis - Geomechanical Modeling, Breakout Prediction, and Trajectory Optimization

Wellbore Stability Analysis - Geomechanical Modeling, Breakout Prediction, and Trajectory Optimization

Wellbore instability is the single largest source of non-productive time (NPT) in drilling operations worldwide, accounting for 10-15% of total well time on average and up to 30-40% in deep offshore and geologically complex wells. The instability mechanisms - shear failure producing breakout cavings, tensile failure producing lost circulation, and reactive formation swelling - are all governed by the same geomechanical framework: the in-situ stress state of the rock, the pore pressure, the rock strength, and the mud pressure applied at the wellbore wall. An engineer who understands this framework can predict where instability will occur before the well is drilled, select a well trajectory that minimizes instability risk, and design the mud weight window that keeps the wellbore in the safe zone between shear failure (too light) and tensile failure (too heavy). This guide provides the complete geomechanical workflow: stress tensor characterization, the Mohr-Coulomb failure criterion, breakout width calculation, and trajectory optimization for stability.

1. The In-Situ Stress State

1.1 The Three Principal Stresses and Their Measurement

The in-situ stress state at any subsurface point is described by three principal stresses that act perpendicular to each other. These stresses, combined with pore pressure, determine whether a given mud weight will stabilize or destabilize the wellbore:

Stress	Definition	Measurement Method	Typical Gradient (ppg)
Vertical stress (Sv)	Weight of overlying rock column. Also called overburden stress.	Integration of bulk density log from surface to depth: Sv = integral(rho_b x g x dz)	18-22 ppg (1.0 psi/ft average)
Minimum horizontal stress (Sh_min)	Smallest horizontal principal stress. The direction perpendicular to Sh_min is where drilling-induced tensile fractures propagate.	Leak-off test (LOT), extended LOT (XLOT), or minifrac test. Most reliable direct measurement.	13-17 ppg (varies widely)
Maximum horizontal stress (Sh_max)	Largest horizontal principal stress. Acts in the direction perpendicular to Sh_min. Controls where compressive failure (breakout) occurs.	Indirect: derived from breakout width on FMI log + Mohr-Coulomb model, or from borehole deformation analysis.	15-20 ppg (most uncertain)

1.2 Stress Regimes and Their Impact on Wellbore Stability

Three tectonic stress regimes (Anderson classification):

Normal faulting (Sv > Sh_max > Sh_min):
Typical: Sv = 19 ppg, Sh_max = 15 ppg, Sh_min = 13 ppg
Vertical wells most stable - Sv is the intermediate stress. Deviated wells drilling toward Sh_min direction become increasingly unstable.

Strike-slip faulting (Sh_max > Sv > Sh_min):
Typical: Sh_max = 20 ppg, Sv = 18 ppg, Sh_min = 14 ppg
Wellbore stability depends strongly on trajectory. Horizontal well drilled parallel to Sh_min most unstable.

Reverse faulting (Sh_max > Sh_min > Sv):
Typical: Sh_max = 22 ppg, Sh_min = 19 ppg, Sv = 17 ppg
Vertical wells very difficult to stabilize. Compressive stresses exceed overburden → high risk of breakout. Often requires mud weights approaching fracture gradient.

Example calculation - effective stresses at wellbore wall (vertical well):
At 9,000 ft TVD: Sv = 18.5 ppg, Sh_max = 17.2 ppg, Sh_min = 15.0 ppg, Pp = 12.5 ppg, MW = 13.0 ppg

Effective stresses (sigma_eff = sigma_total - Pp):
Sv_eff = (18.5 - 12.5) x 0.052 x 9,000 = 6.0 x 468 = 2,808 psi
Sh_max_eff = (17.2 - 12.5) x 468 = 2,200 psi
Sh_min_eff = (15.0 - 12.5) x 468 = 1,170 psi

Stress concentration at wellbore wall (Kirsch equations for vertical well):
At the point of maximum compressive stress (breakout location):
sigma_theta_max = 3 x Sh_max - Sh_min - Pp + dP_mud
= 3 x (17.2 x 0.052 x 9,000) - (15.0 x 0.052 x 9,000) - (12.5 x 0.052 x 9,000) + (13.0 x 0.052 x 9,000)
= 3 x 8,050 - 7,020 - 5,850 + 6,084
= 24,150 - 7,020 - 5,850 + 6,084 = 17,364 psi concentrated compressive stress at wellbore wall

2. Failure Criteria - When Does the Wellbore Fail?

2.1 Mohr-Coulomb Criterion for Shear Failure (Breakout)

Mohr-Coulomb failure criterion:
tau = C + sigma_n x tan(phi)

Or equivalently in terms of principal stresses:
sigma_1 = UCS + q x sigma_3

Where:
UCS = Unconfined Compressive Strength (psi)
q = (1 + sin(phi)) / (1 - sin(phi)) = tan^2(45° + phi/2) where phi = internal friction angle
sigma_1 = maximum principal stress at failure
sigma_3 = minimum principal stress (confining stress)

Breakout occurs when concentrated stress at wellbore wall exceeds UCS:
sigma_theta_max > UCS + q x (MW - Pp_gradient) x 0.052 x TVD

Example continued: sigma_theta_max = 17,364 psi
Rock properties from core: UCS = 6,500 psi, phi = 32° → q = (1 + sin32)/(1-sin32) = 1.530/0.470 = 3.255

Effective mud support at wellbore: sigma_3 = (MW - Pp) x 0.052 x TVD = (13.0 - 12.5) x 0.052 x 9,000 = 234 psi

Failure threshold = UCS + q x sigma_3 = 6,500 + 3.255 x 234 = 6,500 + 762 = 7,262 psi

17,364 psi concentrated stress >> 7,262 psi failure threshold → BREAKOUT WILL OCCUR at this mud weight

Required mud weight to prevent breakout:
Solve for MW when sigma_theta_max = failure threshold:
3 x Sh_max - Sh_min - Pp + MW x 0.052 x TVD = UCS + q x (MW - Pp) x 0.052 x TVD
(3 x 8,050 - 7,020 - 5,850) + MW x 468 = 6,500 + 3.255 x (MW x 468 - 5,850)
9,280 + 468 x MW = 6,500 + 1,524 x MW - 19,040
9,280 + 468 x MW = -12,540 + 1,524 x MW
21,820 = 1,056 x MW
MW = 21,820 / 1,056 = 20.66 ppg minimum MW to prevent breakout

But fracture gradient = 17.5 ppg → 20.66 ppg would fracture the formation before preventing breakout
This is a geomechanically unstable window. Some breakout is unavoidable. Design to minimize it rather than eliminate it.

2.2 Tensile Failure - Lost Circulation

Minimum stress at wellbore wall (breakout direction is at 90° to maximum stress concentration):
sigma_theta_min = 3 x Sh_min - Sh_max - Pp + MW_gradient x 0.052 x TVD

Tensile failure (lost circulation) occurs when:
sigma_theta_min < -T (tensile strength of rock, typically 200-500 psi for sedimentary rocks)

Or simplified (MW too high):
MW_max (ppg) = (3 x Sh_min - Sh_max + Pp + T) / (0.052 x TVD)

At 9,000 ft: Sh_min = 15.0 ppg, Sh_max = 17.2 ppg, Pp = 12.5 ppg, T = 350 psi = 350/(0.052 x 9,000) = 0.748 ppg:
MW_max = 3 x 15.0 - 17.2 + 12.5 + 0.748 = 45 - 17.2 + 13.248 = 41.048 ppg

Wait - dimensional inconsistency. Using consistent pressure units:
sigma_theta_min = 3 x Sh_min_psi - Sh_max_psi - Pp_psi + MW_psi
= 3 x 7,020 - 8,050 - 5,850 + MW_psi
= 21,060 - 8,050 - 5,850 + MW_psi = 7,160 + MW_psi

For tensile failure: 7,160 + MW_psi < -350 → MW_psi < -7,510 → negative = impossible

In this stress state, tensile failure cannot occur at this depth in a vertical well (minimum stress is always compressive). Fracture gradient from LOT (15.0 ppg) is the operative constraint, not tensile failure.

In normal and strike-slip faulting environments with lower horizontal stresses, tensile failure (sigma_theta_min < -T) can occur at moderate MW, causing drilling-induced fractures and lost circulation.

3. Breakout Width - Quantifying Acceptable Instability

3.1 Calculating Breakout Width from Mud Weight

When the mud weight cannot be raised to prevent all breakout (because it would exceed fracture gradient), the engineer must determine what breakout width is acceptable before hole integrity is compromised. Breakout angle wider than 90° typically results in an undrillable hole:

Breakout half-angle (degrees from the Sh_max direction):
cos(2 x Wbo/2) = (UCS + q x MW_eff - (Sh_max_eff + Sh_min_eff)(1+q)) / ((Sh_max_eff - Sh_min_eff)(1-q))

Simplified empirical approach using Zoback formulation:
At MW = 13.5 ppg (trying to stay below fracture gradient 17.5 ppg):
MW_eff = (13.5 - 12.5) x 0.052 x 9,000 = 468 psi
Sh_max_eff = 8,050 - 5,850 = 2,200 psi
Sh_min_eff = 7,020 - 5,850 = 1,170 psi

Failure threshold at MW=13.5: UCS + q x MW_eff = 6,500 + 3.255 x 468 = 6,500 + 1,523 = 8,023 psi
Concentrated stress at wall: 3 x Sh_max_eff - Sh_min_eff + MW_eff = 3 x 2,200 - 1,170 + 468 = 5,898 psi

5,898 < 8,023 → No failure at 13.5 ppg!

Hmm - this contradicts earlier calculation. Rechecking with absolute pressures vs effective:
This demonstrates why consistent units throughout is essential. Using effective stresses (sigma - Pp):
sigma_theta_eff = 3 x Sh_min_eff - Sh_max_eff - MW_eff (effective stresses at well wall in direction of Sh_max)
= 3 x 1,170 - 2,200 - 468 = 3,510 - 2,200 - 468 = 842 psi

vs UCS = 6,500 psi → 842 < 6,500 → no shear failure in this direction

In direction of Sh_min (breakout location):
sigma_theta_eff = 3 x Sh_max_eff - Sh_min_eff - MW_eff = 3 x 2,200 - 1,170 - 468 = 6,600 - 1,170 - 468 = 4,962 psi
vs UCS + q x sigma_3 (confining) = 6,500 + 3.255 x 0 = 6,500 psi (no confining in breakout direction)

4,962 < 6,500 → No breakout at 13.5 ppg mud weight in this particular stress state

Critical MW below which breakout begins:
sigma_theta_eff_breakout = 3 x Sh_max_eff - Sh_min_eff - MW_eff = UCS
3 x 2,200 - 1,170 - MW_eff = 6,500
5,430 - MW_eff = 6,500
MW_eff = -1,070 psi → negative → breakout impossible with any positive MW

Conclusion: In this stress state (Sh_max=17.2, Sh_min=15.0, Pp=12.5, UCS=6,500 psi), the rock is strong enough that no breakout occurs at any reasonable mud weight. The well can be drilled safely at MW = Pp + 0.5 ppg margin = 13.0 ppg.

4. Trajectory Optimization for Wellbore Stability

4.1 How Well Trajectory Affects the Stress Concentration

The Kirsch equations derived above apply to vertical wells. For deviated wells, the stress tensor must be rotated from the geographic frame (Sv, Sh_max, Sh_min) into the wellbore frame (axial, circumferential, radial) using coordinate rotation. The result is that some trajectories are significantly more stable than others in the same stress field:

Well Trajectory	Stress Concentration at Wall	Stability in Normal Faulting (Sv > Sh_max > Sh_min)	Required MW Change vs Vertical
Vertical	sigma_theta = 3Sh_max - Sh_min (standard Kirsch)	Most stable in normal faulting - Sv is intermediate stress	Baseline
Deviated 45° in Sh_max direction (azimuth along Sh_max)	Sv contributes - concentration increases as Sv component rotates into the stress matrix	Moderately less stable than vertical	+0.3 to +0.8 ppg
Horizontal in Sh_max direction (azimuth = Sh_max)	sigma_theta = 3Sv - Sh_min (Sv now acts as maximum stress at wall)	Significantly less stable - Sv (highest stress) concentrates at wall	+1.0 to +2.5 ppg
Horizontal in Sh_min direction (azimuth = Sh_min)	sigma_theta = 3Sv - Sh_max (better than Sh_max azimuth because Sh_max reduces concentration)	Most stable horizontal orientation in normal faulting	+0.5 to +1.5 ppg

4.2 Geomechanical Stability Map - Optimal Trajectory Selection

Minimum mud weight required as function of inclination and azimuth (ppg):
(Derived from full stress tensor rotation and Mohr-Coulomb criterion)

Normal faulting regime: Sv = 19 ppg, Sh_max = 16 ppg, Sh_min = 13.5 ppg, Pp = 11.0 ppg, UCS = 4,500 psi, phi = 30°

Vertical well (0° inclination): MW_min = 11.5 ppg
30° inclination, azimuth = Sh_min direction: MW_min = 11.8 ppg
30° inclination, azimuth = Sh_max direction: MW_min = 12.3 ppg
60° inclination, azimuth = Sh_min direction: MW_min = 12.5 ppg
60° inclination, azimuth = Sh_max direction: MW_min = 13.8 ppg
90° inclination (horizontal), azimuth = Sh_min direction: MW_min = 13.2 ppg
90° inclination (horizontal), azimuth = Sh_max direction: MW_min = 15.4 ppg

Fracture gradient = 14.5 ppg

Conclusion: Horizontal well drilled in the Sh_max azimuth requires MW = 15.4 ppg but fractures at 14.5 ppg → UNDRILLLABLE in this trajectory. Redirecting to Sh_min azimuth requires MW = 13.2 ppg → 1.3 ppg below fracture gradient → DRILLABLE.

This is why FMI imaging to determine Sh_max azimuth before designing horizontal well trajectory is not optional - it determines whether the well trajectory is physically possible.

Conclusion

The trajectory comparison in this article - horizontal well requiring 15.4 ppg MW in the Sh_max azimuth versus 13.2 ppg in the Sh_min azimuth, in a formation that fractures at 14.5 ppg - is the most direct demonstration of trajectory optimization value. The 15.4 ppg requirement is not a design challenge to be solved with better mud or casing design: it is a physical impossibility because fracturing occurs at 14.5 ppg. The only engineering solution is to change the well azimuth. A geomechanical stability map, generated before well planning, identifies this impossibility in pre-drill simulation rather than at 8,000 ft depth when the horizontal section is already committed to a specific azimuth and bit weight.

The breakout analysis that iterates to show no breakout in the specific stress state (Sh_max = 17.2 ppg, UCS = 6,500 psi) demonstrates that breakout is not automatic when horizontal stresses differ - rock strength also matters. A strong rock (high UCS) can remain stable under stress concentrations that would destroy weaker rock. This is why geomechanical wellbore stability analysis requires both stress measurement (from LOT and borehole imaging) and rock strength characterization (from core or scratch tests) - neither alone is sufficient to predict whether a given mud weight will stabilize or destabilize the wellbore.

Want to access our wellbore stability calculator with Kirsch stress equations, Mohr-Coulomb failure analysis, mud weight window, and trajectory stability map, or discuss geomechanical design for a specific well? Join our Telegram group for drilling geomechanics discussions, or visit our YouTube channel for step-by-step tutorials on wellbore stability analysis and trajectory optimization.