A 3D microseismic data-driven damage model for jointed rock mass under hydro-mechanical coupling conditions and its application

Rock mass is a fractured porous medium usually subjected to complex geostress and fluid pressure simultaneously.Moreover,the properties of rock mass change in time and space due to mining-induced fractures.Therefore,it is always challenging to accurately measure rock mass properties.In this study,a three-dimensional(3D)microseismic(MS)data-driven damage model for jointed rock mass under hydro-mechanical coupling conditions is proposed.It is a 3D finite element model that takes seepage,damage and stress field effects into account jointly.Multiple factors(i.e.joints,water and microseismicity)are used to optimize the rock mass mechanical parameters at different scales.The model is applied in Shirengou iron mine to study the damage evolution of rock mass and assess the crown pillar stability during the transition from open-pit to underground mining.It is found that the damage pattern is mostly controlled by the structure,water and rock mass parameters.The damage pattern is evidently different from the two-dimensional result and is more consistent with the field observations.This difference is caused by the MS-derived damage acting on the rock mass.MS data are responsible for gradually correcting the damage zone,changing the direction in which it expands,and promoting it to evolve close to reality.For the crown pillar,the proposed model yields a more trustworthy safety factor.In order to guarantee the stability of the pillar,it is suggested to take waterproof and reinforcement measures in areas with a high degree of damage.

Vibration reduction technology and the mechanisms of surrounding rock damage from blasting in neighborhood tunnels with small clearance

Interlaid rock is an important component in the construction of neighborhood tunnels that supports and reinforces the area between two tunnels.However,the blasting load during excavation can sometimes damage the interlaid rock and threaten the stability of a tunnel’s structure.This paper presents a case study of the small clearance section of the Liantang highway tunnel project in Shenzhen,China,where the minimum distance between the two tunnels involved is only 0.5 m.To reduce the damage to the interlaid rock caused by blasting loads,we proposed a four-part excavation method with a vibrationcushioning rock layer in the following tunnel of neighborhood tunnels.Numerical simulation was used to model the damage prevention mechanism of the vibration-cushioning rock layer and to better understand the propagation of cracks in the interlaid rock.Furthermore,based on the simulation results,combined microseismic controlled-blasting technology was implemented,using innovative blasting patterns combined with different charge structures and blasting equipment designed according to the varying thickness of the interlaid rock.Finally,this implementation succeeded in protecting interlaid rock during blasting operations.

Preliminary results on a near-real-time rock slope damage monitoring system based on relative velocity changes following the September 5,2022 M_(S) 6.8 Luding,China earthquake

Relative seismic velocity change(dv/v)is important for monitoring changes in subsurface material properties and evaluating earthquake-induced rock slope damage in a geological disaster-prone region.In this paper,we present a rapid damage assessment on three slow-moving rock slopes by measuring dv/v decrease caused by the 2022 M_(S) 6.8 Luding earthquake in Southwest China.By applying the stretching method to the cross-correlated seismic wavefields between sensors installed on each slope,we obtain earthquake-induced dv/v decreases of~2.1%,~0.5%,and~0.2%on three slopes at distances ranging from~86 to~370 km to the epicenter,respectively.Moreover,based on seismic data recorded by 16 sensors deployed on the rock slope at a distance of~370 km away from the epicenter,a localized dv/v decease region was observed at the crest of the slope by calculating the spatial dv/v images before and after the earthquake.We also derive an empirical in situ stress sensitivity of -7.29×10^(-8)/Pa by relating the dv/v change to the measured peak dynamic stresses.Our results indicate that a rapid dv/v assessment not only can help facilitate on-site emergency response to earthquakeinduced secondary geological disasters but also can provide a better understanding of the subsurface geological risks under diverse seismic loadings.

Effects of earthquake on damage evolution and failure mechanism of key rock pillars in underground engineering

Rock pillar is the key supporting component in underground engineering.During an earthquake,the key rock pillar must bear both the seismic load and the load transferred from other damaged pillars.This paper attempts to reveal the influence of the mainshock on damage evolution and failure characteristic of the key rock pillar during aftershocks by cyclic loading test of marble.Four levels of pre-damage stress(i.e.,10,30,50 and 70 MPa)in the first cycle were used to simulate the mainshock damage,and then cyclic stress with the same amplitude(namely 10 MPa)was conducted in the subsequent cycles to simulate the aftershock until rock failure.The results indicate that the presence of pre-damage has an obvious weakening effect on the bearing capacity and deformation resistance of rock materials during the aftershock process.Besides,the increase of pre-damage significantly changes the final failure pattern of the key rock pillar,and leads to an increase in the proportion of small-scale rock fragments.This study may contribute to understanding the seismic capacity of the unreinforced rock pillar during mainshock-aftershock seismic sequences and to optimizing the design of the key rock pillar in underground engineering.

Thermal damage constitutive model for rock considering damage threshold and residual strength

With the gradual depletion of mineral resources in the shallow part of the earth,resource exploitation continues to move deeper into the earth,it becomes a hot topic to simulate the whole process of rock strain softening,deformation and failure in deep environment,especially under high temperature and high pressure.On the basis of Lemaitre’s strain-equivalent principle,combined with statistics and damage theory,a statistical constitutive model of rock thermal damage under triaxial compression condition is established.At the same time,taking into account the existing damage model is difficult to reflect residual strength after rock failure,the residual strength is considered in this paper by introducing correction factor of damage variable,the model rationality is also verified by experiments.Analysis of results indicates that the damage evolution curve reflects the whole process of rock micro-cracks enclosure,initiation,expansion,penetration,and the formation of macro-cracks under coupled effect of temperature and confining pressure.Rock thermal damage shows logistic growth function with the increase of temperature.Under the same strain condition,rock total damage decreases with the rise of confining pressure.By studying the electron microscope images(SEM)of rock fracture,it is inferred that 35.40 MPa is the critical confining pressure of brittle to plastic transition for this granite.The model parameter F reflects the average strength of rock,and M reflects the morphological characteristics of rock stress–strain curves.The physical meanings of model parameters are clear and the model is suitable for complex stress states,which provides valuable references for the study of rock deformation and stability in deep engineering.

Rock Damage Structure of the South Longmen-Shan Fault in the 2008 M8 Wenchuan Earthquake Viewed with Fault-Zone Trapped Waves and Scientific Drilling

This article is to review results from scientific drilling and fault-zone trapped waves （FZTWs） at the south Longman-Shan fault （LSF） zone that ruptured in the 2008 May 12 M8 Wenchuan earthquake in Sichuan,China.Immediately after the mainshock,two Wenchuan Fault Scientific Drilling （WFSD） boreholes were drilled at WFSD-1 and WFSD-2 sites approximately 400 m and 1 km west of the surface rupture along the Yinxiu-Beichuan fault （YBF）,the middle fault strand of the south LSF zone.Two boreholes met the principal slip of Wenchuan earthquake along the YBF at depths of 589-m and 1230-m,respectively.The slip is accompanied with a 100-200-m-wide zone consisting of fault gouge,breccia,cataclasite and fractures.Close to WFSD-1 site,the nearly-vertical slip of ～4.3-m with a 190-m wide zone of highly fractured rocks restricted to the hanging wall of the YBF was found at the ground surface after the Wenchuan earthquake.A dense linear seismic array was deployed across the surface rupture at this venue to record FZTWs generated by aftershocks.Observations and 3-D finite-difference simulations of FZTWs recorded at this cross-fault array and network stations close to the YBF show a distinct low-velocity zone composed by severely damaged rocks along the south LSF at seismogenic depths.The zone is several hundred meters wide along the principal slip,within which seismic velocities are reduced by ～30-55％ from wall-rock velocities and with the maximum velocity reduction in the ～200-m-wide rupture core zone at shallow depth.The FZTW-inferred geometry and physical properties of the south LSF rupture zone at shallow depth are in general consistent with the results from petrological and structural analyses of cores and well log at WFSD boreholes.We interpret this remarkable low-velocity zone as being a break-down zone during dynamic rupture in the 2008 M8 earthquake.We examined the FZTWS generated by similar earthquakes before and after the 2008 mainshock and observed that seismic velocities within fault core zone was reduced by ～10％ due to severe damage of fault rocks during the M8 mainshock.Scientific drilling and locations of aftershocks generating prominent FZTWs also indicate rupture bifurcation along the YBF and the Anxian-Guangxian fault （AGF）,two strands of the south LSF at shallow depth.A combination of seismic,petrologic and geologic study at the south LSF leads to further understand the relationship between the fault-zone structure and rupture dynamics,and the amplification of ground shaking strength along the low-velocity fault zone due to its waveguide effect.

Simulation experiment of AE-based localization damage and deformation characteristic on covering rock in mined-out area

The localization damage and deformation characteristics of multiple rock in mined-out area under thick soil layer were investigated by acoustic emission （AE）-based physical similar simulation experiment and damage test combined with the AE theory on rock damage. During different mining courses, the degree of wall rock suffered from mining disturbed process in mined-out area was different, so did the AE density of initial damage and fracture. Some characteristic parameters, including large events, total events, and energy rate, presented fluctuations with temporal and spatial change and take on a certain extent statistic relation. Subsequently, the percent of damage degree was defined and divided into 4 stages, and the localization damage and deformation characteristics of surrounding rock in mined-out area and coal pillars were analyzed quantitatively.

Responses of jointed rock masses subjected to impact loading

Impact-induced damage to jointed rock masses has important consequences in various mining and civil engineering applications. This paper reports a numerical investigation to address the responses of jointed rock masses subjected to impact loading. It also focuses on the static and dynamic properties of an intact rock derived from a series of laboratory tests on meta-sandstone samples from a quarry in Nova Scotia, Canada. A distinct element code(PFC2D) was used to generate a bonded particle model(BPM) to simulate both the static and dynamic properties of the intact rock. The calibrated BPM was then used to construct large-scale jointed rock mass samples by incorporating discrete joint networks of multiple joint intensities into the intact rock matrix represented by the BPM. Finally, the impact-induced damage inflicted by a rigid projectile particle on the jointed rock mass samples was determined through the use of the numerical model. The simulation results show that joints play an important role in the impactinduced rock mass damage where higher joint intensity results in more damage to the rock mass. This is mainly attributed to variations of stress wave propagation in jointed rock masses as compared to intact rock devoid of joints.

A Study of the Rock Breaking Mechanism during Swirling Water Jet Drilling

Based on an analysis of the factors affecting rock breaking and the coupling between rock and fluid during water jet drilling, the rock damage model and the damage-coupling model suitable for the whole rock breaking process under the water jet is established with continuous damage mechanics and micro-damage mechanics. The evolvement of rock damage during swirling water jet drilling is simulated on a nonlinear FEM and dynamic rock damage model, and a decoupled method is used to analyze the rock damage. The numerical results agree with the test results to a high degree, which shows the rock breaking ability of the swirling water jet is strong. This is because the jet particle velocity of the swirling water jet is three-dimensional, and its rock-breaking manner mainly has a slopping impact. Thus, the interference from returning fluid is less. All these aspects make it easy to draw and shear the rock surface. The rock breaking process is to break out an annular on the rock surface first, and then the annular develops quickly in both the radial and axial directions, the last part of the rock broken hole bottom is a protruding awl. The advantage of the swirling water jet breaking rock is the heavy breaking efficiency,large breaking area and less energy used to break rock per unite volume, so the swirling water jet can drill in a hole of a large diameter.

Numerical investigation of effect of eccentric decoupled charge structure on blasting-induced rock damage

Eccentric decoupling blasting is commonly used in underground excavation.Determination of perimeter hole parameters(such as the blasthole diameter,spacing,and burden)based on an eccentric charge structure is vital for achieving an excellent smooth blasting effect.In this paper,the Riedel-Hiermaier-Thoma(RHT)model was employed to study rock mass damage under smooth blasting.Firstly,the parameters of the RHT model were calibrated by using the existing SHPB experiment,which were then verified by the existing blasting experiment results.Secondly,the influence of different charge structures on the blasting effect was investigated using the RHT model.The simulation results indicated that eccentric charge blasting has an obvious pressure eccentricity effect.Finally,to improve the blasting effect,the smooth blasting parameters were optimized based on an eccentric charge structure.The overbreak and underbreak phenomena were effectively controlled,and a good blasting effect was achieved with the optimized blasting parameters.

Meso-damage evolution and mechanical characteristics of low-porosity sedimentary rocks under uniaxial compression

Low pore sedimentary rocks(from Guangxi, China) were subjected to uniaxial compression loading experiment under different initial stresses. The rock samples were investigated by nuclear magnetic resonance before and after the loading. The relationships between the mesoscopic rock damage and macroscopic mechanical parameters were established, and the initial damage stress of the low-porosity sedimentary rock was determined. The results showed that this type of rock has the initial stress of damage. When the initial loading stress is lower than the initial stress of damage, the T2 spectrum area of the rock sample gradually decreases, and the primary pores of the rock are further closed under the stress. The range of the initial stress of damage for this type of rock is 8-16 MPa. When the loading stress exceeds the initial stress of damage, the T2 spectrum area gradually increases, indicating that the porosity of the rock increases and microscopic damage of the rock appears. The rock damage degree is defined, and the nonlinear function between the rock damage degree and the initial loading stress is established.

Effects of loading waveforms on rock damage using particle simulation method

The particle simulation method is used to study the effects of loading waveforms （i.e. square, sinusoidal and triangle waveforms） on rock damage at mesoscopic scale. Then some influencing factors on rock damage at the mesoscopic scale, such as loading frequency, stress amplitude, mean stress, confining pressure and loading sequence, are also investigated with sinusoidal waveform in detail. The related numerical results have demonstrated that： 1） the loading waveform has a certain effect on rock failure processes. The square waveform has the most damage within these waveforms, while the triangle waveform has less damage than sinusoidal waveform. In each cycle, the number of microscopic cracks increases in the loading stage, while it keeps nearly constant in the unloading stage. 2） The loading frequency, stress amplitude, mean stress, confining pressure and loading sequence have considerable effects on rock damage subjected to cyclic loading. The higher the loading frequency, stress amplitude and mean stress, the greater the damage the rock accumulated; in contrast, the lower the confining pressure, the greater the damage the rock has accumulated. 3） There is a threshold value of mean stress and stress amplitude, below which no further damage accumulated after the first few cycle loadings. 4） The high-to-low loading sequence has more damage than the low-to-high loading sequence, suggesting that the rock damage is loading-path dependent.

Simulations of explosion-induced damage to underground rock chambers

A numerical approach is presented to study the explosion-induced pressure load on an underground rock chamber wall and its resultant damage to the rock chamber.Numerical simulations are carried out by using a modified version of the commercial software AUTODYN.Three different criteria,i.e.a peak particle velocity (PPV) criterion,an effective strain (ES) criterion,and a damage criterion,are employed to examine the explosion-induced damaged zones of the underground rock chamber.The results show that the charge chamber geometry,coupling condition and charge configuration affect significantly the dynamic pressure exerted on the rock chamber wall.Thus the chamber is damaged.An inaccurate approximation of pressure boundary ignoring the influences of these factors would result in an erroneous prediction of damaged area and damage intensity of the charge chamber.The PPV criterion yields the largest damaged zone while the ES criterion gives the smallest one.The presented numerical simulation method is superior in consideration of the chamber geometry,loading density,coupling condition and rock quality.The predicted damage intensity of rock mass can be categorized quantitatively by an isotropic damage scalar.Safe separation distance of adjacent chambers for a specific charge weight is also estimated.

Vibration Effect and Damage Evolution Characteristics of Tunnel Surrounding Rock Under Cyclic Blasting Loading

Model test studies based on the similarity theory were conducted to investigate vibration effect and damage evolution characteristics of tunnel surrounding rock under push-type cyclic blasting excavation.The model was constructed with a ratio of 1∶15.By simulating the tunnel excavation of push-type cyclic blasting,the influence of the blasting parameter change on vibration effect was explored.The damage degree of tunnel surrounding rock was evaluated by the change of the acoustic wave velocity at the same measuring point after blasting.The relationship between the damage evolution of surrounding rock and blasting times was established.The research results show that：（1）In the same geological environment,the number of delay initiation is larger,the main vibration frequency of blasting seismic wave is higher,and the attenuation of high frequency signal in the rock and soil is faster.The influence of number of delay initiation on blasting vibration effect cannot be ignored;（2）Under push-type cyclic blasting excavation,there were great differences in the decreasing rates of acoustic wave velocity of the measuring points which have the same distance to the blasting region at the same depth,and the blasting damage ranges of surrounding rock were typically anisotropic at both depth and breadth;（3）When blasting parameters were basically kept as the same,the growth trend of the cumulative acoustic wave velocity decreasing rate at the measuring point was nonlinear under different cycle blasting excavations;（4）There were nonlinear evolution characteristics between the blasting cumulative damage（D）of surrounding rock and blasting times（n）under push-type cyclic blasting loading,and different measuring points had corresponding blasting cumulative damage propagation models,respectively.The closer the measuring point was away from the explosion source,the faster the cumulative damage extension.Blasting cumulative damage effect of surrounding rock had typically nonlinear evolution properties and anisotropic characteristics.

Discussion on blasting vibration monitoring for rock damage control in rock slope excavation

Drill and blast is a commonly used method for rock slope excavation in hydropower engineering.During blasting excavation of rock slopes,far-field vibration monitoring on the first upper berm for statutory compliance is usually performed to control the blast-induced rock damage to the final slope face.In this study,for the rock slope excavation in the Jinping-I hydropower station,the field vibration monitoring and acoustic testing are presented to investigate the vibration characteristics on the first upper berm and the damage depth in the current bench.The relationship between the PPV on the first upper berm and the PPV damage threshold on the damage zone boundary is also studied through three-dimensional FEM simulations.The results show that on the first upper berm,the maximum vibration velocity component occurs in the vertical direction.While on the blasting damage zone boundary,the horizontal radial vibration velocity is the maximum component.For the Jinping-I slope with a bench height of 30 m,the radial PPV on the inner side of the first upper berm is 2.06%of the PPV threshold on the damage zone boundary.This ratio is increased as the bench height decreases.Therefore,the bench height of the rock slope is an important factor that cannot be ignored in determining the allowable vibration velocity for rock damage control.

An overview of the thermography-based experimental studies on roadway excavation in stratified rock masses at CUMTB

This paper presents an overview of experimental investigations conducted at China University of Mining and Technology Beijing(CUMTB) on roadway excavation using large-scale geomechanical model tests.The simulated sedimentary rocks are composed by alternating layers of sandstone, mudstone and coal seam inclined at varied angles with respect to the horizontal including 0°, 45°, 60°, and 90°. During the excavation, infrared thermography was employed to detect the thermal response of the surrounding rocks under excavation. The obtained raw thermograms were processed using denoising algorithm, data reduction procedure and Fourier analysis. The infrared temperature(IRT) characterizes the overall rock response; the processed thermal images represent the structural behavior, and the Fourier spectrum describes damage development in the frequency domain. Deeper understanding was achieved by the comparative analyses of excavation in differently inclined rock masses using the image features of IRTs, thermal images and Fourier spectra.

NUMERICAL SIMULATION STUDY ON ROCK BREAKING MECHANISM AND PROCESS UNDER HIGH PRESSURE WATER JET

The numerical simulation method to study rock breaking process and mechanism under high pressure water jet was developed with the continuous mechanics and the FEM theory. The rock damage model and the damage-coupling model suited to analyze the whole process of water jet breaking rock were established with continuum damage mechanics and micro damage mechanics. The numerical results show the dynamic response of rock under water jet and the evolvement of hydrodynamic characteristic of jet during rock breaking is close to reality, and indicates that the body of rock damage and breakage under the general continual jet occurs within several milliseconds, the main damage form is tensile damage caused by rock unload and jet impact, and the evolvement of rock damage shows a step-change trend. On the whole, the numerical results can agree with experimental conclusions, which manifest that the analytical method is feasible and can be applied to guide the research and application of jet breaking rock theory.

Damage of reservoir rock induced by CO_(2) injection

The study of reservoir rock damage induced by gas injection is of great significance to the design of reservoir stimulation and the improvement of oil and gas recovery. Based on an example horizontal well in the Hudson Oilfield of the Tarim Basin and considering the multi-physics coupling effects among highpressure fluid, rock deformation, and damage propagation during CO_(2) injection, a three-dimensional finite element model for CO_(2) injection in deep reservoir considering seepage-stress-damage coupling was developed. The evolution of reservoir rock damage under different CO_(2) injection conditions was systematically investigated. The results show that tensile damage and shear damage are concentrated in the vertical direction and the horizontal maximum compressive principal stress direction, respectively,and the tensile damage is the main damage mode. At higher CO_(2) injection rate and pressure, the damaged areas near the wellbore are mainly distributed in the direction of the maximum compressive principal stress, and the development of the damaged area near the wellbore will be inhibited by the formation and evolution of far-field damage. CO_(2) injection aggravates the extension of tensile damage,but inhibits the initiation of shear damage, and eventually leads to the gradual transition from shear damage to tensile damage. Under the same injection conditions, CO_(2) injection has superior performance in creating rock damage compared with the injection of nitrogen and water. The results in this study provide guidance for enhanced oil recovery in deep oil and gas reservoirs with CO_(2) injection.

A CREEP DAMAGE EQUATION FOR ROCKS

In this paper,based on many experimental results,a creep damage equation for rocks with a creep modulus that can describe the history of damage has been deduced.According to the equation,a creep damage model for rock materials is constructed,and it is shown that the extent of their damage,or damage degree,during creep can be determined by means of simple experiments.

Analysis on damage and rheological characteristics of deep surrounding rock of shaft engineering

According to the characteristics of deep engineering surrounding rock main shaft of No.3 mining district in Jinchuan, electron microscope scanning and rock mechanics test were adopted to analyze the damage features of rock. The software of FLAG3D and Burgers body （Kelvin-Maxwell model） were used to research on rheological theory, and rheological model was modified. The results indicate that the damage of rock mass is very serious, and the rheological characteristics also outstanding; rheological behavior of deep surrounding rocks of the shaft can be taken as superposition of transient and stable rheology; and there exist the most dangerous zone on 100 m higher than 1 063 m level, so it is necessity that works of monitor and corresponding reinforcement should strengthen.