A novel combined intelligent algorithm prediction model for the risk of the coal and gas outburst

The mechanism of coal and gas outburst disasters is perplexing, and the evaluation methods of outburst disasters based on various sensitive indicators often have some imprecision and fuzziness. With the concept of accurate and intelligent mining in coal mines proposed in China, selecting quantifiable parameters for machine learning risk prediction can avoid the deviation caused by human subjectivity, and improve the accuracy of coal and gas outburst prediction. Aiming at the shortcomings of the support vector machine (SVM) such as low noise resistance and being prone to be influenced by parameters easily, this research proposed a prediction method based on a grey wolf optimizer to optimize the support vector machine (GWO-SVM). To coordinate the global and local optimization ability of the GWO, Tent Chaotic Mapping and DLH strategies were introduced to improve the optimization ability of the GWO and reduce the local optimal probability. The improved prediction model IGWO-SVM was used to predict the coal and gas outburst. The results showed that this model has faster training speed and higher classification prediction accuracy than the SVM and GWO-SVM models, which the accuracy rate reaching 100%. Finally, to obtain the correlation between the parameters of the coal and gas outburst prediction parameters, the random forest algorithm was used for training, and the three parameters with the highest feature importance were selected to rebuild the data set for machine learning. The accuracy of the IGWO-SVM outburst prediction model based on Random Forest was still 100%. Therefore, even if some prediction parameters are missing, the outburst can still be effectively predicted by using the RF-IGWO-SVM model, which is beneficial for the model application and underground safety management.

In 2022, the total output of coal mines in China is 4.45 billion tons, and coal consumption accounts for 56.2% of primary energy consumption 1 .In the future, coal will still be the most important basic energy in China, and it will play an indispensable role in economic and social development.According to the accident statistics published by the websites of the Ministry of Emergency Management of the People's Republic of China annually, with the exploitation of coal resources, coal mine accidents can never be completely avoided, and the number of coal mine accidents in 2022 increased compared with that in 2021, which can be seen from Fig. 1a.During the decade from 2012 to 2022, 68 coal and gas outburst accidents occurred, resulting in 429 deaths as shown in Fig. 1b.In Fig. 1c and d, the number of coal and gas outburst accidents and deaths showed a downward trend as a whole, while fluctuating and rising in recent years, indicating that major accidents were effectively controlled as a whole.The coal and gas outburst is a dynamic disaster in which the gas-bearing coal in the mine suddenly moves rapidly from the coal seam to the mining space, accompanied by a large amount of gas emission, which seriously threatens the life and property safety of coal mine workers.China, Russia, the United States, Poland, Australia, and other major coal-mining countries have continuously invested abundant resources to actively solve the grievous disaster.China is one of the countries with the most severe dynamic disasters in the world.
At present, there are still about 840 high-gas mines and 719 coal and gas outburst mines in China.Besides, with the increasing mining depth, some low-gas mines have evolved into high-outburst mines, and gas burning and explosion accidents and gas overrun in the working face still occur from time to time 2 .Gas content is often used to evaluate outburst risk in coal mines, but due to the change in mining conditions, gas pressure, coefficient of coal firmness, and permeability should also be used to evaluate outburst risk in coal mines, and appropriate outburst threshold limits and controls should be determined 3 .While more and more coal mines are exploiting areas of high gas content and low permeability, the combination of these two factors often reduces the efficiency and effectiveness of gas drainage.It is difficult to reduce the gas content in the coal seam to below the outburst threshold.To control the risk, a large number of experts and scholars have explored the mechanism of outbursts and made great progress.Among them, the comprehensive action hypothesis of outburst is widely recognized, which reveals that outburst is the result of the joint action of the physical and mechanical properties of stress, gas, and coal 4 .Gas affects the effective stress of coal through the adsorption expansion effect and pore pressure.The adsorption expansion effect changes the contact between coal particles and affects the overall bearing of stress, while free gas in pores shares part of the external load 5 .The description of the mechanical behavior of coal-containing gas is mostly based on Terzaghi's effective stress theory, and the elastic-plastic model is used to explain the change, elastic deformation, and failure of the stress path.Experimental evidence suggests that the higher the gas pressure, the greater the expansion deformation, the lower the porosity, and the worse coal seam permeability 6 .It is assumed that under the comprehensive influence of stress, pore pressure, and temperature, the permeability of the coal seam in front of the working face may suddenly increase, causing coal and gas outburst disasters by studying the mechanical properties under various loading and unloading conditions 7 .Xue et al. and Wang et al. through the mechanical test of gas-bearing coal and coal-rock combination under triaxial conditions, the results show that the deformation of coal adsorbed gas has a significant impact on permeability 8,9 .Wang et al. suggested that gas desorption weakens the strength of coal, which is considered to be due to the reduced effective stress and the impact of the internal gas release on coal 10 .Sobczyk and Skoczylas studied the influence of adsorption and desorption processes on outbursts, discussed the danger of outburst of gas-bearing coal with different strength, then concluded that stress and mechanical properties of coal were the main factors determining outbursts strengths 11,12 .Ma and Xue quantitatively studied the energy intensity threshold of coal and gas outbursts under different gas pressures through laboratory tests 13,14 .Compared with gas in free pore gas energy, energy from desorption gas accounts for about 90% of outburst effective energy.Therefore, the laboratory experiments result show that coal and gas outburst is driven not only by gas but also by stress, and coal seam conditions.The multi-factor theory provides a solid theoretical basis for the risk assessment and prediction of coal and gas outburst, and the research on the risk assessment of coal and gas outburst based on the comprehensive action hypothesis combined with various indicators is plentiful [15][16][17][18][19] .The mechanism of coal and gas outburst disasters is complex, and the evaluation methods of outburst disasters based on various sensitive parameters often have some inaccuracies and fuzziness.Many studies have been conducted based on prediction algorithms and models in coal mines [20][21][22][23][24][25] .With the concept of accurate and intelligent mining in coal mines put forward, researchers pay more and more attention to machine learning to predict coal and gas outbursts [26][27][28][29] .Such as, Particle Swarm Optimization (PSO) is used to optimize the prediction model of the back propagation (BP) algorithm 26 , random forest model is used to predict coal and gas outbursts 27 , and genetic algorithm is used to optimize the support vector machine model 28 .Given the lack of sample data, Zheng et al. used the data mining Multiple Imputation (MI) method to fill the missing data, and used the support vector machine (SVM) to predict coal and gas outburst 29 .The above models have significantly improved the prediction precision of coal and gas outburst, but there are still some improvements: for example, the optimization function of the GA algorithm needs to adjust more parameters, and the convergence speed of the GA needs to be improved, while the PSO algorithm has a fast convergence speed, but easily falls into local optimization, SVM, BP and other algorithms are weak at generalizing, which reduces the predicting power to some extent.The grey wolf optimization algorithm 30 is a new type of swarm intelligence optimization algorithm.This algorithm optimizes the search by learning the leadership mode of the grey wolf population and the process of hunting prey.It has the characteristics of a simple principle, few parameter settings, and strong global search capability, but it has the drawbacks of strong development capability and weak exploration capability, and it is easy to fall into local optimization in the later stages.To solve this problem, the grey wolf optimizer algorithm can be improved by extracting features of the predictive parameters, using Tent Chaotic Mapping and DLH strategy 31 .
Given the aforementioned analysis, the mechanism and development process of coal and gas outbursts were investigated, and the outburst prediction parameters reasonably were determined.The data samples of Pansan Mine in Huainan Basin were tested, and the core parameters of the support vector machine were optimized by IGWO.The prediction accuracy of coal and gas outburst was compared between SVM and IGWO-SVM, which showed the IGWO-SVM model could accurately predict coal and gas outburst, and the effectiveness and superiority of the model were verified.Finally, the prediction parameters were extracted by Random Forest and trained by the models, the prediction results showed that the accuracy of the models was not reduced when three prediction parameters were selected.Therefore, Random Forest can effectively reduce the problems of diversified data types and large data volumes.The main contributions of this study are as follows: (1) According to the multi-factor outburst theory and energy conversion theory, the prediction parameters were reasonably determined, and the prediction parameters did not need complex or subjective data processing; (2) The grey wolf algorithm was improved and its effectiveness was verified; (3) Coal and gas outburst can be accurately predicted using the IGWO-SVM model; (4) The random forest algorithm was used to extract the features of the parameters, and the accuracy of the prediction models was not reduced by using the main parameters for training; (5) The superiority of the prediction model was verified by field and experimental data.It is helpful to accurately predict coal and gas outburst in coal mines for effectively preventing outburst disasters, formulating reasonable outburst prevention measures, improving the safety management level, and promoting the safety production capacity of coal mines.

Mechanism of the coal and gas outburst
Multi-factors theory of coal and gas outburst Multi-factor coal and gas theory is based on long-term human exploration of outburst and includes factors such as mine disturbance, effective stress, gas flow, and physical and mechanical properties of coal and rock.Currently, scientists have reached the same understanding that coal is damaged under the action of effective stress and the gas in coal can be quickly desorbed and the damaged coal can be discarded.Therefore, the outburst is the comprehensive effect of the dynamic distribution of geo-stress, gas in the coal seam, and the physical and mechanical properties of the coal seam itself.This has an impact on the gas extraction rate, desorption rate, gas pressure gradient, and mining efficiency.Coal and gas outburst is a mechanical failure process, and the outburst process can be divided into preparing, forming, developing, and ceasing processes as shown in Fig. 2. Li et al. stated that coal is subjected to tensile and shear failure under the action of in-situ stress and gas pressure 32 .In front of the working face, the coal seam is further crushed and destroyed under the action of the leading stress peak in Fig. 2a, which provides the necessary solid material base for the subsequent development of the outburst.The flow velocity of the gas in the pulverized coal is accelerated and the gas can accumulate, providing the necessary kinetic energy base for subsequent development.The dual-pore medium model of gas migration in coal is shown in Fig. 2b.Man-made disturbances such as blasting, tunneling, and pneumatic picking of coal cause a sharp change in the stress state of the coal and rock mass, and the pore pressure increases.This promotes the desorption of adsorbed gas and diffusion into the pores.As shown in Fig. 2c, when the coal seam in the outburst structure is suddenly exposed, a large amount of coal instantly rushes to the working face under the load of the gas, as shown in Fig. 2d.Therefore, the direct contribution of geological tectonic movement to outburst is to form a tectonic coal seam, create a high-stress environment conducive to outburst, and provide a geological structural environment conducive to gas preservation and outburst initiation.The conditions at the outburst site vary widely, but the outburst site must have a geological structural environment conducive to energy accumulation and sudden release.The coal seam in the working face underwent a gradual increase from the initial original rock stress state to the peak stress state.With the occurrence of the unloading confining pressure, the peak stress exceeded the strength limit of the gas-bearing coal seam failure and seepage characteristics under stress perturbation, and the stress decreased to the residual stress state.In a gas-rich coal seam, the gas flow process in the coal is directly related to the stress state.In particular, the permeability is sensitive to the change in stress.As shown in Fig. 2a, the stress state of the coal seam went through a spatiotemporal evolution process, and the permeability of the coal seam changed correspondingly, and the gas flow showed a spatiotemporal characteristic as well.Coal permeabilities are influenced by gas pressure, adsorptive constants, porosity, and stress conditions, which are represented by changes in permeability gradients, adsorptive expansion characteristics of coal, and porosity loading.The effective stress formula is based on the Terzaghi principle and is defined as: where: σ ij is effective stress; σ ij ′ is external stress; p is pore pressure; α is effective stress coefficient.
The effective stress of the gas-bearing coal, taking into account the adsorption expansion of the coal, is expressed as: where: σ s is gas stress; α ′ is the pore pressure coefficient.
Gas adsorption on coal reduces the surface energy of the coal, reducing the force between coal molecules and causing coal volume expansion, while desorption is the reverse process.The adsorption process affects the pore structure and the desorption process influences the pore structure.Macroscopically, coal exhibits adsorption expansion and desorption contraction effects.The expansion stress of coal adsorbing gas is given as: where: a and b are Langmuir adsorption constants respectively; R is the universal constant of gas, 8.3143; T is the absolute temperature, K; µ is Poisson's ratio, p is adsorption pressure; V m is the molar volume, 2.24 × 10 -3 m 3 /mol.Substituting (3) into formula (2) yields: (1)

Tunneling Workface
Primitive With the redistribution of coal stress, the mechanical properties of coal and rock change, the increase or decrease in permeability will cause local variations.This will result in uneven distribution of the pressure gradient and uneven gas release during gas seepage.The permeability increases exponentially, the gas in the coal seam is rapidly released and a very high gas pressure gradient is created within the coal seam.This inevitably promotes coal and gas outburst.

Energy conversion theory of coal and gas outburst
Although the duration of an outburst is short, it still has its phase.From an energetic perspective, it generally has four stages as shown in Fig. 3. Phase I: The preparation stage of the outburst mainly completes the energy accumulation (elastic strain energy E e , gas expansion energy E p , and internal energy of coal U ), and the direct cause of the outburst is related to the energy mutation.Phase II: The formation stage is the sudden instability, destruction, and ejection of the main structural coal seam, resulting in the exposure of the outburst coal seam.At this stage, it often manifests itself as sonic portents, such as the sound of machine guns, muffled thunder, and the sudden rise and fall of gas emissions.Phase III: The developing stage is the process from the exposure of the outburst coal to the end of the outburst.When it encounters a soft-hard coal junction or a new structure formed due to coal dust accumulation and hole clogging, the phenomenon of outburst deceleration, suspension, reignition, and initiation occurs.Phase VI: The ceasing stage is a point in time at the end of the outburst development.When the elastic strain energy, gas expansion energy, and internal energy of the coal and rock mass are less than the energy required for the outburst process (crushing work W 1 , throwing work W 2 , and energy loss of the gas flow field E ), it will stop because the mechanical and energetic conditions are no longer favorable.
The outburst process is subject to the law of conservation of energy.During the occurrence and development of coal and gas outburst, the energy sources are elastic strain energy, gas expansion energy and coal internal energy, and crushing work, throwing work and gas flow field loss energy are the energy required for outburst.Therefore, the physical formula of the outburst process can be expressed as follows: where: E e is the elastic strain energy stored in the coal when stressed and is calculated as follows: (4) Mining operation changes the mechanical properties and stress state of coal and rock mass.
In-situ stress and mining stress are superimposed, forming a large plastic zone.The pressure in the hole is not enough to push the broken coal and rock, and the hole tends to be stable.where: σ 1 , σ 2 , σ 1 are stress in three directions, E is elastic modulus, MPa; ρ c is the density of coal, kg/m 3 .E P is the expansion energy of the gas, and the calculation formula is: where: p gas refers to the gas pressure under standard conditions, which is atmospheric pressure, 0.1 MPa; V 0 is the gas volume involved in the outburst process, including desorbed gas content and free gas content, m 3 /t; γ is a multinomial index.Judging from most outburst cases, the outburst process is close to the adiabatic process, and the approximate calculation of γ can be 1.25.
Coal internal energy is coal seam internal energy, and coal seam internal energy change can be calculated by formula (10) when an outburst occurs and coal seam temperature decreases: where: ϑ is the non-uniformity coefficient; c is the specific heat capacity of coal, 0.79 kJ/(kg•K); m is the quality of coal, kg/m 3 ; T is the falling temperature.
The formulae for calculating the crushing and throwing work per unit mass, based on the principle of increasing energy consumption per unit area, are as follows: where: A is the energy consumed to increase the unit surface area is 505 J/m 2 ; d 0 is the average particle size of the initial coal sample; d i is the average particle size of coal at a specific position; ω is the non-uniformity coef- ficient, 1.2 ~ 1.7.
The energy lost in the gas flow field is the kinetic energy of the gas, which corresponds to the energy dissipated after the gaseous and solid phases of the coal have separated.Theoretically, the formula is the following: where: ρ ′ and v ′ is the density and velocity of gas flow when the gas phase and solid phase are separated.
Through the simulation test of coal and gas outburst in the laboratory, the calculation results of the energy in the outburst process are obtained.Of the total energy, 31% is used for the crushing work and 17% for the throwing work, and the loss of the gas flow field can reach 52%.A large amount of gas is released immediately.The energy change in the coal is used to promote the desorption behavior of the gas.The process of gas expansion and work transforms the internal energy of coal, so it can be considered that gas expansion energy E P and coal internal energy U are gas-related energy, accounting for over 90% of total energy 13 .The gas-related energy is much greater than the elastic energy E e of coal, and the outburst is dominated by the energy associated with gas.

Selection of predictive quantitative indicators
According to the mechanism of coal and gas outburst, eight quantitative parameters of outburst prediction can be summarized as shown in Fig. 4, including gas content, gas pressure, the initial velocity of gas emission from boreholes, index of the initial velocity of diffusion of coal gas, coefficient of firmness of coal, the thickness of the coal seam, maximum drilling cuttings volume, and distance from the geological structure.The relationship of each parameter to coal and gas outburst is explained.From Eq. ( 4), it is clear that when the coal seam is exposed or the coal seam in the working face is exposed, the gas pressure is reduced and the gas in the coal seam is desorbed and liberated, the expansion stress and pore pressure are reduced and the coal seam contracts, which also causes the effective stress to varying significantly.As the porosity of coal is low and nanoscale pores are the main ones, the influence of gas on coal mechanical properties is mainly the combined effect of expansion stress and pore pressure caused by micro-surface adsorption of a large amount of gas.Therefore, the influence of effective stress on the micropore structure of coal is an important factor in causing coal destruction.The magnitude and variation of the initial gas emission velocity and the initial gas emission velocity can characterize the internal stress variation of the coal seam.Gas pressure and gas content are important parameters affecting the internal energy and expansion energy of the gas.The coefficient of firmness of coal, the thickness of the coal seam, and maximum drilling cuttings volume can describe the internal energy of coal seam, while the distance from the geological structure is a parameter that characterizes the stress concentration area.Therefore, selecting a reasonable set of eight indicators related to prominence is important for subsequent machine learning.( 6) Vol.:(0123456789)

Theoretical foundations of the algorithms
This section may be divided by subheadings.It should provide a concise and precise description of the experimental results, their interpretation, as well as the experimental conclusions that can be drawn.

Support vector machine
Support Vector Machine (SVM) was proposed by Vapnik in 1995.It is based on statistical learning to solve the problems of small sample size, nonlinearity and high dimensional identification, and is an optimization algorithm of convex quadratic programming 33 .Constructing an optimal hyperplane with maximum interval to achieve optimal classification between data is the core idea of SVM 25 .For example, a known training set, Constructing a hyperplane to best divide the training set into the maximum interval.The hyperplane is called the prediction hyperplane and is expressed as follows: where: y is the prediction function; ω is the weight; b is biased.However, some data cannot be hyperplane-predicted and the linearly separable SVM model fails due to the uneven distribution of the sample data sets.To solve the above problems, relaxation variables ξ i are introduced as boundary conditions to reflect the situation that some sample points are misclassified.At the same time, to avoid a too strong influence on the classification results, a penalty factor C is introduced to compensate for this influence.ξ i n i=1 ξ i is defined as the misclassified part of the data set, and the linear constraint equation and constraint conditions are: The Lagrange function is introduced to transform the equation into its dual form, and the solution of formula ( 12) can be obtained as follows: where: K(x, x i ) is the kernel function; K(x, x i ) = exp −g|x − x i | 2 , the Gaussian radial basis function with a wider application range is selected as the kernel function.The parameter g of kernel function determines the influence degree of kernel function on the support vector machine.
Therefore, the penalty factor C and the parameter g of the kernel function are the main parameters that affect the prediction performance of the SVM model.

Grey wolf optimizer
Grey Wolf Optimizer (GWO) is inspired by the grey wolf.It stimulates the leadership and hunting mechanism of grey wolves in the wild.Among them, α, β, δ, ω represent the leader wolf, the auxiliary wolf, and the subordinate wolves of the third and fourth classes respectively.In addition, the three main steps of the hunt are realized: finding the prey, surrounding the prey, and attacking the prey.
Gas content, X 1 Gas pressure, X 2 Index of initial velocity of diffusion of coal gas, X 3 Thickness of coal seam, X 6 Strain Initial velocity of gas emission from boreholes, X 8 Maximum drilling cuttings volume, X 7 Distance from geological structure, www.nature.com/scientificreports/Finding the prey In the process of hunting, the distance between the grey wolf and the prey is expressed as: Formula (15) represents the distance between the individual and the prey, and formula ( 15) is the position updating formula of the grey wolf.Where: t is the current iterative algebra; A and C are coefficient vectors; X p (t) and X(t) are the position vectors of the prey and the grey wolf respectively.The calculation formulas for A and C are as follows: where: a is the convergence factor, which linearly decreases from 2 to 0 with the number of iterations; r 1 and r 2 are a random vector, taking a random number between [0,1].

Surrounding the prey
The grey wolf can identify the position of the prey and surround them.When the grey wolf identifies the position of prey, it can guide the wolves to surround the prey under the guidance wolf.The mathematical model of grey wolf tracking prey position is described as follows: where: D α ,D β , and D δ are the distance between α , β and ω with other wolves respectively; and are the optimal solution positions of three wolves respectively; X α (t) , X β (t) , and X δ (t) are the optimal solution positions of three wolves respectively.;A 1 ,A 2 , A 3 , C 1 , C 2 , and C 3 are random vectors, X 1 (t) , X 2 (t) , and X 3 (t) are the moving direction and step size of grey wolves, X(t + 1) is the average distance between wolves.

Attacking the prey
When the prey stops moving, the grey wolf completes the hunt by attacking.For the simulation of approaching prey, the value is gradually reduced so that the fluctuation range is also reduced.In other words, as the value is linearly reduced from 2 to 0, the corresponding value also changes in the interval [− a, a].As shown in Fig. 5, the next position of the grey wolf can be anywhere between the current position and the prey position if the value of A is in the interval.If |A| < 1, the grey wolf will attack its prey (falling into the local optimum).If |A| > 1, the grey wolf separates from the prey in the hope of finding more suitable prey (global optimum).

Improved GWO
In GWO, α, β, and δ guide ω wolves into the search space to find the optimal solution domain, which would cause the GWO algorithm to converge slowly.Furthermore, since the convergence factor A decreases linearly with increasing iteration times, too small a value will gradually move the wolves away from the optimal solution position, leading to local optimization of the algorithm.To overcome the above problems, the tent chaotic map and DLH strategy are introduced.Chaotic mapping produces a chaotic sequence instead of a pseudorandom number generator.In the GWO algorithm, whether the initial population is uniformly distributed has a great influence on the performance of the GWO algorithm.In the case of a uniformly distributed population, after many iterations, it can achieve a larger search area and faster convergence speed than random distribution.The expression of tent chaotic mapping is as follows: where: x k and x k+1 are the chaotic sequences generated by the chaotic mapping, k = 0, 1, 2, … n.
The dimension of the learning-based hunting strategy (DLH) is to allow the wolves to learn from the wolves around them as they hunt.

Initialization stage
In this phase, wolves with a total number of n are randomly distributed in a given search area [ l j , u j ] according to Eq. ( 9), where D is the dimension of the problem, and the wolves form a matrix with N rows and D columns, which X i (t) is computed using the fitness function.

Movement stage
Each wolf learns the strategies of its neighbors.It becomes another candidate for a new position.In GWO, α, β, and δ are the three best wolves.With the linear decrease of the coefficient and the positions of X α , X β and X δ , the position of the surrounding prey is determined.Finally, the first candidate who moves to the new position of the grey wolf X i (t) is named X i−GWO (t + 1) .DLH strategy generates another candidate position X i−DLH,d (t + 1) for wolves, and the calculation formula is as follows: where: X n,d (t) is a wolf nearby; X r,d (t) is a random wolf in the population.
The distance calculation formula is as follows: where: R i (t) is the Euclidean distance between X i (t) and X i−GWO (t + 1) ; M i (t) is the distance of R i (t) of wolves adjacent to X i (t) ; D i is the Euclidean distance between X i (t) and X j (t).

Selection and update stage
In this stage, the better candidate is selected by comparing the fitness values of the sum of the two candidates.
To update the new position X i (t + 1) , if the selected candidate's fitness is less than X i (t) , it is updated with the selected candidate's position.Otherwise, it remains unchanged.Finally, after the process has been carried out on all the individuals, the number of iterations is increased by 1 until the iterative search has reached the pre-defined number of iterations.

Random forest
Random Forest is a commonly used method in machine learning for classifying data.Its model uses ensemble learning theory and the classification model uses a decision tree algorithm.Currently, an important problem in machine learning is overfitting.Random forest is an ensemble learning algorithm that belongs to the bagging type.By combining several weak classifiers, the final result is voted or averaged.This gives the result of the whole model a high accuracy and generalization performance.It can achieve good results mainly due to the "randomness" and "forest", the former making it anti-fitting and the latter making it more accurate.The specific steps of the random forest algorithm are illustrated in Fig. 6.

Prediction model of the coal and outburst building Highlight the forecast parameter set
To construct a vector X = (X 1 , X 2 , X 3 , X 4 , X 5 , X 6 , X 7 , X 8 ), select eight forecast indicators.The elements in the vector represent: X 1 is the gas pressure, X 2 is the gas content, X 3 is the index of the initial velocity of diffusion of coal gas, X 4 is the coefficient of firmness of coal, X 5 is the distance from geological structure, X 6 is the thickness of coal seams, X 7 is the maximum drilling cuttings volume, and X 8 is the initial velocity of gas emission from boreholes.Among X = (X 1 , X 2 , X 3 , X 4 , X 5 , X 6 , X 7 , X 8 ), X 1 is the gas pressure, X 2 is the gas content, X 3 is the index of the initial velocity of diffusion of coal gas, and X 8 is the initial velocity of gas emission from boreholes are variates; X 4 is the coefficient of firmness of coal, X 5 is the distance from geological structure, X 6 is the thickness of coal seams, and X 7 is the maximum drilling cuttings volume are constant in a specific period.The data used

Prediction model building
The prediction model is an improved Grey Wolf optimizer coupled with SVM and Random Forest dimensionality reduction, as shown in Fig. 7.The basic process of prediction is as follows: ( From the prediction results of the SVM model without parameter optimization, it can be seen that the prediction accuracy of the training set is 93.65% and the accuracy of the test set is only 80.95%.The prediction of the test sets in groups 5, 7, 13, and 20 is wrong.The findings show that after repeated training of the SVM, a good fitting effect is achieved on the training set data, but the accuracy of the test set is average, indicating that the non-optimized SVM model does not have sufficient predictive ability for unknown data and still needs to be optimized.At the same time, it can be seen that the model performs well on the training set, but poorly on the test set.The results are quite different, indicating that the model's prediction results have been over-fitted to some extent.Therefore, to improve the predictive ability of the model for unknown data, the GWO grey wolf optimization algorithm is introduced to optimize the C and g parameters and improve the model performance.

GWO-SVM
The original datasets of the training and test sets are left unchanged, with the number of grey wolf populations set to 20, the maximum number of iterations set to 30, and the range of variation of the parameters C and g set to [0.1, 100].Following the specific steps of the GWO, the corresponding mathematical model is constructed in MATLAB, and the predictive results from the operation are shown in Fig. 9.
From Fig. 9, the prediction accuracy of the optimized SVM model in the training set reached 96.83%.The prediction accuracy of the testing set reached 90.48%.The data of the seventh and eighth groups were wrong.It can be seen from the population fitness curve that in the process of population evolution iteration in Fig. 9c, the population fitness gradually increased from generation 0 to generation 7, and fluctuated after generation 8.The fitness decreased and then increased.The final average fitness value is stable between 90 and 95%, and the best fitness value is 95.24%.

IGWO-SVM
The GWO-SVM algorithm has some overfitting due to the large difference between the results of the training and test sets.Moreover, without feature dimensionality reduction and algorithm optimization, the highest accuracy of the test set obtained by the two models is only 90.48%, so there is still room for improvement in the ability to predict unknown data.Tent chaotic mapping and DLH strategy are introduced to improve the GWO algorithm.Through the improved GWO model, the optimal C and g parameters for SVM for outbreak prediction are found, and the results are shown in Fig. 10.
From the results, it can be seen that IGWO has a good effect in optimizing the SVM model, and the accuracy of the training set is 100%, and the accuracy of the test set is 100%, which not only improves the accuracy of the test set, but also narrows the gap between the prediction accuracy of the test set and the training set, solves the over-fitting problem, and improves the accuracy of the test set.From the fitness curve, the fitness of the population gradually increased from Generation 0 to Generation 10.There was no fluctuation in that the fitness suddenly decreased and then increased until the evolution to generation 30.The final average fitness was stable between 0.95 and 1, and the best fitness was 1.Therefore, the effectiveness of the improved GWO by the Tent Chaotic Mapping and the DLH strategy is verified.The initial velocity of gas emission from boreholes, L/min 0.0801 (2) The Grey Wolf Optimizer Algorithm (GWO) can effectively improve the prediction performance of the SVM model.GWO-SVM model has improved the accuracy of the test set to 90.48% from the original non-optimized 80.95%, which shows that the Grey Wolf optimization algorithm can better optimize the parameters C and g of SVM, and the introduction of this optimization algorithm can improve the performance of SVM greatly.(3) The prediction performance of SVM can be further improved by the Improved Grey Wolf Optimize Algorithm (IGWO).It can simplify the data and achieve the same prediction results after introducing the dimension reduction method of Random Forest.The IGWO is introduced to enhance the SVM model, and the prediction result is improved, reaching 100%, which proves the prediction performance of the IGWO-SVM model, and the optimized model can improve the original model's overfitting problem.(4) After the introduction of the Random Forest to reduce the dimension of the features, only three parameters with high importance are kept, which makes the data composition easier.When the data of the remaining features are missing, the model can still achieve better prediction performance.This can better meet the needs of predicting coal mine gas outburst areas and has certain practicability.This algorithm is based on an intelligent algorithm for risk identification, monitoring, and early warning of coal mine gas outburst disaster can provide a technical way for online monitoring, identifying, predicting, and early warning of hidden dangers of coal mine gas outburst disaster.
There are still some areas where our research needs to be improved.Due to the limitation of statistical ability, we can't guarantee that all the causes of coal and gas explosion risks are included.For example, some stress-related indicators such as energy amplitude and frequency, duration, and frequency band energy intensity obtained by micro-seismic monitoring technology, acoustic emission technology, and electromagnetic radiation technology are not included in the attribute list of this study.This may affect the results of risk prediction.In the future, the performance of the model will be further improved by adding more scientific dynamic parameter data to the research (Supplementary information).

Figure 1 .
Figure 1.Statistical analysis of coal mine accidents and coal and gas outburst accidents in China from 2012 to 2022.(a) the data of Chinese coal mine accidents and million-ton mortality; (b) the data of gas outburst accidents and annual total coal output; (c) the variations in the number of kinds of coal and gas outburst accidents; (d) the variations in the death toll of kinds of coal and gas outburst accidents.

Figure 2 .
Figure 2. Mechanism of coal and gas outburst.
and rock mass are destroyed, forming initial holes and developing to the deep part of plastic zone.

Figure 3 .
Figure 3. Description of the mechanical action process of coal and gas outburst.

Figure 6 .
Figure 6.Schematic diagram of the Random Forest algorithm.

Figure 7 .
Figure 7. Diagram of processes of the prediction models.

Figure 8 .
Figure 8. Predictive results of the SVM model.(a) Predictive results of the training set samples; (b) predictive results of the test set samples.

Figure 9 .
Figure 9. Predictive results of the GWO-SVM model.(a) Predictive results of the training set samples; (b) predictive results of the test set samples; (c) fitness curve of the GWO-SVM model.

Figure 10 .
Figure 10.Predictive results of the IGWO-SVM model.(a) Predictive results of the training set samples; (b) predictive results of the test set samples; (c) fitness curve of the IGWO-SVM model.

Figure 12 .
Figure 12.Diagram of results of the prediction models with 3 predictive parameters.

Table 1 .
Normalized partial data set. is from the measured data of the 8#, 11#, and 13# coal seams in the Pansan Mining Area of Huainan Basin, China.In addition, each group of data is divided into outburst danger areas according to the "Prediction Conclusion of Outburst Danger in Outburst Coal Seam Areas in Pansan Mining Area", where 1 represents an area with outburst risk area and 0 represents an area without outburst danger.The data sets are divided into two categories to form an 84 × 9 data matrix.Each data set represents a different underground site as shown in Table1, including 8 types of predictive index values measured at the site and the blast hazard in the area where the site is located。

Table 2 .
The results of the feature important value.

Table 3 .
Comparison of results of multiple models.