Vibration Signal Analysis ofWater Seal Blasting Based onWavelet Threshold Denoising and HHT Transformation

,e blasting vibration signal obtained from tunnel construction monitoring is affected by the external environment, which contains a lot of noise that causes distortion during signal processing. To analyse the blasting vibration signal and determine the appropriate water seal blasting charge structure for construction, combined with wavelet threshold denoising method and HHT transformation, the blasting vibration signals of the four charge structures of conventional charge, water interval charge at both ends, water interval charge at the orifice, and water interval charge at the hole bottom are denoised and HHT is analysed. ,e results show that the wavelet threshold method can effectively eliminate high-frequency noise in the blasting vibration signals and retain information carried by the vibration signal itself. ,e frequency and energy of the blasting vibration signals of the water interval charge at both ends are densely distributed in the range of 0 s to 0.9 s and below 100Hz. ,e frequency and energy of the blasting signals of the other three charging structures are reduced within the same range, sparse areas appear, and the instantaneous total energy is smaller than that with a water interval charge at both ends, which shows that the water interval charge at both ends can effectively apply explosive energy to the surrounding rock and reduce energy loss in the explosive. ,e blasting vibration signal energy of the water interval charge at both ends is mainly concentrated in components IMF2 to IMF5, and the corresponding frequencies are concentrated at 6Hz to 11Hz and 20Hz to 70Hz, while the blasting vibration signal energy of other three charge structures is mainly distributed in components IMF2 to IMF4, corresponding frequencies are concentrated within 20Hz to 70Hz, and the distribution at low frequencies is not obvious. ,erefore, when using the water interval charge at both ends, it is necessary to increase the main vibration frequency of the original vibration signals by reducing the single section charge and using frequency shift technology to avoid the natural frequency of the structure and reduce resonance-induced damage.


Introduction
In the tunnel blasting process, in addition to crushing the rock, the explosive energy propagates within the surrounding rock to form a vibration signal in the form of stress waves. Blasting vibration signals not only reflect the dynamic mechanical characteristics of the surrounding rock but also contain much useful feature information [1][2][3]. erefore, the monitoring and analysis of blasting vibration velocity can provide a reference for the safe construction of tunnels, which is of great importance [4][5][6].
ere are many ways to analyse vibration signals: Liu et al. used MATLAB ™ combined with wavelet analysis to analyse the energy of blasting vibration signals monitored in practice. According to the frequency domain energy of the main vibration frequency of the blasting vibration signals and the frequency domain energy of the controlled structure self-vibration frequency, a new method of blasting hazard evaluation was proposed [7]. Cao et al. used the HHT method to analyse the measured blasting vibration signals from the perspective of their time-frequency and energy distribution characteristics and then assessed the safety of a slope and surrounding buildings [8]. Wang et al. used the HHT method to obtain the Hilbert spectrum, marginal spectrum, and instantaneous energy diagram of blasting vibration signals and studied the effect of blasting vibration on nearby tunnels [9]. Ning et al. use the wavelet packet analysis method to decompose vibration signals collected in the field and assess the attenuation of vibration waves after stope blasting [10]. Li et al. used multifractal analysis to analyse the singular spectrum of vibration signals and then used a multifractal spectrum derived by use of the wavelet transform [11]. Among them, the HHTmethod is a relatively new time-frequency analysis method, which has the advantages of self-adaptability, no need to fix an a priori basis, and suitability for the analysis of nonstationary signals. e blasting vibration signals monitored in the field often contain noise, which affects the effectiveness of signal analysis. erefore, people use various means to denoise such signals. e common denoising methods of blasting vibration include wavelet technology (wavelet threshold denoising, modulus maximum wavelet denoising, wavelet packet threshold denoising, translation invariant wavelet denoising, etc.), empirical mode decomposition technology (EMD denoising, EEMD denoising, etc.), and combinations of the two methods (EMD-wavelet threshold denoising). Yang et al. used EMD-FSWT and EEMD-SPWV methods to analyse blasting vibration signals and discussed the noise reduction effect of EMD-FSWT and the multifractal characteristics of blasting vibration signals in the scale-free zone [12,13]. Shao adopts a low-pass denoising method based on CEEMD (Complementary Ensemble Empirical Mode Decomposition), which can not only eliminate the irregular noise in the tunnel blasting vibration signals but also preserve the authenticity and integrity of the signals [14]. To overcome the fact that EEMD denoising cannot deal with IMFs, Wei et al. combine the "3σ" rule and Singular Value Decomposition (SVD) to improve EEMD denoising [15]. According to the characteristics of decomposition and reconstruction of the EMD method and the flexibility and debuggability of wavelet threshold methods, Rao et al. proposed the EMD-wavelet threshold denoising method [16]. Among them, wavelet threshold denoising is widely used because of its fast calculation and reasonable denoising capabilities.
For environmental protection, water seal blasting has attracted much attention because of its good blasting effect and reducing blasting dust emissions; however, there are few reports on the differences in vibration signals between the charge structure of water-sealed blasting and conventional blasting. erefore, we take the Wenquan Tunnel in Beijing Yanchong Expressway as the engineering background and combine the wavelet threshold denoising and HHT methods to compare and analyse the measured vibration signals of water seal blasting and conventional blasting signals. According to the analytical results, the appropriate water seal blasting charge structure is selected for tunnel construction.

The Principle of Wavelet Threshold
Denoising and HHT 2.1. Wavelet reshold Denoising Principle. Donoho and Johnstone proposed the wavelet threshold shrinkage denoising method in 1994: therein, they mentioned the concept of wavelet threshold denoising [17]. Wavelet threshold denoising is a method that maps the signals to the wavelet domain. According to the different properties and mechanisms of the wavelet coefficients of noise and noise on different scales, the wavelet coefficients of noisy signals are processed. e blasting vibration signals are decomposed by wavelet transformation to generate low-frequency wavelet coefficients and high-frequency wavelet coefficients. e real signal is a low-frequency signal and the noise is a highfrequency signal. When increasing the number of decomposition layers, the correlation coefficients between lowfrequency wavelet coefficients and real signals will change from low to high and then to low. When the correlation coefficient reaches its highest value, it can be set as a threshold value, the wavelet coefficients containing the noise are filtered and inverted, and then the denoised signals can be obtained (Figure 1). Wavelet threshold denoising is one of the most widely used methods of wavelet denoising because of its advantages such as approximate optimal estimation of the original signal, fast computing speed, and wide adaptability [16,[18][19][20]. In addition, wavelet threshold denoising can almost completely suppress the noise, and the peak points reflecting the original signal characteristics can be well preserved.

Selection of reshold.
In the process of threshold selection, choosing the appropriate threshold has an important influence on the effect of noise reduction. If the threshold value is too small, the signal noise reduction is not complete; if the threshold value is too large, the useful signal will be deleted and the accuracy will be affected. e main methods for selecting the threshold are discussed in the following.

General
reshold Method (Sqtwolog). General threshold method is the conclusion and the best threshold for the joint distribution of multidimensional independent normal variables when the dimensionality tends to infinity and the minimum maximum estimation limit. When this method is used, the same threshold value is applied to wavelet decomposition to get all levels of wavelet coefficients, and the estimation results are relatively smooth. e formula for calculating the threshold of the general threshold method is as follows: where N is the sampling length of the signal and σ is the variance of the noise signal. σ � median (W HH )/0.6475, W HH are the orthogonal wavelet coefficients of noise in the highfrequency subband, and median denotes taking the median average.

Stein Unbiased Risk reshold Method (Rigrsure).
Stein unbiased risk threshold method is an adaptive threshold selection based on the principle of Stein unbiased likelihood estimation. Given a threshold value to its likelihood estimation, the selected threshold can be obtained by minimizing the nonlikelihood threshold value. Assuming P � [p 1 , p 2 , . . ., p n ] and p 1 < p 2 < · · · < p n , the elements contained in P are the squares of the wavelet coefficients 2 Advances in Civil Engineering arranged in ascending order. If R is defined as a risk trace, the element r i therein is When the minimum value r min of the element in R is selected as the risk value, the threshold calculation formula is as follows:

Heuristic reshold Method (Heursure)
. e heuristic threshold method selects the optimal prediction variable threshold. In practical processing, when the signal-to-noise ratio is small, the error of the Stein unbiased risk threshold method in calculating the threshold is not conducive to accurate signal denoising analysis. e heuristic threshold method combines the advantages of the general threshold and Stein unbiased risk threshold. When the signal-to-noise ratio is low, the general threshold method is used to calculate the threshold. When the signal-to-noise ratio is high, the Stein unbiased risk threshold method is used to calculate the threshold.

Minimax reshold Method (Minimax).
e minimax threshold method uses the minimax principle to generate the threshold value and uses the minimum mean square error as the objective function to generate an extreme value. is threshold estimator minimizes the maximum mean square error in a given function. e principle of the minimax threshold method is to minimise the maximum risk of prediction. e calculation method of threshold is as follows: Because the threshold selection rules of rigrsure and minimax are conservative, only a part of the coefficients are set to zero. erefore, these two thresholds are very useful to extract small and weak signals when a small part of the highfrequency information is in the noise range. However, the selection rules of heursure threshold and sqtwolog threshold are more effective in removing noise, but the high-frequency part of useful signal may be removed as noise signal. e above four thresholds are used for noise reduction of noisy signals, and the standard deviation is calculated [18]. e calculation results are shown in Table 1. It can be seen from Table 1 that using the sqtwolog threshold noise reduction can get better results.

Selection of reshold Function.
In addition to the aforementioned threshold selection, the threshold function is an important factor affecting wavelet denoising; at present, there are two kinds of commonly used threshold functions: hard threshold functions and soft threshold functions. e two functions are expressed as follows: (1) Hard threshold function (2) Soft threshold function where d j,k are the original wavelet coefficients; d ∧ j,k are the estimated wavelet coefficients after calculating the threshold function; T is the threshold. Figure 2 shows the difference between the two threshold functions.
By comparing mean square error of the two functions, the hard threshold function is better than the soft threshold function; however, due to the influence of discontinuity in the function, the estimated signal from the hard threshold function will produce oscillations, while the wavelet coefficients from the soft threshold function have better overall continuity, so that the estimated signal will not be accompanied by additional oscillations, and the signal, after noise reduction, will be smoother [18]. erefore, soft threshold function is selected for analysis.

Basic Principles of HHT.
Huang et al. proposed the HHT method at the end of the last century [21]. e principle thereof is to change the unstable random signal into a stationary signal. First, empirical mode analysis (EMD) is

Empirical Modal Analysis and Hilbert Transform.
To study transient and nonstationary phenomena, the frequency must be a function related to time; however, no consensus as to the definition of instantaneous frequency has been reached. To obtain an effective and representative instantaneous frequency, Huang and others propose the necessary conditions for satisfying a meaningful instantaneous frequency (IF) in terms of physical characteristics: functions must be symmetric with respect to local zero mean and have the same extremum and zero-crossing points. en the intrinsic mode function (IMF) is defined as the function satisfying the following two conditions that, when using EMD to analyse vibration signals, must be satisfied: (1) In the whole data sequence, the number of extreme points and the number of zero-crossing points must be equal or differ at most by one. (2) At any point on the signal, the mean value of the envelope determined by the local maximum point and the envelope determined by the local minimum point are both zero. at is, the signal is locally symmetric with respect to the time axis.
e Hilbert transform is applied to multiple IMF signals obtained by the EMD method. e instantaneous spectrum of each IMF component can be obtained.
e Hilbert spectrum can be obtained by synthesising the instantaneous spectrum of all IMF components. e Hilbert transform is a type of linear transformation, which emphasises the local properties of the signal. e instantaneous frequency obtained by use of the Hilbert transform is the best definition: it avoids many spurious high-and low-frequency components produced by the Fourier transform.

Engineering Context.
e Beijing Yanchong Expressway (Beijing section) project is a key supporting project jointly organised by Beijing and Zhangjiakou City for the 2022 Winter Olympics. In this test, the Wenquan Tunnel on the Yanchong Expressway (Beijing section) is taken as the engineering background against which to conduct conventional blasting and water seal blasting tests with different charging structures. e excavation of the hot spring tunnel involves drilling through rock (both hard rock and associated secondary rock). e lithology of the tunnel section is weak, the joint fissures are developed, the rock is broken, and the stability thereof is poor. According to the comprehensive consideration of tunnel section size, surrounding rock stability and construction continuity, the millisecond delay blasting method is adopted for upper and lower steps. e  blast-hole on the excavation section of the upper bench is detonated according to the partition arrangement and sequence of the cutting hole, the auxiliary hole, and the peripheral hole and gradually expands to complete a blasting excavation. Horizontal hole slotting blasting is adopted in the lower bench to reduce damage to the lining and the surface of the upper bench.

Test Method
is test is based on the upper bench blasting, and its cyclical footage is controlled to 3 m. e length of auxiliary hole and peripheral hole is 3.3 m, the length of the cut hole is 3.5 m, and the diameter of the borehole is 42 mm. e arrangement of the blast-hole and initiation sequence are shown in Figure 3; relevant parameters are listed in Table 2.

Charge Structure.
To highlight the characteristics of water-sealed blasting, water-sealed blasting is compared with conventional blasting. In this test, four charging structures are used: (a) conventional charge structure, (b) water interval charge at both ends, (c) water interval charge at the orifice, and (d) water interval charge at the hole bottom. e charge structure is shown in Figure 4.
To reduce the influence of external conditions, the quantity of charge in the same type of blast-hole is the same in the four charge structures. e stemming length of the blast-hole is 0.9 m. e water length is 0.6 m (the water length of two sections in the water interval charge at both ends is 0.3 m), and the air-gap length is 0.6 m.

Blasting Vibration Monitoring.
To monitor the vibration signals of blasting effectively, a series of tests was conducted for conventional blasting and water-sealed blasting with different charging structures. According to the actual situation of the site and under the premise of ensuring safety, vibration monitoring point is arranged 50 m from the centreline of the upper step of the tunnel to the tunnel face ( Figure 5).

Monitoring Instrument.
TC-4850 blasting signal recorder is selected as the monitoring instrument, as shown in Figure 6. e recorder is equipped with embedded computer module and LCD, which can directly debug parameters on the instrument and preview the information of vibration speed peak value, frequency, and waveform immediately after sampling. e range of the explosion vibration recorder is set to 0.001 cm/s∼35.4 cm/s, the sampling frequency is 8000 Hz, the sampling time is 5 s, and the triggering level is 0.2 V. e instrument has 128 M storage space, which can continuously trigger the recording and storage. Blasting Safety Regulation (GB6722-2014)" [22], the vibration speed can be used as the standard for blasting vibration damage and protection. By comparing and analysing the in situ monitoring data, the results show that the vertical vibration velocity Z (vertical pavement) is greater than the lateral X (parallel tunnel face) and the longitudinal Y (vertical tunnel face) velocities. erefore, the vertical vibration velocity is selected as the criterion for comparative analysis.

Signal Denoising and Analysis
In the process of field blasting, the blasting vibration signals contain a large amount of blasting information. Time-frequency analysis is an important object of blasting research; however, the environment is often complicated and subject to the effects of interference from reflection and refraction of blasting signals in the propagating medium, external magnetic fields, and so forth, which cause the monitored signals to contain noise that affects the accuracy of signal analysis. It is necessary to denoise the signal to improve the accuracy of signal analysis. In this paper, the general threshold (sqtwolog) combined with wavelet threshold denoising method of soft threshold function is used to process blasting vibration signals of four charge structures [23,24]. e effect of noise reduction is shown in Figure 7. Figure 7 shows that there are few differences between the original signals and the denoised signals after the blasting signals are denoised by the wavelet threshold method. It is shown that the wavelet threshold method can effectively preserve the information carried by the blasting signals while denoising the signals. erefore, we take the vibration signals after noise reduction as the analysis object and compare the vibration signals generated after the four charge structures are detonated. It is found that the vibration signals can reflect the time difference between the different stages of explosive detonation, and the peak vibration velocity is concentrated between 0 and 0.15 s. At this time, the corresponding blast-hole is the cutting hole. is is because there is only one free surface in the cut hole blasting, which is most affected by the surrounding rock, and the energy transferred into the surrounding rock after the explosion is greater. In addition, to provide a good free surface for auxiliary hole blasting, the amount of charge in the cut hole is larger than that in other blast-holes, so the peak vibration velocity in the vibration signal often occurs in the cut hole blasting stage.

Comparative Analysis of Peak Vibration Velocity.
We compare the peak vibration speeds after denoising the blasting vibration signals of the four charge structures ( Table 3).
As can be seen from Table 3, the peak vibration velocities of four charge structures after explosion are different for a given amount of charge. e peak vibration velocities (in descending order) are as follows: the peak vibration velocity of a water interval charge at both ends is 4.14 cm s − 1 , the peak vibration velocity of a water interval charge at the orifice is 3.71 cm s − 1 , the peak velocity of conventional charge structure is 3.52 cm s − 1 cm/s, and the peak vibration velocity Advances in Civil Engineering 5 of water interval charge at the bottom of hole is 3.19 cm s − 1 cm/s. e blasting vibration signal is an intuitive expression of the propagation of the blasting stress wave in the surrounding rock. Its peak value reflects the kinetic energy of the surrounding rock vibration, that is, the energy acting on the surrounding rock after detonation of the explosive. From the above, the energy utilisation of the explosives of the four charge structures can be preliminarily determined. When the water interval charge at both ends is used for construction, the water medium at the bottom of the hole acts like an explosive because of the incompressibility of water [25]. It can apply the most explosive energy to the surrounding rock, so its energy efficiency is highest among the four charge structures tested here. When the water interval charge at the hole bottom is used, as all of the water is at the    bottom of the hole, the stemming length of the hole is slightly shorter than that of the other three kinds of charges and explosive near the orifice, which easily leads to the loss of energy at the orifice after explosion and the minimum utilisation rate of explosive energy. e blasting energy utilisation rate of the water interval charge at the orifice is higher than that of the conventional charge structure. Its energy utilisation rate is between the water interval charge at both ends and the water interval charge at the hole bottom. e energy utilisation ratio of four charge structures is (in descending order) as follows: water interval charge at both ends, water interval charge at the orifice, conventional charge structure, and then the water interval charge at the hole bottom.

Signal Analysis Based on HHT
To analyse the blasting vibration signals more accurately, the HHT method is used to process the signals. HHT includes two parts: empirical mode decomposition (EMD) and the Hilbert transform. First, the signals are decomposed into multiple intrinsic mode functions (IMFs) by EMD, and then a Hilbert transform is applied to each IMF component thus generated. Finally, the energy distributions of IMF components, Hilbert energy spectrum, marginal spectrum, and instantaneous energy spectrum of the signals are obtained.

IMF Layer Energy Distribution Map.
e energy distribution characteristics of blasting vibration signals reflect the relationship between the energy and frequency of blasting signals. Figure 8 shows the energy distribution of each IMF layer after HHT transformation of blasting vibration signals of four charge structures before and after noise reduction. e magnitude of the amplitude corresponds to the energy level of the IMF component. Each IMF component has a different amplitude and frequency: IMF � 1 has the highest component frequency and relatively little energy, which represents the high-frequency noise in the blasting vibration signals and is the main object of noise reduction; 2 ≤ IMF ≤ 5 contains most of the energy of the signals and belongs to the dominant frequency band of the original signals; 6 ≤ IMF ≤ 9 contains less energy after signal decomposition, which is inherent to the signal itself. Comparing the IMF stratification energy before and after noise reduction, it can be seen that the wavelet threshold can effectively remove the noise at IMF � 1. e amplitude of IMF stratified energy before and after structural noise reduction is listed in Table 4.
After noise reduction, the maximum amplitude of the conventional charge structure is 1244.93 cm 2 /s 2 , the maximum amplitude of the water interval charge at both ends is 2733.52 cm 2 /s 2 , the maximum amplitude of the water interval charge at the orifice is 1560.04 cm 2 /s 2 , and the maximum amplitude of the water interval charge at the hole bottom is 1081.49 cm 2 /s 2 . e energy of the conventional charge structure, water interval charge at the orifice, and water interval charge at the hole bottom are mainly concentrated within IMF2 to IMF4, while that of the water interval charge at both ends is mainly concentrated within IMF2 to IMF5. Compared with the other three charge structures, the blasting vibration energy is also distributed in the low-frequency band when the water interval charge at both ends is adopted. It is necessary to pay attention to the influence on the structure to avoid resonance and the damage associated therewith. Figure 9 shows Hilbert energy spectra of blasting vibration signals of the four charge structures: the more dense the points in the figure, the greater the energy. It can be seen that the fluctuations in energy of the blasting vibration signals generated by the four   8 Advances in Civil Engineering the worst in the above range, and a sparse area appears therein. By analysing the Hilbert energy spectrum, it can be explained that the explosive energy can be fully utilised when a water interval charge at both ends is used. e energy utilisation rate of a water interval charge at the orifice is slightly lower than that of a water interval charge at both ends, and use of a conventional charge structure or water interval charge at the hole bottom leads to a loss of explosive energy.

Hilbert Marginal Spectrum.
e marginal spectrum can be obtained by integrating the corresponding amplitude of the frequency over the whole time period. It is the sum of the amplitudes of the frequency in time, reflecting the total energy and energy concentration of frequency in the whole process of time. It can be seen from Figure 10 that the energy of the four charge structures is mainly concentrated in the range of 0 Hz to 100 Hz, and the energy gradually decreases as the frequency increases, which is consistent with the conclusion reached when analysing the energy spectra; however, compared with the other three charging structures, the vibration signal energy of the water interval charged at both ends is not only concentrated in the range of 20 Hz to 70 Hz but is also concentrated in the low-frequency range (6 Hz to 11 Hz), corresponding to IMF5 in the energy distribution of the IMF. Because the natural frequency of the general buildings is from 3 Hz to 9 Hz, when using the water interval charge at both ends for blasting, it is necessary to increase the main vibration frequency of the original vibration signals by reducing the amount of charge in the single section charge and adopt frequency shift technology [26,27] to avoid the natural frequency of the structure and reduce resonance-induced damage thereof.

Instantaneous Energy Distribution Map.
Instantaneous energy is the accumulation of blasting vibration energy over the whole time period and reflects the changes therein. Figure 11 shows a comparison between the vibration signals and the instantaneous energy spectrum after denoising the signals from the four charge structures. It can be seen that the instantaneous energy is consistent with the blasting vibration signals. When the vibration signals change abruptly, the instantaneous energy also changes abruptly and reaches its peak value. We compare the instantaneous peak energies of each explosive charge configuration and the total energy in each section of the four charge structures after detonation (Table 5). Because there is almost no instantaneous energy after the explosion of the 15 segment explosive, it is not counted. It can be seen from Table 5 that, after the explosion of each segment of the blast-hole, the instantaneous total energy is (in descending order) as follows: the water interval charge at both ends, the water interval charge at the orifice, the conventional charging structure, and the water interval charge at the hole bottom. e corresponding maximum instantaneous total energy is 50.94 cm 2 Hz/s 2 , 22.8 cm 2 Hz/s 2 , 18.46 cm 2 Hz/s 2 , and 11.61 cm 2 Hz/s 2 , respectively. is shows that, for a given amount of charge, the water interval charge at both ends can better transfer the explosive energy to the surrounding rock for breaking, which is consistent with the conclusion arising from our analysis of the energy spectra of the blasting vibration signals of the four charge structures.

Blasting Effect.
e blasting effects of the four charging structures are, respectively, calculated. e statistical contents include unit explosive consumption, blast-holes utilisation rate, and rock fragmentation. e results are shown in Table 6.
It can be seen that, with the same blasting volume, the water interval charge at both ends is superior to the other three charge structures in terms of unit explosive consumption, blast-holes utilisation rate, and rock fragmentation.

Dust Suppression Effect.
In the construction process, the dust concentration of the four charging structures after blasting is monitored by the method of filter membrane weighing. In order to ensure the accuracy of monitoring, the ventilation shall be suspended after blasting, and the dust collector shall be turned on. e acquisition parameter is 120 L/min, and the acquisition time is 10 min. See Table 7 for the results.     12 Advances in Civil Engineering Compared with the above four charge structures, the dust concentration of conventional charging structure is 158.75 mg/m 3 , which is the highest among the four charge structures, indicating that water seal blasting can reduce the blasting dust. e dust concentration of the water interval charge at the hole bottom is 140.83 mg/m 3 , which is the highest in the water-sealed charge structure. e dust concentration of the water interval charge at both ends after blasting is slightly higher than that of the water interval charge at the orifice, but the difference is only 1.67 mg/m 3 . Combined with the results obtained in the first half of the paper, the water interval charge at both ends has the advantages of reducing unit explosive consumption, improving blast-holes utilisation rate and rock fragmentation, so the advantages of using the water interval charge at both ends are significant.

Conclusion
(1) e wavelet threshold method can effectively eliminate high-frequency noise in the blasting vibration signals and retain information carried by the vibration signal itself. e frequency and energy of the blasting vibration signals of the water interval charge at both ends are densely distributed in the range of 0 s to 0.9 s and below 100 Hz. e frequency and energy of the blasting signals of the other three charging structures are reduced within the same range, sparse areas appear, and the instantaneous total energy is smaller than that with a water interval charge at both ends, which shows that the water interval charge at both ends can effectively apply explosive energy to the surrounding rock and reduce energy loss in the explosive.
(2) e blasting vibration signal energy of the water interval charge at both ends is mainly concentrated in components IMF2 to IMF5, and the corresponding frequencies are concentrated at 6 Hz to 11 Hz and 20 Hz to 70 Hz, while the blasting vibration signal energy of three other charge structures is mainly distributed in components IMF2 to IMF4, corresponding frequencies are concentrated within 20 Hz to 70 Hz, and the distribution at low frequencies is not obvious. (3) It is the best way to use the water interval charge at both ends for tunnel excavation blasting, but it is necessary to increase the main vibration frequency of the original vibration signals by reducing the single section charge and using frequency shift technology to avoid the natural frequency of the structure and reduce resonanceinduced damage.

Data Availability
e data used to support the findings of this study are included within the article.

Conflicts of Interest
e authors declare that there are no conflicts of interest regarding the publication of this paper.