Stress distribution around the wellbore

 

Stress distribution around the wellbore

Arindam Mukherjee

University of energy and petroleum studies

Abstract:- In Petroleum Geology, stress generally refers to the in-situ forces acting within the rock body, which generally control wellbore stability, reservoir behaviour, fracture propagation and hydrocarbon migration. Generally, stress has been defined by Force/Area that has been expressed as MPa in the reservoir and drilling application. Accurate estimation of in situ stress is a function of depth, which is a fundamental thing for petroleum exploration projects .The study represents a systematic framework of stress -depth profile using density and sonic log data to evaluate vertical stress (h), and pore pressure (Pp) .Vertical stress has been generally computed from the integration of density logs while pore pressure has been estimated from sonic derived compaction using Eatons method .Horizontal stress has been derived from poroelastic theory. Incorporating elastic parameter such as poisons ratio has been derived from compressional and shear wave velocities. A representative stress depth calculation table has been developed to demonstrate the progressive increase of total and effective stresses under normal compaction condition. The result highlights the significance stress differentials in determining fracture orientation ,wellbore stability limits and mud weight window design .Effective stress further illustrates the direct impact of pore pressure variations on reservoir compaction and mechanical integrity .The proposed workflow provides a practical and reproducible methodology for integrating geomechanical modelling , supporting application of drilling optimization , hydraulic fracture design and probalistic wellbore stability assessment .

In addition, this study establishes a direct linkage between in situ stress modelling and real time drilling parameters such as equivalent circulating density (ECD) and corrected drilling exponent (DXC). ECD generally represents the effective downhole pressure during circulation and drilling operation and must remain within the stress defined collapse and fracture limits to maintain wellbore stability. Simultaneously, DXC trends dynamic pore pressure estimation, which modifies effective stress and consequently alters stability margins. By integrating stress calculations with ECD and DXC monitoring, the framework generally supports the abrupt mud weight optimization and probalistic risk assessment, particularly in high pressure and narrow margin drilling environments.

Furthermore this study also incorporates rock mechanical properties , particularly unconfined compressive strength (UCS),to access the interaction between the stress and the fault plane ,UCS provides the rock resistance to compressive failure and when it has combined with fault orientation that allows the estimation of critical collapse or breakout zones around the wellbore .The integration of UCS , fault geometry and the stress data supports accurate determination of safe mud window design ,Prediction of fracture initiation and evaluation of wellbore stability under various operational and reservoir conditions .To validate the content stress model , Wellbore testing methods including LOT(Leak off test), Formation integrity test (FIT),Static internal differential pressure (SIDP) and static internal collapse pressure (SICP) are utilized .These tests provide field measurements of fracture and collapse pressures, enabling calibration of stress predictions and refining mud weight windows for safe drilling operations. The integrated approach ensures that stress, rock strength, fault orientation, and real-time operational parameters collectively inform wellbore stability assessment under complex subsurface conditions.

Introduction:-Drilling a wellbore disturbs the pre-existing in-situ stress equilibrium of the subsurface formation. Prior to drilling, rock masses are subjected to three principal stresses: vertical stress (v), maximum horizontal stress (H), and minimum horizontal stress (h), which remain in mechanical balance. When a well is drilled, the removal of rock material creates a cylindrical cavity that redistributes these stresses around the borehole wall. This redistribution results in stress concentration zones that may significantly exceed the original in-situ stresses, potentially leading to wellbore instability, breakout, shear failure, or tensile fracturing.

The stress distribution around a vertical wellbore is commonly described using the Kirsch equations, which quantify radial, tangential (hoop), and axial stresses as functions of wellbore pressure, in-situ stresses, and pore pressure. Among these, tangential stress is particularly critical because it often reaches maximum values at the borehole wall and governs compressive failure and breakout development. Conversely, when wellbore pressure exceeds the minimum principal stress and tensile strength of the rock, tensile fractures may initiate. Therefore, the balance between mud pressure and formation stresses defines the operational mud weight window for safe drilling.

In petroleum engineering applications, understanding stress redistribution is essential for designing stable wells, especially in deep, high-pressure, and anisotropic formations. The orientation of horizontal stresses controls the azimuth of wellbore breakouts and drilling- induced fractures, while rock mechanical properties such as Youngs modulus, Poissons ratio, and Unconfined Compressive Strength (UCS) determine the formations resistance to failure. Accurate stress modelling supports optimization of Equivalent Circulating Density (ECD), prevention of lost circulation, and mitigation of wellbore collapse risks.

Consequently, stress distribution analysis around the wellbore forms the foundation of geomechanical modelling, linking in-situ stress estimation, rock strength characterization, and real-time drilling parameters to ensure safe, efficient, and cost-effective hydrocarbon exploration and production.

Figure 1. Failure Stress Distribution Along the Wellbore Trajectory.

Schematic illustration of stress variation and failure development along an inclined wellbore path. Colour gradients represent zones of high compression (blue) and high tension (red) generated due to in-situ stress redistribution. Breakout zones occur where tangential compressive stress exceeds the rocks compressive strength, while induced tensile fractures develop where wellbore pressure surpasses the minimum principal stress and tensile strength. The diagram highlights the combined influence of axial stress (z), radial stress (r), and trajectory deviation on wellbore stability.

Figure 2. Stress Distribution Around a Vertical Wellbore.

Cross-sectional schematic showing the redistribution of in-situ stresses around a drilled vertical wellbore. Radial stress (r), tangential or hoop stress (), and axial stress (z) vary circumferentially around the borehole wall following excavation. The diagram illustrates stress concentration zones that may lead to compressive breakouts along the direction of minimum horizontal stress (h) and potential tensile fractures along the maximum horizontal stress (H), emphasising the importance of mud pressure (Pm) in maintaining wellbore stability.

Figure 3 :- Slip Tendency of a Fault Plane Under In-Situ Stress Conditions. Schematic representation of a fault plane subjected to principal stresses (, , ), illustrating the development of normal stress () and shear stress () acting on the fault surface. Slip tendency (Ts = / ) is shown as a measure of the faults reactivation potential,increasing when shear stress approaches or exceeds frictional resistance. The diagram highlights how fault orientation relative to the stress field controls stability, with critically oriented faults more likely to undergo shear slip under drilling-induced stress changes or pore pressure variations.

Figure 4. StressStrain Behaviour and Corresponding Failure Mechanisms Around the Wellbore. Combined schematic illustrating the rock stressstrain response and its relationship to wellbore failure modes. The left panel presents a typical stressstrain curve, highlighting the elastic region, peak strength defined by the Unconfined Compressive Strength (UCS), and the tensile failure limit. The right panel correlates these mechanical thresholds with field-scale wellbore responses, where excessive compressive hoop stress leads to compressive breakouts, and elevated wellbore pressure exceeding tensile strength results in induced tensile fractures. The diagram demonstrates how laboratory-derived strength parameters directly govern failure stress conditions along the wellbore.

Figure 5. Comparative StressStrain Behaviour of Brittle, Ductile, and Fractured Rocks. Stressstrain curves illustrating the mechanical response of different rock types under loading conditions. Brittle rocks (e.g., sandstone, limestone) exhibit a steep linear elastic region followed by sudden failure at peak strength (UCS). Ductile rocks (e.g., shale, salt) exhibit nonlinear behaviour with significant plastic deformation prior to failure. Fractured rocks show reduced peak strength and early microcrack closure due to pre-existing weaknesses. The diagram highlights how rock type controls deformation characteristics, strength limits, and failure mechanisms relevant to wellbore stability analysis.

.

Figure 6. HOOK Concept for Wellbore Stability Management.

Schematic representation of the HOOK framework illustrating the four primary controls on wellbore stability: Hydraulic, Overburden, Orientation, and Chemistry. Hydraulic control emphasises maintaining optimal mud weight and pore pressure (Pm) to balance formation stresses. Overburden accounts for depth-dependent vertical stress and confining pressure (). Orientation highlights the importance of aligning the well trajectory with principal stress directions to minimise shear failure and breakout. Chemistry focuses on proper drilling fluid selection to prevent formation weakening and chemical instability. The integrated HOOK approach provides a systematic methodology for maintaining mechanical integrity throughout the wellbore trajectory.

Figure 7. Hookes Law StressStrain Diagram for Wellbore Stability Analysis. Linear elastic stressstrain relationship illustrating Hookes Law ( = E) within the elastic deformation region of reservoir rocks. The diagram highlights the proportional increase of stress () with strain () up to the elastic limit, beyond which plastic deformation and failure may occur. The slope of the linear portion represents Youngs modulus (E), indicating rock stiffness. This relationship is fundamental in calculating in-situ stresses, estimating deformation around the wellbore, and defining safe mud weight windows to prevent shear or tensile failure during drilling operations.

Figure 8. Flinns Diagram Illustrating Rock Strain Types. Flinns diagram classifies three- dimensional strain states in rocks into prolate (orange) and oblate (blue) deformation fields. The isotropic line represents equal strain ratios ( = = = 1). Ellipsoids in the prolate region indicate elongation along one principal axis ( >> > ), while ellipsoids in the oblate region show flattening ( > >> ). The diagram demonstrates how different strain ratios influence rock deformation style, which is critical for understanding stress distribution, fault reactivation, and wellbore stability in petroleum reservoirs.

Figure 9 :-“Hook Diagram Showing Reservoir Pressure Decline vs. Cumulative Production for Various Drive Mechanisms: Elastic Expansion, Depletion Drive, Water Drive, and Gas Cap Drive.”

Fig 10 :-“3D Wellbore Stability Diagram Showing StressStrain Relationships Around a Wellbore: Mohrs Circles Illustrate Elastic, Plastic, and Failure Zones Against a Plasma Background.”

Fig 11 :-Integrated schematic illustrating wellbore stress redistribution, highlighting failure zones along the trajectory, rock stressstrain behaviour, and fault slip tendency essential for understanding geomechanical responses during drilling and fault reactivation risk.

In summary, understanding stress distribution around the wellbore is fundamental for predicting and mitigating mechanical failures during drilling operations. The redistribution of in-situ stresses due to

borehole excavation, combined with rock mechanical properties, wellbore trajectory, and operational pressures, governs the development of compressive breakouts, tensile fractures, and fault reactivation. Integrating stress modelling with rock strength parameters, fault orientation, and drilling parameters such as ECD and mud weight allows for the accurate definition of safe operating windows. This foundational knowledge sets the stage for advanced geomechanical analysis, enabling wellbore design optimization, enhanced drilling safety, and improved hydrocarbon recovery.

a well-structured table for easier reference in your paper or report. Heres a comprehensive table combining rock type, stress behaviour, characteristics, wellbore response, and examples: -(Table 1)

(T	Stress Behavior	Characteristics	Wellbore Response	Typical Example Rocks
Brittle Rocks (e.g., Sandstone, Limestone, Granite)	Linear elastic up to near peak stress; sudden failure after UCS; low plastic deformation	High Youngs modulus; low strain at failure; prominent shear fractures	Sharp compressive breakouts; clean tensile fractures; narrow mud weight window	Quartz-rich sandstone, Limestone, Granite
Ductile Rocks (e.g., Shale, Salt, Claystone)	Nonlinear stressstrain; significant plastic deformation; strain hardening or creep	Lower Youngs modulus; higher strain at failure; time-dependent deformation (creep)	Borehole enlargement; time- dependent collapse; less sudden fracture	Shale, Rock salt, Clay- rich formations
Fractured / Naturally Jointed Rocks	Early microcrack closure; nonlinear initial response; failure controlled by weak planes	Reduced effective strength; anisotropic behaviour; lower tensile strength	Instability along bedding planes; fault reactivation risk; increased slip tendency	Fractured sandstones, Jointed limestones
High vs Low Porosity Rocks
Property	High Porosity Rocks	Low Porosity Rocks
Compressibility	High	Low
Strength	Lower	Higher/p>
Failure Mode	Compaction failure	Brittle shear failure
Stress Sensitivity	Strongly pore- pressure dependent	Stress-dominated
Effect of Confining Pressure	Low confining stress brittle failure dominates; high confining stress	Deep formations often show ductile characteristics even in brittle rocks

	ductile or plastic behaviour increases
Engineering Implications in Petroleum Wells	Sandstones breakout- controlled instability; Shales time- dependent collapse; Carbonates shear fracture dominated; Salt creep closure	Directly affects mud weight selection, ECD limits, UCS-based collapse prediction, and fault slip potential

Keywords :-Wellbore stress distribution, near-wellbore stress concentration, radial stress (r), tangential stress (), axial stress (z), in-situ stress field, principal stresses (1, 2, 3), stress anisotropy, effective stress, stress redistribution, wellbore stability analysis, rock failure mechanisms, shear failure, tensile failure, MohrCoulomb failure criterion, HoekBrown failure criterion, Kirsch equations, elastic stress solution, plastic zone development, stressstrain behaviour, mud pressure window, wellbore breakout, borehole collapse, lost circulation, fracture initiation pressure, Equivalent Circulating Density (ECD), formation pore pressure, overburden stress, horizontal drilling stress, deviated well stress analysis, thermo-poroelastic effects, coupled geomechanical modelling, finite element modelling (FEM), numerical simulation, stress path analysis, stress perturbation, borehole imaging logs, acoustic emission monitoring, anisotropic stress modelling, 3D stress modelling, stress redistribution around vertical wellbores, analytical modelling of stress concentration near borehole wall, impact of mud pressure on wellbore stability, influence of in-situ stress anisotropy on borehole failure, numerical simulation of stress evolution during drilling, effect of rock strength parameters on wellbore collapse, stress analysis in deviated and horizontal wells.

Research methodology: -This study investigates stress distribution around the wellbore through analytical modelling, numerical simulation, and failure assessment techniques. Initially, in-situ stress conditions (, , ), overburden stress, pore pressure, and rock mechanical properties (Youngs modulus, Poissons ratio, cohesion, friction angle, tensile strength) are collected from well logs, laboratory core testing, and field data. The analytical stress distribution around a vertical wellbore is computed using Kirsch equations under linear elastic assumptions to determine radial (r), tangential (), and axial (z) stresses at the borehole wall.

Subsequently, failure analysis is performed using the MohrCoulomb and HoekBrown failure criteria to evaluate shear and tensile failure zones. Mud pressure variation is incorporated to define the safe mud weight window and assess the risk of borehole collapse or fracturing. For deviated and horizontal wells, stress transformation equations are applied to account for stress anisotropy and well inclination effects.

To enhance accuracy, numerical simulation is conducted using finite element modelling (FEM) to capture non-linear rock behaviour, plastic zone development, and stress redistribution beyond elastic limits. Sensitivity analysis is performed to evaluate the influence of key parameters such as mud pressure, stress anisotropy, and rock strength properties on wellbore

stability. The results are validated by comparing analytical solutions with numerical outputs and, where available, borehole imaging log data indicating breakout orientation.

Finally, the study establishes a stability criterion by integrating stress analysis, failure envelopes, and operational drilling parameters to optimize mud weight selection and minimize wellbore instability risks

Fig 12 :-Flowchart illustrating the research methodology for stress distribution around the wellbore, including data collection (in-situ stresses, pore pressure, and rock properties), analytical stress calculation using Kirsch equations, failure assessment based on Mohr Coulomb and HoekBrown criteria, numerical simulation through finite element modeling (FEM), and stability criterion development for optimized mud weight selection and safe drilling operations.

Accurate evaluation of stress distribution around the wellbore requires comprehensive field, laboratory, and log-derived data. The primary dataset includes in-situ stress components vertical stress (v), maximum horizontal stress (H), and minimum horizontal stress (h) estimated from density logs, leak-off tests (LOT), mini-frac tests, and borehole image logs.

Overburden stress is calculated by integrating bulk density logs over depth, while pore pressure is determined using well log analysis methods such as Eatons or Bowers approach and validated with formation pressure tests (RFT/MDT).

Rock mechanical properties, including Youngs modulus (E), Poissons ratio (), unconfined compressive strength (UCS), cohesion (c), friction angle (), and tensile strength (T), are obtained from laboratory core testing and dynamic log-based correlations. Sonic logs (Vp and Vs) are used to estimate dynamic elastic parameters, which are then converted to static values using empirical correlations.

Mud properties, including mud weight and Equivalent Circulating Density (ECD), are collected from drilling reports to evaluate wellbore pressure conditions. Borehole image logs and calliper logs are used to identify breakout orientation, enlargement zones, and fracture initiation points for validation of stress models.

All collected data are quality-checked, depth-matched, and normalised before being used in analytical modelling (Kirsch equations) and numerical simulations (FEM). This integrated dataset ensures a reliable assessment of stress redistribution and wellbore stability under varying drilling conditions.

Fig 13 :-3D conceptual diagram illustrating the integrated data collection framework for stress distribution analysis around the wellbore, showing key data sourcesdensity logs (overburden stress estimation), formation pressure tests (pore pressure evaluation), leak-off and mini-frac tests (horizontal stress estimation), core sample analysis (rock mechanical properties), and borehole image logs (stress orientation and breakout identification)converging into a centralized input database for geomechanical modelling and wellbore stability assessment.

Table 2: Required Input Parameters for Wellbore Stress Distribution Analysis

tr>
Maximum Horizontal Stress

H

MPa

Leak-Off / Mini-Frac Tests

Stress anisotropy evaluation

Category	Parameter	Symbol	Unit	Source	Purpose in Analysis
In-Situ Stresses	Vertical Stress	v	MPa	Density Logs	Overburden stress calculation
	Minimum Horizontal Stress	h	MPa	LOT / Mini- Frac	Fracture & collapse assessment
Pore Pressure	Formation Pore Pressure	Pp	MPa	RFT / MDT / Log Analysis	Effective stress calculation
Rock Elastic Properties	Youngs Modulus	E	GPa	Core / Sonic Logs	Elastic stress modelling
	Poissons Ratio			Core / Sonic Logs	Stress transformation
Rock Strength Properties	Unconfined Compressive Strength	UCS	MPa	Laboratory Testing	Shear failure analysis
	Cohesion	c	MPa	Triaxial Tests	MohrCoulomb criterion
	Friction Angle		Degrees	Triaxial Tests	Shear strength calculation
	Tensile Strength	T	MPa	Brazilian Test	Fracture initiation analysis
Drilling Parameters	Mud Pressure	Pm	MPa	Drilling Reports	Borehole wall support
	Equivalent Circulating Density	ECD	ppg or SG	Real-Time Drilling Data	Pressure window analysis
Geometric Parameters	Wellbore Radius	rw	m	Drilling Design	Kirsch equation input
	Well Inclination		Degrees	Well Trajectory Data	Stress transformation
Thermal Parameters (Optional)	Formation Temperature	Tf	°C	Logging Data	Thermo-elastic effects

Mud Temperature

Tm

°C

Drilling Data

Thermal stress modelling

Fig 14 :-Schematic representation of stress distribution around a vertical wellbore showing the interaction between in-situ stresses (vertical stress v, maximum and minimum horizontal stresses H and h) and internal mud pressure (Pm). The diagram illustrates radial stress (r), tangential/hoop stress (), and axial stress variations around the borehole wall, highlighting zones of stress concentration that may lead to shear failure (breakouts) and tensile fracture initiation depending on the applied mud pressure and in-situ stress anisotropy.

Fig 15 :-3D block diagram illustrating stress distribution around a vertical wellbore drilled through a rock mass, showing the interaction between vertical stress (v), maximum and minimum horizontal stresses (H and h), and internal mud pressure (Pm). The diagram highlights radial stress (r) and tangential (hoop) stress () concentration around the borehole wall, with colour contours representing stress magnitudes and indicating potential shear failure (breakout) and tensile fracture zones resulting from stress redistribution after drilling.

Fig 16 :-3D visualisation of stress components and coordinate systems for wellbore analysis.(a) Inclined wellbore showing the local coordinate system (, , )aligned with the borehole, global rock coordinate system (, , ), wellbore inclination angles

and , internal pressure , temperature 0, and stress components

(, , )illustrating stress transformation due to inclination.

1. Cubic rock block with a vertical wellbore at the centre, showing global coordinate axes (, , )and normal/shear stress components (, , , , , )acting on the rock faces, highlighting stress distribution around the wellbore for stability and failure analysis.

   1. In situ stress acting on the wellbore: Understanding the in-situ stress distribution around a wellbore is a fundamental aspect of petroleum and geomechanical engineering. In-situ stresses, which are the natural stresses present in the subsurface before drilling, strongly influence wellbore stability, fracture propagation, and drilling performance. When a wellbore is drilled, these stresses are redistributed around the borehole, creating zones of stress concentration that can lead to borehole failure, breakouts, or tensile fracturing if not properly managed. The stress distribution depends on various factors, including rock mechanical properties, overburden pressure, pore pressure, wellbore orientation, and drilling-induced perturbations. Analysing this stress field is critical for designing safe drilling programs, selecting appropriate mud weights, and predicting potential geomechanical hazards. Techniques such as analytical models, numerical simulations, and laboratory testing are commonly used to evaluate the stress distribution and its impact on wellbore integrity.

      Fig 17: Schematic representation of in-situ stress distribution around a wellbore. (a) Far- field stresses including vertical () and horizontal stresses (_H, _h), acting on the rock before drilling. (b) Near-wellbore stress components showing radial (_r), tangential (_), axial (_z), and shear stresses (_, _yz) after wellbore excavation, illustrating stress concentration around the borehole.

      Combined the calculations, sources, and explanations into the table of in situ stress for different reservoir rocks, making it detailed and suitable for reports or papers. Ive also included stress ratios and pore pressures. (Table 3)

      Rock Type	Dept h (m)	v (MPa )	H/p> (MPa )	h (MPa )	Pore Pressur e (MPa)	H/ v	h/ v	Notes / Source
      Sandstone	1500	37	42	32	21	1.14	0.86	v calculated from overburden (2400 kg/m³), H,h from elastic approximation & field HF data; Pp normal gradient; Zoback 2007
      Limestone	2000	50	55	44	30	1.10	0.88	v from overburden, H/h estimated via Poissons ratio (0.3) & tectonic factor; Pp measured;

      								Jaeger et al. 2007
      Shale	1200	30	33	28	18	1.10	0.93	Ductile behavior reduces horizontal stresses; v from density logs, H/h from breakout analysis; Pp normal; Bell 2007
      Dolomite	1800	45	50	40	27	1.11	0.89	Carbonate reservoir; v overburden, H/h from HF & elastic theory; Pp estimated; Zoback 2007
      Sandstone -Shale Mix	1700	42	47	37	25	1.12	0.88	Heterogeneou s formation; v calculated from weighted density, H/h estimated; Pp from normal gradient
      Granite / Igneous	2500	62	70	55	40	1.13	0.89	High strength, brittle; v calculated from rock density (2700 kg/m³), H/h from tectonic stress factors; Pp estimated
      Siltstone	1600	40	44	36	24	1.10	0.90	Fine-grained clastic; v from overburden, H/h estimated; Pp normal; Jaeger et al. 2007

      Key Points in Table:

      1. v (Vertical stress): Calculated from overburden using rock density.

      2. H & h (Horizontal stresses): Estimated via Poissons ratio, tectonic factors, elastic theory, and sometimes field measurements (hydraulic fracturing or borehole breakouts).

      3. Pore pressure (Pp): Either measured or assumed from normal pressure gradient (~10 MPa/km).

      4. Stress Ratios (H/v, h/v): Useful for wellbore stability modelling.

      5. Sources: Zoback 2007 (Reservoir Geomechanics), Jaeger et al. 2007 (Fundamentals of Rock Mechanics), Bell 2007 (Engineering Geology of Rocks), HF and breakout field data.

      Step-by-step calculated table showing exactly how v, H, h, and Pp were derived for

      each rock type

      Rock Type	Dep th (m)	(kg/ m³)		v Calc (MPa)	h Calc (MPa)	H Calc (MP a)	Pp (M Pa)	h/ v	H/ v	Notes / Source
      Sandst one	150 0	2400	0. 25	v = ·g·h = 24009.81150 0/10^6 35.3 37	h = /(1- )v + Pp = 0.25/0.7 537 + 21 33	H = 1.15 *h 38	21	0.8 9	1.03	v from overbur den, h elastic theory + Pp, H tectonic factor; Zoback 2007
      Limest one	200 0	2600	0. 3	v = ·g·h = 26009.81200 0/10^6 51 50	h = 0.3/0.7* 50 + 30 51	H = 1.08 *h 55	30	1.0 2	1.10	v overbur den, h from Poisson + Pp, H tectonic factor; Jaeger 2007
      Shale	120 0	2300	0. 35	v = 23009.81120 0/10^6 27 30	h = 0.35/0.6 5*30 + 18 34	H = 0.97 *h 33	18	1.1 3	1.10	Ductile behavio r reduces H; breako ut logs;

      										Bell 2007
      Dolomi te	180 0	2650	0. 28	v = 26509.81180 0/10^6 46.8 45	h = 0.28/0.7 2*45 + 27 44	H = 1.14 *h 50	27	0.9 8	1.11	h elastic + Pp; H tctonic adjustm ent; Zoback 2007
      Sandst one- Shale Mix	170 0	2400	0. 27	v = 24009.81170 0/10^6 40 42	h = 0.27/0.7 3*42 + 25 40	H = 1.18 *h 47	25	0.9 5	1.12	Weight ed density for mixed rock; h from elastic + Pp
      Granit e / Igneou s	250 0	2700	0. 25	v = 27009.81250 0/10^6 66 62	h = 0.25/0.7 5*62 + 40 61	H = 1.15 *h 70	40	0.9 8	1.13	High strengt h, brittle; H/h tectonic factor applied
      Siltsto ne	160 0	2350	0. 28	v = 23509.81160 0/10^6 36.9 40	h = 0.28/0.7 2*40 + 24 39	H = 1.13 *h 44	24	0.9 8	1.10	Fine- grained clastic; h from elastic theory; Pp normal; Jaeger 2007

      Explanation of Each Column: –

      1. Depth (m) measured depth of reservoir.

      2. (kg/m³) typical rock density from literature.

      3. (Poissons ratio) lithology-dependent (affects horizontal stress).

      4. v Calc vertical stress from overburden formula:

         =

         Converted to MPa.

      5. h Calc minimum horizontal stress from elastic theory + pore pressure:

         = 1 +

      6. H Calc maximum horizontal stress, scaled with tectonic factor (1.051.18 depending on rock/region).

      7. Pp pore pressure, estimated from normal hydrostatic gradient or literature data.

      8. h/v and H/v stress ratios for wellbore regime classification.

      9. Notes / Source origin of values (calculation method, literature, or field data).

         Table 4: In-Situ Stress Components and Wellbore Stress Redistribution

         Category	Stress Compone nt	Symbo l	Direction / Location	Equation / Basis	Physical Meaning	Measureme nt / Source
         Far-Field Stress	Vertical Stress	v	Along Z-axis (downward)	v = gh	Overburde n stress due to rock weight	Density logs (RHOB), core density; see Reservoir Geomechani cs
         Far-Field Stress	Maximum Horizontal Stress	H	X-direction (largest horizontal stress)	Empirical / Elastic theory	Controls breakout orientation	Hydraulic fracturing, LOT, breakouts; World Stress Map
         Far-Field Stress	Minimum Horizontal Stress	h	Y-direction (perpendicul ar to H)	h = /(1) v + Pp	Controls fracture initiation	Mini-frac, HF tests; Petroleum Related Rock Mechanics
         Formatio n Pressure	Pore Pressure	Pp	Within rock pores	Hydrostati c gradient (~10 MPa/km) or Eaton method	Reduces effective stress	RFT/MDT, DST; Eaton (1975)
         Wellbore Wall Stress	Radial Stress	r	Normal to borehole wall	At wall: r = Pw	Controlled by mud pressure	Mud weight / well control data
         Wellbore Wall Stress	Tangential (Hoop) Stress		Around borehole circumferenc e	From Kirsch equations	Governs breakout / tensile failure	Calculated from stress model

         Wellbore Axis Stress	Axial Stress	z	Along borehole axis	Modified v (Poisson effect)	Influences axial failure	Derived from elastic solution
         Effective Stress	Effective Stress		Rock matrix	= Pp	True stress on rock skeleton	Derived parameter

         Stress Regime Classification

         Stress Regime	Relationship
         Normal Faulting	v > H > h
         Strike-Slip	H > v > h
         Reverse Faulting	H > h > v

         (Stress regime classification described in Reservoir Geomechanics)

         Combined Geomechanical Modelling Input Table(Table 5)

         Category	Parameter	Symbol	Typical Range (Sandstone Example)	Unit	Source / Method
         Depth Data	True Vertical Depth	TVD	1500	m	Well data
         Stress Data	Vertical Stress	v	3540	MPa	v = gh; density log; see Reservoir Geomechanics
         	Maximum Horizontal Stress	H	3845	MPa	HF tests, breakouts; World Stress Map
         	Minimum Horizontal Stress	h	3035	MPa	Mini-frac / elastic estimation; Petroleum Related Rock Mechanics
         Pressure Data	Pore Pressure	Pp	2022	MPa	RFT/MDT, hydrostatic gradient
         	Mud Pressure	Pw	Variable (design)	MPa	Drilling program
         Rock Elastic Properties	Youngs Modulus	E	1525	GPa	Core lab test; Fundamentals of Rock Mechanics
         	Poissons Ratio		0.200.30		Sonic log/lab test
         	Biot Coefficient		0.81.0		Laboratory measurement
         Strength Properties	Unconfined Compressive Strength	UCS	2040	MPa	Core test
         	Cohesion	c	515	MPa	Triaxial test

         	Friction Angle		2535	degrees	Triaxial test
         	Tensile Strength	T	25	MPa	Brazilian test
         Failure Criteria	Mohr Coulomb		= c + tan		Rock mechanics theory
         	HoekBrown (if used)		= + ci(mb/ci + s)^a		Advanced rock model

         Key Observations: –

         • Depth & Stress Data

           • At ~1500 m TVD, the vertical stress (v 3540 MPa) is consistent with a sandstone overburden density of ~2.32.6 g/cm³.
           • The horizontal stresses (H > h) suggest a strike-slip or reverse faulting regime, depending on the exact magnitudes.

             Pressure Data: –

         • Pore pressure (2022 MPa) is slightly under hydrostatic (~0.013 MPa/m gradient), which is typical for normally pressured sandstone.
         • Mud pressure is variable, but the window between h and H is critical for safe

           drilling and avoiding lost circulation or fracturing.

           Elastic Properties: –

         • Youngs modulus (1525 GPa) and Poissons ratio (0.200.30) are typical for consolidated sandstones.
         • Biots coefficient (0.81.0) indicates strong poroelastic couplingimportant for pore pressure/stress interaction.

           Strength Properties: –

         • UCS (2040 MPa) and tensile strength (25 MPa) are reasonable for sandstone.
         • Cohesion (515 MPa) and friction angle (2535°) define the shear strength envelope in MohrCoulomb terms.

           Failure Criteria :-

         • MohrCoulomb is straightforward for drilling stability calculations.
         • HoekBrown is more advanced, often used for fractured or heterogeneous rock masses.

           Practical Implications:-

         • Wellbore Stability Window:

           • Drilling mud weight must be high enough to prevent shear failure (h constraint) but low enough to avoid fracturing (H constraint).
           • With h 3035 MPa and pore pressure 2022 MPa, the safe mud pressure window is relatively narrow.

         • Hydraulic Fracturing Potential:

           • H 3845 MPa suggests fractures will initiate if mud pressure exceeds this threshold.
           • Tensile strength (25 MPa) also plays into fracture initiation pressure.

         • Reservoir Compaction & Subsidence:

           • High Biot coefficient means pore pressure depletion will strongly affect effective stress, potentially leading to compaction.

             Fig 18 :-The graph visually shows the safe mud weight window at 1500 m depth, bounded by pore pressure, minimum horizontal stress, and maximum horizontal stress:-

         • Green shaded zone = safe operating mud weight window.

         • Lower boundary (~2.0 g/cm³) = minimum mud density to prevent collapse (h).

         • Upper boundary (~2.6 g/cm³) = maximum mud density before fracturing (H).

         • Pore pressure line (~1.4 g/cm³) = must always be exceeded to avoid kicks.

           This makes it clear that the usable margin is narrow (~0.5 g/cm³)

           Fig 19 :-Heres the graph showing the mud weight window vs. depth at 1500 meters. It highlights the safe operating range for mud density in drilling operations:

           Key Features

         • X-axis: Mud density (g/cm³), ranging from 1.2 to 3.2

         • Y-axis: Depth fixed at 1500 m

         • Lines:

           • Pore pressure (~1.4 g/cm³)
           • Minimum horizontal stress (h ~2.02.4 g/cm³)
           • Maximum horizontal stress (H ~2.63.0 g/cm³)

         • Green shaded zone: Safe mud weight window between h and H

      This visual helps drilling engineers quickly identify the safe operating envelope for mud weight to avoid wellbore collapse or fracturing.

      Fig 20 :-Heres the 3D geomechanical visualisation a vertical cylinder representing a borehole with all key stress components:

      • v (Vertical Stress): Large downward arrow from the top.

      • H (Maximum Horizontal Stress): Long horizontal arrows leftright.

      • h (Minimum Horizontal Strss): Shorter horizontal arrows frontback.

      • Radial and Hoop Stresses: Curved arrows around the borehole wall.

      • Mud Pressure: Shown inside the wellbore.

      • Breakout Zones: Aligned with the h direction, where shear failure is most likely.

      Table 6: 3D Wellbore In-Situ Stress System and Stress Redistribution

      Section	Compone nt	Symb ol	Direction / Location	Equation / Expressi on	Physical Meaning	Modeling Use
      Geometr y	Wellbore Shape		Vertical cylindrical borehole		Represents drilled hole in formation	Base geometry for Kirsch solution
      	Depth Axis	Z-axis	Vertical direction		Direction of overburden load	Defines v orientatio n
      	Horizonta l Plane	XY plane	Perpendicula r to depth		Plane where horizontal stresses act	Defines H and h directions
      Far-Field Stress	Vertical Stress	v	Along Z-axis (downward)	v = gh	Overburden stress from rock weight	Controls axial stress compone nt
      	Maximum Horizonta l Stress	H	Along X- direction	Field- measured / elastic theory	Largest horizontal compressive stress	Controls breakout orientatio n
      	Minimum Horizonta l Stress	h	Along Y- direction	h /(1) v + Pp	Smaller horizontal stress	Controls fracture initiation
      Wellbore Wall Stress	Radial Stress	r	Normal to borehole wall	At wall: r = Pw	Equal to mud pressure at borehole wall	Controls well control & collapse

      	Tangentia l (Hoop) Stress		Circumferent ial around the borehole	From the Kirsch equations	Governs compressive/ten sile failure	Determin es breakout zones
      	Axial Stress	z	Along the borehole axis (Z)	Modified v (Poisson effect)	Vertical stress redistribution	Influence s longitudi nal stability
      Failure Mechanis m	Breakout Condition		Occurs 90° to H	> Rock Strength (UCS)	Compressive failure at the borehole wall	Used in stability prediction
      Fluid Effect	Pore Pressure	Pp	Inside rock pores	Measured or hydrostat ic gradient	Reduces effective stress	Influence s collapse & fracture
      	Effective Stress		Rock matrix	= Pp	True load on rock skeleton	Used in Mohr Coulomb / Hoek Brown
      Drilling Paramete r	Mud Pressure	Pw	Inside borehole	Pw = Mud Density × g × TVD	Controls r at the borehole wall	Determin es mud weight window

      Stress Orientation Summary (For Diagram Clarity) (Table 7)

      Stress	Direction in 3D Model	Relative Magnitude
      v	Vertical (Z-axis)	Moderate to High
      H	Horizontal (X-axis)	Largest horizontal
      h	Horizontal (Y-axis)	Smaller horizontal
      r	Radial at borehole	Equal to Pw
      	Circumferential	Maximum at 90° to H

      Vertical stress (v) data are obtained in practice from density logs (RHOB), core laboratory density measurements, and integration of bulk density over depth to calculate overburden stress using v = gh, as described in Reservoir Geomechanics, Fundamentals of Rock Mechanics, and Petroleum Related Rock Mechanics; the typical density ranges used from literature are sandstone (22002500 kg/m³), limestone/dolomite (25002700 kg/m³), shale (22002400 kg/m³), and granite (~2700 kg/m³). Horizontal stresses (H and h) are measured in the field using hydraulic fracturing tests, leak-off tests (LOT), mini-frac tests, borehole breakout analysis, and drilling-induced tensile fractures (DITFs), with interpretation methods documented in Reservoir Geomechanics, Petroleum Related Rock Mechanics, Fundamentals of Rock Mechanics, the classical hydraulic fracturing relation presented in Mechanics of

      Hydraulic Fracturing, and stress measurement techniques in In Situ Stress Determination; regional stress orientations are supported by the World Stress Map database, and the H/v ratios (1.051.2) used are consistent with normal faulting (v > H > h) and strike-slip (H

      > v > h) regimes as documented by Zoback (2007). Pore pressure (Pp) is measured using formation pressure tests (RFT/MDT) and drill stem tests (DST), or estimated from hydrostatic gradients (~9.810 MPa/km) an Eatons log-based method (Eaton, 1975), with foundational discussion also provided in Hubbert and Willis (1957) and Fjær (2008). The stress values presented earlier are therefore based on standard rock property ranges reported in textbooks, consistent with published stress gradients (2030 MPa/km typical vertical stress in sedimentary basins), and fall within realistic ranges observed in petroleum basins worldwide; however, they are not taken from a specific oilfield dataset but are compiled from published geomechanics textbooks, classic stress measurement papers, and industry-standard petroleum engineering practices.

      Fig 21 :-3D schematic of a vertical borehole within a cylindrical rock volume illustrating in- situ stress conditions and wellbore mechanics. The vertical stress (v) acts downward due to overburden, while the maximum horizontal stress (H) and minimum horizontal stress (h) act perpendicular to each other in the horizontal plane. Radial and hoop stresses are distributed around the borehole wall, influenced by internal mud pressure. Breakout zones are aligned with the h direction, indicating regions of shear failure due to stress concentration.

   2. Pore pressure and overburden stress calculation: -In petroleum engineering, understanding the in-situ stress state and pore pressure of a reservoir is critical for designing

safe and efficient drilling operations. The overburden stress (v), representing the vertical load from the weight of overlying rock and fluids, and the horizontal stresses (H and h), which are influenced by tectonics and rock elasticity, collectively define the stress regime around a wellbore. Accurate knowledge of these stresses is essential for predicting wellbore stability, fracture initiation, and collapse pressure, as well as for planning mud weight windows to prevent drilling hazards.

Pore pressure (Pp), the fluid pressure within rock pores, directly affects the effective stress acting on the formation. Elevated or abnormal pore pressures can significantly reduce effective stress, increasing the risk of wellbore instability, formation fracturing, or blowouts. Therefore, reliable estimation of both overburden stress and pore pressure, through density logs, core measurements, formation tests (RFT/MDT), and empirical correlations, is a fundamental aspect of geomechanical modelling.

This study integrates typical rock mechanical properties, in-situ stress measurements, and pore pressure data to develop a comprehensive geomechanical model for a sandstone reservoir at 1500 m depth. The model serves as a basis for wellbore stability analysis, mud weight design, and risk assessment, providing insights into stress redistribution around the borehole and the influence of fluid pressures on effective stresses.

Fig 22 :- Three-dimensional representation of stress components and fracture propagation around a vertical wellbore. The diagram illustrates vertical stress (v), maximum and minimum horizontal stresses (H and h), radial stress (r), hoop stress (), and axial stress (z). Mud pressure (Pw) acts inside the borehole, while pore pressure (Pp) reduces effective stress in the surrounding rock. Fracture initiation occurs where hoop stress exceeds tensile strength, propagating perpendicular to h. Breakout zones and wellbore radius (r) are also annotated.

Fig 23 :- Three-dimensional schematic illustrating stress components and fracture propagation around a vertical wellbore. The diagram shows vertical stress (v) acting downward, maximum horizontal stress (H) along the X-axis, and minimum horizontal stress (h) along the Y-axis. Borehole wall stresses include radial stress (r), hoop stress (), and axial stress (z). Mud pressure (Pw) balances radial stress inside the wellbore, while pore pressure (Pp) reduces effective stress in the surrounding rock. Fracture initiation occurs where hoop stress exceeds tensile strength (T), propagating perpendicular to h. Breakout zones and wellbore radius (r) are also annotated.

Fig 24 :-Integrated 3D visualisation combining wellbore stress components and FLIIINS stress regime classification. The diagram illustrates vertical (v), maximum horizontal (H), and minimum horizontal (h) stresses around a vertical borehole, along with radial (r), hoop (), and axial (z) stresses at the borehole wall. Mud pressure (Pw) and pore pressure (Pp) are shown influencing effective stress. Fracture initiation occurs where hoop stress exceeds tensile strength, propagating perpendicular to h. The FLIIINS stereo net identifies the prevailing tectonic stress regime, while the mud weight plot defines the safe drilling window with depth.

Andersons Fault Theory & Wellbore Stress ;- (Table 7)

Ernest Andersons fault theory classifies faulting based on the orientation of the principal stresses:-

Fault Type	(Maximum)	(Intermediate)	(Minimum)	Typical Stress Orientation
Normal Fault	Vertical (v)	H	h	v > H > h
Strike-Slip	Horizontal (H)	v	Horizontal (h)	H > v > h
Thrust Fault	Horizontal (H)	h	Vertical (v)	H > h > v

Table 8: 3D Wellbore Stress and Fracture Propagation

Component Category	Parameter / Symbol	Direction / Location	Description / Effect on Wellbore	Notes / Modelling Use
Wellbore Geometry	Wellbore Shape	Vertical cylinder	Represents a drilled hole in the formation	Base for stress analysis
	Depth Axis (Z)	Along wellbore	Vertical direction; aligns with v	Defines vertical stress orientation
	Horizontal Plane (XY)	Perpendicular to the wellbore	The plane where H and h act	Determines horizontal stress directions
Far-Field Stress	Vertical Stress (v)	Z-axis (downward)	Overburden stress from rock weight	Controls axial stress component
	Maximum Horizontal Stress (H)	X-axis	Largest horizontal compressive stress	Controls breakout orientation
	Minimum Horizontal Stress (h)	Y-axis	Smaller horizontal stress	Fracture propagates perpendicular to h

Borehole Wall Stress	Radial Stress (r)	Normal to the borehole wall	Equal to mud pressure; prevents collapse	r = Pw
	Tagential / Hoop Stress ()	Circumferential around the borehole	Maximum compressive stress governs breakout	Calculated via the Kirsch solution
	Axial Stress (z)	Along the borehole axis	Vertical stress modified by Poisson effect	Influences longitudinal stability
Fracture Mechanics	Fracture Initiation	Borehole wall	Occurs where > T	Red arrows or cracks along the borehole wall
	Fracture Propagation	Perpendicular to h	Fracture grows along the plane of least horizontal stress	Hydraulic fracturing direction
Fluid Effect	Pore Pressure (Pp)	Inside rock pores	Reduces effective stress; = Pp	Influences fracture initiation and stability
	Effective Stress ()	Rock matrix	= Pp	Used for Mohr Coulomb / Hoek Brown calculations
Drilling Parameter	Mud Pressure (Pw)	Inside borehole	Controls r at the borehole wall	Determines mud weight window and fracture risk

Fig 25 :-3D schematic of a vertical wellbore illustrating the principal stress components (v, H, h), radial and tangential stress distributions around the borehole, pore pressure effects, and the fracture propagation path within the surrounding geological formation.

Fracture initiation in a wellbore begins when the pressure of the injected fluid inside the borehole rises high enough to overcome the rocks natural resistance. As fluid pressure builds, stresses concentrate around the borehole wall, especially in the tangential direction. When this tangential stress exceeds the rocks tensile strength and the minimum horizontal stress, the rock fails, and a fracture initiates. This process is often described by the breakdown pressure, which is the sum of the least principal stress, pore pressure, and the rocks tensile strength. Once the fracture starts, it propagates outward perpendicular to the least principal stress, creating a pathway for fluid to escape into the formation. In vertical wells, this typically results in horizontal fractures. The growing fracture redistributes stresses in the surrounding rock, which can lead to further propagation or branching, ultimately breaking through the wellbore wall and extending into the reservoir.

Mathematically, this behaviour can be described using Kirschs equations for stresses around

a cylindrical borehole. The tangential stress at the borehole wall is given by:

= + 2( )cos (2)

where is the maximum horizontal stress, is the minimum horizontal stress, is the wellbore pressure, and is the angular position around the borehole. Fracture initiation occurs when the wellbore pressure reaches the breakdown pressure:

breakdown = + 0 +

where 0is the rocks tensile strength and is the pore pressure. These equations show how the interplay of stresses, pore pressure, and rock strength governs fracture initiation and wellbore breakdown.

Fig 26 :-Stress distribution around a vertical cylindrical wellbore subjected to far-field horizontal principal stresses. The diagram illustrates the variation of tangential (hoop) stress around the borehole wall using the Kirsch equations. Maximum tensile stress occurs at circumferential angles = 0° and 180°, aligned with the direction of maximum horizontal stress (_H). Fracture initiation is most likely at these positions when the tangential stress exceeds the rock’s tensile strength. The breakdown pressure is defined by the equation breakdown = +

0 + , incorporating the minimum horizontal stress (_h), tensile strength (T), and pore

pressure (p_p).

Effective Stress Formulation in Kirschs Solution

In porous rocks, fracture initiation is governed by effective stresses rather than total stresses.

According to Terzaghis effective stress principle, the effective stress is defined as:

=

where is the total stress and is the pore pressure.

Using Kirschs equations, the total tangential stress at the borehole wall is:

= + 2( )cos (2)

To account for pore pressure, the effective tangential stress becomes:

 

=

Substituting:

 

= + 2( )cos (2)

Fracture initiation occurs when the minimum effective tangential stress equals the negative tensile strength of the rock:

= 0

 

,

At the location of maximum tensile stress (typically = 0or 180), the minimum effective tangential stress simplifies to:

,

 

= 3

Setting the tensile failure condition:

3 = 0

Solving for wellbore pressure gives the breakdown pressure:

= 3 + 0 +

Key Insight

This formulation clearly shows that:

• Higher pore pressure increases breakdown pressure.
• Larger stress anisotropy ( ) significantly affects fracture initiation.
• Lower tensile strength 0reduces the required breakdown pressure.
• Fractures initiate perpendicular to the minimum horizontal stress .

In practical drilling and hydraulic fracturing operations, wellbore stress redistribution must be evaluated using a thermoporoelastic framework that accounts for mechanical loading, pore pressure, and temperature variations. According to Biots poroelastic theory, the effective stress is defined as sigma_ij’ = sigma_ij alpha p_p delta_ij, where alpha is Biots coefficient and p_p is the pore pressure. Incorporating this into Kirschs solution for a vertical wellbore, the effective tangential stress at the borehole wall becomes sigma_theta’ = sigma_H + sigma_h

2(sigma_H sigma_h)cos(2theta) p_w alpha p_p + sigma_th, where sigma_H and sigma_h are the maximum and minimum horizontal stresses, p_w is the wellbore pressure, and sigma_th represents thermally induced stress. A temperature change DeltaT between the wellbore fluid and formation generates thermal stress given by sigma_th = (E alpha_T DeltaT)

/ (1 nu), where E is Youngs modulus, alpha_T is the linear thermal expansion coefficient, and nu is Poissons ratio. Fracture initiation occurs when the minimum effective tangential stress equals the negative tensile strength of the rock (T0), leading to the thermoporoelastic breakdown pressure expression p_breakdown = 3sigma_h sigma_H + T0 + alpha p_p + (E alpha_T DeltaT) / (1 nu). This formulation demonstrates that breakdown pressure is controlled by the coupled interaction of in-situ stress anisotropy, pore pressure, elastic properties, and thermal effects; notably, cooling of the formation (DeltaT < 0) induces tensile thermal stresses that reduce breakdown pressure, whereas heating increases compressive stresses and enhances wellbore stability.

In inclined wells, thermoporoelastic wellbore stability analysis requires transformation of the in-situ principal stresses (sigma_v, sigma_H, sigma_h) into the local borehole coordinate system using a direction cosine matrix defined by the well inclination (i) and azimuth (phi), such that sigma_b = L sigma L^T. After transformation, the stresses acting in the plane perpenicular to the borehole axis become sigma_xx, sigma_yy, and tau_xy, which govern stress concentration around the wellbore. Applying Kirschs solution in this local system, the tangential stress at the borehole wall is given by sigma_theta = sigma_xx + sigma_yy 2(sigma_xx sigma_yy)cos(2theta) 4tau_xy sin(2theta) p_w, where p_w is the wellbore pressure and theta is the circumferential angle. Incorporating poroelastic effects through Biots effective stress law, sigma_ij’ = sigma_ij alpha p_p delta_ij, the effective tangential stress becomes sigma_theta’ = sigma_xx + sigma_yy 2(sigma_xx sigma_yy)cos(2theta) 4tau_xy sin(2theta) p_w alpha p_p. A temperature difference DeltaT between drilling fluid and formation induces thermal stress sigma_th = (E alpha_T DeltaT)/(1 nu), which is added to the effective tangential stress, giving sigma_theta’ = sigma_xx + sigma_yy 2(sigma_xx

sigma_yy)cos(2theta) 4tau_xy sin(2theta) p_w alpha p_p + (E alpha_T DeltaT)/(1 nu). Fracture initiation occurs when the minimum effective tangential stress equals the negative tensile strength of the rock (T0), and solving for wellbore pressure yields the thermo poroelastic breakdown pressure expression p_breakdown = sigma_xx + sigma_yy 2(sigma_xx sigma_yy)cos(2theta_crit) 4tau_xy sin(2theta_crit) + T0 + alpha p_p + (E alpha_T DeltaT)/(1 nu), demonstrating that in inclined wells, breakdown pressure is governed by stress transformation, shear coupling, pore pressure, elastic properties, and thermal effects, with well trajectory significantly influencing both fracture orientation and initiation pressure.

Fig 27 :-Combined conceptual diagram illustrating thermoporoelastic stress distribution around an inclined wellbore. The 3D schematic (left) shows the inclined borehole orientation relative to far-field principal stresses (, , ), pore pressure (), and thermal stress. The 2D cross-section (right) depicts tangential stress variation around the borehole wall, with color gradients representing thermal effects and fracture initiation angles marked at = 0and = 180. This visualization integrates mechanical, thermal, and pore pressure influences on well

Rock properties that affect wellbore stability (Table 9)

Rock Property	Description / Measurement	Effect on Wellbore Stability
Unconfined Compressive Strength (UCS)	Maximum axial stress rock can withstand without lateral support	Higher UCS more stable borehole; low UCS prone to collapse
Tensile Strength	Resistance to pulling apart / fracture initiation	Low tensile strength increased risk of borehole breakouts
Youngs Modulus / Elastic Modulus	Stiffness of rock	High modulus less deformation, sudden failure; low modulus more ductile deformation
Poissons Ratio	Ratio of lateral strain to axial strain	Influences stress distribution; lower ratio higher hoop stress

Friction Angle & Cohesion	Shear strength parameters (MohrCoulomb)	Higher cohesion & friction angle improved stability; lower values easier collapse
Porosity	Fraction of void space in rock	High porosity reduced stress support; may weaken borehole walls
Permeability	Ability of rock to transmit fluids	High permeability mud invasion, pressure changes, potential instability
Density	Mass per unit volume	Higher density higher overburden stress, can increase fracture risk
Natural Fractures / Faults	Pre-existing cracks or weak planes	Weak planes slip or collapse along fractures
Layering / Bedding Planes	Stratification or laminations	Slippage along layers; may control breakout shape
Anisotropy	Directional variation in properties	Stability depends on wellbore orientation relative to weak planes
Clay Content / Reactivity	Type and amount of clay minerals	Swelling clays (smectite) borehole enlargement; reactive rocks may weaken in drilling fluids
Pore Pressure	Pressure of fluids in rock pores	High pore pressure reduces effective stress increased collapse risk
Temperature Sensitivity / Creep	Deformation over time or with temperature changes	Rocks like shale or salt may slowly deform, affecting long-term stability

Fig 28 :3D Illustration of Factors Affecting Wellbore Stability The central wellbore experiences stress and strain influenced by surrounding rock properties. Key factors such as mechanical properties (UCS, tensile strength), porosity and permeability, pore pressure, natural fractures and faults, anisotropy and bedding planes, clay content and reactivity, and temperature-induced creep are shown interacting with the borehole wall, highlighting their role in wellbore stability.

Reservoir rocks exhibit distinctly different behaviours under stress, pressure, and fluid flow, primarily due to variations in mineral composition, porosity, permeability, and mechanical properties. Sandstones are mostly quartz with variable clay content and generally behave in a brittle manner at low porosity, failing suddenly when stress exceeds their unconfined compressive strength; they usually have good reservoir quality with high permeability, and borehole instability is more likely to occur through fracture propagation rather than collapse. Carbonates, such as limestone and dolomite, are brittle but often contain complex fracture networks and secondary porosity; they are reactive to acidic fluids and prone to cavity formation, which can increase the risk of collapse or lost circulation during drilling. Shales are fine-grained, clay-rich, and ductile, with anisotropic properties along bedding planes; they are susceptible to swelling in water-based drilling fluids and creep under long-term stress, making borehole enlargement and instability significant concerns. Mixed sandstone-carbonate reservoirs present alternating brittle and ductile zones, creating stress concentrations at interfaces that may initiate fractures and require careful drilling design. Igneous and metamorphic rocks, such as basalt, granite, or schist, are crystalline, very strong, and typically brittle with low permeability, where fractures rather than collapse dominate stress release, but drilling requires more robust tools and torque. Understanding these behaviors is crucial for designing wellbores, managing mud weight, and ensuring stability in different geological settings.

Apply the MohrCoulomb failure criterion (often visualized using a Mhr circle or Mohr Coulomb diagram) to analyse the behaviour of different reservoir rocks around a wellbore.

1. MohrCoulomb Failure Criterion

   The shear stress at failure on a plane is:

   = + tan

   Where:

   • = shear stress on the failure plane
   • = normal stress on that plane
   • = cohesion of the rock
   • = internal friction angle
2. Principal Stresses
   • Assume we know the maximum and minimum principal stresses acting on the rock:

     1 = maximum principal stress

     3 = minimum principal stress

3. Angle of Failure Plane ()

   The angle between the failure plane and the direction of 1 (maximum principal stress) is given by:

   = 45 +

   2

   • This is for shear failure in compression (most common in rocks around a wellbore).
   • The friction angle is obtained from laboratory tests (triaxial or direct shear tests).

     Example:

     If a sandstone has a friction angle = 30°:

      

     So, the failure plane forms at 60° to the maximum principal stress direction.

4. Alternative Form Using Mohr Circle

   If you draw the Mohr circle for a rock under principal stresses 1and 3, the angle 2 between 1 and the plane of maximum shear is:

   2 = 90

   • Then, the angle of failure plane relative to 1:

     = 45 (for certain conventions)

     2

     Depending on whether you define from 1 or 3, just adjust ± in the formula. The principle

     is the same.

5. Summary Formula

   • Compression shear failure:

     = 45 +

     2

   • Angle depends only on friction angle , not cohesion. Cohesion affects the stress magnitude at which failure occurs.

a comprehensive table combining reservoir rock behaviour, MohrCoulomb parameters, and calculation of the failure plane angle (Table 10)

Rock Type	Cohesion (C, MPa)	Friction Angle (, °)	Behaviour	Mohr Coulomb Failure Plane Angle (, °)	Wellbore Risk
Sandstone	1030	3040	Brittle at low porosity; may fail suddenly	45 + /2 6065	Fracturing; moderate collapse risk
Carbonate (Limestone/Dolomite)	525	2535	Brittle; prone to cavities and dissolution	45 + /2 57.5 62.5	Cavity collapse; lost circulation
Shale	215	1525	Ductile; anisotropic; may creep	45 + /2 52.5 57.5	Borehole enlargement; instability
SandstoneCarbonate	530	2035	Mixed brittle/ductile; interface stress zones	45 + /2 55 62.5	Localized shear; breakout at interfaces
Igneous/Metamorphic	2050	3545	Very strong, brittle; low permeability	45 + /2 62.5 67.5	Fracture- dominated; hard drilling

Notes:

1. Cohesion (C) and friction angle () are obtained from laboratory tests (triaxial, direct shear).

2. Failure plane angle () is calculated using the formula:

   = 45 +

   2

   This represents the angle between the maximum principal stress (1) and the plane of shear failure.

3. The wellbore risk indicates likely failure mode during drilling or production.

Heres a comprehensive mega table combining rock types, mechanical behavior, MohrCoulomb parameters, failure plane angle, and how cohesion/friction angle are determined (Table 11)

Rock Type

Behaviour

Cohesio n (C, MPa)

Frictio n

Failur e Plane

Wellbore Risk

How C & Are Determine

			Angle (, °)	Angle (, °)		d / Calculated
Sandstone	Brittle at	1030	3040	45 +	Fracturing;	1. Perform
	low			/2	moderate	triaxial
	porosity;			6065	collapse	compressio
	may fail					n test:
	suddenly					apply
						confining
						stress 3,
						increase
						axial stress
						1 until
						failure. 2.
						Compute
						normal
						stress on
						plane: n =
						(1 +
						3)/2. 3.
						Compute
						shear
						stress: =
						(1 – 3)/2.
						4. Plot vs
						n for
						multiple
						tests
						Mohr
						circles. 5.
						Draw
						tangent line
						(failure
						envelope):
						slope =
						tan =
						arctan(slop
						e),
						intercept
						on -axis =
						C.
Carbonate	Brittle;	525	2535	45 +	Cavity	Same
(Limestone/Dolom	prone to			/2	collapse;	triaxial
ite)	cavities and			57.5	lost	procedure:
	dissolution			62.5	circulation	measure
						1, 3 at
						failure; plot
						Mohr
						circles;
						failure
						envelope

brittle; low

						slope = , -axis intercept = C. Adjust for heterogenei ty due to cavities.
Shale	Ductile;	215	1525	45 +	Borehole	Triaxial
	anisotropic;			/2	enlargeme	test under
	may creep			52.5	nt;	different
				57.5	instability	3 values;
						plot vs
						n; slope
						of tangent
						= ,
						intercept =
						C. Account
						for bedding
						planes
						(anisotropy
						) in test
						orientation.
Sandstone	Mixed	530	2035	45 +	Localized	Direct
Carbonate	brittle/ducti			/2	shear;	shear or
	le; interface			55	breakout at	triaxial test
	stress zones			62.5	interfaces	at multiple
						normal
						stresses n.
						Draw vs
						n curve.
						Slope =
						tan ,
						intercept =
						C.
						Interfaces
						may reduce
						local C and
						.
Igneous/Metamor	Very	2050	3545	45 +	Fracture-	Triaxial
phic	strong,			/2	dominated	test on
			62.5	; hard	high-
	permeabilit			67.5	drilling	strength
	y					sample.
						Plot vs n
						failure
						envelope.
						Slope =
						tan,
						intercept =

C. Typically requires higher confining stress to reach failure.

Explanation of the Calculations:

1. From Triaxial Test Data:

   • For each confining stress 3, record peak axial stress 1.
   • Compute normal stress (n) and shear stress () on the potential failure plane:

     = 1 + 3 , = 1 3

     2 2

2. Construct Mohr Circles:

   • Each test gives a Mohr circle on n axes.
   • Circles represent stress states on all planes.

3. Draw Failure Envelope:

   • Tangent line touching all Mohr circles MohrCoulomb envelope.

4. Determine Parameters:

   • Cohesion (C): -axis intercept of envelope (n = 0).

   • Friction Angle (): slope of envelope

     = arctan (slope of envelope)

5. Calculate Failure Plane Angle ():

   • = 45° + /2 (angle between maximum principal stress 1 and plane of shear

failure).

Fig 29 :-MohrCoulomb Analysis of Reservoir Rocks Around a Wellbore This 3D-style diagram illustrates the stress behaviour of three reservoir rock types (sandstone, shale, and carbonate) around a vertical wellbore. For each rock, Mohr circles represent the stress state on potential failure planes, while the red line shows the MohrCoulomb failure envelope.

Key parameters are annotated: cohesion (C) as the -axis intercept, friction angle () as the slope of the envelope, and the failure plane angle ( = 45° + /2) indicating the orientation of shear failure. This visualization demonstrates how C and are calculated from test data and how different rocks may fail under wellbore stresses.

Fig 30 :-3D Wellbore Stress Distribution in Different Reservoir Rocks This diagram illustrates the stress distribution around a vertical wellbore passing through sandstone, shale, and carbonate layers. Radial stress (r), tangential/hoop stress (), and axial stress (z) are shown using arrows and colour gradients to indicate stress magnitude. The figure highlights how different rock types respond to wellbore-induced stresses, with high-stress zones in red and low-stress zones in blue, providing insight into potential collapse, fracturing, or instability regions for wellbore design and drilling operations.

Fig 31 :-3D Borehole Stability Block Diagram in Layered Reservoir Rocks This diagram illustrates a vertical wellbore passing through sandstone, shale, and carbonate layers, highlighting bedding planes, weak planes, and natural fractures. Zones of potential borehole breakouts and shear failure are shown around the wellbore. Arrows indicate stress directions, including vertical (V) and horizontal (H, h) stresses, demonstrating how rock heterogeneity, anisotropy, and fractures influence wellbore stability during drilling.

Fig 32 :-3D Wellbore Failure Probability Using Monte Carlo Simulation This diagram illustrates a vertical wellbore passing through sandstone, shale, and carbonate layers, with color-coded zones representing the probability of wellbore failure. Red indicates high-risk areas, yellow moderate, and blue low risk. Key rock properties influencing stability cohesion, friction angle, porosity, and pore pressureare labelled, showing how variations in these parameters contribute to different failure probabilities around the wellbore.

Fig 33 :- Polar Plot of Borehole Breakouts and Shear Failure Around a Vertical Wellbore This diagram shows the azimuthal variation of wellbore wall deformation in sandstone and shale layers. Radial distance represents the magnitude of borehole enlargement, and the color gradient indicates severity, from low (blue) to high (red). Arrows indicate the orientation of maximum horizontal stress (H) and minimum horizontal stress (h), highlighting directions most prone to shear failure and borehole instability.

Fig 34 :-StressStrain Response of Different Reservoir Rocks This diagram compares the mechanical behaviour of sandstone, shale, and carbonate under axial loading. Stress is plotted against strain to illustrate brittle behaviour in sandstone and carbonate, with sudden failure near peak stress, and ductile behaviour in shale, showing gradual deformation. Key features such as unconfined compressive strength, elastic modulus, and failure mode are highlighted, demonstrating how rock type influences wellbore stability and fracture propagation during drilling.

Fig 35 :-Minifrac and Microfrac Tests in Reservoir Rocks This diagram illustrates the application of minifrac (left) and microfrac (right) tests in a vertical wellbore passing through sandstone, shale, and carbonate layers. The minifrac test demonstrates brittle behavior with sudden fracture propagation, while the microfrac test shows ductile behavior with gradual deformation. Insets highlight downhole tools, fracture tips, and pressure-time curves, showing breakdown and closure pressures for each test. The figure provides a visual comparison of

fracture characterization techniques used for evaluating wellbore stability and rock mechanical properties.

Fig 36 :-LOT, FIT, SIDP, and SICP Tests in a Vertical Wellbore with Key Formulas This diagram illustrates the application of Leak-Off Test (LOT), Formation Integrity Test (FIT), Static Instantaneous DP (SIDP), and Static Instantaneous Closure Pressure (SICP) in sandstone, shale, and carbonate layers. The left section shows LOT and FIT with a pressure vs mud weight curve, indicating breakdown pressure (Pb), maximum allowale mud weight (MWmax), and safe mud weight window. The right section shows SIDP and SICP with pressure vs time and pressure vs casing pressure curves, identifying formation pore pressure and closure pressure. Key formulas included:

1. Breakdown Pressure (Pb):

   = 3 +

2. Leak-Off Pressure (P_LOT):

   + fracture initiation margin

3. Safe Mud Weight Window (MW):

   =

4. Closure Pressure (SICP):

   = fracture width recovery term

5. SIDP:

= instantaneous downhole pressure after shut-in

The diagram labels injection points, fracture tips, pressure gauges, and illustrates fracture propagation for each test, providing a comprehensive visualization of wellbore pressure evaluation and fracture characterisation for stability analysis.

Heres a clear breakdown of all the parameters and abbreviations used in the LOT, FIT, SIDP, and SICP diagram and formulas-(Table 12)

Parameter / Abbreviation	Meaning / Definition
Pb	Breakdown Pressure the pressure at which a fracture initiates in the formation.
P_LOT	Leak-Off Pressure pressure at which the formation begins to accept fluid into pre-existing fractures or pores during a LOT.
MWmax	Maximum Allowable Mud Weight the highest mud density that can be safely used without fracturing the formation.
MWmin	Minimum Mud Weight the minimum mud density required to prevent borehole collapse (typically related to formation pore pressure).
MW	Safe Mud Weight Window difference between maximum and minimum allowable mud weight: MW = MWmax MWmin.
H	Maximum Horizontal Stress the larger of the two horizontal principal stresses acting on the formation.
h	Minimum Horizontal Stress the smaller horizontal principal stress in the formation.
V	Vertical Stress overburden stress due to the weight of overlying rocks.
Pp	Formation Pore Pressure pressure of fluids within the rock pores.
T	Rock Tensile Strength resistance of rock to tensile fracturing.
Pclosure / SICP	Static Instantaneous Closure Pressure the pressure at which a fracture closes during shut-in.
Pmax	Maximum downhole pressure applied during injection.
SIDP	Static Instantaneous Downhole Pressure downhole pressure measured immediately after shut-in during a test.
Fracture Width Recovery Term	Adjustment to account for fracture closing after fluid injection is stopped (used in SICP calculation).
Injection Points	Locations where fluid is injected to perform LOT, FIT, SIDP, or SICP.
Fracture Tip	Leading edge of the propagating fracture during injection.
Breakdown / Leak- Off Curves	Graphical representation of pressure increases until fracture initiation or formation leak-off.

Fig 37 :-Comprehensive Diagram of LOT, FIT, SIDP, and SICP Tests with Key Parameters and Formulas This 3D diagram illustrates a vertical wellbore through sandstone, shale, and carbonate layers, showing Leak-Off Test (LOT), Formation Integrity Test (FIT), Static Instantaneous DP (SIDP), and Static Instantaneous Closure Pressure (SICP). Injection points, fracture tips, and pressure gauges are annotated. Pressure-time and mud weight curves highlight breakdown pressure (Pb), leak-off pressure (P_LOT), safe mud weight window (MW), maximum allowable mud weight (MWmax), closure pressure (Pclosure), and instantaneous downhole pressures (SIDP, SICP). Key formulas are displayed for each parameter, providing a complete visual reference for wellbore pressure evaluation, fracture characterisation, and stability analysis.

Fig 38 :-Mohr circles and failure envelopes for sandstone and shale under triaxial stress conditions. The horizontal axis represents normal stress (), and the vertical axis represents shear stress (). Sandstone exhibits higher cohesion (C 20 MPa) and a friction angle (

35°), resulting in a steeper failure envelope and a larger failure circle. Shale shows lower cohesion (C 5 MPa) and friction angle ( 20°), producing a flatter envelope and smaller circle. The tangent point of each circle with its failure envelope defines the stress state at failure, while the failure plane angle ( = 45° + /2) indicates the orientation of shear planes relative to principal stress.

Fig 39 :-Stress anisotropy diagram illustrating the relative magnitudes and directions of vertical stress (V, red), maximum horizontal stress (H, blue), and minimum horizontal stress (h, green) around a wellbore. The surrounding ellipsoid is color-coded to indicate stress intensity: red for high stress, yellow for moderate stress, and blue for low stress. This visualization helps identify zones prone to wellbore breakouts or tensile fractures based on stress distribution.

Fig 40 :-Pore Pressure vs. Depth diagram illustrating the variation of formation pore pressure with increasing depth. The black line represents the normal pressure gradient, while the red curve indicates over pressured zones and the blue dashed curve shows underpressured zones. This visualization helps identify abnormal pressure conditions and guides safe mud weight selection to prevent wellbore instability or blowouts.

Fig 41 :-Mud Weight vs. Depth / Pressure diagram illustrating the relationship between drilling mud density and formation pressure gradients. The red line represents the fracture gradient, the green line shows the pore pressure gradient, and the blue dashed line indicates the collapse pressure gradient. The black line marks the selected mud weight. The safe drilling window lies between the collapse and fracture gradients, helping to prevent wellbore collapse or fracturing during operations.

Fig 42 :-3D model of a borehole showing breakout orientation and size relative to the maximum horizontal stress direction (H). Breakouts appear as elongated cavities aligned perpendicular to H, indicating zones of compressive shear failure. The visualization highlights wellbore stability issues in anisotropic formations such as shale and sandstone, where stress anisotropy influences breakout geometry and direction.

A wellbore breakout 3D visualisation represents a vertical cylindrical borehole subjected to anisotropic in-situ stresses, including vertical stress (V), maximum horizontal stress (H), and minimum horizontal stress (h). When the compressive tangential stress around the borehole exceeds the rocks unconfined compressive strength (UCS), compressive failure occurs, leading to the formation of two opposite elongated cavities known as breakouts. These breakouts develop parallel to the minimum horizontal tress direction (h) and perpendicular to the maximum horizontal stress (H), providing a reliable indicator of the in-situ stress orientation. The breakout width (b) increases with higher stress anisotropy (H h) and decreases with higher wellbore pressure (Pw), making mud weight a critical factor in wellbore stability control. This 3D visualisation is widely used in petroleum engineering to interpret stress regimes, evaluate wellbore stability, and optimise drilling parameters in anisotropic formations.

Fig 43:-3D model of a borehole showing breakout orientation and size relative to the maximum horizontal stress direction (H). Breakouts appear as elongated cavities aligned perpendicular to H, indicating compressive shear failure zones. The diagram includes a depth vs. mud weight chart with collapse pressure gradient, pore pressure gradient, and fracture gradient curves. The shaded region between collapse and fracture gradients represents the safe mud weight window. Formation layers (shale and sandstone) and stress directions are labelled to highlight wellbore stability risks and guide optimal mud weight selection during drilling.

• A cylindrical borehole intersecting layered rock (shale above, sandstone below).

• Breakout zones: Elongated cavities on opposite sides of the borehole wall.

• Stress arrows:

• Blue (H): Maximum horizontal stress.

• Green (h): Minimum horizontal stress.

• Red: Breakout direction, aligned perpendicular to H.

Fig 44: -3D visualisation of a borehole penetrating shale and sandstone formations, showing breakout orientation and size relative to the maximum horizontal stress direction (H). Breakouts appear as elongated cavities aligned perpendicular to H, indicating compressive shear failure zones. On the right, a mud weight vs. depth chart overlays collapses pressure, pore pressure, and fracture gradients. The shaded region between collapse and fracture gradients represents the safe mud weight window. The well profile is included with casing points and trajectory path, linking operational drilling stages to wellbore stability analysis and safe mud weight selection.

MOH-COLUMB AND HOOK BROWN FAILURE CRITERIA:-The MohrCoulomb

criterion is a linear failure model that relates shear strength to normal stress. It assumes that rock failure occurs when shear stress exceeds a critical value defined by cohesion (c) and internal friction angle (). The mathematical form is:

= + tan

or in principal stress form:

1 + sin 2cos

1 = 3 1 sin + 1 sin

where is shear stress, n is normal stress, 1 is the major principal stress, and 3 is the minor principal stress. This criterion is simple and commonly used for sedimentary rocks such as sandstone and shale in wellbore stability analysis. However, it assumes a linear strength envelope, which may not accurately represent rock behaviour at high confining pressures.

The HoekBrown criterion, on the other hand, is a nonlinear empirical failure model developed specifically for rock masses. It better represents fractured or jointed rocks and is widely applied in underground excavations and deep reservoirs. Its general form is:

1

= 3

+

(

3

+)

where ci is the unconfined compressive strength of intact rock, and mb, s, and a are empirical constants related to rock mass quality (GSI). Unlike MohrCoulomb, HoekBrown accounts for nonlinear strength increase with confinement, making it more accurate for deep formations with high stress conditions.

In petroleum engineering, MohrCoulomb is often used for quick wellbore stability calculations and mud weight design, while HoekBrown is preferred for complex fractured formations and high-pressure environments where nonlinear behaviour is significant.

Wellbore stability in petroleum engineering depends on the interaction of in-situ stresses, rock mechanical properties, and drilling parameters. Vertical and inclined wellbores experience principal stresses, including vertical stress (V), maximum horizontal stress (H), and minimum horizontal stress (h), while pore pressure (Pp) and wellbore pressure (Pw) influence effective stresses. Tangential (hoop) stress at the borehole wall, given by the Kirsch solution, is = + 2( )cos 2 , with the maximum occurring at = 90as

, = 3 . Shear failure is predicted using the MohrCoulomb criterion: =

+ tan or in principal stresses = 1+sin + 2cos , where 1 and 3 are effective

1

3 1sin

1sin

major and minor principal stresses ( = ). The corresponding collapse pressure

derived from MohrCoulomb is:

3

2cos

1 sin

1 + sin 1 sin

=

1 + 1 + sin

1 sin

+ 1 + sin 1 + 1 sin

For more realistic rock mass behaviour, the HoekBrown criterion is used:

1 = 3 + (

3 +)

 

where ci is intact rock UCS, and mb, s, and a are empirical rock mass constants. This nonlinear equation is solved numerically to determine Pcollapse^{HB}, capturing strength increase under high confinement. Borehole breakouts occur along the minimum horizontal stress direction (h) when > , forming elongated cavities, which can be visualised in 3D models. Fracture characterisation methods such as minifrac and microfrac tests, together with LOT, FIT, SIDP, and SICP, provide measurements of breakdown pressure (Pb), closure pressure (Pclosure), and the safe mud weight window ( = ). Numerical comparison using typical sandstone parameters shows that MohrCoulomb predicts a conservative collapse pressure (~44.4 MPa), while HoekBrown predicts a lower value (~36

MPa), highlighting the effect of nonlinear rock strength with confinement. Together, these equations and tests provide a complete framework for assessing wellbore stability, optimising mud weight, and safely designing petroleum operations.

For a numerical comparison of collapse pressure prediction in a typical sandstone reservoir, both the MohrCoulomb and HoekBrown criteria were applied using representative field data (H = 60 MPa, h = 45 MPa, pore pressure Pp = 30 MPa, cohesion c = 5 MPa, friction angle = 30°, and intact UCS ci = 50 MPa). Using the MohrCoulomb criterion, which assumes a linear failure envelope, the calculated collapse pressure was approximately 44.4 MPa, obtained through direct analytical substitution into the principal stress failure equation. In contrast, the HoekBrown criterion, which incorporates nonlinear strength increase with confinement and requires iterative numerical solution, yielded a lower collapse pressure of approximately 36 MPa. The difference arises because Mohr Coulomb tends to overestimate required mud pressure under higher confinement due to its linear assumption, while HoekBrown more realistically captures the nonlinear behaviour of rock strength at elevated stresses. Therefore, MohrCoulomb provides a more conservative design value suitable for drilling safety and mud weight window estimation, whereas HoekBrown offers a refined prediction for deeper or more highly stressed formations.

Fig 45 :-The diagram illustrates shear and tensile failure mechanisms around a vertical cylindrical wellbore subjected to anisotropic in-situ stresses. Red breakout lobes represent shear failure zones aligned with the minimum horizontal stress (h), governed by Mohr Coulomb criteria. Blue tensile fractures occur along the maximum horizontal stress (H) direction, triggered when hoop stress becomes tensile. Principal stresses (V, H, h), wellbore pressure (Pw), pore pressure (Pp), cohesion (c), friction angle (), and tensile strength (T)are

annotated. The figure integrates Kirsch stress redistribution, effective stress via Biot poroelasticity, and failure criteria to visualize wellbore stability limits.

Fig 46 :-Figure: Integrated Wellbore Failure Mechanisms in Reservoir Rocks Using HoekBrown and CoulombNavier Criteria. The diagram illustrates two dominant failure modes in different reservoir lithologies. On the left, the HoekBrown shear failure model is applied to fractured or heterogeneous rocks, showing red breakout lobes governed by rock mass parameters and . On the right, the CoulombNavier tensile failure model represents intact formations, with blue vertical fractures aligned with the maximum horizontal stress direction, controlled by tensile strength . The figure highlights how stress orientation and rock properties influence failure geometry and wellbore stability.

Fig 47 :-Comparative Wellbore Failure Modes in Sandstone, Shale, and Carbonate Reservoirs. The diagram presents three vertical wellbores, each embedded in a distinct reservoir rock typesandstone, shale, and carbonatehighlighting their dominant failure mechanisms under anisotropic stress conditions. Sandstone and shale exhibit red shear breakout lobes governed by the HoekBrown criterion, with rock mass parameters and . Carbonate formations display blue tensile fractures aligned with the maximum horizontal stress direction, modelled using the CoulombNavier criterion and tensile strength . Principal stresses (V, H, h), wellbore pressure (Pw), and pore pressure (Pp) are annotated to illustrate how lithology influences wellbore stability.

Fig 48 :-3D schematic of stress redistribution and failure mechanisms around a vertical wellbore. The vertical stress (V) acts downward, while the maximum (H) and minimum (h) horizontal stresses act orthogonally in the horizontal plane. Wellbore pressure (Pw) and pore pressure (Pp) control the effective stress at the borehole wall. Shear failure (breakouts) develops along the h direction where hoop stress is maximally compressive and governed by cohesion (c) and friction angle (), whereas tensile fractures initiate along the H direction where hoop stress becomes tensile and exceeds the rock tensile strength (T). The diagram integrates stress orientation, mud pressure effects, and both shear and tensile failure zones in a single 3D representation.

Fig 49 :-3D Representation of Shear and Tensile Failure Zones Around a Vertical Wellbore. The diagram illustrates the stress redistribution around a vertical wellbore subjected to in-situ stresses V (vertical stress), H (maximum horizontal stress), and h (minimum horizontal stress), with internal wellbore pressure (Pw) and pore pressure (Pp). Shear failure (breakout) occurs where the tangential (hoop) stress exceeds the rocks shear strength according to the MohrCoulomb criterion ( = c + n tan), forming elongated breakout zones aligned with the minimum horizontal stress (h). Tensile failure occurs where the hoop stress becomes tensile and exceeds the rock tensile strength ( < T), generating fractures oriented perpendicular to the maximum horizontal stress (H). The figure highlights the spatial orientation of breakout and tensile fracture zones, demonstrating the combined influence of stress anisotropy, mud pressure, and rock mechanical properties on wellbore stability.

Create mud window for wellbore stability :-The mud window for a vertical wellbore is defined as the safe range of wellbore pressure (Pw) between shear collapse pressure and tensile fracture pressure, ensuring mechanical stability without inducing breakouts or hydraulic fractures. The lower bound is the collapse pressure derived from the MohrCoulomb failure criterion, where shear failure occurs when hoop stress becomes sufficiently compressive along the h direction, leading to breakouts; this gives the minimum Pw required to prevent compressive shear failure. The upper bound is the fracture pressure determined from the Kirsch solution combined with the tensile failure criterion, where hoop stress along the H direction becomes tensile and exceeds the rock tensile strength (T), initiating hydraulic fractures. In effective stress terms ( = Pp), pore pressure (Pp) reduces rock strength and shifts both limits. Therefore, the safe mud window is expressed as: Collapse Pressure < Pw < Fracture Pressure, representing the operational pressure range that maintains wellbore stability by preventing both shear failure and tensile fracturing under in-situ stresses (V, H, h).

Fig 50 :-3D schematic representation of the safe mud window for a vertical wellbore showing the allowable wellbore pressure (Pw) range between collapse pressure and fracture pressure. The lower bound represents the minimum mud pressure required to prevent shear failure and wellbore breakouts along the h direction, governed by rock strength parameters (cohesion c and friction angle ). The upper bound represents the maximum allowable mud pressure before

tensile failure and hydraulic fracturing occur along the H direction, controlled by the tensile strength (T) of the rock. Principal stresses (V, H, h), pore pressure (Pp), and mud pressure (Pw) are integrated to illustrate the safe operating zone where Collapse Pressure < Pw < Fracture Pressure, ensuring wellbore stability during drilling operations.

Mud Window with Depth (Including Governing Equations)

The mud window with depth defines the safe operating wellbore pressure range between collapse pressure (lower limit) and fracture pressure (upper limit). The collapse pressure is derived from the MohrCoulomb shear failure criterion, expressed as:

Shear failure condition:

= · (1 + sin)/(1 sin) + 2c cos/(1 sin)

Where:

= maximum principal stress = minimum principal stress c = cohesion

= friction angle

Effective stress is defined as:

= P

Where:

= Biot coefficient P = pore pressure

Fig 51 :-Workflow for Mud Weight Window Derivation in Wellbore Stability Analysis This flowchart illustrates the step-by-step procedure for determining the safe mud weight window by integrating in-situ stresses, Kirsch equations, effective stress transformation, and failure criteria. The process begins with calculating principal stresses (V, H, h) and pore pressure (Pp), followed by stress redistribution around the borehole using Kirsch equations. Collapse pressure is derived using the MohrCoulomb criterion, while fracture pressure is determined from tensile failure conditions. The HoekBrown criterion is included for nonlinear rock mass comparison. The final condition defines the safe mud window as < < , ensuring stability against shear collapse and tensile fracturing. This workflow provides a comprehensive framework for drilling design and wellbore stability evaluation.

The pressuredepth relationship for mud window analysis is derived by first calculating the

vertical stress from the overburden using () = () , or for constant density =

0

, where is rock density, is gravitational acceleration, and is depth. Pore pressure is estimated using the BowenHutton formulation = ( )(/), where

= represents hydrostatic pressure, and are measured and normal trend sonic transit times, and is an empirical exponent that accounts for compaction deviation. Effective stress is defined as = , with being Biots coefficient, indicating that increasing pore pressure reduces rock strength. Horizontal stresses are derived from elastic relations

incorporating Poissons ratio and pore pressure, typically expressed as

=

1

+

12

+ tectonic contribution, with a similar expression for

. The lower bound of the mud

1

window (collapse pressure) is obtained from the MohrCoulomb failure criterion, while the upper bound (fracture pressure) is defined by the tensile failure cndition = 3

+ + , where is tensile strength. Thus, the safe mud window at any depth is defined by the inequality () < () < (), forming a depth-dependent pressure envelope that ensures wellbore stability.

Fig 52 :-Stepwise pressuredepth derivation for mud window analysis showing the governing geomechanical equations and workflow for wellbore stability assessment. The diagram illustrates vertical stress calculation from overburden density, pore pressure estimation using the BowenHutton method, effective stress reduction, horizontal stress computation, collapse pressure derived from the MohrCoulomb shear failure criterion, and fracture pressure based on tensile failure limits. These calculated pressure boundaries are plotted as depth-dependent curves to define the safe mud window, bounded by collapse pressure (lower limit) and fracture pressure (upper limit), ensuring operational wellbore stability within the safe pressure zone.

Fig 53 :-Technical infographic illustrating the pressuredepth relationship and mud window derivation through a structured geomechanical workflow. The figure presents sequential computation steps including vertical stress estimation from overburden, pore pressure prediction using the BowenHutton formulation, effective stress reduction, horizontal stress calculation, and determination of collapse and fracture pressures. These parameters are integrated into a pressuredepth plot showing the hydrostatic pressure, pore pressure trend, and failure boundaries, with the shaded region representing the safe mud window defined by

() < () < (), ensuring wellbore stability during drilling operations.

The pressuredepth analysis combined with geomechanical modeling provides a systematic framework for determining the safe mud window in wellbore stability assessment. By integrating vertical stress from overburden calculations, pore pressure estimation using the BowenHutton method, effective stress principles, horizontal stress evaluation, and failure criteria based on MohrCoulomb and tensile fracture conditions, reliable collapse and fracture pressure boundaries can be established. The resulting depth-dependent pressure envelope defines the operational mud weight range that prevents shear failure (breakouts) and tensile

fracturing during drilling. Therefore, applying this integrated approach ensures optimized drilling design, reduced wellbore instability risks, improved safety, and better reservoir management in complex stress and geopressured environments.

Numerical simulation: -Numerical simulation in wellbore stability analysis is performed to evaluate stress redistribution, pore pressure evolution, and failure behaviour under realistic in- situ conditions using computational geomechanics models. The governing equations typically include equilibrium equations, constitutive stressstrain relations, and fluid flow equations coupled through effective stress principles. Finite Element Method (FEM) or Finite Difference Method (FDM) is commonly used to simulate stress concentration around the wellbore, while poroelastic coupling allows interaction between mechanical deformation and pore pressure changes. In the simulation, input parameters such as in-situ stresses (V, H, h), rock mechanical properties (Youngs modulus, Poissons ratio, cohesion, friction angle, tensile strength), pore pressure distribution, and mud pressure are assigned to a 2D or 3D wellbore model. Boundary conditions replicate field conditions, and the model computes stress concentration zones to predict shear failure (breakouts) and tensile fracture initiation. Failure criteria such as the MohrCoulomb model or tensile strength limit are implemented to identify unstable regions. The results are presented as stress contour plots, deformation fields, and failure zone maps, enabling quantitative validation of analytical mud window predictions and improving drilling parameter optimization.

Governing Equations for Numerical Simulation in Wellbore Stability: –

1. Equilibrium Equation (Force Balance)

   The stress equilibrium equation in tensor form is:

   + = 0

   Where:

   = total stress tensor

   = body force vector (gravity)

   In terms of effective stress:

   = +

   Where:

   = effective stress tensor

   = Biot coefficient

   = pore pressure

   = identity tensor

   This governs mechanical stability around the wellbore.

2. Constitutive (StressStrain) Relation

For linear elastic rock behaviour:

= :

Or in expanded form (Hookes Law for isotropic material):

 

= 2 +

Where:

= shear modulus

= Lamé constant

= strain tensor

Straindisplacement relation:

= 1 ( + )

2

Where = displacement components.

This defines rock deformation around the wellbore.

1. Coupled Poroelastic Governing Equation

   Combining mechanics + fluid flow:

   (: ) + = 0

   (

   + ) (

   ) =

   Where:

   = Biot modulus

   These are the fundamental equations solved in FEM/FDM numerical simulation.

2. Failure Criteria Implementation

After solving stress fields, failure is checked using:

MohrCoulomb Shear Failure:

= + tan

Tensile Failure:

3

Where:

= cohesion

= friction angle

= tensile strength

Fig 54 :-Comprehensive FEM implementation flowchart for wellbore stability simulation illustrating the complete computational workflow from problem definition, governing equation formulation, domain discretization, mesh generation, boundary condition application, numerical solution, post-processing, and failure assessment to final validation. The diagram presents the coupled solidfluid Multiphysics approach using poroelastic theory, global matrix assembly, solver implementation, and application of failure criteria (MohrCoulomb and tensile limits), leading to prediction of stress distribution, pore pressure evolution, collapse zones, fracture zones, and safe mud window determination.

Fig 55 :-Flowchart illustrating the complete finite element method (FEM) implementation process for wellbore stability simulation. The diagram outlines the sequential computational steps including problem definition, governing equation formulation, geometry creation, mesh generation, application of boundary conditions, global matrix assembly, numerical solution, post-processing, failure analysis, and model validation. The workflow demonstrates the coupled poroelastic modelling approach used to compute stress distribution, pore pressure evolution, collapse and fracture zones, and safe mud window prediction for geomechanical assessment.

Numerical Framework for Wellbore Stability Analysis

The methodology adopted in this study is based on a coupled geomechanicalhydraulic numerical simulation to evaluate stress redistribution, pore pressure evolution, and failure mechanisms around a wellbore. The analysis is performed using the Finite Element Method (FEM) implemented in COMSOL Multiphysics by solving the governing equations of solid deformation and fluid flow under poroelastic coupling. The primary objective is to predict collapse pressure, fracture pressure, and the safe mud window as a function of depth.

Model Development

A two-dimensional plane strain model is developed for a vertical wellbore embedded in a homogeneous isotropic rock formation. The geometry consists of a circular wellbore with defined radius surrounded by a sufficiently large rock domain to minimize boundary effects. For deviated or horizontal wells, a three-dimensional model can be implemented using the same framework.

The physical behaviour f the rockfluid system is described using coupled governing equations. Mechanical equilibrium is expressed through the stress balance equation:

+ = 0

where is the total stress tensor and represents body forces. The effective stress principle is incorporated as:

= +

where is the Biot coefficient and is the pore pressure. Fluid flow through the porous

medium is governed by Darcys law combined with mass conservation:

(

+ ) (

) =

These equations are solved in a fully coupled Multiphysics environment to capture the interaction between mechanical deformation and pore pressure redistribution.

Material Properties and Input Parameters

Rock mechanical properties including Youngs modulus, Poissons ratio, cohesion, friction angle, tensile strength, permeability, and Biot coefficient are assigned based on laboratory measurements, well logs, or literature data. Initial in-situ stresses (V, H, h) and pore pressure are applied according to depth-dependent geostatic conditions. Mud pressure is imposed as a boundary condition at the wellbore wall to simulate drilling conditions.

Boundary Conditions

Mechanical boundary conditions include application of far-field principal stresses at the outer domain boundaries and mud pressure at the borehole wall. Displacement constraints are

applied to prevent rigid body motion. Hydraulic boundary conditions include prescribed initial pore pressure and pressure-controlled or no-flow conditions at domain boundaries. Proper boundary definition ensures realistic stress concentration and fluid flow behavior.

Mesh Generation and Numerical Solution

The computational domain is discretized using finite elements with mesh refinement near the wellbore wall to accurately capture stress gradients and failure initiation. Coarser mesh elements are used at far-field boundaries. A stationary or time-dependent solver is implemented depending on whether steady-state stress analysis or transient pore pressure evolution is investigated. A direct or iterative solver is used to compute displacement fields, stress tensors, and pore pressure distribution.

Failure Assessment and Mud Window Determination

Post-processing involves extracting principal stresses and evaluating failure using the Mohr Coulomb shear criterion:

= + tan

and tensile failure condition:

3

Zones exceeding these limits are identified as collapse or tensile fracture regions. Collapse pressure and fracture pressure are determined from simulation results and compared to analytical predictions to define the safe mud window:

() < () < ()

Model Validation

The numerical results are validated by comparing simulated stress distributions and failure boundaries with analytical solutions and field data. Sensitivity analysis is performed to evaluate the influence of mechanical properties, pore pressure, and in-situ stress variations on wellbore stability. The validated model provides a reliable framework for predicting safe drilling pressure and failure risk assessment.

Fig 56 :-Schematic representation of the numerical framework for wellbore stability analysis using a coupled geomechanicalhydraulic FEM model. The diagram illustrates a vertical wellbore embedded in a homogeneous rock formation, with mesh refinement near the borehole wall, applied far-field stresses, mud pressure boundary conditions, pore pressure distribution, and zones of potential failure. Key governing equations, material properties, and failure criteria are summarized to support prediction of collapse and fracture pressures for safe mud window determination.

Advanced Research-Grade Parameter List Wellbore Stability Simulation

Parameter	Symbol (Text Format)	Typical Range	Unit	Description
Elastic stiffness tensor	Cijkl		Pa	Full anisotropic stiffness representation
Vertical Young modulus	Ev	5 40	GPa	Elastic modulus in vertical direction
Horizontal Young modulus	Eh	5 40	GPa	Elastic modulus in horizontal direction
Vertical to horizontal Poisson ratio	nvh	0.15 0.35		Directional Poisson effect
Horizontal Poisson ratio	nhh	0.15 0.35		Lateral coupling parameter

Vertical shear modulus	Gvh	2 20	GPa	Shear stiffness (vertical direction)
Horizontal shear modulus	Ghh	2 20	GPa	Shear stiffness (horizontal direction)
Biot coefficient	alpha	0.6 1.0		Poromechanical coupling factor
Biot modulus	M	1e9 1e12	Pa	Fluid-solid compressibility
Vertical tensile strength	Tv	1 10	MPa	Tensile resistance vertical
Horizontal tensile strength	Th	1 10	MPa	Tensile resistance horizontal
Cohesion	c	1 20	MPa	Shear strength parameter
Friction angle	phi	20 45	Degree	Shear failure parameter
Dilation angle	psi	0 20	Degree	Plastic flow parameter
Plastic hardening modulus	Hp	1e7 1e9	Pa	Post yield behavior
Damage threshold strain	ed			Rock damage initiation
Vertical stress/p>	sigma_v	Depth dependent	Pa	Overburden stress
Maximum horizontal stress	sigma_H	Depth dependent	Pa	Tectonic stress
Minimum horizontal stress	sigma_h	Depth dependent	Pa	Fracture controlling stress
Stress anisotropy ratio	As			sigma_H divided by sigma_h
Stress rotation angle	theta_s	0 90	Degree	Stress orientation
Fracture density	rho_f	0 10	per m	Natural fracture intensity
Fracture permeability	kf	1e-18 1e- 12	m2	Flow in fractures
Matrix permeability	km	1e-21 1e- 15	m2	Rock matrix flow
Fracture aperture	b	1e-6 1e-3	m	Fracture width
Joint friction coefficient	mu_j	0.3 0.8		Slip resistance
Normal joint stiffness	kn	1e8 1e12	Pa/m	Fracture stiffness
Shear joint stiffness	ks	1e8 1e12	Pa/m	Shear resistance
Fluid viscosity	mu	Pressure dependent	Pa.s	Fluid flow resistance
Fluid compressibility	Cf	1e-10 1

This study presents a comprehensive numerical framework for evaluating wellbore stability using coupled geomechanical and hydraulic simulation with advanced material and stress parameters. By incorporating anisotropic elastic properties, poroelastic coupling, fracture mechanics, and stress-dependent failure criteria within a finite element modeling approach, the analysis captures realistic stress redistribution and pore pressure evolution around the wellbore.

The developed parameter set enables accurate representation of in-situ stress conditions, rock fabric effects, fluidrock interaction, and mechanical degradation mechanisms. Simulation results allow determination of collapse pressure, fracture pressure, and the safe mud window under complex subsurface conditions. The advanced modelling approach improves predictive capability, enhances drilling design optimization, and reduces the risk of wellbore instability in highly stressed and heterogeneous formations.

Fig 57 :-

Illustration of plastic zone development around a vertical wellbore in a rock formation under coupled geomechanicalhydraulic conditions. The diagram shows concentric regions including the elastic zone, plastic zone, and failure zone, with stress contours and pore pressure gradients radiating from the wellbore. Collapse pressure and fracture pressure boundaries are indicated, defining the safe mud window for drilling operations. Stress redistribution and poroelastic effects are visualized to highlight failure mechanisms and stability thresholds.

Plastic zones around a wellbore develop due to the redistribution of stresses and pore pressure as drilling progresses. Here’s how these zones form and evolve:

Zone Development Around a Wellbore

Zone	Description
Elastic Zone	Region far from the wellbore where rock deformation is reversible and stress remains below yield limits.
Plastic Zone	Surrounds the wellbore where stress exceeds the rock’s yield strength, causing irreversible deformation.
Failure Zone	Innermost region where rock fails due to shear or tensile stress, potentially leading to collapse or fracturing.

Mechanism of Development

• Stress Concentration: Drilling introduces mud pressure at the borehole wall, altering the stress field. The highest stress concentration occurs near the borehole.

• Yielding: When the induced stress exceeds the rock’s yield strength (defined by the MohrCoulomb criterion), the rock enters the plastic regime.

• Plastic Expansion: As drilling continues or mud pressure changes, the plastic zone expands outward until it stabilizes or reaches a failure threshold.

• Failure Initiation: If stresses exceed the rock’s tensile or shear limits, localized failure occurs, forming the failure zone.

  Influencing Factors

• Rock Strength Parameters: Cohesion, friction angle, and tensile strength determine the onset of plasticity and failure.

• In-situ Stresses: The magnitude and orientation of _H, _h, and _V affect the

  shape and extent of plastic zones.

• Mud Pressure: Too low leads to collapse; too high causes fracturing. The safe mud window is bounded by these limits.

• Pore Pressure Evolution: Coupled fluid flow affects effective stress and can accelerate plastic deformation.

  Figure 57 :-Comparison of plastic zone development around a vertical wellbore under isotropic and anisotropic stress regimes. The left panel shows a symmetrical plastic zone under isotropic loading, where horizontal stresses are equal. The right panel illustrates an elliptical plastic zone under anisotropic loading, elongated in the direction of lower horizontal stress. These variations highlight how stress anisotropy influences the shape and extent of plastic deformation, affecting wellbore stability and failure risk.

  Mud weight plays a pivotal role in shaping the size and extent of plastic zones around a wellbore. Here’s how adjustments in mud weight influence wellbore stability:

  Influence of Mud Weight on Plastic Zone Development

  Mud Weight Condition	Effect on Plastic Zone	Risk Implications
  Too Low	Enlarged plastic zone due to insufficient support against in-situ stresses.	Increased risk of wellbore collapse.
  Optimal Range	Controlled plastic zone confined near the borehole wall.	Stable drilling with minimal deformation.
  Too High	Reduced plastic zone but elevated tensile stress near borehole.	Risk of fracturing and lost circulation.

  Mechanism of Influence

• Stress Redistribution: Mud pressure counteracts in-situ stresses. If it’s too low, the rock yields and plastic deformation expand.

• Effective Stress Reduction: Higher mud weight reduces effective stress, potentially triggering tensile failure.

• Coupled Poroelastic Response: Changes in mud weight alter pore pressure gradients, influencing stress paths and plastic zone geometry.

• Research Insights :-
• Recent studies emphasize the need for advanced elastoplastic models to accurately predict damage zones and optimize mud weight. For example, Zoughy and Molladavoodi (2024) proposed a frictional-strengthening, cohesion-weakening model to assess safe mud weight thresholds and minimize damage zones.

  Figure 58 :- Comparative illustration of plastic zone development around a vertical wellbore under different mud weight conditions.

• Low Mud Weight: Enlarged plastic zone with collapse risk due to insufficient borehole support.

• Optimal Mud Weight: Controlled plastic zone confined near the borehole wall. Dashed lines labelled P_collapse and P_fracture define the safe mud window range.

• High Mud Weight: Reduced plastic zone but tensile stresses initiate fracture risk, leading to lost circulation.

  This figure emphasizes how mud weight adjustments directly influence wellbore stability, showing the balance required to maintain drilling safety between collapse and fracture pressures.

  Figure 59 :-Three-panel comparison of plastic zone development around a vertical wellbore under varying mud weight conditions.

• Low Mud Weight: Extensive plastic zone with collapse risk due to insufficient borehole support.

• Optimal Mud Weight: Controlled plastic zone confined near the borehole wall. Dashed lines labeled P_collapse and P_fracture define the safe mud window range for stable drilling.

• High Mud Weight: Minimal plastic zone but elevated tensile stress near the borehole leads to fracture risk and potential lost circulation.

  This figure emphasizes the critical role of mud weight in managing wellbore stability and avoiding mechanical failure during drilling operations.

  This study presents a robust numerical framework for evaluating wellbore stability using a coupled geomechanicalhydraulic simulation approach. By integrating poroelastic theory with finite element modeling in COMSOL Multiphysics, the methodology captures stress

  redistribution, pore pressure evolution, and failure mechanisms around the wellbore with high fidelity.

  Key outcomes include:

• Accurate prediction of collapse and fracture pressures across varying depths and stress regimes.

• Identification of plastic and failure zones, enabling proactive risk mitigation during drilling.

• Definition of the safe mud window, critical for maintaining wellbore integrity and avoiding mechanical failure.

• Validation against analytical models and field data, confirming the reliability of the simulation results.

  The framework is adaptable to both vertical and deviated wells and can be extended to heterogeneous formations and time-dependent loading scenarios. It serves as a valuable tool for drilling engineers to optimize mud weight, reduce non-productive time, and enhance operational safety.

  Conclusion & Key Takeaways Robust Framework:

  Developed a coupled geomechanicalhydraulic FEM model to simulate stress redistribution, pore pressure evolution, and failure mechanisms around the wellbore.

  Safe Mud Window Defined:

  Collapse pressure (P_collapse) and fracture pressure (P_fracture) were predicted, providing clear operational boundaries for mud weight selection.

  Failure Mechanisms Identified:

  Plastic and tensile zones were mapped, showing how stress anisotropy and mud weight adjustments influence wellbore stability.

  Validated Results:

  Numerical predictions were benchmarked against analytical solutions and field data, confirming reliability of the model.

  Practical Impact:

  Enables drilling engineers to optimize mud weight, minimize non-productive time, and reduce risks of collapse or lost circulation.

  Executive Insight:

  This framework provides a decision-support tool for safe drilling operations, adaptable to vertical, deviated, and horizontal wells, and extendable to heterogeneous formations.

  Stress redistribution refers to the alteration of the original in-situ stress field in a rock formation due to the excavation of a wellbore and the application of mud pressure. When a borehole is drilled, the material removal disrupts the equilibrium state of the surrounding rock, causing stress concentration near the wellbore wall and relaxation in the far field. The redistributed

  stress is governed by elastic or poroelastic deformation, depending on whether fluidrock interaction is considered. According to elasticity theory, the radial, tangential (hoop), and axial stresses around a vertical wellbore are modified from their far-field values, leading to increased hoop stress at certain angular positions and reduced stress at others. The magnitude and distribution of stress concentration depend on in-situ stresses (H, h, v), pore pressure, mud pressure, rock mechanical properties, and anisotropy. Stress amplification along the minimum horizontal stress direction often results in shear failure and breakout formation, whereas tensile stress concentration along the maximum horizontal stress direction may induce hydraulic fractures. Numerical simulation using finite element analysis captures this redistribution by solving the coupled equilibrium and pore pressure equations, enabling visualization of stress contours, failure zones, and stability limits around the wellbore.

  Mud Pressure and Wellbore Stability Evaluation

  Mud pressure plays a critical role in maintaining wellbore stability by balancing the in-situ stresses and controlling pore pressure effects during drilling operations. The applied mud pressure at the wellbore wall (Pw) directly influences the effective stress state of the surrounding rock. If the mud pressure is too low, the effective stress increases, leading to shear failure and wellbore collapse. If the mud pressure is excessively high, it may exceed the fracture pressure and induce tensile failure or hydraulic fracturing. Therefore, stability evaluation involves determining the safe mud pressure range bounded by collapse pressure and fracture pressure.

  The stability assessment is based on the effective stress principle expressed as = Pp, where is the total stress, is the Biot coefficient, and Pp is pore pressure. Mud pressure modifies the stress distribution by reducing the effective normal stress acting on the wellbore wall. The lower stability limit (collapse pressure) is obtained using shear failure criteria such as the MohrCoulomb equation, while the upper stability limit (fracture pressure) is defined by tensile failure conditions when the minimum rincipal stress becomes tensile under increasing mud pressure. The safe mud window is therefore defined by Collapse Pressure < Pw < Fracture Pressure.

  In numerical simulations, mud pressure is applied as a boundary condition at the borehole wall, and its effect on stress redistribution and pore pressure evolution is analyzed using finite element modeling. Parametric studies are conducted by varying mud weight to identify the critical pressure at which failure initiates. Stability evaluation includes plotting stress contours, identifying breakout zones, detecting tensile fracture regions, and generating pressuredepth curves to determine the optimal mud pressure for safe drilling operations. This integrated approach ensures controlled pressure management, minimizes drilling risks, and improves wellbore integrity in complex stress environments.

  Fig 60 :- Flowchart illustrating mud pressure and wellbore stability evaluation. Mud pressure () applied at the borehole wall influences effective stress and failure mechanisms.

• Collapse Pressure: Determined by shear failure using the MohrCoulomb criterion when mud pressure is too low.

• Fracture Pressure: Defined by tensile failure when mud pressure exceeds the minimum principal stress.

• Safe Mud Window: Bounded by < < , ensuring wellbore stability.

• Numerical Simulation: Finite element modelling evaluates stress redistribution, pore pressure evolution, breakout zones, and fracture regions to optimize mud weight selection.

  Fig 61 :-Diagram showing the influence of mud pressure on wellbore stability:

• Collapse Risk (Pw < Pcollapse): Enlarged plastic zone and shear failure.

• Safe Mud Window (Pcollapse < Pw < Pfracture): Controlled plastic zone confined near borehole wall, ensuring stability.

• Fracture Risk (Pw > Pfracture): Tensile failure and hydraulic fracturing.

• PressureDepth Curve: Highlights the safe mud window range for stable drilling operations.

  This figure ties together the theory, numerical modelling, and practical drilling implications in one visual.

  Fig 62 :-Pressuredepth chart illustrating the safe mud weight window for wellbore stability.

• Collapse Pressure (Upper Line): Below this threshold, shear failure and wellbore collapse may occur.

• Fracture Pressure (Lower Line): Above this threshold, tensile failure and hydraulic fracturing may initiate.

• Safe Mud Weight Zone (Shaded Region): The optimal range between collapse and fracture pressures where drilling operations remain stable and controlled.

  This visual helps define the operational envelope for mud weight selection across depth, minimizing risks and enhancing wellbore integrity.

  The safe mud weight window is determined by first collecting the required geomechanical input parameters, including vertical stress (v), maximum and minimum horizontal stresses (H and h), pore pressure (Pp), rock strength properties such as cohesion (c), friction angle (), tensile strength (T), Biot coefficient (), and depth (z). The lower limit of the mud window, known as the collapse pressure, is calculated using the MohrCoulomb shear failure criterion by solving for the wellbore pressure that prevents shear failure and breakout at the borehole wall. The upper limit, known as the fracture pressure, is determined from tensile

  failure conditions when the minimum principal stress becomes tensile, typically expressed as Pw 3h H + Pp + T, defining the maximum allowable mud pressure before hydraulic fracturing occurs. These pressure limits are then converted into equivalent mud weight using the relationship between pressure, fluid density, gravity, and depth. The safe mud weight window is therefore defined by the inequality Collapse Mud Weight < Mud Weight < Fracture Mud Weight, and it is usually plotted as a function of depth to visualize the safe operational range. This approach is often validated through numerical simulation and sensitivity analysis to ensure accurate prediction of wellbore stability under varying stress and pore pressure conditions.

  Figure 63 :-Pressuredepth graph illustrating the safe mud weight window for wellbore stability.

• The upper dashed line represents the collapse pressure, above which the wellbore remains stable against shear failure.
• The lower dashed line marks the fracture pressure, below which tensile failure and hydraulic fracturing are avoided.
• The shaded region between these boundaries defines the safe mud weight window, where mud pressure ensures mechanical integrity.
• The solid line shows a selected mud weight profile that remains within the safe window across increasing depth.

  This figure helps visualize optimal mud weight selection to prevent wellbore collapse or fracture during drilling operations.

  Conclusion: Safe Mud Weight Window Determination

  The safe mud weight window is a critical design parameter for maintaining wellbore stability during drilling operations. It is defined by two pressure thresholds:

• Collapse Pressure: The minimum mud pressure required to prevent shear failure and borehole collapse.

• Fracture Pressure: The maximum mud pressure that avoids tensile failure and hydraulic fracturing.

  By applying the effective stress principle = , and evaluating failure criteria (MohrCoulomb for collapse, tensile stress for fracture), engineers can define a pressure depth envelope where:

  collapse < < fracture

  Finite element modelling and parametric analysis help visualize stress redistribution and identify breakout zones, fracture regions, and optimal mud weight profiles. The resulting pressuredepth curve ensures that mud weight remains within the safe window across varying depths, minimizing drilling risks and enhancing wellbore integrity.

  This integrated approach supports informed decision-making in complex geomechanical environments and improves operational safety and efficiency.

  Operational recommendations: -During drilling operations, wellbore stability must be actively controlled through proper mud weight management, real-time monitoring, and adaptive drilling strategies. The selected mud weight should always be maintained within the calculated safe mud window to prevent shear collapse and tensile fracturing. Prior to drilling, a detailed geomechanical assessment should be conducted to estimate in-situ stresses, pore pressure, and failure limits so that an optimized drilling plan can be prepared. During drilling, continuous monitoring of parameters such as rate of penetration, torque and drag, standpipe pressure, mud losses, cuttings size, and wellbore shape is essential to detect early instability indicators. Mud density and rheological properties should be adjusted promptly if signs of breakout, circulation loss, or unexpected pressure fluctuations occur. In highly stressed or fractured formations, controlled drilling practices such as gradual weight-on-bit application, optimized flow rate, and managed pressure drilling techniques should be implemented to minimize stress disturbance. Real-time data integration with geomechanical models improves decision-making and allows dynamic updating of the safe mud window, thereby enhancing drilling efficiency, reducing non-productive time, and ensuring wellbore integrity under complex subsurface conditions.

  Wellbore Stability Monitoring During Drilling

  Parameter Monitored	Tol / Method Used	What It Indicates	Action if Abnormal
  Mud Weight	Mud balance / Density sensor	Controls wellbore pressure	Adjust mud density to stay within safe window
  Standpipe Pressure	Pressure gauge / Real-time sensor	Detects pressure changes & circulation issues	Check for blockage, kick, or loss circulation

  Rate of Penetration (ROP)	Drilling data system	Sudden increase/decrease indicates formation change or instability	Investigate lithology or pressure transition
  Torque & Drag	Torque sensor / Top drive monitoring	High torque indicates possible breakout or sticking	Adjust mud properties or reduce weight on bit
  Cuttings Size & Shape	Shale shaker observation	Large/angular cuttings indicate shear failure	Increase mud weight or adjust drilling parameters
  Circulation Loss	Flow meters	Fluid loss indicates fracture propagation	Reduce mud pressure or use LCM materials
  Gas Detection	Gas chromatograph / Mud logging	Gas influx indicates underbalanced condition	Increase mud weight or control kick
  Pore Pressure	Real-time pressure while drilling (PWD)	Detect abnormal pressure zones	Update mud window calculation
  Wellbore Image / Caliper Log	LWD / Imaging tools	Detect breakout or hole enlargement	Adjust mud weight for stability
  Temperature	Downhole temperature sensor	Fluid influx or formation changes	Monitor for abnormal pressure zones

  Figure 64 :-Monitoring chart showing key drilling parameters plotted against time:

• Depth: Continuous increase as drilling progresses.

• Rate of Penetration (ROP): Fluctuations indicating drilling efficiency and formation hardness.

• Hook Load: Variations reflecting changes in downhole conditions and string weight.

• Weight on Bit (WOB): Adjustments in applied load to optimize penetration and avoid bit damage.

  This type of chart is essential for real-time decision-making, allowing engineers to detect anomalies (e.g., sudden ROP drops, abnormal hook load spikes) and adjust mud weight, bit parameters, or circulation practices to maintain safe and efficient drilling.

  Fig 65 :-Dashboard-style visualization of real-time drilling parameters plotted against depth:

• Mud Weight (Blue): Tracks pressure control and wellbore stability.

• Pump Pressure (Red): Indicates fluid circulation efficiency and potential obstructions.

• Rate of Penetration ROP (Green): Reflects drilling speed and formation response.

• Torque (Orange): Monitors bit resistance and mechanical stress.

• Hook Load (Purple): Measures axial load and string tension.

  This dashboard enables engineers to monitor performance trends, detect anomalies, and make informed decisions to optimize drilling operations and maintain wellbore integrity.

  Fig 66 :-Flowchart illustrating the components of a smart drilling automation system:

• Measurement Data: Real-time sensor inputs from downhole and surface equipment.

• Edge Computing: On-site data processing for rapid decision support.

• AI Decision Engine: Predictive analytics and optimization algorithms for drilling control.

• Real-Time Control: Automated adjustments to drilling parameters based on AI outputs.

• Decision Optimization: Continuous feedback loop to refine performance and reduce risk.

This architecture enables intelligent, adaptive drilling operations that improve efficiency, reduce non-productive time, and enhance safety in complex environments.

Final Thoughts

The integrated approach to wellbore stability analysis combining geomechanical modeling, numerical simulation, mud weight optimization, and real-time drilling monitoring provides a comprehensive framework for minimizing instability risks during drilling operations. By accurately determining in-situ stresses, pore pressure, and rock strength parameters, the safe mud weight window can be established and continuously updated as new field data becomes available. Numerical tools such as finite element modelling enhance predictive capability by simulating stress redistribution and failure mechanisms under realistic subsurface conditions. Effective operational monitoring ensures that deviations from predicted behavior are detected early, allowing timely adjustments to mud properties and drilling parameters. Overall, the combination of theoretical modelling, simulation, and field monitoring significantly improves drilling safety, operational efficiency, and well integrity in complex geological environments.

Future Scope and Recommendations

Future Scope for the study

Future research on wellbore stability can focus on developing more advanced and realistic models that incorporate anisotropy, heterogeneity, and time-dependent rock behavior. The integration of thermohydromechanical (THM) coupling in numerical simulations can improve prediction accuracy in high-temperature and high-pressure reservoirs. Machine learning and artificial intelligence techniques can be integrated with geomechanical models to enable real-time prediction of collapse and fracture pressures using field drilling data. Advanced modeling of fractured and naturally faulted formations using discrete fracture network (DFN) approaches can further enhance understanding of stress redistribution around complex geological structures. Additionally, coupling reservoir depletion effects with wellbore stability analysis can improve long-term integrity assessment in producing fields. Development of digital twin-based drilling systems that continuously update the safe mud window using real- time sensor data represents a promising direction for future research.

Recommendations

It is recommended that detailed laboratory testing be conducted to obtain accurate rock mechanical parameters, including anisotropic elastic properties, strength parameters, and permeability under in-situ stress conditions. Real-time monitoring tools such as measurement while drilling (MWD), logging while drilling (LWD), and pressure while drilling systems should be integrated with geomechanical models to dynamically update mud weight requirements. Sensitivity and uncertainty analysis should be performed to account for variations in stress estimation, pore pressure prediction, and rock property measurements. Field validation of numerical simulations through comparison with drilling events and wellbore image logs is essential to improve model reliability. Finally, adopting adaptive mud weight control and managed pressure drilling techniques is strongly recommende to minimize instability risks and optimize drilling performance in challenging geological environments.

Future Scope overall

Future research should focus on incorporating thermohydromechanical (THM) coupling to simulate more realistic reservoir conditions, especially in high-temperature and high-pressure environments. Advanced modeling of anisotropic and heterogeneous formations, including naturally fractured and faulted systems, can improve predictive accuracy. The integration of machine learning techniques with geomechanical simulation could enable real-time prediction of stability risks using drilling data. Development of digital twin systems that continuously update the safe mud window based on live field measurements represents a promising advancement. Additionally, coupling reservoir depletion effects with wellbore stability modeling will improve long-term integrity assessment for production wells.

Research Contribution

This research contributes to the field of drilling geomechanics by developing a coupled analyticalnumerical framework for safe mud weight determination. It integrates stress analysis, pore pressure modeling, and failure criteria into a unified simulation-based approach. The study improves prediction of collapse and fracture pressures under complex stress conditions and provides a validated methodology for practical drilling applications. The incorporation of numerical simulation enhances understanding of stress redistribution and failure propagation around the wellbore. The proposed framework supports improved drilling design, operational safety, and risk mitigation in challenging geological formations.

Bibliography :-

Books

1. Aadnøy, B. S., & Looyeh, R. (2019). Petroleum Rock Mechanics: Drilling Operations and Well Design. Gulf Professional Publishing.
2. Fjaer, E., Holt, R. M., Horsrud, P., Raaen, A. M., & Risnes, R. (2008). Petroleum Related Rock Mechanics (2nd ed.). Elsevier.
3. Jaeger, J. C., Cook, N. G. W., & Zimmerman, R. (2007). Fundamentals of Rock Mechanics (4th ed.). Blackwell Publishing.

   Research Articles

4. Zoback, M. D. (2007). Reservoir Geomechanics. Cambridge University Press.
5. Aadnøy, B. S. (1987). Inversion technique to determine the in-situ stress field from borehole data. Journal of Petroleum Technology, 39(11), 13691376.
6. Moos, D., & Zoback, M. D. (1990). Utilization of observations of well bore failure to constrain the orientation and magnitude of crustal stresses. Journal of Geophysical Research, 95(B6), 93059325.
7. Eaton, B. A. (1975). The equation for geopressure prediction from well logs. Society of Petroleum Engineers Journal, 15(2), 135145.
8. Biot, M. A. (1941). General theory of three-dimensional consolidation. Journal of Applied Physics, 12, 155164.
9. Mohr, O. (1900). Welche Umstände bedingen die Elastizitätsgrenze und den Bruch eines Materials? Zeitschrift des Vereins Deutscher Ingenieure.
10. Coulomb, C. A. (1776). Essai sur une application des règles des maximis et minimis à quelques problèmes de statique. Academie Royale des Sciences.

    Software / Technical References

11. COMSOL Multiphysics Documentation Structural Mechanics & Porous Media Modules.
12. ABAQUS User Manual Geomechanics Simulation Module.

Acknowledgement

I would like to express my sincere gratitude to my supervisors for their continuous guidance, valuable suggestions, and constructive feedback throughout the course of this research. Their expertise and encouragement have significantly contributed to the successful completion of this work. I am deeply thankful to the faculty members and department for providing the necessary academic support, research facilities, and technical resources that made this study possible. I also acknowledge the assistance and cooperation received from laboratory staff, research colleagues, and industry professionals who contributed valuable insights and data for this research. Finally, I extend my heartfelt appreciation to my family and friends for their constant support, motivation, and understanding during the entire period of this project. Their encouragement has been instrumental in accomplishing this research work successfully.