Accurate prediction of hydraulic fracture propagation is vital for Enhanced Geothermal System (EGS) design. We study the first hydraulic fracturing job at the GR1 well in the Gonghe Basin using field data, where the overall direction of hydraulic fractures does not show a delineated shape parallel to the maximum principal stress orientation. A field-scale numerical model based on the distinct element method is set up to carry out a fully coupled hydromechanical simulation, with the explicit representation of natural fractures via the discrete fracture network (DFN) approach. The effects of injection parameters and in situ stress on hydraulic fracture patterns are then quantitatively assessed. The study reveals that shear-induced deformation primarily governs the fracturing morphology in the GR1 well, driven by smaller injection rates and viscosities that promote massive activation of natural fractures, ultimately dominating the direction of hydraulic fracturing. Furthermore, the increase of in situ differential stress may promote shear damage of natural fracture surfaces, with the exact influence pattern depending on the combination of specific discontinuity properties and in situ stress state. Finally, we provide recommendations for EGS fracturing based on the influence characteristics of multiple parameters. This study can serve as an effective basis and reference for the design and optimization of EGS in the Gonghe basin and other sites. (c) 2025 Institute of Rock and Soil Mechanics, Chinese Academy of Sciences. Published by Elsevier B.V. This is an open access article under the CC BY-NC-ND license (http://creativecommons.org/licenses/by-nc-nd/ 4.0/).

Wellbore breakout is one of the critical issues in drilling due to the fact that the related problems result in additional costs and impact the drilling scheme severely. However, the majority of such wellbore breakout analyses were based on continuum mechanics. In addition to failure in intact rocks, wellbore breakouts can also be initiated along natural discontinuities, e.g. weak planes and fractures. Furthermore, the conventional models in wellbore breakouts with uniform distribution fractures could not reflect the real drilling situation. This paper presents a fully coupled hydro-mechanical model of the SB-X well in the Tarim Basin, China for evaluating wellbore breakouts in heavily fractured rocks under anisotropic stress states using the distinct element method (DEM) and the discrete fracture network (DFN). The developed model was validated against caliper log measurement, and its stability study was carried out by stress and displacement analyses. A parametric study was performed to investigate the effects of the characteristics of fracture distribution (orientation and length) on borehole stability by sensitivity studies. Simulation results demonstrate that the increase of the standard deviation of orientation when the fracture direction aligns parallel or perpendicular to the principal stress direction aggravates borehole instability. Moreover, an elevation in the average fracture length causes the borehole failure to change from the direction of the minimum in-situ horizontal principal stress (i.e. the direction of wellbore breakouts) towards alternative directions, ultimately leading to the whole wellbore failure. These findings provide theoretical insights for predicting wellbore breakouts in heavily fractured rocks. (c) 2025 Institute of Rock and Soil Mechanics, Chinese Academy of Sciences. Published by Elsevier B.V. This is an open access article under the CC BY-NC-ND license (http://creativecommons.org/licenses/by-nc-nd/ 4.0/).

Borehole instability in naturally fractured rocks poses significant challenges to drilling. Drilling mud invades the surrounding formations through natural fractures under the difference between the wellbore pressure (Pw) and pore pressure (Pp) during drilling, which may cause wellbore instability. However, the weakening of fracture strength due to mud intrusion is not considered in most existing borehole stability analyses, which may yield significant errors and misleading predictions. In addition, only limited factors were analyzed, and the fracture distribution was oversimplified. In this paper, the impacts of mud intrusion and associated fracture strength weakening on borehole stability in fractured rocks under both isotropic and anisotropic stress states are investigated using a coupled DEM (distinct element method) and DFN (discrete fracture network) method. It provides estimates of the effect of fracture strength weakening, wellbore pressure, in situ stresses, and sealing efficiency on borehole stability. The results show that mud intrusion and weakening of fracture strength can damage the borehole. This is demonstrated by the large displacement around the borehole, shear displacement on natural fractures, and the generation of fracture at shear limit. Mud intrusion reduces the shear strength of the fracture surface and leads to shear failure, which explains that the increase in mud weight may worsen borehole stability during overbalanced drilling in fractured formations. A higher in situ stress anisotropy exerts a significant influence on the mechanism of shear failure distribution around the wellbore. Moreover, the effect of sealing natural fractures on maintaining borehole stability is verified in this study, and the increase in sealing efficiency reduces the radial invasion distance of drilling mud. This study provides a directly quantitative prediction method of borehole instability in naturally fractured formations, which can consider the discrete fracture network, mud intrusion, and associated weakening of fracture strength. The information provided by the numerical approach (e.g. displacement around the borehole, shear displacement on fracture, and fracture at shear limit) is helpful for managing wellbore stability and designing wellbore-strengthening operations. (c) 2024 Institute of Rock and Soil Mechanics, Chinese Academy of Sciences. Production and hosting by Elsevier B.V. This is an open access article under the CC BY-NC-ND license (http://creativecommons.org/ licenses/by-nc-nd/4.0/).