Probabilistic analysis of tunnel face seismic stability in layered rock masses using Polynomial Chaos Kriging metamodel

Face stability is an essential issue in tunnel design and construction.Layered rock masses are typical and ubiquitous;uncertainties in rock properties always exist.In view of this,a comprehensive method,which combines the Upper bound Limit analysis of Tunnel face stability,the Polynomial Chaos Kriging,the Monte-Carlo Simulation and Analysis of Covariance method(ULT-PCK-MA),is proposed to investigate the seismic stability of tunnel faces.A two-dimensional analytical model of ULT is developed to evaluate the virtual support force based on the upper bound limit analysis.An efficient probabilistic analysis method PCK-MA based on the adaptive Polynomial Chaos Kriging metamodel is then implemented to investigate the parameter uncertainty effects.Ten input parameters,including geological strength indices,uniaxial compressive strengths and constants for three rock formations,and the horizontal seismic coefficients,are treated as random variables.The effects of these parameter uncertainties on the failure probability and sensitivity indices are discussed.In addition,the effects of weak layer position,the middle layer thickness and quality,the tunnel diameter,the parameters correlation,and the seismic loadings are investigated,respectively.The results show that the layer distributions significantly influence the tunnel face probabilistic stability,particularly when the weak rock is present in the bottom layer.The efficiency of the proposed ULT-PCK-MA is validated,which is expected to facilitate the engineering design and construction.

Evaluation of slope stability through rock mass classification and kinematic analysis of some major slopes along NH-1A from Ramban to Banihal, North Western Himalayas

The network of Himalayan roadways and highways connects some remote regions of valleys or hill slopes,which is vital for India’s socio-economic growth.Due to natural and artificial factors,frequency of slope instabilities along the networks has been increasing over last few decades.Assessment of stability of natural and artificial slopes due to construction of these connecting road networks is significant in safely executing these roads throughout the year.Several rock mass classification methods are generally used to assess the strength and deformability of rock mass.This study assesses slope stability along the NH-1A of Ramban district of North Western Himalayas.Various structurally and non-structurally controlled rock mass classification systems have been applied to assess the stability conditions of 14 slopes.For evaluating the stability of these slopes,kinematic analysis was performed along with geological strength index(GSI),rock mass rating(RMR),continuous slope mass rating(CoSMR),slope mass rating(SMR),and Q-slope in the present study.The SMR gives three slopes as completely unstable while CoSMR suggests four slopes as completely unstable.The stability of all slopes was also analyzed using a design chart under dynamic and static conditions by slope stability rating(SSR)for the factor of safety(FoS)of 1.2 and 1 respectively.Q-slope with probability of failure(PoF)1%gives two slopes as stable slopes.Stable slope angle has been determined based on the Q-slope safe angle equation and SSR design chart based on the FoS.The value ranges given by different empirical classifications were RMR(37-74),GSI(27.3-58.5),SMR(11-59),and CoSMR(3.39-74.56).Good relationship was found among RMR&SSR and RMR&GSI with correlation coefficient(R 2)value of 0.815 and 0.6866,respectively.Lastly,a comparative stability of all these slopes based on the above classification has been performed to identify the most critical slope along this road.

Borehole stability in naturally fractured rocks with drilling mud intrusion and associated fracture strength weakening:A coupled DFN-DEM approach

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(P w)and pore pressure(P p)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.

Heterogeneous information phase space reconstruction and stability prediction of filling body–surrounding rock combination

Traditional research believes that the filling body can effectively control stress concentration while ignoring the problems of unknown stability and the complex and changeable stress distribution of the filling body–surrounding rock combination under high-stress conditions.Current monitoring data processing methods cannot fully consider the complexity of monitoring objects,the diversity of monitoring methods,and the dynamics of monitoring data.To solve this problem,this paper proposes a phase space reconstruction and stability prediction method to process heterogeneous information of backfill–surrounding rock combinations.The three-dimensional monitoring system of a large-area filling body–surrounding rock combination in Longshou Mine was constructed by using drilling stress,multipoint displacement meter,and inclinometer.Varied information,such as the stress and displacement of the filling body–surrounding rock combination,was continuously obtained.Combined with the average mutual information method and the false nearest neighbor point method,the phase space of the heterogeneous information of the filling body–surrounding rock combination was then constructed.In this paper,the distance between the phase point and its nearest point was used as the index evaluation distance to evaluate the stability of the filling body–surrounding rock combination.The evaluated distances(ED)revealed a high sensitivity to the stability of the filling body–surrounding rock combination.The new method was then applied to calculate the time series of historically ED for 12 measuring points located at Longshou Mine.The moments of mutation in these time series were at least 3 months ahead of the roadway return dates.In the ED prediction experiments,the autoregressive integrated moving average model showed a higher prediction accuracy than the deep learning models(long short-term memory and Transformer).Furthermore,the root-mean-square error distribution of the prediction results peaked at 0.26,thus outperforming the no-prediction method in 70%of the cases.

A stability evaluation method for deep-seated toppling in the upper Lancang river,Southwestern China

Deep-seated toppling in the upper reaches of the Lancang River,southwest China involves deformations exceeding 100 m in depth.The slope deformation is initiated by river downcutting and evolves distinctive characteristics with a depth of river incision.In this study,we propose a system for evaluating the stability of deep-seated toppled slopes in different evolutionary stages.This system contains identification criteria for each evolutionary stage and provides the corresponding stability evaluation methods.Based on the mechanical and kinematic analysis of slope blocks,the specific stage of slope movement can be identified in the field through outcrop mapping,in situ tests,surface displacement monitoring,and adit and borehole explorations.The stability evaluation methods are established based on the limiting equilibrium theory and the strain compatibility between the undisturbed zone and the toppled zone.Finally,several sample slopes in different evolution stages have been investigated to verify the applicability and accuracy of the proposed stability evaluation system.The results indicate that intense tectonic activity and rapid river incision lead to a maximum principal stress ratio exceeding 10 near the slope surface,thus triggering widespread toppling deformations along the river valley.When considering the losses of joint cohesion during the further rotation process,the safety factor of the slope drops by 7%e28%.The self-stabilization of toppling deformation can be recognized by the layer symmetry configuration after the free rotation of the deflected layers.Intensely toppled rock blocks mainly suffer sliding failures beyond the layer symmetry condition.The factor of safety of the K73 rockslide decreased from 1.17 to 0.87 by considering the development of the potential sliding surface and the toesaturated zone.

Automatic identification of rock discontinuity and stability analysis of tunnel rock blocks using terrestrial laser scanning

Local geometric information and discontinuity features are key aspects of the analysis of the evolution and failure mechanisms of unstable rock blocks in rock tunnels.This study demonstrates the integration of terrestrial laser scanning(TLS)with distinct element method for rock mass characterization and stability analysis in tunnels.TLS records detailed geometric information of the surrounding rock mass by scanning and collecting the positions of millions of rock surface points without contact.By conducting a fuzzy K-means method,a discontinuity automatic identification algorithm was developed,and a method for obtaining the geometric parameters of discontinuities was proposed.This method permits the user to visually identify each discontinuity and acquire its spatial distribution features(e.g.occurrences,spac-ings,trace lengths)in great detail.Compared with hand mapping in conventional geotechnical surveys,the geometric information of discontinuities obtained by this approach is more accurate and the iden-tification is more efficient.Then,a discrete fracture network with the same statistical characteristics as the actual discontinuities was generated with the distinct element method,and a representative nu-merical model of the jointed surrounding rock mass was established.By means of numerical simulation,potential unstable rock blocks were assessed,and failure mechanisms were analyzed.This method was applied to detection and assessment of unstable rock blocks in the spillway and sand flushing tunnel of the Hongshiyan hydropower project after a collapse.The results show that the noncontact detection of blocks was more labor-saving with lower safety risks compared with manual surveys,and the stability assessment was more reliable since the numerical model built by this method was more consistent with the distribution characteristics of actual joints.This study can provide a reference for geological survey and unstable rock block hazard mitigation in tunnels subjected to complex geology and active rockfalls.

Stability prediction of hard rock pillar using support vector machine optimized by three metaheuristic algorithms

Hard rock pillar is one of the important structures in engineering design and excavation in underground mines.Accurate and convenient prediction of pillar stability is of great significance for underground space safety.This paper aims to develop hybrid support vector machine(SVM)models improved by three metaheuristic algorithms known as grey wolf optimizer(GWO),whale optimization algorithm(WOA)and sparrow search algorithm(SSA)for predicting the hard rock pillar stability.An integrated dataset containing 306 hard rock pillars was established to generate hybrid SVM models.Five parameters including pillar height,pillar width,ratio of pillar width to height,uniaxial compressive strength and pillar stress were set as input parameters.Two global indices,three local indices and the receiver operating characteristic(ROC)curve with the area under the ROC curve(AUC)were utilized to evaluate all hybrid models’performance.The results confirmed that the SSA-SVM model is the best prediction model with the highest values of all global indices and local indices.Nevertheless,the performance of the SSASVM model for predicting the unstable pillar(AUC:0.899)is not as good as those for stable(AUC:0.975)and failed pillars(AUC:0.990).To verify the effectiveness of the proposed models,5 field cases were investigated in a metal mine and other 5 cases were collected from several published works.The validation results indicated that the SSA-SVM model obtained a considerable accuracy,which means that the combination of SVM and metaheuristic algorithms is a feasible approach to predict the pillar stability.

Classification and rating of disintegrated dolomite strata for slope stability analysis

Although disintegrated dolomite,widely distributed across the globe,has conventionally been a focus of research in underground engineering,the issue of slope stability issues in disintegrated dolomite strata is gaining increasing prominence.This is primarily due to their unique properties,including low strength and loose structure.Current methods for evaluating slope stability,such as basic quality(BQ)and slope stability probability classification(SSPC),do not adequately account for the poor integrity and structural fragmentation characteristic of disintegrated dolomite.To address this challenge,an analysis of the applicability of the limit equilibrium method(LEM),BQ,and SSPC methods was conducted on eight disintegrated dolomite slopes located in Baoshan,Southwest China.However,conflicting results were obtained.Therefore,this paper introduces a novel method,SMRDDS,to provide rapid and accurate assessment of disintegrated dolomite slope stability.This method incorporates parameters such as disintegrated grade,joint state,groundwater conditions,and excavation methods.The findings reveal that six slopes exhibit stability,while two are considered partially unstable.Notably,the proposed method demonstrates a closer match with the actual conditions and is more time-efficient compared with the BQ and SSPC methods.However,due to the limited research on disintegrated dolomite slopes,the results of the SMRDDS method tend to be conservative as a safety precaution.In conclusion,the SMRDDS method can quickly evaluate the current situation of disintegrated dolomite slopes in the field.This contributes significantly to disaster risk reduction for disintegrated dolomite slopes.

Effect of surface retaining elements on rock stability:laboratory investigation with sand powder 3D printing

This study aims to investigate the benefcial efects of surface retaining elements (SREs) on the mechanical behaviors of bolted rock and roadway stability. 3D printing (3DP) technology is utilized to create rock analogue prismatic specimens for conducting this investigation. Uniaxial compression tests with acoustic emission (AE) and digital image correlation techniques have been conducted on 3DP specimens bolted with diferent SREs. The results demonstrate that the strength and modulus of elasticity of the bolted specimens show a positive correlation with the area of the SRE;the AE characteristics of the bolted specimens are higher than those of the unbolted specimen, but they decrease with an increase in SRE area, thus further improving the integrity of the bolted specimens. The reinforcement efect of SREs on the surrounding rock of roadways is further analyzed using numerical modelling and feld test. The results provide a better understanding of the role of SREs in rock bolting and the optimization of rock bolting design. Furthermore, they verify the feasibility of 3DP for rock analogues in rock mechanics tests.

Probabilistic assessment of rock slopes stability using the response surface approach——A case study

This paper examines the stability condition of a jointed rock slope in the south western region of Saudi Arabia using deterministic and probabilistic approaches,under both dry and wet conditions.The study area is characterized by complex geology in rugged terrains.The stability analysis is carried out using the code FLAC3D to generate a 3-dimensional,ubiquitous joint model,to determine the influence of the dominant,unfavourable discontinuity orientation with respect to the slope face.The deterministic analysis is first implemented using the mean values of the selected random variables,namely the dip,dip direction and friction angle of the dominant discontinuity set,and the stability condition is assessed with a factor of safety based on the classical frictional joint constitutive model.A Box-Behnken design(BBD)approach is then adopted to create the surface response function as a second order polynomial for the factor of safety.To do so,fifteen FLAC3D models are generated in accordance with the BBD.Based on this,10,000 simulations of different slope realizations are carried out using Monte-Carlo simulation technique,and the probability of unsatisfactory of performance of the rock slope is assessed.It is shown that the probabilistic approach provides more insight and confidence in the stability condition of the rock slope,both under dry and steady state heavy rainfall conditions.A discussion is presented on the significance of accepting lower safety factors when heavy rainfall conditions are encountered.

Application of rock mass classification systems to rock slope stability assessment:A case study

The stability of rock slopes is considered crucial to public safety in highways passing through rock cuts, as well as to personnel and equipment safety in open pit mines. Slope instability and failures occur due to many factors such as adverse slope geometries, geological discontinuities, weak or weathered slope materials as well as severe weather conditions. External loads like heavy precipitation and seismicity could play a significant role in slope failure. In this paper, several rock mass classification systems developed for rock slope stability assessment are evaluated against known rock slope conditions in a region of Saudi Arabia, where slopes located in rugged terrains with complex geometry serve as highway road cuts. Selected empirical methods have been applied to 22 rock cuts that are selected based on their failure mechanisms and slope materials. The stability conditions are identified, and the results of each rock slope classification system are compared. The paper also highlights the limitations of the empirical classification methods used in the study and proposes future research directions.

Microseismic monitoring and numerical simulation on the stability of high-steep rock slopes in hydropower engineering

For high-steep slopes in hydropower engineering, damage can be induced or accumulated due to a seriesof human or natural activities, including excavation, dam construction, earthquake, rainstorm, rapid riseor drop of water level in the service lifetime of slopes. According to the concept that the progressivedamage （microseismicity） of rock slope is the essence of the precursor of slope instability, a microseismicmonitoring system for high-steep rock slopes is established. Positioning accuracy of the monitoringsystem is tested by fixed-position blasting method. Based on waveform and cluster analyses of microseismicevents recorded during test, the tempo-spatial distribution of microseismic events is analyzed.The deformation zone in the deep rock masses induced by the microseismic events is preliminarilydelimited. Based on the physical information measured by in situ microseismic monitoring, an evaluationmethod for the dynamic stability of rock slopes is proposed and preliminarily implemented bycombining microseismic monitoring and numerical modeling. Based on the rock mass damage modelobtained by back analysis of microseismic information, the rock mass elements within the microseismicdamage zone are automatically searched by finite element program. Then the stiffness and strengthreductions are performed on these damaged elements accordingly. Attempts are made to establish thecorrelation between microseismic event, strength deterioration and slope dynamic instability, so as toquantitatively evaluate the dynamic stability of slope. The case studies about two practical slopes indicatethat the proposed method can reflect the factor of safety of rock slope more objectively. Numericalanalysis can help to understand the characteristics and modes of the monitored microseismic events inrock slopes. Microseismic monitoring data and simulation results can be used to mutually modify thesensitive rock parameters and calibrate the model. Combination of microseismic monitoring and numericalsimulation provides a more objective basis for the numerical model and parameters and a solidmechanical foundation for the microseismic monitoring.

Deformation Mechanism and Stability of a Rocky Slope

A high slope is located on the side of the spillway at a hydropower station in Southwest China, which has some weak inter-layers inclining outwards. Parts of the slope show heavy weathering and unloading. There appeared deformation and tensile crack either on the surface or on the afteredge of the slope during excavation, and under a platform （elev. 488 m）, two levels of slopes collapsed on the downriver side. Based on the investigation in situ and the analysis of the geological structure, the conceptual model of deformation and failure mechanism was erected for this slope. Furthermore, the deformation characteristics were studied with FLAC^3D numerical simulation. Comprehensive analysis shows that the whole deformation of the slope is unloading rebound in certain depth scope and the whole body does not slide along any weak interlayer. In addition, two parts with prominent local deformation in the shallow layer of the slope show the models of ＂creep sfiding-tensile cracking＂ and ＂slidlng-tensile cracking＂, respectively. Based on the above analysis, the corresponding project of support and reinforcement is proposed to make the slope more stable.

3D Identification and Stability Analysis of Key Surface Blocks of Rock Slope

Complicated geological structures make it difficult to analyze the stability of rock slopes, such as faults, weak intercalated layers or joint fissures. Based on 3D geological modeling and surface block identifying methods, an integrated methodology framework was proposed and realized to analyze the stability of surface blocks in rock slopes. The surface blocks cut by geological structures, fissures or free faces could be identified subjected to the four principles of closure, completeness, uniqueness and validity. The factor of safety(FOS)of single key block was calculated by the limit equilibrium method. If there were two or more connected blocks, they were defined as a block-group. The FOS of a block-group was computed by the Sarma method. The proposed approach was applied to an actual rock slope of a hydropower project, and some possible instable blocks were demonstrated and analyzed visually. The obtained results on the key blocks or block-groups provide essential information for determining potential instable region of rock slopes and designing effective support scheme in advance.

Model test and numerical simulation on the dynamic stability of the bedding rock slope under frequent microseisms

Shake table testing was performed to investigate the dynamic stability of a mid-dip bedding rock slope under frequent earthquakes. Then, numerical modelling was established to further study the slope dynamic stability under purely microseisms and the influence of five factors, including seismic amplitude, slope height, slope angle, strata inclination and strata thickness, were considered. The experimental results show that the natural frequency of the slope decreases and damping ratio increases as the earthquake loading times increase. The dynamic strength reduction method is adopted for the stability evaluation of the bedding rock slope in numerical simulation, and the slope stability decreases with the increase of seismic amplitude, increase of slope height, reduction of strata thickness and increase of slope angle. The failure mode of a mid-dip bedding rock slope in the shaking table test is integral slipping along the bedding surface with dipping tensile cracks at the slope rear edge going through the bedding surfaces. In the numerical simulation, the long-term stability of a mid-dip bedding slope is worst under frequent microseisms and the slope is at risk of integral sliding instability, whereas the slope rock mass is more broken than shown in the shaking table test. The research results are of practical significance to better understand the formation mechanism of reservoir landslides and prevent future landslide disasters.

Assessment of rock slope stability by probabilistic-based Slope Stability Probability Classification method along highway cut slopes in Adilcevaz-Bitlis(Turkey)

Rock slope stability is of great concern along highway routes as stability problems on cut slopes may cause fatal events as well as loss of property.In rock slope engineering,stability evaluations are commonly performed by means of analytical or numerical analyses,principally considering the factor of safety concept.As a matter of fact,the probabilistic assessment of slope stability is progressively getting popularity due to difficulties in assigning the most appropriate values to design parameters in analytical or numerical methods.Additionally,the effect of heterogeneities in rock masses and discontinuities on the analysis results is minimized through the probabilistic concept.In this study,slope stability of high and steep sedimentary rock cut slopes along a state highway in AdilcevazBitlis(Turkey) was evaluated on the basis of probabilistic approach using the Slope Stability Probability Classification(SSPC) system.The probabilistic assessment indicates major slope stability problems because of discontinuity controlled and discontinuity orientation independent mass movements.Almost all studied cut slopes suffer from orientation-independent stability problems with very low stability probabilities.Additionally,the probability of planar and toppling failures issignificantly high with respect to the SSPC system.The stability problems along the investigated rock slopes were also verified by field reconnaissance.Remedial measures such as slope re-design and reinforcement at the studied locations should be taken to prevent hazardous events along the highway.On the other hand,the probabilistic approach may be a useful tool during rock slope engineering to overcome numerous uncertainties when probabilistic and analytic results are compared.

Seismic stability of jointed rock slopes under obliquely incident earthquake waves

Seismic stability of slopes has been traditionally analyzed with vertically propagated earthquake waves.However,for rock slopes,the earthquake waves might approach the outcrop still with a evidently oblique direction.To investigate the impact of obliquely incident earthquake excitations,the input method for SV and P waves with arbitrary incident angles is conducted,respectively,by adopting the equivalent nodal force method together with a viscous-spring boundary.Then,the input method is introduced within the framework of ABAQUS software and verified by a numerical example.Both SV and P waves input are considered herein for a 2 D jointed rock slope.For the jointed rock mass,the jointed material model in ABAQUS software is employed to simulate its behavior as a continuum.Results of the study show that the earthquake incident angles have significance on the seismic stability of jointed rock slopes.The larger the incident angle,the greater the risk of slope instability.Furthermore,the stability of the jointed rock slopes also is affected by wave types of earthquakes heavily.P waves induce weaker responses and SV waves are shown to be more critical.

Rock slope stability evaluation using kinematic and kinetic methods along the Kamyaran-Marivan road,west of Iran

Interest in rock slope stability in mountainous regions has increased greatly in recent years.This issue has become a topic of major interest for geoscientists and engineering professionals,as well as for private citizens and local administrators,in many parts of the world.This paper evaluates the stability of seven rock slopes along the KamyaranMarivan tourist road,Kurdistan province,Iran,using various methods.The two main reasons for performing this research were to determine whether different methods of stability analysis provide the same results,and to determine how different factors such as the presence of water,tension cracks,and seismic forces affect the stability of these rock slopes.Firstly,field investigations were performed to obtain the engineering characteristics of the rock masses,discontinuities,and intact rocks of the slopes.Secondly,laboratory tests were carried out on rock samples obtained from the slopes,to determine the engineering properties of the intact rocks.Then for each rock slope,the contour diagram of discontinuities and slope face was drawn in the Dips v.5.1 software environment,and the failure mechanism was determined based on the kinematic or stereographic method.Next,the factors of safety of the rock slopes were calculated using the limit equilibrium method,based on the failure mechanisms resulting from the kinematic method.The accuracy of the results obtained by these two methods was investigated using SWedge v.4.0 software.The results indicated that four rock slopes have a potential for plane,wedge,and toppling failure,and three others are stable.Also,it was found that the stability of the studied rock slopes decreases greatly in the presence of water,tension cracks,and seismic forces.

Numerical investigations of rock bridge effect on open pit slope stability

In this study,the effect of rock bridges on rock slope stability was investigated by incorporating nonpersistent joint networks in numerical models,and the critical profiles of an open pit mine were analysed.Parallel deterministic networks of infinite and finite lengths,ubiquitous joint network model and Veneziano joint network model were used in order to simulate the rock fractures.Materials were modelled based on the generalised Hoek-Brown and equivalent Mohr-Coulomb failure criteria.The parallel deterministic infinite and the ubiquitous joint network models produced lower safety factors.The introduction of rock bridges along discontinuity planes in the parallel deterministic network and Veneziano joint network models significantly contributed to the stability and strain distribution,which should be considered in stability analysis of rock mass in open pit by rock slope practitioners.The results show the significance of joints in hard rock behaviour and the joints should be included in order to attain practical and realistic simulations.

Geostructural stability assessment of cave using rock surface discontinuity extracted from terrestrial laser scanning point cloud

The use of terrestrial laser scanning（TLS） in the caves has been growing drastically over the last decade.However, TLS application to cave stability assessment has not received much attention of researchers.This study attempted to utilize rock surface orientations obtained from TLS point cloud collected along cave passages to（1） investigate the influence of rock geostructure on cave passage development, and（2）assess cave stability by determining areas susceptible to different failure types. The TLS point cloud was divided into six parts（Entry hall, Chamber, Main hall, Shaft 1, Shaft 2 and Shaft 3）, each representing different segments of the cave passages. Furthermore, the surface orientation information was extracted and grouped into surface discontinuity joint sets. The computed global mean and best-fit planes of the entire cave show that the outcrop dips 290° with a major north-south strike. But at individual level, the passages with dip angle between 26° and 80° are featured with dip direction of 75°-322°. Kinematic tests reveal the potential for various failure modes of rock slope. Our findings show that toppling is the dominant failure type accounting for high-risk rockfall in the cave, with probabilities of 75.26%, 43.07%and 24.82% in the Entry hall, Main hall and Shaft 2, respectively. Unlike Shaft 2 characterized by high risk of the three failure types（32.49%, 24.82% and 50%）, the chamber and Shaft 3 passages are not suffering from slope failure. The results also show that the characteristics of rock geostructure considerably influence the development of the cave passages, and four sections of the cave are susceptible to different slope failure types, at varying degrees of risk.