The creep phenomenon of inelastic deformation of surrounding rock may occur under the action of deep geological stress for a long period of time, potentially resulting in large-scale deformations or even instability failure of the underground engineering. Accurate characterization of the creep behavior of the surrounding rock is essential for evaluating the long-term stability and safety of high-level radioactive waste (HLW) disposal repositories. Although the laboratory creep tests of brittle undamaged rocks, such as granite, have been extensively performed, the creep characteristics of fractured surrounding rock under the multi-field coupling environment still require attention. In this study, a series of creep experiments was conducted on Beishan granite, which was identified as the optimal candidate surrounding rock for the disposal repository in China. The effects of various factors, including inclination angle of fractures, stress conditions, temperatures, and water content, were investigated. The experimental results show that the axial total strain increases linearly with increasing stress level, while the lateral total strain, axial and lateral creep strain rates increase exponentially. The failure time of saturated specimens fractured at 45 degrees and 60 degrees is approximately 1.05 parts per thousand and 0.84 parts per thousand of that of dry specimens, respectively. The effect of temperature, ranging from room temperature to 120 degrees C, is minimal, compared to the substantial variations in strain and creep rates caused by stress and water content. The creep failure of specimens fractured at 30 degrees is dominated by rock material failure, whereas the creep failure of specimens fractured at 60 degrees is dominated by pre-existing fracture slip. At a 45 degrees fracture angle, a composite failure mechanism is observed that includes both rock material failure and pre-existing fracture slip. (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/).

Asa calculation method based on the Galerkin variation, the numerical manifold method (NMM) adopts a double covering system, which can easily deal with discontinuous deformation problems and has a high calculation accuracy. Aiming at the thermo-mechanical (TM) coupling problem of fractured rock masses, this study uses the NMM to simulate the processes of crack initiation and propagation in a rock mass under the in fluence of temperature field, deduces related system equations, and proposes a penalty function method to deal with boundary conditions. Numerical examples are employed to con firm the effectiveness and high accuracy of this method. By the thermal stress analysis of a thick-walled cylinder (TWC), the simulation of cracking in the TWC under heating and cooling conditions, and the simulation of thermal cracking of the Swedish & Auml;sp & ouml; Pillar Stability Experiment (APSE) rock column, the thermal stress, and TM coupling are obtained. The numerical simulation results are in good agreement with the test data and other numerical results, thus verifying the effectiveness of the NMM in dealing with thermal stress and crack propagation problems of fractured rock masses. (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/).

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/).

Freeze-thaw failure of frozen rock slope often occurs during engineering construction and mining in cold area, which poses a great threat to engineering construction and people's life safety. The properties of rock mass in cold region will change with the periodic change of temperature, which makes it difficult to accurately evaluate the stability of slope under the action of freeze-thaw cycle by conventional methods. Based on field investigation and literature review, this paper discusses the characteristics of frozen rock mass and the failure mechanism of frozen rock slope, and gives the types and failure modes of frozen rock slope. Then, the research status of frozen rock slope is analyzed. It is pointed out that the failure of frozen rock slope is the result of thermo-hydro-mechanical (THM) coupling. It is considered that freeze-thaw cycle, rainfall infiltration and fracture propagation have significant effects on the stability of frozen rock slope, and numerical simulation is used to demonstrate. The research shows that the safety factor of frozen rock slope changes dynamically with the surface temperature, and the safety factor of slope decreases year by year with the increase of freeze-thaw cycles, and the fracture expansion will significantly reduce the safety factor. Based on the above knowledge, a time-varying evaluation method of frozen rock slope stability based on THM coupling theory is proposed. This paper can deepen scholars' understanding of rock fracture slope in cold area and promote related research work.

A 3D high-resolution subsurface characteristic (HSC) numerical model to assess migration and distribution of subsurface DNAPLs was developed. Diverse field data, including lithologic, hydrogeologic, petrophysical, and fracture information from both in situ observations and laboratory experiments were utilized for realistic model representation. For the first time, the model integrates hydrogeologic characteristics of both porous (unconsolidated soil (US) and weathered rock (WR)) and fractured rock (FR) media distinctly affecting DNAPLs migration. This allowed for capturing DNAPLs behavior within US, WR, and FR as well as at the boundary between the media, simultaneously. In the 3D HSC model, hypothetical 100-year DNAPLs contamination was simulated, quantitatively analyzing its spatiotemporal distributions by momentum analyses. Twelve sensitivity scenarios examined the impact of WR and FR characteristics on DNAPLs migration, delineating significant roles of WR. DNAPLs primarily resided in WR due to low permeability and limited penetration into FR through sparse inlet fractures. The permeability anisotropy in WR was most influential to determine the DNAPLs fate, surpassing the impacts of FR characteristics, including rock matrix permeability, fracture aperture size, and fracture + rock mean porosity. This study first attempted to apply the field-data-based multiple geological media concept in the DNAPLs prediction model. Consequently, the field-scale effects of WR and media transitions, which have been often overlooked in evaluating DNAPLs contamination, were underscored.