Predicting excavation damage zone depths in brittle rocks

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.

Seismic damage zone and width-throw scaling along the strike-slip faults in the Ordovician carbonates in the Tarim Basin

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.

Calibration of coupled hydro-mechanical properties of grain-based model for simulating fracture process and associated pore pressure evolution in excavation damage zone around deep tunnels

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.

Distribution characteristics and the evolution law of excavation damage zone in the large-span transition section of high-speed railway tunnel based on microseismic monitoring

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).

Evaluation of excavation damaged zones(EDZs)in Horonobe Underground Research Laboratory(URL)

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.

Oxidation damage zone formed in creep fatigue crack growth of GH4169 alloy at 650℃

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℃.

Predicting the Thickness of an Excavation Damaged Zone around the Roadway Using the DA-RF Hybrid Model

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.

THE DAMAGE PROCESS ZONE CHARACTERISTICS AT CRACK TIP IN CONCRETE

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.

Predicting the excavation damaged zone within brittle surrounding rock masses of deep underground caverns using a comprehensive approach integrating in situ measurements and numerical analysis

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.

In-situ observations of damage-fracture evolution in surrounding rock upon unloading in 2400-m-deep tunnels

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.

Stresses and Shear Fracture Zone of Jinshazhou Tunnel Surrounding Rock in Rich Water Region

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.

BRITTLE-DAMAGE ANALYSIS OF CRACKS UNDER CONDITIONS OF PLANE-STRAIN TENSILE LOADING

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.

Numerical modelling of spatially and temporally distributed on‑fault induced seismicity:implication for seismic hazards

Induced seismicity is strongly related to various engineering projects that cause anthropogenic in-situ stress change at a great depth.Hence,there is a need to estimate and mitigate the associated risks.In the past,various simulation methods have been developed and applied to induced seismicity analysis,but there is still a fundamental diference between simulation results and feld observations in terms of the spatial distribution of seismic events and its frequency.The present study aims to develop a method to simulate spatially distributed on-fault seismicity whilst reproducing a complex stress state in the fault zone.Hence,an equivalent continuum model is constructed,based on a discrete fracture network within a fault damage zone,by employing the crack tensor theory.A fault core is simulated at the center of the model as a discontinuous plane.Using the model,a heterogeneous stress state with stress anomalies in the fault zone is frst simulated by applying tractions on the model outer boundaries.Subsequently,the efective normal stress on the fault plane is decreased in a stepwise manner to induce slip.The simulation result is validated in terms of the b-value and other seismic source parameters,hence demonstrating that the model can reproduce spatially and temporally distributed on-fault seismicity.Further analysis on the parameters shows the variation of frequency-magnitude distribution before the occurrence of large seismic events.This variation is found to be consistent with feld observations,thus suggesting the potential use of this simulation method in evaluating the risk for seismic hazards in various engineering projects.

Numerical simulations for describing generation of excavation damaged zone: Important case study at Horonobe underground research laboratory

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。

First-arrival traveltime tomography for monitoring the excavation damaged zone in the Horonobe Underground Research Laboratory

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.

EXPERIMENTAL GH SINGULARITY FIELD AROUND GROWING CRACK-TIP

The displacement fields u_x,u_y at growing crack tip of LYI2-M specimens with double edge cracks are measured using moire method.The experimental singularity fields are compared with GH theoretical field [12-14].The size and shape of the experimental GH singularity fields are obtained.The error in both the experimental and theoretical evaluations is controlled within ±10%.The experiments show that there is (In A/r)^(+1)singularity dominant around a growing crack tip.The shape of this dominant region ranges from butterfly wing to oblate and circular.Inside GH-field,there is a 3-D deformed damage zone where no GH sin- gularity exists.

Experimental and numerical investigation into the non-explosive excavation of tunnels

The use of explosives is restricted on some important holidays,and the handling of unexploded charge is very dangerous.Therefore,an innovative non-explosive technology called instantaneous expansion(IE)was developed for tunneling.IE,whose components are derived from solid wastes such as coal gangue and straw conduces to realizing the reuse of waste.Moreover,its cost is lower than explosives.Blind guns of IE are easy to treat with water.The IE tunneling method is classified into two categories,i.e.IE with a single fracture(IESF)and IE with multiple fractures(IEMF),which are used to form the tunnel crosssection directionally cross-section and to fragment the rocks inside the cross-section,respectively.In this study,the principle of IE tunneling was elaborated first.Then,tunneling experiments and numerical simulations were performed on IE,conventional blasting(CB)and shaped charge blasting(SCB)in comparison.The experimental and numerical results show that IE achieved the best performance of directional rock breaking and corresponded to the most minor excavation-induced damage zone of the surrounding rock.Besides,the tunnel cross-section created by IE was flat and smooth.Comparing IE with CB and SCB,the over/under-excavation area decreased by 64%and 17%,and the excavation-induced damage zone fell by 26%and 11%,respectively.The range of the loose circle is reduced,which is conducive to improving the long-term stability of the roadway.The research provides a safe and economical tunneling method with excellent application prospects.

Laboratory tests and numerical simulations of brittle marble and squeezing schist at Jinping II hydropower station,China

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.

TIMODAZ:A successful international cooperation project to investigate the thermal impact on the EDZ around a radioactive waste disposal in clay host rocks

Disposal of spent nuclear fuel and long lived radioactive waste in deep clay geological formations is one of the promising options worldwide. In this concept of the geological disposal system, the host clay formation is considered as a principal barrier on which the fulfillment of key safety functions rests. Between 2006 and 2010, the European Commission project TIMODAZ, which gathered 15 partners from 8 countries, has investigated the coupled thermo-hydro-mechanical （THM） effects on clay formations for geological disposal of radioactive waste, and specific attention was paid to investigating the thermal effect on the evolution of the damaged zone （DZ）. Three types of potential host clay formations were investigated： the Boom Clay （Belgium）, the Opalinus Clay （Switzerland） and the Callovo-Oxfordian argillite （France）. Intensive experimental （laboratory and in situ in underground research laboratories） and numerical studies have been performed. Multi-scale approach was used in the course of the project. High degree of similarities between the failure modes, sealing process, stress paths, deformation, etc., observed in laboratories and in situ has been obtained, which increased the confidence in the applicability of laboratory test results and up-scaling perspective. The results of the laboratory and in situ tests obtained allowed the parameters for numerical models at various scales to be derived and provided the basis for the simplified performance assessment models that are used to assess the long-term safety of a repository. The good cooperation between the numerical modeler and experimenters has allowed an in-depth analysis of the experimental results and thus better understanding the underlying processes, and consequently increased the capabilities to model the THM effects in claystones. This paper presents the main achievements obtained by TIMODAZ project and shows how a European scientific community investigates a problem of concern in a collaborative way and how the obtained main results are applied to the performance assessment of a geological repository.

Fracture development around deep underground excavations: Insights from FDEM modelling

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.