Rock Mass Research Articles

Diffusion mechanism of quick-setting slurry in water-rich fractured rock mass based on circle-outburst diffusion model

Grouting is the most commonly used methods in dealing with water inrush issue in mine and tunnel engineering. In order to better predict the grouting effect, research on slurry diffusion mechanism became a hotspot for scholars. At present, the mainstream theoretical model used to study the slurry diffusion mechanism are the circle diffusion model and the modified ones in plane plate fracture. However, these models cannot explain the “contradiction” between the rapid setting characteristics of quick-setting slurry and the continuous long-term injection of slurry in grouting engineering, which cannot accurately predict the grouting pressure or grouting diffusion range. Therefore, a “circle-outburst diffusion model” which can explain the above “contradiction” was proposed in this paper. Based on the new model, stepwise algorithms were developed to predict the grouting pressure and slurry diffusion range. By means of conducting fracture grouting simulation test and collecting grouting pressure data in an actual grouting project, the new model was verified and the slurry diffusion mechanism was studied. Comparison results indicate that the new model can reveal the intrinsic reason for the irregular diffusion phenomenon of the slurry and forecast the outburst moment accurately. The circle diffusion stage and outburst diffusion stage elaborated in the new “circle-outburst diffusion model” were consistent with the staged characteristics of the grouting process presented in the grouting simulating test. Differences between theoretical and experimental pressure value were within 10 %, indicating a high degree of consistency. The number and distribution of outburst diffusion points determine slurry diffusion form and significantly influence the diffusion range. The grouting pressure and the length of slurry flow path increase nonlinearly with time in each diffusion stage. It is hoped that the new model can provide a theoretical basis for the study of diffusion mechanism of quick-setting slurry.

Multifactor-coupled study on freeze-thaw forces of rocks in cold regions

This study delves into the mechanical properties of steep and rocky slopes subjected to long-term freeze-thaw actions. Considering the unique climatic conditions in cold regions, especially the significant impact of seasonal and diurnal temperature variations on slope excavation, the research focuses on a high-cold region iron ore mine. Four types of rocks commonly found in the mining area are thoroughly examined, taking into account the hydrogeological conditions of the mining area. The study systematically analyzes the mechanisms of various factors such as weathering, freeze-thaw cycles, and ice-water phase changes on the stability of cold region fractured rock masses. The research reveals that under prolonged freeze-thaw actions, crack water within the rock continuously undergoes ice-water phase changes, generating substantial freeze expansion forces that result in structural damage to the rock mass. This damage is evident not only in the development of existing microcracks but also leads to the generation of new fractures, ultimately causing deterioration in the rock mass structure. The study of the evolution patterns of freeze-thaw forces contributes to a better understanding of slope stability issues in cold region mineral resource extraction, offering crucial insights for the design, construction, and operation of related engineering projects.

Study on the evolution of mechanical properties of hot dry rocks after supercritical CO2 injection

Characterizing the evolution of mechanical properties of hot dry rock (HDR) after supercritical CO2 (CO2(sc)) injection is crucial for assessing the heat extraction rate and reservoir security of CO2 based enhanced geothermal systems. This study designed the experiments of triaxial seepage and mechanical properties considering no CO2(sc) injection, CO2(sc) injection, and alternating injection of water-CO2(sc) (AIWC) in granite at 150–300 ℃. The experiments can reveal the mechanical properties of HDR in single-phase CO2 zone, CO2-water two-phase zone and dissolved CO2 liquid phase zone in HDR reservoir. The results indicate that the failure mode of the rock samples primarily exhibits sudden instability after no CO2(sc) injection and AIWC, whereas it predominantly manifests progressive instability after CO2(sc) injection. Compared with 25 ℃, the uniaxial compressive strength (UCS) after no CO2(sc) injection at 150–300 ℃ decreased by 13.86%–32.92%. After CO2(sc) injection, the UCS decreased by 40.79%–59.60%. After AIWC, the UCS decreased by 27.74–40.48%. This shows that the strength of rock mass in the single-phase CO2 zone is lower than that in the other two zones, and this weakening phenomenon increases with the increase of temperature difference. At the same temperature, the elasticity modulus after AIWC was greater than that after no CO2(sc) injection and CO2(sc) injection. With no CO2(sc) injection, when the temperature was increased to 200 ℃ and 300 ℃, intergranular cracks and transgranular appeared respectively. After AIWC, mineral crystals such as calcite were precipitated on the surfaces of the connected large cracks, accompanied by kaolinite clay minerals. This increases the frictional contact of the mineral particles and enhances the stability of the HDR reservoir.

Failure mechanisms and reinforcement support of soft rock roadway in deep extra-thick coal seam: A case study

The soft rock roadway in deep extra-thick coal seam is prone to large deformation disaster. In order to solve this problem effectively, taking the W3101 coalface of Xiao’kang coal mine as the engineering background, the failure mechanisms and the corresponding control technology were investigated by the methods of laboratory test, numerical simulation and field measurement. First, a series of tests were carried out on the basic mechanical properties of rock mass, in-situ stress and loose zone, and it was found that the failure mechanisms of soft rock roadway in deep extra-thick coal seam are poor lithology of surrounding rock, high in-situ stress, large loose circle and unreasonable original support. Then, based on the supporting mechanisms of grouting cable and concrete-filled steel tube (CFST) of bottom angle, a new collaborative reinforcement support (CRS) technology combining the two was proposed. Next, the numerical simulation shows that the CRS technology has good applicability in controlling large deformation of soft rock roadway in deep extra-thick coal seam. Finally, to exhibit the engineering application effect of the proposed CRS scheme, it was also applied in Xiao’kang coal mine. Comparing to the controlling effect of original support scheme, it can be seen that the average deformation of soft rock roadway reduces by 73.2 % under the proposed CRS scheme, and the steel arch frame of bottom angle is slightly deformed. Meanwhile, it is found that the grouts diffusion radius is nearly 1.21 m, and the deepest diffusion depth is 4.25 m, which also indicates again that the feasibility of the proposed CRS scheme in controlling large deformation of soft rock roadway in deep extra-thick coal seam.

Study on mechanical properties and damage characteristics of acid corrosion in granite based on NMR technology

In acidic environments, rock masses are frequently subjected to severe chemical corrosion, resulting in the initiation of numerous geological engineering disasters. This study aimed to collect physical and mechanical parameters of granite exposed to prolonged acid corrosion and analyze fracture characteristics using acoustic emission (AE) techniques. Additionally, it examined the evolution of pore structure and damage mechanisms through the use of low-field nuclear magnetic resonance (NMR) and fractal theory. The results demonstrate a monotonic decrease in mass, volume, density, P-wave velocity, and S-wave velocity of granite with increasing corrosion time. Particularly notable is the phased reduction observed in uniaxial compressive strength and elastic modulus. The transition from brittle to ductile failure in corroded granite is accompanied by a gradual decrease in internal fracture strength. The trend in the correlation dimension reveals the relationship between the formation time of the main fracture surface and the pore structure. Additionally, total porosity and macropores (D, Da) exhibit significant fractal characteristics. The fractal dimension correlates positively with the damage variable and inversely with uniaxial compressive strength and elastic modulus. This indicates that more severe pore structure damage leads to a higher fractal dimension and lower mechanical performance. Among these, Da demonstrates higher sensitivity in characterizing rock mechanical properties. These findings provide important basis for evaluating the stability of granite geotechnical engineering in acidic environments.

Damage deformation properties and acoustic emission characteristics of hard-brittle rock under constant amplitude cyclic loading

In order to study the deformation and damage characteristics of the limestone specimens with high strength and brittleness under constant amplitude cyclic loading, the deformation and the acoustic emission (AE) characteristics were analysed, and the relationship between them was sought. The damage variables under different amplitude cyclic loading were defined by AE counts. The results showed that the radial deformation of the limestone specimens was more sensitive and unstable than the axial deformation. The concept of apparent residual strain was proposed to describe the specimen deformation characteristics, and it resulted that the radial apparent residual strain produced at higher stress state would recover at lower stress state. The limestone specimens showed obvious Kaiser effect and Felicity effect under cyclic loading. When the upper limit of the cyclic loading was close to the peak stress of the specimen, the AE counts generated in unloading sections were almost the same as that in the loading sections. The damage was increased as the amplitude and the stress level increased and the unloading process at higher stress level would also lead to the aggravation of damages. Specimens would absorb more energy under cyclic loading than under uniaxial loading. Reasonable driving parameters should be controlled in underground excavation practice, to ensure that the stress level of surrounding rock mass in a periodic stress state is located before peak stress and such that to limit the occurrence of rock burst to a certain extent.

Experimental and meshless numerical simulation on the crack propagation processes of marble SCB specimens

In order to correctly understand the rock multi-fissure interactions and crack propagation laws, three-point bending fracture tests of the marble semi-circular bend samples (SCB samples) are carried out, and the crack propagation processes are simulated by the Smoothed Particle Hydrodynamics (SPH) method. Results show that: Cracks initiate from the guide fissure tip and propagate vertically in the SCB sample with no obstacle fissures; The overlap positions of crack1 and the obstacle fissure are biased towards the lower left end of the obstacle fissure with the increase of the inclination angle α; With the increase of the eccentric distance d, the overlap positions of crack1 and the obstacle fissure are biased towards the upper right end of the obstacle fissure; The peak load of the samples decreases with the increase of inclination angle α, while increases with the increase of eccentric distance d. The “failure coefficient” η is introduced into the conventional SPH momentum equations, which can simulate the multi-fissure propagations of the SCB samples. The simulation results are consistent with the experimental results, and the mechanisms of the crack propagation morphologies between the guide fissure and the obstacle fissure are discussed. The research results can provide some references for correctly understanding the crack propagation mechanisms of multi-fissure rock masses and the applications of the SPH method into simulating rock multi-fissure interactions.

Crack coalescence prediction and load-bearing mechanism of defective specimen based on computer vision recognition model

The coalescence of cracks in rock, triggered by external stress or geological activities, plays a pivotal role in determining the mechanical properties and stability of rock structures. Consequently, the prediction method of crack coalescence become critical tools in ensuring the safety and stability of geotechnical engineering facilities. In this paper, based on the strain field data obtained by digital image correlation (DIC), a set of programs was developed to automatically identify the values and percentage changes of different strain intervals in the strain field of the specimen. Then, a computer vision recognition (CVR) model is established to study the process and prediction of crack coalescence in sandstone specimens with open flaws. This method overcomes the limitations, subjectivity and unpredictability of the traditional method of identifying cracks through artificial vision. The results show that the increase of the flaw width causes the displacement trend to be squeezed and deflected along the width direction, thereby changing the angle of crack initiation and exhibiting different forms of crack coalescence. The increase in the inclination angle of the flaw leads to an increase in the effective compressive area (ECA) of sandstone, and the magnitude of ECA is positively correlated with the uniaxial compressive strength (UCS) of the sandstone. Additionally, the CVR model identifies two types of numerical fluctuation signals before crack coalescence, namely the plastic characteristic signal (PCS) in the early stage of the experiment and the early-warning signal (EWS) near the period of crack coalescence, where EWS can be used as a predictive signal for crack coalescence. The research results provide a new algorithm technical support for the process analysis and prediction of crack coalescence, and provide a basis for early warning of rock mass engineering disasters.

Quantifying uncertainty in rock mass properties: Implications for GSI, RMi, and RMR assessments

Probability-based empirical methods were employed as an alternative approach to predicting uncertainties associated with rock mass properties. The focus was on developing probabilistic spreadsheets to forecast rock mass classification indexes. Histograms were constructed to describe the best distribution in predicting rock mass properties. The developed models also offer utility in predicting the impact of discontinuities within the rock mass on rock strength and rock mass classification systems. Statistical analyses identified volumetric joint count, joint spacing, joint frequency, and rock strength as the most influential parameters. Moreover, the statistical analysis revealed varying degrees of correlation among different rock mass properties. While some properties demonstrated significant correlations suitable for modelling, others did not align well with any correlation model. The results highlight the need for a comprehensive approach to rock mass characterization, considering multiple factors beyond volumetric joint count. Geological complexities, including tectonic activity and weathering processes, may obscure direct correlations. These results emphasize the importance of empirical modelling and detailed site investigations for accurate assessment of rock mass quality and stability in the Himalaya.

Fractal evolution characteristics of fracture meso-damage in uniaxial compression rock masses using bonded block model

With regard to deep mining in metal mines, an investigation into the failure mode of deep fractured rock masses and their corresponding acoustic emission signal characteristics is conducted via uniaxial compression tests. Subsequently, a fractal damage renormalization group mechanical model is developed to explain the behavior of those fractured rock masses. Employing the bonded block model (BBM) numerical simulation method, fracture process in synthetic rock samples is analyzed, thereby validating the efficacy of the mechanical model. The numerical simulations highlight the critical role of fractures expansion in underlying the deterioration of rock mass strength. As the peak load decreases, the fracture fractal dimension increases, leading to a significant 14.2% reduction in compressive strength accompanied by an approximate 8.7% rise in average fracture fractal dimension. A comparative analysis of tetrahedral and voronoi block synthetic rock samples reveals the tetrahedral block samples exhibit a superior ability to depict the fracture behavior of fractured rock masses. Specifically, they offer a more accurate simulation of acoustic emission characteristics and failure modes. Furthermore, variations in the fracture fractal dimension with respect to the hole defect’s position are observed, with the maximum value occurring along the vertical axis of the hole defect. This observation underscores the potential utility of visually monitoring deep rock fracture dynamics as an effective mean for quantitatively evaluating fracture damage and strength degradation in deep rock formations.

A novel tool for seismic response analysis of tunnel in multilayered media based on kinematic earthquake source

Earthquake damage of underground structure indicates that the seismic response of tunnels is affected by the earthquake source, propagation medium, and local tunnel-surrounding rock system. Current research on the seismic response of tunnels primarily focuses on local tunnel systems, posing a challenge in comprehensively assessing the impact of earthquakes on tunnel group. This study introduces a theoretical seismological approach for the seismic response analysis of tunnels and presents a novel methodology for assessing the response. The proposed approach uses a kinematic seismic source as the earthquake source, assuming the geological strata as a layered medium and a tunnel as a point within a semi-infinite space. The wavenumber-integral method was used to solve for the seismic wave field and strain. To reflect the followability of underground structure to the surrounding rock deformation, we improved it and derived an analytical formula for the strain. The 2022 Menyuan earthquake was considered as a case study; the seismic damage to tunnels along the Lan-Xin Railway was analyzed using the proposed approach. The results confirm that there is a strong correlation between tunnel damage and free-field strain induced by earthquakes. The theoretical influence range of earthquake on tunnel was investigated using the diameter deformation ratio of the tunnel as an indicator. For tunnels embedded in a grade III rock mass, the theoretical range was approximately 7.0 km on both sides of the fault. For tunnels embedded in a grade V rock mass, the theoretical range can reach up to 25.9 km. The results also revealed that tunnels on the hanging-wall side of the fault had a more significant impact than those on the foot-wall side, although their affected areas were smaller.

Magmatic and sedimentological arguments for an Ediacaran active margin in the Bayankhongor Zone in western Mongolia, Central Asian Orogenic Belt

Geochronological and geochemical investigation of the magmatic and sedimentary rocks from the south-eastern tip of the Bayankhongor Zone (central Mongolia) constrains the Neoproterozoic evolution of the northern margin of the Baidrag Block. There, ultramafic to felsic igneous rocks of the Khan-Uul Massif intruded the volcano-sedimentary sequence of the Ulziit Gol Unit. U–Pb zircon ages show that both the plutonic and volcanic rocks were coeval products of the same Ediacaran (598–564 Ma) magmatism. Most of the samples have low contents of high field strength elements typical of arc-related magmas; however, some mafic rock samples display N-MORB-like chemistry. Magmatic rocks in the Khan-Uul Massif and Ulziit Gol Unit yielded exclusively positive initial εHf in zircon values (+4.9 to +13.8), implying derivation from depleted-mantle sources. Furthermore, heterogeneity of the whole-rock initial εNd values (–4.0 to +2.1) documents that fractional crystallization was accompanied by variable crustal contamination. The detrital zircon age patterns of the sandstone and tuffites in the Ulziit Gol Unit indicate that the sequence was filled dominantly by Ediacaran magmatic detritus derived from the Khan-Uul Massif and its volcanic equivalents, while the older cratonic material from the Baidrag basement represented only a subordinate component. Combined magmatic and sedimentary records imply that both the Khan-Uul Massif and the Ulziit Gol Unit may have formed in the same back-arc basin environment. This challenges the previous view that the entire Bayankhongor Zone was ophiolitic in nature. Distribution of coeval Ediacaran supra-subduction systems on the scale of the Mongolian Collage indicates the closure of multiple oceanic basins rimming the Siberian Continent. Following Rodinia break-up, this peri-Siberian realm was presumably formed by sub-parallel continental ribbons and oceanic basins. It is proposed that the amalgamation of these blocks and closure of intervening oceans reflected the Ediacaran advancing mode of the Palaeo-Pacific subduction.

Seepage-heat accumulation characteristics and hydrothermal transport calculation model in mining-induced fractured rock mass: Implications for geothermal water control

The aim of this paper is to explore hydrothermal transport in surrounding rocks during deep mining, focusing on the mechanisms and key factors influencing groundwater seepage and heat accumulation disasters. Stress–seepage coupling tests were conducted on fractured rock mass with varying roughness under constant-pressure loading/unloading conditions. Also, a permeability model considering joint roughness coefficient and confining pressure was established. The study further examined the evolution of heat transfer in rough fractured rock mass and variations in convective heat transfer characteristics. The heat transfer process was observed to progress through stages of surge, attenuation, and stabilization. Two calculation models were developed using multivariate linear regression and random forest regression to quantify how different factors influenced water–rock convective heat transfer in mining-induced fractured rock mass. The average correlation coefficient R2 of the proposed models exceeded 0.95. The importance scores of osmotic pressure, inlet water temperature, rock temperature, and permeability were 0.417, 0.296, 0.223, and 0.065, respectively. The flow rate and water–rock temperature difference were determined to be the key factors influencing hydrothermal transport. This study enhanced the understanding of the hydrothermal migration mechanisms in surrounding rocks during mining and provided insights into geothermal water control.

Dynamic impact mechanical properties of red sandstone based on digital image correlation method

ABSTRACT In the process of underground blasting excavation, the safety and stability of rock mass are significantly influenced by dynamic load. To study the full-field strain and displacement behaviour of rock materials subjected to dynamic compressive force, the impact tests were done on red sandstone specimens using an electromagnetically driven Hopkinson pressure bar. At the same time, the data were collected by a high-speed camera, and non-contact measurements of the red sandstone were carried out with the digital image correlation method. The resulting images were screened and analysed, based on which the stress-strain curve, full-field strain cloud diagram and displacement cloud diagram were obtained. Results show that the dynamic peak stress increases with the increase of impact velocity. During the impact process, substantial internal uneven deformation occurs in the red sandstone close to the transmission bar’s end face. Under the impact load, the cracks initiate from the upper and lower sides near the transmission bar’s end face, followed by formation of cracks in the middle part. These cracks expand throughout the rock to form through-cracks, eventually resulting in complete failure. For rock materials, the evolution of their strain field during the failure process can reflect the initiation and propagation of their internal cracks. Therefore, the use of digital image correlation method proves effective for studying the deformation and failure mechanisms of rock. The conclusions obtained in this study provide a valuable reference for the macro-meso deformation and failure mechanisms of geotechnical materials.

Beyond the empirical pillar design method: The strain criterion and the pillar load inversion concepts

Pillar design is a crucial aspect of underground mining engineering, as it directly impacts the safety, stability, and overall effectiveness of mining operations. In this paper we present a method to calculate pillar stability based on the strain criterion and pillar load inversion concept, which complements the empirical pillar strength formulae. Those formulae do not account for certain parameters that can affect the stability of the pillars, e.g. weak layers, changes in stress orientations due to mining, orebody dip, mining layout complexity, and the influence of regional stabilizing pillars. The approach presented here is not intended to replace the empirical method; rather, we suggest using it as the initial step in the pillar design process. This should be followed by iterative numerical modelling to study pillar behaviour, define pillar stability, and optimize the pillar design by considering site-specific and representative rock mass properties and criteria.

Probabilistic evaluation method for the stability of large underground cavern considering the uncertainty of rock mass mechanical parameters: A case study of Baihetan underground powerhouse project

For deep large underground cavern projects under complex geological conditions, the determination of rock mass mechanical parameters and stability analysis is fraught with a great deal of uncertainty. This paper proposes a set of methods for estimating rock mass mechanical parameters and probabilistic stability assessment suitable for construction characteristics of large underground caverns. On the one hand, according to the Hoek-Brown criterion and Monte Carlo simulation, a dynamic estimation method for rock mass mechanical parameters is proposed by using Rock Mass Rating (RMR), the uniaxial compressive strength (UCS) of intact rock, and the material constant (mb) as input parameters. On the other hand, considering the uncertainty of rock mass mechanical parameters, a probabilistic evaluation method applicable to the unloading response characteristics of rock mass (including displacement and fracturing zones) around the large cavern is put forward, combining the point estimation method and numerical simulation. The proposed methods were applied to an underground powerhouse, with an excavation span of 34 m and height of 88.7 m, at the Baihetan hydropower station in the southwest of China. Based on extensive field investigations and tests, the probability distributions of key rock mass mechanical parameters were obtained. Furthermore, the high sensitivity of input parameters in the H-B criterion for estimating rock mass mechanical parameters were revealed. Through dynamic simulation, the probability distributions of the displacement and fracturing zones in surrounding rock during the layer-by-layer excavation of the cavern were presented. Comprehensive on-site tests showed that the simulated results were in basic agreement with the actual unloading response behavior during the cavern excavation, validating the correctness and applicability of the proposed methods. The study provides a relatively comprehensive approach for estimating rock mass mechanical parameters and stability evaluation in similar underground engineering projects, and also has guiding significance for predicting and preventing engineering rock mass hazards.

An analytical modelling of the rotation of in-situ principal stresses orientation in the vicinity of an underground opening

ABSTRACT An underground excavation would lead to the redistribution of not only stress magnitude but also stress direction in the host rock mass. In this paper, an equation was derived to quantify the principal stress orientation rotation angle along the roadway radial direction. Besides, the principal stress orientation rotation gradient was analytically modelled followed by a sensitivity analysis of various influential factors on the principal stress orientation rotation. It was found that the principal stress orientation rotation angle tends to decrease as moving further radially into the rock mass from the excavation surface, where the principal stress rotation angle reaches its maximum value of 90°. Four zones have been classified according to the principal stress orientation rotation gradient including high-frequency , medium-frequency, low-frequency and no rotation zone. Finally, the numerical analysis confirmed the redistribution of the principal stress difference and the rotation angle and provided guidance for optimising the bolts setting angle.

Study on dynamic mechanical properties and microscopic damage mechanisms of granite after dynamic triaxial compression and thermal treatment

As underground mining operations gradually extend deeper, the conditions for orebody occurrence become increasingly complex, and various geological disasters occur frequently. Rock masses are prone to different degrees and types of damage, making it impractical to continue using intact rock as a reference. To study the dynamic mechanical properties of damaged rock under actual conditions, this study subjected granite samples to impact and high-temperature damage. Detailed observations were made of the samples' surface morphology and microstructure before and after damage, and the patterns of damage changes were investigated. Subsequently, uniaxial compression tests at different loading rates were conducted on the damaged samples. By calculating the loading rate effect sensitivity, it was found that as the damage increased, the rate effect gradually diminished. In addition, this study also summarized the influence of damage and loading rate on the macroscopic failure characteristics of the samples. The novelty of this study lies in focusing on damaged rock, which more closely resembles the actual rock conditions encountered in most underground mining operations today. This research can provide a reference for stability assessment and safe construction in deep mining rock engineering and offers important insights into the feasibility of non-explosive extraction of damaged rock.

The Influence of Grain Size Gradation on Deformation and the Void Structure Evolution Mechanism of Broken Rock Mass in the Goaf

The Influence of Grain Size Gradation on Deformation and the Void Structure Evolution Mechanism of Broken Rock Mass in the Goaf

Thermo-mechanical behavior of sandstone joints: Findings from direct shear tests

Rock masses typically contain joints and defects affecting their safety and stability under compression, shear, or tension. The shear strength of these discontinuities, influenced by environmental factors like temperature, is critical in deep rock engineering. This study conducted direct shear tests on sandstone joints at room temperature (RT) and after thermal treatment at 250 °C and 500 °C. The behavior under compressive and shear forces was investigated by using the direct shear test. The results showed that compressive and shear properties improved at 250 °C but deteriorated following thermal treatment at 500 °C. Peak shear strength before and after thermal treatment closely followed Patton's bilinear failure criterion. As intact asperities were sheared under normal stress of 0.75MPa, the residual shear strength was less affected by temperature than peak shear strength. The shear strength of intact asperities increased by 8.8% following thermal treatment at 250 °C and decreased by 4.3% at 500 °C. The joint closure, normal stress, and stiffness were consistent with a modified Bandis et al. (1983) criterion. Increasing the thermal treatment temperature from RT to 250 °C and 500 °C changed maximum joint closure from 0.201 to 0.125and 0.301 mm, and subsequently, the initial normal stiffness from 8 to 19.6 and 6.62MPa/mm, respectively. These findings offer valuable insights into the temperature-dependent behavior of sandstone joints, aiding the design and maintenance of deep underground structures. This research is particularly relevant to engineering geology, where understanding the thermo-mechanical behavior of rock joints is essential for projects such as geothermal energy extraction and the storage of nuclear waste. The results can improve the safety and efficiency of geological engineering practices in thermally active environments.