A methodology for damage evaluation of underground tunnels subjected to static loading using numerical modeling

We have proposed a methodology to assess the robustness of underground tunnels against potential failure.This involves developing vulnerability functions for various qualities of rock mass and static loading intensities.To account for these variations,we utilized a Monte Carlo Simulation(MCS)technique coupled with the finite difference code FLAC^(3D),to conduct two thousand seven hundred numerical simulations of a horseshoe tunnel located within a rock mass with different geological strength index system(GSIs)and subjected to different states of static loading.To quantify the severity of damage within the rock mass,we selected one stress-based(brittle shear ratio(BSR))and one strain-based failure criterion(plastic damage index(PDI)).Based on these criteria,we then developed fragility curves.Additionally,we used mathematical approximation techniques to produce vulnerability functions that relate the probabilities of various damage states to loading intensities for different quality classes of blocky rock mass.The results indicated that the fragility curves we obtained could accurately depict the evolution of the inner and outer shell damage around the tunnel.Therefore,we have provided engineers with a tool that can predict levels of damages associated with different failure mechanisms based on variations in rock mass quality and in situ stress state.Our method is a numerically developed,multi-variate approach that can aid engineers in making informed decisions about the robustness of underground tunnels.

Determination of uncertainties of geomechanical parameters of metamorphic rocks using petrographic analyses

Geomechanical parameters of intact metamorphic rocks determined from laboratory testing remain highly uncertain because of the great intrinsic variability associated with the degrees of metamorphism.The aim of this paper is to develop a proper methodology to analyze the uncertainties of geomechanical characteristics by focusing on three domains,i.e.data treatment process,schistosity angle,and mineralogy.First,the variabilities of the geomechanical laboratory data of Westwood Mine(Quebec,Canada)were examined statistically by applying different data treatment techniques,through which the most suitable outlier methods were selected for each parameter using multiple decision-making criteria and engineering judgment.Results indicated that some methods exhibited better performance in identifying the possible outliers,although several others were unsuccessful because of their limitation in large sample size.The well-known boxplot method might not be the best outlier method for most geomechanical parameters because its calculated confidence range was not acceptable according to engineering judgment.However,several approaches,including adjusted boxplot,2MADe,and 2SD,worked very well in the detection of true outliers.Also,the statistical tests indicate that the best-fitting probability distribution function for geomechanical intact parameters might not be the normal distribution,unlike what is assumed in most geomechanical studies.Moreover,the negative effects of schistosity angle on the uniaxial compressive strength(UCS)variabilities were reduced by excluding the samples within a specific angle range where the UCS data present the highest variation.Finally,a petrographic analysis was conducted to assess the associated uncertainties such that a logical link was found between the dispersion and the variabilities of hard and soft minerals.

Evaluation of the effect of geometrical parameters on stope probability of failure in the open stoping method using numerical modeling

Stress-induced failure is among the most common causes of instability in Canadian deep underground mines.Open stoping is the most widely practiced underground excavation method in these mines,and creates large stopes which are subjected to stress-induced failure.The probability of failure(POF)depends on many factors,of which the geometry of an open stope is especially important.In this study,a methodology is proposed to assess the effect of stope geometrical parameters on the POF,using numerical modelling.Different ranges for each input parameter are defined according to previous surveys on open stope geometry in a number of Canadian underground mines.A Monte-Carlo simulation technique is combined with the finite difference code FLAC3D,to generate model realizations containing stopes with different geometrical features.The probability of failure(POF)for different categories of stope geometry,is calculated by considering two modes of failure;relaxation-related gravity driven(tensile)failure and rock mass brittle failure.The individual and interactive effects of stope geometrical parameters on the POF,are analyzed using a general multi-level factorial design.Finally,mathematical optimization techniques are employed to estimate the most stable stope conditions,by determining the optimal ranges for each stope’s geometrical parameter.