Rockburst Prediction on the Superimposed Effect of Excavation Accumulation Energy and Blasting Vibration Energy in Deep Roadway

In deep mining, much elastic energy is stored in rock mass due to the high geostress. Rockburst will be induced by accumulated energy during excavation.Meanwhile, because of blasting vibration energy in the host rock, there will be an obvious superimposed effect on the probability and intensity of rockburst. To explore the most reasonable and effective method for understanding rockburst problem under blasting, a deep roadway of Sanshandao gold mine was studied. On the basis of in situ geostress data, the accumulated energy of three-centered arch roadway after quasi-static excavation was derived. *en, a series of in situ blasting vibration were monitored, and the blasting vibration energy was calculated by employing the equivalent theory of elastic vibration boundary. Finally, the tendentiousness of rockburst was evaluated qualitatively with the superimposed energy. *e results indicated that the disaster-driven energy was increased by 45.1% and 28.2% on different places of roadway. Also, the probability and intensity of rockburst would be raised.


Introduction
Nowadays, resource exploration has entered deep mining. Much elastic energy is stored in rock mass, and deeper buried depth leads to more energy. Many scholars have paid attention to energy storage characteristics of rocks [1][2][3][4]. Affected by excavation, the stress and strain will redistribute, resulting in rock damage and destruction. Kong et al. gave a new method to estimate the pressure arch formation above underground excavation in rock mass [5]. Also, the distribution characteristics of stored energy will change. Some energy will release to damage the host rock. Gao et al. studied the rock damage based on the energy principles [6]. Xiao et al. raised a new method for calculating energy release rate in tunnel excavation subjected to high in situ stress [7]. Sometimes, the intense energy releasing process will cause dynamic disasters, especially rockburst. Recently, energy theory has been regarded as the most effective approach to explain the mechanism of deep dynamic disasters and has attracted attentions of scholars around the world. Chen et al. studied the rock burst intensity classification based on the radiated energy with damage intensity at Jinping II Hydropower Station, China [8]. Xu et al. established a new energy index for evaluating the tendency of rockburst and its engineering application also in Jinping II hydropower station [9]. Wang et al. did the research on the energy criterion for rockbursts induced by broken hard and thick rock strata [10]. Miao et al. predicted the rockburst based on in situ stress and energy accumulation theory [11].
Meanwhile, drill and blast is still the main excavation way for the rock mass tunneling [12][13][14]. Under the action of blasting wave, there are many negative effects, such as host rock damage [15][16][17], high initial stress release [18], intense stress and strain [19,20], and blasting vibration of surrounding rock mass [21,22], to name a few. Also, due to blasting effect, the energy accumulation degree is much higher [23], so the possibility and intensity of rockburst and other dynamic disasters will increase accordingly. Also, some scholars have noticed the relationship between rockburst intensity and rock blasting. Fan et al. studied the influence of tunneling methods on the strainburst characteristics during the excavation of deep rock mass [24]. Liu et al. put forward a method for dynamic risk assessment and management of rockbursts in drill and blast tunnels [25]. However, these methods mentioned above are mainly based on theoretical derivations or numerical simulations. e superposition effect and degree cannot be quantitatively given. Hence, based on the in situ geostress data and energy calculation method of circular tunnels, the accumulated energy of three-centered arch roadway in Sanshandao gold mine after quasi-static excavation was derived. en, a series of modified SHPB dynamic tests were conducted after employing the equivalent theory of elastic vibration boundary and the monitoring data of near field blasting impact velocity. So, the energy attenuation in the long granite specimens (400 mm in length and 50 mm in diameter) and blasting vibration energy dispersion with different incident bars were confirmed. Finally, according to the superimposed energy, the tendentiousness of rockburst was evaluated qualitatively. So, the effect of blasting vibration on rockburst intensity was quantified, which provides theoretical guidance and supporting data for dynamic disasters of deep rock under blasting vibration conditions.

Engineering Cases and Geostress
Measurement. e largest exploitation depth in Sanshandao gold mine has already exceeded 1400 m. Rockburst and other dynamic disasters are becoming more and more serious. To study the rockburst tendency, a 1000 m depth roadway was selected as the research cases. e cross-section diagram of excavated roadway is displayed in Figure 1. e dimensions of the roadway are 4.2 m in height and 4.5 m in width. As shown, region I is the roof and the circular arc area of side wall, with the arc radius 2.25 m; region II is the transition region of side wall and the floor, with the arc radius 0.94 m; and region III is the circular arc area of floor, with the arc radius 4.6 m.
A large number of studies have shown that deep high geostress is the essential cause of dynamic disasters; to obtain accurate geostress measurement data, the stress relief method was employed to measure the in situ geostress, and the adopted in situ stress probe is shown in Figure 2.
In order to acquire the accurate distribution law of geostress, in situ geostress measurements of 11 groups were carried out. After measuring, it is found that all curves show similar stress characteristics. As relief length increases, the measured strain values increase constantly and finally become stable. Typical measuring data curves of geostress are shown in Figure 3. In the legend of Figure 3, A, B, and C are the three sets of strain rosette and 0, 45, 90, and 135 are the direction of strain gage in a strain rosette. e stress relief method is to reverse the stress value according to the measured value of strain recovery value with the basis of generalized Hooke's law, so the elastic modulus and Poisson's ratio of rock cores are necessary for data processing; they can be acquired by the confining pressure rating test (shown in Figure 4(a)), and the following formulas are obtained: where E is the elastic modulus; v is the Poisson ratio; P c is the confining pressure; ε θ and ε z are radial and axial strain, respectively; and L and l are outer and inner diameter of the rock core, respectively. Considering the influence of colloid, formula (1) can be modified as follows: where K 1 is the compensation factor. When E ≥ 50 MPa, K 1 is 1.12.
After in situ strain monitoring, confining pressure calibration tests were carried out on the drilled cores; typical confining pressure rating curves are displayed in Figure 4(b). After analyzing the measured data of 11 points, the value and direction of geostress of Sanshandao gold mine are listed in Table 1, and the fitted results are illustrated in Figure 5.
As listed in Table 1, for the geostresses in three directions, the dip angles of maximum principal stress (σ 1 ) are between −10°and 4°; this principal stress is close to the horizontal direction; it is named as the maximum horizontal principal stress (σ H ). For the intermediate principal stress (σ 2 ) and minimum principal stress (σ 3 ), one of them is in the near horizontal state and the other one is in the near vertical state by the dip angle analysis. e geostress with near horizontal state is named as the minimum horizontal principal stress (σ h ), and the geostress with vertical dip angle is the vertical principal stress (σ v ).
ese three principal stresses all increase linearly with depth, and the maximum horizontal principal stress increases much more rapidly. e relationships between three different principal stress and depth are as follows: where H is the depth.

Strain gauge rosette
Controlling system and data logging system Rehargeable battery

Analysis of Energy Accumulation in Excavation.
To avoid the effect of tectonic stress, most of the underground mine main roadway will be designed and excavated along with the direction of the maximum horizontal stress. So, the maximum horizontal principal stress will be released during excavation. When the roadway is excavated along the maximum horizontal principal stress, the vertical principal stress is very close to the minimum horizontal principal stress based on the geostress monitoring results. So, the excavated roadway can be equivalent to an infinite threecentered arch roadway, and the far field geostress effect can be approximate to hydrostatic pressure, which means P 0 � σ v ≈ σ h . Using the plane strain problem method, the confining pressure of circular roadway can be expressed as follows [26]: where σ 1 , σ 2 , and σ 3 are the first, second, and third principal stress, respectively; R is the radius of roadway; and r is the distance between the surrounding rock to the roadway center. en, the strain energy density of surrounding rock is Under hydrostatic pressure conditions, the stresses in the far field have the rule as Substituting equations (5) and (7) into equation (6), the strain energy density of surrounding rock before and after excavation can be calculated as follows: where U 0 and U w are the strain energy density of surrounding rock before and after excavation, respectively. For any point N(r, θ) in the three-centered arch roadway, according to the region partition in Figure 4, when θ ∈ (−α, α), which means N(r, θ) is in the region I, the strain energy density of surrounding rock after excavation is When N(r, θ) is in the region II,   Shock and Vibration When N(r, θ) is in the region III,

Shock and Vibration
Substituting in situ geostress data, elastic modulus, and Poisson's ratio into formulas (9)- (14), the quasi static energy distribution of host rock mass can be fitted, as shown in Figure 6.
It can be seen from Figure 6 that after excavation, the maximum accumulated energy on the top and both sides of the roadway is 89.7 kJ/m 3 . For the bottom of the roadway, the value is 143.6 kJ/m 3 . e maximum values all appear in the host rock at the location 5 m from the free face.

Confirmation of Equivalent Elastic Vibration Zone.
Blasting in the rock mass is a very complicated process; the energy produced by the explosion of an explosive can be roughly divided into two parts: one is for the excavation of the target rock mass and the other can produce the rock mass vibration at a certain distance. Under these combination effects, different damage zones are formed. Based on the damage types and degrees, the recognized zones are blasting crushed zone, blasting fragmentation zone, and elastic vibration zone, as shown in Figure 7, As seen in Figure 7, some of the energy of the explosion will be consumed in blasting crushed zone and blasting fragmentation zone to excavate target rock mass. e remaining energy will propagate outward in the form of a vibrational wave; this part energy will have a superimposed effect on the accumulated energy of rock mass. erefore, for the analysis of energy superposition effect by blasting vibration, an equivalent elastic vibration boundary (region 4 in Figure 7) could be adopted. In this region, the blasting wave can be recognized as the elastic wave gradually decays along the radius [27]. When single-hole blasting is used for infinite rock, the failure mode of rocks in blasting crushed zone is mainly shearing failure. While for the rocks in blasting fragmentation zone, it is tensile failure. e radius can be calculated by the following equation [28]: where r 0 is the radius of blasting drill hole; r 1 is the radius of blasting crushed zone; r 2 is the radius of blasting fragmentation zone; c p is the velocity of P wave in the rock; σ c is the dynamic compressive strength; σ t is the dynamic tensile strength; σ * is the static compressive strength; P is the blasting loading; α is the attenuation coefficient of vibration wave; and ρ r is the density of rock.
Calculations suggest that the radius of single-hole blasting crushed zone is 4 times the radius of blasting hole. For the blasting fragmentation zone, the value is ten times. erefore, in this paper, 10 times of blasting hole radius can be equivalent to elastic vibration zone [29].

In Situ Blasting Vibration Wave
Monitoring. For the in situ blasting, the drilled blasting hole is 40 mm; according to the equivalent elastic vibration zone model, the range of blasting crush zone in the rock mass is 160 mm around blasting hole and the range of blasting fragmentation zone is about 400 mm around blasting hole. It also means that the blasting energy will propagate outward in the form of vibration waves out of the range of 0.4 m. To acquire accurate velocity of blasting vibration wave, Ubox-5016 blasting vibration monitor (see Figure 8) was used in in situ blasting monitoring.
Because it is close to the roadway working face, it is difficult to conduct in situ monitoring. Here, 7 groups of near-field blasting vibration monitoring were applied, and the monitoring location and monitoring results are shown in Table 2.
After data analysis, the attenuation law of vibration velocity with blasting distance could be fitted as shown in Figure 9.
As shown in Figure 9, the peak velocity of blasting decreases gradually with the increase of the propagation distance in rock mass, and the total vibration velocity will attenuate by the formula in Figure 9.

Rockburst Tendentiousness Analysis
According to the value of elastic energy, the rockburst intensity can be divided into four degrees, as shown in Table 3.
As calculated in Section 2.2, without considering the explosion vibration effect, the maximum accumulated energy on the top and both sides of the roadway is 89.7 kJ/m 3 after excavation. For the bottom of the roadway, the value is 143.6 kJ/m 3 .
ese two values indicate that the rockburst intensity is slight and medium, respectively. e maximum values all appear at about 5 m from the free face.
At 5 m position in the host rock, the velocity caused by blasting is 3.86 m/s which is calculated by the formula in Figure 9. Taking the theorem of kinetic energy, the energy density per unit volume of rock mass could be calculated by considering the granite density 2712 kg/m 3 , and the blasting vibration energy density is 40.41 kJ/m 3 . en, the accumulated energy has risen to 130.11 kJ/m 3 on the top and both sides of the roadway and 184.01 kJ/m 3 on the bottom of the roadway. e accumulated energy increase rates are 45.1% and 28.2%, respectively, and the rockburst intensity has risen to medium. It can be seen that blasting vibration obviously increases the probability and intensity of the rockburst, which increases the potential threat of dynamic disasters at a certain extent. 6 Shock and Vibration

Conclusions
Considering the superimposed effect of excavation accumulation energy and blasting vibration energy in deep roadway, (1) Based on the in situ geostress data and energy calculation method of circular tunnels, the accumulated energy of three-centered arch roadway after quasistatic excavation was derived. e maximum accumulated energy on the top and both sides of the roadway is 89.7 kJ/m 3 after excavation. For the bottom of the roadway, the value is 143.6 kJ/m 3 . e maximum values all appear at about 5 m from the free face. (2) Based on near-field monitored data, the attenuation law blasting vibration velocity could be fitted, and the result shows that the peak velocity of blasting will decrease gradually with the increase of the propagation distance.
(3) After considering blasting vibration, the accumulated energy has risen to 130.11 kJ/m 3 on the top and both sides of the roadway and 184.01 kJ/m 3 on the bottom of the roadway. ere is an increase of 45.1% and 28.2%, respectively. It can be seen that blasting vibration obviously increases the probability and intensity of the rockburst, which increases the potential threat of dynamic disasters at a certain extent.

Data Availability
e datasets used in the present study are available from the corresponding author upon reasonable request.