Research on the Method of Methane Emission Prediction Using Improved Grey Radial Basis Function Neural Network Model

: E ﬀ ectively avoiding methane accidents is vital to the security of manufacturing minerals. Coal mine methane accidents are often caused by a methane concentration overrun, and accurately predicting methane emission quantity in a coal mine is key to solving this problem. To maintain the concentration of methane in a secure range, grey theory and neural network model are increasingly used to critically forecasting methane emission quantity in coal mines. A limitation of the grey neural network model is that researchers have merely combined the conventional neural network and grey theory. To enhance the accuracy of prediction, a modiﬁed grey GM (1,1) and radial basis function (RBF) neural network model is proposed, which combines the amended grey GM (1,1) model and RBF neural network model. In this article, the proposed model is put into a simulation experiment, which is built based on Matlab software (MathWorks.Inc, Natick, Masezius, U.S). Ultimately, the conclusion of the simulation experiment veriﬁed that the modiﬁed grey GM (1,1) and RBF neural network model not only boosts the precision of prediction, but also restricts relative error in a minimum range. This shows that the modiﬁed grey GM (1,1) and RBF neural network model can make more e ﬀ ective and precise predict the predicts, compared to the grey GM (1,1) model and RBF neural network model.


Introduction
Throughout the years, coal mines have played a crucial role in energy resources around the world. Ranking first in coal production, China is the largest consumer of coal all over the world-consuming roughly half of the world's coal consumption. Coal exploitation has made an important contribution to developing China's economy by providing employment opportunities and generating purchasing power. Coal still plays a dominant role in the total primary energy, even though coal production is declining in China [1,2]. Nevertheless, the production system in underground coal mines is a peculiarly complex disaster system, which causes coal mine accidents to frequently occur [3].
Safe production in underground coal mines is threatened by a wide variety of coal mine disasters. In China, high rates of coal mine accidents cause huge losses to the Chinese economy, and are responsible for thousands of fatalities among coal miners every year [4]. Among various types of coal mine accidents, methane accidents are normally perceived as the most hazardous [5]. After CO 2 , methane is the second-biggest source of global greenhouse methane emission; and a main source of methane emission is coal mines [6][7][8]. The energy contained in methane is in direct proportion to its risk degree. For this reason, it is becoming more popular in the energy sector [19]. Indeed, abundant methane released during mining may create security problems, but coal mine enterprises can also obtain obvious economic benefits if they use coal mine methane rationally [20]. On account of the low concentration of methane, some researchers focused on technologies of treating and recovering methane effectively [21]. Meanwhile, in contrast to traditional fossil fuel energy, methane is a cleaner energy resource.
China is the largest coal producer all over the world, in excess of 95% of the national coal production from the underground coal mine. Furthermore, comparing surface mines, underground coal mines have higher methane emissions. An enormous amount of methane is released during the mining process, especially from state-owned coal mines, which have high-methane content. An Between 1950-2015, extremely serious methane explosion disasters have frequently occurred in China. In this time, the number of extremely serious methane explosion disasters reached 126 [18]. Hence, it is necessary to carry out an investigation that can accurately predict methane outbursts, and provide guides for coal mine support.
The energy contained in methane is in direct proportion to its risk degree. For this reason, it is becoming more popular in the energy sector [19]. Indeed, abundant methane released during mining may create security problems, but coal mine enterprises can also obtain obvious economic benefits if they use coal mine methane rationally [20]. On account of the low concentration of methane, some researchers focused on technologies of treating and recovering methane effectively [21]. Meanwhile, in contrast to traditional fossil fuel energy, methane is a cleaner energy resource.
China is the largest coal producer all over the world, in excess of 95% of the national coal production from the underground coal mine. Furthermore, comparing surface mines, underground coal mines have higher methane emissions. An enormous amount of methane is released during the mining process, especially from state-owned coal mines, which have high-methane content. An estimated amount of 19 billion cubic meters of methane discharged from China's coal mine every year, leading the world, but a few pockets of methane were utilized [11,[22][23][24][25]. Methane is categorized as an important greenhouse gas. To protect the atmosphere, it is important that methane is extracted effectively [26]. However, reality dictates that during the process of mining, a considerable amount of coal mine methane just be discharged into the atmosphere through the ventilation system, meaning some of the unused extracted methane is also released [9]. This leads to global warming and energy waste. Previous research has shown that methane concentration that exceeds the limit is commonly triggered by high methane emission. In addition, excessive methane concentration will lead to many serious consequences, such as fire and explosion [27]. Coal mine accidents that are caused by methane are extremely dangerous and increase of depth of the exploited seam-resulting in an increase of methane content. Predicting methane emission is crucial to imposing restrictions on transfinite methane concentration [28].
In spite of methane being present in nearly all coal mines in China, it is unrealistic to use other energy in lieu of coal in recent years. Mining excavations not only affect economic development, but are also inseparable from the people's livelihood. According to statistics, coal account for 70% of the gross energy supply in China [29]. To guarantee a regular national economy, declining the methane emission and improving the utility ratio of extracted methane are appropriate.
Methane emission is admittedly considered as a significant parameter of methane accidents, and has dramatically impacted the mining process and the mine ventilation design [30]. Hence, to achieve safety production, it is necessary to analyze data in terms of methane emission, and forecast its change tendency in advance.
For the sake of processing complex and enormous data connected with methane emission, various forecasting approaches are proposed. Computational fluid dynamic (CFD) analysis is used to simulate methane emission in underground coal mines [31].
As early as the 1996s, researchers have suggested methods to forecast methane emissions that is associated with disturbed and undisturbed longwall faces [32]. With the rapid development of computer technology in recent years, numerical simulations were widely used to predict methane emissions.
In order to predict methane emission quantity precisely, a time series analysis on the strength of the Gaussian process regression model has been used for methane emission prediction [33]. Moreover, study [34] considered the steeply inclined and extremely thick coal seams, a used new method (that employed numerical simulation), which was applied to forecasting methane emission quantity. The issue related to methane emissions in a multifarious geospatial context was analyzed in Reference [35].
Some researchers have used several kinds of mathematical models to predict and analyze methane emissions simultaneously. Wei et al. focused on methane emission prediction based on a grey prediction model, new information model, and a metabolism model, and compared the results of these three prediction models [36]. Jing et al. simultaneously used model GM (1,1) of grey, and one element linear regression forecasting methane emission, and concluded that the former has a greater precision by contrast [30].
In this study, one improved radial basis function (RBF) neural network model of methane emissions from underground coal mine of China's Shanxi Province is applied to forecast. Our analysis includes six main affecting factors of methane emissions (coal seam thickness, coal seam spacing, coal seam depth, coal seam methane content, daily production, and daily progress). By processing these factors, a tool used in prediction (associated with methane emission model) was built by the artificial neural network. Predictions of methane emission are presented in this paper, and the prediction accuracy is also discussed. Moreover, this paper compares predicted values with actual values.

RBF Neural Network
An RBF neural network is a model based on Cover which can use primary function approximating arbitrary function.
where X = (x 1 , x 2 , . . . , x n ) T ∈ R n , y i and i represent input vector, output vector and NO. i node, respectively. c i is the centre and i express NO. i primary function. w ij represent the weight between the hidden layer and output layer, and i, as well as j, express NO. i of hidden layer, and NO. j node of output layer. · is Euclidean norm. N is the number of centers and hidden layer node. φ(·) is primary function of hidden layer. To solve output vector y i , parameters in the Equation (1) should be confirmed which are c i , w ij and φ(·).
The center c i of RBF neural network can be determined by means of K-means algorithm. It is considering that the elect center can simplify network. Thus, the data that are selected as RBF neural network center should immensely reduce number of elements in the hidden layer.
Confirming radius δ i should guarantee the space of training sample, which can be contained by the set of acceptance region of RBF unit. The general method for determining c i and δ i is done by employing a K-means algorithm, as shown in Equation (2).
Then, we used a linear least square method adjust weight w ij between the hidden layer and output layer, through Equations (3) and (4).
It is a common practice to consider Gaussian function as φ(·), as shown in Equation (5).
where X = (x 1 , x 2 , . . . , x n ) T represent input vector and σ i is NO. i perceptive variate which defines width. In this equation c i express center and subscript i represent NO. i primary function. N is the number of nodes in hidden layer. Eventually, the output vector can be ciphered out though Equation (6)

Grey Prediction Model
The grey prediction model is to process the seemingly irregular data column into a regular data column, using a certain method. Therefore, the grey prediction model generates a new sequence and sequence model. Modified grey forecasting model GM (1,1) is a valid model of bating the randomness of sample data and enhancing accuracy in forecasting methane emission quantity. The first step is to acquire data of factors that can influence methane emission quantity. Using a computing logarithmic [10] is an eminent way to make data smoother, and improves the accuracy of prediction. This algorithm is shown below. There are x groups training samples and y groups predicted samples, compute logarithmic of these x + y groups of samples. Take methane emission quantity Q (0) , for example: Then, obtaining x + y groups samples, respectively, are (6) correspond six influence factors about methane emission quantity, including thickness of coal seam, interval of coal seam, depth of coal seam, concentration of coal seam, daily output, daily progress. Moreover, methane emission quantity homologous sample is x (0) (7).
Secondly, a grey predicted model will be built. The most familiar method to establish grey system model is accumulation.
Adding up x groups of training samples receive a list of data: . Then, we used this data as an input vector for the modified RBF neural network model.
By using grey forecasting model GM (1,1) Then, to solve the parameters in this equation, we utilize the least square method: where B and y n are built by Finally, to obtainx (1) (6),x (1) (7), the objective vector of modified RBF neural network model, we establish the GM (1,1) model:x The calculated objective vector and input vector were then input into the modified network model for training.

Modified RBF Neural Network
The improved RBF Neural Network Structure is shown in Figure 2. (1) Determine the width of the Gaussian function, δ, set a variable A (1) used to store the sum of heterogeneous outputs, and set a variable B (1) used to count the heterogeneous samples, where 1 is the number of classes. (2) With x 1 , y 1 ) , as the first sequence, a network is established, and there is only one hidden neuron node in the network. Let c 1 = x 1 , A(1) = y 1 , B(1) = 1, c 1 is the center of the hidden layer node, and the weight of the hidden layer node to the output layer is: (3) For the second sequence, x 2 , y 2 ) , calculate the distance between x 2 and c 1 : c 2 is the hidden node center. A neuron node is added to the hidden layer of step (2), and the weight from the node to the output layer is w 2 = A(2)/B(2). (4) The k sequence: (x k , y k ), k = 3, 4, · · ·, P, the number of nodes in the hidden layer of the network is M, in which the center of the node is c 1 , c 2 , · · ·, c M in order, and the distance between x k and c 1 , c 2 , · · ·, c M is calculated, in turn,  6) Add incentives x k = x k 0 , x k 1 , · · ·, x k n−1 , x k j as the input of node j at time k. (7) Calculation output. Let the output of the i node of the hidden layer be: The network output is: (8) Adjustment weight. The objective function is: Adjust the weights w, µ, and θ to obtain the increment factor by the gradient descent method, where: The iteration formula to get weights and biases is: Energies 2020, 13, 6112 7 of 15 Among them: θ = [θ 1 , θ 2 · ··, θ m ] T , hidden layer output R = [R 1 , R 2 , · · ·, R H ] T , w i represents column i of w, µ j represents column j of µ, d = y k − f (x k ), n, H, m, in turn, indicate the number of input layer, hidden layer and output layer nodes, η represents the learning rate, 0 < η < 1; α is the momentum factor and 0 < α < 1; k is the number of iterations. (9) Go to the next sample and repeat steps (6) to (9), until the error E = 1 2 k y k − f (x k ) 2 reaches the specified accuracy.
Energies 2020, 13, x FOR PEER REVIEW 7 of 15 Among them:  

Modeling and Prediction of Improved Grey RBF Neural Network Model
The modeling and prediction steps, and the structure of the improved grey RBF neural network model are shown in Figures 3 and 4. The modeling prediction steps are as follows: (1) Group x is the training sample and group y is the prediction sample. Obtaining by evaluating the logarithm of the sample x + y, where ( 0 ) (7 ) x in the y group is not processed, ( 0 ) , , (1), ( 2 ), , (1), ( 2 ), , , this data is listed as the input vector of the improved RBF neural network model; (3) The data column ( 0 ) , ,

Results and Discussion
In order to test the prediction effect of the model, this paper uses 300 sets of methane data from a coal mine in Shanxi, and 15 groups are listed in Table 1. Among them, represent the coal seam depth, coal seam thickness, coal seam methane content, coal seam spacing, daily progress, daily output, and methane emission, respectively. In order to ensure the prediction accuracy, a sensitivity analysis of the number of training samples is made, and the results show that the reasonable number of training samples is 275. The first 275 groups are taken as training samples, and the last 25 groups, i.e., 276 to 300 groups, are used as prediction samples. Using MATLAB 7.8, we wrote programs to build RBF neural network models, grey RBF neural network models, and improved grey RBF neural network models. These three models were used to make predictions, respectively. The predicted values were compared with the actual values, and the relative errors of the three models are calculated for analysis and comparison.   4 x /(m) 5 x /(m/d) 6 x /(t/d) 7 x /(m 3 /min)

Results and Discussion
In order to test the prediction effect of the model, this paper uses 300 sets of methane data from a coal mine in Shanxi, and 15 groups are listed in Table 1. Among them, x 1 , x 2 , x 3 , x 4 , x 5 , x 6 , x 7 , represent the coal seam depth, coal seam thickness, coal seam methane content, coal seam spacing, daily progress, daily output, and methane emission, respectively.
In order to ensure the prediction accuracy, a sensitivity analysis of the number of training samples is made, and the results show that the reasonable number of training samples is 275. The first 275 groups are taken as training samples, and the last 25 groups, i.e., 276 to 300 groups, are used as prediction samples. Using MATLAB 7.8, we wrote programs to build RBF neural network models, grey RBF neural network models, and improved grey RBF neural network models. These three models were used to make predictions, respectively. The predicted values were compared with the actual values, and the relative errors of the three models are calculated for analysis and comparison. Table 1. Methane sample data.

Sample
Serial Number

RBF Neural Network Model Prediction
The input vector of the RBF neural network model is the six main factors that affect the methane emission amount, and the output vector is the methane emission amount. The modeling and prediction steps are as follows: (1) Take the first 275 sets of training samples x 1 , x 2 , x 3 , x 4 , x 5 , x 6 as the input vector and x 7 as the target vector; (2) Input the input vector and the target vector into the RBF neural network model to learn and train; (3) The last 25 sets of prediction samples x 1 , x 2 , x 3 , x 4 , x 5 , x 6 are input into the trained RBF neural network model, and the prediction value is calculated; (4) Finally, we can calculate the relative error. Figure 5 is the comparison between the predicted and actual values of the model, and Figure 6 is the relative error of the model.

RBF Neural Network Model Prediction
The input vector of the RBF neural network model is the six main factors that affect the methane emission amount, and the output vector is the methane emission amount. The modeling and prediction steps are as follows: (1) Take the first 275 sets of training samples   6  5  4  3  2  1   ,  ,  ,  , , as the input vector and 7 x as the target vector; (2) Input the input vector and the target vector into the RBF neural network model to learn and train; (3) The last 25 sets of prediction samples are input into the trained RBF neural network model, and the prediction value is calculated; (4) Finally, we can calculate the relative error. Figure 5 is the comparison between the predicted and actual values of the model, and Figure 6 is the relative error of the model.

RBF Neural Network Model Prediction
The input vector of the RBF neural network model is the six main factors that affect the methane emission amount, and the output vector is the methane emission amount. The modeling and prediction steps are as follows: (1) Take the first 275 sets of training samples  6  5  4  3  2  1   ,  ,  ,  ,  ,  x  x  x  x  x  x are input into the trained RBF neural network model, and the prediction value is calculated; (4) Finally, we can calculate the relative error. Figure 5 is the comparison between the predicted and actual values of the model, and Figure 6 is the relative error of the model.

Grey RBF Neural Network Model Prediction
The grey RBF neural network model combines the GM (1,1) model with the RBF neural network model. The input layer and the output layer are six and two neuron nodes, respectively. The modeling and prediction steps are as follows: (1) The first 275 sets of training samples x 1 , x 2 , x 3 , x 4 , x 5 , x 6 accumulate to generate a sequence x (1) (1), x (1) (2), . . . , x (1) (6), which is used as the input vector of the RBF neural network model; (2) The first 275 sets of training samples x 1 , x 2 , x 3 , x 4 , x 5 , x 6 , x 7 are calculated using the GM (1,1) model to obtainx (1) (6),x (1) (7), and usex (1) (6),x (1) (7) as the target vectors of the RBF neural network model; (3) Input the input vector and the target vector into the RBF neural network model to learn and train; (4) A cumulative operation of the last 25 sets of prediction samples x 1 , x 2 , x 3 , x 4 , x 5 , x 6 to obtain the generated sequence; (5) The generated sequence in (4) is input into the RBF neural network model that has been trained, and the output values are:x (1) (6),x (1) (7); is the final predicted value; (7) Find the relative error. Figure 7 is the comparison between the predicted value and the actual value of the grey RBF neural network model, and Figure 8 is the relative error of the model. model. The input layer and the output layer are six and two neuron nodes, respectively. The modeling and prediction steps are as follows: (1) The first 275 sets of training samples accumulate to generate a sequence ) 6 ( ),..., 2 ( ), 1 ( , which is used as the input vector of the RBF neural network model; (2) The first 275 sets of training samples

Improved Grey RBF Neural Network Model Prediction
Follow the steps to build an improved grey RBF neural network model and make predictions. The modeling and prediction steps are as follows: (1) The first 275 groups of training samples

Improved Grey RBF Neural Network Model Prediction
Follow the steps to build an improved grey RBF neural network model and make predictions. The modeling and prediction steps are as follows: