Discriminant Model of Coal Mining Microseismic and Blasting Signals Based on Waveform Characteristics

Generally, there are two important types of microseismic (MS) signals caused by mining and blasting activities at coal mines. The waveform characteristics of MS signals using FFT, STA/LTA method, and envelope analysis were studied to distinguish these two types of MS signals. The main results are as follows: the dominant frequency and duration of two types of signals are significantly different. The following peak envelope curves of two types of MS signals fit a power function.The power exponent was obtained to describe the attenuated speed of theMS signals.The attenuation of the coalminingMS signals is slower andmore fluctuant than that of the blasting signal. Waveform characteristics consisting of the dominant frequency, duration, and attenuation coefficient were extracted as the discriminating parameters. The discriminating performance of these parameters was compared and discussed. Based on the waveform characteristics, a discriminant model for coal mining MS and blasting signals was established by using Fisher linear discriminant method and its performance was checked.The results show that the accuracy of the discriminant model is more than 85%, which can meet the requirements of MS monitoring at coal mines.


Introduction
The rockburst is a sudden, violent fracture of rock mass in tunnels and mines, generally caused by failure of highly stressed rock and the rapid instantaneous release of accumulated stain energy [1][2][3].The rockburst often causes coal mass failure, roadways collapse, casualties, mining equipment, ground building damage, and so on.The occurrence of the rockburst is accompanied by many small energy microseismic activities which are of great importance for rockburst research [4].A MS monitoring system installed in a mine can monitor the state of the rock mass by analyzing the MS events caused by the mining activities.The core parts of the MS monitoring system include the layout of the monitoring network, data processing, MS source location, and MS early warning [5].One of the challenges of the data processing is to identify the MS signals caused by coal mining activities among the many interference signals [6][7][8].In these complex waveforms, the blasting signals are dominant and can be easily confused with the coal mining MS signals.Thus, distinguishing between the coal mining MS signals and blasting signals becomes the main task of data processing and worth further study.
Recently, MS monitoring systems were widely installed in coal mines in China, and considerable MS data has been recorded.However, little attention has been paid to the discrimination of coal mining MS signals and blasting signals and there is also lack of related literature based on waveform analysis.By contrast, the discrimination of earthquake signals against those of nuclear explosions has been studied extensively.The widely used methods of discrimination can be divided into two categories.One is based on source or waveform characteristics which includes the source location, source depth, origin time, direction of first motion of P wave, ratio of surface wave magnitude to body wave magnitude, ratio of P wave first arrival amplitude to maximum amplitude, ratio of Love wave amplitude to Rayleigh wave amplitude, and ratio of P wave spectral amplitude to S wave spectral amplitude [9][10][11].The other is based on statistical identification which includes Fisher discriminant classifier, logistic regression, unascertained measurement, and neural networks [12][13][14].
Studies on MS signals were also conducted in non-coal mines.Dong analyzed the source parameters of blasting and MS events of three mines in Australia and Canada and established three statistical discriminant models by the Fisher classifier, the naive Bayesian classifier, and the logistic Regression.The classification performances and discriminant precision of three statistical techniques were discussed and compared [8,15].Ma et al. analyzed seven parameters to distinguish between mining MS signals and blasting signals in a phosphorite mine, based on which, two statistical discriminant models were proposed by the Bayesian classifier and the Fisher classifier.Their models were more accurate than traditional methods [16,17].Based on the works of Booker and Taylor, Malovichko applied a multivariate maximum-likelihood Gaussian classifier technique to discriminate between mining MS and blasting events.Analysis of a large number of events showed that approximately 20% (1431 out of 7053) of MS mining events were reclassified as blasting events [7].Vallejos and Mckinnon applied logical regression and neural network modeling to classify mine earthquake and blasting events.The location error and 13 seismic parameters provided by full-waveform system of ESG were used to calibrate the models.The logical regression and neural network approaches can effectively classify different event categories under nonlinear conditions [18].Zhao et al. decompose the signals in time and frequency domains, and time-frequency characteristics of two kinds of waveforms are analyzed.The results showed that the energy of MS signals was mainly concentrated in the lower frequency band than blast vibration signals [19].Pan et al. proposed a comprehensive analysis method to identify and calibrate microseismic events by using STA/LTA method and waveform information [20].Zhu et al. analyzed the fractal dimension box of specified frequency bands through the wavelet analysis and the fractal theory and obtained 23dimensional values of pattern recognition feature vector.Finally, a SVM network model was established to recognize MS events, but the efficiency of identification still needs further improvement [21].Jiang et al. presented a strategy for classifying local multichannels MS waveforms, triggered by a single event.However, the method required a lot of computation and the efficiency of identification needs further improvement [22].Li et al. analyzed the characteristics of coal mining MS and blasting signals using the Hilbert-Huang Transform (HHT) method, but, no clear line of demarcation between these characteristics of the two types of signals was obtained [23].
As the geological conditions in coal mine areas are complex and the transmission medium is not homogenous, the MS signals in coal mines are more complicated than those in non-coal mines and natural earthquakes.Thus, it is difficult to differentiate the two kinds of signals with high precision.Moreover, because most MS monitoring sensors used in coal mines are single component and the monitoring region is small, the difference between the propagation time of the P waves and S waves produced by the rupture of coal and rock is small.The S waves are likely to overlap the P coda waves, making it difficult to detect them.Therefore, the recognition method, utilizing the P/S-wave energy ratio which is widely used to identify nuclear explosion signals and natural earthquake signals, is not suitable for coal mines.In this paper we analyzed three significant characteristics of these two types of signals using the Fourier transform (FFT), short-term averaging/long-term averaging (STA/LTA), and envelope analysis.Based on the characteristics, a discriminant model for coal mining MS signals and blasting signals, with the help of Fisher linear discriminant method, is established, which are useful to improve the efficiency of automatic discrimination and reduce the workload of artificial discrimination.

Outline of the Mine and MS
Monitoring System  the Qianqiu coal mine is equipped with 16 MS sensors, S1, S2, . .., S16.The layout of the MS monitoring network is shown in Figure 1.

Sample Database.
The MS events recorded by the ARAMIS MS monitoring system can be divided into two main categories: (1) coal mining MS events, which refer to the vibration caused by mining; they contain abundant information on the fracturing of the coal and rock; (2) blasting events, which refer to the vibration caused by blasting.The type of recorded event was manually confirmed and marked with different symbols in the database of ARAMIS system.To establish discriminant model of coal mining MS and blasting signals, fifty coal mining MS events and fifty blasting events were randomly selected as the training sample; the other fifty coal mining MS events and fifty blasting events were selected as the testing sample.As the distance between sensors and hypocenter increases, the attenuation and changes of signals become stronger.Therefore, the earliest available signal, received by the nearest sensor from the source, should be selected as the representative signal of the event.Finally, a total of one hundred training sample signals and one hundred testing sample signals are obtained.

Characteristics of the Coal Mining MS and
Blasting Signals   3).We define  start as the arrival time and  end as the terminate time of a wave.For the sample frequency   , the duration is given by The arrival time and termination time can be determined by the short-term average/long-term average (STA/LTA) method as [27][28][29] where  is the sampling moment, which represents the common position of short time window and long time window. ranges in each of short time and long time window, and it is used to calculate average values of signal amplitude in the windows.  is the length of the long time window, and   is the length of the short time window.CF() represents the characterization function at point .In practice, CF() is often replaced with the signal energy or the absolute value of amplitude.When the ratio   >  ( is the trigger threshold), an event is detected.Thus, automatic recognition and arrival picking can be implemented.In this study, the values for  start are   = 50 (sampling points),   = 20 (sampling points), and  = 2 and the values for  end are   = 50 (sampling points),   = 20 (sampling points), and  = 1.5.

Duration Statistics of the Coal Mining MS and Blasting
Signals.Taking one of the coal mining MS signals as an example, the STA/LTA time series of the arrival time and terminate time are shown in Figure 3.The windows move from the first point to the last point for picking  start , and the windows move from the last point to the first point for picking  end .The coal mining MS signal starts at the 920th data point and ends at the 1398th data point in Figure 3.So the duration is  = (1398 − 920)/500 = 0.956 s according to (2).

Attenuation Characteristics of the Coal Mining MS
and Blasting Signals amplitude) to the corresponding point of terminate time as the attenuation of the MS signal.The following peak attenuation curve is proposed for quantitative analysis of the MS signal attenuation.In Figure 5,  high is the time of peak point,  end is the termination time, the envelope fitting curve from  high to  end is defined as the following peak attenuation curve.
To investigate the following attenuation curve, the envelope curve should be obtained first.The steps to obtain the envelope curve are as follows: (1) Picking all the peaks of the waveform (2) Interpolating the discrete peaks using cubic spline interpolation to obtain a complete envelope curve.
The cubic spline interpolation [30,31] is based on the following method: the interval [ ] is divided as  ≤  0 <  1 < ⋅ ⋅ ⋅ <   ≤ .For the point   ,   = (  ), if the function () is cubic polynomial in every small interval [  ,  +1 ] and also has second-order continuous derivatives in [, ], () is regarded as a cubic spline function in after derivation, the expression of () is To study the function form corresponding to the following peak envelope curve of the coal mining MS and blasting signals recorded in the mine and then to describe the signal attenuation process quantitatively, it is necessary to fit the following peak envelope curve of the signals.Gao et al. found that the shock vibration energy in rock-soil medium diminished in a power relation with respect to the travel distance [32].This power function is adopted here to investigate the attenuation of the recorded signals: where  is the amplitude of the signal,  is the sampled point, and  and  are the fitting parameters corresponding to the peak value and attenuation rate, respectively.Generally speaking, the energy of the signal decreases more rapidly when the absolute of  (||) increases; thus, || is called the attenuation coefficient.The fitting precision is represented by the adjusted -squared value, which is a calibration correlation coefficient and is widely used in Origin and MATLAB.A coefficient value close to 1 indicates a highfitting precision and stable attenuation of the signal's energy.

Statistics of Attenuation Coefficient and Fitting Precision.
Taking the coal mining MS signal in Figure 3 as an example, the envelope curve of the waveform is given by the cubic spline interpolation and the termination time point is determined by the STA/LTA method.The fitted envelope curve stretching from the maximum amplitude point to the termination time point is shown in Figure 6.The attenuation coefficient of the signal is 27.73 with a fitting precision (Adj.-square) of 0.96, indicating that the following peak envelope curve satisfies the power function.
The fitting precision and attenuation coefficient of 100 training sample signals are obtained as shown in Figure 7. Figure 7(a) shows the distribution of the fitting precision of the following peak attenuation curve.Of the 100 sampled signals, 93% show a fitting precision above 0.8, indicating that the attenuation curve of both types of signals satisfy the power function.The fitting precision of the following peak attenuation curve of the coal mining MS signals is mainly within 0.8∼1, with 42% between 0.8 and 0.9 and 48% between 0.9 and 1.The fitting precision of the following peak attenuation curve of the blasting signal is higher, with 84% of the signals falling in the range of 0.9∼1.Therefore, in most   The basic principle of Fisher linear discriminant analysis is to project multidimensional and multiclass data to the same direction and separate two or more classes [15,17].In fact, it is a dimension reduction process.Suppose the number of collectivity is , that is,  1 ,  2 , . . .,   .The corresponding mean vectors and covariance matrices are  (1) ,  (2) , . . .,  () and  (1) ,  (2) , . . .,  () , respectively.The number of   samples is selected from the total samples, as follows [33,34]:

Discriminant Model of the Coal Mining MS and Blasting Signals
()  is the projection of  ()  on the axis, which is defined as where  () and  are the intergroup average and the total average, respectively.The intergroup difference is as follows: where  is the maximum eigenvalue 's corresponding eigenvector and   is the deviation matrix of   samples  ()  ( = 1, 2, . . .,   ) in   samples.The intergroup difference is as follows: To let the discrimination function better distinguish samples from different population samples, the difference between two groups from two population samples is expected to be bigger while the intergroup difference is expected to be smaller.Above all, For this, the Lagrange's multiplier is applied and the equation of  =    − (   − 1) can be obtained.The partial differential result for  is 0: By calculating the above formula, the following result could be obtained.
Then the coefficients of the discrimination function could be calculated. is the ratio of sum of squares of deviation within the group to the sum of squares of deviation between groups.By solving the equations, -2 discrimination functions could be obtained ( = min{ − 1, }).Generally, most of the information of the sample could be explained by the first equation.If the first discrimination function is difficult to distinguish, the following discrimination functions could be applied for comprehensive consideration.
Writing the function  =    and plugging  () into  =   ,  () could be obtained.Then values of  () are ranged from small to large; for example, if  (1) <  (2) < ⋅ ⋅ ⋅  () , the threshold value between   and  +1 is as follows: Take the training sample data of the known type into the above formula and predict its type.For the given sample , if  =    is between   ( − 1, ) and   (,  + 1),  belongs to   .Importing the characteristic matrix and categorical data of samples into SPSS to establish the model, partial results are as follows.

Discriminant Model. To use Fisher linear discriminant analysis to get the discriminant function of coal mining MS and blasting signals and establish discriminant model, it is
(1) Tests of Equality of Group Means.Table 1 shows the tests of equality of group means.It reflects these two types of signals' mean differences in different indexes.The probability  (Sig.'s corresponding data in Table 1) is all less than 0.001 in the detection of the coal mining MS signals and blasting signals.
It is indicated that there are significant differences between these two types of signals in each index, so the discriminant analysis could be carried out.
(2) Structure Matrix.The structure matrix in Table 2 shows signals' dominant frequency, duration, and attenuation coefficient and the final discriminant functions' correlations between groups.It could be seen that the bigger the correlation coefficient's absolute value, the greater the effects of the variable on these two types of signals' discriminant functions.As shown in the table, the dominant frequency has the greatest effects on the discriminant functions followed by duration.The attenuation coefficient's impact is minimal, but the correlation coefficient still reaches 0.448.
(3) Central Values of the Final Discriminant Function in the Coal Mining MS and Blasting Signals.  the coal mining MS signal.When the discriminant score is close to 1.442, the signal will be considered as the blasting signal.Due to the only discriminant function, the obtained scores change in the 1D line, and the values of the coal mining MS signal and blasting signal are symmetric about the origin.
Then the demarcation point of the scores could be considered as the origin and the discriminant process could be simplified.
When the signal's discriminant score is less than zero, the signal will be determined as the coal mining MS signal, whereas the signal is considered as the blasting signal.The discriminant result is more reliable with the score larger or smaller than zero.
(4) Discriminant Function Coefficients of Coal Mining MS and Blasting Signals.Table 4 shows the discriminant function coefficients of coal mining MS and blasting signals through modeling.With   representing the signal's dominant frequency,   representing the duration, and   referring to the attenuation coefficient, the discriminant function formula of the coal mining MS and blasting discriminant model could be expressed as follows: where (  ,   ,   ) refers to the discriminant score of the signal.If (  ,   ,   ) < 0, the signal is determined as the coal mining MS signal; if (  ,   ,   ) > 0, the signal will belong to the blasting signal.The discriminant result is more reliable with (  ,   ,   ) larger or smaller than zero.

Performance of the Discriminant Model.
Table 5 shows the eigenvalues of the discriminant analysis.The eigenvalues represent the amount of information, so the larger the eigenvalue, the higher the discriminating performance of the function.To further check the performance of the discriminant model, 100 testing sample signals were identified and classified by using this model and then compared with the real category.Besides, two other classifiers (Bayes and Random forest (RF)) were applied to establish discriminant models using 100 training sample signals and check the accuracy of the discriminant models using 100 testing sample signals.The discriminant results of three classifiers for 100 testing sample signals were compared with the discriminant results using single characteristic.As shown in Figure 8, the discriminant results using multiple characteristic indexes are much better than the single indexes.The performance of the Fisher, Bayes, and RF will have little difference as long as the proper characteristic indexes are used, which indicates that selecting proper characteristic indexes is one of the key steps to distinguish different types of MS signals before establishing discriminant model.Take the discriminant results of Fisher as an example, among 50 groups of coal mining MS signals, 7 groups of signals are wrongly recognized as the blasting signals.Among 50 groups of blasting signals, 5 groups are wrongly recognized as the coal mining MS signals.It could be obtained that the correct recognition rates of the model for the coal mining MS signals and blasting signals are 86% and 90%, respectively.The average accuracy of Fisher for two types of signals is 88%, which is higher than the other two classifiers (the average accuracy are both 87%).Thus, this Fisher discriminant model has relatively high accuracy and stability and its application could be promoted.

Application of the Discriminant Model.
According to the analysis results of field MS monitoring at Qianqiu coal mine, a typical destructive mining MS event with energy more than 10 7 J and a roof-break blasting test occurred at number 21141 longwall panel, respectively.Figures 9 and 10 are the waveforms of these MS events recorded by the ARAMIS MS monitoring system.As shown in Figure 9, there are 15 waveforms, among which, the waveform of S16 is not complete and its zero line drifts.Therefore, there are 14 waveforms involved in the source location.As shown in Figure 10, there are 16 waveforms, among which, S13 and S15 have low signal-noise ratio (SNR) and the arrival time and the terminal time could not be picked up.Therefore, the other 13 waveforms are involved in the source location.
The sensor with the earliest arrival time usually has the smallest travel distance from the source, and the recorded waveform has high SNR and is suitable for distinguishing the type of MS signals.For the destructive mining MS and blasting events, sensor S7 and sensor S3 had the earliest arrival time, and the signals of S7 and S3 were selected as the analysis object, respectively.
The analysis results of the destructive mining MS event are shown in Figure 11.The arrival time and the terminal time of sensor S7 are obtained by the STA/LTA method (the 1152th sampling point and the 2906th sampling point, resp.).Therefore the signal duration is (2906 − 1152)/500 = 3.51 s. Figure 11 The analysis results of the blasting event are shown in Figure 12.The signal duration is (1686 − 990)/500 = 1.39 s (Figure 12

Conclusions
Three significant characteristics of the coal mining MS and blasting signals were studied by using FFT, STA/LTA, and envelope analysis.Based on these characteristics and Fisher method, a discriminant model for these two types of signals was established.The conclusions are as follows.(2) Waveform characteristics consisting of the dominant frequency, the duration, and the attenuation coefficient were extracted as the discriminating parameters.And the boundaries between the two types of signals on the three characteristics were obtained.At the same time, the discriminating performance of the three characteristics was compared and discussed.The dominant frequency is most significant in distinguishing between two types of signals followed by duration.The attenuation coefficient's impact is relatively minimal.
(3) Based on the waveform characteristics, a discriminant model of coal mining MS and blasting signals was established using Fisher method.Its performance was checked and compared with other methods by using the real data from Qianqiu coal mine.Results show that the correct discriminant rate of the coal mining MS and blasting signals in the testing samples is 86 percent and 90 percent, which can basically meet the requirements of the MS monitoring data processing at the coal mine site.Besides, results also demonstrated that the identification results using multiple characteristic indexes are much better than the single indexes.The performance of different recognition methods such as the Fisher method, Bayes method, and Random Forest will have little difference as long as the proper characteristic indexes are used, which indicates that selecting proper characteristic indexes is one of the key steps to distinguish different types of MS signals before establishing discriminant model.

Figure 1 :
Figure 1: Number 21141 longwall panel and layout of the ARAMIS M/E system sensors network.

Figure 2 :Figure 3 :
Figure 2: Distribution of the dominant frequency of coal mining MS and blasting signals.

Figure 4 :
Figure 4: Duration distribution of coal mining MS and blasting signals recorded in the coal mine.

Figure 4
shows the durations of 100 training sample signals.The duration distribution of all the signals is in the range of 0.5∼6 s; most of the coal mining MS signals last longer than the blasting signals.With 2 s as the dividing line, 44 signals lasted 2∼6 s, making up 88% of the coal mining signals and 35 signals lasted 0.5∼2 s, contributing 70% of the blasting signals.Therefore, the duration of the signal can be another important index to differentiate between the coal mining MS and blasting signals in coal mines.

Figure 6 :Figure 7 :
Figure 6: Envelope curve and following peak attenuation curve of a coal mining MS signal.

Figure 7 (
b) shows the distribution of the attenuation coefficient for the coal mining MS and blasting signals, which is generally within the range of 2∼25.The attenuation coefficients of the coal mining MS signals are smaller than those of the blasting signals, which demonstrates the slower damping of the coal mining MS signals.Among all the coal mining MS signals, the attenuation coefficients of 39 signals (78%) are located in the range 2∼8.Among all the blasting signals, the attenuation coefficients of 35 signals (70%) are located in the 8∼25 range.This demonstrates that the attenuation coefficient can reflect the difference between the two kinds of signals and can be used as an important index to distinguish between the two types of signals.

Figure 8 :
Figure 8: Comparison of discriminant results of three classifiers and single characteristic.

Figure 9 :
Figure 9: Waveforms of the destructive mining MS event.
(b) is the spectrum of the signal, which was obtained by FFT method.It is found that the frequency distribution is dispersed.There are wideband signals with 10-250 Hz frequency.The dominant frequency is 37.76 Hz and the subdominant frequency and the third frequency are 86.7 Hz and 123.1 Hz.As shown in Figure 11(c), the following peak envelope curve (red line) corresponds to the section from the 1344th sampling point to the 2906th sampling point.The declining process of the signal is not steady, which is the typical characteristic of the coal mining MS signals.The attenuation coefficient (||) is 6.14 and the fitting precision (Adj.-square) is 0.85.The discriminant score is −2.26, which means this is a coal mining MS event. (a)).
Figure 12(b) is the spectrum of the signal, and the frequency distribution is relatively dispersed.There are obvious signals in 20-150 Hz frequency range.The dominant frequency is 135.97Hz. Figure 12(c) is the following peak envelope curve and attenuation curve of the signal.The following peak envelope curve (red line) corresponds to the section from the 1110th sampling point to the 1686th sampling point.The declining process of the signal is relatively steady and there are no obvious fluctuations.The attenuation coefficient (||) is 9.04 and the fitting precision (Adj.-square) is 0.96.The discriminant score is 2.19, indicating that the signal belongs to blasting event.

( 1 )
The following peak attenuation curves of the mining MS and blasting signals satisfy the power function, which makes it reasonable to quantify the attenuation and stability of the signals utilizing the corresponding fitting coefficient and fitting precision values.Generally, the fitting precision and the attenuation coefficient of the coal mining MS are smaller than that of the blasting signals, which means that the attenuation of the coal mining MS signals is slower and more fluctuant than that of the blasting signals.

( 4 )
The results are useful to improve the efficiency of automatic recognition for coal mining MS signals and blasting signals and reduce the workload of artificial recognition.Although the established discriminant model in this paper is based on the monitoring data in Qianqiu coal mine, the research results could be
As shown in Figure 2, the dominant frequencies lie in the band of below 150 Hz.The majority of the coal mining MS signals fall within the low-frequency range while most of the blasting signals are in the high-frequency range.Below 80 Hz there are 39 signals, which are 78% of all the coal mining MS signals.Above 80 Hz there are 42 signals, comprising 84% of all the blasting signals.The results demonstrate that the dominant frequency of the coal mining MS signals is generally lower than that of the blasting signals.Thus the dominant frequency can be an important index for distinguishing the two types of signals.
) 3.1.2.Dominant Frequency Statistics for the Coal Mining MS and Blasting Signals.For the signal processing in this work, the MATLAB command FFT was used to perform the frequency spectrum analysis.The dominant frequencies of 100 training sample signals (50 coal mining signals and 50 blasting signals) were obtained and analyzed statistically.3.2.Duration of the Coal Mining MS and Blasting Signals3.2.1.Calculating the Signal Duration.The duration of the signals is the time from the arrival time to the end of undulations (Figure mining MS signals and blasting signals and the data are not subject to the classical statistics, so the pattern recognition has become an effective research method.The pattern recognition methods include Bayes method, minimum distance method, ICHAM, Fisher, SVM, and BEG, in which Fisher linear discriminant analysis has been widely used in seismology.Therefore, this paper introduces the Fisher linear discriminant analysis for the discrimination of coal mining MS signals and blasting signals.
[12][13][14]fully used in various fields.For instance, it has been applied to analyze nuclear explosion and earthquakes and has achieved some results[12][13][14].However, it is seldom applied in the discrimination of coal mining MS signals and blasting signals.Considering that there is a lack of the corresponding geophysical models and equations in the discrimination of coal

Table 1 :
Tests of equality of group means.
frequency, duration, and attenuation coefficient.Meanwhile, the modeling needs one column of data to show the real type of signals in each group.Thus the coal mining MS signals are marked as 1 and the blasting signals are marked as 2.

Table 3
shows the central values of the final discriminant functions in the coal mining MS signals and blasting signals.The result shows only one discriminant function, with a discriminant score of −1.442 and 1.442 in the center of coal mining MS signals and blasting signals, respectively.Therefore, when the signal's discriminant score is near −1.442, the signal will belong to

Table 3 :
Central values of the final discriminant function in the coal mining MS signals and blasting signals.

Table 4 :
Discriminant function coefficients of two types of signals.
The canonical correlation coefficient refers to

Table 5 :
Eigenvalue of discriminant analysis.

Table 6 :
Identification results of training sample.

Table 6 .
Among 50 groups of coal mining MS signals, 45 groups are recognized correctly, and 5 groups are recognized as the blasting signals mistakenly.Among 50 groups of blasting signals, 46 groups are recognized correctly, and 4 groups are wrongly recognized as the coal mining MS signals.Therefore, the correct recognition rates of the discriminant model for the coal mining MS signals and blasting signals are 90% and 92%, respectively, which has relatively higher recognition rate than the single criterion.