Understanding the scaling relation of damage zone width with displacement of faults is important for predicting subsurface faulting mechanisms and fluid flow processes. The understanding of this scaling relationship i...Understanding the scaling relation of damage zone width with displacement of faults is important for predicting subsurface faulting mechanisms and fluid flow processes. The understanding of this scaling relationship is influenced by the accuracy of the methods and types of data utilized to investigate faults. In this study, seismic reflection data are used to investigate the throw and damage zone width of five strike-slip faults a ecting Ordovician carbonates of the Tarim intracraton basin,NW China. The results indicate that fault slips with a throw less than 200 m had formed wide damage zones up to 3000 m in width. Also, damage zone width is found to have both a positive correlation and a power-law relation with throw of two orders of magnitude, with a ratio of these values varying in a range of 2–15. The relationship between throw and damage zone width is not a simple power-law and changes its slope from small to larger size faults. The results indicate that throw scales well with damage zone width for the studied faults, and hence these can be used to predict fault geometries in the Tarim Basin. The study of the wide carbonate damage zones presented here provides new insights into scaling of large-size faults, which involve multiple faulting stages.展开更多
During the construction of an underground excavation, damage occurs in the surrounding rock mass due in large part to stress changes. While the predicted damage extent impacts profile selection and support design, the...During the construction of an underground excavation, damage occurs in the surrounding rock mass due in large part to stress changes. While the predicted damage extent impacts profile selection and support design, the depth of damage is a critical aspect for the design of permeability sensitive excavations, such as a deep geological repository(DGR) for nuclear waste. Review of literature regarding the depth of excavation damage zones(EDZs) indicates three zones are common and typically related to stress induced damage. Based on past developments related to brittle damage prediction using continuum modelling, the depth of the EDZs has been examined numerically. One method to capture stress induced damage in conventional engineering software is the damage initiation and spalling limit(DISL) approach. The variability of depths predicted using the DISL approach has been evaluated and guidelines are suggested for determining the depth of the EDZs around circular excavations in brittle rock masses. Of the inputs evaluated, it was found that the tensile strength produces the greatest variation in the depth of the EDZs. The results were evaluated statistically to determine the best fit relation between the model inputs and the depth of the EDZs. The best correlation and least variation were found for the outer EDZ and the highly damaged zone(HDZ) showed the greatest variation. Predictive equations for different EDZs have been suggested and the maximum numerical EDZ depths, represented by the 68% prediction interval, agreed well with the empirical evidence. This suggests that the numerical limits can be used for preliminary depth prediction of the EDZs in brittle rock for circular excavations.展开更多
The objective of this paper is to develop a methodology for calibration of a discrete element grain-based model(GBM)to replicate the hydro-mechanical properties of a brittle rock measured in the laboratory,and to appl...The objective of this paper is to develop a methodology for calibration of a discrete element grain-based model(GBM)to replicate the hydro-mechanical properties of a brittle rock measured in the laboratory,and to apply the calibrated model to simulating the formation of excavation damage zone(EDZ)around underground excavations.Firstly,a new cohesive crack model is implemented into the universal distinct element code(UDEC)to control the fracturing behaviour of materials under various loading modes.Next,a methodology for calibration of the components of the UDEC-Voronoi model is discussed.The role of connectivity of induced microcracks on increasing the permeability of laboratory-scale samples is investigated.The calibrated samples are used to investigate the influence of pore fluid pressure on weakening the drained strength of the laboratory-scale rock.The validity of the Terzaghi’s effective stress law for the drained peak strength of low-porosity rock is tested by performing a series of biaxial compression test simulations.Finally,the evolution of damage and pore pressure around two unsupported circular tunnels in crystalline granitic rock is studied.展开更多
In the creep fatigue crack growth of GH4169 alloy,oxidation is a prominent damage source,which is mainly manifested as the oxidation damage zone in front of crack tip.In order to investigate the property of the oxidat...In the creep fatigue crack growth of GH4169 alloy,oxidation is a prominent damage source,which is mainly manifested as the oxidation damage zone in front of crack tip.In order to investigate the property of the oxidation damage zone formed in the creep fatigue crack growth,crack growth tests of directly aged GH4169 alloy were conducted at 650℃ in air under various load conditions.Interrupted tests were performed to observe the damage characteristics at crack tip.Block tests were systematically executed to quantify the dependency of oxidation damage zone size on load and holding time.The crack propagation of the GH4169 alloy has a close relationship with grain boundary oxidation at 650℃.An oxidation damage zone in front of crack tip includes intergranular microcracks and oxidised but uncracked grain boundaries.Its size has been calculated from transient crack growth rate and described as a function of maximum stress intensity factor and holding time.Based on oxidation damage zone size,a novel model has been developed to predict the creep fatigue crack growth rate of the GH4169 alloy at 650℃.展开更多
Purpose–The microseismic monitoring technique has great advantages on identifying the location,extent and the mechanism of damage process occurring in rock mass.This study aims to analyze distribution characteristics...Purpose–The microseismic monitoring technique has great advantages on identifying the location,extent and the mechanism of damage process occurring in rock mass.This study aims to analyze distribution characteristics and the evolution law of excavation damage zone of surrounding rock based on microseismic monitoring data.Design/methodology/approach–In situ test using microseismic monitoring technique is carried out in the large-span transition tunnel of Badaling Great Wall Station of Beijing-Zhangjiakou high-speed railway.An intelligent microseismic monitoring system is built with symmetry monitoring point layout both on the mountain surface and inside the tunnel to achieve three-dimensional and all-round monitoring results.Findings–Microseismic events can be divided into high density area,medium density area and low density area according to the density distribution of microseismic events.The positions where the cumulative distribution frequencies of microseismic events are 60 and 80%are identified as the boundaries between high and medium density areas and between medium and low density areas,respectively.The high density area of microseismic events is regarded as the high excavation damage zone of surrounding rock,which is affected by the grade of surrounding rock and the span of tunnel.The prediction formulas for the depth of high excavation damage zone of surrounding rock at different tunnel positions are given considering these two parameters.The scale of the average moment magnitude parameters of microseismic events is adopted to describe the damage degree of surrounding rock.The strong positive correlation and multistage characteristics between the depth of excavation damage zone and deformation of surrounding rock are revealed.Based on the depth of high excavation damage zone of surrounding rock,the prestressed anchor cable(rod)is designed,and the safety of anchor cable(rod)design parameters is verified by the deformation results of surrounding rock.Originality/value–The research provides a new method to predict the surrounding rock damage zone of large-span tunnel and also provides a reference basis for design parameters of prestressed anchor cable(rod).展开更多
Excavation of underground caverns,such as mountain tunnels and energy-storage caverns,may cause the damages to the surrounding rock as a result of the stress redistribution.In this influenced zone,new cracks and disco...Excavation of underground caverns,such as mountain tunnels and energy-storage caverns,may cause the damages to the surrounding rock as a result of the stress redistribution.In this influenced zone,new cracks and discontinuities are created or propagate in the rock mass.Therefore,it is effective to measure and evaluate the acoustic emission(AE)events generated by the rocks,which is a small elastic vibration,and permeability change.The authors have developed a long-term measurement device that incorporates an optical AE(O-AE)sensor,an optical pore pressure sensor,and an optical temperature sensor in a single multi-optical measurement probe(MOP).Japan Atomic Energy Agency has been conducting R&D activities to enhance the reliability of high-level radioactive waste(HLW)deep geological disposal technology.In a high-level radioactive disposal project,one of the challenges is the development of methods for long-term monitoring of rock mass behavior.Therefore,in January 2014,the long-term measurements of the hydro-mechanical behavior of the rock mass were launched using the developed MOP in the vicinity of 350 m below the surface at the Horonobe Underground Research Center.The measurement results show that AEs occur frequently up to 1.5 m from the wall during excavation.In addition,hydraulic conductivity increased by 2e4 orders of magnitude.Elastoplastic analysis revealed that the hydraulic behavior of the rock mass affected the pore pressure fluctuations and caused micro-fractures.Based on this,a conceptual model is developed to represent the excavation damaged zone(EDZ),which contributes to the safe geological disposal of radioactive waste.展开更多
Underground excavation can lead to stress redistribution and result in an excavation damaged zone(EDZ),which is an important factor affecting the excavation stability and support design.Accurately estimating the thick...Underground excavation can lead to stress redistribution and result in an excavation damaged zone(EDZ),which is an important factor affecting the excavation stability and support design.Accurately estimating the thickness of EDZ is essential to ensure the safety of the underground excavation.In this study,four novel hybrid ensemble learning models were developed by optimizing the extreme gradient boosting(XGBoost)and random forest(RF)algorithms through simulated annealing(SA)and Bayesian optimization(BO)approaches,namely SA-XGBoost,SA-RF,BO XGBoost and BO-RF models.A total of 210 cases were collected from Xiangxi Gold Mine in Hunan Province and Fankou Lead-zinc Mine in Guangdong Province,China,including seven input indicators:embedding depth,drift span,uniaxial compressive strength of rock,rock mass rating,unit weight of rock,lateral pressure coefficient of roadway and unit consumption of blasting explosive.The performance of the proposed models was evaluated by the coefficient of determination,root mean squared error,mean absolute error and variance accounted for.The results indicated that the SA-XGBoost model performed best.The Shapley additive explanations method revealed that the embedding depth was the most important indicator.Moreover,the convergence curves suggested that the SA-XGBoost model can reduce the generalization error and avoid overfitting.展开更多
Deep underground excavations within hard rocks can result in damage to the surrounding rock mass mostly due to redistribution of stresses.Especially within rock masses with non-persistent joints,the role of the pre-ex...Deep underground excavations within hard rocks can result in damage to the surrounding rock mass mostly due to redistribution of stresses.Especially within rock masses with non-persistent joints,the role of the pre-existing joints in the damage evolution around the underground opening is of critical importance as they govern the fracturing mechanisms and influence the brittle responses of these hard rock masses under highly anisotropic in situ stresses.In this study,the main focus is the impact of joint network geometry,joint strength and applied field stresses on the rock mass behaviours and the evolution of excavation induced damage due to the loss of confinement as a tunnel face advances.Analysis of such a phenomenon was conducted using the finite-discrete element method(FDEM).The numerical model is initially calibrated in order to match the behaviour of the fracture-free,massive Lac du Bonnet granite during the excavation of the Underground Research Laboratory(URL)Test Tunnel,Canada.The influence of the pre-existing joints on the rock mass response during excavation is investigated by integrating discrete fracture networks(DFNs)of various characteristics into the numerical models under varying in situ stresses.The numerical results obtained highlight the significance of the pre-existing joints on the reduction of in situ rock mass strength and its capacity for extension with both factors controlling the brittle response of the material.Furthermore,the impact of spatial distribution of natural joints on the stability of an underground excavation is discussed,as well as the potentially minor influence of joint strength on the stress induced damage within joint systems of a non-persistent nature under specific conditions.Additionally,the in situ stress-joint network interaction is examined,revealing the complex fracturing mechanisms that may lead to uncontrolled fracture propagation that compromises the overall stability of an underground excavation.展开更多
This paper presents a comprehensive derivation of fracture process zone size which closely parallels similar work in fracture of metals and anisotropic solid, but is adapted to conrete. Some nonlinear mechanics models...This paper presents a comprehensive derivation of fracture process zone size which closely parallels similar work in fracture of metals and anisotropic solid, but is adapted to conrete. Some nonlinear mechanics models of concrete materials will be discussed by using uniaxial stress assumptions. For uniaxial stress assumption, energy model and fracture model will be presented for nonlinear softening models. Finally, we make a comparison with those models.展开更多
Excavation Damaged Zone(EDZ)scope is important for optimizing excavation and support schemes in deep underground caverns.However,accurately predicting the full EDZ scope within the surrounding rock masses of deep unde...Excavation Damaged Zone(EDZ)scope is important for optimizing excavation and support schemes in deep underground caverns.However,accurately predicting the full EDZ scope within the surrounding rock masses of deep underground caverns during excavation remains a pressing problem.This study presents a comprehensive EDZ scope prediction approach(CESPA)for the brittle surrounding rock masses of deep underground caverns by coupling numerical simulation with quantitative analysis of borehole wall images and ultrasonic test results.First,the changes in both P-velocity(V_(p))and joint distribution of the surrounding rock masses before and after excavation damage are captured using ultrasonic tests and borehole digital cameras.Second,the quality Q-parameters of the surrounding rock mass before and after excavation damage are preliminarily rated with the rock mass descriptions provided by borehole wall images,and the rock mass V_(p)-parameter values are determined according to the V_(p)-borehole depth curves.Third,the Q-parameter ratings are further finely adjusted by updating the related Q-values to be similar with the Q-values estimated by V_(p)-parameter values.Fourth,the initial and residual mechanical parameters for the rock mass deterioration model(RDM)are estimated by the adjusted Q-parameter ratings based on the modified Q-based relations,and the elastic modulus deterioration index(EDI)threshold to describe the EDZ boundary is determined with the V_(p)-parameter values.Finally,EDZ scope is predicted using the elastoplastic numerical simulation with RDM and EDI based on the mechanical parameter estimates and EDI threshold.Analyses of applications in Sub-lab D1 in Jinping II project show that CESPA can provide a reliable and operable solution for predicting full EDZ scopes within the brittle surrounding rock masses of deep underground caverns.展开更多
Research on the permeability and pressure distribution characteristics of the roadway surrounding rock in the excavation damaged zone(EDZ) is beneficial for the development of gas control technology. In this study, an...Research on the permeability and pressure distribution characteristics of the roadway surrounding rock in the excavation damaged zone(EDZ) is beneficial for the development of gas control technology. In this study, analytical solutions of stress and strain of the roadway surrounding rock were obtained, in which the creep deformation and strain softening were considered. Using the MTS815 rock mechanics testing system and a gas permeability testing system, permeability tests were conducted in the complete stress-strain process, and the evolution characteristics of permeability and strain were studied over the whole loading process. Based on the analytical solutions of stress and strain and the governing equation of gas seepage flow, this paper proposes a hydro-mechanical(HM) model, which considers three different zones around the roadway. Then the gas flow process in the roadway surrounding rock in three different zones was simulated according to the engineering geological conditions, thus obtaining the permeability and pressure distribution characteristics of the roadway surrounding rock in three different zones. These results show that the surrounding rock around the roadway can be divided into four regions-the full flow zone(FFZ), flow-shielding zone(FSZ), transitive flow zone(TFZ), and in-situ rock flow zone(IRFZ). These results could provide theoretical guidance for the improvement of gas extraction and gas control technology.展开更多
Surrounding rock deterioration and large deformation have always been a significant difficulty in designing and constructing tunnels in soft rock.The key lies in real-time perception and quantitative assessment of the...Surrounding rock deterioration and large deformation have always been a significant difficulty in designing and constructing tunnels in soft rock.The key lies in real-time perception and quantitative assessment of the damaged area around the tunnel.An in situ microseismic(MS)monitoring system is established in the plateau soft tock tunnel.This technique facilitates spatiotemporal monitoring of the rock mass's fracturing expansion and squeezing deformation,which agree well with field convergence deformation results.The formation mechanisms of progressive failure evolution of soft rock tunnels were discussed and analyzed with MS data and numerical results.The results demonstrate that:(1)Localized stress concentration and layered rock result in significant asymmetry in micro-fractures propagation in the tunnel radial section.As excavation continues,the fracture extension area extends into the deep surrounding rockmass on the east side affected by the weak bedding;(2)Tunnel excavation and long-term deformation can induce tensile shear action on the rock mass,vertical tension fractures(account for 45%)exist in deep rockmass,which play a crucial role in controlling the macroscopic failure of surrounding rock;and(3)Based on the radiated MS energy,a three-dimensional model was created to visualize the damage zone of the tunnel surrounding rock.The model depicted varying degrees of damage,and three high damage zones were identified.Generally,the depth of high damage zone ranged from 4 m to 12 m.This study may be a valuable reference for the warning and controlling of large deformations in similar projects.展开更多
The strike-slip fault system in the central Tarim Craton controls a complex petroleum system with estimated reserves exceeding 1×10^(9)t,the fault-related fractures are important for hydrocarbon accumulation.In t...The strike-slip fault system in the central Tarim Craton controls a complex petroleum system with estimated reserves exceeding 1×10^(9)t,the fault-related fractures are important for hydrocarbon accumulation.In this paper,the basic parameters such as density and width of fractures are counted and classified,and the effects of fractures on reservoirs are analyzed.The results show that:(1)Structural fractures and stylolite were widely developed in Halahatang area and experienced at least three stages of activity based on the infilling materials and crosscutting relationship.(2)Fracture density,width,aperture,and dip angle vary in different wells,but the relationship between the above parameters and the distance to the fault core indicates the fracture differences in the fault damage zone and further provides a method to divide the inner units in the fault damage zone.In addition,oil and gas wells with high production mainly concentrate in the inner unit.(3)The infilling materials and degree of fractures vary.Fractures formed in the early stage are more filled and less open,while the fractures formed in the late stage are relatively less filled and more open.(4)Fractures improve porosity to a certain extent but greatly increase permeability,especially in the inner zone of fault damage zone with large quantity,multiple inclinations,less filling and large width.These features contribute to the formation of a higher-quality reservoir,further improving oil and gas production.This paper provides a quantitative characterization method for the study of strike-slip fault-related fracture-caved reservoirs,and points out that fault damage zone,especially the inner zone of the fault damage zone,is the potential goal for oil and gas exploration.展开更多
Field evidence has shown that large-scale and unstable discontinuous planes in the rock mass surrounding tunnels in rich water region are probably generated after excavation. The tunnel surrounding rock was divided in...Field evidence has shown that large-scale and unstable discontinuous planes in the rock mass surrounding tunnels in rich water region are probably generated after excavation. The tunnel surrounding rock was divided into three zones, including elastic zone, plastic damage zone and shear fracture zone for assessing the stability of the tunnel surrounding rock. By local hydrogeology, the stresses of surrounding rock of Jinshazhou circular tunnel was analyzed and the stress solutions on the elastic and plastic damage zones were obtained by applying the theories of fluid-solid coupling and elasto-plastic damage mechanics. The shear fracture zone generated by joints was studied and its range was determined by using Molar-Coulomb strength criterion. Finally, the correctness of the theoretical results was validated by comparing the scopes of shear fracture zones calculated in this paper with those from literature.展开更多
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.展开更多
After the excavation of the roadway,the original stress balance is destroyed,resulting in the redistribution of stress and the formation of an excavation damaged zone(EDZ)around the roadway.The thickness of EDZ is the...After the excavation of the roadway,the original stress balance is destroyed,resulting in the redistribution of stress and the formation of an excavation damaged zone(EDZ)around the roadway.The thickness of EDZ is the key basis for roadway stability discrimination and support structure design,and it is of great engineering significance to accurately predict the thickness of EDZ.Considering the advantages of machine learning(ML)in dealing with high-dimensional,nonlinear problems,a hybrid prediction model based on the random forest(RF)algorithm is developed in this paper.The model used the dragonfly algorithm(DA)to optimize two hyperparameters in RF,namely mtry and ntree,and used mean absolute error(MAE),rootmean square error(RMSE),determination coefficient(R^(2)),and variance accounted for(VAF)to evaluatemodel prediction performance.A database containing 217 sets of data was collected,with embedding depth(ED),drift span(DS),surrounding rock mass strength(RMS),joint index(JI)as input variables,and the excavation damaged zone thickness(EDZT)as output variable.In addition,four classic models,back propagation neural network(BPNN),extreme learning machine(ELM),radial basis function network(RBF),and RF were compared with the DA-RF model.The results showed that the DARF mold had the best prediction performance(training set:MAE=0.1036,RMSE=0.1514,R^(2)=0.9577,VAF=94.2645;test set:MAE=0.1115,RMSE=0.1417,R^(2)=0.9423,VAF=94.0836).The results of the sensitivity analysis showed that the relative importance of each input variable was DS,ED,RMS,and JI from low to high.展开更多
An abrupt damage model, taking full account of finite geometry changes, is used to study both the shapes of damage zones and the stress strain fields for a plane-strain tensile crack under the small-scale yielding con...An abrupt damage model, taking full account of finite geometry changes, is used to study both the shapes of damage zones and the stress strain fields for a plane-strain tensile crack under the small-scale yielding condition. Two typical crack-tip damage profiles are simulated by the element vanish technique. The fracture toughness increment due to damage dissipation is evaluated in terms of the energy release rate.展开更多
Degradation of rock mass produced by rock blasting,stress relief,and other causes is an important factor in the assessment of rock strength.Quantified as a disturbance factor,such degradation varies depending on blast...Degradation of rock mass produced by rock blasting,stress relief,and other causes is an important factor in the assessment of rock strength.Quantified as a disturbance factor,such degradation varies depending on blasting control,stress state and stress relief,and rock mass quality.This study focuses on the impact of disturbance on the safety of slopes.The disturbance in the rock mass is characterized by the geometry of the disturbed zone,its size,the magnitude,and the decaying rate with the distance away from the slope surface.A method accounting for decay of rock disturbance is presented.A study of the impact of rock disturbance characteristics on the quantitative stability measures of slopes was carried out.These characteristics included disturbed zone geometry,its thickness,the maximum magnitude of the disturbance factor,and the rate of disturbance decaying.The thickness of the disturbed zone and the maximum factor of disturbance were found to have the greatest impact.For example,the factor of safety for a 45slope in low-quality rock mass can decrease from 1.96 to 1.09 as the thickness of the disturbed zone increases from 1/4 of slope height H to the double of H and the maximum disturbance factor increases from 0.5 to 1.Uniform thickness of a disturbed zone was found to yield more conservative outcomes than the triangular zones did.The critical failure surfaces were found to be shallow for high rates of disturbance decay,and they were the deepest for spatially uniform disturbance factors.展开更多
Using pressure-preserved coring technique to determine in-situ gas content provides a more precise assessment of gas resource reserves and safeguard of mining safety in coal seams. How coring technique and depth affec...Using pressure-preserved coring technique to determine in-situ gas content provides a more precise assessment of gas resource reserves and safeguard of mining safety in coal seams. How coring technique and depth affect the determination of gas content is unclear due to borehole zoning rupture caused by roadway excavation and drilling disturbance. To this end, a proposed coupling model of stress distribution and gas migration was simulated and validated by FLAC^(3D) and COMSOL Multiphysics considering superposition effects of roadway excavation and drilling disturbance. The findings indicate that the roadway surrounding rock displays distinct zoning features including stress relief zone, stress concentration zone that is composed of plastic zone, elastic zone, and original stress zone;and the broken situations depending on the borehole peeping are consistent with the corresponding simulation results.On this basis, this study proposes a set of drilling coring depth calculation and prediction model for the gas desorption affected area under engineering disturbance. Optimal depth of coring drilling is not only approach to the in-situ coal bulk, but also can get the balance of the drilling workload and cost controlling. According to the typical mine site geological conditions and the numerical simulation results in this study, if the roadway excavation time is ~1 year, it is recommended that the pressure-preserved coring depth should be greater than 17 m.展开更多
In drill and blast method,due to uncontrolled blasting operations,the blast-induced damaged zone(BIDZ)is formed,whose mechanical and hydraulic properties are altered.This zone affects the behavior of the rock mass suc...In drill and blast method,due to uncontrolled blasting operations,the blast-induced damaged zone(BIDZ)is formed,whose mechanical and hydraulic properties are altered.This zone affects the behavior of the rock mass such that it reduces the strength of surrounding mass,and stability of the excavation.On the other hand,the groundwater is also effected by damaged zone induced stress and displacement,leading to a change in hydraulic flow around the tunnel which subsequently could produce new stress and displacement fields.In this research,an analytical solution for evaluating the stress and displacement of deep circular tunnels in elastoplastic rock mass is proposed,assuming the presence of BIDZ and hydraulic flow around the tunnel.The tunnel is subjected to in situ hydrostatic stresses,under radial hydraulic flow,and the damaged zone is supposed cylindrical shaped surrounding the cavity.Four different scenarios are predicted for stress evolution around the cavity considering the seepage zone,damage zone and plastic zone spread for elastic brittle-plastic behavior of surrounding mass.The analytical solution is validated using FLAC software,which shows excellent agreement.Examples are given to investigate the effect of BIDZ on the stress and displacement fields around the tunnel in both drained and undrained condition.The results show a significant impact on tunnel wall displacement especially for small magnitude of the ratio of seepage zone to damage zone radii,indicating its great significance in tunnel practice in terms of support and ground control.展开更多
基金partly supported by National Natural Science Foundation of China(Grant No.41472103)Technology Major Project(2016ZX05004001)
文摘Understanding the scaling relation of damage zone width with displacement of faults is important for predicting subsurface faulting mechanisms and fluid flow processes. The understanding of this scaling relationship is influenced by the accuracy of the methods and types of data utilized to investigate faults. In this study, seismic reflection data are used to investigate the throw and damage zone width of five strike-slip faults a ecting Ordovician carbonates of the Tarim intracraton basin,NW China. The results indicate that fault slips with a throw less than 200 m had formed wide damage zones up to 3000 m in width. Also, damage zone width is found to have both a positive correlation and a power-law relation with throw of two orders of magnitude, with a ratio of these values varying in a range of 2–15. The relationship between throw and damage zone width is not a simple power-law and changes its slope from small to larger size faults. The results indicate that throw scales well with damage zone width for the studied faults, and hence these can be used to predict fault geometries in the Tarim Basin. The study of the wide carbonate damage zones presented here provides new insights into scaling of large-size faults, which involve multiple faulting stages.
基金funded by the Natural Sciences and Engineering Research Council of Canadaby the Nuclear Waste Management Organization(NWMO)of Canada
文摘During the construction of an underground excavation, damage occurs in the surrounding rock mass due in large part to stress changes. While the predicted damage extent impacts profile selection and support design, the depth of damage is a critical aspect for the design of permeability sensitive excavations, such as a deep geological repository(DGR) for nuclear waste. Review of literature regarding the depth of excavation damage zones(EDZs) indicates three zones are common and typically related to stress induced damage. Based on past developments related to brittle damage prediction using continuum modelling, the depth of the EDZs has been examined numerically. One method to capture stress induced damage in conventional engineering software is the damage initiation and spalling limit(DISL) approach. The variability of depths predicted using the DISL approach has been evaluated and guidelines are suggested for determining the depth of the EDZs around circular excavations in brittle rock masses. Of the inputs evaluated, it was found that the tensile strength produces the greatest variation in the depth of the EDZs. The results were evaluated statistically to determine the best fit relation between the model inputs and the depth of the EDZs. The best correlation and least variation were found for the outer EDZ and the highly damaged zone(HDZ) showed the greatest variation. Predictive equations for different EDZs have been suggested and the maximum numerical EDZ depths, represented by the 68% prediction interval, agreed well with the empirical evidence. This suggests that the numerical limits can be used for preliminary depth prediction of the EDZs in brittle rock for circular excavations.
文摘The objective of this paper is to develop a methodology for calibration of a discrete element grain-based model(GBM)to replicate the hydro-mechanical properties of a brittle rock measured in the laboratory,and to apply the calibrated model to simulating the formation of excavation damage zone(EDZ)around underground excavations.Firstly,a new cohesive crack model is implemented into the universal distinct element code(UDEC)to control the fracturing behaviour of materials under various loading modes.Next,a methodology for calibration of the components of the UDEC-Voronoi model is discussed.The role of connectivity of induced microcracks on increasing the permeability of laboratory-scale samples is investigated.The calibrated samples are used to investigate the influence of pore fluid pressure on weakening the drained strength of the laboratory-scale rock.The validity of the Terzaghi’s effective stress law for the drained peak strength of low-porosity rock is tested by performing a series of biaxial compression test simulations.Finally,the evolution of damage and pore pressure around two unsupported circular tunnels in crystalline granitic rock is studied.
基金supported by the National Key R&D Program of China(No.2022YFF0609300)the National Major Science and Technology Projects of China(J2019-VI-0021-0137).
文摘In the creep fatigue crack growth of GH4169 alloy,oxidation is a prominent damage source,which is mainly manifested as the oxidation damage zone in front of crack tip.In order to investigate the property of the oxidation damage zone formed in the creep fatigue crack growth,crack growth tests of directly aged GH4169 alloy were conducted at 650℃ in air under various load conditions.Interrupted tests were performed to observe the damage characteristics at crack tip.Block tests were systematically executed to quantify the dependency of oxidation damage zone size on load and holding time.The crack propagation of the GH4169 alloy has a close relationship with grain boundary oxidation at 650℃.An oxidation damage zone in front of crack tip includes intergranular microcracks and oxidised but uncracked grain boundaries.Its size has been calculated from transient crack growth rate and described as a function of maximum stress intensity factor and holding time.Based on oxidation damage zone size,a novel model has been developed to predict the creep fatigue crack growth rate of the GH4169 alloy at 650℃.
基金supported by the Fundamental Research Funds for Chinese National Natural Science Foundation under Grant 51678035National Key Research and Development Programs of China under Grant 2017YFC0805401China Railway Corporation Research and Development Program of Science and Technology under Grant 2014004-C.
文摘Purpose–The microseismic monitoring technique has great advantages on identifying the location,extent and the mechanism of damage process occurring in rock mass.This study aims to analyze distribution characteristics and the evolution law of excavation damage zone of surrounding rock based on microseismic monitoring data.Design/methodology/approach–In situ test using microseismic monitoring technique is carried out in the large-span transition tunnel of Badaling Great Wall Station of Beijing-Zhangjiakou high-speed railway.An intelligent microseismic monitoring system is built with symmetry monitoring point layout both on the mountain surface and inside the tunnel to achieve three-dimensional and all-round monitoring results.Findings–Microseismic events can be divided into high density area,medium density area and low density area according to the density distribution of microseismic events.The positions where the cumulative distribution frequencies of microseismic events are 60 and 80%are identified as the boundaries between high and medium density areas and between medium and low density areas,respectively.The high density area of microseismic events is regarded as the high excavation damage zone of surrounding rock,which is affected by the grade of surrounding rock and the span of tunnel.The prediction formulas for the depth of high excavation damage zone of surrounding rock at different tunnel positions are given considering these two parameters.The scale of the average moment magnitude parameters of microseismic events is adopted to describe the damage degree of surrounding rock.The strong positive correlation and multistage characteristics between the depth of excavation damage zone and deformation of surrounding rock are revealed.Based on the depth of high excavation damage zone of surrounding rock,the prestressed anchor cable(rod)is designed,and the safety of anchor cable(rod)design parameters is verified by the deformation results of surrounding rock.Originality/value–The research provides a new method to predict the surrounding rock damage zone of large-span tunnel and also provides a reference basis for design parameters of prestressed anchor cable(rod).
文摘Excavation of underground caverns,such as mountain tunnels and energy-storage caverns,may cause the damages to the surrounding rock as a result of the stress redistribution.In this influenced zone,new cracks and discontinuities are created or propagate in the rock mass.Therefore,it is effective to measure and evaluate the acoustic emission(AE)events generated by the rocks,which is a small elastic vibration,and permeability change.The authors have developed a long-term measurement device that incorporates an optical AE(O-AE)sensor,an optical pore pressure sensor,and an optical temperature sensor in a single multi-optical measurement probe(MOP).Japan Atomic Energy Agency has been conducting R&D activities to enhance the reliability of high-level radioactive waste(HLW)deep geological disposal technology.In a high-level radioactive disposal project,one of the challenges is the development of methods for long-term monitoring of rock mass behavior.Therefore,in January 2014,the long-term measurements of the hydro-mechanical behavior of the rock mass were launched using the developed MOP in the vicinity of 350 m below the surface at the Horonobe Underground Research Center.The measurement results show that AEs occur frequently up to 1.5 m from the wall during excavation.In addition,hydraulic conductivity increased by 2e4 orders of magnitude.Elastoplastic analysis revealed that the hydraulic behavior of the rock mass affected the pore pressure fluctuations and caused micro-fractures.Based on this,a conceptual model is developed to represent the excavation damaged zone(EDZ),which contributes to the safe geological disposal of radioactive waste.
基金Project(52204117)supported by the National Natural Science Foundation of ChinaProject(2022JJ40601)supported by the Natural Science Foundation of Hunan Province,China。
文摘Underground excavation can lead to stress redistribution and result in an excavation damaged zone(EDZ),which is an important factor affecting the excavation stability and support design.Accurately estimating the thickness of EDZ is essential to ensure the safety of the underground excavation.In this study,four novel hybrid ensemble learning models were developed by optimizing the extreme gradient boosting(XGBoost)and random forest(RF)algorithms through simulated annealing(SA)and Bayesian optimization(BO)approaches,namely SA-XGBoost,SA-RF,BO XGBoost and BO-RF models.A total of 210 cases were collected from Xiangxi Gold Mine in Hunan Province and Fankou Lead-zinc Mine in Guangdong Province,China,including seven input indicators:embedding depth,drift span,uniaxial compressive strength of rock,rock mass rating,unit weight of rock,lateral pressure coefficient of roadway and unit consumption of blasting explosive.The performance of the proposed models was evaluated by the coefficient of determination,root mean squared error,mean absolute error and variance accounted for.The results indicated that the SA-XGBoost model performed best.The Shapley additive explanations method revealed that the embedding depth was the most important indicator.Moreover,the convergence curves suggested that the SA-XGBoost model can reduce the generalization error and avoid overfitting.
基金the Natural Sciences and Engineering Research Council of Canadathe Ministry of National Defensethe RMC Green Team for providing the funding and the resources
文摘Deep underground excavations within hard rocks can result in damage to the surrounding rock mass mostly due to redistribution of stresses.Especially within rock masses with non-persistent joints,the role of the pre-existing joints in the damage evolution around the underground opening is of critical importance as they govern the fracturing mechanisms and influence the brittle responses of these hard rock masses under highly anisotropic in situ stresses.In this study,the main focus is the impact of joint network geometry,joint strength and applied field stresses on the rock mass behaviours and the evolution of excavation induced damage due to the loss of confinement as a tunnel face advances.Analysis of such a phenomenon was conducted using the finite-discrete element method(FDEM).The numerical model is initially calibrated in order to match the behaviour of the fracture-free,massive Lac du Bonnet granite during the excavation of the Underground Research Laboratory(URL)Test Tunnel,Canada.The influence of the pre-existing joints on the rock mass response during excavation is investigated by integrating discrete fracture networks(DFNs)of various characteristics into the numerical models under varying in situ stresses.The numerical results obtained highlight the significance of the pre-existing joints on the reduction of in situ rock mass strength and its capacity for extension with both factors controlling the brittle response of the material.Furthermore,the impact of spatial distribution of natural joints on the stability of an underground excavation is discussed,as well as the potentially minor influence of joint strength on the stress induced damage within joint systems of a non-persistent nature under specific conditions.Additionally,the in situ stress-joint network interaction is examined,revealing the complex fracturing mechanisms that may lead to uncontrolled fracture propagation that compromises the overall stability of an underground excavation.
文摘This paper presents a comprehensive derivation of fracture process zone size which closely parallels similar work in fracture of metals and anisotropic solid, but is adapted to conrete. Some nonlinear mechanics models of concrete materials will be discussed by using uniaxial stress assumptions. For uniaxial stress assumption, energy model and fracture model will be presented for nonlinear softening models. Finally, we make a comparison with those models.
基金funded by the National Natural Science Foundation of China under Grant Nos.51979268,U1765206,41877256。
文摘Excavation Damaged Zone(EDZ)scope is important for optimizing excavation and support schemes in deep underground caverns.However,accurately predicting the full EDZ scope within the surrounding rock masses of deep underground caverns during excavation remains a pressing problem.This study presents a comprehensive EDZ scope prediction approach(CESPA)for the brittle surrounding rock masses of deep underground caverns by coupling numerical simulation with quantitative analysis of borehole wall images and ultrasonic test results.First,the changes in both P-velocity(V_(p))and joint distribution of the surrounding rock masses before and after excavation damage are captured using ultrasonic tests and borehole digital cameras.Second,the quality Q-parameters of the surrounding rock mass before and after excavation damage are preliminarily rated with the rock mass descriptions provided by borehole wall images,and the rock mass V_(p)-parameter values are determined according to the V_(p)-borehole depth curves.Third,the Q-parameter ratings are further finely adjusted by updating the related Q-values to be similar with the Q-values estimated by V_(p)-parameter values.Fourth,the initial and residual mechanical parameters for the rock mass deterioration model(RDM)are estimated by the adjusted Q-parameter ratings based on the modified Q-based relations,and the elastic modulus deterioration index(EDI)threshold to describe the EDZ boundary is determined with the V_(p)-parameter values.Finally,EDZ scope is predicted using the elastoplastic numerical simulation with RDM and EDI based on the mechanical parameter estimates and EDI threshold.Analyses of applications in Sub-lab D1 in Jinping II project show that CESPA can provide a reliable and operable solution for predicting full EDZ scopes within the brittle surrounding rock masses of deep underground caverns.
基金financially supported by the Natural Science Foundation of Jiangsu Province,China(No.BK20140189)the Postdoctoral Science Foundation of China(No.2014M550315)
文摘Research on the permeability and pressure distribution characteristics of the roadway surrounding rock in the excavation damaged zone(EDZ) is beneficial for the development of gas control technology. In this study, analytical solutions of stress and strain of the roadway surrounding rock were obtained, in which the creep deformation and strain softening were considered. Using the MTS815 rock mechanics testing system and a gas permeability testing system, permeability tests were conducted in the complete stress-strain process, and the evolution characteristics of permeability and strain were studied over the whole loading process. Based on the analytical solutions of stress and strain and the governing equation of gas seepage flow, this paper proposes a hydro-mechanical(HM) model, which considers three different zones around the roadway. Then the gas flow process in the roadway surrounding rock in three different zones was simulated according to the engineering geological conditions, thus obtaining the permeability and pressure distribution characteristics of the roadway surrounding rock in three different zones. These results show that the surrounding rock around the roadway can be divided into four regions-the full flow zone(FFZ), flow-shielding zone(FSZ), transitive flow zone(TFZ), and in-situ rock flow zone(IRFZ). These results could provide theoretical guidance for the improvement of gas extraction and gas control technology.
基金the funding support from the National Natural Science Foundation of China(Grant Nos.U23A2060,42177143 and 42277461).
文摘Surrounding rock deterioration and large deformation have always been a significant difficulty in designing and constructing tunnels in soft rock.The key lies in real-time perception and quantitative assessment of the damaged area around the tunnel.An in situ microseismic(MS)monitoring system is established in the plateau soft tock tunnel.This technique facilitates spatiotemporal monitoring of the rock mass's fracturing expansion and squeezing deformation,which agree well with field convergence deformation results.The formation mechanisms of progressive failure evolution of soft rock tunnels were discussed and analyzed with MS data and numerical results.The results demonstrate that:(1)Localized stress concentration and layered rock result in significant asymmetry in micro-fractures propagation in the tunnel radial section.As excavation continues,the fracture extension area extends into the deep surrounding rockmass on the east side affected by the weak bedding;(2)Tunnel excavation and long-term deformation can induce tensile shear action on the rock mass,vertical tension fractures(account for 45%)exist in deep rockmass,which play a crucial role in controlling the macroscopic failure of surrounding rock;and(3)Based on the radiated MS energy,a three-dimensional model was created to visualize the damage zone of the tunnel surrounding rock.The model depicted varying degrees of damage,and three high damage zones were identified.Generally,the depth of high damage zone ranged from 4 m to 12 m.This study may be a valuable reference for the warning and controlling of large deformations in similar projects.
基金supported by the Natural Science Foundation of China-Youth Foundation(42402163)Natural Science Foundation of Sichuan Province of China(2024NSFSC0814)Science and Technology Cooperation Project of the CNPC-SWPU Innovation Alliance(2020CX010101).
文摘The strike-slip fault system in the central Tarim Craton controls a complex petroleum system with estimated reserves exceeding 1×10^(9)t,the fault-related fractures are important for hydrocarbon accumulation.In this paper,the basic parameters such as density and width of fractures are counted and classified,and the effects of fractures on reservoirs are analyzed.The results show that:(1)Structural fractures and stylolite were widely developed in Halahatang area and experienced at least three stages of activity based on the infilling materials and crosscutting relationship.(2)Fracture density,width,aperture,and dip angle vary in different wells,but the relationship between the above parameters and the distance to the fault core indicates the fracture differences in the fault damage zone and further provides a method to divide the inner units in the fault damage zone.In addition,oil and gas wells with high production mainly concentrate in the inner unit.(3)The infilling materials and degree of fractures vary.Fractures formed in the early stage are more filled and less open,while the fractures formed in the late stage are relatively less filled and more open.(4)Fractures improve porosity to a certain extent but greatly increase permeability,especially in the inner zone of fault damage zone with large quantity,multiple inclinations,less filling and large width.These features contribute to the formation of a higher-quality reservoir,further improving oil and gas production.This paper provides a quantitative characterization method for the study of strike-slip fault-related fracture-caved reservoirs,and points out that fault damage zone,especially the inner zone of the fault damage zone,is the potential goal for oil and gas exploration.
基金The Program for New Century Excellent Talents in University ( No. NCET-06-0649)the Natural Science Foundation of Hubei Province (No.2005ABA303)
文摘Field evidence has shown that large-scale and unstable discontinuous planes in the rock mass surrounding tunnels in rich water region are probably generated after excavation. The tunnel surrounding rock was divided into three zones, including elastic zone, plastic damage zone and shear fracture zone for assessing the stability of the tunnel surrounding rock. By local hydrogeology, the stresses of surrounding rock of Jinshazhou circular tunnel was analyzed and the stress solutions on the elastic and plastic damage zones were obtained by applying the theories of fluid-solid coupling and elasto-plastic damage mechanics. The shear fracture zone generated by joints was studied and its range was determined by using Molar-Coulomb strength criterion. Finally, the correctness of the theoretical results was validated by comparing the scopes of shear fracture zones calculated in this paper with those from literature.
基金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.
基金funded by the National Science Foundation of China(42177164)the Distinguished Youth Science Foundation of Hunan Province of China(2022JJ10073)the Innovation-Driven Project of Central South University(2020CX040).
文摘After the excavation of the roadway,the original stress balance is destroyed,resulting in the redistribution of stress and the formation of an excavation damaged zone(EDZ)around the roadway.The thickness of EDZ is the key basis for roadway stability discrimination and support structure design,and it is of great engineering significance to accurately predict the thickness of EDZ.Considering the advantages of machine learning(ML)in dealing with high-dimensional,nonlinear problems,a hybrid prediction model based on the random forest(RF)algorithm is developed in this paper.The model used the dragonfly algorithm(DA)to optimize two hyperparameters in RF,namely mtry and ntree,and used mean absolute error(MAE),rootmean square error(RMSE),determination coefficient(R^(2)),and variance accounted for(VAF)to evaluatemodel prediction performance.A database containing 217 sets of data was collected,with embedding depth(ED),drift span(DS),surrounding rock mass strength(RMS),joint index(JI)as input variables,and the excavation damaged zone thickness(EDZT)as output variable.In addition,four classic models,back propagation neural network(BPNN),extreme learning machine(ELM),radial basis function network(RBF),and RF were compared with the DA-RF model.The results showed that the DARF mold had the best prediction performance(training set:MAE=0.1036,RMSE=0.1514,R^(2)=0.9577,VAF=94.2645;test set:MAE=0.1115,RMSE=0.1417,R^(2)=0.9423,VAF=94.0836).The results of the sensitivity analysis showed that the relative importance of each input variable was DS,ED,RMS,and JI from low to high.
文摘An abrupt damage model, taking full account of finite geometry changes, is used to study both the shapes of damage zones and the stress strain fields for a plane-strain tensile crack under the small-scale yielding condition. Two typical crack-tip damage profiles are simulated by the element vanish technique. The fracture toughness increment due to damage dissipation is evaluated in terms of the energy release rate.
基金supported by the National Science Foundation(Grant No.CMMI-1901582)supported by the Nuclear Research and Development Program of the National Research Foundation of Korea(Grant No.2024-M2E3A2007963)the Korea Electric Power Corporation(Grant No.R22XO05-05).
文摘Degradation of rock mass produced by rock blasting,stress relief,and other causes is an important factor in the assessment of rock strength.Quantified as a disturbance factor,such degradation varies depending on blasting control,stress state and stress relief,and rock mass quality.This study focuses on the impact of disturbance on the safety of slopes.The disturbance in the rock mass is characterized by the geometry of the disturbed zone,its size,the magnitude,and the decaying rate with the distance away from the slope surface.A method accounting for decay of rock disturbance is presented.A study of the impact of rock disturbance characteristics on the quantitative stability measures of slopes was carried out.These characteristics included disturbed zone geometry,its thickness,the maximum magnitude of the disturbance factor,and the rate of disturbance decaying.The thickness of the disturbed zone and the maximum factor of disturbance were found to have the greatest impact.For example,the factor of safety for a 45slope in low-quality rock mass can decrease from 1.96 to 1.09 as the thickness of the disturbed zone increases from 1/4 of slope height H to the double of H and the maximum disturbance factor increases from 0.5 to 1.Uniform thickness of a disturbed zone was found to yield more conservative outcomes than the triangular zones did.The critical failure surfaces were found to be shallow for high rates of disturbance decay,and they were the deepest for spatially uniform disturbance factors.
基金supported by National Natural Science Foundation of China (Nos. 51827901, 52104096)the Shenzhen National Science Fund for Distinguished Young Scholars (No. RCJC20210706091948015)+2 种基金Open Fund of State Key Laboratory of Hydraulics and Mountain River Engineering, Sichuan University (No. SKHL2216)National Key R & D Program of China (No. 2022YFB3706605)Youth Foundation of Sichuan Natural Science Foundation (No. 2023NSFSC0780)。
文摘Using pressure-preserved coring technique to determine in-situ gas content provides a more precise assessment of gas resource reserves and safeguard of mining safety in coal seams. How coring technique and depth affect the determination of gas content is unclear due to borehole zoning rupture caused by roadway excavation and drilling disturbance. To this end, a proposed coupling model of stress distribution and gas migration was simulated and validated by FLAC^(3D) and COMSOL Multiphysics considering superposition effects of roadway excavation and drilling disturbance. The findings indicate that the roadway surrounding rock displays distinct zoning features including stress relief zone, stress concentration zone that is composed of plastic zone, elastic zone, and original stress zone;and the broken situations depending on the borehole peeping are consistent with the corresponding simulation results.On this basis, this study proposes a set of drilling coring depth calculation and prediction model for the gas desorption affected area under engineering disturbance. Optimal depth of coring drilling is not only approach to the in-situ coal bulk, but also can get the balance of the drilling workload and cost controlling. According to the typical mine site geological conditions and the numerical simulation results in this study, if the roadway excavation time is ~1 year, it is recommended that the pressure-preserved coring depth should be greater than 17 m.
文摘In drill and blast method,due to uncontrolled blasting operations,the blast-induced damaged zone(BIDZ)is formed,whose mechanical and hydraulic properties are altered.This zone affects the behavior of the rock mass such that it reduces the strength of surrounding mass,and stability of the excavation.On the other hand,the groundwater is also effected by damaged zone induced stress and displacement,leading to a change in hydraulic flow around the tunnel which subsequently could produce new stress and displacement fields.In this research,an analytical solution for evaluating the stress and displacement of deep circular tunnels in elastoplastic rock mass is proposed,assuming the presence of BIDZ and hydraulic flow around the tunnel.The tunnel is subjected to in situ hydrostatic stresses,under radial hydraulic flow,and the damaged zone is supposed cylindrical shaped surrounding the cavity.Four different scenarios are predicted for stress evolution around the cavity considering the seepage zone,damage zone and plastic zone spread for elastic brittle-plastic behavior of surrounding mass.The analytical solution is validated using FLAC software,which shows excellent agreement.Examples are given to investigate the effect of BIDZ on the stress and displacement fields around the tunnel in both drained and undrained condition.The results show a significant impact on tunnel wall displacement especially for small magnitude of the ratio of seepage zone to damage zone radii,indicating its great significance in tunnel practice in terms of support and ground control.