# Research on Identification Technology of Explosive Vibration Based on EEMD Energy Entropy and Multiclassification SVM

Authors: Huayuan Ma; Xinghua Li; Qiang Liu; Xie Xingbo; Chong Ji; Changxiao Zhao
Journal: Shock and Vibration
 (2020)

---

## Abstract

In this study, the authors introduced energy entropy as a reference feature into the field of blast vibration recognition classification and achieved good results. On the basis of the previous experimental database, 4 kinds of typical vibration signals were selected to form the sample group (building collapse vibration, surface rock blast vibration, underground tunnel blast vibration, and natural gas pipeline explosion vibration). EEMD (ensemble empirical mode decomposition) algorithm was used to calculate the energy entropy of each signal. Taking eigenvector composed of CEE (components of energy entropy) as input, multiclassification SVM algorithm was used for training and prediction. Prediction accuracy was more than 80%. Compared with BP (backpropagation) neural network algorithm, SVM (support vector machine) algorithm has higher training efficiency. The research results can be used in urban vibration monitoring, identify the nature of vibration source in time, and provide technical support for rapid response of emergency rescue.


---

## Body

## 1. Introduction

The blasting vibration signal is a typical short-time nonstationary random signal. The characteristics of the blasting vibration signal can reflect the blasting type, blasting parameters, and site media information. The feature extraction of blasting vibration signal is the basis of studying blasting vibration prediction, structural response under blasting vibration excitation, and reducing blasting vibration. At present, the characteristics of blasting vibration signals mainly include peak characteristics, spectral features, time-holding characteristics, and energy characteristics. A large number of scholars have carried out in-depth research on the influence of characteristics and blasting parameters on signal characteristics. Min et al. [1] studied the strain law of buried pipeline under collapse vibration load. Through range analysis and variance significance analysis, it was found that the dynamic strain of buried pipeline under impact was more significant than that of collapse height and self-weight. Lu [2] obtained the energy distribution of each frequency band after multiscale lifting wavelet packet decomposition of the measured blasting vibration signal and summarized the energy distribution characteristics of the frequency band of the blasting vibration signal. Based on the dimension analysis theory, Zhong [3] discussed the influence factors of blasting vibration holding time and deduced the prediction formula of blasting vibration holding time.With the deepening of research, a new branch of blasting vibration characteristics has appeared: energy entropy. Shannon referred to the concept of thermodynamics and called the average amount of information excluding redundancy “energy entropy” [4]. Energy entropy is an important concept in information theory and a description of the degree of system uncertainty. It can reflect the degree of random change of signals. In recent years, it has been widely used in mechanical fault diagnosis [5–7], speech signal processing [8], mine micro-earthquake signal recognition [9, 10], and other fields [11–13]. In this paper, EEMD algorithm was used to decompose the energy entropy of vibration signal, and the energy characteristics of each IMF are obtained as the input vector of SVM. Through the analysis and training of many kinds of vibration signals, the accurate prediction of vibration types was realized. The research results can provide technical reserves for urban vibration monitoring system and quick and accurate information of explosion source types for decision-makers of disaster prevention and mitigation command departments.


## 2. EEMD Energy Entropy

The EEMD (ensemble empirical mode decomposition) algorithm is a mature nonstationary random signal decomposition algorithm based on the EMD (empirical mode decomposition) algorithm. The original intention is to decompose any nonstationary random signal into multiple IMFs (intrinsic mode function), so as to provide feasibility for further Hilbert transform. EMD algorithm is different from the traditional wavelet and Fourier transform, its decomposition process has no fixed base, it has adaptive characteristics, and so it is famous for its efficiency and accuracy. But, the frequency aliasing effect of EMD algorithm has been perplexing scholars. The improved idea of EEMD is to add white noise at the beginning of EMD decomposition cycle, so as to solve the problem of frequency aliasing in traditional EMD algorithm.The signal was decomposed by the EEMD algorithm to obtain a plurality of IMFs, which represent the characteristic distribution of the signal energy. Due to the difference of explosion source, initiation mode, and propagation medium of different types of explosion vibration signals, the energy distribution characteristics of IMFs are different.Here we introduce the concept of EEMD energy entropy:(1)HEn=−∑i=1nEiElogEiE,where E is the total energy of the signal and Ei is the energy of each IMF. In this paper, −Ei/ElogEi/E is called CEE (components of energy entropy).Several vibration signals from different sample groups are selected for decomposition, and their respective energy entropies are calculated, as shown in Table1. The signal acquisition parameters are shown in Table 2. A simple statistical analysis of the energy entropy data of four kinds of vibration shows that their statistical characteristics are significantly different. Therefore, the vibration energy entropy can be studied in depth as a characteristic parameter of vibration signal classification.Table 1
Comparison of energy entropy of typical vibration signals.
Signal typePipeline explosion vibrationUnderground tunnel blast vibrationDemolition collapse vibrationSurface rock blast vibrationEnergy entropySignal 12.16561.92901.55781.5580Signal 22.24542.32161.08541.9108Signal 32.06222.27221.68361.7539Signal 42.43682.38131.61462.0222Expectation2.22752.22601.48541.8112Variance0.13720.17580.23520.1745Table 2
Parameters of recorder.
Type of recorderType of sensorNumber of channelsFrequency response range/HzSampling rate/SPSTC-4850Piezoelectric35∼5008000Blast-UMPiezoelectric35∼30010000


## 3. Multiclassification SVM

Support vector machine is a binary linear classifier with supervised learning mode [14]. As shown in the Figure 1, for a given input, the sample set S is linearly separable. Each sample contains multiple features, thus forming a feature space (2-dimensional space shown in the figure). The goal of the algorithm is to find the hyper plane with the largest margin ωX + b = 0 [15]. The decision boundary satisfying the condition actually constructs two parallel hyperplanes as the interval boundary to distinguish the classification of samples [16]:(2)ωX+b≥1,⇒y=1,ωX+b≤1,⟹y=−1.Figure 1
Schematic diagram of SVM.All samples above the upper interval boundary belong to a positive class and those below the lower interval boundary belong to a negative class. The distance2/ω between two interval boundaries is defined as the margin [17]. The samples of positive and negative classes on the boundary of interval are support vectors.Support vector machine uses the hinge loss function to calculate the empirical risk. It adds the regularization module to optimize the structural risk and has good robustness [18]. For the nonlinear classification problem, SVM can learn through kernel function method, and the classification effect is still excellent [19].As mentioned above, SVM is initially a linear classifier. With the in-depth study of scholars, SVM algorithm has been successfully extended to the field of nonlinear separability [20]. The principle is to map the nonlinear separable problem from the original feature space to a higher dimensional Hilbert space by using nonlinear functions [21]. There are hyperplanes [22] in space that can separate positive classes from negative classes as shown in Figure 2.Figure 2
Nonlinear separable case.
(a)(b)The standard SVM computing process was used to construct multiple decision boundaries in order to achieve multiclassification of samples [23]. The commonly used structures are “one-to-all (OAA-SVM)” and “one-to-one (OAO-SVM)”. OAA-SVM establishes M decision boundaries for M classifications, and each decision boundary determines the attribution of one classification to all other classifications. OAO-SVM used in this paper is a voting method [24]. The calculation process is to establish decision boundaries for any two of the M classifications, that is, there are MM−1/2 decision boundaries in total. The final sample category is determined by the score order of decision boundary.


## 4. Application of Experimental Data Analysis

In this paper, vibration data under four typical explosion scenarios are selected.Scene 1: pipeline explosion vibration: at present, the global natural gas pipeline construction mileage has exceeded 370000 km [25]. The pressure of pipeline operation keeps increasing. Due to various reasons, natural gas pipeline explosion accidents occur frequently, which brings hidden dangers to the production safety of human society. Therefore, the author has carried out many full-scale high-pressure natural gas pipeline explosion tests, as shown in Figure 3, and recorded the pipeline explosion vibration data as a kind of sample for this study.Figure 3
Pipeline explosion vibration. (a) OD1422-12Mpa explosion experiment of natural gas pipeline [26]. (b) Vibration signal.
(a)(b)Scene 2: building collapse vibration: with the expansion of the city scale and the exhaustion of the design life of the building itself, the blasting demolition engineering of the high-rise building in the city is increasing day by day. The high-rise building disintegrates in the process of collapse, which causes violent vibration when hitting the ground, a serious threat to surrounding buildings. The research team collected a number of groups of building collapse vibration data as a kind of sample for this study, as shown in Figure4.Figure 4
Demolition collapse vibration. (a) Blasting project of investment Plaza building in Jingzhou, Hubei Province. (b) Vibration signal.
(a)(b)Scene 3: surface rock blast vibration: rock blasting is a key method for the rapid and high-quality construction of large-scale infrastructure projects such as nuclear power plants, airports, tunnels, and ports. The high standard, high difficulty, and high strength of large-scale engineering construction put forward new challenges to the application of blasting technology. The impact of blasting vibration of subsequent projects on the projects that have been put into operation cannot be ignored. The research team obtained a large number of surface rock blast vibration data from the foundation pit excavation project of Fujian Fuqing Nuclear Power Plant for this study, as shown in Figure5. The blast parameters are shown in Table 3.Figure 5
Surface rock blast vibration. (a) “Hualong No.1” foundation pit blasting project. (b) Vibration signal.
(a)(b)Table 3
Parameters of rock blasting.
Hole spacing (m)Row spacing (m)Hole diameter (mm)Hole depth (m)Explosive mass (kg)Resistance line (m)Explosive typeDelay time (ms)5.53.211513852.2Emulsion explosive φ9050Scene 4: blasting vibration in tunneling: because of its strong adaptability to geological conditions and low excavation cost, DBM (drilling and blasting method) is especially suitable for the construction of hard rock tunnel, broken rock tunnel, and a large number of short tunnels [27]. The vibration effect of DBM on the surrounding buildings, especially the subway lines which have been built and in operation, cannot be ignored. The research team had been stationed at the construction site of Nanjing Metro Tunnel for a long time, as shown in Figure 6, and collected the blasting vibration data of metro tunnel as a kind of sample for this study. The blast parameters are shown in Table 4.Figure 6
Underground tunnel blast vibration. (a) Blasting site of Nanjing metro line 3 tunnel. (b) Vibration signal.
(a)(b)Table 4
Parameters of tunnel blasting.
Blasting methodHole spacing (cm)Hole diameter (mm)Hole depth (m)Specific consumption (kg·m−3)Explosive typeDelay time (ms)Cut blasting50401.20.6∼1Emulsion explosive φ3250The above data are collected by TC-4850 and Blast-UM vibration recorders, and the specific sampling parameters are shown in the table.As mentioned above, the data samples in this study were divided into four categories, each of which takes 10 groups of data, 40 samples in total. In each training-prediction cycle, 30 groups of data were randomly selected as training sets, and the remaining 10 groups were detection sets. EEMD decomposition was carried out for each data. According to statistics, the number of IMF of 40 sample data was between 10 and 13. Considering that the main energy components were concentrated in the first several IMF layers, the CEE of the first nine IMF layers was taken and normalized to obtain the input matrix as shown in Table5.Table 5
Input matrix.
SignalSerial numberEigenvectorsCEE1CEE2CEE3CEE4CEE5CEE6CEE7CEE8CEE9Pipeline explosion vibration10.02790.00750.08190.24890.02170.26530.48510.53050.496820.03030.01630.00590.00620.01250.02900.02620.51660.485830.08280.00270.04200.09610.46830.31960.51110.44340.096340.00130.02820.29930.53030.48250.39420.35090.28980.060350.00110.01620.09890.41060.51620.50190.32010.10500.059760.00220.02210.01750.40690.51100.49410.34280.06980.108870.00060.00470.04930.07000.28780.37460.47240.53050.455480.02390.00360.02420.07330.40830.49340.46620.50400.391590.02360.09070.04120.07050.40060.53060.40530.49420.2492100.00490.01720.16970.33500.36980.49490.47250.09110.0287Underground tunnel blast vibration110.53070.53000.30880.39420.08730.04390.01700.00970.0076120.34780.34380.52610.51940.39420.13730.02810.01330.0115130.38440.23990.47190.52650.50420.08530.03370.01180.0143140.39900.42140.53070.48710.23070.11770.10230.06850.0239150.22090.27400.53070.53000.37470.10710.05610.02260.0148160.41890.41150.51880.52220.23830.09620.09810.03320.0099170.53040.50260.29670.33650.36510.23690.11030.03150.0452180.33320.34040.51730.49700.45970.20480.08670.02470.0322190.36100.50300.44220.52870.29820.10740.03770.04630.0178200.42280.53060.41530.37570.39520.18260.08600.04010.0257Demolition collapse vibration210.00310.00400.00160.00040.00160.14780.46800.48960.4417220.00060.00010.00120.00240.10700.40570.28470.21100.0726230.00160.00040.00010.00030.43430.30620.41340.16870.3587240.00100.00030.00020.00010.00140.49670.20370.52470.3866250.00010.00010.00000.00580.37940.49420.49500.33460.0972260.00070.00030.00010.00010.00020.35540.45910.49190.4404270.00060.00020.00010.00010.40030.33240.37920.28810.1832280.00040.00030.00340.02040.05820.07860.47140.53050.2545290.00060.00020.00010.00010.00690.52590.51420.29970.1634300.00100.00020.00010.00030.15720.52830.53070.43870.3204Surface rock blast vibration310.00080.00030.00020.13900.47070.52690.25770.11720.0453320.00100.00030.00020.53070.51080.51450.24420.07310.0359330.00030.00010.08350.43390.51370.52590.13290.04070.0230340.00030.00010.01250.14990.52500.50540.47300.18490.1713350.00090.00040.00020.18700.51880.52980.37800.15590.0485360.00060.00020.00850.17430.29470.38460.22940.05910.0183370.00070.00050.00020.48270.52770.52700.14830.03480.0167380.00260.00180.00080.08780.51840.50110.50900.07740.0280390.00210.00110.00040.16930.46680.53070.15300.04720.0094400.00540.00180.00110.00070.52020.52900.49170.28580.2168The feature matrix of 40 × 9 was input into OAO-SVM to train and test. The results are shown in Figure7.Figure 7
Comparison between the prediction results of SVM and the real situation. (a) First operation. (b) Second operation. (c) Third operation. (d) Fourth operation.
(a)(b)(c)(d)It can be seen from the figure that the prediction accuracy can be maintained above 80% in multiple training and tests. The prediction accuracy statistics of the four types of signals is shown in Table6.Table 6
Prediction result statistics.
Signal typePipeline explosion vibration (%)Underground tunnel blast vibration (%)Demolition collapse vibration (%)Surface rock blast vibration (%)Prediction accuracy5010088.992.3It can be seen from the statistical results that this method has high prediction accuracy for underground tunnel blast vibration and surface rock blast vibration. On the other hand, this algorithm often misjudges the pipeline explosion vibration. The cause of this phenomenon may be that the sample size was too small.The EEMD-SVM algorithm provides us with a new means for vibration classification. Compared with commonly used BP neural network algorithms, the new method has the characteristics of high accuracy and short training time. The comparison between the two methods is shown in Table7.Table 7
Algorithm effect comparison.
SVM algorithmBP neural network algorithmPrediction accuracyTraining time (s)Prediction accuracyTraining time (s)Example 18013170502Example 29012760478Example 39015270615It can be seen from the table that SVM algorithm has certain advantages in prediction accuracy compared with neural network algorithm, while in training efficiency, it was much better than neural network algorithm.Furthermore, the approximate entropy and permutation entropy are calculated on the basis of the inherent modal components as input eigenvectors, which are compared with the results of the algorithm in this paper, as shown in Table8. It can be seen that, in terms of prediction accuracy and learning efficiency, the three algorithms are basically the same. However, due to the increase of approximate entropy and permutation entropy in the process of entropy calculation, the CPU time used in the construction of input eigenvectors is significantly increased, and it is highly related to the embedding dimension (permutation entropy) and reconstruction dimension (approximate entropy) parameters, so the comprehensive efficiency is at a disadvantage.Table 8
Algorithm effect comparison.
EEMD energy entropyApproximate entropyPermutation entropyPrediction accuracyTraining time (s)Prediction accuracyTraining time (s)Prediction accuracy (%)Training time (s)Example 1801358014290129Example 2901248013690131Example 3901567016190167Example 41001198012180118


## 5. Discussion and Conclusion

In this study, through the cross-research of EEMD algorithm and SVM algorithm, combined with the training and prediction of 4 types of typical explosion vibration signals, the relevant conclusions are obtained as follows:(1)
EEMD energy entropy can be used as a significant feature for explosion vibration signal recognition(2)
Compared with the traditional BP neural network algorithm, the method in this paper has greater advantages in training efficiency and prediction accuracy(3)
Compared with permutation entropy and approximate entropy, the method in this paper has no obvious advantages in training efficiency and prediction accuracy, but the comprehensive efficiency of the algorithm has been improved.The new algorithm can provide reference for urban vibration monitoring and help managers to know the type of accidents in a short time to make timely decisions.Considering that the prediction accuracy of pipeline explosion vibration signal is insufficient, the author will further collect and sort out the data and expand the sample capacity, in order to achieve higher prediction accuracy.



---

*Source: 7893925-2020-08-08.xml*