Dissipation of Impact Stress Waves within the Artificial Blasting Damage Zone in the Surrounding Rocks of Deep Roadway

Artificial explosions are commonly used to prevent rockburst in deep roadways. However, the dissipation of the impact stress wave within the artificial blasting damage zone (ABDZ) of the rocks surrounding a deep roadway has not yet been clarified. The surrounding rocks were divided into the elastic zone, blasting damage zone, plastic zone, and anchorage zone in this research. Meanwhile, the ABDZ was divided into the pulverizing area, fractured area, and cracked area from the inside out. Besides, the model of the normal incidence of the impact stress waves in the ABDZ was established; the attenuation coefficient of the amplitude of the impact stress waves was obtained after it passed through the intact rock mass, and ABDZ, to the anchorage zone. In addition, a numerical simulation was used to study the dynamic response of the vertical stress and impact-induced vibration energy in the surrounding rocks. By doing so, the dissipation of the impact stress waves within the ABDZ of the surrounding rocks was revealed. As demonstrated in the field application, the establishment of the ABDZ in the surrounding rocks reduced the effect of the impactinduced vibration energy on the anchorage support system of the roadway.


Introduction
South Africa, Poland, Russia, Canada, Australia, and others have gradually begun to mine their deeper coal resources.However, while exploiting deep coal resources, rock burst occurs more often and with greater intensity in mine roadways.According to relevant statistics, approximately 80% of rock burst occurs in the roadways of deep coal seams annually [1].Rock burst is usually a sudden outburst, giving rise to vibration and damage in coal and rock masses.Meanwhile, its energy is expected to eject coal and rock masses so as to block the roadway, leading to collapse and damage to supports and equipment in the stope, as well as casualties.The control of rock burst in deep roadways has been an important topic in the past few years [2,3].
When the excavation in a deep roadway becomes stable, the fractured zone, plastic zone, and elastic zone are expected to be generated in the surrounding rocks: a large elastic strain energy usually accumulates at the junction of the elastic and plastic zones [4,5].Disturbed by the mining of the working face, the elastic strain energy is expected to be spread to the roadway surface in the form of stress waves (the impactinduced vibration energy).When the support system of the roadway (the supporting layer constituted by the anchorage zones of the bolt and anchor cable) cannot resist the impactinduced vibration energy generated by the stress waves, a rock burst occurs [6].
When a deep roadway is excavated, the radii of the fractured zone and the plastic zone in the surrounding rocks are affected by many factors, such as the ground stress, the radius of the roadway, and the support resistance of the support.As a matter of fact, the thicknesses of the fractured zone and the plastic zone formed naturally after the excavation usually fail to meet anti-impact requirements.Therefore, pressure relief technology using large diameter drill holes [7], artificial explosions [8], and coal seam water infusion [9] can be used to smash the surrounding rocks.The thickness of the fractured zone and the size of the fractured coal and rock masses are, therefore, regulated artificially, forming an anti-impact system with favourable energy absorption and cushioning performance.The system can absorb the energy released by stress waves and, therefore, provide ideal anti-impact effects.As for artificial explosions used to prevent rock burst, firstly, the drilling rig is usually employed to drill boreholes in the surrounding rocks.Then, explosives are placed in the borehole (usually at the junction of the elastic and plastic zones) before the borehole is sealed.While the artificial blasting damage zone (ABDZ) was established in the surrounding rocks, concentrated elastic strain energy can be transferred to greater depth.Besides, when the elastic strain energy is disturbed and spread outward in the form of stress waves, the fractured coal and rock masses in the ABDZ can dissipate some of the impact stress waves energy [10], as shown in Figure 1.In addition, the support system constituted by flexible energy absorption devices, such as the bolt and anchor cable outside the ABDZ, can resist the rest of energy, thus preventing the rock burst successfully [11].
However, the dissipation of impact stress waves in the fractured coal and rock masses inside the ABDZ is unclear.As a result, the artificial explosions adopted by many mines have difficulty in preventing rock burst.
Based on the fractured structures of the coal and rock masses in the ABDZ, a simplified model of the normal incidence of the impact stress waves was established.In this way, the attenuation coefficient of the amplitude of the impact stress waves was obtained after the waves passed through the ABDZ.Meanwhile, a numerical simulation was used to reveal the dissipation of the impact stress waves in the ABDZ.Thereafter, the model was verified in the field.

Simplified Model of the Normal Incidence of Impact Stress Waves
After the detonation of the spherical explosive cartridge in the surrounding rocks, the pulverizing zone, fractured zone, and crack zone were, respectively, generated from the inside out in the rock mass near the explosive cartridge [12].This zonal structure caused by the explosion of a single explosive  cartridge is called the "spherical rupture element."Moreover, the spherical rupture elements are supposed to constitute the ABDZ (Figure 2).As for a single spherical rupture element, when the elastic strain energy within the surrounding rocks is disturbed by mining, it is supposed to be spread in the coal mass in the form of impact stress waves.In this process, the impact stress wave is expected to pass through the crack zone, the fractured zone, and the pulverizing zone of the spherical rupture element in sequence and then pass through these zones again.At last, it enters into anchorage support zone of the roadway.Owing to the rupture element being spherical and assuming that the impact stress waves vertically penetrate the subareas of each rupture element in the ABDZ during their propagation, the waves finally act on the anchorage zone.

The Attenuation Coefficients of the Impact
Stress Waves within the Intact Rock Mass and the ABDZ in the Surrounding Rocks

The Propagation and Attenuation Coefficient of the Impact
Stress Waves within Intact Rock Mass.After being disturbed by mining, the elastic strain energy in the surrounding rocks of the deep roadway is expected to be propagated in the rock masses in the form of stress waves.Assuming that the intact rock mass within the surrounding rocks follows Kelvin-Voigt viscoelastic body behaviour, the amplitude of the stress waves at each point in the rock mass is expected to decrease with increasing propagation distance and time.
The two attenuation effects are called spatial and temporal attenuation, respectively [13].
The governing equation of the stress waves in a Kelvin-Voigt viscoelastic body is expressed as where   ,   , and  are the elastic modulus, viscosity coefficient, and time, respectively,  is the displacement of the stress waves along the propagation direction , and  0 is the density of the intact rock mass.Suppose that the distance between the vibration source and the particle is , and according to the harmonic wave equation, the amplitude attenuation of the wave at a certain point of the rock can be expressed as follows: where  and  are the amplitude and angular frequency and  is the responsive wave number of the stress waves in the coal and rock mass.By combining (1) and ( 2), the attenuation coefficients (  and   ) of the amplitude of the stress waves in a Kelvin-Voigt viscoelastic body, varying with time and space, are given by Therefore, the amplitude  of the stress waves with propagation distance () and time () in the intact rock mass can be expressed as where  0 represents the initial amplitude of the impact stress waves.

The Attenuation Coefficients of the Impact Stress
Waves in the ABDZ.Since the impact stress waves are projected vertically into the ABDZ, a microcuboidal element is separated from the spherical rupture element in the ABDZ.When the stress waves are propagated to the ABDZ, they are expected to pass through the crack zone, fractured zone, and pulverizing zone of the spherical rupture element in sequence and then pass through these zones again in the opposite order.At last, it acted on the anchorage zone of the roadway.Meanwhile, the stress waves are expected to generate complicated transmission and reflection phenomena between the subareas of the blasting damage zone [14], which will affect the propagation and attenuation of the wave stress in the rock mass (Figure 3).When the stress waves vertically enter the intact rock mass and the interface of the blasting damage zone  1 , part of the stress waves (with amplitude  1 ) will be transmitted through the partitioning interface  1 .Then, these waves are expected to generate the catadioptric phenomenon on interfaces  1 and  2 (the interface between the crack zone and the fractured zone in ABDZ), while some parts of the stress waves are expected to be transmitted through the interface  2 (the amplitudes of the reflected waves and transmitted waves are   2 and  2 , resp.) and then continue to be spread to the interface  3 .Afterwards, the stress waves continue to be propagated in additional zone of ABDZ according to the aforementioned order.Finally, the amplitude of the transmitted waves at the interface  6 is   .
Assume that V 0 , V 1 , . . ., V 4 are the velocities of the stress waves within the intact rock mass, crack zone, fractured zone, and pulverizing zone, respectively, and  0 ,  1 , . . .,  4 are the densities of the intact rock mass, crack zone, fractured zone, and pulverizing zone, separately.When the stress waves are propagated from the intact rock mass to the anchorage zone, its transmission coefficients on interfaces  1 ,  2 ,  3 ,  4 ,  5 , and  6 are  1 ,  2 ,  3 ,  4 ,  5 , and  6 , respectively.The widths of the intact rock mass, crack zone, fractured zone, and pulverizing zone are Δ 0 , Δ 1 , Δ 2 , and Δ 3 , separately.The initial wavelength of the stress waves is .
In accordance with other studies [15], the reflection frequency of the stress waves within the jointed rock mass is determined by the ratio of the wavelength to the width of the joint layer /Δ  ( = 1, 2, 3), as shown in When the stress waves are transmitted through the joint layer  1 for the first time, the attenuation rate of the amplitude is (Δ  ) 1 .After  reflections (assuming that  is an even number) between interfaces  1 and  2 , the amplitude of the stress waves is According to the literature [15], the amplitude of the stress waves in mine rocks can be expressed as where   , , and  are the spatial attenuation index, the propagation distance of the stress waves, and a constant value, respectively.According to (7), we obtain Let (Δ  )  =   ( = 1, 2, 3, 4, 5); thereinto, Δ 4 = Δ 2 and Δ 5 = Δ 1 .
After passing through the crack zone, the total amplitude  1 of the attenuated stress waves is Similarly, after passing through the rock masses in the crack zone, fractured zone, and pulverizing zone and entering the anchorage zone by taking the reverse journey, the amplitude of the stress waves is given by where  1 ,

Crack zone
Fractured zone Pulverizing zone

Fractured zone
Crack zone

Geological and Mining
Conditions in the Field.The 2304 working face in Tangkou coal mine (Zibo Mining Group Co., Ltd., Shandong Province, China) is located in the third coal seam with a burial depth of 1,000 m.In addition, the thickness of the coal seam ranges from 2.5 m to 3.5 m with an average thickness of 3.0 m, the dip angle of the coal seam ranges from 1 ∘ to 8 ∘ with an average dip of 6 ∘ , and the Protodyakonov coefficient of the coal is 0.7.The 2304 working face is located at the north flank of the western return airway and its south is the 2303 working face which is mining currently.There is no excavation at its north area.At present, the 2304 working face has not yet been mined, and its ventilation roadway is under excavation (Figure 4).Fully mechanised mining of the inclined retreating longwall is used: the strike and inclination lengths being 200 m and 1,700 m, respectively.The 2304 ventilation roadway is excavated in a rectangular section, which is 4 m long, 3 m wide, and supported permanently by cable anchors (Figure 5).Besides, thread steel resin bolts containing no longitudinal reinforcement with diameters and lengths of 20 mm and 2,400 mm, respectively, are used in the roof.The bolts, each of which is fixed using two CK2850 resin anchor agents, are arranged with row and line spacings on an 800 × 800 mm grid.The anchor cables, which are 17.8 mm in diameter and 9,000 mm long, are set with the intervals in the rows and line spacings being 1,200 mm and 1,600 mm, respectively.Moreover, three CK2350 resin anchor agents are placed into each hole for end anchorage.Thread steel resin bolts without longitudinal reinforcement are used on the walls of the roadway.The bolts have a diameter and length of 20 mm and 2,200 mm, respectively, and are fixed at 800 mm intervals (both row and line spacings).In addition, the 100 × 100 mm reinforcing mesh is made of reinforcing bars with a diameter of 6.5 mm.The length, width, and thickness of the steel strip are 3,600 mm, 80 mm, and 5 mm, respectively.
As shown in Figure 4, the 2304 ventilation roadway has been excavated for a length of 1,300 m.During roadway excavation, dynamic accidents such as the projection of coal and rock masses and the collapse of the roof and the wall of roadway took place frequently.
The TDS-6 microseismic acquisition system was used to monitor and collect the signal from the impact stress waves in the heading face of the 2304 ventilation roadway every 0.2 ms.After these collected signals were processed, the velocity-time curves of the signals of the impact stress waves were obtained (Figure 6(a)).To simplify the study, the harmonic stress wave with the maximum vibration velocity and shortest period was selected from the curves obtained above (Figure 6(b)).Then, the selected harmonic stress curve was applied as the impact stress wave acting on the surrounding rocks of the roadway for two periods of the numerical simulation.Rock cores were obtained in the field by drilling.Meanwhile, laboratory testing to find the rock mechanical parameters was carried out according to the Standard for Tests Method of Engineering Rock Masses (GB/T 50266-2013).The lithology and mechanical parameters of the roof and floor are summarised in Table 1.

Establishment of the Model and Simulation Schemes.
Based on the geological and mining conditions of the 2304 ventilation roadway, the dissipation of the impact-induced vibration energy in the ABDZ was studied.Meanwhile, FLAC 3D was used to establish a numerical calculation model, whose length, width, and height were 52 m, 10 m, and 34 m, respectively.The width, height, and length of the roadway were 4 m, 3 m, and 10 m, respectively.According to in situ stress measurements on the 2304 ventilation roadway, the rock stress was mainly geostatic at 25 MPa in the vertical direction and 32.5 MPa in the horizontal direction (the lateral pressure coefficient being 1.25).Therefore, horizontal restraint was applied to the and -directions of the model, while the bottom of -direction was fixed [16].Besides, the vertical stress ℎ of 25 MPa generated by the self-weight of the overlying strata was imposed on the top boundary.When the calculation reaches a balance after the excavation process, the impact stress waves were imposed.During the simulation, the static boundary was set and an observation line established   along the vertical direction at the middle of the roadway roof to monitor the change of elastic strain energy.The simulation schemes are as follows.
Scheme 1.After the static calculation reaches equilibrium, impact stress waves are imposed on the boundary between the elastic region and the plastic region at the midline above the roadway roof for subsequent dynamic calculation.This scheme aims to simulate the propagation of impact energy in the non-blast-damage zone.Scheme 2. After reaching equilibrium, the ABDZ, to a range of 2 m, is set on the boundary between the elastic region and the plastic region of the roadway roof according to published data [17] (Figure 7).When the static calculation reaches equilibrium for a second time, the impact stress waves are applied on the boundary between the elastic region and the ABDZ at the midline above the roadway roof for subsequent dynamic calculation.
The changes in vertical stress and impact-induced vibration energy in the surrounding rocks were observed at the midline above the roadway roof, so as to analyse the influence of the ABDZ on the impact stress waves and determine the dissipation of the impact-induced vibration energy therein.on the roof of the roadway.Then, according to Scheme 1, the impact stress waves are applied to the zone (the interface between the elastic and the plastic regions) at a distance of 9 m from the surface of the roadway roof.According to Scheme 2, the blasting damage zone lying within a range of 2 m was set in the boundary between the elastic and the plastic regions at the midline above the roadway roof: the impact stress waves were then imposed on the boundary between the elastic region and the blasting damage zone (Figure 7).The variations in vertical stress and impact-induced vibration energy of the unit body in the surrounding rocks illustrated in Schemes 1 and 2 are shown in Figures 8-12.
(1) The Change of Vertical Stress in the Surrounding Rocks under the Action of the Impact Stress Waves.Under the influence of the impact stress waves for the 3 ms, the change in vertical stress in the surrounding rocks is shown in Figures 8-11.
The vertical stress in the roadway roof is large at  = 2 ms.It can be explained by the fact that the compression and tension waves are contained within one period of the impact stress waves.The vertical stress in the surrounding rocks is large under the action of the compression wave, while it is expected to decrease when affected by the tension wave.At  = 2 ms, according to Scheme 1, the vertical stress, at approximately 9 m above the roadway roof, reaches 70 MPa, with a stress concentration factor of 2.8, while the vertical stress decreases gradually from a point 9 m above the roadway roof to 2 m above it, with the decreasing amplitude per unit length being 6.14 MPa.The decreasing amplitude is approximately 10 MPa⋅m −1 from 2 m above the roadway roof to the roadway surface (Figure 9).In contrast, the maximum vertical stress at the boundary between the ABDZ and the elastic zone in the roadway surrounding rocks (11 m above the roof at  = 2 ms) is about 70 MPa in Scheme 2. However, Figure 11 shows the decreasing amplitude of the vertical stress in the ABDZ of the surrounding rocks to be approximately 15 MPa⋅m −1 .Furthermore, the decreasing amplitudes of the vertical stress are 2.85 MPa⋅m −1 from 9 m to 2 m above the roof and 10 MPa⋅m −1 approximately from 2 m above the roof to the roadway surface, separately.
(2) The Dynamic Response of the Impact-Induced Vibration Energy in the Surrounding Rocks under the Action of the Impact Stress Waves. Figure 12 shows the dynamic response of the impact-induced vibration energy at the medial axis of the roadway roof under the effect of the impact stress waves.For Scheme 1, the impact-induced vibration energy is large (248.6 kJ at 9 m above the roof), while it decreases gradually from 9 m to 2 m above the roof at 20 kJ⋅m −1 .The impactinduced vibration energy close to the roadway surface is approximately 30 kJ, where a rock burst is likely.By contrast, in Scheme 2, the impact-induced vibration energy shows a maximum decreasing amplitude in the ABDZ, being approximately 40 kJ⋅m −1 .The decreasing amplitude of the impactinduced vibration from 9 m to 2 m above the roof is basically the same as that in Scheme 1.The impact-induced vibration energy near the roadway is approximately 2 kJ, which is far smaller than the threshold value (25 kJ) for the occurrence of a rock burst [18].
According to the comparative analysis of the dynamic responses of the vertical stress and the impact-induced vibration energy in the roadway surrounding rocks with, and   without, the ABDZ, under the action of impact stress waves, the following conclusions can be obtained: under the action of the impact stress waves, the vertical stress and the impactinduced vibration energy within the roadway surrounding rocks are increased significantly (Figures 8-12).When the impact stress waves are propagated to the roadway surface through the ABDZ of the surrounding rocks, most of the impact-induced energy is dissipated in the form of heat energy, surface energy, and plastic energy [19].This is because the ABDZ is a fractured structure consisting of multiple joints and fractures, and the impact-induced vibration energy leads to the dislocation, extrusion, and repeated damage of the fractured rock masses [20,21].As a result, the impact-induced vibration energy is weakened substantially  to generate a reduced stress zone in the ABDZ.Meanwhile, the impact-induced vibration energy within the anchorage zone of the surrounding rocks is reduced accordingly (the area situated within 9 m to 11 m above the roadway roof in Figure 12).Therefore, the proposed Scheme 2 can effectively reduce the impact-induced vibration energy.4).A high-power pneumatic drill was used to drill five boreholes in the two walls and roof of the roadway (Figure 13): the length and diameter of the boreholes were 20 m and 42 mm, respectively.The method of drill bit analysis was used to determine the energy concentration zone in the surrounding rocks.Since the surrounding rocks from 1 m to 3 m depth into the roadway sides are already damaged, the cuttings obtained from depths of 1 to 3 m in the boreholes were abandoned in the field.Samples were therefore cut from no less than 4 m deep.The test results are shown in Figure 14.

Application Example
As shown in Figure 14, at borehole depths of 9 to 10 m, the cuttings index of the boreholes reaches a maximum (4.1 kg/m).This indicates that strong elastic strain energy is contained herein; however, if the elastic strain energy is disturbed by the mining of the excavated working face, it is likely to generate intensive impact stress waves and even give rise to a rock burst.This is because the anchorage support system and the roadway surrounding rocks themselves are not enough to resist the impact-induced vibration energy.
To reduce the possibility of a rock burst in the 2304 ventilation roadway, a drilling rig was used to drill a borehole (length, 9 m and radius, 0.75 m) in the surrounding rocks.Soon afterwards, a spherical explosive cartridge was used to conduct coupling blasting, on the one hand, to transfer the concentration zone of the elastic strain energy to the deeper surrounding rocks.On the other hand, it aims to form an ABDZ, in which the fractured coal and rock masses can dissipate some of the impact-induced vibration energy (Figure 15).Immediately after the explosion, lengthened anchor cables, 15 m long, were installed in the two walls and roof of the roadway (Figure 16).

Impact Resistance of the 2304 Ventilation Roadway.
To verify the dissipation effect of the ABDZ on the impactinduced vibration energy, the method of drill bit analysis was performed on the aforementioned observation stations (1, 2, and 3 in Figure 4) to determine the cuttings index of the borehole in roadway surrounding rocks (Figure 17).By establishing the ABDZ in the 2304 ventilation roadway, the cuttings  index of the boreholes in the surrounding rocks is reduced, with the maximum decrease being approximately 33%.That is, the elastic strain energy in the surrounding rocks of the deep roadway is reduced.Meanwhile, by increasing the number of anchorage cables, the resistance of the anchorage support system in the roadway, to prevent impact-induced vibration energy, is strengthened, which gives rise to a significantly improved impact resistance of the deep roadway.As demonstrated in the field monitoring, no rock burst was found in the 2304 ventilation roadway from November, 2012, to November, 2013, and there was little deformation of the surrounding rocks.Therefore, the 2304 ventilation roadway satisfies the mine ventilation, human entry, and transportation requirements.

Conclusions
(1) While establishing the artificial blasting damage zone in the surrounding rocks of the deep roadway, the pulverizing zone, fractured zone, and crack zone are expected to be generated in the blasting damage zone to form the spherical rupture element.When the elastic strain energy within the surrounding rocks of the deep roadway is disturbed, it is expected to be propagated through the rock mass in the form of impact stress waves.Thereafter, the stress waves are supposed to pass through the crack zone, fractured zone, and pulverizing zone of the spherical rupture element in sequence and then pass through these zones again in the opposite order, last to enter the anchorage zone of the roadway.
(2) A numerical simulation was used to study the changes in vertical stress and elastic strain energy in the surrounding rocks under the action of the impact stress waves.As demonstrated by the simulation results, while establishing the blasting damage zone in the rocks surrounding this deep roadway, the decreasing amplitudes of the vertical stress and the impactinduced vibration energy in the blasting damage zone are 15 MPa⋅m −1 and 40 KJ⋅m −1 , respectively.
(3) As demonstrated by the case study, the establishment of the ABDZ in the rocks surrounding this deep roadway and the enhancement of the anchorage support system have strengthened the ability to prevent impact-induced vibration energy.Meanwhile, the effect of the impact-induced vibration energy on the anchorage support system of the roadway is reduced, which gives rise to the significantly increased impact resistance of this deep roadway.

Figure 1 :
Figure 1: Stress distribution in the surrounding rocks with the ABDZ of the deep roadway.

Figure 2 :
Figure 2: Simplified model of the normal incidence of the impact stress waves.

Figure 3 :
Figure 3: The incidence and reflection of a stress wave on the interface of the joint layer.

Figure 4 :Figure 5 :
Figure 4: Plan layout of the working face.

Figure 9 :
Figure 9: The change in vertical stress in the surrounding rocks: Scheme 1.
Action time  = 3 ms

Figure 10 :
Figure 10: Dynamic response of the vertical stress in the surrounding rocks: Scheme 2.

Figure 11 :
Figure 11: Change in vertical stress in the surrounding rocks: Scheme 2.

Figure 12 :
Figure 12: Dynamic response of the impact-induced vibration energy of the strata in the roadway roof.

Table 1 :
The lithology and mechanical parameters of the roadway roof and floor.