# The Vibration Characteristics of Ground of Rock Blasting in Silt-Rock Strata

Authors: Tao Yin; Chuanbo Zhou; Changqing Zheng; Jugen Fu; Ziru Guo
Journal: Shock and Vibration
 (2021)

---

## Abstract

As the shield section passes through the silt-rock strata, the rock stratum of the tunnel section has to be blasted into blocks in advance, and the diameter of the blocks should be less than 30 cm after breaking, and then, the blocks could smoothly enter the soil cabin through the opening of the tunnel boring machine (TBM) cutter head and finally be discharged through the screw machine. The geology of rock blasting in silt-rock strata is complex, and the vibrations caused by blasting threaten the safety of buildings around the blasting area. According to the measured data of blasting vibrations at the sites, the waveform duration of vibration acceleration and the distribution characteristics of dominant frequency of vibration velocity were analyzed, the energy characteristics of vibration velocity were researched by wavelet analysis, and the attenuation laws of vibration velocity were studied by dimensional analysis (DA). The dominant frequency bands of vibration energy of ground are in the range of 0–15.625 Hz, and the distribution characteristics of frequency bands of vibration energy in different directions of the ground are similar to each other, but the energy magnitude is different from each other. The research results could provide a reference for the safe blasting distance of buildings under similar geological conditions.


---

## Body

## 1. Introduction

With the fast economic growth in China, the demand for urban underground space engineering is growing rapidly. By 2020, there have been 33 cities with 177 subway lines in China [1]. In the construction of metro tunnel shield method, the silt-rock strata often encounter in the range of section, such as Shenzhen Metro Line 11 [2], Xiamen Metro Line 2 [3], Shenzhen Metro Line 5 [4], and Fuzhou Metro Line 1 [5], and when the shield machine passes through the silt-rock strata, the rock stratum in the lower part is not easy to break, while the silt stratum in the upper part is easy to enter into the soil cabin, especially the silt stratum, which is rich in groundwater. Once the soil pressure in the sealed soil cabin is unbalanced, it is easy to cause soil loss, large settlement, and collapse accident. Many scientists had researched to solve this problem [1–5]. They concluded that when the rock stratum of the tunnel section was blasted into blocks in advance, and the diameter of the blocks should be less than 30 cm after breaking, then the blocks could smoothly enter the soil cabin through the opening of the TBM cutter head and finally be discharged through the screw machine.There are many studies about blasting in rock, soil, air, and water, but few researches about rock blasting in the silt-rock strata. The removing and replacing silt layer by explosion has been investigated comprehensively in China. For example, Yang et al. [6] studied the propagation laws of blast waves in silt, which were similar to that in water and satisfied the geometrical similarity law. Li et al. [7] researched the vibration effect of removing and replacing silt layer by explosion in deep silt and obtained that the high tide level could reduce the negative effect of blasting vibrations and the low tide level had a positive impact on the peak velocity of blasting vibrations. Xu et al. [8] applied the regression analysis by Sadovski's formula and obtained the attenuation laws of peak particle velocity (PPV) with distance in the removing and replacing silt layer by explosion. W. Wei [9] adopted many blasting tests in soft marine silt clay at the sites of Ningbo and obtained the time history of blast wave pressure, excess pore water pressure, and temperature of soft soil under different buried depths of explosive. Some foreign researchers had conducted a lot of experiments to compact soil by blasting. C. WA et al. [10] studied that the blasting vibrations could cause the liquefaction of saturated soil; however, the liquefaction caused by blasting load was different from that caused by seismic load. Mikulyak [11] modeled the blasting by computer in geological conditions of soft soil and investigated the formation characteristics of two kinds of spherical cavities by different blasting schemes.Some scientists studied the blasting vibration waveform by wavelet analysis. D. Huang [12] analyzed the time-frequency characteristics of measured vibration signal of the mountain tunnel by the Fourier transform and the wavelet packet transform. Y. Yang [13] studied the blasting vibration signals of different elevations and explosion distances of an open-pit mine based on the wavelet packet theory. Józef Pyra [14] used various methods of time and frequency to analyze the vibration induced by blasting works carried out in open-cast mines.Other scientists have studied the propagation law of blasting vibration by DA. M. Khandelwal [15] performed a DA technique to propose a new formula for the prediction of the PPV. Dehghani [16] evaluated and calculated the blast-induced ground vibration by artificial neural networks and DA techniques. Olarewaju [17] used the finite element method to study the dynamic dimensionless behaviors of underground pipes due to blast loads.The analysis of vibration characteristics of ground of rock blasting in silt-rock strata by wavelet analysis and DA has been less reported at present. This study relied on the rock blasting engineering in silt-rock strata in the shield section of Zhuhai Metro, the waveform and the distribution characteristics of frequency of blasting vibrations were analyzed, the energy characteristics of blasting vibrations in different frequency bands were studied by wavelet analysis, and the attenuation laws of vibration velocity were fitted by DA. The dominant frequency bands of vibration energy of ground are in the range of 0–15.625 Hz, and the distribution characteristics of frequency bands of vibration energy in different directions of the ground are similar to each other, but the energy magnitude is different from each other. The research results could provide a reference for the safe blasting distance of buildings under similar geological conditions.


## 2. Project Overview

The shield section between Hengqinbei Railway Station and Hengqin Railway Station of Zhuhai-Zhuhai Airport Intercity Railway is located in the Economic Development Zone of Hengqin Island, Zhuhai City, as shown in Figure1. It runs along the coastline of South China Sea, with a total length of 2.656 km, all of which is constructed by the shield method.Figure 1
Zhuhai-Zhuhai Airport Intercity Railway.As the shield section passes through the silt-rock strata, the rock stratum of the tunnel section has to be blasted into blocks in advance, and the diameter of the blocks should be less than 30 cm after breaking, and then, the blocks could smoothly enter the soil cabin through the opening of the cutter head and finally be discharged through the screw machine. The diagram of charging on field is shown in Figure2. The rock cores before blasting and after blasting are shown in Figure 3.Figure 2
Diagram of charging on field.Figure 3
Rock cores (a) before blasting and (b) after blasting.
(a)(b)The blasting area is located in the artificial reclamation area, and there are three strata from the ground surface to the subway tunnel. (1) The first stratum is artificial fill. It is yellowish-brown, slightly wet, mainly in saturated loose state, and slightly compacted, and the artificial fill is mostly pavement concrete with a thickness of 20–35 cm, and the layer thickness is in the range of 0.30–15.10 m, and the average layer thickness is 2.82 m. (2) The second stratum is silt. It is from gray to dark gray, saturated, and in flow plastic state, with uniform and fine soil, locally enriching with shell debris and sand. The layer thickness is in the range of 0.50–30.30 m. (3) The third stratum is rock. It is weakly weathered granite. It is brown-yellow, gray-brown, gray-white, and other colors, medium-coarse grain structure, and massive structure. The joint fissures are developed relatively. The rock quality is hard. The buried depth of the layer fluctuates greatly, and the distribution is discontinuous. The physical and mechanical parameters of artificial fill, silt, and rock are listed in Tables1 and 2, respectively.Table 1
Physical and mechanical parameters of artificial fill.
Densityρ (g/cm3)Poisson’s ratiouCohesion (KPa)Shear modulus (MPa)Angle of internal friction (rad)1.980.351001500.31Table 2
Physical and mechanical parameters of rock and silt.
Densityρ (g·cm−3)Elastic modulusE (GPa)Poisson’s ratiouP wave velocity (m/s)Yield strength (MPa)Tangent modulus (GPa)Hardening parameterStrain rate parameterCRock2.50.740.2553573.0e-50.0421.00.06Silt1.633.3e-50.4512504.0e-90.0——The length of the blasting area is 960 cm, and the width of the blasting area is 160 cm. The casing drilling method is used for drilling blastholes, and the PVC casing with a diameter of 7.5 cm is installed immediately after the blasthole is formed for protecting the blasthole wall. The spacing and row spacing of blastholes are both 80 cm, the diameter of blastholes is 14.6 cm, and the depth of blastholes is 3400 cm. Three rows of blastholes are detonated each time, 13 blastholes are arranged in a row, and the millisecond blasting network in the blasthole is adopted. The millisecond sections of three rows of blastholes are 6, 7, and 8 in turn, and the delay time is 50 ms. The layouts of blastholes are shown in Figure4.Figure 4
Sketch map of blastholes (cm). (a) Planar graph, (b) cross section of A-A, and (c) cross section of B-B.
(a)(b)(c)After drilling the blastholes, the blastholes are filled with underground water immediately, and the blastholes are water coupling charge. The diameter of the explosive roll is 6 cm, and the density of the explosive is 1 061 kg/m3. The continuous charge is adopted, the charge length is the depth of the rock at the tunnel section, and the thickness of the rock to be blasted is in the range of 1–8 m. The blastholes are stemmed with gravels of 1–1.5 cm diameter, and the stemming length is 1 000 cm.


## 3. Test Points of Blasting Vibrations

The TC-4850 self-recording instrument (developed and manufactured by Zhongke (Chengdu) Instruments Company Limited) was used as the blasting vibrations measuring instrument during the blasting, each instrument had three channels, and it could be collected in parallel. The arrangement of blasting vibration meters is shown in Figure5.Figure 5
Arrangement of blasting vibration meters.Three acceleration measuring points and six velocity measuring points were arranged in total. The measuring points were arranged along the ray with the explosion source as the center. The schematic diagram of explosion source and measuring points is shown in Figure6.Figure 6
Schematic diagram of explosion source and measuring pointsm.


## 4. Analysis of the Characteristics of Blasting Vibrations of Ground

### 4.1. Analysis of the Characteristics of Vibration Acceleration

The first blasting was carried out at 10 : 00 a.m. on March 31, 2016, and the second blasting was carried out at 3 : 00 p.m. on March 31, 2016, with an interval of 10 days, and the third blasting was carried out at 10 : 00 a.m. on April 10, 2016. The vibration acceleration test points were all 10 meters away from the explosion source, and the results of three tests of blasting are shown in Table3.Table 3
Results of three tests of blasting.
The total explosives (kg)The maximum charge of one section (kg)SectionsThe dominant frequency (Hz)The first blasting43.234.85,6,7X47.06Y55.56Z36.04The second blasting432.0150.05,6,7X10.34Y8.32Z14.60The third blasting331.2115.25,7,8X36.70Y37.38Z66.67The waveform duration of vibration accelerations of three blasting is shown in Figure7.Figure 7
Acceleration waveform of blasting vibrations, (a) acceleration waveform of the first blasting, (b) acceleration waveform of the second blasting, and (c) acceleration waveform of the third blasting.
(a)(b)(c)As shown in the acceleration waveform of blasting vibrations from Figure7, the peak particle acceleration of the three blasting vibrations has little difference to each other; the duration of acceleration waveform is about 0.4 s at 10 : 00 a.m. on March 31, 2016, in the first blasting, and about 0.8 s at 3 : 00 p.m. on March 31, 2016, in the second blasting, which is 0.4 s longer than the duration of the first blasting, with an interval of 10 days, and in the third blasting at 10 : 00 a.m. on April 10, 2016, the duration of acceleration waveform is about 0.5 s, which is close to the duration of the first blasting.The above data show that the silt has a good structural strength during the first blasting, but after the disturbance of the blasting load, the silt has been thixotropic liquefied [18]. During the second blasting with a short time interval, the explosive gas generates a cavity in the silt after the explosion. When the pressure of the explosive gas is less than the pressure of the silt outside the cavity, the silt around the cavity begins to compress into the cavity, and then, the cavity comes into being. The cavity is overcompressed, and then, the cavity expands again to work on the silt. In the process of pulsation, the water and stemming materials in the blast hole are continuously pushed out of the blast hole until they reach the orifice and escape when they come into contact with the atmosphere, which results in the extension of the waveform duration of acceleration in the second blasting. During the third blasting, the duration of acceleration waveform is close to that of the first blasting, which indicates that after a period of sedimentation, the silt has recovered a certain structural strength.


### 4.2. Analysis of the Distribution Characteristics of the Dominant Frequency of Blasting Vibration Velocity

Based on the analysis of the dominant frequency of vibration velocity of the three blasting tests in Table3, it is concluded that the dominant frequency of vibration velocity of the first blasting is between 30 and 60 Hz, which indicates that the silt has good structural strength; during the second blasting, the silt undergoes thixotropic liquefaction after the disturbance of the first blasting load. The dominant vibration frequency of dynamic velocity is reduced to 5–15 Hz, which is close to the natural vibration frequency of the building, and it is easy to cause the building resonance. At 10 : 00 am on April 10, 2016, the dominant frequency of the blasting vibration velocity is between 30 and 70 Hz during the third blasting, and the analysis shows that the structural strength of the silt has been restored due to 10 days from the last blasting. Therefore, the continuous blasting operation should be avoided in the field blasting. Each blasting operation should be separated for a period of time, and the next blasting will be beneficial to reduce the vibration damage until the silt recovers structural strength.


### 4.3. Energy Percentage of Each Frequency Bands of Blasting Vibrations

Taking a typical test data of blasting vibrations as an example, the field monitoring data of blasting vibrations are shown in Table4, and the time history curves of velocity of measured blasting vibration signals are shown in Figure 8.Table 4
Field monitoring data of blasting vibrations.
The total explosives (kg)The maximum charge of one section (kg)SectionsThe distance from explosion source (m)PPV (cm/s)XYZ4321505,6,71300.150.130.10Note. X was the horizontal direction pointing to the explosion source, Y was the horizontal direction perpendicular to X, and Z was the vertical direction.Figure 8
Time history curves of velocity of measured blasting vibrations.


#### 4.3.1. Reconstruction Error

The signals of blasting vibrations belonged to random signals. At present, the db8 wavelet basis function was widely used to analyze the nonstationary vibration signal [19–21]. The error of wavelet reconstruction signal based on db8 is shown in Figure 9.Figure 9
Wavelet reconstruction and error based on db8. (a) Original blasting vibration signals and (b) error of reconstructed signals.
(a)(b)It could be seen from Figure9 that using db8 wavelet basis function to decompose the vibration signals in 8 layers, the relative error magnitude between the reconstructed signal and the original signal was 10−11, which met the analysis requirements.


#### 4.3.2. Decomposition and Reconstruction

Using db8 wavelet basis function to decompose the signal of blasting vibrations with the scale of 8, the wavelet coefficients of 9 frequency bands were obtained [22]. During blasting vibration monitoring, the sampling frequency of the signal was set to 8000 Hz (2 s sampling 16000 points). According to the Shannon sampling theorem [21], the Nyquist frequency of sampling was 4000 Hz. According to the principle of wavelet decomposition [23], the bandwidth of nine wavelets was a8: 0–15.625 Hz, d8: 15.625–31.25 Hz, d7: 31.25–62.5 Hz, d6: 62.5–125 Hz, d5: 125–250 Hz, d4: 250–500 Hz, d3: 500–1000 Hz, d2: 1000–2000 Hz, and d1: 2000–4000 Hz. The blasting vibration components of different wavelet frequency bands are shown in Figure 10.Figure 10
Stratified reconstruction of signals by wavelet decomposition.


#### 4.3.3. Percentage of Energy per Frequency Band

The statistics of energy percentage of each wavelet band of blasting vibration signals in different directions are shown in Table5.Table 5
Energy percentage of each frequency band in different directions (%).
Frequency bandsa8d8d7d6d5d4d3d2d1X43.809.6924.9421.120.380.010.010.020.04Y67.021.8111.2519.410.460.010.010.010.02Z88.914.934.131.800.150.030.010.020.03During blasting rock in silt-rock strata, the frequency bands of ground vibration energy are wide, but the frequency bands are mainly concentrated in 0–15.625 Hz, 31.25–62.5 Hz, and 62.5–125 Hz, which are accounting for more than 90% of the total signal energy. The dominant frequency bands of blasting vibration energy are in 0–15.625 Hz, which are close to the natural frequency of buildings and have a great impact on the vibration damage of surface buildings.The high-frequency signals above 125 Hz continue the whole monitoring period, which are represented by noise signal, with relatively small energy, and limited impact on surface buildings.


### 4.4. The Energy of Each Frequency Bands of Blasting Vibrations

Taking the test data of the first blasting vibrations on the morning of March 31, 2016, as an example, the distance between the test point and the explosion source was 60 m, and the energy distributions of each wavelet band in different directions are shown in Figure11.Figure 11
Energy of each frequency bands in different directions (cm/s)2.It can be seen from Figure11 that the distribution characteristics of the frequency bands of blasting vibration energy in different directions of the ground are similar to each other, but the energy magnitudes are different from each other.The wavelet energy analysis was carried out on the six test data of blasting vibrations of rock blasting in silt-rock strata, and the directions of the maximum energy are shown in Table6.Table 6
Directions of the maximum energy.
The distance from explosion source (m)The first blastingThe second blastingThe third blastingThe forth blastingThe fifth blastingThe sixth blasting30XZZZXZ60ZXZZXZ90ZZZZXZ120YYXXZZ150The instrument was not triggeredZZZYZ180ZZXThe instrument was not triggeredZZNote. X was the horizontal direction pointing to the explosion source, Y was the horizontal direction perpendicular to X, and Z was the vertical direction.According to Table6, when the distance of the test points from the explosion source is in the range of 30–60 m, the number of times of maximum energy in Z direction is 8, and the number of times of maximum energy in X direction is 4.According to the analysis, because of the end effect of the cylindrical charge, there is air above the ground, without clamping effect, and there are more times of maximum energy in theZ direction; when the amount of cylindrical charge is less, the end effect is not obvious, and the explosion stress waves propagate outward in the form of cylinder, 13 blastholes in a row produce a cylinder wave, and the cylinder wave would overlap at the measuring point, so there are 4 times in the X direction.When the distance of the test point from the explosion source is in the range of 90–180 m, the number of times of maximum energy is 15 in theZ direction, 4 in the X direction, and 3 in the Y direction. It is considered that the maximum number of times of energy in the X direction and Y direction is the most due to the influence of pinch effect on explosion stress wave in the Z direction.The times of the maximum energy in theZ direction are the most, followed by the X direction when the distance of the test point from the explosion source is in the range of 30–60 m. Because the vibration energy in the Z direction is the largest and easy to cause up and down vibration damage to the building, the blasting vibration energy in the Z direction should be focused on analysis, followed by the X direction.The vibration energy in theZ direction is mainly analyzed when the distance of the test point from the explosion source is in the range of 90–180 m.


### 4.5. Attenuation Laws of Vibration Velocity of Rock Blasting in Silt-Rock Strata

To study the attenuation laws of vibration velocity of rock blasting in the silt-rock strata, four vibration velocity tests were carried out in the blasting area. Six measuring points of vibration velocity were arranged for each test. The test results of blasting vibration velocity are shown in Table7.Table 7
Testing results of blasting vibration velocity.
Test timesThe maximum charge of one sectionQ (kg)Vibration velocityNo. of measuring pointThe distance from explosion source (m)PPVυ (cm/s)TheX directionTheY directionTheZ directionResultantNo. 134.8ZC1301.671.321.931.94ZC2600.370.390.490.54ZC3900.180.170.330.34ZC41300.150.130.10.2ZC51700.040.060.090.1No. 2150ZC1304.683.28.18.58ZC2606.973.874.227.6ZC3900.460.470.971.01ZC41201.111.140.811.36ZC51500.340.390.620.63ZC61800.340.140.480.49No. 3165.6ZC1304.953.115.315.36ZC2602.371.758.058.25ZC3901.220.842.092.09ZC41200.730.360.830.87ZC51500.010.0170.00570.017No. 454ZC1302.130.680.630.63ZC2600.570.680.381.41ZC3901.60.820.51.61ZC41200.370.262.12.12ZC51500.220.181.021.03ZC61800.180.120.630.63


#### 4.5.1. Prediction Model of Blasting Vibration Velocity Based on DA

Natural phenomena and engineering problems could be described by a series of physical quantities. To find the laws between these physical quantities, firstly, we could list the physical quantities that involved in these problems, ignore some physical quantities with less influence, and then classify the physical quantities involved according to their attributes; secondly, we could find out the interrelations between different physical quantities; and finally, we could find out the causal relationship between some physical quantities and the other physical quantities.There were relations among all kinds of physical quantities, which showed that their structure must be composed of several unified basic components. This basic component of physical quantity was called dimension. Because physics studied the evolution and movement of matter in space and time, all quantitative problems could not be separated from the three basic quantities of mass (M), time (T), and length (L). Therefore, M, T, and L were the most suitable dimensions for these three basic quantities. The dimensions of all other derived quantities could be expressed as the combination of these three basic quantities according to the definition or objective laws.The attenuation of seismic wave was affected by the factors of explosion source, geological conditions, distance of explosion source, and depth of explosion source when the blasting seismic wave propagated in silt-rock strata [24, 25]. There were 9 physical quantities involved in propagation velocity in silt-rock strata [26, 27], as shown in Table 8.Table 8
Physical quantities involved in blasting vibration velocity in silt-rock strata.
VariablesDimensionsThe dependent variablesParticle vibration displacementSLPPVυLT−1Peak particle accelerationαLT−2Frequency of particle vibrationfT−1The independent variablesQuality of explosivesQMHorizontal distance between the test point and the explosion sourcerLDepth of explosion sourceHLDensity of rockρML−3Propagation velocity of vibration wavecLT−1Detonation timetTNote. L, T, and M were dimensions of length, time, and mass, respectively.Based on the DA of Buckingham's theorem (π theorem), the PPV v of the ground during the rock blasting in the silt-rock strata could be expressed as follows:(1)v=ΦS,ρ,H,a,f,t,Q,r,c.According to theπ theorem, where independent dimensions were Q, r, and c, and π represents dimensionless quantities, then,(2)π=vQαrβcγ,π1=SQα1rβ1cγ1,π2=ρQα2rβ2cγ2,π3=HQα3rβ3cγ3,π4=aQα4rβ4cγ4,π5=fQα5rβ5cγ5,π6=tQα6rβ6cγ6.where α, β, and γ were undetermined coefficients. According to the dimensional homogeneous theorem, there were as follows:(3)π=vc,π1=Sr,π2=ρQr−3,π3=Hr,π4=ar−1c2,π5=fr−1c,π6=trc−1.Substituting (3) into (1) could obtain the following:(4)vc=ΦSr,ρQr−3,Hr,ar−1c2,fr−1c,trc−1.Since the product and power of different dimensionless numbersπ were still dimensionless, π2 and π3 are combined as follows to get a new dimensionless number π7:(5)π7=π21/3β1π3β2=ρ3rQ3β1Hrβ2,where β1 and β2 were the indexes of π2 and π3, respectively.At the same site,ρ and c could be approximately constant. Therefore, from (5), it could be considered that v and Q3/rβ1H/rβ2 had a functional relationship.To sum up, this functional relationship could be expressed as follows:(6)lnv=α1+β1lnQ3r+α2+β2lnHr.Iflnv0=α1+β1lnQ3/R., there were as follows:(7)lnv0=α1+β1lnQ3−β1lnr.In equation (6), α1 and α2 were the coefficients given in the process of function transformation; −β1ln r was the attenuation of blasting vibration velocity with distance r; β1 was the attenuation index, which mainly reflected the effect of site medium conditions; and α1 + (β1lnQ)/3 comprehensively reflected the contribution of site medium conditions and explosive quantity to the PPV of ground.Iflnk1=lnα1, there were as follows:(8)v0=k1Q3rβ1.Equation (8) was Sadovski's equation. without the effect of the depth of explosion source. Substituting (8) into (6), we could get the following:(9)lnv=lnv0+α2+β2lnHr.Iflnk1=lnα1, then (9) could be changed into the following:(10)v=k1k2Q3rβ1Hrβ2.Ifk=k1k2, the mathematical model of the attenuation laws of blasting vibration velocity in silt-rock strata considering the effect of burial depth was established as follows:(11)v=kQ3rβ1Hrβ2,where k was the site effect coefficient, β1 was the attenuation coefficient of blasting vibrations, and β2 was the effect coefficient of explosion source depth.


#### 4.5.2. Comparative Analysis of Prediction Results of Blasting Vibration Velocity

The explosion source was below the ground surface during blasting rock in silt-rock strata, while measuring points of blasting vibration velocity were often located on the ground surface, which was also the key content of the analysis of attenuation laws of blasting vibrations. Therefore, this study selected the ground data for analysis and established the prediction model of vibration velocity of rock blasting in silt-rock strata. According to Z. C. Zhang’s research [28], the attenuation laws of vibration velocity near blasting zone were different from the attenuation laws of vibration velocity far blasting zone. R¯=r/Q1/3 was the scaled distance, r was the horizontal distance between the test point and the explosion source, and Q was the quality of explosives. The R¯<5.0 was in the near blasting zone. The R¯ of monitoring data in this paper is greater than 5.0, and the monitoring data were not within the range of near blast area. The monitoring data could be used to analyze the attenuation law. Equation (11) in Section 4.5.1 is used for regression fitting analysis of the blasting vibration data listed in Table 7 to obtain the prediction model of vibration velocity of rock blasting in the silt-rock strata.In regression analysis, the correlation coefficient was an indicator to measure the linear correlation between variables. The larger the correlation coefficient was, the smaller the error was, the higher the linear correlation between variables was, the higher the reliability of the prediction model was, and vice versa. To determine the accuracy of the established mathematical model, the established mathematical model (11) was compared with the traditional Sadovski's formula (8) and the prediction accuracy of the two eqns. for blasting vibration velocity was evaluated according to the correlation coefficient of the fitting curve. The analysis results are shown in Table 9.Table 9
Prediction model of attenuation laws of blasting vibration velocity.
DirectionEquation (8)Correlation coefficientEquation (11)Correlation coefficientXv=152.93Q3/r1.890.60v=237.46Q3/r2.13H/r−0.280.62Yv=63.43Q3/r1.670.63v=92.76Q3/r1.88H/r−0.240.65Zv=107.77Q3/r1.660.38v=340.36Q3/r2.30H/r−0.730.42Resultantv=114.43Q3/r1.590.43v=333.62Q3/r2.19H/r−0.680.47It can be concluded from Table9:(1)
According to the obtained prediction model of vibration velocity considering the explosion depth, using the control variable method, it can be found that with the increase in horizontal distance, the PPV decreases; with the increase in explosive quantity, the PPV increases; and with the increase in explosion depth, the PPV decreases; this is consistent with the results of the paper analysis and common sense. It is shown that the prediction model of vibration velocity of rock blasting in silt-rock strata established is effective in predicting PPV of ground.(2)
The correlation coefficient obtained by fitting 4.5 section blasting vibration velocity data by using the vibration velocity prediction model (equation (11)) considering explosion depth is greater than that of traditional Sadovski's formula (equation (8)), which indicates that the attenuation laws of ground blasting vibrations in silt-rock strata are more complex, and the mathematical prediction model established in this study can better reflect the effect of rock blasting vibrations on the ground in silt-rock strata.


## 4.1. Analysis of the Characteristics of Vibration Acceleration

The first blasting was carried out at 10 : 00 a.m. on March 31, 2016, and the second blasting was carried out at 3 : 00 p.m. on March 31, 2016, with an interval of 10 days, and the third blasting was carried out at 10 : 00 a.m. on April 10, 2016. The vibration acceleration test points were all 10 meters away from the explosion source, and the results of three tests of blasting are shown in Table3.Table 3
Results of three tests of blasting.
The total explosives (kg)The maximum charge of one section (kg)SectionsThe dominant frequency (Hz)The first blasting43.234.85,6,7X47.06Y55.56Z36.04The second blasting432.0150.05,6,7X10.34Y8.32Z14.60The third blasting331.2115.25,7,8X36.70Y37.38Z66.67The waveform duration of vibration accelerations of three blasting is shown in Figure7.Figure 7
Acceleration waveform of blasting vibrations, (a) acceleration waveform of the first blasting, (b) acceleration waveform of the second blasting, and (c) acceleration waveform of the third blasting.
(a)(b)(c)As shown in the acceleration waveform of blasting vibrations from Figure7, the peak particle acceleration of the three blasting vibrations has little difference to each other; the duration of acceleration waveform is about 0.4 s at 10 : 00 a.m. on March 31, 2016, in the first blasting, and about 0.8 s at 3 : 00 p.m. on March 31, 2016, in the second blasting, which is 0.4 s longer than the duration of the first blasting, with an interval of 10 days, and in the third blasting at 10 : 00 a.m. on April 10, 2016, the duration of acceleration waveform is about 0.5 s, which is close to the duration of the first blasting.The above data show that the silt has a good structural strength during the first blasting, but after the disturbance of the blasting load, the silt has been thixotropic liquefied [18]. During the second blasting with a short time interval, the explosive gas generates a cavity in the silt after the explosion. When the pressure of the explosive gas is less than the pressure of the silt outside the cavity, the silt around the cavity begins to compress into the cavity, and then, the cavity comes into being. The cavity is overcompressed, and then, the cavity expands again to work on the silt. In the process of pulsation, the water and stemming materials in the blast hole are continuously pushed out of the blast hole until they reach the orifice and escape when they come into contact with the atmosphere, which results in the extension of the waveform duration of acceleration in the second blasting. During the third blasting, the duration of acceleration waveform is close to that of the first blasting, which indicates that after a period of sedimentation, the silt has recovered a certain structural strength.


## 4.2. Analysis of the Distribution Characteristics of the Dominant Frequency of Blasting Vibration Velocity

Based on the analysis of the dominant frequency of vibration velocity of the three blasting tests in Table3, it is concluded that the dominant frequency of vibration velocity of the first blasting is between 30 and 60 Hz, which indicates that the silt has good structural strength; during the second blasting, the silt undergoes thixotropic liquefaction after the disturbance of the first blasting load. The dominant vibration frequency of dynamic velocity is reduced to 5–15 Hz, which is close to the natural vibration frequency of the building, and it is easy to cause the building resonance. At 10 : 00 am on April 10, 2016, the dominant frequency of the blasting vibration velocity is between 30 and 70 Hz during the third blasting, and the analysis shows that the structural strength of the silt has been restored due to 10 days from the last blasting. Therefore, the continuous blasting operation should be avoided in the field blasting. Each blasting operation should be separated for a period of time, and the next blasting will be beneficial to reduce the vibration damage until the silt recovers structural strength.


## 4.3. Energy Percentage of Each Frequency Bands of Blasting Vibrations

Taking a typical test data of blasting vibrations as an example, the field monitoring data of blasting vibrations are shown in Table4, and the time history curves of velocity of measured blasting vibration signals are shown in Figure 8.Table 4
Field monitoring data of blasting vibrations.
The total explosives (kg)The maximum charge of one section (kg)SectionsThe distance from explosion source (m)PPV (cm/s)XYZ4321505,6,71300.150.130.10Note. X was the horizontal direction pointing to the explosion source, Y was the horizontal direction perpendicular to X, and Z was the vertical direction.Figure 8
Time history curves of velocity of measured blasting vibrations.


### 4.3.1. Reconstruction Error

The signals of blasting vibrations belonged to random signals. At present, the db8 wavelet basis function was widely used to analyze the nonstationary vibration signal [19–21]. The error of wavelet reconstruction signal based on db8 is shown in Figure 9.Figure 9
Wavelet reconstruction and error based on db8. (a) Original blasting vibration signals and (b) error of reconstructed signals.
(a)(b)It could be seen from Figure9 that using db8 wavelet basis function to decompose the vibration signals in 8 layers, the relative error magnitude between the reconstructed signal and the original signal was 10−11, which met the analysis requirements.


### 4.3.2. Decomposition and Reconstruction

Using db8 wavelet basis function to decompose the signal of blasting vibrations with the scale of 8, the wavelet coefficients of 9 frequency bands were obtained [22]. During blasting vibration monitoring, the sampling frequency of the signal was set to 8000 Hz (2 s sampling 16000 points). According to the Shannon sampling theorem [21], the Nyquist frequency of sampling was 4000 Hz. According to the principle of wavelet decomposition [23], the bandwidth of nine wavelets was a8: 0–15.625 Hz, d8: 15.625–31.25 Hz, d7: 31.25–62.5 Hz, d6: 62.5–125 Hz, d5: 125–250 Hz, d4: 250–500 Hz, d3: 500–1000 Hz, d2: 1000–2000 Hz, and d1: 2000–4000 Hz. The blasting vibration components of different wavelet frequency bands are shown in Figure 10.Figure 10
Stratified reconstruction of signals by wavelet decomposition.


### 4.3.3. Percentage of Energy per Frequency Band

The statistics of energy percentage of each wavelet band of blasting vibration signals in different directions are shown in Table5.Table 5
Energy percentage of each frequency band in different directions (%).
Frequency bandsa8d8d7d6d5d4d3d2d1X43.809.6924.9421.120.380.010.010.020.04Y67.021.8111.2519.410.460.010.010.010.02Z88.914.934.131.800.150.030.010.020.03During blasting rock in silt-rock strata, the frequency bands of ground vibration energy are wide, but the frequency bands are mainly concentrated in 0–15.625 Hz, 31.25–62.5 Hz, and 62.5–125 Hz, which are accounting for more than 90% of the total signal energy. The dominant frequency bands of blasting vibration energy are in 0–15.625 Hz, which are close to the natural frequency of buildings and have a great impact on the vibration damage of surface buildings.The high-frequency signals above 125 Hz continue the whole monitoring period, which are represented by noise signal, with relatively small energy, and limited impact on surface buildings.


## 4.3.1. Reconstruction Error

The signals of blasting vibrations belonged to random signals. At present, the db8 wavelet basis function was widely used to analyze the nonstationary vibration signal [19–21]. The error of wavelet reconstruction signal based on db8 is shown in Figure 9.Figure 9
Wavelet reconstruction and error based on db8. (a) Original blasting vibration signals and (b) error of reconstructed signals.
(a)(b)It could be seen from Figure9 that using db8 wavelet basis function to decompose the vibration signals in 8 layers, the relative error magnitude between the reconstructed signal and the original signal was 10−11, which met the analysis requirements.


## 4.3.2. Decomposition and Reconstruction

Using db8 wavelet basis function to decompose the signal of blasting vibrations with the scale of 8, the wavelet coefficients of 9 frequency bands were obtained [22]. During blasting vibration monitoring, the sampling frequency of the signal was set to 8000 Hz (2 s sampling 16000 points). According to the Shannon sampling theorem [21], the Nyquist frequency of sampling was 4000 Hz. According to the principle of wavelet decomposition [23], the bandwidth of nine wavelets was a8: 0–15.625 Hz, d8: 15.625–31.25 Hz, d7: 31.25–62.5 Hz, d6: 62.5–125 Hz, d5: 125–250 Hz, d4: 250–500 Hz, d3: 500–1000 Hz, d2: 1000–2000 Hz, and d1: 2000–4000 Hz. The blasting vibration components of different wavelet frequency bands are shown in Figure 10.Figure 10
Stratified reconstruction of signals by wavelet decomposition.


## 4.3.3. Percentage of Energy per Frequency Band

The statistics of energy percentage of each wavelet band of blasting vibration signals in different directions are shown in Table5.Table 5
Energy percentage of each frequency band in different directions (%).
Frequency bandsa8d8d7d6d5d4d3d2d1X43.809.6924.9421.120.380.010.010.020.04Y67.021.8111.2519.410.460.010.010.010.02Z88.914.934.131.800.150.030.010.020.03During blasting rock in silt-rock strata, the frequency bands of ground vibration energy are wide, but the frequency bands are mainly concentrated in 0–15.625 Hz, 31.25–62.5 Hz, and 62.5–125 Hz, which are accounting for more than 90% of the total signal energy. The dominant frequency bands of blasting vibration energy are in 0–15.625 Hz, which are close to the natural frequency of buildings and have a great impact on the vibration damage of surface buildings.The high-frequency signals above 125 Hz continue the whole monitoring period, which are represented by noise signal, with relatively small energy, and limited impact on surface buildings.


## 4.4. The Energy of Each Frequency Bands of Blasting Vibrations

Taking the test data of the first blasting vibrations on the morning of March 31, 2016, as an example, the distance between the test point and the explosion source was 60 m, and the energy distributions of each wavelet band in different directions are shown in Figure11.Figure 11
Energy of each frequency bands in different directions (cm/s)2.It can be seen from Figure11 that the distribution characteristics of the frequency bands of blasting vibration energy in different directions of the ground are similar to each other, but the energy magnitudes are different from each other.The wavelet energy analysis was carried out on the six test data of blasting vibrations of rock blasting in silt-rock strata, and the directions of the maximum energy are shown in Table6.Table 6
Directions of the maximum energy.
The distance from explosion source (m)The first blastingThe second blastingThe third blastingThe forth blastingThe fifth blastingThe sixth blasting30XZZZXZ60ZXZZXZ90ZZZZXZ120YYXXZZ150The instrument was not triggeredZZZYZ180ZZXThe instrument was not triggeredZZNote. X was the horizontal direction pointing to the explosion source, Y was the horizontal direction perpendicular to X, and Z was the vertical direction.According to Table6, when the distance of the test points from the explosion source is in the range of 30–60 m, the number of times of maximum energy in Z direction is 8, and the number of times of maximum energy in X direction is 4.According to the analysis, because of the end effect of the cylindrical charge, there is air above the ground, without clamping effect, and there are more times of maximum energy in theZ direction; when the amount of cylindrical charge is less, the end effect is not obvious, and the explosion stress waves propagate outward in the form of cylinder, 13 blastholes in a row produce a cylinder wave, and the cylinder wave would overlap at the measuring point, so there are 4 times in the X direction.When the distance of the test point from the explosion source is in the range of 90–180 m, the number of times of maximum energy is 15 in theZ direction, 4 in the X direction, and 3 in the Y direction. It is considered that the maximum number of times of energy in the X direction and Y direction is the most due to the influence of pinch effect on explosion stress wave in the Z direction.The times of the maximum energy in theZ direction are the most, followed by the X direction when the distance of the test point from the explosion source is in the range of 30–60 m. Because the vibration energy in the Z direction is the largest and easy to cause up and down vibration damage to the building, the blasting vibration energy in the Z direction should be focused on analysis, followed by the X direction.The vibration energy in theZ direction is mainly analyzed when the distance of the test point from the explosion source is in the range of 90–180 m.


## 4.5. Attenuation Laws of Vibration Velocity of Rock Blasting in Silt-Rock Strata

To study the attenuation laws of vibration velocity of rock blasting in the silt-rock strata, four vibration velocity tests were carried out in the blasting area. Six measuring points of vibration velocity were arranged for each test. The test results of blasting vibration velocity are shown in Table7.Table 7
Testing results of blasting vibration velocity.
Test timesThe maximum charge of one sectionQ (kg)Vibration velocityNo. of measuring pointThe distance from explosion source (m)PPVυ (cm/s)TheX directionTheY directionTheZ directionResultantNo. 134.8ZC1301.671.321.931.94ZC2600.370.390.490.54ZC3900.180.170.330.34ZC41300.150.130.10.2ZC51700.040.060.090.1No. 2150ZC1304.683.28.18.58ZC2606.973.874.227.6ZC3900.460.470.971.01ZC41201.111.140.811.36ZC51500.340.390.620.63ZC61800.340.140.480.49No. 3165.6ZC1304.953.115.315.36ZC2602.371.758.058.25ZC3901.220.842.092.09ZC41200.730.360.830.87ZC51500.010.0170.00570.017No. 454ZC1302.130.680.630.63ZC2600.570.680.381.41ZC3901.60.820.51.61ZC41200.370.262.12.12ZC51500.220.181.021.03ZC61800.180.120.630.63


### 4.5.1. Prediction Model of Blasting Vibration Velocity Based on DA

Natural phenomena and engineering problems could be described by a series of physical quantities. To find the laws between these physical quantities, firstly, we could list the physical quantities that involved in these problems, ignore some physical quantities with less influence, and then classify the physical quantities involved according to their attributes; secondly, we could find out the interrelations between different physical quantities; and finally, we could find out the causal relationship between some physical quantities and the other physical quantities.There were relations among all kinds of physical quantities, which showed that their structure must be composed of several unified basic components. This basic component of physical quantity was called dimension. Because physics studied the evolution and movement of matter in space and time, all quantitative problems could not be separated from the three basic quantities of mass (M), time (T), and length (L). Therefore, M, T, and L were the most suitable dimensions for these three basic quantities. The dimensions of all other derived quantities could be expressed as the combination of these three basic quantities according to the definition or objective laws.The attenuation of seismic wave was affected by the factors of explosion source, geological conditions, distance of explosion source, and depth of explosion source when the blasting seismic wave propagated in silt-rock strata [24, 25]. There were 9 physical quantities involved in propagation velocity in silt-rock strata [26, 27], as shown in Table 8.Table 8
Physical quantities involved in blasting vibration velocity in silt-rock strata.
VariablesDimensionsThe dependent variablesParticle vibration displacementSLPPVυLT−1Peak particle accelerationαLT−2Frequency of particle vibrationfT−1The independent variablesQuality of explosivesQMHorizontal distance between the test point and the explosion sourcerLDepth of explosion sourceHLDensity of rockρML−3Propagation velocity of vibration wavecLT−1Detonation timetTNote. L, T, and M were dimensions of length, time, and mass, respectively.Based on the DA of Buckingham's theorem (π theorem), the PPV v of the ground during the rock blasting in the silt-rock strata could be expressed as follows:(1)v=ΦS,ρ,H,a,f,t,Q,r,c.According to theπ theorem, where independent dimensions were Q, r, and c, and π represents dimensionless quantities, then,(2)π=vQαrβcγ,π1=SQα1rβ1cγ1,π2=ρQα2rβ2cγ2,π3=HQα3rβ3cγ3,π4=aQα4rβ4cγ4,π5=fQα5rβ5cγ5,π6=tQα6rβ6cγ6.where α, β, and γ were undetermined coefficients. According to the dimensional homogeneous theorem, there were as follows:(3)π=vc,π1=Sr,π2=ρQr−3,π3=Hr,π4=ar−1c2,π5=fr−1c,π6=trc−1.Substituting (3) into (1) could obtain the following:(4)vc=ΦSr,ρQr−3,Hr,ar−1c2,fr−1c,trc−1.Since the product and power of different dimensionless numbersπ were still dimensionless, π2 and π3 are combined as follows to get a new dimensionless number π7:(5)π7=π21/3β1π3β2=ρ3rQ3β1Hrβ2,where β1 and β2 were the indexes of π2 and π3, respectively.At the same site,ρ and c could be approximately constant. Therefore, from (5), it could be considered that v and Q3/rβ1H/rβ2 had a functional relationship.To sum up, this functional relationship could be expressed as follows:(6)lnv=α1+β1lnQ3r+α2+β2lnHr.Iflnv0=α1+β1lnQ3/R., there were as follows:(7)lnv0=α1+β1lnQ3−β1lnr.In equation (6), α1 and α2 were the coefficients given in the process of function transformation; −β1ln r was the attenuation of blasting vibration velocity with distance r; β1 was the attenuation index, which mainly reflected the effect of site medium conditions; and α1 + (β1lnQ)/3 comprehensively reflected the contribution of site medium conditions and explosive quantity to the PPV of ground.Iflnk1=lnα1, there were as follows:(8)v0=k1Q3rβ1.Equation (8) was Sadovski's equation. without the effect of the depth of explosion source. Substituting (8) into (6), we could get the following:(9)lnv=lnv0+α2+β2lnHr.Iflnk1=lnα1, then (9) could be changed into the following:(10)v=k1k2Q3rβ1Hrβ2.Ifk=k1k2, the mathematical model of the attenuation laws of blasting vibration velocity in silt-rock strata considering the effect of burial depth was established as follows:(11)v=kQ3rβ1Hrβ2,where k was the site effect coefficient, β1 was the attenuation coefficient of blasting vibrations, and β2 was the effect coefficient of explosion source depth.


### 4.5.2. Comparative Analysis of Prediction Results of Blasting Vibration Velocity

The explosion source was below the ground surface during blasting rock in silt-rock strata, while measuring points of blasting vibration velocity were often located on the ground surface, which was also the key content of the analysis of attenuation laws of blasting vibrations. Therefore, this study selected the ground data for analysis and established the prediction model of vibration velocity of rock blasting in silt-rock strata. According to Z. C. Zhang’s research [28], the attenuation laws of vibration velocity near blasting zone were different from the attenuation laws of vibration velocity far blasting zone. R¯=r/Q1/3 was the scaled distance, r was the horizontal distance between the test point and the explosion source, and Q was the quality of explosives. The R¯<5.0 was in the near blasting zone. The R¯ of monitoring data in this paper is greater than 5.0, and the monitoring data were not within the range of near blast area. The monitoring data could be used to analyze the attenuation law. Equation (11) in Section 4.5.1 is used for regression fitting analysis of the blasting vibration data listed in Table 7 to obtain the prediction model of vibration velocity of rock blasting in the silt-rock strata.In regression analysis, the correlation coefficient was an indicator to measure the linear correlation between variables. The larger the correlation coefficient was, the smaller the error was, the higher the linear correlation between variables was, the higher the reliability of the prediction model was, and vice versa. To determine the accuracy of the established mathematical model, the established mathematical model (11) was compared with the traditional Sadovski's formula (8) and the prediction accuracy of the two eqns. for blasting vibration velocity was evaluated according to the correlation coefficient of the fitting curve. The analysis results are shown in Table 9.Table 9
Prediction model of attenuation laws of blasting vibration velocity.
DirectionEquation (8)Correlation coefficientEquation (11)Correlation coefficientXv=152.93Q3/r1.890.60v=237.46Q3/r2.13H/r−0.280.62Yv=63.43Q3/r1.670.63v=92.76Q3/r1.88H/r−0.240.65Zv=107.77Q3/r1.660.38v=340.36Q3/r2.30H/r−0.730.42Resultantv=114.43Q3/r1.590.43v=333.62Q3/r2.19H/r−0.680.47It can be concluded from Table9:(1)
According to the obtained prediction model of vibration velocity considering the explosion depth, using the control variable method, it can be found that with the increase in horizontal distance, the PPV decreases; with the increase in explosive quantity, the PPV increases; and with the increase in explosion depth, the PPV decreases; this is consistent with the results of the paper analysis and common sense. It is shown that the prediction model of vibration velocity of rock blasting in silt-rock strata established is effective in predicting PPV of ground.(2)
The correlation coefficient obtained by fitting 4.5 section blasting vibration velocity data by using the vibration velocity prediction model (equation (11)) considering explosion depth is greater than that of traditional Sadovski's formula (equation (8)), which indicates that the attenuation laws of ground blasting vibrations in silt-rock strata are more complex, and the mathematical prediction model established in this study can better reflect the effect of rock blasting vibrations on the ground in silt-rock strata.


## 4.5.1. Prediction Model of Blasting Vibration Velocity Based on DA

Natural phenomena and engineering problems could be described by a series of physical quantities. To find the laws between these physical quantities, firstly, we could list the physical quantities that involved in these problems, ignore some physical quantities with less influence, and then classify the physical quantities involved according to their attributes; secondly, we could find out the interrelations between different physical quantities; and finally, we could find out the causal relationship between some physical quantities and the other physical quantities.There were relations among all kinds of physical quantities, which showed that their structure must be composed of several unified basic components. This basic component of physical quantity was called dimension. Because physics studied the evolution and movement of matter in space and time, all quantitative problems could not be separated from the three basic quantities of mass (M), time (T), and length (L). Therefore, M, T, and L were the most suitable dimensions for these three basic quantities. The dimensions of all other derived quantities could be expressed as the combination of these three basic quantities according to the definition or objective laws.The attenuation of seismic wave was affected by the factors of explosion source, geological conditions, distance of explosion source, and depth of explosion source when the blasting seismic wave propagated in silt-rock strata [24, 25]. There were 9 physical quantities involved in propagation velocity in silt-rock strata [26, 27], as shown in Table 8.Table 8
Physical quantities involved in blasting vibration velocity in silt-rock strata.
VariablesDimensionsThe dependent variablesParticle vibration displacementSLPPVυLT−1Peak particle accelerationαLT−2Frequency of particle vibrationfT−1The independent variablesQuality of explosivesQMHorizontal distance between the test point and the explosion sourcerLDepth of explosion sourceHLDensity of rockρML−3Propagation velocity of vibration wavecLT−1Detonation timetTNote. L, T, and M were dimensions of length, time, and mass, respectively.Based on the DA of Buckingham's theorem (π theorem), the PPV v of the ground during the rock blasting in the silt-rock strata could be expressed as follows:(1)v=ΦS,ρ,H,a,f,t,Q,r,c.According to theπ theorem, where independent dimensions were Q, r, and c, and π represents dimensionless quantities, then,(2)π=vQαrβcγ,π1=SQα1rβ1cγ1,π2=ρQα2rβ2cγ2,π3=HQα3rβ3cγ3,π4=aQα4rβ4cγ4,π5=fQα5rβ5cγ5,π6=tQα6rβ6cγ6.where α, β, and γ were undetermined coefficients. According to the dimensional homogeneous theorem, there were as follows:(3)π=vc,π1=Sr,π2=ρQr−3,π3=Hr,π4=ar−1c2,π5=fr−1c,π6=trc−1.Substituting (3) into (1) could obtain the following:(4)vc=ΦSr,ρQr−3,Hr,ar−1c2,fr−1c,trc−1.Since the product and power of different dimensionless numbersπ were still dimensionless, π2 and π3 are combined as follows to get a new dimensionless number π7:(5)π7=π21/3β1π3β2=ρ3rQ3β1Hrβ2,where β1 and β2 were the indexes of π2 and π3, respectively.At the same site,ρ and c could be approximately constant. Therefore, from (5), it could be considered that v and Q3/rβ1H/rβ2 had a functional relationship.To sum up, this functional relationship could be expressed as follows:(6)lnv=α1+β1lnQ3r+α2+β2lnHr.Iflnv0=α1+β1lnQ3/R., there were as follows:(7)lnv0=α1+β1lnQ3−β1lnr.In equation (6), α1 and α2 were the coefficients given in the process of function transformation; −β1ln r was the attenuation of blasting vibration velocity with distance r; β1 was the attenuation index, which mainly reflected the effect of site medium conditions; and α1 + (β1lnQ)/3 comprehensively reflected the contribution of site medium conditions and explosive quantity to the PPV of ground.Iflnk1=lnα1, there were as follows:(8)v0=k1Q3rβ1.Equation (8) was Sadovski's equation. without the effect of the depth of explosion source. Substituting (8) into (6), we could get the following:(9)lnv=lnv0+α2+β2lnHr.Iflnk1=lnα1, then (9) could be changed into the following:(10)v=k1k2Q3rβ1Hrβ2.Ifk=k1k2, the mathematical model of the attenuation laws of blasting vibration velocity in silt-rock strata considering the effect of burial depth was established as follows:(11)v=kQ3rβ1Hrβ2,where k was the site effect coefficient, β1 was the attenuation coefficient of blasting vibrations, and β2 was the effect coefficient of explosion source depth.


## 4.5.2. Comparative Analysis of Prediction Results of Blasting Vibration Velocity

The explosion source was below the ground surface during blasting rock in silt-rock strata, while measuring points of blasting vibration velocity were often located on the ground surface, which was also the key content of the analysis of attenuation laws of blasting vibrations. Therefore, this study selected the ground data for analysis and established the prediction model of vibration velocity of rock blasting in silt-rock strata. According to Z. C. Zhang’s research [28], the attenuation laws of vibration velocity near blasting zone were different from the attenuation laws of vibration velocity far blasting zone. R¯=r/Q1/3 was the scaled distance, r was the horizontal distance between the test point and the explosion source, and Q was the quality of explosives. The R¯<5.0 was in the near blasting zone. The R¯ of monitoring data in this paper is greater than 5.0, and the monitoring data were not within the range of near blast area. The monitoring data could be used to analyze the attenuation law. Equation (11) in Section 4.5.1 is used for regression fitting analysis of the blasting vibration data listed in Table 7 to obtain the prediction model of vibration velocity of rock blasting in the silt-rock strata.In regression analysis, the correlation coefficient was an indicator to measure the linear correlation between variables. The larger the correlation coefficient was, the smaller the error was, the higher the linear correlation between variables was, the higher the reliability of the prediction model was, and vice versa. To determine the accuracy of the established mathematical model, the established mathematical model (11) was compared with the traditional Sadovski's formula (8) and the prediction accuracy of the two eqns. for blasting vibration velocity was evaluated according to the correlation coefficient of the fitting curve. The analysis results are shown in Table 9.Table 9
Prediction model of attenuation laws of blasting vibration velocity.
DirectionEquation (8)Correlation coefficientEquation (11)Correlation coefficientXv=152.93Q3/r1.890.60v=237.46Q3/r2.13H/r−0.280.62Yv=63.43Q3/r1.670.63v=92.76Q3/r1.88H/r−0.240.65Zv=107.77Q3/r1.660.38v=340.36Q3/r2.30H/r−0.730.42Resultantv=114.43Q3/r1.590.43v=333.62Q3/r2.19H/r−0.680.47It can be concluded from Table9:(1)
According to the obtained prediction model of vibration velocity considering the explosion depth, using the control variable method, it can be found that with the increase in horizontal distance, the PPV decreases; with the increase in explosive quantity, the PPV increases; and with the increase in explosion depth, the PPV decreases; this is consistent with the results of the paper analysis and common sense. It is shown that the prediction model of vibration velocity of rock blasting in silt-rock strata established is effective in predicting PPV of ground.(2)
The correlation coefficient obtained by fitting 4.5 section blasting vibration velocity data by using the vibration velocity prediction model (equation (11)) considering explosion depth is greater than that of traditional Sadovski's formula (equation (8)), which indicates that the attenuation laws of ground blasting vibrations in silt-rock strata are more complex, and the mathematical prediction model established in this study can better reflect the effect of rock blasting vibrations on the ground in silt-rock strata.


## 5. Conclusions

Based on the analysis of vibration characteristics of rock blasting in silt-rock strata, the following conclusions are drawn:(1)
In the first blasting, the silt has good structural strength, and in the second blasting with an interval of 5 hours, the silt is subject to thixotropic liquefaction after the disturbance of blasting load. When the third blasting is carried out with an interval of 10 days, the waveform duration of acceleration is close to that of the first blasting waveform, and the dominant frequency of blasting vibration velocity is close to that of the first blasting too. It is considered that after a period of sedimentation, the silt recovered to a certain structural strength, and the residual oscillation of pulsation is weakened.(2)
The distribution of frequency bands of vibration energy is wide, but they are mainly concentrated in the range of 0–15.625 Hz, 31.25–62.5 Hz, and 62.5–125 Hz, and they are accounting for more than 90% of the total energy of signals during blasting rock. The dominant frequency bands of blasting vibration energy are in the range of 0–15.625 Hz.(3)
The distribution characteristics of frequency bands of blasting vibration energy in different directions are similar to each other during blasting rock, but the energy magnitudes are different from each other. The times of maximum energy inZ direction are the most, followed by X direction when the distance from the explosion source is in the range of 30–60 m. The blasting vibration energy in Z direction should be selectively analyzed, followed by X direction. The blasting vibration energy in Z direction should be mainly analyzed when the distance from the explosion source is in the range of 90–180 m.(4)
The correlation coefficients of the prediction model based on the DA are greater than that of the traditional Sadovski's formula through fitting the measured blasting vibration data.



---

*Source: 3318965-2021-11-30.xml*