# GA-Based Early Warning Method for Rock Burst with Microseismic and Acoustic Emission in Steeply Inclined Coal Seam

Authors: Yarong Xue; Dazhao Song; Zhenlei Li; Jianqiang Chen; Xueqiu He; Shengquan He; Donghui Li
Journal: Shock and Vibration
 (2020)

---

## Abstract

Aiming at problem of low efficacy of early warning of rock burst in coal mine, a multisystem and multiparameter integrated early warning method based on genetic algorithm (GA) is proposed. In this method, firstly, the temporal-spatial-intensity information of energy incubation process of rock burst is deeply mined, and the multidimensional precursory characteristic parameter system of rock burst is constructed. Secondly, the genetic algorithm is used to train the historical monitoring data to obtain the optimal critical value and fitness value of each precursory characteristic parameter, and then the early warning indexWC of each monitoring system is calculated. Finally, the integrated rock burst early warning index IC is obtained by synthesizing the early warning index WC of each system. The value of IC corresponds to the specific rock burst risk level of the mine. This method is applied to Wudong coal mine in Xinjiang, China. Based on the actual situation of the mine, a multidimensional precursory characteristic parameter system of rock burst is constructed, which includes energy deviation (DE), frequency ratio (Fr), frequency deviation (DF), degree of dispersion (DS), and total high value of energy deviation (DH). After analyzing the rock burst danger status and risk level in the monitoring area, the early warning capability of this method is found to reach 0.896. Combining with the specific prevention and control measures corresponding to different rock burst risk levels, it can provide effective guidance for the field work.


---

## Body

## 1. Introduction

Rock burst is one of the main dynamic disasters that affect the underground production of coal mine, which often causes serious consequences such as casualties and major property losses [1–4]. There are 139 mines threatened by rock burst in China, among which 26 are mines with a mining depth of more than one kilometer [5]. With the continuous extension of mining depth and intensity, rock burst disasters continue to escalate, and the number of mines threatened by rock burst is also increasing, which seriously affects the lives and property safety of coal miners.In order to solve this problem, many scholars have made a lot of valuable research from the aspects of rock burst mechanism [6–9], monitoring method [10–12], disaster reduction measures [13, 14], protection means [15], and so on. Among them, the accurate monitoring and early warning method is the premise of the prevention and control of rock burst disaster.At present, the monitoring and early warning methods for rock burst are mainly divided into two categories: traditional rock mechanics methods (e.g., drilling cutting method, coal and rock mass stress measurement, and core discing method) [16–21] and geophysical methods (e.g., electromagnetic radiation (EMR), microseismic (MS), and acoustic emissions (AE)) [11, 22–27]. In recent years, geophysics method has been widely used in underground because of its advantages of noncontact, data continuity, and automatic monitoring. At the same time, in order to reduce the risk of rock burst, coal production enterprises will generally choose the above-mentioned means to improve the accuracy of early warning. However, the monitoring principle and range of different monitoring technology and equipment are quite different, which may lead to the inconsistency of early warning results, thus affecting the judgment of relevant personnel on the actual dangerous situation and the measures taken and reducing the efficiency of the system.Recognizing the problem of insufficient accuracy of single monitoring system, some excellent works have been carried out to establish multisystem monitoring methods to improve the early warning accuracy. He et al. [28] and Li et al. [29] combined MS and EMR and found that the energy and frequency of MS events had a “quiet period” within 1–3 days before the occurrence of rock burst, and the “quiet period” was corresponding to the peak intensity of EMR. Using this law, a coupling evaluation system was established to improve the prediction accuracy of rock burst. He et al. [30] studied the precursory information law of MS and AE monitoring system for rock burst events and proposed characteristic parameters. It found that it has obvious response law to the rock burst events of steeply inclined and extremely thick coal seam, which provided a solid foundation for the safe mining of such coal seam. Qin et al. [31] proposed an improved drill cutting method to make up for the defects of traditional drilling cutting method where the depth cannot reach the stress concentration area and cannot accurately reflect the internal stress variations of coal body. After the combination of improved drilling cutting method and the MS monitoring system, it overcomes the limitations of a single method and has achieved successful application in the underground. Lou et al. [32] carried out the uniaxial compression tests to the burst-prone coal and rock samples, and the results show that AE has a high correlation with the EME spectrum, and EME is usually accompanied by stress drop and high-intensity AE. This feature can be used to combine the two monitoring methods to improve the rock burst early warning accuracy of the working surface range. Tan et al. [33] used the energy release before the rock burst as a precursor to establish a multi-index risk warning system based on MS and EMR monitoring methods. Combined with the drilling bits method to determine the burst region, the targeted stress relief measures were taken to ensure the safe and efficient production of the mine. Xia et al. [34] extracted the energy deviation value and other parameters from the MS and AE monitoring data. After the capability test of each parameter with R-value scoring method, it is found that its early warning capability is good, and the prediction of rock burst can be more accurate after multisystem integration. Many applications of the multisystem monitoring methods [35–38] can also be found in various aspects of mining engineering.All of the above researches have made beneficial exploration of the methods and indexes of monitoring and early warning of rock burst and achieved certain results in the field practice. However, in practice, it relies on artificially determining the influencing factors of coal rock stress and their weights, which has a certain degree of subjectivity. At the same time, the historical data of each monitoring system contains insufficient information mining and low data analysis efficiency, and the timeliness and accuracy of early warning have room to improve.In recent years, the intelligent and precise mining of coal mine has become a hot issue. Using machine learning method to analyze and mine a large amount of historical data of various monitoring systems has become a new way to realize the underground intelligent monitoring and early warning. Faradonbeh et al. [39] used two robust clustering techniques of self-organizing map (SOM) and fuzzy c-mean (FCM) to find the latent relationship between the rock burst-related parameters based on the compiled data samples from deep underground projects. And the result shows that the maximum tangential stress around the openings has a crucial role in rock burst clustering and has the most influence on the occurrence of strong and moderate rock burst types, which can be used as a method to predict the rock burst. Xue et al. [40] established a PSO-ELM rock burst prediction model using a database of 344 collected worldwide rock burst cases, which greatly increase the accuracy of rock burst prediction. Wu et al. [41] extracted six quantitatively characteristic parameters including the maximum tangential stress, the uniaxial compressive strength, and the uniaxial tensile strength from historical rock burst cases and established a rock burst prediction probability model based on the Least Squares Support Vector Machine (LSSVM) which was optimized by particle swarm optimization algorithm. Li et al. [42] established a rock burst risk evaluation model based on particle swarm optimization and BP neural network, which improved the accuracy of the standard BP model by 15% and has higher adaptability. Zhou et al. [43] applied the genetic algorithm and particle swarm optimization algorithm to the parameter optimization of support vector machine. The accuracy of the prediction results is higher than that of the measured samples, and it becomes an effective tool for monitoring and early warning of rock burst. Adoko et al. [44] used a large amount of on-site monitoring data to establish a rock burst prediction model based on fuzzy inference system (FIS) and adaptive neurofuzzy inference systems (ANFIS). The practical results show that the model has a high prediction accuracy for rock burst. Lin et al. [45] and Zhou et al. [46] used cloud model to classify and predict rock burst, using uniaxial compressive strength, tensile strength, tangential stress, rock brittleness coefficient, stress coefficient, and elastic energy index to establish the evaluation index system. The results show that the cloud model has higher discriminatory ability than Bayes, KNN, and RF methods, which can predict the rock burst event more accurately and reasonably. Li et al. [47] proposed a new Multiresolution Feature Fusion SVM (MRFF-SVM) recognition approach to compute a comprehensive feature vector for coal rock burst risk recognition and forecasting, and the risk of rock burst can be effectively monitored to reduce casualties and property losses. The above researches have made positive contribution to improving the early warning accuracy of rock burst and effectively improved the efficiency of data analysis.Genetic algorithm is an adaptive global optimization probability search algorithm, which has strong robustness. Research shows that it has good application efficiency in dealing with multipeak problems, which means that it is suitable for solving complex mine monitoring data [48–51]. Therefore, we apply it to the optimization of rock burst danger critical value. Based on this, this paper puts forward a new multiparameter integrated early warning method for rock burst based on genetic algorithm and combines it with the measured data of No. B3 + 6 mining face of + 450 mining level in Wudong coal mine (WCM), Xinjiang, China, carries out simulation research, and verifies its feasibility in MATLAB R2014a environment. The research results play an important role in the construction of continuous, real-time, and accurate monitoring and early warning system of rock burst and can effectively guide the on-site production and disaster relief work.


## 2. Methods

### 2.1. Genetic Algorithm Optimization of Rock Burst Danger Critical Value

#### 2.1.1. Principle of Genetic Algorithm

Genetic algorithm (GA) is a kind of directed process search optimal solution algorithm that mimics the evolutionary mechanism of natural organisms in the target space. It can optimize the multiparameter or multicombination simultaneously. The essence is to use the mechanism of “survival of the fittest” in nature and finally select the samples with high environmental adaptability through iterative methods by genetic manipulation such as selection, crossover, and mutation [52, 53]. The calculation process of GA is shown in Figure 1, and the main parameter values are shown in Table 1.Figure 1
Genetic algorithm calculation process.Table 1
Main parameter values of GA.
OptionsDescriptionsValuesPopulationSizeSize of the population30CrossoverFractionThe fraction of the population at the next generation, not including elite children0.4StallGenLimitPositive integer. The algorithm stops if the average relative change in the best fitness function value over StallGenLimit generations is less than or equal to TolFun (default value is 1E − 06)20MigrationFractionScalar between 0 and 1 specifying the fraction of individuals in each subpopulation that migrates to a different subpopulation0.2EliteCountPositive integer specifying how many individuals in the current generation are guaranteed to survive to the next generation6The indexes of the commonly used rock burst monitoring system, such as MS and AE, are discontinuous, multipeak, and noisy, which makes the critical value easy to fall into the local optimal solution in the process of optimization, and the search for the optimal critical value is closely related to the accurate early warning of rock burst. The GA naturally has excellent ability to find the global optimal solution, which makes it converge to the global optimal solution with a large probability in this complex environment, and the application performance is good. Specific steps are as follows:(a)
Combine the precursory characteristic parameters of each monitoring system into pairs and input them into the GA model as training samples.(b)
Automatically generate initial population at random.(c)
Calculate the fitness function of the initial population. The greater the fitness value is, the faster the fitness function converges.(d)
Determine whether the termination condition is satisfied, and if so, terminate the operation and output an optimal critical value of each characteristic parameter and a corresponding fitness value, and if not, perform step (e).(e)
Perform selection, crossover, and mutation operators, retain the individuals with high fitness value, eliminate the individuals with low fitness, and generate new populations. Return to step (d) to continue the operation until the termination condition is met.


#### 2.1.2. Fitness Function

The fitness function of genetic algorithm has no fixed requirement for form, so it has good adaptability and scalability. In this paper, theR-value scoring method of earthquake prediction capability is used as the fitness function of GA. The R-value scoring is a method proposed in the 1970s to evaluate the capability of earthquake prediction and has been widely used for decades with good application effect [36, 54, 55]. The idea is clear and easy to operate; neither false warnings nor missing warnings are encouraged. It can be calculated as follows:(1)R=a−d=n11N1−n01N0,where n11 represents the accurate prediction number of large energy tremors and rock bursts, N1 represents the number of large energy tremors and rock bursts, n01 represents the number of false warnings, N0 represents the number of MS events that are neither rock burst nor large energy tremor, a=n11/N1 represents the early warning accuracy, and d=n01/N0 represents the false warning rate. The closer the R-value is to 1, the higher the accuracy of early warning is and the lower the rate of false or missing warning is.


### 2.2. Extraction of Precursory Characteristic Parameters of Rock Burst

Accurate and clear physical meaning of the precursory characteristic parameters of rock burst is the basis of the effective use of GA for early warning. At present, the monitoring information of each monitoring system is not fully utilized. Most of the mines still use a single critical value of energy or frequency to early warning rock burst [36], which does not consider the large amount of internal information such as amplitude and growth rate in the monitoring data, resulting in the low accuracy of early warning. Therefore, it is necessary to establish a precursory characteristic parameter system that can reflect the evolution features of underground tremor.


#### 2.2.1. Temporal Characteristic Parameter

In general, the commonly used geophysical monitoring methods of rock burst include MS, AE, and EMR, which automatically record the variation of stress activities in coal and rock mass in time sequence. For example, the MS monitoring system will automatically collect and filter the tremor signal and calculate its occurrence time, three-dimensional coordinate position, energy, and so on. The greater the energy of the tremor occurred, the larger the fracture inside the coal rock mass is, which may induce a dynamic disaster, and the risk is evaluated accordingly.However, the raw data of each monitoring system is mostly multipeak and noisy, and a single energy or frequency critical value is likely to cause false or missing warning. Therefore, based on the real-time data of monitoring system, the characteristic parameters including information such as growth rate and amplitude should be extracted to improve the efficiency of information use.


#### 2.2.2. Spatial Characteristic Parameter

The tremor may occur in the far-field of rock stratum, near-field of surrounding rock, or local dangerous point of roadway, and it is in the dynamic change with the coal mining [56]. The smaller the positional dispersion of the tremor, the more concentrated the stress inside the rock mass. And the probability of dynamic disaster in the stress concentration area is significantly higher than that in other areas, so the description of the location information of tremors is also very important. For example, the MS monitoring system can be used to record the three-dimensional coordinates of the MS events in order to establish spatial characteristic parameters.


#### 2.2.3. Intensity Characteristic Parameter

There are two kinds of tremors caused by mining activities in underground mine: one is low-frequency, high-energy, which occurs in the high-level rock stratum and belongs to the category of MS; the other is high-frequency, low-energy, which occurs in the low-level surrounding rock and belongs to the category of AE. Studies have shown that the intensity of such tremors can reflect the change trend of the internal stress and rock burst risk of rock mass, and the commonly used early warning indicators includeb value, A(b) value, and P(b) value [57].In summary, the multidimensional precursory characteristic parameters with clear physical significance and sensitivity extracted from a large amount of historical data of each monitoring system provide effective supports for the accurate early warning of rock burst danger.


### 2.3. Early Warning Classification Criteria and Construction of Integrated Early Warning Index

#### 2.3.1. Early Warning Index of Single Monitoring System,WC

While using GA to optimize the critical value of precursory characteristic parameter, firstly, combine the characteristic parameters from different dimensions of each system in pairs and then input different pairs of characteristic parameters into the GA to train the optimal critical value recorded asV1 and V2, and at the same time, the fitness value Rj corresponding to V1 and V2 calculated by the fitness function is obtained, where j represents the different characteristic parameter combination.Compare the fitness values of different characteristic parameter combinations, take the largest one or two combinations as the characteristic parameters of the monitoring system, record thei-th of them as ai and bi, and then compare them, respectively, with V1 and V2 to obtain the warning degree Wi. The calculation method is shown as follows:(2)Wi=0,ai<V1&bi<V2,1,ai<V1orbi<V2,2,ai>V1&bi>V2.Then, the warning degreeWi at the i-th moment and the fitness value Rj are substituted into equation (3) to obtain the early warning index of single monitoring system WCi:(3)WCi=∑WimaxWi×Rj∑Rj,where Wi represents the warning degree of different characteristic parameter combinations of a single system or different characteristic parameters of different sensors of a single system; Rj represents the fitness value of different characteristic parameter combinations of a single system or different characteristic parameters of different sensors of a single system.


#### 2.3.2. Integrated Early Warning Index of the Multiple Monitoring System,IC

In order to synthesize the early warning results of the precursory characteristic parameters of each system and form the real-time rock burst risk early warning index with unified early warning criteria and critical values, combined with the synthetical anomaly index method [58, 59], the integrated early warning index of rock burst risk is constructed and can be calculated as follows:(4)ICi=∑e−e1−WCie−1×RCi∑RCi,where WCi represents the early warning index of a monitoring system and RCi represents the early warning capability of the corresponding monitoring system calculated with R-value scoring method.The integrated early warning indexIC can be obtained by the monitoring data of a certain period of time; in order to timely and clearly determine the current risk of rock burst, IC is divided into four levels, corresponding to no rock burst risk, weak rock burst risk, medium rock burst risk, and strong rock burst risk, respectively. As shown in Table 2, different prevention and control measures can be adopted for different risk levels.Table 2
Rock burst hazard classification.
Integrated early warning index of rock burst,ICLevels of rock burst riskState of rock burst risk0 ≤IC < 0.25INo rock burst risk0.25 ≤IC < 0.5IIWeak rock burst risk0.5 ≤IC < 0.75IIIMedium rock burst risk0.75 ≤IC ≤ 1IVStrong rock burst riskA multisystem and multiparameter integrated early warning method for rock burst danger is established as above. This method can adjust the weight coefficient of each monitoring system in real time according to the calculation result. The higher the early warning capability, the greater the weight of the system, which can effectively improve the accuracy of early warning. The flowchart is shown in Figure2.Figure 2
The flowchart of the multiparameter integrated early warning method for rock burst danger.


## 2.1. Genetic Algorithm Optimization of Rock Burst Danger Critical Value

### 2.1.1. Principle of Genetic Algorithm

Genetic algorithm (GA) is a kind of directed process search optimal solution algorithm that mimics the evolutionary mechanism of natural organisms in the target space. It can optimize the multiparameter or multicombination simultaneously. The essence is to use the mechanism of “survival of the fittest” in nature and finally select the samples with high environmental adaptability through iterative methods by genetic manipulation such as selection, crossover, and mutation [52, 53]. The calculation process of GA is shown in Figure 1, and the main parameter values are shown in Table 1.Figure 1
Genetic algorithm calculation process.Table 1
Main parameter values of GA.
OptionsDescriptionsValuesPopulationSizeSize of the population30CrossoverFractionThe fraction of the population at the next generation, not including elite children0.4StallGenLimitPositive integer. The algorithm stops if the average relative change in the best fitness function value over StallGenLimit generations is less than or equal to TolFun (default value is 1E − 06)20MigrationFractionScalar between 0 and 1 specifying the fraction of individuals in each subpopulation that migrates to a different subpopulation0.2EliteCountPositive integer specifying how many individuals in the current generation are guaranteed to survive to the next generation6The indexes of the commonly used rock burst monitoring system, such as MS and AE, are discontinuous, multipeak, and noisy, which makes the critical value easy to fall into the local optimal solution in the process of optimization, and the search for the optimal critical value is closely related to the accurate early warning of rock burst. The GA naturally has excellent ability to find the global optimal solution, which makes it converge to the global optimal solution with a large probability in this complex environment, and the application performance is good. Specific steps are as follows:(a)
Combine the precursory characteristic parameters of each monitoring system into pairs and input them into the GA model as training samples.(b)
Automatically generate initial population at random.(c)
Calculate the fitness function of the initial population. The greater the fitness value is, the faster the fitness function converges.(d)
Determine whether the termination condition is satisfied, and if so, terminate the operation and output an optimal critical value of each characteristic parameter and a corresponding fitness value, and if not, perform step (e).(e)
Perform selection, crossover, and mutation operators, retain the individuals with high fitness value, eliminate the individuals with low fitness, and generate new populations. Return to step (d) to continue the operation until the termination condition is met.


### 2.1.2. Fitness Function

The fitness function of genetic algorithm has no fixed requirement for form, so it has good adaptability and scalability. In this paper, theR-value scoring method of earthquake prediction capability is used as the fitness function of GA. The R-value scoring is a method proposed in the 1970s to evaluate the capability of earthquake prediction and has been widely used for decades with good application effect [36, 54, 55]. The idea is clear and easy to operate; neither false warnings nor missing warnings are encouraged. It can be calculated as follows:(1)R=a−d=n11N1−n01N0,where n11 represents the accurate prediction number of large energy tremors and rock bursts, N1 represents the number of large energy tremors and rock bursts, n01 represents the number of false warnings, N0 represents the number of MS events that are neither rock burst nor large energy tremor, a=n11/N1 represents the early warning accuracy, and d=n01/N0 represents the false warning rate. The closer the R-value is to 1, the higher the accuracy of early warning is and the lower the rate of false or missing warning is.


## 2.1.1. Principle of Genetic Algorithm

Genetic algorithm (GA) is a kind of directed process search optimal solution algorithm that mimics the evolutionary mechanism of natural organisms in the target space. It can optimize the multiparameter or multicombination simultaneously. The essence is to use the mechanism of “survival of the fittest” in nature and finally select the samples with high environmental adaptability through iterative methods by genetic manipulation such as selection, crossover, and mutation [52, 53]. The calculation process of GA is shown in Figure 1, and the main parameter values are shown in Table 1.Figure 1
Genetic algorithm calculation process.Table 1
Main parameter values of GA.
OptionsDescriptionsValuesPopulationSizeSize of the population30CrossoverFractionThe fraction of the population at the next generation, not including elite children0.4StallGenLimitPositive integer. The algorithm stops if the average relative change in the best fitness function value over StallGenLimit generations is less than or equal to TolFun (default value is 1E − 06)20MigrationFractionScalar between 0 and 1 specifying the fraction of individuals in each subpopulation that migrates to a different subpopulation0.2EliteCountPositive integer specifying how many individuals in the current generation are guaranteed to survive to the next generation6The indexes of the commonly used rock burst monitoring system, such as MS and AE, are discontinuous, multipeak, and noisy, which makes the critical value easy to fall into the local optimal solution in the process of optimization, and the search for the optimal critical value is closely related to the accurate early warning of rock burst. The GA naturally has excellent ability to find the global optimal solution, which makes it converge to the global optimal solution with a large probability in this complex environment, and the application performance is good. Specific steps are as follows:(a)
Combine the precursory characteristic parameters of each monitoring system into pairs and input them into the GA model as training samples.(b)
Automatically generate initial population at random.(c)
Calculate the fitness function of the initial population. The greater the fitness value is, the faster the fitness function converges.(d)
Determine whether the termination condition is satisfied, and if so, terminate the operation and output an optimal critical value of each characteristic parameter and a corresponding fitness value, and if not, perform step (e).(e)
Perform selection, crossover, and mutation operators, retain the individuals with high fitness value, eliminate the individuals with low fitness, and generate new populations. Return to step (d) to continue the operation until the termination condition is met.


## 2.1.2. Fitness Function

The fitness function of genetic algorithm has no fixed requirement for form, so it has good adaptability and scalability. In this paper, theR-value scoring method of earthquake prediction capability is used as the fitness function of GA. The R-value scoring is a method proposed in the 1970s to evaluate the capability of earthquake prediction and has been widely used for decades with good application effect [36, 54, 55]. The idea is clear and easy to operate; neither false warnings nor missing warnings are encouraged. It can be calculated as follows:(1)R=a−d=n11N1−n01N0,where n11 represents the accurate prediction number of large energy tremors and rock bursts, N1 represents the number of large energy tremors and rock bursts, n01 represents the number of false warnings, N0 represents the number of MS events that are neither rock burst nor large energy tremor, a=n11/N1 represents the early warning accuracy, and d=n01/N0 represents the false warning rate. The closer the R-value is to 1, the higher the accuracy of early warning is and the lower the rate of false or missing warning is.


## 2.2. Extraction of Precursory Characteristic Parameters of Rock Burst

Accurate and clear physical meaning of the precursory characteristic parameters of rock burst is the basis of the effective use of GA for early warning. At present, the monitoring information of each monitoring system is not fully utilized. Most of the mines still use a single critical value of energy or frequency to early warning rock burst [36], which does not consider the large amount of internal information such as amplitude and growth rate in the monitoring data, resulting in the low accuracy of early warning. Therefore, it is necessary to establish a precursory characteristic parameter system that can reflect the evolution features of underground tremor.


### 2.2.1. Temporal Characteristic Parameter

In general, the commonly used geophysical monitoring methods of rock burst include MS, AE, and EMR, which automatically record the variation of stress activities in coal and rock mass in time sequence. For example, the MS monitoring system will automatically collect and filter the tremor signal and calculate its occurrence time, three-dimensional coordinate position, energy, and so on. The greater the energy of the tremor occurred, the larger the fracture inside the coal rock mass is, which may induce a dynamic disaster, and the risk is evaluated accordingly.However, the raw data of each monitoring system is mostly multipeak and noisy, and a single energy or frequency critical value is likely to cause false or missing warning. Therefore, based on the real-time data of monitoring system, the characteristic parameters including information such as growth rate and amplitude should be extracted to improve the efficiency of information use.


### 2.2.2. Spatial Characteristic Parameter

The tremor may occur in the far-field of rock stratum, near-field of surrounding rock, or local dangerous point of roadway, and it is in the dynamic change with the coal mining [56]. The smaller the positional dispersion of the tremor, the more concentrated the stress inside the rock mass. And the probability of dynamic disaster in the stress concentration area is significantly higher than that in other areas, so the description of the location information of tremors is also very important. For example, the MS monitoring system can be used to record the three-dimensional coordinates of the MS events in order to establish spatial characteristic parameters.


### 2.2.3. Intensity Characteristic Parameter

There are two kinds of tremors caused by mining activities in underground mine: one is low-frequency, high-energy, which occurs in the high-level rock stratum and belongs to the category of MS; the other is high-frequency, low-energy, which occurs in the low-level surrounding rock and belongs to the category of AE. Studies have shown that the intensity of such tremors can reflect the change trend of the internal stress and rock burst risk of rock mass, and the commonly used early warning indicators includeb value, A(b) value, and P(b) value [57].In summary, the multidimensional precursory characteristic parameters with clear physical significance and sensitivity extracted from a large amount of historical data of each monitoring system provide effective supports for the accurate early warning of rock burst danger.


## 2.2.1. Temporal Characteristic Parameter

In general, the commonly used geophysical monitoring methods of rock burst include MS, AE, and EMR, which automatically record the variation of stress activities in coal and rock mass in time sequence. For example, the MS monitoring system will automatically collect and filter the tremor signal and calculate its occurrence time, three-dimensional coordinate position, energy, and so on. The greater the energy of the tremor occurred, the larger the fracture inside the coal rock mass is, which may induce a dynamic disaster, and the risk is evaluated accordingly.However, the raw data of each monitoring system is mostly multipeak and noisy, and a single energy or frequency critical value is likely to cause false or missing warning. Therefore, based on the real-time data of monitoring system, the characteristic parameters including information such as growth rate and amplitude should be extracted to improve the efficiency of information use.


## 2.2.2. Spatial Characteristic Parameter

The tremor may occur in the far-field of rock stratum, near-field of surrounding rock, or local dangerous point of roadway, and it is in the dynamic change with the coal mining [56]. The smaller the positional dispersion of the tremor, the more concentrated the stress inside the rock mass. And the probability of dynamic disaster in the stress concentration area is significantly higher than that in other areas, so the description of the location information of tremors is also very important. For example, the MS monitoring system can be used to record the three-dimensional coordinates of the MS events in order to establish spatial characteristic parameters.


## 2.2.3. Intensity Characteristic Parameter

There are two kinds of tremors caused by mining activities in underground mine: one is low-frequency, high-energy, which occurs in the high-level rock stratum and belongs to the category of MS; the other is high-frequency, low-energy, which occurs in the low-level surrounding rock and belongs to the category of AE. Studies have shown that the intensity of such tremors can reflect the change trend of the internal stress and rock burst risk of rock mass, and the commonly used early warning indicators includeb value, A(b) value, and P(b) value [57].In summary, the multidimensional precursory characteristic parameters with clear physical significance and sensitivity extracted from a large amount of historical data of each monitoring system provide effective supports for the accurate early warning of rock burst danger.


## 2.3. Early Warning Classification Criteria and Construction of Integrated Early Warning Index

### 2.3.1. Early Warning Index of Single Monitoring System,WC

While using GA to optimize the critical value of precursory characteristic parameter, firstly, combine the characteristic parameters from different dimensions of each system in pairs and then input different pairs of characteristic parameters into the GA to train the optimal critical value recorded asV1 and V2, and at the same time, the fitness value Rj corresponding to V1 and V2 calculated by the fitness function is obtained, where j represents the different characteristic parameter combination.Compare the fitness values of different characteristic parameter combinations, take the largest one or two combinations as the characteristic parameters of the monitoring system, record thei-th of them as ai and bi, and then compare them, respectively, with V1 and V2 to obtain the warning degree Wi. The calculation method is shown as follows:(2)Wi=0,ai<V1&bi<V2,1,ai<V1orbi<V2,2,ai>V1&bi>V2.Then, the warning degreeWi at the i-th moment and the fitness value Rj are substituted into equation (3) to obtain the early warning index of single monitoring system WCi:(3)WCi=∑WimaxWi×Rj∑Rj,where Wi represents the warning degree of different characteristic parameter combinations of a single system or different characteristic parameters of different sensors of a single system; Rj represents the fitness value of different characteristic parameter combinations of a single system or different characteristic parameters of different sensors of a single system.


### 2.3.2. Integrated Early Warning Index of the Multiple Monitoring System,IC

In order to synthesize the early warning results of the precursory characteristic parameters of each system and form the real-time rock burst risk early warning index with unified early warning criteria and critical values, combined with the synthetical anomaly index method [58, 59], the integrated early warning index of rock burst risk is constructed and can be calculated as follows:(4)ICi=∑e−e1−WCie−1×RCi∑RCi,where WCi represents the early warning index of a monitoring system and RCi represents the early warning capability of the corresponding monitoring system calculated with R-value scoring method.The integrated early warning indexIC can be obtained by the monitoring data of a certain period of time; in order to timely and clearly determine the current risk of rock burst, IC is divided into four levels, corresponding to no rock burst risk, weak rock burst risk, medium rock burst risk, and strong rock burst risk, respectively. As shown in Table 2, different prevention and control measures can be adopted for different risk levels.Table 2
Rock burst hazard classification.
Integrated early warning index of rock burst,ICLevels of rock burst riskState of rock burst risk0 ≤IC < 0.25INo rock burst risk0.25 ≤IC < 0.5IIWeak rock burst risk0.5 ≤IC < 0.75IIIMedium rock burst risk0.75 ≤IC ≤ 1IVStrong rock burst riskA multisystem and multiparameter integrated early warning method for rock burst danger is established as above. This method can adjust the weight coefficient of each monitoring system in real time according to the calculation result. The higher the early warning capability, the greater the weight of the system, which can effectively improve the accuracy of early warning. The flowchart is shown in Figure2.Figure 2
The flowchart of the multiparameter integrated early warning method for rock burst danger.


## 2.3.1. Early Warning Index of Single Monitoring System,WC

While using GA to optimize the critical value of precursory characteristic parameter, firstly, combine the characteristic parameters from different dimensions of each system in pairs and then input different pairs of characteristic parameters into the GA to train the optimal critical value recorded asV1 and V2, and at the same time, the fitness value Rj corresponding to V1 and V2 calculated by the fitness function is obtained, where j represents the different characteristic parameter combination.Compare the fitness values of different characteristic parameter combinations, take the largest one or two combinations as the characteristic parameters of the monitoring system, record thei-th of them as ai and bi, and then compare them, respectively, with V1 and V2 to obtain the warning degree Wi. The calculation method is shown as follows:(2)Wi=0,ai<V1&bi<V2,1,ai<V1orbi<V2,2,ai>V1&bi>V2.Then, the warning degreeWi at the i-th moment and the fitness value Rj are substituted into equation (3) to obtain the early warning index of single monitoring system WCi:(3)WCi=∑WimaxWi×Rj∑Rj,where Wi represents the warning degree of different characteristic parameter combinations of a single system or different characteristic parameters of different sensors of a single system; Rj represents the fitness value of different characteristic parameter combinations of a single system or different characteristic parameters of different sensors of a single system.


## 2.3.2. Integrated Early Warning Index of the Multiple Monitoring System,IC

In order to synthesize the early warning results of the precursory characteristic parameters of each system and form the real-time rock burst risk early warning index with unified early warning criteria and critical values, combined with the synthetical anomaly index method [58, 59], the integrated early warning index of rock burst risk is constructed and can be calculated as follows:(4)ICi=∑e−e1−WCie−1×RCi∑RCi,where WCi represents the early warning index of a monitoring system and RCi represents the early warning capability of the corresponding monitoring system calculated with R-value scoring method.The integrated early warning indexIC can be obtained by the monitoring data of a certain period of time; in order to timely and clearly determine the current risk of rock burst, IC is divided into four levels, corresponding to no rock burst risk, weak rock burst risk, medium rock burst risk, and strong rock burst risk, respectively. As shown in Table 2, different prevention and control measures can be adopted for different risk levels.Table 2
Rock burst hazard classification.
Integrated early warning index of rock burst,ICLevels of rock burst riskState of rock burst risk0 ≤IC < 0.25INo rock burst risk0.25 ≤IC < 0.5IIWeak rock burst risk0.5 ≤IC < 0.75IIIMedium rock burst risk0.75 ≤IC ≤ 1IVStrong rock burst riskA multisystem and multiparameter integrated early warning method for rock burst danger is established as above. This method can adjust the weight coefficient of each monitoring system in real time according to the calculation result. The higher the early warning capability, the greater the weight of the system, which can effectively improve the accuracy of early warning. The flowchart is shown in Figure2.Figure 2
The flowchart of the multiparameter integrated early warning method for rock burst danger.


## 3. Case Study

### 3.1. Basic Conditions of the No.B3 + 6 Working Face, WCM

WCM is located in Urumqi, Xinjiang, China. The overall geological section of the north and south mining area is shown in Figure3, and up to now, the current mining situation of WCM is at No. B3 + 6 + 425 m working face. Both the No. B3 + 6 and No. B1 + 2 are classified as steeply inclined and extremely thick coal seam, and there is a rock pillar with 70–101 m thickness between them. The average inclination of the coal seam is 87°. The horizontal section full-mechanized top-coal caving mining technology is being adopted for mining, the ratio of mining to caving height is 1 : 7, and the horizontal section height is 25 m. Due to the special occurrence conditions of the coal seam, with the increase of the mining depth, a towering rock pillar is formed between the two coal seams. The rock pillar is affected by the horizontal tectonic stress and self-heavy stress, which causes the tilting effect of the intersection between the lower side of the rock column and the coal seam. Under the disturbance of mining, the elastic energy of rock mass accumulation is released through the coal body, which is easy to produce different degrees of dynamic disasters. There are six times of rock burst in the process of No. B3 + 6 coal seam mining.Figure 3
The overall geological section of WCM.At present, the mine has installed various monitoring systems such as ARAMIS M/E MS monitoring system and KJ623 AE monitoring system. The ARAMIS M/E is a particular MS monitoring system, which can automatically collect and filter seismic signals, and then calculates the occurrence time, released energy, and three-dimensional coordinates of both tremors and rock bursts. The south mining area of the WCM was surrounded by a total of 12 sensors, each sensor covered a range of 2 km, and the layout of the sensors until 26 April 2017 when the last rock burst occurred is shown in Figure4 (please refer to Li et al. 2014 for more details). The KJ623 AE monitoring system contains a total of 11 sensors that was installed on the coal and rock mass to collect AE events within a frequency range of 300–2000 Hz. The sensitivity of the AE transducer was 65 ± 3 dB (V/m/s), the bolt was used as the AE conductor, and the sensor was installed on the bolt. The sensors are located at 30 m in front of the working face, and the spacing between each sensor is 30–40 m. As the working face moves, the sensor layout moves accordingly.Figure 4
Network layout of MS and AE monitoring system in No.B3 + 6 coal seam +450 horizontal working face of WCM during the research period.According to the actual situation of WCM, the large energy mine tremors are defined as MS event with energy exceeding 1 × 105J. Taking the No. B3 + 6 coal seam + 450 horizontal working face of the coal mine as the research object, the MS and AE monitoring data of 2016/06/06–2017/12/28 time period are selected (see Tables S1-S2 in the Supplementary Material for complete field monitoring data of MS and AE). During this period, there are three rock burst events and ten large energy mine tremors in total, and the distribution of rock burst events is shown in Figure 5. Table 3 shows the specific manifestation records of these three rock burst events.Figure 5
The distribution map of MS events with energy exceeding 1 × 103J in No. B3 + 6 coal seam + 450 horizontal working face WCM during 2016/06/06–2017/12/28.Table 3
On-site rock burst manifestation records.
DateLocation of working face (m)Location of MS sourceMS energy (J)Damage form and degree2016/11/242026Located in theB2-B3 rock wall, 12 m away from the B3 bottom plate9.5 × 106Upper wall moved 0.3 m, lower wall moved 0.5 m, floor heaved by 0.4 m, damage to hydraulic props2017/02/011824Located at 9 m from theB6 roof slate2.1 × 108Upper wall moved 0.3–0.7 m, lower wall moved 0.5–1.3 m, damage to belt conveyor2017/04/261568Located in theB2-B3 rock pillar, 35 m away from the B3 bottom plate2.2 × 106Lower wall moved 0.3–0.5 m, floor heaved by 0.2–0.3 m, damage to hydraulic props and belt conveyor, reversed loader change direction


### 3.2. Application of the Integrated Early Warning Method in WCM

#### 3.2.1. Construction of Precursory Characteristic Parameter System

(1) Temporal Characteristic Parameter(a)
Frequency deviation (DF)The information of the mine tremor monitored by the MS system can reflect the internal stress of the rock mass. Prior to the occurrence of rock burst, the source of the MS increased and clustered anomaly compared to the previous normal mining period, which shows that the stress in the study area is close to the strength limit of coal and rock mass. Therefore, the frequency deviation value (DF) of a certain time window and sliding step length in the monitoring area is calculated. It can reflect the trend of the MS events over a period of time; the denser the MS tremor is in time sequence, the larger the value is. It can be calculated as follows:(5)DF=PF−PF¯PF¯,wherePF represents the number of MS tremors within the first 24 hours and PF¯ represents the mean value of daily tremors during the period from the time of the last large energy mine tremor to the current time.(b)
Frequency ratio (Fr)The higher the frequency of MS tremor is, the more frequent the stress activities in the rock mass are and the more likely the rock burst is to be induced; otherwise, the less likely it is. Therefore, the frequency ratio (Fr) of the characteristic parameter is defined, and the calculation method is as follows:(6)Fr=Ft1¯Ft2¯,whereFt1¯ represents the average number of the MS tremors in the previous t1 time period and Ft2¯ represents the average number of the MS tremors in the previous t2 time period, in which t1 < t2. According to the actual situation of WCM, let t1 = 8 h and t2 = 48 h.(2) Spatial Characteristic Parameter. The studies show that the smaller the dispersion of the source position of the MS tremor in a period of time, the greater the possibility of a large energy tremor after that; the greater the energy of the MS tremor, the greater the possibility of a large energy tremor occurring after the same degree of dispersion. Therefore, degree of dispersion (DS) is defined by combining the two dimensions of the dispersion and energy, which is a quantitative description of the risk of rock burst. The calculation method is as follows.First, calculate the weighted average distance (ds) of all MS tremors in the time window by(7)ds=∑i=1nwixi−x¯2+yi−y¯2+zi−z¯2n,where wi represents the energy of the i-th MS tremor, xi/yi/zi represent the coordinates of the i-th MS tremor, x¯/y¯/z¯ represent the average coordinates of the MS tremors within the time window, and n represents the number of MS tremors within the time window.Second, calculate the characteristic energy (Er) by(8)Er=Emax+Emedian,where Emax/Emedian represents the maximum/median energy of MS tremors within the time window.Third, sortds of all MS tremors in the time window from small to large, and find the MS tremors which located in 1.2, 7.0, 27.5, 80.8, 88.5, 95.0, and 98.0 quantiles. Average the energy of the MS tremors within 1 m around the above seven MS tremors to get the critical characteristic energy (Ec), and then use interpolation to connect the above seven points to obtain the curve of Ec which can be seen in Figure 6. Each Er point can find its corresponding Ec on the curve, and then DS can be calculated as follows:(9)Ds=lgErEc+1.Figure 6
The critical characteristic energy curve.(3) Intensity Characteristic Parameter(a)
Energy deviation (DE)The research shows that the intensity of mine tremors’ energy caused by mining activities can reflect the variation trend of rock stress and rock burst risk. Prior to the occurrence of rock burst, microcracks propagated and connected, eventually forming macrofractures, leading to an abnormal increase in MS/AE energy, andDE was higher than before. Therefore, the energy deviation (DE) of the MS and AE monitoring systems is defined as precursor to warn rock burst hazard. It can be calculated as follows:(10)DE=PE−PE¯PE¯,wherePE represents the total energy of MS/AE monitoring system within the first 24 hours and PE¯ represents the mean value of daily MS/AE monitoring system’s energy during the period from the time of the last large energy mine tremor to the current time.(b)
Total high value of energy deviation (DH)The higherDE indicates that the stress variation in the rock and the risk of rock burst are more severe than in the past. The total high value of energy deviation (DH) is extracted as the precursory characteristic parameter to improve the accuracy of rock burst risk early warning. The calculation method is shown as follows:(11)Ni=1,DEi≥C,0,DEi<C,whereDEi represents the i-th DE and Ni represents the judgment result corresponding to the i-th DE with the critical value C. One has(12)DHi=∑i=1nNi,wheren represents the length of time window and DHi represents the total number of high values exceeding the critical value C corresponding to the i-th DE. According to the actual situation of WCM, let C = 20.In summary, the multidimensional precursory characteristic parameter system of rock burst in WCM shown in Figure7 can be obtained.Figure 7
The multidimensional precursory characteristic parameter system of rock burst in WCM.The time window and sliding step length of the precursory characteristic parameters should be set according to the actual situation of WCM. If they are too long, the detailed characteristics of the data will be lost; if they are too short, the credibility of the data will be reduced (a large number of invalid data will be produced). Therefore, using the 24-hour time window and 1-hour sliding step length and taking the monitoring data within 1.5 months before and after the “02.01” rock burst as an example, plot the time series curve of precursory characteristic parameters as shown in Figure8. For convenience of explanation, the MS tremors with energy larger than 1 × 105 J are numbered according to the sequence of ①–⑤.Figure 8
The time series curve of each precursory characteristic parameters, (a) energy of MS event; (b) time series curve ofDE of MS system; (c) time series curve of DF of MS system; (d) time series curve of DS of MS system; (e) time series curve of Fr of MS system; (f) time series curve of DE of AE system; (g) time series curve of DH of AE system ① 2016/12/23 5 : 50 : 37, 1.10 E + 07 J; ② 2017/1/5 14 : 45 : 25, 1.50 E + 06 J; ③ 2017/1/8 14 : 23 : 01, 1.60 E + 05 J; ④ 2017/1/10 20 : 07 : 07, 7.10 E + 05 J; ⑤ 2017/2/1 21 : 51 : 21, 2.10 E + 08 J.
(a)(b)(c)(d)(e)(f)(g)In the case of a single precursory characteristic parameter, it is considered as accurate warning that there are significant fluctuations within two days before the occurrence of a rock burst or large energy mine tremor (hereinafter referred to as “event”) and exceeds the corresponding critical value; otherwise, it is false or missing warning. It can be seen from Figure8 that the time series variations of these six characteristic parameters have a good correlation with events and can reflect the precursory early warning information of events from different dimensions of time, space, and intensity.At the same time, it is found that the warning results of the events between different characteristic parameters are not the same; Figures8(c)–8(e) are accurate for the ④ event and Figures 8(f) and 8(g) are not. After analyzing the reasons, Figures 8(c)–8(e) are characteristic parameters of the MS monitoring system, and Figures 8(f) and 8(g) are characteristic parameters of the AE monitoring system. Because different monitoring systems have different principles and ranges, that is, the MS system monitors the fracture of high-field coal and rock mass, its vibration frequency is less than 150 Hz, the monitoring range is 10 km, the AE system monitors the crack propagation of low-field coal and rock mass, its vibration frequency is between 150 and 3000 Hz, and the monitoring range is 200 m. The two systems can separately monitor the whole process of rock burst incubation and development, and the monitoring results complement each other, so there are certain differences. Apart from this, due to the different network layout of each monitoring system, most of the AE sensors are arranged not far from the front of the working face, while the MS sensors are arranged in a wide range, which can reach the whole mining area. Because of the complex underground environment, the signals collected by different sensors may be interfered, so the results may vary from system to system. By using this feature, the characteristic parameters extracted from different monitoring systems can be integrated synthetically, which can make them complement each other’s advantages and more efficiently and accurately early warning the risk of rock burst.


#### 3.2.2. Critical Value Optimization and Early Warning Index Classification

The characteristic parameters of different monitoring systems are combined into pairs and input into the GA model, and the training is carried out in the computer simulation environment. The paper uses MATLAB R2014a to compile the code and the operation process is shown in Figure9. The code is run under the following environment: processor: AMD Ryzen 5 4600U with Radeon Graphics 2.10 GHz, RAM: 16 GB, 64 bits, and Windows 10.Figure 9
The operation process of GA in MATLAB R2014a. (a) The change of the best fitness value and the mean penalty value during the operation. (b) The vector value of the best fitness individual in each generation.
(a)(b)According to the calculation results shown in Table4, one or two pairs of characteristic parameters with the highest fitness value are selected as the precursory characteristic parameters of each monitoring system, and Wi of different monitoring systems can be obtained by substituting into equation (2), and then use equation (3) to get WCi of each monitoring system. Finally, according to the above results, IC can be obtained, and the different values correspond to four levels of no, weak, medium, and strong rock burst risk. The flowchart of multisystem and multiparameter integrated early warning method of rock burst risk in WCM is shown in Figure 10.Table 4
Optimal critical value and fitness of precursory characteristic parameters.
Name of the monitoring systemOptimal critical valueV1Optimal critical valueV2Fitness valueRARAMIS M/E MS monitoring systemDsFr0.609Optimal combination0.3641.594Suboptimal combinationDEDF0.40639.08740.815KJ623 AE monitoring systemDEDH1# sensor107.06112.6390.6712# sensor26.5857.6440.4233# sensor115.99714.8480.7154# sensor105.23712.3670.515Figure 10
The flowchart of multisystem and multiparameter integrated early warning method of rock burst risk in WCM.


### 3.3. Early Warning Results and Analysis

Due to the complex and changeable underground mine conditions, the sensors of each monitoring system are interfered by multiple factors, which may lead to problems such as the discontinuity of the monitoring data or the instability of the monitoring time. When integrating multisystem monitoring data, it is necessary to ensure that the monitoring data of each sensor is in the same time period. Therefore, in order to use as many sensors as possible to provide sufficient precursory information of rock burst risk, the monitoring data from 2017/03/07 to 2017/12/28 is selected for integration.There are ten events from 2017/03/07 to 2017/12/28, including “04.26” rock burst event. In order to avoid systematic errors, it is considered that the MS events with energy larger than 1 × 105 J, which occurred at a very close time and geographical coordinate, belong to the same large energy event. For example, the four consecutive MS events which occurred on September 21 are the same event; therefore, it is finally considered that six events have occurred in this time period, numbered according to the sequence of ①–⑥. The R-value scoring method (equation (1)) is used to calculate the early warning capability of each monitoring system in this time period. The basis for judging the success or failure of early warning is as follows: if IC exceeds 0.5 within three days before the event, the early warning is successful; otherwise, it has failed. The results are shown in Table 5.Table 5
Early warning capability of different monitoring systems.
Name of the monitoring systemn11n01N1N0adRARAMIS M/E MS monitoring system661686711910.2370.763KJ623 AE monitoring system4632371190.6670.0450.621Using equation (4) to calculate the integrated early warning index of rock burst risk IC, the result is shown in Figure 11.Figure 11
Multisystem and multiparameter integrated early warning index of rock burst risk.As shown in the red boxes in Figure11, IC has exceeded 0.5 within three days of the events. According to Table 2, it is judged as level III or IV rock burst risk, indicating that the variation of the index has an obvious correlation with the occurrence of the event, which can provide effective early warning for the event. The index synthesizes the multidimensional characteristic parameters of multisystem and provides precursory information for rock burst danger. According to the R-value scoring method, the early warning capability is R = 6/6–735/7119 = 0.896, which is significantly higher than the early warning capability of MS single system of 0.763 and the early warning capability of AE single system of 0.621. This is because when the multisystem early warning results are integrated, different monitoring systems provide each other with early warning criterion and multidimensional characteristic parameters’ early warning results verify each other, eliminating the influence of single system accidental factors, thus improving the early warning accuracy and finally making its early warning capability higher than that of a single monitoring system.The integrated early warning indexIC can effectively evaluate the risk degree of internal dynamic disaster of coal and rock mass, and it is a single index integrating multiple systems and parameters, which is easy to read and can provide early warning information for mines in real time and quantitatively. More efficient and accurate prevention and control of underground dynamic disasters, thus avoiding major casualties and property losses, are of great significance to the production field.


## 3.1. Basic Conditions of the No.B3 + 6 Working Face, WCM

WCM is located in Urumqi, Xinjiang, China. The overall geological section of the north and south mining area is shown in Figure3, and up to now, the current mining situation of WCM is at No. B3 + 6 + 425 m working face. Both the No. B3 + 6 and No. B1 + 2 are classified as steeply inclined and extremely thick coal seam, and there is a rock pillar with 70–101 m thickness between them. The average inclination of the coal seam is 87°. The horizontal section full-mechanized top-coal caving mining technology is being adopted for mining, the ratio of mining to caving height is 1 : 7, and the horizontal section height is 25 m. Due to the special occurrence conditions of the coal seam, with the increase of the mining depth, a towering rock pillar is formed between the two coal seams. The rock pillar is affected by the horizontal tectonic stress and self-heavy stress, which causes the tilting effect of the intersection between the lower side of the rock column and the coal seam. Under the disturbance of mining, the elastic energy of rock mass accumulation is released through the coal body, which is easy to produce different degrees of dynamic disasters. There are six times of rock burst in the process of No. B3 + 6 coal seam mining.Figure 3
The overall geological section of WCM.At present, the mine has installed various monitoring systems such as ARAMIS M/E MS monitoring system and KJ623 AE monitoring system. The ARAMIS M/E is a particular MS monitoring system, which can automatically collect and filter seismic signals, and then calculates the occurrence time, released energy, and three-dimensional coordinates of both tremors and rock bursts. The south mining area of the WCM was surrounded by a total of 12 sensors, each sensor covered a range of 2 km, and the layout of the sensors until 26 April 2017 when the last rock burst occurred is shown in Figure4 (please refer to Li et al. 2014 for more details). The KJ623 AE monitoring system contains a total of 11 sensors that was installed on the coal and rock mass to collect AE events within a frequency range of 300–2000 Hz. The sensitivity of the AE transducer was 65 ± 3 dB (V/m/s), the bolt was used as the AE conductor, and the sensor was installed on the bolt. The sensors are located at 30 m in front of the working face, and the spacing between each sensor is 30–40 m. As the working face moves, the sensor layout moves accordingly.Figure 4
Network layout of MS and AE monitoring system in No.B3 + 6 coal seam +450 horizontal working face of WCM during the research period.According to the actual situation of WCM, the large energy mine tremors are defined as MS event with energy exceeding 1 × 105J. Taking the No. B3 + 6 coal seam + 450 horizontal working face of the coal mine as the research object, the MS and AE monitoring data of 2016/06/06–2017/12/28 time period are selected (see Tables S1-S2 in the Supplementary Material for complete field monitoring data of MS and AE). During this period, there are three rock burst events and ten large energy mine tremors in total, and the distribution of rock burst events is shown in Figure 5. Table 3 shows the specific manifestation records of these three rock burst events.Figure 5
The distribution map of MS events with energy exceeding 1 × 103J in No. B3 + 6 coal seam + 450 horizontal working face WCM during 2016/06/06–2017/12/28.Table 3
On-site rock burst manifestation records.
DateLocation of working face (m)Location of MS sourceMS energy (J)Damage form and degree2016/11/242026Located in theB2-B3 rock wall, 12 m away from the B3 bottom plate9.5 × 106Upper wall moved 0.3 m, lower wall moved 0.5 m, floor heaved by 0.4 m, damage to hydraulic props2017/02/011824Located at 9 m from theB6 roof slate2.1 × 108Upper wall moved 0.3–0.7 m, lower wall moved 0.5–1.3 m, damage to belt conveyor2017/04/261568Located in theB2-B3 rock pillar, 35 m away from the B3 bottom plate2.2 × 106Lower wall moved 0.3–0.5 m, floor heaved by 0.2–0.3 m, damage to hydraulic props and belt conveyor, reversed loader change direction


## 3.2. Application of the Integrated Early Warning Method in WCM

### 3.2.1. Construction of Precursory Characteristic Parameter System

(1) Temporal Characteristic Parameter(a)
Frequency deviation (DF)The information of the mine tremor monitored by the MS system can reflect the internal stress of the rock mass. Prior to the occurrence of rock burst, the source of the MS increased and clustered anomaly compared to the previous normal mining period, which shows that the stress in the study area is close to the strength limit of coal and rock mass. Therefore, the frequency deviation value (DF) of a certain time window and sliding step length in the monitoring area is calculated. It can reflect the trend of the MS events over a period of time; the denser the MS tremor is in time sequence, the larger the value is. It can be calculated as follows:(5)DF=PF−PF¯PF¯,wherePF represents the number of MS tremors within the first 24 hours and PF¯ represents the mean value of daily tremors during the period from the time of the last large energy mine tremor to the current time.(b)
Frequency ratio (Fr)The higher the frequency of MS tremor is, the more frequent the stress activities in the rock mass are and the more likely the rock burst is to be induced; otherwise, the less likely it is. Therefore, the frequency ratio (Fr) of the characteristic parameter is defined, and the calculation method is as follows:(6)Fr=Ft1¯Ft2¯,whereFt1¯ represents the average number of the MS tremors in the previous t1 time period and Ft2¯ represents the average number of the MS tremors in the previous t2 time period, in which t1 < t2. According to the actual situation of WCM, let t1 = 8 h and t2 = 48 h.(2) Spatial Characteristic Parameter. The studies show that the smaller the dispersion of the source position of the MS tremor in a period of time, the greater the possibility of a large energy tremor after that; the greater the energy of the MS tremor, the greater the possibility of a large energy tremor occurring after the same degree of dispersion. Therefore, degree of dispersion (DS) is defined by combining the two dimensions of the dispersion and energy, which is a quantitative description of the risk of rock burst. The calculation method is as follows.First, calculate the weighted average distance (ds) of all MS tremors in the time window by(7)ds=∑i=1nwixi−x¯2+yi−y¯2+zi−z¯2n,where wi represents the energy of the i-th MS tremor, xi/yi/zi represent the coordinates of the i-th MS tremor, x¯/y¯/z¯ represent the average coordinates of the MS tremors within the time window, and n represents the number of MS tremors within the time window.Second, calculate the characteristic energy (Er) by(8)Er=Emax+Emedian,where Emax/Emedian represents the maximum/median energy of MS tremors within the time window.Third, sortds of all MS tremors in the time window from small to large, and find the MS tremors which located in 1.2, 7.0, 27.5, 80.8, 88.5, 95.0, and 98.0 quantiles. Average the energy of the MS tremors within 1 m around the above seven MS tremors to get the critical characteristic energy (Ec), and then use interpolation to connect the above seven points to obtain the curve of Ec which can be seen in Figure 6. Each Er point can find its corresponding Ec on the curve, and then DS can be calculated as follows:(9)Ds=lgErEc+1.Figure 6
The critical characteristic energy curve.(3) Intensity Characteristic Parameter(a)
Energy deviation (DE)The research shows that the intensity of mine tremors’ energy caused by mining activities can reflect the variation trend of rock stress and rock burst risk. Prior to the occurrence of rock burst, microcracks propagated and connected, eventually forming macrofractures, leading to an abnormal increase in MS/AE energy, andDE was higher than before. Therefore, the energy deviation (DE) of the MS and AE monitoring systems is defined as precursor to warn rock burst hazard. It can be calculated as follows:(10)DE=PE−PE¯PE¯,wherePE represents the total energy of MS/AE monitoring system within the first 24 hours and PE¯ represents the mean value of daily MS/AE monitoring system’s energy during the period from the time of the last large energy mine tremor to the current time.(b)
Total high value of energy deviation (DH)The higherDE indicates that the stress variation in the rock and the risk of rock burst are more severe than in the past. The total high value of energy deviation (DH) is extracted as the precursory characteristic parameter to improve the accuracy of rock burst risk early warning. The calculation method is shown as follows:(11)Ni=1,DEi≥C,0,DEi<C,whereDEi represents the i-th DE and Ni represents the judgment result corresponding to the i-th DE with the critical value C. One has(12)DHi=∑i=1nNi,wheren represents the length of time window and DHi represents the total number of high values exceeding the critical value C corresponding to the i-th DE. According to the actual situation of WCM, let C = 20.In summary, the multidimensional precursory characteristic parameter system of rock burst in WCM shown in Figure7 can be obtained.Figure 7
The multidimensional precursory characteristic parameter system of rock burst in WCM.The time window and sliding step length of the precursory characteristic parameters should be set according to the actual situation of WCM. If they are too long, the detailed characteristics of the data will be lost; if they are too short, the credibility of the data will be reduced (a large number of invalid data will be produced). Therefore, using the 24-hour time window and 1-hour sliding step length and taking the monitoring data within 1.5 months before and after the “02.01” rock burst as an example, plot the time series curve of precursory characteristic parameters as shown in Figure8. For convenience of explanation, the MS tremors with energy larger than 1 × 105 J are numbered according to the sequence of ①–⑤.Figure 8
The time series curve of each precursory characteristic parameters, (a) energy of MS event; (b) time series curve ofDE of MS system; (c) time series curve of DF of MS system; (d) time series curve of DS of MS system; (e) time series curve of Fr of MS system; (f) time series curve of DE of AE system; (g) time series curve of DH of AE system ① 2016/12/23 5 : 50 : 37, 1.10 E + 07 J; ② 2017/1/5 14 : 45 : 25, 1.50 E + 06 J; ③ 2017/1/8 14 : 23 : 01, 1.60 E + 05 J; ④ 2017/1/10 20 : 07 : 07, 7.10 E + 05 J; ⑤ 2017/2/1 21 : 51 : 21, 2.10 E + 08 J.
(a)(b)(c)(d)(e)(f)(g)In the case of a single precursory characteristic parameter, it is considered as accurate warning that there are significant fluctuations within two days before the occurrence of a rock burst or large energy mine tremor (hereinafter referred to as “event”) and exceeds the corresponding critical value; otherwise, it is false or missing warning. It can be seen from Figure8 that the time series variations of these six characteristic parameters have a good correlation with events and can reflect the precursory early warning information of events from different dimensions of time, space, and intensity.At the same time, it is found that the warning results of the events between different characteristic parameters are not the same; Figures8(c)–8(e) are accurate for the ④ event and Figures 8(f) and 8(g) are not. After analyzing the reasons, Figures 8(c)–8(e) are characteristic parameters of the MS monitoring system, and Figures 8(f) and 8(g) are characteristic parameters of the AE monitoring system. Because different monitoring systems have different principles and ranges, that is, the MS system monitors the fracture of high-field coal and rock mass, its vibration frequency is less than 150 Hz, the monitoring range is 10 km, the AE system monitors the crack propagation of low-field coal and rock mass, its vibration frequency is between 150 and 3000 Hz, and the monitoring range is 200 m. The two systems can separately monitor the whole process of rock burst incubation and development, and the monitoring results complement each other, so there are certain differences. Apart from this, due to the different network layout of each monitoring system, most of the AE sensors are arranged not far from the front of the working face, while the MS sensors are arranged in a wide range, which can reach the whole mining area. Because of the complex underground environment, the signals collected by different sensors may be interfered, so the results may vary from system to system. By using this feature, the characteristic parameters extracted from different monitoring systems can be integrated synthetically, which can make them complement each other’s advantages and more efficiently and accurately early warning the risk of rock burst.


### 3.2.2. Critical Value Optimization and Early Warning Index Classification

The characteristic parameters of different monitoring systems are combined into pairs and input into the GA model, and the training is carried out in the computer simulation environment. The paper uses MATLAB R2014a to compile the code and the operation process is shown in Figure9. The code is run under the following environment: processor: AMD Ryzen 5 4600U with Radeon Graphics 2.10 GHz, RAM: 16 GB, 64 bits, and Windows 10.Figure 9
The operation process of GA in MATLAB R2014a. (a) The change of the best fitness value and the mean penalty value during the operation. (b) The vector value of the best fitness individual in each generation.
(a)(b)According to the calculation results shown in Table4, one or two pairs of characteristic parameters with the highest fitness value are selected as the precursory characteristic parameters of each monitoring system, and Wi of different monitoring systems can be obtained by substituting into equation (2), and then use equation (3) to get WCi of each monitoring system. Finally, according to the above results, IC can be obtained, and the different values correspond to four levels of no, weak, medium, and strong rock burst risk. The flowchart of multisystem and multiparameter integrated early warning method of rock burst risk in WCM is shown in Figure 10.Table 4
Optimal critical value and fitness of precursory characteristic parameters.
Name of the monitoring systemOptimal critical valueV1Optimal critical valueV2Fitness valueRARAMIS M/E MS monitoring systemDsFr0.609Optimal combination0.3641.594Suboptimal combinationDEDF0.40639.08740.815KJ623 AE monitoring systemDEDH1# sensor107.06112.6390.6712# sensor26.5857.6440.4233# sensor115.99714.8480.7154# sensor105.23712.3670.515Figure 10
The flowchart of multisystem and multiparameter integrated early warning method of rock burst risk in WCM.


## 3.2.1. Construction of Precursory Characteristic Parameter System

(1) Temporal Characteristic Parameter(a)
Frequency deviation (DF)The information of the mine tremor monitored by the MS system can reflect the internal stress of the rock mass. Prior to the occurrence of rock burst, the source of the MS increased and clustered anomaly compared to the previous normal mining period, which shows that the stress in the study area is close to the strength limit of coal and rock mass. Therefore, the frequency deviation value (DF) of a certain time window and sliding step length in the monitoring area is calculated. It can reflect the trend of the MS events over a period of time; the denser the MS tremor is in time sequence, the larger the value is. It can be calculated as follows:(5)DF=PF−PF¯PF¯,wherePF represents the number of MS tremors within the first 24 hours and PF¯ represents the mean value of daily tremors during the period from the time of the last large energy mine tremor to the current time.(b)
Frequency ratio (Fr)The higher the frequency of MS tremor is, the more frequent the stress activities in the rock mass are and the more likely the rock burst is to be induced; otherwise, the less likely it is. Therefore, the frequency ratio (Fr) of the characteristic parameter is defined, and the calculation method is as follows:(6)Fr=Ft1¯Ft2¯,whereFt1¯ represents the average number of the MS tremors in the previous t1 time period and Ft2¯ represents the average number of the MS tremors in the previous t2 time period, in which t1 < t2. According to the actual situation of WCM, let t1 = 8 h and t2 = 48 h.(2) Spatial Characteristic Parameter. The studies show that the smaller the dispersion of the source position of the MS tremor in a period of time, the greater the possibility of a large energy tremor after that; the greater the energy of the MS tremor, the greater the possibility of a large energy tremor occurring after the same degree of dispersion. Therefore, degree of dispersion (DS) is defined by combining the two dimensions of the dispersion and energy, which is a quantitative description of the risk of rock burst. The calculation method is as follows.First, calculate the weighted average distance (ds) of all MS tremors in the time window by(7)ds=∑i=1nwixi−x¯2+yi−y¯2+zi−z¯2n,where wi represents the energy of the i-th MS tremor, xi/yi/zi represent the coordinates of the i-th MS tremor, x¯/y¯/z¯ represent the average coordinates of the MS tremors within the time window, and n represents the number of MS tremors within the time window.Second, calculate the characteristic energy (Er) by(8)Er=Emax+Emedian,where Emax/Emedian represents the maximum/median energy of MS tremors within the time window.Third, sortds of all MS tremors in the time window from small to large, and find the MS tremors which located in 1.2, 7.0, 27.5, 80.8, 88.5, 95.0, and 98.0 quantiles. Average the energy of the MS tremors within 1 m around the above seven MS tremors to get the critical characteristic energy (Ec), and then use interpolation to connect the above seven points to obtain the curve of Ec which can be seen in Figure 6. Each Er point can find its corresponding Ec on the curve, and then DS can be calculated as follows:(9)Ds=lgErEc+1.Figure 6
The critical characteristic energy curve.(3) Intensity Characteristic Parameter(a)
Energy deviation (DE)The research shows that the intensity of mine tremors’ energy caused by mining activities can reflect the variation trend of rock stress and rock burst risk. Prior to the occurrence of rock burst, microcracks propagated and connected, eventually forming macrofractures, leading to an abnormal increase in MS/AE energy, andDE was higher than before. Therefore, the energy deviation (DE) of the MS and AE monitoring systems is defined as precursor to warn rock burst hazard. It can be calculated as follows:(10)DE=PE−PE¯PE¯,wherePE represents the total energy of MS/AE monitoring system within the first 24 hours and PE¯ represents the mean value of daily MS/AE monitoring system’s energy during the period from the time of the last large energy mine tremor to the current time.(b)
Total high value of energy deviation (DH)The higherDE indicates that the stress variation in the rock and the risk of rock burst are more severe than in the past. The total high value of energy deviation (DH) is extracted as the precursory characteristic parameter to improve the accuracy of rock burst risk early warning. The calculation method is shown as follows:(11)Ni=1,DEi≥C,0,DEi<C,whereDEi represents the i-th DE and Ni represents the judgment result corresponding to the i-th DE with the critical value C. One has(12)DHi=∑i=1nNi,wheren represents the length of time window and DHi represents the total number of high values exceeding the critical value C corresponding to the i-th DE. According to the actual situation of WCM, let C = 20.In summary, the multidimensional precursory characteristic parameter system of rock burst in WCM shown in Figure7 can be obtained.Figure 7
The multidimensional precursory characteristic parameter system of rock burst in WCM.The time window and sliding step length of the precursory characteristic parameters should be set according to the actual situation of WCM. If they are too long, the detailed characteristics of the data will be lost; if they are too short, the credibility of the data will be reduced (a large number of invalid data will be produced). Therefore, using the 24-hour time window and 1-hour sliding step length and taking the monitoring data within 1.5 months before and after the “02.01” rock burst as an example, plot the time series curve of precursory characteristic parameters as shown in Figure8. For convenience of explanation, the MS tremors with energy larger than 1 × 105 J are numbered according to the sequence of ①–⑤.Figure 8
The time series curve of each precursory characteristic parameters, (a) energy of MS event; (b) time series curve ofDE of MS system; (c) time series curve of DF of MS system; (d) time series curve of DS of MS system; (e) time series curve of Fr of MS system; (f) time series curve of DE of AE system; (g) time series curve of DH of AE system ① 2016/12/23 5 : 50 : 37, 1.10 E + 07 J; ② 2017/1/5 14 : 45 : 25, 1.50 E + 06 J; ③ 2017/1/8 14 : 23 : 01, 1.60 E + 05 J; ④ 2017/1/10 20 : 07 : 07, 7.10 E + 05 J; ⑤ 2017/2/1 21 : 51 : 21, 2.10 E + 08 J.
(a)(b)(c)(d)(e)(f)(g)In the case of a single precursory characteristic parameter, it is considered as accurate warning that there are significant fluctuations within two days before the occurrence of a rock burst or large energy mine tremor (hereinafter referred to as “event”) and exceeds the corresponding critical value; otherwise, it is false or missing warning. It can be seen from Figure8 that the time series variations of these six characteristic parameters have a good correlation with events and can reflect the precursory early warning information of events from different dimensions of time, space, and intensity.At the same time, it is found that the warning results of the events between different characteristic parameters are not the same; Figures8(c)–8(e) are accurate for the ④ event and Figures 8(f) and 8(g) are not. After analyzing the reasons, Figures 8(c)–8(e) are characteristic parameters of the MS monitoring system, and Figures 8(f) and 8(g) are characteristic parameters of the AE monitoring system. Because different monitoring systems have different principles and ranges, that is, the MS system monitors the fracture of high-field coal and rock mass, its vibration frequency is less than 150 Hz, the monitoring range is 10 km, the AE system monitors the crack propagation of low-field coal and rock mass, its vibration frequency is between 150 and 3000 Hz, and the monitoring range is 200 m. The two systems can separately monitor the whole process of rock burst incubation and development, and the monitoring results complement each other, so there are certain differences. Apart from this, due to the different network layout of each monitoring system, most of the AE sensors are arranged not far from the front of the working face, while the MS sensors are arranged in a wide range, which can reach the whole mining area. Because of the complex underground environment, the signals collected by different sensors may be interfered, so the results may vary from system to system. By using this feature, the characteristic parameters extracted from different monitoring systems can be integrated synthetically, which can make them complement each other’s advantages and more efficiently and accurately early warning the risk of rock burst.


## 3.2.2. Critical Value Optimization and Early Warning Index Classification

The characteristic parameters of different monitoring systems are combined into pairs and input into the GA model, and the training is carried out in the computer simulation environment. The paper uses MATLAB R2014a to compile the code and the operation process is shown in Figure9. The code is run under the following environment: processor: AMD Ryzen 5 4600U with Radeon Graphics 2.10 GHz, RAM: 16 GB, 64 bits, and Windows 10.Figure 9
The operation process of GA in MATLAB R2014a. (a) The change of the best fitness value and the mean penalty value during the operation. (b) The vector value of the best fitness individual in each generation.
(a)(b)According to the calculation results shown in Table4, one or two pairs of characteristic parameters with the highest fitness value are selected as the precursory characteristic parameters of each monitoring system, and Wi of different monitoring systems can be obtained by substituting into equation (2), and then use equation (3) to get WCi of each monitoring system. Finally, according to the above results, IC can be obtained, and the different values correspond to four levels of no, weak, medium, and strong rock burst risk. The flowchart of multisystem and multiparameter integrated early warning method of rock burst risk in WCM is shown in Figure 10.Table 4
Optimal critical value and fitness of precursory characteristic parameters.
Name of the monitoring systemOptimal critical valueV1Optimal critical valueV2Fitness valueRARAMIS M/E MS monitoring systemDsFr0.609Optimal combination0.3641.594Suboptimal combinationDEDF0.40639.08740.815KJ623 AE monitoring systemDEDH1# sensor107.06112.6390.6712# sensor26.5857.6440.4233# sensor115.99714.8480.7154# sensor105.23712.3670.515Figure 10
The flowchart of multisystem and multiparameter integrated early warning method of rock burst risk in WCM.


## 3.3. Early Warning Results and Analysis

Due to the complex and changeable underground mine conditions, the sensors of each monitoring system are interfered by multiple factors, which may lead to problems such as the discontinuity of the monitoring data or the instability of the monitoring time. When integrating multisystem monitoring data, it is necessary to ensure that the monitoring data of each sensor is in the same time period. Therefore, in order to use as many sensors as possible to provide sufficient precursory information of rock burst risk, the monitoring data from 2017/03/07 to 2017/12/28 is selected for integration.There are ten events from 2017/03/07 to 2017/12/28, including “04.26” rock burst event. In order to avoid systematic errors, it is considered that the MS events with energy larger than 1 × 105 J, which occurred at a very close time and geographical coordinate, belong to the same large energy event. For example, the four consecutive MS events which occurred on September 21 are the same event; therefore, it is finally considered that six events have occurred in this time period, numbered according to the sequence of ①–⑥. The R-value scoring method (equation (1)) is used to calculate the early warning capability of each monitoring system in this time period. The basis for judging the success or failure of early warning is as follows: if IC exceeds 0.5 within three days before the event, the early warning is successful; otherwise, it has failed. The results are shown in Table 5.Table 5
Early warning capability of different monitoring systems.
Name of the monitoring systemn11n01N1N0adRARAMIS M/E MS monitoring system661686711910.2370.763KJ623 AE monitoring system4632371190.6670.0450.621Using equation (4) to calculate the integrated early warning index of rock burst risk IC, the result is shown in Figure 11.Figure 11
Multisystem and multiparameter integrated early warning index of rock burst risk.As shown in the red boxes in Figure11, IC has exceeded 0.5 within three days of the events. According to Table 2, it is judged as level III or IV rock burst risk, indicating that the variation of the index has an obvious correlation with the occurrence of the event, which can provide effective early warning for the event. The index synthesizes the multidimensional characteristic parameters of multisystem and provides precursory information for rock burst danger. According to the R-value scoring method, the early warning capability is R = 6/6–735/7119 = 0.896, which is significantly higher than the early warning capability of MS single system of 0.763 and the early warning capability of AE single system of 0.621. This is because when the multisystem early warning results are integrated, different monitoring systems provide each other with early warning criterion and multidimensional characteristic parameters’ early warning results verify each other, eliminating the influence of single system accidental factors, thus improving the early warning accuracy and finally making its early warning capability higher than that of a single monitoring system.The integrated early warning indexIC can effectively evaluate the risk degree of internal dynamic disaster of coal and rock mass, and it is a single index integrating multiple systems and parameters, which is easy to read and can provide early warning information for mines in real time and quantitatively. More efficient and accurate prevention and control of underground dynamic disasters, thus avoiding major casualties and property losses, are of great significance to the production field.


## 4. Early Warning Verification

In practical applications, for newly installed monitoring systems or mines, mining new coal seams, the data to be monitored shall be at least two months, so as to ensure that there is sufficient data for training to avoid overfitting; if different monitoring systems have been installed and the mining coal seams have not been replaced in a certain period of time, the method described in Section2 can be directly applied to calculate the integrated early warning index IC.The specific calculation method is as follows: take 24 hours as the time window, 1 hour as the sliding step to calculate the characteristic parameters, and one month as the drawing time window. Taking the monitoring data of the No.B3 + 6 coal seam + 450 horizontal working face in January 2018 as an example, there were four MS tremors with energy over 1 × 105 J. Similarly, the two MS tremors that occurred on 01/06 belong to the same event. For the convenience of explanation, all events are numbered according to the sequence of ①–③. The time series curve of each precursor characteristic parameter for this time period is shown in Figure 12.Figure 12
Time series curve of precursory characteristic parameters of each system in January 2018, (a) energy of MS event; (b) time series curve ofDE of MS system; (c) time series curve of DF of MS system; (d) time series curve of DS of MS system; (e) time series curve of Fr of MS system; (f) time series curve of DE of AE system; (g) time series curve of DH of AE system ① 2018/1/6 04 : 30 : 48, 1.9 E + 05 J; 2018/1/6 04 : 30 : 31, 9.8 E + 05 J; ② 2018/1/16 01 : 25 : 45, 4.50 E + 05 J; ③ 2018/1/21 04 : 07 : 24, 6.6 E + 06 J.
(a)(b)(c)(d)(e)(f)(g)It can be seen from Figure12 that the variation of each characteristic parameter has a significant correlation with the energy of MS tremor, which basically exceeds the optimal critical value for a period of time before the occurrence of the event and does not fluctuate much most of the time. Except that DS and Fr are more sensitive than other parameters, it is low affected by noise and has high early warning capability.Combined with the optimal critical value and fitness value of the same coal seam historical monitoring data, the real-timeIC is calculated, and the calculation results are shown in Figure 13.Figure 13
Integrated early warning indexIC of rock burst risk in January 2018 ①2018/1/6 04 : 30 : 48, 1.9 E + 05J; 2018/1/6 04 : 30 : 31, 9.8 E + 05J; ② 2018/1/16 01 : 25 : 45, 4.50 E + 05J; ③ 2018/1/21 04 : 07 : 24, 6.6 E + 06J.It can be seen from Figure13 that the index in the red box has reached more than 0.5 within three days before the occurrence of the event, which is in the warning level of III or IV, and the number of false warnings is less, which greatly reduces the impact of early warning on normal production activities of the mine and achieves the balance between safe production and efficient production.The verification results show that the method can reflect the rock burst risk status, specific risk areas, and levels of the current monitoring area in real time and quantitatively, eliminate the influence of subjective factors in the early warning process, and improve the accuracy of the early warning of the rock burst risk. Combined with the different risk levels in Table2, it can provide useful guidance for mining enterprises to prevent coal and rock dynamic disasters.


## 5. Conclusion

Rock burst often causes huge casualties and property losses, and accurate monitoring and early warning is the key to disaster prevention. Therefore, this paper constructs a multidimensional precursory characteristic parameter system, and proposes a multisystem and multiparameter integrated rock burst risk early warning method based on GA, which forms a highly reliable risk evaluation method with unified criterion, index, and critical value. The conclusions are as follows:(1)
Combining the GA with theR-value scoring method to optimize the critical value of the precursory characteristic parameters of each monitoring system, the optimal critical value and the corresponding early warning capability are obtained. The influence of missing or false warning is considered comprehensively, and the weight coefficients of each monitoring system can be adjusted according to the real-time calculation results. Under the long-term monitoring data, it still has high accuracy and can exclude the influence of subjective factors on the results in the early warning process.(2)
The multidimensional precursory characteristic parameter system of WCM is constructed, including energy deviation (DE), frequency ratio (Fr), frequency deviation (DF), degree of dispersion (DS), and total high value of energy deviation (DH). The system combines the information of temporal-spatial-intensity rock burst incubation process with high accuracy and clear physical significance.(3)
The example shows that the proposed multisystem and multiparameter integrated rock burst risk early warning method has an early warning capability of 0.896, which is higher than 0.763 of the MS monitoring system and 0.621 of the AE monitoring system. It realizes the complementary advantages of multiple monitoring systems, improves the early warning accuracy, and further complements the early warning method of the underground dynamic disaster.(4)
The integrated early warning indexIC of rock burst risk is proposed, which can reflect the change of rock burst risk in coal and rock in real time, and its quantification is divided into four levels. Combined with the specific measures corresponding to different risk levels, it can provide effective guidance for the disaster prevention work.



---

*Source: 8865654-2020-09-03.xml*