Identification of Microseismic Signals Based on Multiscale Singular Spectrum Entropy

-e accurate identification of effectivemicroseismic events has great significance in themonitoring, early warning, and forecasting of rockburst hazards. However, the conventional identification methods have displayed difficulties in achieving satisfactory results. A microseismic signal identification method which combines variational mode decomposition (VMD) and multiscale singular spectrum entropy was proposed in this paper. -e original signal was firstly broken down into a given number K variational mode components, which are ranked by frequency in descending order. -en, the characteristic pattern matrix was constructed according to the mode component signals, and the identification model of the microseismic signals based on the support vector machine was built by performing a multiscale singular spectrum entropy calculation of the collected vibration signals, constructing eigenvectors of signals. Finally, a comparative analysis of the microseismic events and blasting vibration signals in the experiment proved that the different characteristics of the two kinds of signals can be fully expressed by using multiscale singular spectrum entropy. Experimental results further confirmed the effective identification performance of this proposed method.


Introduction
In recent years, microseismic monitoring has become an advanced and effective means for monitoring rock and coal dynamical disasters. It can be carried out in real time for the continuous and online monitoring of the microseismic activities of rock fractures and thereby has become a useful source of microseismic monitoring data. It is well known that the mining environment is very complex due to the large amounts of noise signals such as the background noise and blasting vibrations. ese noise signals generally cause the inaccurate identification of the microseismic signals in microseismic monitoring systems. In addition, technical personnel are required to manually identify effective microseismic events, which seriously affect the efficiency of the microseismic monitoring systems. Considering the fact that blasting operations often occur in coal mines, and the waveforms of rock fractures and blasting vibrations are very similar, the artificial identification method often makes processing errors, which is difficult to effectively identify.
Presently, the identification methods used for the microseismic signal mainly include the first arrival times of the seismic phases, feature recognitions, and parameter discriminations. e most commonly used seismic phase first arrival methods include a STA/LTA (short-time average/long-time average) method and AIC (Akaike information criterion) criterion. e STA/LTA method [1][2][3] dynamically reflects the amplitude changes of the signals using the energy average of the short-term window of the time domain, along with the energy average of the long-term window. is method is able to automatically detect microseismic events by identifying their first arrival times. e STA/LTA method is known to be simple and fast. However, the error rates have been considered to be slightly high, especially when the signal-to-noise ratio is low. Also, the automatic detection function of the method has been found to be unreliable, and the antinoise performances of the algorithm have been observed to be poor. e AIC method was proposed to distinguish the first arrivals of the seismic phases based on the statistical differences in the wave data before and after the arrivals of seismic waves [4]. It has been observed that the AIC function is sharp in the local minimum peak at the beginning of earthquake phases, and the minimum points of the AIC values are the demarcation points of two stable time series. However, due to the fact that the actual microseismic signals are nonstationary signals, the identification effects of the AIC method tend to be greatly influenced by the window selections and signal-to-noise ratios. e feature recognition method mainly uses Fourier transform, wavelet transform, empirical mode decomposition (EMD), and other time-frequency analysis methods to extract signal features and identify signals. Lu et al. [5] used the Fourier transform to analyze the power spectrum and amplitude frequency characteristics of roof-pressure relief-blasting microseismic signals and coal-seam pressure-blasting microseismic signals and preliminarily identified the different types of microseismic signals in mines. However, the Fourier transform is traditionally used to analyze periodic stationary signals and has not been found to be effective for random nonstationary microseismic signals with spikes and mutations [6,7]. Wavelet analysis can simultaneously perform time and frequency analyses [8]. However, this type of analysis requires the choosing of a suitable wavelet base to achieve better decomposition effects [9,10]. Presently, many researchers have applied wavelet analysis to the waveform characteristics of natural earthquakes [11,12] and energy variation of seismomagnetic signals [13], as well as the analyses of microseismic signals in mines [14], and so on, and have achieved better application results. Empirical mode decomposition, which was proposed by Huang et al. [15], is an algorithmic method used to detect and decompose a signal into principal "modes." is method has been proven to be suitable for handling random nonstationary signals. Scholars have applied it to the noise reduction of mine microseismic signals [16], feature extraction [17], classification identification [18,19], and other fields. However, there have been boundary effects and a mode-aliasing phenomenon observed in the EMD method [20,21], which causes the decomposition results to experience problems, such as instability and nonuniqueness. ese defects in the EMD make it difficult to effectively identify signals. As a novel method for signal identification, the combination of deep learning technology and feature extraction has attracted the interest of many researchers [22,23]. e parameter discrimination used a linear regression fitting method [24,25] and extracted the two slope values of the starting-up trend line and the two coordinate values of the first and maximum peaks as the characteristic parameters. en, it established an identification model by applying a Fisher discriminant analysis [26]. However, this method required a higher signal waveform, which also resulted in difficulties in the identification of signals.
Variational mode decomposition (VMD) is an entirely nonrecursive decomposition model in which the bandlimited intrinsic mode functions (BLIMFs) or modes are extracted concurrently [27]. e model searches for mode ensembles, as well as their respective center frequencies, in which the modes collectively reproduce the input signals, while each becomes smooth after demodulation into a baseband.
e VMD method has a strong theoretical foundation which overcomes the shortcomings of the boundary effects and mode-aliasing phenomenon of the EMD and other recursive decomposition algorithms. Also, the VMD method avoids the difficulties of selecting the wavelet base functions of wavelet or wavelet packet analyses. At the present time, the VMD method has been successfully applied to the functional coupling analyses of electroencephalograms and electromyograms [28], fault diagnoses [29], and feature extractions [30,31], as well as other fields. Singular spectrum entropy has been widely used to analyze a signal's sequence complexity analysis and perform nonlinearity analysis because of its high robustness and antinoise. Multiscale singular spectrum entropy (MSSE) is a combination of multiscale and singular spectrum entropy and is more effective for analyzing the information of a time sequence. It has been applied to the isolated island detection in power grids [32], fault diagnoses [33][34][35], and other fields. e support vector machine is based on a statistical theory that follows the principle of structural risk minimization and considers the empirical risk as well as the confidence risk. In a small sample, nonlinear and high-dimensional pattern recognition problems are used and show many unique advantages [36]. In the case of wired data samples, the data sample features can be well described and have been widely used in fault diagnosis, pattern recognition, and other fields [37,38]. e least-squares support-vector machine (LS-SVM) is based on the regularization theory [39], which solves the quadratic programming problem in the classical SVM algorithm to solve linear equations, improves the convergence speed of the algorithm, and reduces the time needed to perform the algorithm. e above analysis shows that MSSE can reflect the essential characteristics of a onedimensional time series from different aspects, is sensitive to the complexity of signals at different scales, and can fully express the different characteristics of microseismic signals and blasting vibration signals. A LS-SVM classifier can be trained by learning microseismic signal samples and blasting vibration signal samples. A LS-SVM classifier can accurately identify microseismic signals, so this paper introduces multiscale singular spectrum entropy to measure microseismic signal complexity, and the LS-SVM model was used to identify microseismic events.
In view of the shortcomings of the abovementioned multiscale decomposition methods and the limitations of feature extraction methods, in this paper a combination of MSSE and SVM, named MSSE-SVM, is proposed to identify microseismic signals. e signal is firstly broken down into the K modes through VMD, and then the singular spectrum entropy of each mode is estimated using MSSE; finally, the singular spectrum entropy eigenvector is input into the trained classifier to identify the signal type. In order to verify the identifying effectiveness, this proposed MSSE-SVM method is applied to a large number of historical sample data. Compared to EMD-based techniques and energy characteristic methods, e MSSE characteristic can distinguish microseismic signals from blasting vibration signals, and has better accuracy and stability.

2
Shock and Vibration

Variational Mode Decomposition.
e VMD algorithm is a new completely nonrecursive signal decomposition method [27], it can adaptively decompose any signal x into the K discrete number of the subsignals (modes) u k (k � 1, 2, . . . , K), where each mode is a limited bandwidth and compact around their respective center frequencies ω k (k � 1, 2, . . . , K). e constraint condition is that the sum of the bandwidth of all modes is the smallest, and the sum of all the modes is equal to the input signal x. e VMD algorithm suggests to making use of both a quadratic penalty term α and Lagrangian multipliers λ to address the constrained problem, in order to render the problem unconstrained. An alternate direction method of multipliers algorithm was used to obtain the optimal solution of the constrained variational problem, and the input signal is decomposed into each mode and its central frequency adaptively.
e concrete steps of VMD were as follows: Input: the original signal x and the decomposition number K. Output: K modes.
Step 2: n � n + 1, and update u k and ω k according to the following equations: Step 3: update λ according to the following equation: where τ represents the update step parameter of the Lagrangian multiplier.
Step 4: for a given discriminant accuracy ε > 0, if equation (4) is satisfied, then stop the iteration and output the K modes; otherwise, go back to Step 2: where ε is the convergence tolerance level, which can control the relative error effectively. In VMD algorithm, the default value of ε is 1e − 6. In this paper, the value of ε is 1e − 7. Figure 1 depicts as an example of the decomposition of the noisy triharmonic signal with length 1000. It is affected by noise: where η ∼ N(0, 0.1) represents the Gaussian additive noise. Figure 1(b) shows the five IMFs from EMD. Figure 1(c) shows the BLIMFs extracted by VMD, where the number K is set to 5 to coincide with EMD. As shown in Figure 1, the EMD produces 5 estimated modes. e first IMF contains the highest-frequency harmonic and amounts of noise but is severely distorted. e last two IMFs contain the middle and the lowest-frequency harmonics; however, serious mode aliasing has been appeared in these two IMFs. In the 5 BLIMFs extracted by VMD, the first two modes perfectly extracted the lowestfrequency and the middle harmonics. e forth mode contains the highest-frequency harmonic and small amount of noise. Figures 1(d) and 1(e) show that the frequency of BLIMFs is more compact and clear than that of IMFS. ese demonstrate that the VMD decomposition performance is better than EMD.
In general, VMD is a powerful signal decomposition method. It has a solid foundation of mathematical theory and can solve the mode-aliasing problem better. Its operation is carried out in the frequency domain, and the operation efficiency is high. is paper uses these advantages of the VMD method to decompose the two kinds of vibration signals, so as to achieve the purpose of accurate calculation of the signal frequency band singular spectrum entropy and the extraction of multiscale singular spectrum entropy eigenvectors.

Multiscale Singular Spectrum
Entropy. Singular spectrum analysis is an effective time-domain analysis method. It transforms the embedded space into an equivalent orthogonal coordinate system and obtains the signal trajectory in the subspace with the minimum embedding dimension. It eliminates the linear dependence and artificial symmetry between the delay coordinates, enhances the signal-to-noise ratio, and sharpens the singularity of the signal. However, as a time domain analysis method, traditional singular spectrum analysis is not conducive to multiscale monitoring and feature extraction of signal singularity. On this basis, a novel microseismic identification method that combines multiscale analysis and singular spectral entropy is proposed in this paper. By calculating the singular spectrum entropy of the components in different frequency bands, the uncertainty of the distribution of a singular characteristic in the different frequency space of the analyzed signal is accurately reflected.
Given the one-dimensional time series x � x 1 , x 2 , . . . , x N , the number of the total sampling points was represented as N, x was decomposed using the VMD method, and the K modes were denoted as Shock and Vibration u 1 , . . . , u k , . . . , u K . For each mode, a m × n pattern matrix [40] was constructed as follows: e SVD was used to obtain singular values δ k1 ≥ δ k2 ≥ · · · ≥ δ kl of the R k , and the singular value spectrum of the k-th mode u k was formulated. e singular spectrum entropy of mode u k can be defined as follows: where q ki is the proportion of the i-th singular value in the entire singular value spectrum as follows: Applying the above method, a K-dimensional vector can be constructed from a multiscale perspective to describe the singular spectral characteristics of a signal.  Shock and Vibration

Extraction of the Eigenvector.
e VMD method has the ability to adaptively realize the effective separation of the frequency domain and the components of the signal, along with preserving the local characteristics of the signal, which can effectively avoid the mode-aliasing problem. In this study, based on the VMD, the singular value spectrum was obtained by the SVD of each mode, and the corresponding singular spectrum entropy was calculated. en, the extraction of the eigenvector could be accurately realized.
e main steps of this study's signal eigenvector extraction based on the VMD method and multiscale singular spectrum entropy were as follows: Step 1: . , x N was decomposed using the VMD method, and the K variational mode components with different frequency bands were obtained u k � u 1 , u 2 , . . . , u K .
Step 2: e m-dimensional phase space of each mode u k was reconstructed to form the pattern matrix R k , and the singular value spectrum δ k1 , δ k2 , . . . , δ kl of each mode component signal was obtained by SVD. en, the singular spectrum entropy S k was calculated successfully.
Step 3: e singular spectrum entropy eigenvector T � (S 1 , S 2 , . . . , S K ) of the original signal x was obtained by combining the singular spectrum entropy of each mode. T is used as the eigenvector to identify the rock fracture microseismic signal and the blasting vibration signal. e feature extraction method quantificationally described the signal's frequency band information in the form of singular spectral entropy and also presented the determinations of the singular complexity of the original signal in each frequency band. e eigenvectors extracted by this method are repeatable and stable and can distinguish the two kinds of vibration signals.

Support Vector Machine.
e support-vector machine (SVM) is a type of small sample learning algorithm which is based on statistical learning. It is effective in solving the problems related to small numbers of samples and nonlinear and high-dimensional pattern recognition. It has been widely applied in image recognition, fault diagnoses, etc. In this study, a least-squares support-vector machine (LS-SVM) was used to identify microseismic events. e LS-SVM is a model of the SVM which is used to solve KKT (Karush-Kuhn-Tucker) optimization problems. e LS-SVM has the ability to significantly reduce computational costs and improve the convergence speed on the basis of maintaining the standard SVM advantages [41,42]. e basic principle is as follows: For a given training set Q � (x i , y i ) | i � 1, 2, . . . , d , x i is the i-th input data, and y i is the output class of the i-th input data. e optimization problem of the LS-SVM can be expressed as follows: where w is the weight vector, ϑ is the adjustable constant, and b denotes the threshold value. A Lagrange method was used to solve the optimization problem and was obtained as follows: In the formula a � [a 1 , a 2 , . . . , a d ], a i is the Lagrange multiplier, and a i ≥ 0. en, by optimizing equation (10), the partial derivative could be determined as follows: In order to avoid the problem of dimensionality caused by a direct solution, a Gauss radial basis function was introduced as the kernel function as follows: en, after equation (12) was calculated, optimal classification equation (13) could be used to determine whether the waveform was a microseismic signal as follows:

Proposed Identification Model
First, select an equal number of microseismic samples and blasting vibration samples, and then, perform the multiscale singular spectrum entropy calculation on the sample data and extract the eigenvector of each signal to form the training set. e LS-SVM software package is used to train the training set of the sample data to form the LS-SVM classifier. With the help of the LS-SVM classifier, the test samples can be classified. e signal identification flowchart is shown in Figure 2.

Shock and Vibration 5
In this study, the identification algorithm of the microseismic signals based on the VMD, multiscale singular spectrum entropy, and LS-SVM (hereinafter referred to the MSSE-SVM algorithm) was as follows: Input: the sample set Q, regularization parameter ϑ, and kernel function parameter σ 2 . Output: the classification results.
Step 1: Sample selection was done. e appropriate data were selected to form a sample set Q.
Step 2: Eigenvector extraction was done. e singular spectrum entropy eigenvector of each sample signal was extracted from the sample set Q, and the characteristics of the sample signals were represented by the multiscale singular spectrum entropy T � (S 1 , S 2 , . . . , S K ).
Step 3: In order to determine the LS-SVM parameters, a W-fold cross validation method was used to determine the regularization parameter ϑ and kernel function parameter σ 2 . e basic process included the following: dividing the sample set Q into W parts; taking the W − 1 parts as the training set QT; the remaining part served as the verification set QV; the classifier was trained with the QT; and the performance of the classifier was evaluated by the verification set QV. en, the parameters were adjusted and the W times of the above process were repeated, in order to train and test all of the samples in the sample set Q. e parameters with the highest classification accuracy on the W verification set were considered to be the optimal parameters.
Step 4: Output the identification result. e trained LS-SVM classifier is used to classify the input signals and output the types of signals. e W-fold cross validation method is the most common application and has been observed to achieve good results [43]. Finally, by taking into account the computational efficiency and references in other related research results [44], this study selected W � 10, ϑ � 1.8125, and σ 2 � 1.2075. e purpose of developing the MSSE-SVM algorithm is to correctly identify the microseismic signal from the microseismic detection system with high noise and blasting vibration and effectively facilitate the analysis of the microseismic signal such as source location. Compared with the similar algorithm based on EMD method, it has the merits of high noise robustness, simple operation, and suitability for small sample.

Experimental Data Sources and Preprocessing.
e experimental data used in this study were derived from a coal mine microseismic monitoring system located in western China, which had been mainly used for monitoring the rockburst disasters in mines. e seismic detector adopted a single component geophone sensor with a frequency response range between 3 Hz and 2 kHz. e signal sampling frequency was 1 kHz. An artificial pickup is performed to obtain a total number of 100 rock fracture microseismic signals and 100 blasting vibration signals and is recorded as the sample set Q. For the ease of analysis, data form 200 waveforms are interpreted as time sequences of equal length, and each sequence has 3000 sampling points.

Example of the Analysis Process.
In this study, in order to explain the process of the feature extraction based on the VMD, singular spectrum entropy, and SVM classification, typical rock fracture microseismic and blasting vibration signals were selected from groups A and B, respectively. en, the two signals were decomposed using the VMD method, and six variational mode components were obtained. e waveforms of the original signal, along with its variational mode components, are shown in Figure 3. e number of the sampling points N was 3000. e phase space dimension m of the pattern matrix was 300. e singular spectrum entropy of the two signals was calculated, and the results are shown in Table 1. As can be seen from Table 1, among the eigenvectors of the microseismic signal, the singular spectrum entropy of u 3 is the largest and that of u 4 is the smallest. In the eigenvectors of blasting vibration signals, the singular spectrum entropy values of u 3 and u 6 are smaller. In summary, the singular spectrum entropy could be used to characterize the singular state of the signal on the different modes.
In this study, the sample set Q was divided into two groups: A and B. In group A, there were 100 rock fracture microseismic signals observed with clear take-off waveforms, and 100 blasting vibration signals in group B. From groups A and B, 15 microseismic signals and 15 blasting vibration signals were randomly selected, respectively. en, these signals were decomposed by VMD, and the singular spectrum entropy of each component was calculated. e results are shown in Table 2.
e average values of the singular spectrum entropy were calculated, as shown in Figure 4. It was confirmed that after the decomposition of the rock fracture microseismic signals, the singular spectrum entropy values of the two modes u 4 and u 6 were smaller, and also, mode u 3 was the maximum. In regard to the blasting vibration signals, the singular spectral entropy values of the two modes u 3 and u 6 were observed to be smaller, and the two modes u 1 and u 5 were larger. As detailed in Table 2, the singular spectrum entropy of the mode u 6 in the first five sets of blasting vibration signals was approximately 1.5. e results showed that the five sets of blasting vibration signals almost do not include vibration waves of the mode u 6 , so their singular complexity is low and the singularity value is small.
From the above analysis, it was determined that the eigenvectors based on the multiscale singular spectrum entropy were stable in the same types of signals. Moreover, there were obvious differences observed in the singular spectrum entropy eigenvectors of the microseismic and blasting vibration signals, respectively. e singular spectrum entropy was found to be beneficial to the classification and identification of the two kinds of signals. Shock and Vibration 7

Performance of SVM Classifier and Comparative
Examination. In order to test the performance of this study's algorithm, 70 groups of vibration data were selected from the groups A and B, respectively, and the training set QT was composed. e remaining 60 groups of vibration data in the two groups constituted the test set QV. e final classification results for the test set are shown in Table 3.
As can be seen from Table 3, the singular spectrum entropy could effectively characterize the features of the rock fracture microseismic signals. Among the test set QV, 30 rock microseismic signals were 100% correctly classified, and two groups of 30 blasting vibration signals were found to be misclassified, resulting in an overall classification preparation rate of 93.33%. en, using the 60 signals in the test set QV, a comparison between the results of this method and other methods was completed, as detailed in Table 4. It can be seen in the table that this study's method can take on higher robust and better classification comparing with the other algorithms. Although the proposed method and the energy feature method have the same classification accuracy in this test set, when we have more samples, the SVM based identification model has better adaptability and accuracy than the hard threshold method. erefore, the proposed method was confirmed to be suitable for the classification and recognition of rock fracture microseismic signals.

Influence of Singular Spectrum Entropy under Different
Numbers of Sample Points. In this study, the durations of the microseismic events were relatively short, generally lasting for only approximately several seconds. When intercepting the microseismic data, the selection of the number of sampling points was determined to affect the values of the singular spectrum entropy, as well as the accuracy of the final recognition.
In order to illustrate the influences of the number of sampling points on the singular spectral entropy, a typical rock fracture microseismic signal was selected as an example (Figure 3(a)) and the number of the sample points N was 2,000, 3,000, 4,000, and 5,000, respectively. e time series signals of the four groups of different durations of the microseismic signals were referred to as RF1A, RF1B, RF1C, and RF1D, respectively, as shown in Figure 5.
As detailed in Figure 5, when the number of sampling points was 2,000, the tail of the microseismic event was truncated, and the microseismic waveform data were slightly deficient. However, when the number of sampling points was 3,000, 4,000, and 5,000, the microseismic waveform data were observed to be complete. e singular spectrum entropy of each of the modes of the four signals was calculated, and the results are shown in Table 5.
e singular spectrum entropy of each of the modes of the four rock fracture microseismic signals was projected to the histogram for comparison, as shown in Figure 6. In the figure, the minimum amplitudes of the two adjacent signals   8 Shock and Vibration in the same mode on the singular spectrum entropy growth are denoted by the blue digital display, and the maximum amplitude is denoted by the red digital display. As can be seen in Figure 4, the singular spectrum entropy of each mode increased with the increases in the number of sampling points. When the number of signal sampling points changed from 3,000 to 4,000, the maximum increased amplitude difference of the singular spectrum entropy of each mode was 0.1, and the singular spectrum entropy of the modes maintained the same increase difference. When the number of signal sampling points changed from 4,000 to 5,000, the maximum increased amplitude difference of the singular spectrum entropy of each mode was 0.07, and the singular spectrum entropy of the modes also maintained the same increase difference. However, when the number of signal sampling points changed from 2,000 to 3,000, the maximum increase in the amplitude difference of the singular spectrum entropy of each mode was 0.513, and the singular spectrum entropy increases between the modes were observed to be quite different.
In this study, following the training of the LS-SVM classifier, it is usually necessary to normalize each component value of an eigenvector to a [0, 1] interval in order to form a standardized eigenvector. e purpose was to speed up the training speed and improve the classification accuracy. e normalized eigenvectors of the four singular spectral entropy feature vectors detailed in Table 4 are shown in Figure 7. As can be seen in Figure 7, the RF1A was very different from the singular spectrum entropy standard eigenvectors of the other three microseismic signals. However, the singular spectrum entropy standard eigenvectors of RF1B, RF1C, and RF1D were found to be almost the same. It was found that when the number of sampling points was 2,000, the data of RF1A were deficient, and the eigenvector lacked reliability. However, when the number of sampling points was 3,000, 4,000, and 5,000, the data of RF1B, RF1C, and RF1D were observed to be complete, and the eigenvectors were determined to be reliable. erefore, it was necessary to pay attention to the influences of the number of sampling points on the eigenvectors of the singular spectrum entropy. In accordance with the results of the large number of experiments completed in this study, it was suggested that 3,000 sampling points were more suitable for both the rock fracture microseismic signals and the blasting vibration signals.

Investigation of Identification Result under Different K
Values. From the previous discussion, it can be seen that each signal could be used as a six-dimensional singular spectrum entropy vector to express its characteristics. In order to reveal the effect of the K value on identification     results, 100 rock fracture microseismic signals and 100 blasting vibration signals were used in this study, of which K was 2, 4, 6, 8, 10, 12, and 14, respectively. e identification accuracy at the different K values is shown in Figure 8. It is clear in Figure 8 that when the K value was 2, the identification accuracy of the blasting vibration signal was approximately 50%, and the identification accuracy of the rock fracture microseismic signal was approximately 80%. As the K values increased, the identification accuracy was also observed to increase. is was found to be due to the increased K values causing the layer of signal decomposition to increase, allowing for a more accurate expression of the complexity of the signal. erefore, the differences in the rock fracture microseismic and blasting vibration signals were more easily distinguished. However, when the K value was more than 6, the accuracy rate of the identification was observed to be reduced. ese results were due to the fact that when the number of the modes was too high, over decomposition of the signal occurred. When the original signal had only a very small number of signals in some of the modes, the singular spectrum entropy of the mode may have been too high or too low, and the differences in the singular spectral entropy eigenvectors of the two kinds of signals were observed to be affected. In accordance with the results of the experiments conducted in this study, it was suggested that the K value was 6.

Conclusions
In this paper, the identification model of the microseismic signal is established to accurately identify and record the microseismic signal under the influence of blasting vibration signals.
e microseismic signal identification method is proposed with a combination of MSSE and SVM based on the LS-SVM theory, which can effectively distinguish the microseismic signal and blasting vibration signal and has high identification accuracy and stability. We can draw the following conclusions form the results of numerical computation and experiment: (1) Compared with EMD and wavelet analysis, the VMD algorithm captures the relevant center frequencies quite precisely, the resulting signal separation is relatively good, and the frequency of each mode is more compact and clear. e main problem of the application of the VMD algorithm is to determine the value of K, which is related to the signal-to-noise ratio and dominant spectrum distribution of the signal. (2) e number of sampling points of the signals have certain influences on the eigenvector of the singular spectral entropy. It was determined that when the number of sampling points was 3000, the singular spectral entropy eigenvectors of the same types of vibration signals could reach a stable state. (3) e model experiment demonstrated that the multiresolution singular spectrum entropy can be used as a classification feature to classify and identify the rock fracture microseismic and blasting vibration signals. It was verified that the proposed method has high identification accuracy.

Data Availability
All data included in this study are available upon request by contact with the corresponding author.

Conflicts of Interest
e authors declare that they have no conflicts of interest regarding the publication of this paper.  10 Shock and Vibration