Characterization of Spatial Variability in Rock Mass Mechanical Parameters for Slope Stability Assessment: A Comprehensive Case Study

The spatial variability in rock mass mechanical parameters critically affects slope stability assessments. This study investigated the southern slope of the Bayan Obo open-pit mine. A representative elementary volume (REV) with a side length of 14 m was determined through discrete fracture network (DFN) simulations. Based on the rock quality designation (RQD) data from 40 boreholes, a three-dimensional spatial distribution model of the RQD was constructed using Ordinary Kriging interpolation. The RQD values were converted into geological strength index (GSI) values through an empirical correlation, and the generalized Hoek–Brown criterion was applied to develop a spatially heterogeneous equivalent mechanical parameter field. Numerical simulations were performed using FLAC3D, with the slope stability evaluated using the point safety factor (PSF) method. For comparison, three homogeneous benchmark models based on the 5th, 25th, and 50th percentiles produced profile-scale safety factors of 0.96–1.92 and failed to replicate the observed failure geometry. By contrast, the heterogeneous model yielded safety factors of approximately 1.03–1.08 and accurately reproduced the mapped sliding surface. These findings demonstrate that incorporating spatial heterogeneity significantly improves the accuracy of slope stability assessments, providing a robust theoretical basis for targeted monitoring and reinforcement design.