Microseismic (MS) frequency response is an important part of high-efficiency data mining to achieve the aim of coal and gas outburst (CGOB) early warning. Based on the variation pattern of acoustic emission (AE) signal in the coal failure process, the experimental characteristics of MS activity and typical signals CGOB were obtained in this study. First, the AE behavior of coal failure experiment was studied, and an explanation of laws was provided as follows: the fracture behavior of coal sample exhibits certain characteristics of AE response in terms of AE event count, signal amplitude, and frequency; each stage has its own physical meaning during the process of loading test. Based on these laws, CGOB experiments were carried out using a large CGOB physical simulation system with a MS monitoring system. Notching filter and wavelet packet transform technique were used in the denoising and feature extraction of six typical MS events (signals). The features of each stage, including the time-frequency domain, were extracted and quantitatively expressed. We finally arrive at the following conclusions: (1) CGOB exhibits significantly periodic characteristics, and each CGOB stage corresponds to the significant response characteristics of MS. CGOB presents varying characteristics, such as “valley-peaks-valley”. (2) From the incubation stage to happen stage of outburst, the spectrum significantly moved from extremely low frequency (100-200 Hz) to high-frequency band (approach to 1600 Hz). During the residual stage, MS frequency manifested the concentration distribution (50 Hz) and offered the advantage of energy concentration. (3) The phenomenon of signal energy also shows the trend of energy transform low to high and to low modes along with the process. Signals total energy distribution (42.81%, 1,437.5-1,812.5 Hz) in the happen stage are markedly larger than those of events in incubation stage (7.01%) and residual stage (1.44%). The methodology presented in this paper for CGOB signal analysis provides a new method to obtain MS response precursor and predict CGOB disaster. This approach can be useful for rockburst anticipation and control during mining in gas and highly stressed coal mines.
National Natural Science Foundation of China5160411551504096Fundamental Research Funds for the Central Universities3142017002Open Foundation of the Hebei State Key Laboratory of Mine Disaster PreventionKJZH2017K14Technology Project and Self-Financing Foundation of Hebei Province16275512Project of Langfang Key Technology Research and Development Program, China20160131161. Introduction
Coal and gas outburst (CGOB) disasters are characterized by a complicated mechanism, sudden occurrence, and massive harm [1–4]. CGOB disasters not only affect normal operations of mining but also seriously restrict mine safety and efficient production. To date, existing parameters and methods are limited in making accurate predictions and fail to satisfy the requirements of mining safety; Figure 1 shows a typical unforecasted CGOB accident that occurred in Guizhou Province. Therefore, prediction and prevention studies should be related to reduce casualties and economic losses caused by CGOB disasters.
Disaster of CGOB in Guizhou Province, China.
Outburst and piling up seam
Dumping and buried single hydraulic props
CGOBs are controlled by multiple factors, such as stress, gas, physical, and mechanical properties of coal [5–12]. Capturing characteristic information before the disaster occurs is crucial to achieve accurate and timely warning. According to the investigation and analysis of accidents, CGOB is commonly accompanied by precursor information [13, 14], such as drilling jamming, gas abnormal effusing, temperature changes of coal body, gas pressure increase, and hissing noise released from coal body. Many researchers suggested using gas pressure, coal bed stresses, and methane concentration parameters based on these characteristic pieces of information to predict CGOB disasters. The main contactless monitoring and early warning methods to predict the CGOB disasters include electromagnetic radiation (EMR) method [15–18], acoustic emission (AE) method [19, 20], and microseismic (MS) method [21]. Compared with AE and EMR, MS inversion technology exhibits advantages in source inversions, monitoring range, and signal connotation aspects over a large size monitoring range [22–24]. MS method has been widely used for mining disaster prediction, such as CGOB, because onsite monitoring and laboratory research indicate that microfracturing is a typical precursor of dynamic disaster [25–27].
MS signal analysis, processing, and feature mining are key to realize accurate and early warning, which is still an urgent problem. In practice, the characteristics of MS event, signal recognition [28], and precise positioning [29, 30] are the important parts of signal process and analysis. The related theory of signal has been advanced for many years. At present, certain methods such as time-domain analysis, frequency-domain analysis, and time-frequency domain analyses have been conducted to capture signal characteristics. Among these methods, time-frequency analysis method is widely used in blasting, machinery, and the other areas of signal analysis, especially seismic signal processing [31–33]. Typical methods, such as Fourier transform, continuous wavelet transform, S transform, Wigner–Ville distribution, and Hilbert transform (HHT), are based on signal frequency analysis with primary short-time Fourier transform [34], wavelet transform [35], discrete wavelet transform (DWT) [36], and wavelet packet transform [37]. Time frequency analysis is a complex signal processing method that can be used for constituent and transient nonstationary MS signals. Kang (2010) analyzed the time frequency of AE signals from rock damage and obtained time-delay estimation in a single band, which resulted in a damage localization with enhanced accuracy [38]. Zhu (2012) also conducted research on wavelet packet energy to classify blasting and MS signals [39]. Guizhou Coal Mine Design and Research Institute has investigated MS signal analysis and prediction of CGOB in recent years and obtained significant results through a series of laboratory experiments [40, 41].
This paper aims to find an effective feature for describing the process changes of CGOB and explore a new approach to disaster prediction. In general, this study can be primarily divided into two parts: uniaxial compression test (UCT) of coal specimens containing gas (CSCG) and large-scale CGOB test in the laboratory. In previous experiments [40, 41], we found that the main frequency of AE signal changes with the coal specimen failure process. On the basis of this result, we present a novel method to analyze the weak and nonstationary transient characteristics of MS signals during CGOB. Through analysis and processing, we observed that AE or MS signal frequency was regularly changed during CGOB. We describe the problem in detail and demonstrate how to find a solution.
2. AE Behavior of Coal Failure Process2.1. AE Dynamic Response under UCT
Previous studies show that the failure process of coal-rock combined body exhibits apparent AE characteristics that are significant when the specimen is initially damaged, yielded, and destroyed [42–44]. The AE characteristics in UCT with different gas pressures should be understood and interpreted to confirm the coal and rock fracture. Therefore, we performed the UCT study on AE characteristics of coal failure with coal samples stored in a gas pressure environment before the CGOB experiment. The experimental study on the mechanical properties of coal containing methane was performed on raw coal samples in a coal mine located in Guizhou Province, China. AE system choice of high precision equipment from PAC (Physical Acoustic Corporation), the parameters of AE device were set as follows: main amplifier, 40 Hz; threshold value, 35 db; sensor resonant frequency, 20–400 kHz; sample frequency, 1 MHz (more details about this system and parameters can be seen in [43–45]). Figure 2 clearly shows that the AE event count changes with loading stress in different stages. During the laboratory experiment, AE signals occurred along with uniaxial loading at every stage, and AE space-time evolution corresponded to the complete stress-strain curves. Clearly, the fracture behavior of coal sample exhibits certain characteristics, and each stage denotes its own physical meaning during loading test. After a gap, the stage of specimen failure arrives. The count of AE event increases sharply along with the loading stress and then reaches a maximum value. Afterward, the process is completed and enters a quiet period.
AE event count and loading stress for different gas pressure levels.
2.2. Peak Frequency Variation of AE Signal
In addition, several new conclusions can be defined as an indications or precursors for the CSCG to sprout and develop during loading. As shown in Figure 3, the frequency of AE signals shows an apparent periodic pattern. Sporadic AE responses with high frequency (> 200 kHz) appeared in the initial period. The process then developed into a brief period of “blank.” When the signal starts to become unsteady, the signal frequency presented a sharp increase from 20 kHz to 500 kHz and remained at high-frequency levels for a short time until the coal ultimately failed. Three specimens with different gas pressure levels displayed similar frequency response patterns, most of which depicted frequency changes from sporadic, blank to shoot up, and from low to high frequency. Our results gave agreement with the corresponding result reported by JI (2015) [46], in which the basic AE characteristics in the low and high frequencies are similar during failure matching the mechanical properties. Thus, signal frequency regularly changes in the cracking process of rock and coal.
Relationship curves of stress-time-AE signal peak frequency.
The process of UCT exhibits certain regularity of AE response in terms AE event count, signal amplitude, and frequency. This finding is followed by a description of CGOB laboratory and a detailed presentation of how to quantitatively represent the MS signal and apply it in predicting hazards. In the subsequent sections, we describe the experiment on CGOB disaster to verify whether the result is the same.
3. CGOB Experiment
In this section, we introduce the CGOB experimental apparatus and scheme. The experiments were completed on December 25, 2011. Experiment phenomena, especially the MS dynamic response of the experiment, are shown in the following description.
3.1. Experimental Apparatus and Scheme
CGOB experiments were carried out by means of a large CGOB physical simulation system, which was developed independently by Chongqing University (see [7] for more details about this system). The simulation system includes three functional parts (see Figure 4): model steel structure (including sample box), loading device and air ejecting device (including gas compressor), and external monitoring equipment (such as gas pressure and temperature). MS equipment from Engineering Seismology Group, Canada Inc. (ESG), was one of the external equipment used to record the MS vibration in real-time dynamic monitoring during the total CGOB process. The parameters of MS monitoring device were set as follows: sampling frequency, 4,000 Hz; continuous-acquisition buffer length, 15 min with a subsequent STA/LTA for picking up and intercepting the events; speed sensor frequency, 50–5 kHz with a sensitivity of 30 V/g; acquisition frequency, 0–2,000 Hz.
Analog simulation system of coal and gas outburst (ASSCGO): (a) steel experimental structure, (b) gas compressor, (c) sample box, and (d) installation of MS sensors.
The total experiment process can be summarized as follows: coal sample production, coal packing and pressing, box sealing and installation, chamber air tightness test, coal inflatable adsorption, MS equipment installation, and press loading. During each experiment, the signals were recorded by the ESG MS equipment in real-time. The recorded events can be used to analyze the response and internal activities in the tested chamber. The rest of the detailed parameters of the process are referred to in the literature [40].
3.2. Experimental Process and Phenomena
The duration of coal outburst simulation parameters and result are summarized in Table 1, which corresponds to the test parameters and results of an experiment completed on December 25, 2011. The monitoring and analysis results indicate that the duration of the outburst process was approximately 15 s. A total of 20.125 kg coal was ejected (77.5 kg coal with 3.6 kg yellow mud infilling was installed). The key time points and parameters, including beginning time, occurring time, gas pressure, and other parameters, are shown in Table 1. Figure 5 reflects outburst situation at different times. Figure 5(d) displays the ejected coal in the outburst stage.
Experimental parameters and result.
Items
Content
Result
Key parameters
1
Gas content
20.125 m3/t
2
Gas pressure
1.41 MPa
3
Description of thrown body
Shape
Cylinder
4
Size
None recorded
Key moments
1
Beginning time
15:30:50
2
Shutting time
16:29:00
3
Preparation time
17:06:54
4
Occurring time of CGOB
19:04:31
Phenomena of CGOB.
Responses during CGOB process were monitored by the MS system. Figure 6 shows the typical MS waveforms during the CGOB experiment. We can clearly observe the variation in MS signal in different moments, including amplitude and waveforms. For CGOB quantitative characterization, three aspects of the issues are discussed and validated in the next step.
Typical waveforms of MS signals in the CGOB laboratory.
4. Time-Frequency Characteristic Analysis of MS Signals
Mining MS signals possess typical nonstationary random and nonlinear properties. Time-frequency analysis is an important new branch of mathematics that has been recently developed. Time-frequency analysis has become of considerable research interest worldwide and thus has been used to analyze nonstationary MS signals.
4.1. Data Preprocessing
During the experiment, the vibration of loading device can trigger the MS monitoring system. The signal during CGOB was different from microcrack signal caused by coal and rock fracture. As the effect of gas, this kind of signal belonged to high-frequency vibration, and its frequency was higher than the original. In this experiment, the maximum frequency was close to 2,000 Hz (sampling frequency of 4,000 dBm/Hz). Therefore, signal component distributed in the range of 0–100 Hz was not useful for analysis in this study and should be deleted before signal analysis.
Based on the comparison of several denoising methods, the notch filter was eventually elected to denoise MS signals. Figures 7(a1) and 7(b1) exhibit variations in the original waveform (Figure 7(a)) and the denoising waveform (Figure 7(b)) versus time, and Figures 7(a2) and 7(b2) show the signal frequency of original and filter signal. Information of the original signal compared with Figure 7 shows that wavelet packet denoising better retained details and useful wavelet denoising. In addition, the processed signal was smoother than before.
MS signal preprocessing results sampling point versus amplitude and frequency versus (a) original MS signal and (b) filtered MS signal.
Original MS signal
Filter MS signal
4.2. MS Dynamic Response
The response frequency of the experiment was recorded at every 4 min during the burst process to highlight the features of the burst. The peak of the signal frequency (Hz) versus time (min) (solid lines) and gas pressure (MPa) versus time (min) (solid lines with marker “o”) are given in Figure 8. From the signal intensity of this figure, the distribution can be roughly divided into four stages, namely, incubation, inspire, happen, and residue stages, which agreed with the findings of Hu (2008) [10]. With the objective to observe the frequency change for a long period, the raw signal data were sliced to data pieces and data lengths of 4,000 sampling points (equal to 1 second) to mitigate the effects caused by data size. Figure 8 shows that the frequency distribution of MS signals scattered around is low in Stage I (event 1, almost less than 200 Hz).
Variation in gas pressure levels and MS response (frequency) during the experiment.
In Stage II (events 2 and 3), a general but slow increase in MS signal frequency of less than 400 Hz was observed with the gas pressure increase. The highest value of frequency in this stage can reach 632.68 Hz. However, Stage III (events 4 and 5) deviated from Stages I and II. Large numbers of high-frequency signals emerged, and the signal frequency markedly changed from low (0–400 Hz) to high frequency (up to 1,600 Hz). In Stage IV (event 6), calm was restored, and the frequency decreased to 50 Hz (power frequency). The variation curve of gas pressure matches the frequency change of MS signal. The CGOB exhibits significant periodic characteristics, and each CGOB stage corresponds to significant MS signal frequency.
This figure provides further evidence that the process of CGOB changes according to different stages along with a clear frequency change trend; that is, the variation trend in signal frequency moves from low to high.
4.3. Frequency Spectrum Analysis
According to Section 4.2, time-frequency analysis, an important method, was used for feature extraction of MS signal [35, 47]. This method can describe the relationship between time and frequency domains and provides conjoint distribution information of both sides (namely, energy density and intensity). Six typical events numbered 1–6 were detected from the CGOB experimental process and analyzed using short-time Fourier transform (see Figure 9) to mirror the microseismic signal’s features in different stages.
Evolution diagrams of MS signal spectrum under the CGOB test.
Figure 9 shows that the distribution of frequency and energy with the variation in time and frequency can reflect certain characteristics. Events 1 and 2 occurred during Stage I. Event 1 mainly concentrated in 800 Hz with the subjective maximum value of 6,000 dBm/Hz. This event was speculated as a mechanical vibration caused by experimental apparatus. Event 2 differed from 1 and 2 focused on 100–200 Hz and belonged to background noise. After a quiet period, event 3 occurred in Stage II with the increase in frequency (from 300 to 800 Hz). CGOBs occurred in Stage III. Event 4 occurred with a series of MS events. Signal dominant frequency and energy reached 1,350 Hz (in the range 480–1,600 Hz) and 1,800 dBm/Hz in this stage. Frequency was then distributed from 800 to 1,380 Hz and was concentrated at 850 Hz (see 5). After coal ejection, the frequency decreased to 50 Hz, and energy was 700 dBm/Hz at maximum (see 6). Comprehensive analysis shows no apparent changes in the energy and frequency of MS signals. In the outburst stage, the complicated frequency components gradually shifted from low to high frequency. This result validates the assumption in Section 2. The remaining problem was how to quantitatively represent the MS signal and apply it in CGOB hazard prediction.
5. Prominent Signs of Seismic Energy Development5.1. Time-Frequency Analysis of Wavelet Packet Method
Wavelet packet transform (WPT), which can capture detailed information with fine resolution in both low- and high-frequency bands, was proposed by Wickerhauser [48] and Coifman [49] as generalization and improvement of wavelet transform. The WPT roughly divides the frequency space of a signal into a low-frequency and a high-frequency component by using a pair of low pass and high pass filters; the same splitting procedure is repeated for all obtained components until a specified level is reached. Through this process, the total frequency space of a signal can be divided into various narrow bands, allowing improved time-frequency localization and frequency components containing low energy to be readily identified [50]. This method exhibits good time-frequency characteristic, multiple scales, and multiresolution for confirming the features of MS signals. By assuming the MS signal is s(t), wavelet packet decomposition model can be expressed as follows:(1)st=∑n=02k-1snjtwhere snj(t) represents the n-th decomposed signal at j layer scale.
Thus, the wavelet packet time-frequency component spectrum WPSnjτ,f is given as follows:(2)WPSnjτ,f=∫Runj,0t-τsnjte-j2πftdt=2-j/2∫Run2-jt-τsnjte-j2πftdt
The time-frequency component spectrum in k layer and j scale can also be expressed as(3)WPSjτ,f=∑n=02k-1WPSnjτ,f
Meanwhile, the time-frequency component spectrum accords with the conservation of energy. Hence, (3) can be expressed in Fourier transform method:(4)∫Rst2dt=∑n=02k-1∫R∫RWPSnjτ,f2dτdf
With the help of the above equations, analysis of MS signals by using wavelet packet method becomes feasible.
5.2. Frequency Band Energy Calculation Method
Signal s(t) (see (1)) can be decomposed into 2j band components [39], forming s0jt,s1jt,…,snjt, and s2j-1jt because the MS signal is decomposed with wavelet packet, and the signal of each frequency band is reconstructed. The energy Ekj of component skjt is given as(5)Ekj=∫Rskjt2dtwhere k∈{0,1,2,…,2j–1}.
The energy percentage of k-th component in the total original signal can be described as(6)Pnj=EkjE=Ekj∑n=02j-1Enj=∫Rskjt2dt∑n=02j-1∫Rsnjt2dt
5.3. Distribution of Wavelet Packet Energy
Wavelet packet analysis was applied to analyze the MS signal to extract the feature of stage signals. The basis function of wavelet package was “sym6”. One channel signal was divided into 5 layers and 32 subbands, and every subband was 62.5 Hz. Table 2 and Figure 10 show six events, and the distribution of band energy variation in time reflects the characteristics of four stages.
Signal energy distribution percentages of frequency bands in different stages.
FB
Energy distribution of Typical Signal /%
FB
Energy distribution of Typical Signal /%
①
②
③
④
⑤
⑥
①
②
③
④
⑤
⑥
S0
1.15
5.37
1.91
3.51
3.28
74.76
S16
0.02
0.01
0.07
0.18
0.55
0.03
S1
16.58
47.77
8.59
2.26
0.62
13.60
S17
0.03
0.01
0.24
0.41
0.31
0.01
S2
3.15
6.54
32.27
1.80
0.56
0.84
S18
0.03
0.00
0.14
1.31
4.10
0.03
S3
14.06
34.68
5.44
1.01
0.34
5.77
S19
0.01
0.01
0.04
0.18
0.86
0.01
S4
0.20
0.06
1.34
1.46
0.94
0.05
S20
0.05
0.02
0.16
1.98
3.11
0.05
S5
0.23
0.21
3.26
7.56
0.95
0.07
S21
0.10
0.02
0.25
2.83
1.26
0.05
S6
0.97
1.89
2.51
0.45
0.76
0.59
S22
0.01
0.00
0.05
0.63
0.63
0.04
S7
0.35
1.21
1.47
0.80
0.61
0.12
S23
0.07
0.00
0.05
0.07
0.85
0.11
S8
0.63
0.13
1.51
2.69
5.49
0.23
S24
0.07
0.02
1.04
7.33
2.21
0.09
S9
0.84
0.07
0.98
0.91
13.27
0.13
S25
1.28
0.09
2.61
6.34
18.03
0.14
S10
32.13
0.13
3.05
7.69
4.03
0.15
S26
3.34
0.06
0.65
10.59
9.39
0.04
S11
5.31
0.69
1.29
5.10
7.45
0.67
S27
2.09
0.08
3.83
7.39
2.45
0.76
S12
0.49
0.11
1.79
1.63
0.57
0.06
S28
0.30
0.02
0.32
5.11
2.24
0.01
S13
1.27
0.06
6.25
2.53
0.76
0.05
S29
0.12
0.03
1.60
3.14
2.51
0.24
S14
10.48
0.64
12.75
1.84
3.93
0.07
S30
1.93
0.02
1.54
10.19
6.03
0.16
S15
1.93
0.04
2.77
0.18
0.84
0.05
S31
0.75
0.00
0.17
0.95
1.04
0.16
Frequency band energy distribution percentages of MS signals at different stages.
Table 2 and Figure 10 show apparent characteristics of frequency and energy variation, and frequency band energy distributions of MS signal largely differed from different stages during gas outburst. In this stage (before outburst), signal energy of number 2 is mainly distributed in the range of 187.5–250 Hz (S3), which accounts for up to 34.68% of the total signal. The same situation occurred in number 1. The distribution trend changes to high frequency along with the process. In the second stage (event 3), frequency composition is dispersed along a wide frequency range from S1 to S14. Signal energy increased to 12.75% in the S14 band. At the same time, bands S29 and S30 experience a modest increase of accumulated energy. The variation of frequency and energy distribution is apparent in events 4 and 5. The signal energy in this stage (outburst occurrence) distributes in a local area, such as S24–S30. In this frequency band range, total energy is 42.76% and 42.86% (average value 42.81%), which are markedly larger than those of events 1–3 (7.01%, (9.13%, 0.32%, and 11.59%, respectively)) and 6 (1.44%). This phenomenon also shows the trend of signal frequency and energy transfer from low to high modes along with the process. Until the end of outburst (event 6, later period), the signal frequency changes and reverts to low frequency (S0).
For clearly showing the frequency band energy distribution of CGOB, MS signal versus time dependence was generated to obtain the characteristics of frequency band energy. Frequency band energy distribution of CGOB MS signals in different stages was clearly given in Figure 10, in which the high amplitude denotes large concentration and high energy of frequency. Signal change during the total process can be characterized as a progressive procedure involving “low frequency distribution-high-frequency transformation-low frequency rotation;” each step corresponded to “before-during-after” stages of the gas outburst. Therefore, the MS precursor characteristics of CGOB process can be excavated effectively by means of this quantitative method.
According to the above analysis, MS signal characteristics can be effectively obtained by the wavelet packet method. This method can provide detailed information by using wavelet packet decomposition. Moreover, wavelet packet method can provide the frequency spectrum and energy characteristics of CGOB. Therefore, finding the precursor information of CGOB is crucial.
The discussion in this paper is only an initial conclusion, and numerous problems have to be solved until the goal of CGOB forecasting is achieved. Identification and intercept of valid waveform from MS real-time data and warning threshold value setting may be key issues in the future.
6. Results and Discussion
Through the above analysis, we obtain preliminary precursor and regularity of CGOB. This work is only the initial research about CGOB disaster prediction, and much remains to be done. The foreboding information is taken as prediction index for CGOB, and a system of CGOB disaster prediction based on MS monitoring is proposed (Figure 11). The system is established based on the fundamental research about mechanism, influence factor of CGOB disaster, and multiparameter monitoring, including MS, EMR, AE, and gas density and pressure. The principle of MS method uses feature mining and waveform recognition between monitoring data and typical MS signal of CGOB disaster. This research content includes several aspects, such as signal recognition and precise positioning, MS event spatial clustering and statistical analysis, signal feature extraction, quantitative expression, and early warning indicators and system.
Prediction of CGOB disaster based on MS monitoring.
7. Conclusions
Rockburst, especially the mixed disaster of coal and gas rockburst, is a challenging problem in coal mines that imposes serious danger to personnel and investment safety. Owing to its complexity associated with stress, gas, and physical and mechanical properties of coal, CGOB damage is difficult to forecast and control. Onsite monitoring and laboratory research indicate that microfracturing is a typical precursor of dynamic disaster, and MS technology has been widely used for mining disaster prediction. For this purpose, some CGOB experiments were carried out by means of a large CGOB physical simulation system with MS monitoring system.
Based on the MS monitoring system experiments and corresponding signal analysis and extractions, the conclusions of CGOB in this study can be drawn as follows. (1) The MS/AE behavior of coal failure experiment was studied first, followed by an explanation of laws: the fracture behavior of coal sample exhibits certain characteristics of MS/AE response in terms of MS/AE event count, signal amplitude, and frequency, and each stage presents its own physical meaning during loading test. The precursor characteristics of AE signal/event can be a foundation of CGOB laboratory studies. (2) The approach presented in this study can capture dynamic response under CGOB testing. Through data preprocessing of MS signals, the method of wavelet packet transform is preferred in time-frequency, nonstationary outburst MS data, and even characteristic signal extraction. The features of each stage (incubation, inspire, happen, and residue stages), including the time-frequency domain, were extracted and quantitatively expressed. (3) Results of the CGOB MS data analysis showed that (1) the process exhibits its own MS response in different stages and represents using the time-frequency analysis method; (2) the frequency and energy density gradually transformed to high-frequency band (600–1,800 Hz), especially in 1,600 Hz; (3) varying characteristics, such as “valley-peaks-valley,” were observed in the total CGOB process.
Therefore, we assumed that the typical waveforms at different times in CGOB also exhibited different characteristics, especially the significant signal frequency changes before burst. If the assumption is reasonable, the quantization of CGOB signal characteristic, which is the basis of our present study, can be achieved.
The approach presented in this study can capture the dynamic response of a rock mass. The results of MS data analysis and interpretation provide dynamic response characteristics, changing characteristics, and a new idea for disaster prediction.
Conflicts of Interest
The authors declare that there are no conflicts of interest regarding the publication of this paper.
Acknowledgments
The authors thank the people who greatly helped in this work and everyone else who is involved with the project, including the stuff who participated in the laboratory. This study was funded by the National Natural Science Foundation of China (nos. 51604115 and 51504096), the Fundamental Research Funds for the Central Universities (3142017002), the Open Foundation of the Hebei State Key Laboratory of Mine Disaster Prevention (KJZH2017K14), the Technology Project and Self-Financing Foundation of Hebei Province (no. 16275512), and the Project of Langfang Key Technology Research and Development Program, China (Grant no. 2016013116).
BeamishB. B.CrosdaleP. J.Instantaneous outbursts in underground coal mines: an overview and association with coal type1998351–4275510.1016/S0166-5162(97)00036-02-s2.0-0032001681YuanL.Control of coal and gas outbursts in Huainan mines in China: A review2016845595672-s2.0-8499155622610.1016/j.jrmge.2016.01.005FanC.LiS.LuoM.DuW.YangZ.Coal and gas outburst dynamic system201727149552-s2.0-8500756626410.1016/j.ijmst.2016.11.003SobczykJ.A comparison of the influence of adsorbed gases on gas stresses leading to coal and gas outburst20141152882942-s2.0-8488128277410.1016/j.fuel.2013.07.016WoldM. B.ConnellL. D.ChoiS. K.The role of spatial variability in coal seam parameters on gas outburst behaviour during coal mining200875111410.1016/j.coal.2008.01.0062-s2.0-44149125502ValliappanS.WohuaZ.Role of gas energy during coal outbursts19994478758952-s2.0-003961390010.1002/(SICI)1097-0207(19990310)44:7<875::AID-NME527>3.0.CO;2-GZbl0947.76087PengS. J.XuJ.YangH. W.LiuD.Experimental study on the influence mechanism of gas seepage on coal and gas outburst disaster20125048168212-s2.0-8485645149210.1016/j.ssci.2011.08.027ZhuW. C.WeiC. H.LiuJ.QuH. Y.ElsworthD.A model of coal-gas interaction under variable temperatures2011862-32132212-s2.0-7995462549710.1016/j.coal.2011.01.011AlexeevA. D.RevvaV. N.AlyshevN. A.ZhitlyonokD. M.True triaxial loading apparatus and its application to coal outburst prediction20045842452502-s2.0-254253454710.1016/j.coal.2003.09.007HuQ.-T.ZhouS.-N.ZhouX.-Q.Mechanical mechanism of coal and gas outburst process20083312136813722-s2.0-58149196246LuC. P.DouL. M.LiuH.LiuH. S.LiuB.DuB. B.Case study on microseismic effect of coal and gas outburst process20125310111010.1016/j.ijrmms.2012.05.0092-s2.0-84861334228PanY.-S.Integrated study on compound dynamic disaster of coal-gas outburst and rockburst20164111051122-s2.0-84959869042AguadoM. B. D.NiciezaC. G.Control and prevention of gas outbursts in coal mines, Riosa-Olloniego coalfield, Spain200769425326610.1016/j.coal.2006.05.0042-s2.0-33846433667TorañoJ.TornoS.AlvarezE.RiesgoP.Application of outburst risk indices in the underground coal mines by sublevel caving2012509410110.1016/j.ijrmms.2012.01.0052-s2.0-84857652705ZhouH.ZhangR.ChengY.DaiH.GeC.ChenJ.Methane and coal exploitation strategy of highly outburst-prone coal seam configurations20152363692-s2.0-8492153519110.1016/j.jngse.2015.01.018LinB.YanF.ZhuC.ZhouY.ZouQ.GuoC.LiuT.Cross-borehole hydraulic slotting technique for preventing and controlling coal and gas outbursts during coal roadway excavation2015265185252-s2.0-8493685797210.1016/j.jngse.2015.06.035FridV.Electromagnetic radiation method for rock and gas outburst forecast1997382971042-s2.0-003142430410.1016/S0926-9851(97)00017-7XueY.GaoF.GaoY.ChengH.LiuY.HouP.TengT.Quantitative evaluation of stress-relief and permeability-increasing effects of overlying coal seams for coal mine methane drainage in Wulan coal mine2016321221372-s2.0-8498229326910.1016/j.jngse.2016.04.029FridV.Rockburst hazard forecast by electromagnetic radiation excited by rock fracture199730422923610.1007/BF010457192-s2.0-0031251171DouL. M.XqH. E.WangE. Y.Electromagnetic emission technique of monitoring rock burst and its application2004294396399SiG.DurucanS.JamnikarS.LazarJ.AbrahamK.KorreA.ShiJ.-Q.ZavšekS.MutkeG.LurkaA.Seismic monitoring and analysis of excessive gas emissions in heterogeneous coal seams201514941542-s2.0-8493795864810.1016/j.coal.2015.06.016ZhaoY.LiuL.PanY.JiaoB.ZhangC.Experiment study on microseismic, charge induction, self-potential and acoustic emission during fracture process of rocks20173611071232-s2.0-85015261848ZouY. H.ZhaoX. S.LiuS.Research on sound transmitted continued prediction technology for coal and gas outburst20053366164GaoB. B.LiH. G.LiH. M.YuanR. F.LiuY. P.Current situation of the study on acoustic emission and microseismic monitoring of coupling dynamic catastrophe for gas-filled coal-rock2014292689697ShadrinA. V.Geophysical criterion of pre-outburst crack propagation in coal beds20165246706822-s2.0-8501606777410.1134/S106273911604107XRibeiroH.CostaL.AlvesL.SantoroP.PicoliS.LenziE.MendesR.Analogies Between the Cracking Noise of Ethanol-Dampened Charcoal and Earthquakes2015115210.1103/PhysRevLett.115.025503DecJ.PietschK.A possibility of seismic identification of gas accumulations in Zechstein carbonate formations of the fore-sudetic monocline2012283931122-s2.0-84866856390DongL.-J.WesselooJ.PotvinY.LiX.-B.Discriminant models of blasts and seismic events in mine seismology20168628229110.1016/j.ijrmms.2016.04.0212-s2.0-84969277573DongL.SunD.LiX.DuK.Theoretical and Experimental Studies of Localization Methodology for AE and Microseismic Sources Without Pre-Measured Wave Velocity in Mines2017516818168282-s2.0-8502850825010.1109/ACCESS.2017.2743115DongL. J.ZouW.LiX. B.ShuW. W.WangZ. W.Collaborative localization method using analytical and iterative solutions for microseismic/acoustic emission sources in the rockmass structure for underground mining2018132AnF. H.ChengY. P.An explanation of large-scale coal and gas outbursts in underground coal mines: The effect of low-permeability zones on abnormally abundant gas201414847514775WangE.-Y.HeX.-Q.LiuZ.-T.LiZ.-H.Study on frequency spectrum characteristics of acoustic emission in coal or rock deformation and fracture20042932892922-s2.0-3442886474LuC.-P.DouL.-M.ZhangN.XueJ.-H.LiuG.-J.Microseismic and acoustic emission effect on gas outburst hazard triggered by shock wave: a case study2014731715173110.1007/s11069-014-1167-72-s2.0-84897362286ZhaoG. Y.DengQ. L.MaJ.Recognition of mine microseismic signals based on FSWT time-frequency analysis2015372306312YanZ.MiyamotoA.JiangZ.Frequency slice wavelet transform for transient vibration response analysis2009235147414892-s2.0-6374911462210.1016/j.ymssp.2009.01.008LiX.ShangX.WangZ.DongL.WengL.Identifying P-phase arrivals with noise: An improved Kurtosis method based on DWT and STA/LTA201613350612-s2.0-8498270626010.1016/j.jappgeo.2016.07.022KimiaefarR.SiahkoohiH. R.HajianA. R.KalhorA.Seismic random noise attenuation using artificial neural network and wavelet packet analysis2016931112-s2.0-84962499039KangY.ZhuW.BaiQ.XieL.TangC.Time-delay estimation of acoustic emission signals of rock using time-frequency energy analysis based on wavelet transform2010295101010162-s2.0-77953776932ZhuQ.JiangF.YuZ.YinY.LuL.Study on energy distribution characters about blasting vibration and rock fracture microseismic signal20123147237302-s2.0-84862330186LeiW.LiS.ShangP.ZhangY.XuJ.WangZ.LiQ.Coal and gas outburst simulation experiment reacted by microseismic monitoring20143111611662-s2.0-84897756545ZhuQ. J.LiQ. S.LiS. Q.Microseismic dynamic response and characteristic analysis of coal and gas outburst experiment2015314723730ZuoJ.PeiJ.-P.LiuJ.-F.PengR.-D.LiY.-C.Investigation on acoustic emission behavior and its time-space evolution mechanism in failure process of coal-rock combined body2011308156415702-s2.0-80052383523ZhaoH.YinG.LiH.WangZ.Analysis on AE characteristic and its confining pressure effect of outburst coal containing gas201336111011072-s2.0-84891640038ButtS. D.FrempongP. K.MukherjeeC.UpshallJ.Characterization of the permeability and acoustic properties of an outburst-prone sandstone20055811122-s2.0-2254446582410.1016/j.jappgeo.2005.04.002YinG.QinH.HuangG.Experimental study of characteristics of seepage and acoustic emission of gas-filled coal under different stress paths2013327131513202-s2.0-84879724702JiH.LuX.Characteristics of acoustic emission and rock fracture precursors of granite under conventional triaxial201534469470210.13722/j.cnki.jrme.2015.04.0052-s2.0-84938386805LuC. P.DouL. M.WuX. R.Frequency spectrum analysis on microseismic monitoring and signal differentiation of rock material2008277772775WickerhauserM. V.Comparison of picture compression methods: wavelet, wavelet packet, and local cosine transform coding1993Academic Press Professional, Inc.679700CoifmanR.MeyerY.WickerhauserV.Signal processing and compression with wavelet packets1994Springer363379MR1173863LawL.KimJ. H.LiewW. Y. H.LeeS.An approach based on wavelet packet decomposition and HilbertHuang transform (WPDHHT) for spindle bearings condition monitoring2012331972112-s2.0-8486565101110.1016/j.ymssp.2012.06.004