Stability Of Rock Mass Research Articles

Digital tools used on the teaching-learning process in geotechnical engineering

This article presents the incorporation of information and communication technologies on the teaching-learning process of geotechnical engineering to improve the quality of education and provide practical knowledge to civil engineering students. The content of this paper is divided into three main modules, which are: Treatment of laboratory tests results, Slope stability and Rock mass stability. For this, some software were used, including 2D limit equilibrium stability analysis, stability analysis of rock wedges and a dynamic mathematics, to obtain and analyze soil shear strength parameters, among others. Furthermore, a digital tool was developed to analyze the results gathered in compaction, one-dimensional consolidation, and direct soil shear tests, in order to clarify the relationship between theoretical concepts and practical results of the tests and analyses and to help students on doubts, in addition to increase their interest and motivation to perform the complete interpretation of the collected data. The activities were conducted in classes of geotechnical disciplines of the undergraduate course in Civil Engineering at the Federal University of Viçosa, aiming to promote active learning and improve teaching quality. Based on the results of an applied feedback questionnaire, it was observed that most students were satisfied with the resources used in the classroom, demonstrating that the implemented digital tools work as a didactic instrument that facilitates learning and comprehension of practical problems, in addition to enabling the resolution of several geotechnical engineering problems much more quickly and efficiently.

Stability Deterioration Analysis of Toppling Dangerous Rock Mass under the Combined Effect of Vibration and Freeze-Thaw Cycles

To investigate the degree of stability deterioration in toppling dangerous rock mass under different action times, firstly, the vibration load generated by train travel in a single action and its decay law were studied, and the theoretical relational equation of rock strength deterioration under the action of vibration loading was obtained by analyzing the increase in half axle expansion of micro fissures inside the rock. Secondly, the freezing stress of the microfissures was determined according to the rock freezing and swelling theory, and the theoretical relational equation of rock strength deterioration under the action of freeze-thaw cycles was obtained by integrating the rock confinement and moisture migration conditions of the fissures wall. The stability coefficient of the toppling dangerous rock mass exhibited a rapid and then slow decreasing trend with increasing number of actions, which was non-linearly negatively correlated with the train running speed and the initial porosity of the rock, and non-linearly positively correlated with the freezing temperature. Controlling the freezing temperature and initial porosity of rock can significantly reduce the degree of deterioration in the stability of dangerous rock mass.

Estimating dynamic compressive strength of rock subjected to freeze-thaw weathering by data-driven models and non-destructive rock properties

ABSTRACT The dynamic compressive strength (DCS) of frozen-thawed rock influences the stability of rock mass in cold regions, especially when rock masses are possibly disturbed by dynamic loads. Laboratory freeze-thaw weathering treatment is usually time-consuming, and the dynamic strength test is destructive. Therefore, this paper attempts to quickly predict the DCS of frozen-thawed sandstones using data-driven methods, non-destructive rock properties, and basic environmental parameters. The sparrow search algorithm (SSA), gorilla troops optimiser, and dung beetle optimiser were chosen to develop two hyperparameters in the random forest (RF). The classic RF, back propagation neural network, and support vector regression models were taken as the control group. These six models were developed to predict the DCS. Their prediction results were compared. Finally, the sensitivity analysis was carried out to assess the significance of all input variables. The results indicate that the SSA – RF model yields the best prediction result, and three optimised models have better performance than single machine-learning models. Strain rate, dry density, and wave velocity are found to be the three most important parameters in DCS prediction, which further indicates that there is also a strong correlation between the characteristic impedance of the rock and the DCS.

A precise modeling method of three‐dimensional discrete fracture network based on rectangular joint model

AbstractAs the weak structure and main seepage channel of rock mass, the discrete fracture network (DFN) has a significant influence on the physical and mechanical properties of rock mass. In this paper, based on the rectangular joint model, two algorithms are proposed to realize the precise modeling of the DFN in any polyhedron. First, the equivalence of DFN modeling with the rectangular joint model and the elliptical joint model is theoretically deduced. On this basis, the chain splicing algorithm suitable for the rectangular joint model is proposed, which can generate the DFN model according to the actual distribution of joint spacing, trace, and bridge length. Aiming at the concave boundary in the model domain, the segmentation‐stitching operation is put forward, and the continuity of the joint when passing through the interface is ensured. Finally, the effectiveness of the above algorithm is verified by several simple cases, and the relationship between the size of the geometrical representative elementary volume (REV) and discontinuity parameters is discussed. The results show that when the mean values of the joint distribution parameters are the same, even if the distribution rules are different, the number of joints generated in the same spatial region is approximately equal. The size of the geometrical REV is independent of the spacing and bridge length and is approximately three to five times the mean of the trace. The relevant algorithms and conclusions in this paper can be used as the basis for subsequent rock mass stability analysis and seepage calculation.

The Mechanism of Mineral Dissolution on the Development of Red-Bed Landslides in the Wudongde Reservoir Region

The Wudongde reservoir region exhibits a notable prevalence of landslides within the red-bed reservoir stratum. The red bed is a clastic sedimentary rock layer dominated by red continental deposits. It is mainly composed of sandstone, mudstone, and siltstone. The lithology is diverse and uneven. In this study, we delve into the impact of mineral dissolution on the development of red-bed landslides in the reservoir area by utilizing the Xiaochatou landslide as a representative case study. Considering the inherent susceptibility of red-bed formations to erosion, collapse, and softening when exposed to water, an investigation was conducted to examine the consequences of mineral dissolution on landslides occurring in these areas. We conducted a mineral analysis and an identification of rock samples from the Xiaochatou landslide site, revealing alternating layers of sandstone and mudstone. Sandstone and conglomerate specimens were immersed in deionized water, and advanced techniques such as scanning electron microscopy (SEM), ion chromatography (IC), and inductively coupled plasma (ICP) analysis were used to examine the effects of water immersion. We also employed the hydrogeochemical simulation software PHREEQC to understand the dissolution mechanism of gypsum during soaking. Our findings reveal that sandstone and conglomerates harbor a notable quantity of gypsum, which readily dissolves in water. Prolonged immersion leads to erosion cavities within the sandstone, thereby augmenting its permeability. The concentration of SO42− ions in the soaking solution emerges as the highest, followed by Ca2+ and Na+. The notable significance is the dissolution of gypsum, whose intricate mechanism is contingent upon diverse environmental conditions. Variations in ion concentration profoundly influence the saturation index (SI) value, with the pH value playing a crucial role in shifting the reaction equilibrium. Regarding the deformation mode of the landslide, it manifests as a combination of sliding compression and tension cracking. The fracture surface of the landslide assumes a step-like configuration. As the deformation progresses, the mudstone layer takes control over the sliding process, causing the sandstone to develop internal narrow-top and wide-bottom cracks, which propagate upward until the stability of the slope rock mass is compromised, resulting in its rupture. In this manuscript, we delve into the dissolution traits of red-bed soft rock in the Wudongde reservoir area, using a landslide case as a reference. We simulate this rock’s dissolution under environmental water influences, examining its interaction with diverse water types through rigorous experiments and simulations. This study’s importance lies in its potential to shed light on the crucial engineering characteristics of red-bed soft rock.

Numerical Investigation for the Effect of Joint Persistence on Rock Slope Stability Using a Lattice Spring-Based Synthetic Rock Mass Model

This study underscores the profound influence of rock joints, both persistent and non-persistent with rock bridges, on the stability and behavior of rock masses—a critical consideration for sustainable engineering and natural structures, especially in rock slope stability. Leveraging the lattice spring-based synthetic rock mass (LS-SRM) modeling approach, this research aims to understand the impact of persistent and non-persistent joint parameters on rock slope stability. The Slope Model, a Synthetic Rock Mass (SRM) approach-based code, is used to investigate the joint parameters such as dip angle, spacing, rock bridge length, and trace overlapping. The results show that the mobilizing zones in slopes with non-persistent joints were smaller and shallower compared to slopes with fully persistent joints. The joint dip angle was found to heavily influence the failure mode in rock slopes with non-coplanar rock bridges. Shallow joint dip angles led to tensile failures, whereas steeper joint dip angles resulted in shear-tensile failures. Slopes with wider joint spacings exhibited deeper failure zones and a higher factor of safety, while longer rock bridge lengths enhanced slope stability and led to lower failure zones. The overlapping of joint traces has no apparent impact on slope stability and failure mechanism. This comprehensive analysis contributes valuable insights into sustainable rock engineering practices and the design of resilient structures in natural environments.

Mechanical responses and stress distribution of rock-like specimen containing a spherical defect under uniaxial compression

Spherical defect is often seen in rock material, and the size and position of the defect seriously affect the stability of rock mass. In this study, a locating device for preparing spherical defects was designed. A series of uniaxial compression tests and simulations were conducted to study the influence of the defect size and location on the mechanical, micro-cracking, failure characteristics, and stress distribution of the sample containing a spherical defect. The results indicate that the uniaxial compressive strength (UCS) and Young’s modulus both decrease with increasing defect size. The location variation of the same-sized spherical defect has a significant influence on the UCS but hardly influences Young’s modulus. The closer the defect is to the boundary, the lower the UCS is. The acoustic emission (AE) monitoring results indicate that the damage point aggregates and expands around the defect, forming a failure surface in the later stage of specimen loading. The failure modes of the specimen are divided into tensile failure, shear failure, and tensile-shear failure. Laboratory test results match well with the numerical simulation results obtained through the PFC3D program. The numerical stress monitoring results indicate that there is significant stress concentration near the defect. And the concentrated stress at the defect boundary drops before the the model reaches its peak stress, which can be used to predict the failure of the defected rock mass. The research findings could supply a certain reference for stability analysis and control of rock material containing similar types of defects.

DEM-based analysis of water inrush process of underground engineering face with intermittent joints in karst region


 
 
 
 Water inrush disaster of karst tunnel often lead to significant economic losses and serious casualties, which is an urgent engineering roadblock to be solved in the construction of tunnel in karst area. In this paper, three-dimensional discrete element method considering fluid-solid coupling effect and structural characteristics of water-mud resistant rock mass is adopted to systematically study the evolution law of displacement field and seepage field of intermittent joint type water-mud resistant rock mass of tunnel face and its water inrush critical characteristics during the process of sequential excavation of karst tunnel close to the frontal high-pressure water-rich karst cavity. The results show that: With the tunnel face gradually approaching the front-concealed high-pressure water-rich karst cavity, the stability of water-mud resistant rock mass is increasingly affected by high-pressure karst water, and karst water pressure gradually becomes the main control factor. The closer the tunnel face is to the front-concealed high-pressure water-rich karst cavity, the greater the extrusion displacement of karst tunnel face and its increase amplitude, the higher damage degree of water-mud resistant rock mass of face. With the advance of tunnel excavation, the intermittent cracks in the water-mud resistant rock mass of face gradually connect and form a stable hydraulic connection. The flow velocity and seepage pressure of karst water rise significantly at the moment of overall instability of face and the formation of water inrush channel, showing obvious precursor characteristics. The research achievements provide a reference for early warning and prevention and control of water inrush disaster of karst tunnel face.
 
 
 

Exploring the potential of smartphone photogrammetry for field measurement of joint roughness

Rock joints play a crucial role in determining the behavior and stability of rock masses. Despite their importance, there is a lack of devices for easily and accurately measuring joint roughness in field conditions. This research explores the potential of smartphone photogrammetry as a solution for large-scale joint roughness mapping. A large natural rock joint from an open pit mine in Yunnan, China, served as the specimen. Both two dimensional (2D) and three dimensional (3D) roughness parameters were utilized to characterize joint roughness. The study found that the accuracy of smartphone photogrammetry is acceptable when operated with guidance. The measurement errors were 1.7%, 4.7%, 7.8%, and 6.1% for the Z1, Ai, Z2, and PSD-D roughness parameters, respectively. In conclusion, this study demonstrates that using smartphone photogrammetry as a portable and cost-effective method for measuring joint roughness in field conditions has the potential to significantly benefit engineering geological investigations.

A Novel Deformation Analytical Solution and Constitutive Model for Fractured Rock Masses

In order to study the deformation and stability of a fractured rock mass, existing research suggests that fracture deformation, which is usually obtained by evaluating the equivalent deformation modulus, dominates the deformation of a fractured rock mass. However, this parameter is difficult to obtain in theory and practice, which limits the application of rock mass deformation analysis methods. In order to calculate the deformation of fractured rock masses, the mass of the rock is regarded as a sponge-like material, and it is assumed that the deformation of a water-saturated fractured rock mass under external force load is approximately equal to the net flow of fracture water. Based on this assumption, firstly, the relationship between rock mass deformation and fracture flow is studied through a single-fracture rock mass model. The hydraulic properties of the fracture are characterized by the permeability coefficient, and the fracture deformation in the rock mass is equivalent to the fracture flow. The fracture deformation calculation formula is derived from the fracture hydraulics calculation formula, and this is compared with the measured data. The rationality of the calculation formula was verified. On this basis, the calculation formula for rock mass deformation, including multiple groups of fracture surfaces, is proposed, and the stress-strain constitutive relationship of a complex rock mass is established. The correctness of the calculation method was verified by comparing it with other theoretical calculation results.

Implementing close-range remote surveys for the digitally supported rock mass stability analysis

A digitally supported rock mass stability analysis is herein presented for the implementation of non-contact close-range remote surveys against landslide hazard in sheltered environments, such as Marine Protected Areas. An airborne photogrammetric survey was carried out at the Lachea Islet (Sicily), which is part of a small volcanic archipelago with a peculiar geological history, hosting unstable rock masses belonging to the first eruptive stages of Mount Etna. The digital model built for this study allowed extracting the discontinuity orientation data, performing a kinematic analysis along the recognized discontinuity clusters, and quantifying the volume of potentially unstable boulders. Laboratory tests allowed achieving the main physical-mechanical parameters of involved rocks, which were used to perform a quantitative stability analysis under both deterministic and probabilistic approaches, based on the extracted digital data. Finally, a 3D rockfall trajectory analysis was carried out by simulating statistical rockfall path on the digital model built for the study. Results highlighted a poor stability in the studied islet sector, which is worsened by taking into account the high seismicity that characterizes the area. The main tourist elements of the islet must be considered at risk, as they are located within the simulated trajectories of the landslides. The methodological protocol carried out for this study, based on conventional procedures but using remotely derived and digitally processed information, can be considered a useful tool in the frame of nature-based solutions for hydro-meteorological hazards connected to climatic changes.

A 2D energy-conserving contact model for the combined finite-discrete element method (FDEM)

The combined finite-discrete element method (FDEM) has been widely used to simulate the progressive fracturing of brittle materials from continua to discontinua. However, due to the relatively simple definition of contact potential, the contact interaction algorithm utilized in the original FDEM is element size-dependent. As a result, a non-smooth and fluctuating contact force direction can occur when a contact point moves from one area of an element to another, and the contact interaction algorithm may fail to accurately evaluate the tangential contact force. To circumvent these limitations, we adopt an energy-conserving model developed recently to calculate the contact force in 2D FDEM, in which the magnitude and direction of contact force can be directly determined by geometrical features. A suite of numerical benchmarks is conducted to validate the effectiveness and robustness of the proposed method for contact interaction processing between discrete bodies, and the advantages of the proposed method are also demonstrated. As an application, two typical examples are performed to systematically explore the potential of the proposed approach for evaluating the stability of rock mass and rock slope. The proposed contact model may help enhance the applicability and accuracy of FDEM for rock fracturing simulation.

3D digital analysis for geo-structural monitoring and virtual documentation of the saint Michael cave in Minervino Murge, Bari (Italy)

The presence of caves and, in general, of subterranean landforms produced by karst processes is among the most typical geological features of the territory of Apulia, southern Italy. Several examples can be counted throughout the region, especially in the Murge area. The need to perform geo-structural analysis aimed at evaluating the stability conditions is especially significant for those caves open to the public. In the Murge area, the cave of Saint Michael at Minervino Murge is among the most famous. Beside the religious and historical interests, the whole cave needed a detailed analysis of the rock mass stability, both for protection of its architectural and archaeological values, and for safeguard of the visitors. At these aims, we performed various digital surveys by means of laser scanners and drones equipped with high-resolution cameras, in order to understand the geological setting of the entire cave and to identify unstable blocks.

Delineating homogeneous domains of fractured rocks using topological manifolds and deep learning

Determining homogeneous domains statistically is helpful for engineering geological modeling and rock mass stability evaluation. In this text, a technique that can integrate lithology, geotechnical and structural information is proposed to delineate homogeneous domains. This technique is then applied to a high and steep slope along a road. First, geological and geotechnical domains were described based on lithology, faults, and shear zones. Next, topological manifolds were used to eliminate the incompatibility between orientations and other parameters (i.e. trace length and roughness) so that the data concerning various properties of each discontinuity can be matched and characterized in the same Euclidean space. Thus, the influence of implicit combined effect in between parameter sequences on the homogeneous domains could be considered. Deep learning technique was employed to quantify abstract features of the characterization images of discontinuity properties, and to assess the similarity of rock mass structures. The results show that the technique can effectively distinguish structural variations and outperform conventional methods. It can handle multisource engineering geological information and multiple discontinuity parameters. This technique can also minimize the interference of human factors and delineate homogeneous domains based on orientations or multi-parameter with arbitrary distributions to satisfy different engineering requirements.

Determination of the minimum number of specimens required for laboratory testing of the shear strength of rock joints

The shear strength of rock joints is critical in determining rock mass stability. However, the reliability of shear strength determined by testing only a few rock joint specimens is questionable due to the inherent heterogeneity of rock joints. This study aimed to investigate shear strength heterogeneity and determine the required minimum number of specimens (RMN) for obtaining reliable values. Specimens were extracted from various positions within a natural rock joint, and their roughness and shear strength variations were examined. A significant correlation between shear strength heterogeneity and joint roughness was found. Monte Carlo simulation was used to simulate engineering practice sampling procedures. The effectiveness of specimen numbers under different normal stresses was evaluated using relative error and corresponding probability. Results showed that the RMN depended on acceptable relative error, specified probability, and normal stress. For acceptable relative errors of 5% to 10% and corresponding probabilities of 85–95%, the RMN exhibits a decrease with an increase in normal stress, while it demonstrates an increase when considering acceptable relative error and probability requirements. A smaller RMN was needed under high normal stress for the same rock joint, demonstrating a decrease in shear strength heterogeneity with increasing normal stress. However, commonly used specimen numbers like 3 or 5 did not always align with the recommended RMN. The results also emphasize the significance of considering rock material properties when determining the RMN. The proposed method provides a quantitative and reliable tool for determining the RMN, considering the heterogeneity of rock joints.

Computer information technology-based green excavation of tunnels in complex strata and technical decision of deformation control

Abstract Currently, information processing of tunnel engineering has mainly adopted conventional mathematical statistics-based methods. Moreover, some nonlinear processing methods are implemented to derive more insights, even though the degree of research is not deep enough. In the research, the rock mechanics test is carried out by drilling a method and taking samples in situ according to the construction technology of tunnels in complex geological conditions and implementing computer information-based methods. Also, rock mechanics tests are carried out in the excavation test area of the flat tunnel. Based on the tests using physical properties, such as deformation, tensile, uniaxial compression, triaxial compression, and longitudinal wave velocity, the physical and mechanical characteristics of the surrounding rock in the tunnel area are comprehensively evaluated, and the stability of the tunnel rock mass is assessed to devise convenient conditions for the subsequent research of the complex geological tunnels based on green excavation. The particle density of sandy mudstone, the bulk density, the porosity, and the natural water content are represented by 2.67 ± 0.61 g/cm3, 2.56 ± 1.42 g/cm3, 7.45%, and 2.86%, respectively, in terms of physical characteristics. These indicate that the sandy mudstone structure is relatively loose, with relatively large pores, micro-fractures, and a high degree of natural water content. The representative deformation test curve of the rock block shows that the ratio of deformation modulus to the compressive strength of the rock block is 650 on average, and Poisson’s ratio ranges from 0.21 to 0.38. These show that sandy mudstone has deformation properties after compression. The tensile strength of sandy mudstone, the shear strength, and c are represented by 1.25 ± 0.23 MPa, f = 1.32, and = 2.35 MPa, respectively. The stated test results can provide a scientific basis for selecting engineering design and its construction parameters in similar areas. In addition, the measurement results show that the surrounding rock will gradually increase, and the tunnel space will gradually become shorter with the increase of buried depth when the gravity stress field occurs. The linear elastic displacement of soft rock is smaller than that of elastic–plastic analysis, and the deeper the tunnel is buried, the larger the displacement difference would be. Therefore, establishing a stable and orderly monitoring and detection system could fully understand the intrinsic law between surrounding rock stress release and surrounding rock pressure and obtain accurate monitoring and measured data to evaluate the grading management standard of a tunnel at the ultimate displacement. In a word, this research provides a feasible idea to study the decision process of green excavation and deformation control technology of tunnels in complex strata.

Study on the Dynamic Stability and Spectral Characteristics of a Toppling Dangerous Rock Mass under Seismic Excitation

To evaluate the dynamic stability of dangerous rock masses under seismic excitation more reasonably, a mass viscoelasticity model was adopted to simulate the two main controlling surfaces of a toppling dangerous rock mass. Based on the principles of structural dynamics, a dynamic response analysis model and motion equations were established for toppling dangerous rock masses. The Newmark-β method was utilized to establish a calculation method for the dynamic stability coefficient of a toppling dangerous rock mass. This method was applied to the WY2 dangerous rock mass developed in a steep cliff zone in Luoyi Village, and the dynamic stability coefficient time history was calculated. Subsequently, the acceleration response signals of the dangerous rock mass in different directions were analyzed using wavelet packet transform. The results show that the sum of the energy proportions of the first to third frequency bands in the n1 and s2 directions exceeded 95%. This suggests that the n1 and s2 directions of the WY2 dangerous rock mass suffered the initial damage under bidirectional seismic actions. Finally, the marginal spectra variations of the acceleration response signals in different directions were analyzed based on the HHT. The results show that the seismic energy in the n1 and s2 directions of the dangerous rock mass was found to be the most significant under seismic loading, indicating that the rock mass experienced the most severe damage along these two directions. This reveals that the failure mode of the dangerous rock mass is inclined toppling, consistent with the results of wavelet packet analysis.

An Enhanced RIME Optimizer with Horizontal and Vertical Crossover for Discriminating Microseismic and Blasting Signals in Deep Mines.

Real-time monitoring of rock stability during the mining process is critical. This paper first proposed a RIME algorithm (CCRIME) based on vertical and horizontal crossover search strategies to improve the quality of the solutions obtained by the RIME algorithm and further enhance its search capabilities. Then, by constructing a binary version of CCRIME, the key parameters of FKNN were optimized using a binary conversion method. Finally, a discrete CCRIME-based BCCRIME was developed, which uses an S-shaped function transformation approach to address the feature selection issue by converting the search result into a real number that can only be zero or one. The performance of CCRIME was examined in this study from various perspectives, utilizing 30 benchmark functions from IEEE CEC2017. Basic algorithm comparison tests and sophisticated variant algorithm comparison experiments were also carried out. In addition, this paper also used collected microseismic and blasting data for classification prediction to verify the ability of the BCCRIME-FKNN model to process real data. This paper provides new ideas and methods for real-time monitoring of rock mass stability during deep well mineral resource mining.

Diffusion mechanism of solid waste product utilization pulping and fracture network grouting

China's commitment to achieving carbon neutrality by 2060 has highlighted the pressing need for effective solid waste recycling. The reinforcement of mine water hazards has led to increased production costs due to the significant consumption of slurry. This study focuses on the development of mine crack grouting sealing and reinforcement materials using mine slag powder and waste rubber particles, along with cement reinforcement and alkali slag stimulation. The impact of water-cement ratio, rubber particle content, and slag mass fraction on slurry performance is analyzed, and the optimal values for these parameters are determined. A two-phase seepage simulation using COMSOL is conducted to investigate the flow and diffusion of slurry in a fractured network rock mass. The findings indicate that a water-cement ratio of 0.7 and a slag and rubber particle mass fraction of 70% and 30% respectively result in the slurry with favorable fluidity and mechanical properties. The diffusion pattern of the slurry in the fractured rock mass is characterized by vortex-like diffusion. The stability of the fractured rock mass is influenced by both grouting time and pressure, which should be carefully considered during on-site operations. This research has significant implications for the integrated prevention and control of mine water hazards and the large-scale utilization of high-dosage solid waste, thereby ensuring production safety, ecological environment safety, and sustainable development in mining areas.

Study on the evolution of microscopic pore structure of sandstone under freeze-thaw cycles

In cold regions, rocks undergo freeze-thaw (FT) cycles, leading to damage caused by the evolution of pores. This phenomenon significantly influences the rock mass stability. A comprehensive and accurate study of FT damage evolution in rock pores was conducted in this research. The microscopic morphology and pore size distribution (PSD) of sandstone samples subjected to FT cycles were measured utilizing scanning electron microscopy (SEM), nitrogen adsorption/desorption (NAD) method, and mercury intrusion porosimetry (MIP). Multifractal analysis was applied to explore the evolution of PSD. A gradual weakening of the cementation degree of mineral particles under FT cycles was observed based on SEM results, resulting in the formation of larger connected pores. A decrease of 37.3% in the number of pores was recorded after 40 cycles. According to NAD and MIP results, a greater effect on mesopores and large pores were found, with a minimal impact on micropores and transition pores. Multifractal characteristics were identified in the PSD of sandstone samples with different FT cycles. Initially, the connectivity of sandstone pores was found to decrease and then increase gradually after 10 cycles. Uneven PSD was observed to concentrate stress among pores, exacerbate frost heave damage, and drive the PSD towards a more uniform direction. It was concluded that the non-uniform degree of PSD does not consistently increase with FT cycles. The findings of this study can serve as a reference for further research on rock FT damage mechanisms.