Prediction of Vibration Velocity of Bench Blasting Reflecting Negative Elevation Effect

In the field of mine blasting practice, accurate prediction of blasting vibration is considered as a critical task. In accordance with the theory of elastic wave, the reflection of blasting vibration wave at the interface of negative elevation is analyzed in this paper. The negative elevation effect exerted by blasting vibration is interpreted considering the mechanism, and the formula of the blasting vibration prediction step (negative) 
 
 V
 =
 K
 
 
 
 
 
 
 Q
 
 
 1
 /
 3
 
 
 /
 R
 
 
 
 
 α
 
 
 ⋅
 
 
 
 
 1
 /
 H
 
 
 
 
 β
 
 
 
 is derived, reflecting the negative elevation effect. Besides, the formula accuracy is verified by the measured vibration velocity of the mass in the Jinou coal mine. The step (negative) formula acts as a more feasible candidate for the prediction of step blasting vibration.


Introduction
In large-scale deep open pits, the blasting scale is large on the whole, and the blasting operations are relatively frequent. Accordingly, the vibration effect exerted by blasting is a critical factor affecting the safety of blasting and the mass stability of slope rock. Thus, the blasting vibration effect predicted accurately lays a basis for controlling the blasting vibration effect. Such prediction is a vital task in the practice of mine blasting, as well as a crucial problem in the research of the technology for blasting vibration control.
In the blasting for strip mining, the elevation effect attributed to the blasting vibration wave propagation along the slope of the stope is considered in the study of the law of vibration attenuation. In recent years, the elevation effect of blasting vibration has been extensively studied in-depth. According to the studies by Lü et al. [1], the convex surface above the source can significantly amplify the blasting vibration. Besides, the larger the height difference, the more obvious the blasting vibration will be amplified. Shu et al. [2] found that the amplifying effect of elevation is associated with the slope. When the slope ratio is above 1/2, the slope will lead to the amplification of the elevation; otherwise, there will be no amplifying effect. Tang and Li [3] considered that the open pit slope can amplify the blasting vibration of the explosion source underneath the foot of the slope, and the vibration speed is upregulated unevenly with the rise in the height. When the critical height is reached, the amplification factor will reach its peak. When the height rises further, the amplification factor will gradually decrease. Chen et al. [4] considered that the amplifying effect of elevation on the blasting vibration in an open pit is affected by the shape of the slope and the dynamic load of the blasting. When the slope rises, or when the shape varies, the vibration intensity of the upper step may be greater than that of the lower step. Through experimental observation and analysis, Wan et al. [5] considered that the amplifying effect of the step terrain blasting vibration is affected by the whiplash effect and the slope effect. Zhang et al. [6], Hu et al. [7], and Zhou and Li [8] analyzed the physical quantities associated with blasting vibration. Then, they conducted the dimensional analysis to derive the blasting vibration formula reflecting the amplifying effect of elevation and applied it in engineering. In the meantime, the dynamic finite element method was successfully introduced in the blasting dynamic simulation, and it was proved that the vibration velocity of the particle at the top of the slope exerts an amplifying effect [9]. The amplification phenomenon is reflected as a local dynamic response [10][11][12][13][14][15][16][17][18]. However, most of the above studies on elevation effects were conducted by field monitoring and using numerical models. The studies on the mechanism of the elevation effect have been mostly the qualitative description. The blasting vibration prediction formulas using the dimensional analysis are different in form, exhibiting poor theoretical stringency, and the elevation effect cannot be explained from the mechanism. In accordance with the theory of elastic wave, the reflection of blasting vibration wave at the interface of negative elevation is studied in this paper from the angle of blasting vibration elastic wave, and the negative elevation effect of blasting vibration from its mechanism is explained. The blasting vibration step (negative) prediction formula reflecting the negative elevation effect was proposed in this study, and the accuracy of the formula was verified by the actual measured particle vibration velocity in Jinou openpit coal mine.

Mechanism Analysis of Negative Elevation Effect of Blasting Vibration
It is assumed that the underground rock formation is isotropic and the explosion source is spherically symmetric. Blasting seismic waves propagate in the formation in the form of spherical waves. At the same time, the ground is assumed to be a free surface. After the seismic waves reach the ground, they are reflected and superimposed on the ground. By comparing the reflection of the flat ground model and the step model on the free surface, the negative elevation effect of blasting vibration is explained. The main theory used here is the theory of elastic spherical wave propagation in uniform and isotropic media. There is also the reflection theory of elastic waves on free surfaces. In the process of elastic wave reflection analysis on free surface, a representative longitudinal (P) wave is analyzed.

Propagation of Blasting Vibration Elastic
Wave. The seismic wave attributed to blasting refers to an elastic wave propagating in a subterranean formation [19], so the source of the explosion is assumed to be spherically symmetric. When the seismic wave propagates inside the rock mass and does not reach the ground, it can be considered a divergent spherical wave. If the spherical wave consists of a P wave component and a S wave component, set u p , u s as the displacement field, respectively, and φ and ψ still represent the potential of the displacement field, and Take the P wave as an example. It yields The displacement vector field is defined as where A denotes the amplitude of ψ at the source, equal to a constant; w is the circular frequency of the wave; k = w/c = 2π/λ represents the circular wave number; r is the distance from the observation point to the source. When r > >λ, 1/λ ≥ 1/r. Because two of the brackets in Equation (3) are known, 1/r 2 can be omitted. Then, Equation (3) is written as Equation (4) expresses the simple harmonic of the spherical displacement, indicating that the amplitude of the displacement at a distance away from the source is inversely proportional to the distance from the source.
The magnitude of the wave displacement can be calculated by the above equation: n S 2 X ZF igure 1: P wave incidence.

Geofluids
The above equation was used to derive the time, i.e., the particle velocity at the observation point is expressed as 2.2. Reflection of Elastic Waves on Free Surfaces. The free surface acts as a special interface, i.e., the boundary between the medium and the vacuum. Thus, there exists no medium on one side of the interface. In practice, the earth is commonly considered a free surface [20]. Figure 1 shows that taking the x-y plane as a free surface, a P wave is incident from the lower medium onto the free surface. Since there exists no medium on the surface, when the wave reaches the boundary of the medium, it can only fold back into the original medium to propagate without passing through it, i.e., there is only a reflected wave rather than a transmitted wave. When the P wave is incident on the free surface, it will cause the displacement along the normal direction of the surface and along the tangential direction. Thus, the reflected wave consists of both P wave and SV wave. In the meantime, according to the symmetry of the problem and the independence of the P wave, the SV wave, and the SH wave, the reflected wave contains no SH wave component (e.g., the incident wave). The displacement vectors of the incident P wave, the reflected P wave, and the reflected SV wave are expressed as where w denotes the circular frequency of the wave; V P and V S are the speeds of P and S waves, respectively; In accordance with the propagation law of the elastic wave at the interface, on the free surface z = 0, the magnitude of the displacement of each wave is expressed as The above displacement formula was adopted to calculate the time t, and the calculated velocity of the particle at the observation point of each wave is expressed as follows: Based on the geometric relationship, the reflected P wave and the reflected SV wave are superimposed on each other at the free surface, and the velocity in the x-direction and the z-direction is expressed as

Mechanism Analysis of Negative Elevation Effect of
Step Blasting Vibration. To study the negative elevation effect mechanism of step blasting vibration, a geometric model was built, as shown in Figure 2. Assume that the step model exhibits a high steep shape, there is a certain distance of the observation point from the explosion source, and the depth of the explosive is relatively shallow. In the step model, the distance between the blasthole and  3 Geofluids the edge of the step is almost equal to the depth of the explosive. The depth of the medicine package is denoted as h, the step height is H, the distance from the medicine bag to the edge of the step is l, and the horizontal distance of the observation point from the drug pack is x. According to the Huygens-Fresnel principle [21], each point on the wavefront (wave front) can be considered the wave source of the transmitted wavelet in the propagation process. At any subsequent time, the envelope surface of these wavelets can act as a new wavefront. The envelope surface is formed by the interference of the wavelets generated by the points on the wavefront at the observation points. In the meantime, based on the Fermat principle [22], the wave always propagates along the path with the smallest propagation time in the medium. These paths are rays, and the rays are straight lines in a uniform medium. Subsequently, from the explosion source to the observation point E at the same horizontal distance from the explosion source, the seismic wave propagation in the step model According to Equation (10), in the step model, after the elastic waves are superposed on each other at the free surface E, the particle velocity of the point can be, respectively, expressed in the X-direction and the Z-direction as follows: In the flat model, after the elastic waves are superimposed on each other at the free surface E, the velocity of the point particle is expressed, respectively, in the X -direction and the Z-direction as follows: where α 1 , β 1 denote the incident angles and reflection angles of the waves at a certain point in the step model. α 2 , β 2 represent the incident angle and reflection angle of the wave at a certain point in the flat model, respectively. In the step model, compared with the flat model, the vibration velocity magnification in the X-direction and the Z-direction is expressed as Substituting Equations (11) and (12) into the above equations, they are simplified according to Equation (9) as where and SVwave reflection coefficients, respectively.
In the step model, according to Equation (6), the velocity of the particle at the observation point E (where the elastic wave does not reach the critical point of the free surface) is expressed as Subsequently, the amplitude of the particle vibration velocity is expressed as In the flat model, according to Equation (6), the velocity of the particle at the observation point E (where the elastic wave does not reach the critical point of the free surface) is expressed as Next, its amplitude of particle vibration velocity is expressed as Then, The relationship between the velocity ratio of elastic longitudinal wave and transverse wave and Poisson's ratio v is defined as [23] In accordance with Snell's law [24], it yields Set v = 0:25 According to the relationship between Equations (15) and (16) and other parameters, through the MATLAB math software drawing, the relationships between the X-direction and Z-direction magnification and the horizontal distance and vertical distance from the explosion source were plotted, as shown in Figures 3 and 4, respectively. Figure 3 shows that the magnification of the vibration velocity in the X-direction is above 1, and the remaining magnifications are all no less than 1, at a horizontal distance of 50-65 m from the source of the explosion. Accordingly, the existence of the negative elevation of the step model is primarily reflected in the decreasing effect on the vibration velocity. Figure 4 shows the Z-direction particle vibration velocity magnification is extremely small, almost zero. This reveals that the vibration velocity of the particle is significantly downregulated because of the existence of the negative elevation of the step model. In the meantime, Figure 4 suggests that the larger the step negative elevation, the smaller the amplification factor of the vibration speed will be. Figures 3 and 4 indicate that the magnification of the vibration velocity of the step negative elevation is inevitably associated with the horizontal distance and the elevation difference.

Prediction Formula of Blasting Vibration Velocity Reflecting Negative Elevation Effect
From the above, the vibration velocity ratio of the measuring points with equal horizontal distance and different vertical distance is obtained, but it is more complicated and needs to be simplified. Then, combined with the Sadowski formula, the vibration velocity prediction formula reflecting the negative elevation effect was deduced.

Negative Elevation Effect Magnification Analytic
Simplification. In large-scale deep open pits, the general step slope is relatively high and steep. Besides, there is a certain In the step model, α 1 = 90 ∘ , so it yields In the flat model, the burial depth h of the blasthole is smaller than the distance X of the observation point from the center of the explosion source. Assume h < 6h ≤ x.
Since the value of α 2 varies slightly, to simplify the calculation, set α 2 = 86 ∘ . Then, it yields According to the conditions of the initial hypothesis, the step is high, and there is a certain distance of the observation point from the source. Thus, the depth of the explosive is smaller than the distance from the observation point to explosive source, which is negligible. Furthermore, the distance between the depth of the explosive is nearly equal to the length of the blasthole from the edge of the step, which is expressed as Then, Because the step is high steep, and the depth of the explosive is associated with the elevation difference H, which is relatively small, the depth is neglected and further simplified as Substituting Equations (27), (28), (29), and (32) into Equations (15) and (16), it yields Equations (33) and (34) are unified, and the parameter k 1 is introduced, i.e., the amplification factor R cT of the step model negative acceleration vibration speed is expressed as To reduce the error caused by the approximation during the above simplification, the parameter β was introduced, and Equation (35) is written as

Prediction Formula of Vibration Velocity Reflecting
Negative Elevation Effect. The vibration velocity at the blasting observation point of the negative elevation step model is denoted as V 1 . According to the Sadowski formula, the formula for predicting the vibration velocity of the observation point widely used in flat terrain is written as [22] where V is the ground particle peak vibration velocity, cm/s; Q is the explosive quantity (the total charge amount when bursting together, the maximum charge amount when delaying blasting), kg; R is the distance from the observation point to the source, m; k, α are the coefficient and attenuation coefficient associated with the topography and geological conditions between the blasting point and the calculation point.
Subsequently, the vibration speed V 1 of the negative elevation step model was compared with the vibration speed of the flat terrain, and the magnification is expressed as Then, according to Equation (36), Equation (37) is rewritten as Next, Equation (39) can be rewritten as Equation (41) is the prediction formula for the negative elevation blasting vibration velocity of the step model. Since R represents the explosion source distance in this formula, the above equation can also be rewritten as where V is the negative elevation step ground particle peak vibration velocity, cm/s; Q is the amount of explosives (the total charge when exploding together, the maximum charge during delayed blasting), kg; D is the horizontal distance from the observation point to the source, m; H is the vertical distance from the observation point to the source, m; k, α are the coefficient and attenuation coefficient associated with the topography and geological conditions between the blasting point and the calculation point. β is the error coefficient    Table 1.
The data in Table 1 suggests that the distribution law of the vibration velocity with the horizontal distance of the blasting heart can be concluded, as shown in Figures 5 and 6. Figure 5 shows that compared with 1#, 2#, 3#, and 4# measuring points, the vibration velocity of 2# and 4# measuring points rises; compared with 5# measuring points, particle vibration speed of 4# measuring points decreases significantly. Step shape Explosive source 2# Figure 6: Monitoring data on October 16.
9 Geofluids Figure 6 suggests that compared with 5#, 6#, and 7# measuring points, the particle vibration speed decreases slowly and keeps basically flat. The difference in vibration speed is the largest between 2# and 3#. Figures 5 and 6 indicate that the vibration velocity of the particle gradually decreases with the rise in the horizontal distance, whereas the existence of the elevation difference leads to a more significant decrease in the vibration velocity of the particle.

Discussion
From the above, not only the vibration velocity prediction formula reflecting the negative elevation effect was obtained but also the field measurement of blasting vibration was carried out. Next, in order to test the accuracy of the prediction formula, a regression analysis was performed on the field measured data. Then, the prediction formula obtained in this paper is compared with other improved prediction formulas for error analysis.

Regression Analysis of Prediction Formula Reflecting
Negative Elevation Effect Velocity. Some researchers considered that the commonly used Sadowski formula only considers the vibration velocity V, which decreases as the explosion heart distance R rises, without considering the effect of elevation on the blasting vibration effect.
Accordingly, some optimized formulas are proposed, such as the Sadowski formula (space distance formula) [25].
where D denotes the horizontal distance between the observation point and the source (m); H is the vertical distance between the observation point and the source. However, literature [26,27] considers the elevation effect of blasting vibration wave propagation, and the optimized particle vibration velocity prediction formula is written as Zhu et al. [28] considered that the vibration velocity has an amplifying effect along the elevation, so the calculation of the vibration velocity is defined as where R denotes the explosion source distance, i.e., the distance between the center of the explosion source and the measuring point (m), which can be written as Tang and Li [29] optimized a calculation formula reflecting the amplifying effect of elevation by the dimensional analysis: The regression analysis was conducted on the aboveoptimized vibration velocity calculation formula reflecting the elevation effect with Equation (42), and the achieved regression results are listed in Table 2. Table 2 shows that the Sax (space) formula exhibits the lowest fitting correlation coefficient, the largest degree of difference, and the smallest correction coefficient, indicating that the fitting effect of the formula is the worst. However, 10 Geofluids the fitting correlation coefficient, the degree of difference, and the correction determination coefficient of the optimized 2 formula, the optimized 3 formula, and the step (negative) formula are highly consistent with each other.

Comparison of Blasting Vibration Speed Calculation
Accuracy. Based on the regression results achieved in the previous section, the relative error and the average relative error of the monitoring data of Table 1 are compared with the predicted results of the vibration velocity calculation formula reflecting the elevation effect. The comparison results are listed in Table 3 and Figure 7. Figure 7 suggests that the Sax (horizontal distance) formula exhibits the largest prediction error, reaching 57.73%, followed by the Sa's (space distance) formula, with an error of 44.2%. The step formula (negative) exhibits the smallest error of 33.12%.

Conclusion
Aiming at the height effect of bench blasting vibration, this paper uses theoretical research, field measurement, and other methods to carry out research. Established a geometric model of the steps and analyzed the reflection of vibration waves on the free surface of the steps. Conducted on-site monitoring and statistical of blasting vibration. The main conclusions of the paper are as follows: A vibration propagation model of bench blasting was established, which revealed the mechanism of elevation effect of blasting vibration wave propagating in bench terrain, and provided a more effective new idea for the prediction of blasting vibration.
Based on the reflection theory of elastic waves on a free interface, a calculation formula for predicting the peak vibration velocity of surface particles, the step (negative) formula is proposed. It is found that this formula exhibits a significantly lower error of the prediction results than the commonly used Sa's formula at home and abroad and other optimized formulas.
The feasibility and accuracy of the method to predict the peak velocity of the negative elevation surface in the step blasting project were verified in this study. However, in the positive elevation difference surface and underground blasting engineering, the prediction of the peak velocity of the particle requires further verification.

Data Availability
The test data used to support the findings of this study are included within the article. Readers can obtain data supporting the research results from the test data table in the paper.

Conflicts of Interest
There is no conflict of interest to declare.