The weakening of rock mass after immersion substantially affects slope stability. Considering the mudstone in the north slope of the Fushun West Open pit as the research object, we investigated the failure and creep characteristics of mudstone under different immersion times and confining pressures. For this, we used a self-developed experimental apparatus that applies biaxial compression loading in a water immersion environment. Water significantly changed the strength and failure mode of the mudstone. The increase in confining pressure limits the expansion of the crack aperture, reduces the rock strength-weakening coefficient, and delays the rock weakening time. The change in cohesion is primarily responsible for the weakening of shear strength during immersion. Based on the creep test results, the creep curve of the soaking specimen was identical to that of the dry specimen in the initial phase. The weakening effect of water on the specimen increased as the soaking duration increased, due to which the creep curve of the soaking specimen gradually approached that of the saturated specimen. Water entered the microcracks and reacted with the rock interior, accelerating rock failure. Consequently, the soaking specimen had the lowest stress level for creep failure and the shortest failure time. The results of this experiment provide a reference for the long-term stability of the flooded slope.

To evaluate the effect of hydraulic fracturing in a low-permeability reservoir, a propagation path tracing method for hydraulically created fractures was established based on microseismic monitoring data. First, the numerical simulation of the wave propagation process, grid search, and error-weight coefficient method was combined to locate the microseismic source. Then, the moment tensor inversion method was used to determine the tensile angle and source mechanism of hydraulically created fractures. Next, the tensile angle was used as the weight-index to determine the size of the mixed-source mechanism fracture combined with the shear and tensile source size quantization model. Finally, the spatial topological relationship between fractures was determined by comprehensively considering the spatial location and radius of the fractures, to realize the propagation path tracing of hydraulically created fractures. These tracking results can be used as one of the bases for the evaluation of the hydraulic fracturing effect.

This study develops a discrete element model incorporating the water–ice phase transition volume effect to simulate frost damage in saturated granite. The model investigates the damage evolution and mechanical degradation under freeze–thaw cycles. The results show that during freeze–thaw cycles, the model’s temperature field exhibits non-uniform distribution characteristics and geometric dependency, with lower maximum temperature differences in Brazilian disk models versus uniaxial compression specimens. Frost heave damage progresses through three distinct stages: localized bond fractures (1~5 cycles); accelerated crack interconnection and branching (15~20 cycles); and fully interconnected damage zones (25~30 cycles). As the number of freeze–thaw cycles increases, the crack network significantly influences the mechanical behavior of the model under load. The failure mode of the loaded model undergoes a transformation from brittle penetration to ductile fragmentation. Freeze–thaw cycles cause more significant degradation in the tensile strength of granite compared to compressive strength. After 30 freeze–thaw cycles, the uniaxial compressive strength and Brazilian tensile strength decrease by 47.5% and 93.8%, respectively. These findings provide theoretical support for assessing frost heave damage in geotechnical engineering in cold regions.

Previous studies have mainly examined sandstone freeze–thaw degradation under unloaded conditions, whereas the continuous influence of sustained stress on compressive and tensile properties remains unclear. To address this issue, water-saturated specimens were subjected to axial loads of 0–2 MPa and 0–10 freeze–thaw cycles, followed by uniaxial compression and Brazilian splitting tests. A two-variable exponential strength model and a load–freeze–thaw coupled discrete element model were established. Both strengths decreased with increasing freeze–thaw cycles. Under no load, after 10 cycles, the uniaxial compressive strength and tensile strength decreased by 27.85% and 36.67%, respectively, indicating higher freeze–thaw sensitivity of the tensile property. Loading mitigated strength loss: after 10 cycles, the two strengths of the 2 MPa group were 9.59% and 10.53% higher than those of the no-load group. The retention effect on compressive strength first increased and then decreased, whereas that on tensile strength increased with cycle number. The discrete element results showed that frost-heave cracks evolved from local initiation to connection and clustering, with tensile cracks dominating. Loading reduced tensile bond breakage and slowed crack accumulation and subsequent coalescence, linking the observed strength retention to the inhibition of mesoscopic tensile damage.