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.展开更多
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.展开更多
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.展开更多
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.展开更多
The aim of the present research was to establish a case study for the prediction of the unknown EDZ(Excavation Damaged Zone)distribution using a numerical analysis calibrated by replicating the trends in the EDZ obser...The aim of the present research was to establish a case study for the prediction of the unknown EDZ(Excavation Damaged Zone)distribution using a numerical analysis calibrated by replicating the trends in the EDZ observed from one of the representative underground research fields in Japan(Horonobe URL).In this study,a 2D numerical analysis using a damage model,which can determine rock deformation and fracturing simultaneously,is presented.It was calibrated to reproduce the excavation of the gallery at the Horonobe URL at a depth of 350 m.Simulated results show an excellent agreement with the extent of the measured EDZ and capture the failure modes of EDZ fractures suggested by the in-situ observations.Finally,the calibrated numerical analysis was used to realistically estimate the EDZ formation for the geological disposal of high-level radioactive waste(HLW)under the same environment as that of the above-mentioned galley at the Horonobe URL.Consequently,it was shown that the tensile/shear hybrid fractures dominantly constituted the EDZ and propagated to a maximum extent of about 0.3 m from the cavity wall during the cavity excavation for the HLW disposal.Overall,the calibrated numerical analysis and resulting estimations,targeted for the environment at the depth of 350 m at the Horonobe URL,where mudstone is located,should be useful for predicting the trends in the EDZ distribution expected in the implementation of HLW disposal projects under deep geological conditions,such as those that exist in Japan,which are dominated by sedimentary rocks,including mudstone。展开更多
Subsurface excavation results in the formation of a zone called excavation damaged zone(EDZ)around the tunnel wall.An EDZ is a major concern in the field of high-level radioactive waste disposal because it may act as ...Subsurface excavation results in the formation of a zone called excavation damaged zone(EDZ)around the tunnel wall.An EDZ is a major concern in the field of high-level radioactive waste disposal because it may act as a flow path after the closure of a repository.In this study,first-arrival traveltime tomography was repeatedly conducted on the EDZ at a depth of 350 m in the Horonobe Underground Research Laboratory.However,the acquired data was highly affected by the support structure on the drift wall.For proper visualization of the EDZ,information about the structure was incorporated into the inversion by modifying the model constraint.The synthetic study showed that the approach reproduced the EDZ in the model without the artifacts.The method was applied to field data,and the EDZ around the drift was detected.The inversion was extended to a time-lapse inversion to trace the changes in P-wave velocity in the EDZ.The synthetic study demonstrated that temporal changes in the P-wave velocity distribution could be detected.Data obtained from 12 surveys under open-drift conditions were analyzed by time-lapse inversion.The results indicated that the EDZ did not undergo sealing or evolution at the site for approximately seven years.展开更多
Over the past twenty years, there has been a growing interest in the development of numerical modelsthat can realistically capture the progressive failure of rock masses. In particular, the investigation ofdamage deve...Over the past twenty years, there has been a growing interest in the development of numerical modelsthat can realistically capture the progressive failure of rock masses. In particular, the investigation ofdamage development around underground excavations represents a key issue in several rock engineeringapplications, including tunnelling, mining, drilling, hydroelectric power generation, and the deepgeological disposal of nuclear waste. The goal of this paper is to show the effectiveness of a hybrid finitediscreteelement method (FDEM) code to simulate the fracturing mechanisms associated with theexcavation of underground openings in brittle rock formations. A brief review of the current state-of-theartmodelling approaches is initially provided, including the description of selecting continuum- anddiscontinuum-based techniques. Then, the influence of a number of factors, including mechanical and insitu stress anisotropy, as well as excavation geometry, on the simulated damage is analysed for threedifferent geomechanical scenarios. Firstly, the fracture nucleation and growth process under isotropicrock mass conditions is simulated for a circular shaft. Secondly, the influence of mechanical anisotropy onthe development of an excavation damaged zone (EDZ) around a tunnel excavated in a layered rockformation is considered. Finally, the interaction mechanisms between two large caverns of an undergroundhydroelectric power station are investigated, with particular emphasis on the rock mass responsesensitivity to the pillar width and excavation sequence. Overall, the numerical results indicate that FDEMsimulations can provide unique geomechanical insights in cases where an explicit consideration offracture and fragmentation processes is of paramount importance. 2014 Institute of Rock and Soil Mechanics, Chinese Academy of Sciences. Production and hosting byElsevier B.V. All rights reserved.展开更多
Four 16.7 km-long tunnels with diameters ranging from 12.4 to 14.6 m are now under construction at Jinping II hydropower station along the Yalong River.The tunnels pass through Triassic rocks below Jinping Mountain.Th...Four 16.7 km-long tunnels with diameters ranging from 12.4 to 14.6 m are now under construction at Jinping II hydropower station along the Yalong River.The tunnels pass through Triassic rocks below Jinping Mountain.The tunnels are characterized with high overburden,long alignment and complex geological conditions.Brittle failure in marble and squeezing in schist are the primary problems in tunnelling.This paper introduces the studies of laboratory tests on Jinping II marble as well as numerical prediction of excavation damaged zone (EDZ) of tunnel section in brittle marble and determination of reinforced concrete lining thickness for restraining time-dependent deformation in the schist tunnel section.Laboratory tests indicate that Jinping II marble presents a complex brittle-ductile-plastic transition behavior of post-peak response with increasing confining pressure.Such behavior can be described numerically with the Hoek-Brown model.The EDZ was calibrated and predicted using both fast Lagrangian analysis of continua (FLAC) and particle flow code (PFC).The predicted result of EDZ in sections with different qualities of rock mass under various overburden pressures is quite helpful in understanding EDZ characterization and support design.A power-law creep model was used to support the lining design,especially in determining the lining thickness.Field convergence measurement data over 19 months were used to calibrate the creep model properties,followed by a sensibility analysis of reinforced concrete lining thickness that was launched to present the maximum lining compressive stress.展开更多
In this paper, the coupled thermo-mechanical (TM) processes in the AEspoe Pillar Stability Experiment (APSE) carried out by the Swedish Nuclear Fuel and Waste Management Company (SKB) were simulated using both c...In this paper, the coupled thermo-mechanical (TM) processes in the AEspoe Pillar Stability Experiment (APSE) carried out by the Swedish Nuclear Fuel and Waste Management Company (SKB) were simulated using both continuum and discontinuum based numerical methods. Two-dimensional (2D) and three- dimensional (3D) finite element method (FEM) and 2D distinct element method (DEM) with particles were used. The main objective for the large scale in situ experiment is to investigate the yielding strength of crystalline rock and the formation of the excavation disturbed/damaged zone (EDZ) during excavation of two boreholes, pressurizing of one of the boreholes and heating. For the DEM simulations, the heat flow algorithm was newly introduced into the original code. The calculated stress, displacement and temperature distributions were compared with the ones obtained from in situ measurements and FEM simulations. A parametric study for initial microcracks was also performed to reproduce the spalling phenomena observed in the APSE.展开更多
In the present study,a numerical model is first calibrated against the crack networks and pressure attenuation data in laboratory blasting test.Then,based on the calibrated numerical model,two-hole plane models are de...In the present study,a numerical model is first calibrated against the crack networks and pressure attenuation data in laboratory blasting test.Then,based on the calibrated numerical model,two-hole plane models are developed and used to perform a series of sim-ulations of smoothwall blasting in deep tunnelling subjected to in-situ stress.The evolutions of rock fracture and excavation damage zone in the roof/floor and sidewalls under different far-field hydrostatic pressure and anisotropic in-situ stress conditions are numerically investigated.The findings in numerical modelling are also analytically interpreted with the stress distributions around the designed tunnel perimeter and perimeter borehole.The numerical and analytical results show that the variations of rock cracking and excavation dam-aged zone induced by smoothwall blasting with in-situ stress are mainly attributed to the high tangential compressive stress concentration around the remaining rock after inner primary blasts and the tensile stress acting on the wall of perimeter hole,which control the crack propagation and initiation respectively.At last,the implications of findings for practical smoothwall blasting in deep tunnelling are discussed.展开更多
基金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 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.
基金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.
基金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.
基金supported by the JST FOREST Program(Grant no.JPMJFR216Y,Japan)JSPS KAKENHI(Grant nos.20K14826,22H01589,and 22H00229).
文摘The aim of the present research was to establish a case study for the prediction of the unknown EDZ(Excavation Damaged Zone)distribution using a numerical analysis calibrated by replicating the trends in the EDZ observed from one of the representative underground research fields in Japan(Horonobe URL).In this study,a 2D numerical analysis using a damage model,which can determine rock deformation and fracturing simultaneously,is presented.It was calibrated to reproduce the excavation of the gallery at the Horonobe URL at a depth of 350 m.Simulated results show an excellent agreement with the extent of the measured EDZ and capture the failure modes of EDZ fractures suggested by the in-situ observations.Finally,the calibrated numerical analysis was used to realistically estimate the EDZ formation for the geological disposal of high-level radioactive waste(HLW)under the same environment as that of the above-mentioned galley at the Horonobe URL.Consequently,it was shown that the tensile/shear hybrid fractures dominantly constituted the EDZ and propagated to a maximum extent of about 0.3 m from the cavity wall during the cavity excavation for the HLW disposal.Overall,the calibrated numerical analysis and resulting estimations,targeted for the environment at the depth of 350 m at the Horonobe URL,where mudstone is located,should be useful for predicting the trends in the EDZ distribution expected in the implementation of HLW disposal projects under deep geological conditions,such as those that exist in Japan,which are dominated by sedimentary rocks,including mudstone。
基金This study was conducted as part of a commissioned project entitled“The project for validating sealing of the geological repository(2020 and 2021 FY)”funded by the Ministry of Economy,Trade,and Industry of Japan.
文摘Subsurface excavation results in the formation of a zone called excavation damaged zone(EDZ)around the tunnel wall.An EDZ is a major concern in the field of high-level radioactive waste disposal because it may act as a flow path after the closure of a repository.In this study,first-arrival traveltime tomography was repeatedly conducted on the EDZ at a depth of 350 m in the Horonobe Underground Research Laboratory.However,the acquired data was highly affected by the support structure on the drift wall.For proper visualization of the EDZ,information about the structure was incorporated into the inversion by modifying the model constraint.The synthetic study showed that the approach reproduced the EDZ in the model without the artifacts.The method was applied to field data,and the EDZ around the drift was detected.The inversion was extended to a time-lapse inversion to trace the changes in P-wave velocity in the EDZ.The synthetic study demonstrated that temporal changes in the P-wave velocity distribution could be detected.Data obtained from 12 surveys under open-drift conditions were analyzed by time-lapse inversion.The results indicated that the EDZ did not undergo sealing or evolution at the site for approximately seven years.
基金supported by the Natural Science and Engineering Research Council (NSERC) of Canada in the form of discovery grant No. 341275the Swiss National Cooperative for the Disposal of Radioactive Waste (NAGRA)
文摘Over the past twenty years, there has been a growing interest in the development of numerical modelsthat can realistically capture the progressive failure of rock masses. In particular, the investigation ofdamage development around underground excavations represents a key issue in several rock engineeringapplications, including tunnelling, mining, drilling, hydroelectric power generation, and the deepgeological disposal of nuclear waste. The goal of this paper is to show the effectiveness of a hybrid finitediscreteelement method (FDEM) code to simulate the fracturing mechanisms associated with theexcavation of underground openings in brittle rock formations. A brief review of the current state-of-theartmodelling approaches is initially provided, including the description of selecting continuum- anddiscontinuum-based techniques. Then, the influence of a number of factors, including mechanical and insitu stress anisotropy, as well as excavation geometry, on the simulated damage is analysed for threedifferent geomechanical scenarios. Firstly, the fracture nucleation and growth process under isotropicrock mass conditions is simulated for a circular shaft. Secondly, the influence of mechanical anisotropy onthe development of an excavation damaged zone (EDZ) around a tunnel excavated in a layered rockformation is considered. Finally, the interaction mechanisms between two large caverns of an undergroundhydroelectric power station are investigated, with particular emphasis on the rock mass responsesensitivity to the pillar width and excavation sequence. Overall, the numerical results indicate that FDEMsimulations can provide unique geomechanical insights in cases where an explicit consideration offracture and fragmentation processes is of paramount importance. 2014 Institute of Rock and Soil Mechanics, Chinese Academy of Sciences. Production and hosting byElsevier B.V. All rights reserved.
基金Support by the China Postdoctor Fund (20100471757)
文摘Four 16.7 km-long tunnels with diameters ranging from 12.4 to 14.6 m are now under construction at Jinping II hydropower station along the Yalong River.The tunnels pass through Triassic rocks below Jinping Mountain.The tunnels are characterized with high overburden,long alignment and complex geological conditions.Brittle failure in marble and squeezing in schist are the primary problems in tunnelling.This paper introduces the studies of laboratory tests on Jinping II marble as well as numerical prediction of excavation damaged zone (EDZ) of tunnel section in brittle marble and determination of reinforced concrete lining thickness for restraining time-dependent deformation in the schist tunnel section.Laboratory tests indicate that Jinping II marble presents a complex brittle-ductile-plastic transition behavior of post-peak response with increasing confining pressure.Such behavior can be described numerically with the Hoek-Brown model.The EDZ was calibrated and predicted using both fast Lagrangian analysis of continua (FLAC) and particle flow code (PFC).The predicted result of EDZ in sections with different qualities of rock mass under various overburden pressures is quite helpful in understanding EDZ characterization and support design.A power-law creep model was used to support the lining design,especially in determining the lining thickness.Field convergence measurement data over 19 months were used to calibrate the creep model properties,followed by a sensibility analysis of reinforced concrete lining thickness that was launched to present the maximum lining compressive stress.
基金conducted within the context of the international DECOVALEX Project (DEvelopment of COupled models and their VALidation against EXperiments)financed by Japan Atomic Energy Agency (JAEA) who was also one of the Funding Organizations of the projectChrister Anders-son from Swedish Nuclear Fuel and Waste Management Co.(SKB),Sweden
文摘In this paper, the coupled thermo-mechanical (TM) processes in the AEspoe Pillar Stability Experiment (APSE) carried out by the Swedish Nuclear Fuel and Waste Management Company (SKB) were simulated using both continuum and discontinuum based numerical methods. Two-dimensional (2D) and three- dimensional (3D) finite element method (FEM) and 2D distinct element method (DEM) with particles were used. The main objective for the large scale in situ experiment is to investigate the yielding strength of crystalline rock and the formation of the excavation disturbed/damaged zone (EDZ) during excavation of two boreholes, pressurizing of one of the boreholes and heating. For the DEM simulations, the heat flow algorithm was newly introduced into the original code. The calculated stress, displacement and temperature distributions were compared with the ones obtained from in situ measurements and FEM simulations. A parametric study for initial microcracks was also performed to reproduce the spalling phenomena observed in the APSE.
基金the financial support from the National Natural Science Foundation of China(Grant no.51974360)for carrying out this research work.
文摘In the present study,a numerical model is first calibrated against the crack networks and pressure attenuation data in laboratory blasting test.Then,based on the calibrated numerical model,two-hole plane models are developed and used to perform a series of sim-ulations of smoothwall blasting in deep tunnelling subjected to in-situ stress.The evolutions of rock fracture and excavation damage zone in the roof/floor and sidewalls under different far-field hydrostatic pressure and anisotropic in-situ stress conditions are numerically investigated.The findings in numerical modelling are also analytically interpreted with the stress distributions around the designed tunnel perimeter and perimeter borehole.The numerical and analytical results show that the variations of rock cracking and excavation dam-aged zone induced by smoothwall blasting with in-situ stress are mainly attributed to the high tangential compressive stress concentration around the remaining rock after inner primary blasts and the tensile stress acting on the wall of perimeter hole,which control the crack propagation and initiation respectively.At last,the implications of findings for practical smoothwall blasting in deep tunnelling are discussed.