As a calculation method based on the Galerkin variation,the numerical manifold method(NMM)adopts a double covering system,which can easily deal with discontinuous deformation problems and has a high calculation accura...As a calculation method based on the Galerkin variation,the numerical manifold method(NMM)adopts a double covering system,which can easily deal with discontinuous deformation problems and has a high calculation accuracy.Aiming at the thermo-mechanical(TM)coupling problem of fractured rock masses,this study uses the NMM to simulate the processes of crack initiation and propagation in a rock mass under the influence of temperature field,deduces related system equations,and proposes a penalty function method to deal with boundary conditions.Numerical examples are employed to confirm the effectiveness and high accuracy of this method.By the thermal stress analysis of a thick-walled cylinder(TWC),the simulation of cracking in the TWC under heating and cooling conditions,and the simulation of thermal cracking of the SwedishÄspöPillar Stability Experiment(APSE)rock column,the thermal stress,and TM coupling are obtained.The numerical simulation results are in good agreement with the test data and other numerical results,thus verifying the effectiveness of the NMM in dealing with thermal stress and crack propagation problems of fractured rock masses.展开更多
The aperture of natural rock fractures significantly affects the deformation and strength properties of rock masses,as well as the hydrodynamic properties of fractured rock masses.The conventional measurement methods ...The aperture of natural rock fractures significantly affects the deformation and strength properties of rock masses,as well as the hydrodynamic properties of fractured rock masses.The conventional measurement methods are inadequate for collecting data on high-steep rock slopes in complex mountainous regions.This study establishes a high-resolution three-dimensional model of a rock slope using unmanned aerial vehicle(UAV)multi-angle nap-of-the-object photogrammetry to obtain edge feature points of fractures.Fracture opening morphology is characterized using coordinate projection and transformation.Fracture central axis is determined using vertical measuring lines,allowing for the interpretation of aperture of adaptive fracture shape.The feasibility and reliability of the new method are verified at a construction site of a railway in southeast Tibet,China.The study shows that the fracture aperture has a significant interval effect and size effect.The optimal sampling length for fractures is approximately 0.5e1 m,and the optimal aperture interpretation results can be achieved when the measuring line spacing is 1%of the sampling length.Tensile fractures in the study area generally have larger apertures than shear fractures,and their tendency to increase with slope height is also greater than that of shear fractures.The aperture of tensile fractures is generally positively correlated with their trace length,while the correlation between the aperture of shear fractures and their trace length appears to be weak.Fractures of different orientations exhibit certain differences in their distribution of aperture,but generally follow the forms of normal,log-normal,and gamma distributions.This study provides essential data support for rock and slope stability evaluation,which is of significant practical importance.展开更多
Natural slopes usually display complicated exposed rock surfaces that are characterized by complex and substantial terrain undulation and ubiquitous undesirable phenomena such as vegetation cover and rockfalls.This st...Natural slopes usually display complicated exposed rock surfaces that are characterized by complex and substantial terrain undulation and ubiquitous undesirable phenomena such as vegetation cover and rockfalls.This study presents a systematic outcrop research of fracture pattern variations in a complicated rock slope,and the qualitative and quantitative study of the complex phenomena impact on threedimensional(3D)discrete fracture network(DFN)modeling.As the studies of the outcrop fracture pattern have been so far focused on local variations,thus,we put forward a statistical analysis of global variations.The entire outcrop is partitioned into several subzones,and the subzone-scale variability of fracture geometric properties is analyzed(including the orientation,the density,and the trace length).The results reveal significant variations in fracture characteristics(such as the concentrative degree,the average orientation,the density,and the trace length)among different subzones.Moreover,the density of fracture sets,which is approximately parallel to the slope surface,exhibits a notably higher value compared to other fracture sets across all subzones.To improve the accuracy of the DFN modeling,the effects of three common phenomena resulting from vegetation and rockfalls are qualitatively analyzed and the corresponding quantitative data processing solutions are proposed.Subsequently,the 3D fracture geometric parameters are determined for different areas of the high-steep rock slope in terms of the subzone dimensions.The results show significant variations in the same set of 3D fracture parameters across different regions with density differing by up to tenfold and mean trace length exhibiting differences of 3e4 times.The study results present precise geological structural information,improve modeling accuracy,and provide practical solutions for addressing complex outcrop issues.展开更多
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...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.展开更多
Stability analysis of underground constructions requires a model study of rock masses’ long-term performance. Creep tests under different stress conditions was conducted on intact granite and granite samples fracture...Stability analysis of underground constructions requires a model study of rock masses’ long-term performance. Creep tests under different stress conditions was conducted on intact granite and granite samples fractured at 30° and 45° angles. The experimental results indicate that the steady creep strain rates of intact and fractured rock present an exponential increase trend with the increase of stress level. A nonlinear creep model is developed based on the experimental results, in which the initial damage caused by fracture together with the damage caused by constant load have been taken into consideration. The fitting analysis results indicated that the model proposed is more accurate at identifying the full creep regions in fractured granite, especially the accelerated stage of creep deformation. The least-square fit error of the proposed creep model is significantly lower than that of Nishihara model by almost an order of magnitude. An analysis of the effects of elastic modulus, viscosity coefficient, and damage factors on fractured rock strain rate and creep strain is conducted. If no consideration is given to the effects of the damage, the proposed nonlinear creep model can degenerate into to the classical Nishihara model.展开更多
The surrounding rock of roadways exhibits intricate characteristics of discontinuity and heterogeneity.To address these complexities,this study employs non-local Peridynamics(PD)theory and reconstructs the kernel func...The surrounding rock of roadways exhibits intricate characteristics of discontinuity and heterogeneity.To address these complexities,this study employs non-local Peridynamics(PD)theory and reconstructs the kernel function to represent accurately the spatial decline of long-range force.Additionally,modifications to the traditional bondbased PD model are made.By considering the micro-structure of coal-rock materials within a uniform discrete model,heterogeneity characterized by bond random pre-breaking is introduced.This approach facilitates the proposal of a novel model capable of handling the random distribution characteristics of material heterogeneity,rendering the PD model suitable for analyzing the deformation and failure of heterogeneous layered coal-rock mass structures.The established numerical model and simulation method,termed the sub-homogeneous PD model,not only incorporates the support effect but also captures accurately the random heterogeneous micro-structure of roadway surrounding rock.The simulation results obtained using this model show good agreement with field measurements from the Fucun coal mine,effectively validating the model’s capability in accurately reproducing the deformation and failure mode of surrounding rock under bolt-supported(anchor cable).The proposed subhomogeneous PD model presents a valuable and effective simulation tool for studying the deformation and failure of roadway surrounding rock in coal mines,offering new insights and potential advancements.展开更多
Considering the importance of fractured rock aquifers in the hydrogeologic process,this research aimed to analyze the flow regime,internal degree of karstification,and estimate storage volume in fractured rock aquifer...Considering the importance of fractured rock aquifers in the hydrogeologic process,this research aimed to analyze the flow regime,internal degree of karstification,and estimate storage volume in fractured rock aquifers of the Germi Chai Basin in northwest Iran,which is attributed to its active tectonics,erosion,and the lithological diversity.Given the geological setting,the hypothesis is that this basin is characterized by a high degree of karstification and diffuse or intermediate flow regime leading to variation in discharge flow rate.The hydrodynamic and hadrochemical analysis was conducted on 9 well distributed springs across the basin from 2019 to 2020.The maximum flow rate in most of the springs appeared in the early wet season despite their different levels of fluctuations on the monthly discharge time series.Analyzing the spring recession curve form revealed an aquifer containing multiple micro-regimes withαrecession coefficients and a degree of karstification ranging between 0.001 to 0.06 and 0.55 to 2.61,respectively.These findings indicated a dominant diffuse and intermediate flow system resulting from the development of a high density of fractures in this area.The electrical conductivity of the spring changes inversely proportional to the change in flow discharge,indicating the reasonable hydrological response of the aquifer to rainfall events.Hydrograph analysis revealed that the delay time of spring discharge after rainfall events mostly varies between 10 to 30 days.The total dynamic storage volume of the spring for a given period(2019-2020)was estimated to be approximately 1324 million cubic meters reflecting the long-term drainage potential and high perdurability of dynamic storage.Estimating the maximum and minimum ratio revealed that the springs recharging system in Germi Chai Basin comes under the slow aquifers category.This finding provides valuable insight into the hydrogeological properties of fractured rock aquifers contributing to effective water management strategy.展开更多
The fractured surrounding rocks of roadways pose major challenges to safe mining.Grouting has often been used to reinforce the surrounding rocks to mitigate the safety risks associated with fractured rocks.The aim of ...The fractured surrounding rocks of roadways pose major challenges to safe mining.Grouting has often been used to reinforce the surrounding rocks to mitigate the safety risks associated with fractured rocks.The aim of this study is to develop highly efficient composite ultrafine cement(CUC)grouts to reinforce the roadway in fractured surrounding rocks.The materials used are ultrafine cement(UC),ultrafine fly ash(UF),ultrafine slag(US),and additives(superplasticizer[SUP],aluminate ultrafine expansion agent[AUA],gypsum,and retarder).The fluidity,bleeding,shrinkage,setting time,chemical composition,microstructure,degree of hydration,and mechanical property of grouting materials were evaluated in this study.Also,a suitable and effective CUC grout mixture was used to reinforce the roadway in the fractured surrounding rock.The results have shown that the addition of UF and US reduces the plastic viscosity of CUC,and the best fluidity can be obtained by adding 40%UF and 10%US.Since UC and UF particles are small,the pozzolanic effect of UF promotes the hydration reaction,which is conductive to the stability of CUC grouts.In addition,fine particles of UC,UF,and US can effectively fill the pores,while the volumetric expansion of AUA and gypsum decreases the pores and thus affects the microstructure of the solidified grout.The compressive test results have shown that the addition of specific amounts of UF and US can ameliorate the mechanical properties of CUC grouts.Finally,the CUC22‐8 grout was used to reinforce the No.20322 belt roadway.The results of numerical simulation and field monitoring have indicated that grouting can efficaciously reinforce the surrounding rock of the roadway.In this research,high‐performance CUC grouts were developed for surrounding rock reinforcement of underground engineering by utilizing UC and some additives.展开更多
To gain insight into the flow mechanisms and stress sensitivity for fractured-vuggy reservoirs,several core models with different structural characteristics were designed and fabricated to investigate the impact of ef...To gain insight into the flow mechanisms and stress sensitivity for fractured-vuggy reservoirs,several core models with different structural characteristics were designed and fabricated to investigate the impact of effective stress on permeability for carbonate fractured-vuggy rocks(CFVR).It shows that the permeability performance curves under different pore and confining pressures(i.e.altered stress conditions)for the fractured core models and the vuggy core models have similar change patterns.The ranges of permeability variation are significantly wider at high pore pressures,indicating that permeability reduction is the most significant during the early stage of development for fractured-vuggy reservoirs.Since each obtained effective stress coefficient for permeability(ESCP)varies with the changes in confining pressure and pore pressure,the effective stresses for permeability of four representative CFVR show obvious nonlinear characteristics,and the variation ranges of ESCP are all between 0 and 1.Meanwhile,a comprehensive ESCP mathematical model considering triple media,including matrix pores,fractures,and dissolved vugs,was proposed.It is proved theoretically that the ESCP of CFVR generally varies between 0 and 1.Additionally,the regression results showed that the power model ranked highest among the four empirical models mainly applied in stress sensitivity characterization,followed by the logarithmic model,exponential model,and binomial model.The concept of“permeability decline rate”was introduced to better evaluate the stress sensitivity performance for CFVR,in which the one-fracture rock is the strongest,followed by the fracture-vug rock and two-horizontalfracture rock;the through-hole rock is the weakest.In general,this study provides a theoretical basis to guide the design of development and adjustment programs for carbonate fractured-vuggy reservoirs.展开更多
Accurately picking P-and S-wave arrivals of microseismic(MS)signals in real-time directly influences the early warning of rock mass failure.A common contradiction between accuracy and computation exists in the current...Accurately picking P-and S-wave arrivals of microseismic(MS)signals in real-time directly influences the early warning of rock mass failure.A common contradiction between accuracy and computation exists in the current arrival picking methods.Thus,a real-time arrival picking method of MS signals is constructed based on a convolutional-recurrent neural network(CRNN).This method fully utilizes the advantages of convolutional layers and gated recurrent units(GRU)in extracting short-and long-term features,in order to create a precise and lightweight arrival picking structure.Then,the synthetic signals with field noises are used to evaluate the hyperparameters of the CRNN model and obtain an optimal CRNN model.The actual operation on various devices indicates that compared with the U-Net method,the CRNN method achieves faster arrival picking with less performance consumption.An application of large underground caverns in the Yebatan hydropower station(YBT)project shows that compared with the short-term average/long-term average(STA/LTA),Akaike information criterion(AIC)and U-Net methods,the CRNN method has the highest accuracy within four sampling points,which is 87.44%for P-wave and 91.29%for S-wave,respectively.The sum of mean absolute errors(MAESUM)of the CRNN method is 4.22 sampling points,which is lower than that of the other methods.Among the four methods,the MS sources location calculated based on the CRNN method shows the best consistency with the actual failure,which occurs at the junction of the shaft and the second gallery.Thus,the proposed method can pick up P-and S-arrival accurately and rapidly,providing a reference for rock failure analysis and evaluation in engineering applications.展开更多
Rock slope with horizontal-layered fractured structure(HLFS)has high stability in its natural state.However,a strong earthquake can induce rock fissure expansion,ultimately leading to slope failure.In this study,the d...Rock slope with horizontal-layered fractured structure(HLFS)has high stability in its natural state.However,a strong earthquake can induce rock fissure expansion,ultimately leading to slope failure.In this study,the dynamic response,failure mode,and spectral characteristics of rock slope with HLFS under strong earthquake conditions were investigated based on the large-scale shaking table model test.On this basis,multiple sets of numerical calculation models were further established by UDEC discrete element program.Five influencing factors were considered in the parametric study of numerical simulations,including slope height,slope angle,bedding-plane spacing and secondary joint spacing as well as bedrock dip angle.The results showed that the failure process of rock slope with HLFS under earthquake action is mainly divided into four phases,i.e.,the tensile crack of the slope shoulder joints and shear dislocation at the top bedding plane,the extension of vertical joint cracks and increase of shear displacement,the formation of step-through sliding surfaces and the instability,and finally collapse of fractured rock mass.The acceleration response of slopes exhibits elevation amplification effect and surface effect.Numerical simulations indicate that the seismic stability of slopes with HLFS exhibits a negative correlation with slope height and angle,but a positive correlation with bedding-plane spacing,joint spacing,and bedrock dip angle.The results of this study can provide a reference for seismic stability evaluation of weathered rock slopes.展开更多
This study employs similar simulation testing and discrete element simulation coupling to analyze the failure and deformation processes of a model coal seam's roof.The caving area of the overburden rock is divided...This study employs similar simulation testing and discrete element simulation coupling to analyze the failure and deformation processes of a model coal seam's roof.The caving area of the overburden rock is divided into three zones:the delamination fracture zone,broken fracture zone,and compaction zone.The caving and fracture zones'heights are approximately 110 m above the coal seam,with a maximum subsidence of 11 m.The delamination fracture zone's porosity range is between 0.2 and 0.3,while the remainder of the roof predominantly exhibits a porosity of less than 0.1.In addition,the numerical model's stress analysis revealed that the overburden rock's displacement zone forms an'arch-beam'structure starting from 160 m,with the maximum and minimum stress values decreasing as the distance of advancement increases.In the stress beam interval of the overburden rock,the maximum value changes periodically as the advancement distance increases.Based on a comparative analysis between observable data from on-site work and numerical simulation results,the stress data from the numerical simulation are essentially consistent with the actual results detected on-site,indicating the validity of the numerical simulation results.展开更多
Due to the presence of ice and unfrozen water in pores of frozen rock,the rock fracture behaviors are susceptible to temperature.In this study,the potential thawing-induced softening effects on the fracture behaviors ...Due to the presence of ice and unfrozen water in pores of frozen rock,the rock fracture behaviors are susceptible to temperature.In this study,the potential thawing-induced softening effects on the fracture behaviors of frozen rock is evaluated by testing the tension fracture toughness(KIC)of frozen rock at different temperatures(i.e.-20℃,-15℃,-12℃,-10℃,-8℃,-6℃,-4℃,-2℃,and 0℃).Acoustic emission(AE)and digital image correlation(DIC)methods are utilized to analyze the microcrack propagation during fracturing.The melting of pore ice is measured using nuclear magnetic resonance(NMR)method.The results indicate that:(1)The KIC of frozen rock decreases moderately between-20℃ and-4℃,and rapidly between-4℃ and 0℃.(2)At-20℃ to-4℃,the fracturing process,deduced from the DIC results at the notch tip,exhibits three stages:elastic deformation,microcrack propagation and microcrack coalescence.However,at-4℃e0℃,only the latter two stages are observed.(3)At-4℃e0℃,the AE activities during fracturing are less than that at-20℃ to-4℃,while more small events are reported.(4)The NMR results demonstrate a reverse variation trend in pore ice content with increasing temperature,that is,a moderate decrease is followed by a sharp decrease and-4℃ is exactly the critical temperature.Next,we interpret the thawing-induced softening effect by linking the evolution in microscopic structure of frozen rock with its macroscopic fracture behaviors as follow:from-20℃ to-4℃,the thickening of the unfrozen water film diminishes the cementation strength between ice and rock skeleton,leading to the decrease in fracture parameters.From-4℃ to 0℃,the cementation effect of ice almost vanishes,and the filling effect of pore ice is reduced significantly,which facilitates microcrack propagation and thus the easier fracture of frozen rocks.展开更多
P-and SV-wave dispersion and attenuation have been extensively investigated in saturated poroelastic media with aligned fractures.However,there are few existing models that incorporate the multiple wave attenuation me...P-and SV-wave dispersion and attenuation have been extensively investigated in saturated poroelastic media with aligned fractures.However,there are few existing models that incorporate the multiple wave attenuation mechanisms from the microscopic scale to the macroscopic scale.Hence,in this work,we developed a unified model to incorporate the wave attenuation mechanisms at different scales,which includes the microscopic squirt flow between the microcracks and pores,the mesoscopic wave-induced fluid flow between fractures and background(FB-WIFF),and the macroscopic Biot's global flow and elastic scattering(ES)from the fractures.Using Tang's modified Biot's theory and the mixed-boundary conditions,we derived the exact frequency-dependent solutions of the scattering problem for a single penny-shaped fracture with oblique incident P-and SV-waves.We then developed theoretical models for a set of aligned fractures and randomly oriented fractures using the Foldy approximation.The results indicated that microcrack squirt flow considerably influences the dispersion and attenuation of P-and SV-wave velocities.The coupling effects of microcrack squirt flow with the FB-WIFF and ES of fractures cause much higher velocity dispersion and attenuation for P waves than for SV waves.Randomly oriented fractures substantially reduce the attenuation caused by the FB-WIFF and ES,particularly for the ES attenuation of SV waves.Through a comparison with existing models in the limiting cases and previous experimental measurements,we validated our model.展开更多
Grouting is a widely used approach to reinforce broken surrounding rock mass during the construction of underground tunnels in fault fracture zones,and its reinforcement effectiveness is highly affected by geostress.I...Grouting is a widely used approach to reinforce broken surrounding rock mass during the construction of underground tunnels in fault fracture zones,and its reinforcement effectiveness is highly affected by geostress.In this study,a numerical manifold method(NMM)based simulator has been developed to examine the impact of geostress conditions on grouting reinforcement during tunnel excavation.To develop this simulator,a detection technique for identifying slurry migration channels and an improved fluid-solid coupling(FeS)framework,which considers the influence of fracture properties and geostress states,is developed and incorporated into a zero-thickness cohesive element(ZE)based NMM(Co-NMM)for simulating tunnel excavation.Additionally,to simulate coagulation of injected slurry,a bonding repair algorithm is further proposed based on the ZE model.To verify the accuracy of the proposed simulator,a series of simulations about slurry migration in single fractures and fracture networks are numerically reproduced,and the results align well with analytical and laboratory test results.Furthermore,these numerical results show that neglecting the influence of geostress condition can lead to a serious over-estimation of slurry migration range and reinforcement effectiveness.After validations,a series of simulations about tunnel grouting reinforcement and tunnel excavation in fault fracture zones with varying fracture densities under different geostress conditions are conducted.Based on these simula-tions,the influence of geostress conditions and the optimization of grouting schemes are discussed.展开更多
For expedited transportation,vehicular tunnels are often designed as two adjacent tunnels,which frequently experience dynamic stress waves from various orientations during blasting excavation.To analyze the impact of ...For expedited transportation,vehicular tunnels are often designed as two adjacent tunnels,which frequently experience dynamic stress waves from various orientations during blasting excavation.To analyze the impact of dynamic loading orientation on the stability of the twin-tunnel,a split Hopkinson pressure bar(SHPB)apparatus was used to conduct a dynamic test on the twin-tunnel specimens.The two tunnels were rotated around the specimen’s center to consider the effect of dynamic loading orientation.LS-DYNA software was used for numerical simulation to reveal the failure properties and stress wave propagation law of the twin-tunnel specimens.The findings indicate that,for a twin-tunnel exposed to a dynamic load from different orientations,the crack initiation position appears most often at the tunnel corner,tunnel spandrel,and tunnel floor.As the impact direction is created by a certain angle(30°,45°,60°,120°,135°,and 150°),the fractures are produced in the middle of the line between the left tunnel corner and the right tunnel spandrel.As the impact loading angle(a)is 90°,the tunnel sustains minimal damage,and only tensile fractures form in the surrounding rocks.The orientation of the impact load could change the stress distribution in the twin-tunnel,and major fractures are more likely to form in areas where the tensile stress is concentrated.展开更多
The damage-fracture evolution of deep rock mass has obvious particularity,which is revealed in 2400-mdeep tunnels by field tests.The evolution of the excavation damaged zone depth is consistent with that of the fractu...The damage-fracture evolution of deep rock mass has obvious particularity,which is revealed in 2400-mdeep tunnels by field tests.The evolution of the excavation damaged zone depth is consistent with that of the fractured zone depth.The ratio of the excavation damaged zone depth to the excavation fractured zone depth is greater than 2.0 in a rock mass with both high strength and good integrity,but less than1.5 in a rock mass with lower strength or poor integrity.Zonal disintegration in a rock mass with high strength and fair integrity is more likely to occur when it contains more than two groups of primary fractures in damaged zones.Fractures develop outward in zonal disintegration but are totally different from the single-zone fracture,in which the fractures develop inward,and it is the starting position of the fractured zone when the excavation surface of the middle pilot is 7–9 m close to the pre-set borehole and it stops after the excavation surface of the baseplate is 11–14 m away.The most intense evolution occurs around 2–4 m from the pre-set borehole in the sidewall expansion stage.The research results provide a reference for the monitoring scheme and support design of CJPL-Ⅲin its future construction.展开更多
基金supported by the National Natural Science Foundation of China(Grant No.42277165)the Fundamental Research Funds for the Central Universities,China University of Geosciences(Wuhan)(Grant No.CUGCJ1821)the National Overseas Study Fund(Grant No.202106410040).
文摘As a calculation method based on the Galerkin variation,the numerical manifold method(NMM)adopts a double covering system,which can easily deal with discontinuous deformation problems and has a high calculation accuracy.Aiming at the thermo-mechanical(TM)coupling problem of fractured rock masses,this study uses the NMM to simulate the processes of crack initiation and propagation in a rock mass under the influence of temperature field,deduces related system equations,and proposes a penalty function method to deal with boundary conditions.Numerical examples are employed to confirm the effectiveness and high accuracy of this method.By the thermal stress analysis of a thick-walled cylinder(TWC),the simulation of cracking in the TWC under heating and cooling conditions,and the simulation of thermal cracking of the SwedishÄspöPillar Stability Experiment(APSE)rock column,the thermal stress,and TM coupling are obtained.The numerical simulation results are in good agreement with the test data and other numerical results,thus verifying the effectiveness of the NMM in dealing with thermal stress and crack propagation problems of fractured rock masses.
基金This work was supported by the National Nature Science Foundation of China(Grant Nos.42177139 and 41941017)the Natural Science Foundation Project of Jilin Province,China(Grant No.20230101088JC).The authors would like to thank the anonymous reviewers for their comments and suggestions.
文摘The aperture of natural rock fractures significantly affects the deformation and strength properties of rock masses,as well as the hydrodynamic properties of fractured rock masses.The conventional measurement methods are inadequate for collecting data on high-steep rock slopes in complex mountainous regions.This study establishes a high-resolution three-dimensional model of a rock slope using unmanned aerial vehicle(UAV)multi-angle nap-of-the-object photogrammetry to obtain edge feature points of fractures.Fracture opening morphology is characterized using coordinate projection and transformation.Fracture central axis is determined using vertical measuring lines,allowing for the interpretation of aperture of adaptive fracture shape.The feasibility and reliability of the new method are verified at a construction site of a railway in southeast Tibet,China.The study shows that the fracture aperture has a significant interval effect and size effect.The optimal sampling length for fractures is approximately 0.5e1 m,and the optimal aperture interpretation results can be achieved when the measuring line spacing is 1%of the sampling length.Tensile fractures in the study area generally have larger apertures than shear fractures,and their tendency to increase with slope height is also greater than that of shear fractures.The aperture of tensile fractures is generally positively correlated with their trace length,while the correlation between the aperture of shear fractures and their trace length appears to be weak.Fractures of different orientations exhibit certain differences in their distribution of aperture,but generally follow the forms of normal,log-normal,and gamma distributions.This study provides essential data support for rock and slope stability evaluation,which is of significant practical importance.
基金supported by the National Key Research and Development Program of China(Grant No.2022YFC3080200)the National Natural Science Foundation of China(Grant No.42022053)the China Postdoctoral Science Foundation(Grant No.2023M731264).
文摘Natural slopes usually display complicated exposed rock surfaces that are characterized by complex and substantial terrain undulation and ubiquitous undesirable phenomena such as vegetation cover and rockfalls.This study presents a systematic outcrop research of fracture pattern variations in a complicated rock slope,and the qualitative and quantitative study of the complex phenomena impact on threedimensional(3D)discrete fracture network(DFN)modeling.As the studies of the outcrop fracture pattern have been so far focused on local variations,thus,we put forward a statistical analysis of global variations.The entire outcrop is partitioned into several subzones,and the subzone-scale variability of fracture geometric properties is analyzed(including the orientation,the density,and the trace length).The results reveal significant variations in fracture characteristics(such as the concentrative degree,the average orientation,the density,and the trace length)among different subzones.Moreover,the density of fracture sets,which is approximately parallel to the slope surface,exhibits a notably higher value compared to other fracture sets across all subzones.To improve the accuracy of the DFN modeling,the effects of three common phenomena resulting from vegetation and rockfalls are qualitatively analyzed and the corresponding quantitative data processing solutions are proposed.Subsequently,the 3D fracture geometric parameters are determined for different areas of the high-steep rock slope in terms of the subzone dimensions.The results show significant variations in the same set of 3D fracture parameters across different regions with density differing by up to tenfold and mean trace length exhibiting differences of 3e4 times.The study results present precise geological structural information,improve modeling accuracy,and provide practical solutions for addressing complex outcrop issues.
基金financially supported by National Natural Science Foundation of China(Grant Nos.52074312 and 52211530097)CNPC Science and Technology Innovation Foundation(Grant No.2021DQ02-0505).
文摘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.
基金supported by the National Natural Science Foundation of China(No.42307258)the technological research projects in Sichuan Province(No.2022YFSY0007)the China Atomic Energy Authority(CAEA)through the Geological Disposal Program.
文摘Stability analysis of underground constructions requires a model study of rock masses’ long-term performance. Creep tests under different stress conditions was conducted on intact granite and granite samples fractured at 30° and 45° angles. The experimental results indicate that the steady creep strain rates of intact and fractured rock present an exponential increase trend with the increase of stress level. A nonlinear creep model is developed based on the experimental results, in which the initial damage caused by fracture together with the damage caused by constant load have been taken into consideration. The fitting analysis results indicated that the model proposed is more accurate at identifying the full creep regions in fractured granite, especially the accelerated stage of creep deformation. The least-square fit error of the proposed creep model is significantly lower than that of Nishihara model by almost an order of magnitude. An analysis of the effects of elastic modulus, viscosity coefficient, and damage factors on fractured rock strain rate and creep strain is conducted. If no consideration is given to the effects of the damage, the proposed nonlinear creep model can degenerate into to the classical Nishihara model.
基金supported by the National Natural Science Foundation of China(Nos.12302264,52104004,12072170,and 12202225)the Natural Science Foundation of Shandong Province(No.ZR2021QA042)Special Fund for Taishan Scholar Project(No.Tsqn202211180).
文摘The surrounding rock of roadways exhibits intricate characteristics of discontinuity and heterogeneity.To address these complexities,this study employs non-local Peridynamics(PD)theory and reconstructs the kernel function to represent accurately the spatial decline of long-range force.Additionally,modifications to the traditional bondbased PD model are made.By considering the micro-structure of coal-rock materials within a uniform discrete model,heterogeneity characterized by bond random pre-breaking is introduced.This approach facilitates the proposal of a novel model capable of handling the random distribution characteristics of material heterogeneity,rendering the PD model suitable for analyzing the deformation and failure of heterogeneous layered coal-rock mass structures.The established numerical model and simulation method,termed the sub-homogeneous PD model,not only incorporates the support effect but also captures accurately the random heterogeneous micro-structure of roadway surrounding rock.The simulation results obtained using this model show good agreement with field measurements from the Fucun coal mine,effectively validating the model’s capability in accurately reproducing the deformation and failure mode of surrounding rock under bolt-supported(anchor cable).The proposed subhomogeneous PD model presents a valuable and effective simulation tool for studying the deformation and failure of roadway surrounding rock in coal mines,offering new insights and potential advancements.
基金the Water and Wastewater Company of East Azarbaijan Province for providing the funding for this research
文摘Considering the importance of fractured rock aquifers in the hydrogeologic process,this research aimed to analyze the flow regime,internal degree of karstification,and estimate storage volume in fractured rock aquifers of the Germi Chai Basin in northwest Iran,which is attributed to its active tectonics,erosion,and the lithological diversity.Given the geological setting,the hypothesis is that this basin is characterized by a high degree of karstification and diffuse or intermediate flow regime leading to variation in discharge flow rate.The hydrodynamic and hadrochemical analysis was conducted on 9 well distributed springs across the basin from 2019 to 2020.The maximum flow rate in most of the springs appeared in the early wet season despite their different levels of fluctuations on the monthly discharge time series.Analyzing the spring recession curve form revealed an aquifer containing multiple micro-regimes withαrecession coefficients and a degree of karstification ranging between 0.001 to 0.06 and 0.55 to 2.61,respectively.These findings indicated a dominant diffuse and intermediate flow system resulting from the development of a high density of fractures in this area.The electrical conductivity of the spring changes inversely proportional to the change in flow discharge,indicating the reasonable hydrological response of the aquifer to rainfall events.Hydrograph analysis revealed that the delay time of spring discharge after rainfall events mostly varies between 10 to 30 days.The total dynamic storage volume of the spring for a given period(2019-2020)was estimated to be approximately 1324 million cubic meters reflecting the long-term drainage potential and high perdurability of dynamic storage.Estimating the maximum and minimum ratio revealed that the springs recharging system in Germi Chai Basin comes under the slow aquifers category.This finding provides valuable insight into the hydrogeological properties of fractured rock aquifers contributing to effective water management strategy.
基金supported by the National Natural Science Foundation of China(NSFC)(grant No.52074169,No.51704280)the China Postdoctoral Science Foundation(No.2023M732109)the Opening Foundation of Shandong Key Laboratory of Civil Engineering Disaster Prevention and Mitigation(No.CDPM2021FK02).
文摘The fractured surrounding rocks of roadways pose major challenges to safe mining.Grouting has often been used to reinforce the surrounding rocks to mitigate the safety risks associated with fractured rocks.The aim of this study is to develop highly efficient composite ultrafine cement(CUC)grouts to reinforce the roadway in fractured surrounding rocks.The materials used are ultrafine cement(UC),ultrafine fly ash(UF),ultrafine slag(US),and additives(superplasticizer[SUP],aluminate ultrafine expansion agent[AUA],gypsum,and retarder).The fluidity,bleeding,shrinkage,setting time,chemical composition,microstructure,degree of hydration,and mechanical property of grouting materials were evaluated in this study.Also,a suitable and effective CUC grout mixture was used to reinforce the roadway in the fractured surrounding rock.The results have shown that the addition of UF and US reduces the plastic viscosity of CUC,and the best fluidity can be obtained by adding 40%UF and 10%US.Since UC and UF particles are small,the pozzolanic effect of UF promotes the hydration reaction,which is conductive to the stability of CUC grouts.In addition,fine particles of UC,UF,and US can effectively fill the pores,while the volumetric expansion of AUA and gypsum decreases the pores and thus affects the microstructure of the solidified grout.The compressive test results have shown that the addition of specific amounts of UF and US can ameliorate the mechanical properties of CUC grouts.Finally,the CUC22‐8 grout was used to reinforce the No.20322 belt roadway.The results of numerical simulation and field monitoring have indicated that grouting can efficaciously reinforce the surrounding rock of the roadway.In this research,high‐performance CUC grouts were developed for surrounding rock reinforcement of underground engineering by utilizing UC and some additives.
基金This work was supported by the Joint Fund of NSFC for Enterprise Innovation and Development(Grant No.U19B6003-02-06)the National Natural Science Foundation of China(Grant No.51974331)+1 种基金the Natural Science Foundation of Jiangsu Province(Grant No.BK20200525)The authors would like to sincerely acknowledge these funding programs for their financial support.Particularly,the support provided by the China Scholarship Council(CSC)during a visit of Ke Sun(File No.202106440065)to the University of Alberta is also sincerely acknowledged.
文摘To gain insight into the flow mechanisms and stress sensitivity for fractured-vuggy reservoirs,several core models with different structural characteristics were designed and fabricated to investigate the impact of effective stress on permeability for carbonate fractured-vuggy rocks(CFVR).It shows that the permeability performance curves under different pore and confining pressures(i.e.altered stress conditions)for the fractured core models and the vuggy core models have similar change patterns.The ranges of permeability variation are significantly wider at high pore pressures,indicating that permeability reduction is the most significant during the early stage of development for fractured-vuggy reservoirs.Since each obtained effective stress coefficient for permeability(ESCP)varies with the changes in confining pressure and pore pressure,the effective stresses for permeability of four representative CFVR show obvious nonlinear characteristics,and the variation ranges of ESCP are all between 0 and 1.Meanwhile,a comprehensive ESCP mathematical model considering triple media,including matrix pores,fractures,and dissolved vugs,was proposed.It is proved theoretically that the ESCP of CFVR generally varies between 0 and 1.Additionally,the regression results showed that the power model ranked highest among the four empirical models mainly applied in stress sensitivity characterization,followed by the logarithmic model,exponential model,and binomial model.The concept of“permeability decline rate”was introduced to better evaluate the stress sensitivity performance for CFVR,in which the one-fracture rock is the strongest,followed by the fracture-vug rock and two-horizontalfracture rock;the through-hole rock is the weakest.In general,this study provides a theoretical basis to guide the design of development and adjustment programs for carbonate fractured-vuggy reservoirs.
基金We acknowledge the funding support from National Natural Science Foundation of China(Grant No.42077263).
文摘Accurately picking P-and S-wave arrivals of microseismic(MS)signals in real-time directly influences the early warning of rock mass failure.A common contradiction between accuracy and computation exists in the current arrival picking methods.Thus,a real-time arrival picking method of MS signals is constructed based on a convolutional-recurrent neural network(CRNN).This method fully utilizes the advantages of convolutional layers and gated recurrent units(GRU)in extracting short-and long-term features,in order to create a precise and lightweight arrival picking structure.Then,the synthetic signals with field noises are used to evaluate the hyperparameters of the CRNN model and obtain an optimal CRNN model.The actual operation on various devices indicates that compared with the U-Net method,the CRNN method achieves faster arrival picking with less performance consumption.An application of large underground caverns in the Yebatan hydropower station(YBT)project shows that compared with the short-term average/long-term average(STA/LTA),Akaike information criterion(AIC)and U-Net methods,the CRNN method has the highest accuracy within four sampling points,which is 87.44%for P-wave and 91.29%for S-wave,respectively.The sum of mean absolute errors(MAESUM)of the CRNN method is 4.22 sampling points,which is lower than that of the other methods.Among the four methods,the MS sources location calculated based on the CRNN method shows the best consistency with the actual failure,which occurs at the junction of the shaft and the second gallery.Thus,the proposed method can pick up P-and S-arrival accurately and rapidly,providing a reference for rock failure analysis and evaluation in engineering applications.
基金supported by Central Guiding Local Science and Technology Development Special Fund Project(No.ZYYD2023B02)the National Natural Science Foundation of China(Nos.52078432 and 52168066)the Scientific Research Project of China Railway First Survey and Design Institute Group Co.(No.20-06).
文摘Rock slope with horizontal-layered fractured structure(HLFS)has high stability in its natural state.However,a strong earthquake can induce rock fissure expansion,ultimately leading to slope failure.In this study,the dynamic response,failure mode,and spectral characteristics of rock slope with HLFS under strong earthquake conditions were investigated based on the large-scale shaking table model test.On this basis,multiple sets of numerical calculation models were further established by UDEC discrete element program.Five influencing factors were considered in the parametric study of numerical simulations,including slope height,slope angle,bedding-plane spacing and secondary joint spacing as well as bedrock dip angle.The results showed that the failure process of rock slope with HLFS under earthquake action is mainly divided into four phases,i.e.,the tensile crack of the slope shoulder joints and shear dislocation at the top bedding plane,the extension of vertical joint cracks and increase of shear displacement,the formation of step-through sliding surfaces and the instability,and finally collapse of fractured rock mass.The acceleration response of slopes exhibits elevation amplification effect and surface effect.Numerical simulations indicate that the seismic stability of slopes with HLFS exhibits a negative correlation with slope height and angle,but a positive correlation with bedding-plane spacing,joint spacing,and bedrock dip angle.The results of this study can provide a reference for seismic stability evaluation of weathered rock slopes.
基金National Key R&D Program of China(2023YFC3009100,2023YFC3009102)National Natural Science Foundation of China(52304198)Open Fund of the National and Local Joint Engineering Research Center for Safe and Accurate Coal Mining(EC2021016).
文摘This study employs similar simulation testing and discrete element simulation coupling to analyze the failure and deformation processes of a model coal seam's roof.The caving area of the overburden rock is divided into three zones:the delamination fracture zone,broken fracture zone,and compaction zone.The caving and fracture zones'heights are approximately 110 m above the coal seam,with a maximum subsidence of 11 m.The delamination fracture zone's porosity range is between 0.2 and 0.3,while the remainder of the roof predominantly exhibits a porosity of less than 0.1.In addition,the numerical model's stress analysis revealed that the overburden rock's displacement zone forms an'arch-beam'structure starting from 160 m,with the maximum and minimum stress values decreasing as the distance of advancement increases.In the stress beam interval of the overburden rock,the maximum value changes periodically as the advancement distance increases.Based on a comparative analysis between observable data from on-site work and numerical simulation results,the stress data from the numerical simulation are essentially consistent with the actual results detected on-site,indicating the validity of the numerical simulation results.
基金We acknowledge the funding support from the National Natural Science Foundation of China(Grant No.42271148).
文摘Due to the presence of ice and unfrozen water in pores of frozen rock,the rock fracture behaviors are susceptible to temperature.In this study,the potential thawing-induced softening effects on the fracture behaviors of frozen rock is evaluated by testing the tension fracture toughness(KIC)of frozen rock at different temperatures(i.e.-20℃,-15℃,-12℃,-10℃,-8℃,-6℃,-4℃,-2℃,and 0℃).Acoustic emission(AE)and digital image correlation(DIC)methods are utilized to analyze the microcrack propagation during fracturing.The melting of pore ice is measured using nuclear magnetic resonance(NMR)method.The results indicate that:(1)The KIC of frozen rock decreases moderately between-20℃ and-4℃,and rapidly between-4℃ and 0℃.(2)At-20℃ to-4℃,the fracturing process,deduced from the DIC results at the notch tip,exhibits three stages:elastic deformation,microcrack propagation and microcrack coalescence.However,at-4℃e0℃,only the latter two stages are observed.(3)At-4℃e0℃,the AE activities during fracturing are less than that at-20℃ to-4℃,while more small events are reported.(4)The NMR results demonstrate a reverse variation trend in pore ice content with increasing temperature,that is,a moderate decrease is followed by a sharp decrease and-4℃ is exactly the critical temperature.Next,we interpret the thawing-induced softening effect by linking the evolution in microscopic structure of frozen rock with its macroscopic fracture behaviors as follow:from-20℃ to-4℃,the thickening of the unfrozen water film diminishes the cementation strength between ice and rock skeleton,leading to the decrease in fracture parameters.From-4℃ to 0℃,the cementation effect of ice almost vanishes,and the filling effect of pore ice is reduced significantly,which facilitates microcrack propagation and thus the easier fracture of frozen rocks.
基金This work was supported by the Laoshan National Laboratory Science and Technology Innovation Project(No.LSKJ202203407)the National Natural Science Foundation of China(Grant Nos.42174145,41821002,42274146)+1 种基金Guangdong Provincial Key Laboratory of Geophysical High-resolution Imaging Technology(2022B1212010002)Shenzhen Stable Support Plan Program for Higher Education Institutions(20220815110144003).
文摘P-and SV-wave dispersion and attenuation have been extensively investigated in saturated poroelastic media with aligned fractures.However,there are few existing models that incorporate the multiple wave attenuation mechanisms from the microscopic scale to the macroscopic scale.Hence,in this work,we developed a unified model to incorporate the wave attenuation mechanisms at different scales,which includes the microscopic squirt flow between the microcracks and pores,the mesoscopic wave-induced fluid flow between fractures and background(FB-WIFF),and the macroscopic Biot's global flow and elastic scattering(ES)from the fractures.Using Tang's modified Biot's theory and the mixed-boundary conditions,we derived the exact frequency-dependent solutions of the scattering problem for a single penny-shaped fracture with oblique incident P-and SV-waves.We then developed theoretical models for a set of aligned fractures and randomly oriented fractures using the Foldy approximation.The results indicated that microcrack squirt flow considerably influences the dispersion and attenuation of P-and SV-wave velocities.The coupling effects of microcrack squirt flow with the FB-WIFF and ES of fractures cause much higher velocity dispersion and attenuation for P waves than for SV waves.Randomly oriented fractures substantially reduce the attenuation caused by the FB-WIFF and ES,particularly for the ES attenuation of SV waves.Through a comparison with existing models in the limiting cases and previous experimental measurements,we validated our model.
基金This work was supported by the Guangdong Basic and Applied Basic Research Foundation(Grant No.2021A1515110304)the Na-tional Natural Science Foundation of China(Grant Nos.42077246 and 52278412).
文摘Grouting is a widely used approach to reinforce broken surrounding rock mass during the construction of underground tunnels in fault fracture zones,and its reinforcement effectiveness is highly affected by geostress.In this study,a numerical manifold method(NMM)based simulator has been developed to examine the impact of geostress conditions on grouting reinforcement during tunnel excavation.To develop this simulator,a detection technique for identifying slurry migration channels and an improved fluid-solid coupling(FeS)framework,which considers the influence of fracture properties and geostress states,is developed and incorporated into a zero-thickness cohesive element(ZE)based NMM(Co-NMM)for simulating tunnel excavation.Additionally,to simulate coagulation of injected slurry,a bonding repair algorithm is further proposed based on the ZE model.To verify the accuracy of the proposed simulator,a series of simulations about slurry migration in single fractures and fracture networks are numerically reproduced,and the results align well with analytical and laboratory test results.Furthermore,these numerical results show that neglecting the influence of geostress condition can lead to a serious over-estimation of slurry migration range and reinforcement effectiveness.After validations,a series of simulations about tunnel grouting reinforcement and tunnel excavation in fault fracture zones with varying fracture densities under different geostress conditions are conducted.Based on these simula-tions,the influence of geostress conditions and the optimization of grouting schemes are discussed.
基金supported by the National Natural Science Foundation of China(Grant Nos.52204104 and U19A2098)the Science and Technology Department of Sichuan Province,China(Grant No.2023YFH0022).
文摘For expedited transportation,vehicular tunnels are often designed as two adjacent tunnels,which frequently experience dynamic stress waves from various orientations during blasting excavation.To analyze the impact of dynamic loading orientation on the stability of the twin-tunnel,a split Hopkinson pressure bar(SHPB)apparatus was used to conduct a dynamic test on the twin-tunnel specimens.The two tunnels were rotated around the specimen’s center to consider the effect of dynamic loading orientation.LS-DYNA software was used for numerical simulation to reveal the failure properties and stress wave propagation law of the twin-tunnel specimens.The findings indicate that,for a twin-tunnel exposed to a dynamic load from different orientations,the crack initiation position appears most often at the tunnel corner,tunnel spandrel,and tunnel floor.As the impact direction is created by a certain angle(30°,45°,60°,120°,135°,and 150°),the fractures are produced in the middle of the line between the left tunnel corner and the right tunnel spandrel.As the impact loading angle(a)is 90°,the tunnel sustains minimal damage,and only tensile fractures form in the surrounding rocks.The orientation of the impact load could change the stress distribution in the twin-tunnel,and major fractures are more likely to form in areas where the tensile stress is concentrated.
基金supported by the National Natural Science Foundation of China(Nos.51909136 and 42177168)Project of Youth Inno vation Promotion Association of Chinese Academy of Sciences(No.2021326)the Open Research Fund of Key Laboratory of Geological Hazards on Three Gorges Reservoir Area(China Three Gorges University),Ministry of Education(No.2020KDZ03)。
文摘The damage-fracture evolution of deep rock mass has obvious particularity,which is revealed in 2400-mdeep tunnels by field tests.The evolution of the excavation damaged zone depth is consistent with that of the fractured zone depth.The ratio of the excavation damaged zone depth to the excavation fractured zone depth is greater than 2.0 in a rock mass with both high strength and good integrity,but less than1.5 in a rock mass with lower strength or poor integrity.Zonal disintegration in a rock mass with high strength and fair integrity is more likely to occur when it contains more than two groups of primary fractures in damaged zones.Fractures develop outward in zonal disintegration but are totally different from the single-zone fracture,in which the fractures develop inward,and it is the starting position of the fractured zone when the excavation surface of the middle pilot is 7–9 m close to the pre-set borehole and it stops after the excavation surface of the baseplate is 11–14 m away.The most intense evolution occurs around 2–4 m from the pre-set borehole in the sidewall expansion stage.The research results provide a reference for the monitoring scheme and support design of CJPL-Ⅲin its future construction.