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

School of Civil Engineering and Architecture, Anhui University of Science and Technology, No. 168 Taifeng Street, Huainan 232001, Anhui, China Faculty of Engineering, China University of Geosciences, No. 388 Lumo Road, Wuhan 430074, Hubei, China Zhuhai Baopoxin Research & Development Co., Ltd, No. 2099 Fenghuang North Road, Zhuhai 519099, Guangdong, China School of Chemical Engineering, Anhui University of Science and Technology, No. 168 Taifeng Street, Huainan 232001, Anhui, China


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][2][3][4][5]. ey 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.
ere are many studies about blasting in rock, soil, air, and water, but few researches about rock blasting in the siltrock strata. e 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 blastinduced 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. e analysis of vibration characteristics of ground of rock blasting in silt-rock strata by wavelet analysis and DA has been less reported at present. is 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. e 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. e research results could provide a reference for the safe blasting distance of buildings under similar geological conditions.

Project Overview
e 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 Figure 1. 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.
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.
e diagram of charging on field is shown in Figure 2. e rock cores before blasting and after blasting are shown in Figure 3. e blasting area is located in the artificial reclamation area, and there are three strata from the ground surface to the subway tunnel. (1) e 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) e 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. e layer thickness is in the range of 0.50-30.30 m. (3) e third stratum is rock. It is weakly weathered granite. It is brown-yellow, gray-brown, graywhite, and other colors, medium-coarse grain structure, and massive structure. e joint fissures are developed relatively. e rock quality is hard. e buried depth of the layer fluctuates greatly, and the distribution is discontinuous. e physical and mechanical parameters of artificial fill, silt, and rock are listed in Tables 1 and 2, respectively. e length of the blasting area is 960 cm, and the width of the blasting area is 160 cm. e 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. e 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. ree rows of blastholes are detonated each time, 13 blastholes are arranged in a row, and the millisecond blasting network in the blasthole is adopted. e millisecond sections of three rows of blastholes are 6, 7, and 8 in turn, and the delay time is 50 ms. e layouts of blastholes are shown in Figure 4.
After drilling the blastholes, the blastholes are filled with underground water immediately, and the blastholes are water coupling charge. e diameter of the explosive roll is 6 cm, and the density of the explosive is 1 061 kg/m 3 . e 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. e blastholes are stemmed with gravels of 1-1.5 cm diameter, and the stemming length is 1 000 cm.

Test Points of Blasting Vibrations
e 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. e arrangement of blasting vibration meters is shown in Figure 5. ree acceleration measuring points and six velocity measuring points were arranged in total. e measuring points were arranged along the ray with the explosion source

Analysis of the Characteristics of Vibration
Acceleration. e 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. e vibration acceleration test points were all 10 meters away from the explosion source, and the results of three tests of blasting are shown in Table 3. e waveform duration of vibration accelerations of three blasting is shown in Figure 7.
As shown in the acceleration waveform of blasting vibrations from Figure 7, 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.
e 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. e 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.

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 Table 3, 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. e 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. erefore, 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.

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 Table 4, and the time history curves of velocity of measured blasting vibration signals are shown in Figure 8.

Reconstruction Error.
e 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][20][21]. e error of wavelet reconstruction signal based on db8 is shown in Figure 9.
It could be seen from Figure 9 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.

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:

Percentage of Energy per Frequency Band.
e statistics of energy percentage of each wavelet band of blasting vibration signals in different directions are shown in Table 5.
During 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. e 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. e 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.
e 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 Figure 11.
It can be seen from Figure 11 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. e 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 Table 6.
According to Table 6, 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 the Z 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 the Z 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. e times of the maximum energy in the Z 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. e vibration energy in the Z direction is mainly analyzed when the distance of the test point from the explosion source is in the range of 90-180 m.

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 Shock and Vibration 7 measuring points of vibration velocity were arranged for each test. e test results of blasting vibration velocity are shown in Table 7.

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. ere were relations among all kinds of physical quantities, which showed that their structure must be composed of several unified basic components. is 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). erefore, M, T, and L were the most suitable dimensions for these three basic quantities. e dimensions of all other derived quantities could be expressed as the combination of these three basic quantities according to the definition or objective laws. e 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].
ere were 9 physical quantities involved in propagation velocity in silt-rock strata [26,27], as shown in Table 8.
Based on the DA of Buckingham's theorem (π theorem), the PPV v of the ground during the rock blasting in the siltrock strata could be expressed as follows: v � Φ (S, ρ, H, a, f, t, Q, r, c). (1)

Shock and Vibration
According to the π theorem, where independent dimensions were Q, r, and c, and π represents dimensionless quantities, then, , 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 : Note. X was the horizontal direction pointing to the explosion source, Y was the horizontal direction perpendicular to X, and Z was the vertical direction.
10 Shock and Vibration where β 1 and β 2 were the indexes of π 2 and π 3 , respectively. At the same site, ρ and c could be approximately constant. erefore, from (5), it could be considered that v and To sum up, this functional relationship could be expressed as follows: If ln v 0 � α 1 + β 1 ln( �� Q 3 √ /R)., there were as follows: In equation (6), α 1 and α 2 were the coefficients given in the process of function transformation; −β 1 ln 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 + (β 1 lnQ)/3 comprehensively reflected the contribution of site medium conditions and explosive quantity to the PPV of ground.
If ln k 1 � ln α 1 , there were as follows: Equation (8) was Sadovski's equation. without the effect of the depth of explosion source. Substituting (8) into (6), we could get the following: If ln k 1 � ln α 1 , then (9) could be changed into the following: If k � k 1 k 2 , 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: 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.

Comparative Analysis of Prediction Results of Blasting
Vibration Velocity. e 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. erefore, 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/Q 1/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. e R < 5.0 was in the near blasting zone. e 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. e 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. e 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.
e analysis results are shown in Table 9. It can be concluded from Table 9: (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 siltrock strata established is effective in predicting PPV of ground. (2) e 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.

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) e 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 e dominant frequency bands of blasting vibration energy are in the range of 0-15.625 Hz.
(3) e 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. e times of maximum energy in Z direction are the most, followed by X direction when the distance from the explosion source is in the range of 30-60 m. e blasting vibration energy in Z direction should be selectively analyzed, followed by X direction. e 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) e 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.

Data Availability
e processed data required to reproduce these findings cannot be shared at this time as the data also form part of an ongoing study.

Conflicts of Interest
e authors declare that they have no conflicts of interest.