Next Article in Journal
Robust Control Allocation for Space Inertial Sensor under Test Mass Release Phase with Overcritical Conditions
Previous Article in Journal
Optical Imaging, Optical Sensing and Devices
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Combined Prediction Model of Gas Concentration Based on Indicators Dynamic Optimization and Bi-LSTMs

1
School of Civil and Resource Engineering, University of Science and Technology Beijing, Beijing 100083, China
2
Key Laboratory of Ministry of Education for High-Efficient Mining and Safety of Metal, University of Science and Technology Beijing, Beijing 100083, China
3
State Key Laboratory of Coking Coal Exploitation and Comprehensive Utilization, Pingdingshan 467000, China
4
College of Safety Engineering, North China Institute of Science and Technology, Langfang 065201, China
*
Author to whom correspondence should be addressed.
Sensors 2023, 23(6), 2883; https://doi.org/10.3390/s23062883
Submission received: 2 February 2023 / Revised: 17 February 2023 / Accepted: 27 February 2023 / Published: 7 March 2023
(This article belongs to the Section Environmental Sensing)

Abstract

:
In order to accurately predict the gas concentration, find out the gas abnormal emission in advance, and take effective measures to reduce the gas concentration in time, this paper analyzes multivariate monitoring data and proposes a new dynamic combined prediction method of gas concentration. Spearman’s rank correlation coefficient is applied for the dynamic optimization of prediction indicators. The time series and spatial topology features of the optimized indicators are extracted and input into the combined prediction model of gas concentration based on indicators dynamic optimization and Bi-LSTMs (Bi-directional Long Short-term Memory), which can predict the gas concentration for the next 30 min. The results show that the other gas concentration, temperature, and humidity indicators are strongly correlated with the gas concentration to be predicted, and Spearman’s rank correlation coefficient is up to 0.92 at most. The average R2 of predicted value and real value is 0.965, and the average prediction efficiency R for gas abnormal or normal emission is 79.9%. Compared with the other models, the proposed dynamic optimized indicators combined model is more accurate, and the missing alarm of gas abnormal emission is significantly alleviated, which greatly improves the early alarming accuracy. It can assist the safety monitoring personnel in decision making and has certain significance to improve the safety production efficiency of coal mines.

1. Introduction

A gas disaster is one of the most serious disasters in the coal mining process. Gas explosions, coal and gas outbursts, and other gas disasters will not only cause serious loss of life and property but also seriously affect mine safety production [1,2,3,4]. On 27 November 2014, a major gas explosion accident occurred in Songlin Coal Mine in Guizhou Province of China, resulting in 11 deaths and 8 injuries. This accident was caused by the accumulated gas in the roadway encountering frictional sparks from the fan blades running. The keys to preventing such accidents are to accurately predict the gas concentration and take effective measures to reduce the gas concentration in time [5,6].
The safety monitoring system, as one of the six major systems of the mine, provides real-time monitoring of gas concentration, wind speed, and other indicators in the mining working face [7]. The safety monitoring system is generally connected to the mine power supply system, which will automatically cut off the power supply once any gas abnormal emission is detected [8]. However, the start-up and shutdown of the power supply system will seriously affect the productivity of the coal mine. If the gas concentration can be predicted, the relevant measures to reduce the gas concentration can be taken in advance [9]. The multivariate monitoring data in the safety monitoring system provides conditions for the gas concentration prediction. The dynamic and accurate prediction of gas concentration can provide a decision-making basis for safety monitoring personnel, which is also significant to improve the safety production efficiency of coal mines.
In recent years, scholars have conducted in-depth research on the gas concentration prediction and put forward various prediction methods. Traditional prediction methods include gray prediction [10], chaotic time series [11], ARIMA time series prediction [12], etc. These methods make full use of the time series characteristics of gas concentration, but the prediction accuracy needs to be further improved. In latest years, with the gradual application of the neural network algorithm, many scholars have used it in gas concentration prediction, including long short-term memory network [13,14,15,16], recurrent neural network [17,18], and random neural networks [19], etc. As a deep learning technology, neural network has strong adaptability in time series prediction. It can rely on its special structural units to effectively use historical series information to realize deep mining of potential correlations between data, thus improving the accuracy of gas concentration prediction [13,20]. To further improve the prediction accuracy, scholars have made further improvements in terms of input data and model combinations. Cheng et al. [21] proposed an Evolutionary attention-based temporal graph convolutional network (EAT-GCN) to simultaneously capture the spatial and temporal dependences. Wang et al. [22] proposed the LSTM-LightGBM model with a variable weight combination method based on residual assignment. Wu et al. [23] proposed a model that integrates the t-distributed stochastic neighbor embedding algorithm (t-SNE) and the support vector regression (SVR) algorithm. Meng et al. [24] proposed a novel prediction method that combines classical time series analysis with these deep learning models. Dey et al. [25] proposed the t-SNE_VAE_bi-LSTM prediction model that combines the t-SNE, VAE, and bi-LSTM networks. Xu et al. [26] constructed a new IWOA-LSTM-CEEMDAN prediction model based on the improved whale optimization algorithm (IWOA). The above methods have improved the prediction accuracy.
In order to further improve the prediction accuracy of gas concentration, this paper analyzes the multivariate monitoring data and proposes a new combined prediction model of gas concentration based on indicators dynamic optimization and Bi-LSTMs. The proposal of this new prediction model significantly improves the accuracy of gas concentration prediction and also has a certain significance to improve the safety production efficiency of coal mines.

2. Research and Theoretical Background

2.1. Data Sources

The research data in this paper are obtained from IJCRS’15 Data Challenge: Mining Data from Coal Mines, organized at the Knowledge Pit web platform [27,28,29,30]. These data are the multivariate monitoring data near the long wall of the fully mechanized mining face in Poland, mainly including gas concentration, wind speed, temperature, humidity, air pressure, working current of the fully mechanized mining machine, etc. The layout of various sensors is shown in Figure 1. In Figure 1, sensors descriptions are as follows: AN-anemometer [m/s] (including AN311, AN422, and AN423); TP-temperature [°C] (including TP1711 and TP1721, the sensor type of TP1711 is temperature THP, three-component sensor THP2/94, and the sensor type of TP1721 is temperature THP, three-component sensor THP2/93); RH-humidity [%RH] (including RH1712 and RH1722, the sensor type of RH1712 is temperature THP, three-component sensor THP2/94, and the sensor type of RH1722 is temperature THP, three-component sensor THP2/93); BA-barometer [hPa] (including BA1713 and BA1723); MM-methane meter [%CH4] (including MM252, MM261, MM262, MM263, MM264, MM256 and MM211); CM861-high concentration methane meter [%CH4], the range is 0–100; TC862-temperature inside the pipeline [°C], the range is 10–40; WM868-methane delivery calculated according to CM861, CR863, P, TC862 [m3/min], the range is 0–50; CR863-sensor for pressure difference on the methane drainage flange [Pa], the range is 0–250; P-pressure inside the methane drainage pipeline [kPa], the range is 0–110; AMP1-current of the left cutting head of the cutter loader [A]; AMP2-current of the right cutting head of the cutter loader [A]; DMP3-current of the left haulage in the cutter loader [A]; DMP4-current of the right haulage in the cutter loader [A]; AMP5-current of the hydraulic pump engine in the cutter loader [A]; F-driving direction, 1 = left, 0 = right; V-cutter loader speed [Hz], Vmin = 3 Hz, Vmax = 100 Hz. Herz values are then transformed into m/min, 100 Hz equal to about 20 m/min.
The fully mechanized mining machine moves along the long wall between MM262 and MM264. The larger the working current of the coal cutter, the higher efficient the coal cutting will be, which theoretically leads to more gas diffusion into the mined-out area, and the monitoring value of the gas sensor will increase accordingly. The three key sensors, MM263, MM264, and MM256, are located near the upper corner of the fully mechanized mining face, which is prone to gas accumulation. The three sensors can monitor the change of gas concentration at this location in real-time, and the coal cutter will automatically shut down when the gas concentration of any sensors reaches the alarm threshold (the alarm threshold value of the three sensors is 1.0%). If the gas abnormal emissions can be predicted in advance, it is possible to reduce the coal cutter’s cutting speed, i.e., to give more time for the gas to diffuse before the entire line is stopped. Therefore, it is particularly important to predict the gas concentration of MM263, MM264, and MM256.

2.2. Prediction Theories

2.2.1. Spearman’s Rank Correlation Coefficient

In order to accurately predict the gas concentration near the fully mechanized mining machine, the Spearman’s rank correlation coefficient (SRC) is used to dynamically optimize prediction indicators. Spearman’s rank correlation coefficient (SRC) refers to the degree of correlation between indicators. The larger the SRC, the stronger the correlation between the two indicators, indicating that the indicator is useful for predicting the other indicator [31]. SRC is not affected by the overall distribution and sample size of data and is more suitable for non-normal distributed coal mine time series data [32]. Suppose that the two variables are P and Q respectively, and the number of their elements is both N. Pi and Qi represent the ith (1 ≤ iN) value of P and Q, respectively, and the two variables are sorted (both in ascending or descending order) to obtain the two variable ranking sets p and q, where the elements pi and qi represent the ranking of Pi and Qi in P and Q, respectively. A ranked difference set d is obtained by subtracting the corresponding elements in the sets p and q, where di = piqi (1 ≤ iN), and the SRC between the two variables P and Q is the following:
ρ = 1 6 i = 1 N d i 2 N ( N 2 1 )
If the SRC is 0.8~1.0, there is an extremely strong correlation between variables; 0.6~0.8, strong correlation; 0.4~0.6, moderate correlation; 0.2~0.4, weak correlation; 0.0~0.2, very weak correlation or no correlation. When the SRC is a negative, there is negative correlation between variables.

2.2.2. Long Short-Term Memory Network

Bi-directional Long Short-term Memory (Bi-LSTM) and Long Short-Term Memory (LSTM) networks are applied to establish a combined prediction model of gas concentration. Long Short-Term Memory (LSTM) network is an improved recurrent neural network (RNN) [33,34]. Three logic control units, namely, input gate, output gate, and forget gate, are added to the basic structure of RNN to control the input and output of information flow and the state of the cell units by setting the weights at the edges where the neural network memory cells are connected to other parts. LSTM overcomes the shortcomings of RNN, which cannot deal with long-distance dependence and is prone to gradient disappearance and explosion. The basic structure of LSTM is shown in Figure 2.
Ct and ht represent the memory state and the hidden state of LSTM, respectively, and xt represents the input of the model, σ is a sigmoid activation function. Ft, it, and C ˜ t calculation formula are the following:
f t = σ ( W f [ h t 1 , x t ] + b f )
i t = σ ( W i [ h t 1 , x t ] + b i )
C ˜ t = tanh ( W c [ h t 1 , x t ] + b c )
O t = σ ( W o [ h t 1 , x t ] + b o )
h t = O t × tanh ( C t )
where tanh is the hyperbolic tangent activation function; Wf, Wi, Wc, and Wo are the weight matrices connected to the input vector xt and the previous hidden state ht−1 for each layer, respectively; bf, bi, bc and bo are the bias terms of each layer, respectively.

2.2.3. Bi-Directional Long Short-Term Memory Network

Bi-directional Long Short-Term Memory (Bi-LSTM) network [35,36] combines the information of the input sequence in both forward and backward directions based on the LSTM, i.e., the sequence is input into LSTM in forward and reverse order, respectively. Bi-LSTM makes the feature acquired at moment t to have information between the past and the future, which can effectively ensure the accuracy of time sequence prediction [37,38]. For the output at moment t, the forward LSTM layer has the information at moment t in the input sequence and the previous moments, while the backward LSTM layer has the information at moment t in the input sequence and the moments after. The structure of Bi-LSTM is shown in Figure 3.
It can be seen from Figure 3 that Bi-LSTM is composed of two LSTMs. When solving, the hidden vector ht−1 of the front term is calculated to generate a new hidden vector h t , and the hidden vector ht+1 of the rear term is also calculated to generate a new hidden vector h t . Combine the output results of the positive and negative input sequences to obtain the result Yt, and the calculation formula is as follows:
h t = LSTM ( x t , h t 1 )
h t = LSTM ( x t , h t 1 )
Y t = tanh ( W Y [ h t + h t ] + b Y )
where, WY is the weight matrix of each layer connected to the previous hidden state ht−1; bY is the bias term.

3. Methodology

3.1. Prediction Methods

Through the analysis of multivariate monitoring data, this paper proposes a new combined prediction method of gas concentration based on indicators dynamic optimization and Bi-LSTMs, as shown in Figure 4. First, the multivariate monitoring data is processed into time series data, and the SRC between the gas concentration to be predicted and other indicators within the time window length l are analyzed. Considering the SRC changes with time, the prediction indicators are dynamically optimized. Then, the time series and spatial topology features of the optimized indicators are extracted to build the feature matrix X, which is input into the Bi-LSTM single indicator model to obtain the prediction result Y. After that, each optimized single indicator prediction result is input into the LSTM combined model again and obtain the gas concentration prediction value Ypredict. Finally, the model is continuously evaluated and optimized to improve the applicability of the prediction model.

3.2. Prediction Model

3.2.1. Correlation Analysis

The safety monitoring system stores all kinds of monitoring data, but not all indicators are useful for gas concentration prediction, and some may even be counterproductive. Redundant indicators will prevent the model from finding potential laws of data, cause dimensional disasters, and reduce the prediction efficiency of the model [20]. SRC shows the degree of correlation between indicators. Take MM256 as an example, set the length of time window to 1 day, calculate the SRC of MM256 and other indicators, and select the top three indicators of the strongest correlation with MM256 for 7 consecutive days, as shown in Table 1.
As can be seen from Table 1, the top three indicators of the strongest correlation with MM256 are almost different every day. Most of the three indicators contain a gas indicator, that is, the other gas indicators are mostly useful for MM256 prediction. However, the indicator with strongest correlation is not necessarily the gas indicator. For example, on the 26th and 27th days, the indicator with strongest correlation is RH1712, and the SRC is more than 0.7, showing a strong correlation. The SRC between each indicator and MM256 also changes every day. In order to clarify the SRC between other indicators and MM256, the frequency of each indicator in the top three indicators is counted, as shown in Figure 5.
It can be seen from Figure 5 that the frequency of each indicator is different, but the frequency of gas indicators is higher than other indicators on average, which also shows that gas indicators are mostly useful for MM256 prediction. In addition, humidity indicators (RH1212, RH1722) and temperature indicators (TP1721, TP1711, TC862) also appear more frequently, indicating that MM256 is more influenced by temperature and humidity. Among the top three indicators of the strongest correlation, MM263, RH1712, and MM264 have the highest frequency. Therefore, take these three indicators as an example to calculate the SRC with time, as shown in Figure 6.
As can be seen from Figure 6, the SRC of each indicator fluctuates drastically with time, and the highest SRC of three indicators all reach above 0.9, which is an extremely strong correlation. The extreme differences are 1.53, 1.1, and 1.22, respectively, which indicates that each indicator’s prediction effect on MM256 is unstable and the effect changes with time. Therefore, it is difficult to apply fixed one or some indicators to achieve stable and accurate prediction.

3.2.2. Indicators Dynamic Optimization

From the correlation analysis in Section 3.2.1, it can be seen that gas concentration, humidity, and temperature indicators are highly correlated with MM256. However, the top three indicators are not fixed, and the SRC of each indicator fluctuates with time. Analyze the above reasons, the underground environment of the coal mine is complex, and the gas concentration is disturbed by various factors (such as geological, human, equipment factors, etc.). In a certain period of time, the influence of various factors on gas concentration is different. For a certain factor, its influence on the gas concentration is not fixed, so it is difficult to determine its relationship with the gas concentration. Similar results are obtained when analyzing MM263 and MM264 data. It is more important to consider the correlation between gas concentration and other indicators at the prediction moment. Therefore, a dynamic optimization method for prediction indicators based on SRC is proposed, as shown in Figure 7.
For the multivariate time series monitoring data, calculate the SRC of each indicator with the gas indicator to be predicted within the time window length l of moment tc, and obtain ρ = (ρ1, ρ2, …, ρn). According to the SRC, the indicators with moderate correlation and above are optimized, i.e., the SRC ρ′ = (ρ1, ρ2, …, ρm) is higher than 0.4, those indicators form the optimized indicators set. If there are less than 2 indicators in the optimized indicators set, i.e., except for the gas indicators to be predicted, no other indicators’ SRC is higher than 0.4, then the five indicators corresponding to the 5 largest SRC in ρ = (ρ1, ρ2, …, ρn), gas concentration, temperature and humidity indicators form the optimized indicators set. Among them, gas concentration indicators include MM252, MM261, MM262, MM263, MM264, MM256, MM211 and CM861, recorded as MM; temperature indicators include TP1721, TP1711 and TC862, recorded as T; humidity indicators include RH1212 and RH1722, recorded as H. This method not only considers the spatiotemporal relationship of indicators but also the change of correlation between other indicators and gas concentration with time.

3.2.3. Features Extraction

Feature extraction is a very critical step in prediction model design, and the quality of features directly affects the model prediction performance. Time series and spatial topology features are extracted from the optimized indicators.
(1)
Time series features
The monitoring data are typical time series data, and the gas time series is set as follows:
C ( t ) = { C 1 ( t 1 ) , C 2 ( t 2 ) , C 3 ( t 3 ) , , C n ( t n ) }
where, ti is time; Ci(ti) is the gas concentration at moment ti; n is the time length of gas time series.
Calculate the features of the gas time series C(t):
Gas concentration monitoring value Cc(tc) at the moment tc; First order difference value of gas concentration at the moment tc, Dc(tc) = Cc(tc) − Cc−1(tc−1); Statistical features of the gas concentration time series Cl(t) = {Ccl(tcl), …, Cc−2(tc−2), Cc−1(tc−1), Cc(tc)} with time window length l before moment tc, including 10 statistical features of maximum, average, root mean square, variance, standard deviation, dispersion coefficient, peak factor, skewness, kurtosis and range of Cl(t) [39].
(2)
Spatial topology features
During the production process of a coal mine, gas flows into the roadway from fallen coal, coal wall, and mined-out area, and the gas continuously diffuses in the roadway. According to the correlation analysis in Section 3.2.1, there is a strong correlation between different indicators at different monitoring points. Randomly select the time series monitoring data of different sensors with a time length of 1 day, as shown in Figure 8.
As shown in Figure 8, The monitoring values of different sensors are correlated in time and space, and their time series data show a similar trend, but the peak value shows a certain time lag, which is obviously related to sensor layout. Therefore, when predicting the gas concentration at a certain measuring point, in addition to calculating its own time series features, it also calculates the time series features of the optimized indicators and inputs them into the model as spatial topology features.

3.2.4. Gas Concentration Prediction Model

(1)
Bi-LSTM single indicator model
Indicators dynamic optimization is applied to predict gas concentration, and single indicator prediction model is established based on single indicator time series data. Set the single indicator time series as follows:
X ( t ) = { x 1 ( t 1 ) , x 2 ( t 2 ) , x 3 ( t 3 ) , , x n ( t n ) }
where ti is the time; xi(ti) is the monitoring value at the moment ti; n is the time length of the single indicator time series. According to the features extraction method in this paper, calculate the time series features of the sequence Xl(t) = {Xcl(tcl), …, Xc−2(tc−2), Xc−1(tc−1), Xc(tc)} with the following time window length l before moment tc:
X t c = [ x 1 t c , x 2 t c , x 3 t c , , x 12 t c ]
Then the features matrix with lag step of h and prediction step of p is the following:
X = [ X t c h p + p , , X t c p , X t c ] T
i.e.,
X = [ x 1 t c h p + p x 2 t c h p + p x 12 t c h p + p x 1 t c p x 2 t c p x 12 t c p x 1 t c x 2 t c x 12 t c ]
Different features have different dimensions and units, which will affect the model prediction accuracy. Therefore, firstly, the features matrix X is mapped between [0, 1] by using the method of MinMaxScaler_fit_transform to obtain the standardized feature matrix XN. Input XN into Bi-LSTM, and then input its output results into 3 layers DenseNet. Feature reuse is achieved through the connection of features on the channel, and the activation function RELU is used to de-linearize between DenseNet layers. The Dropout layer is set to discard neurons from the network with a probability of 20% to prevent overfitting of the model. Dense connections with weight sharing are used to filter out process noise and interference information. Different features are learned in a supervised manner to output predicted values. Finally, the maximum gas concentration Y within T time after the moment tc is obtained through MinMaxScaler_inverse_transform. T is defined as the prediction length. Bi-LSTM single indicator model network structure is shown in Figure 9.
(2)
LSTM dynamic optimized indicators combined model
Gas disasters are often the result of multiple factors coupling, and combined prediction can excavate more hidden information. The simple combination of indicators does not necessarily improve the prediction accuracy, so it is necessary to optimize and process the information from various indicators. Bi-LSTM single indicator model can predict the maximum value Y of gas concentration in time T after the moment tc. If Y is larger than the threshold, the gas emission is abnormal in time T, and corresponding measures should be taken at once. By applying the proposed dynamic optimization method, m prediction indicators at moment tc are optimized. According to Formulas (11)–(14), the features matrix X = [X1, X2, …, Xm] with the lag step of h and the prediction step of p are calculated respectively. The features matrix of each indicator is input as the Bi-LSTM single indicator model, and obtain the prediction result Y = [Y1, Y2, …, Ym].
Combined prediction is to combine the results of different prediction models to improve the prediction accuracy as much as possible [40]. In this paper, the prediction result Y = [Y1, Y2, …, Ym] is combined and standardized to obtain YN, input YN into the LSTM of the two layers in turn, and the Dropout layer is also set to discard neurons from the network with a 20% probability. DenseNet layer connection and weight sharing are adopted. The final prediction value Ypredict of gas concentration in T time after the moment tc is obtained through MinMaxScaler_inverse_transform. LSTM-optimized indicators combined model network structure are shown in Figure 10.

3.3. Model Evaluation

In order to quantitatively evaluate the performance of the prediction model, evaluation indexes are used to evaluate the prediction results. The errors are divided into longitudinal error and transverse error. Longitudinal error is used to analyze the long-term operation of system in amplitude, while transverse error is utilized to study the prediction performance in time delay [40]. Longitudinal error includes mean absolute error (MAE), root mean square error (RMSE), etc.; transverse error includes determination coefficient (R2), etc. The closer the longitudinal error is to 0, and the closer the transverse error is to 1, the better the prediction model performance is. These evaluation indexes are defined as follows:
M A E = 1 n i = 0 n | y i y ^ i |
R M S E = 1 n i = 0 n ( y i y ^ i ) 2
R 2 = 1 i = 1 n ( y i y ^ i ) 2 i = 1 n ( y i y ¯ ) 2 y ¯ = 1 n i = 1 n y i
where, n is the total number of samples; yi is the real value of gas concentration; y ^ i is the predicted value of gas concentration.
According to the predicted value of gas concentration and the threshold value (1.0%), it can be inferred whether the gas emission is abnormal. The gas emission situations are divided into the following two categories: gas abnormal emission and gas normal emission. Considering the influence of sample imbalance [41], missing alarm and false alarm comprehensively, the false alarm rate (FAR), missing alarm rate (MAR), and prediction efficiency (R) [39,42,43] are applied to objectively evaluate the prediction model for gas abnormal or normal emission. These evaluation indexes are defined as follows:
F A R = F P F P + T N
M A R = F N F N + T P
R = T P T P + F N F P F P + T N
where TP is the number of samples with both predicted and real values exceeding the threshold; FN is the number of samples with real values exceeding the threshold but not the predicted value, i.e., the number of missing alarm samples; FP is the number of samples with predicted value exceeds the threshold but not the real value, i.e., the number of false alarm samples; TN is the number of samples with neither predicted value, nor real value exceeds the threshold.

4. Results

4.1. Data Preprocessing

The data provided by the Knowledge Pit web platform includes a training set and a test set. The training set contains the monitoring values of 28 sensors, which are given in time order, with each type of sensor producing one monitoring value per second, but there is a time overlap in the training set. The test set has the same data format as the training set, and its time does not overlap, but it is generated in random order. Therefore, this paper only analyzes the training set data. The time window length of the training set is 10 min, and the sliding step length is 1 min. The training set overlapping data is processed into time series data, as shown in Figure 11.

4.2. Model Training

The time overlapping test set data is processed through the method as shown in Figure 11, and 36 days of multivariate monitoring data of fully mechanized mining face is obtained. Set the time window length l = 1 d, lag step length h = 5, prediction length T = 30 min, that is, according to the multivariate monitoring data with 1 day before the current moment, calculate the features matrix of optimized indicators with the length of 5, and predict gas abnormal or normal emission within 30 min after the current moment. The prediction step p = 30 min, that is, it is predicted once every 30 min.
For each optimized indicator, according to the feature extraction and sample construction methods in this paper, 1629 sample datasets between features matrix X and Y corresponding one by one are constructed with the help of the Python-Keras framework. The first 80% of the sample datasets are divided into the training set and the last 20% into the test set. Input the training set and test set into the Bi-LSTM single indicator model, respectively, and the model is continuously trained and optimized. In the model training process, the Nadam optimizer is used to train the network, Epochs = 50, batch_size = 32, and lr = 0.001. For the prediction result Y of the Bi-LSTM single indicator model, the sample datasets corresponding to Y and Ypredict are constructed, which are input into the LSTM optimized indicators combined model again. In order to avoid overfitting, 5-fold cross-validation is used to adjust the hyperparameters for automatic optimization. The model training is reflected by cross-entropy loss (Loss), which is defined as follows:
L o s s = 1 n i = 1 n [ y i log y ^ i + ( 1 y i ) log ( 1 y ^ i ) ]
For the three key sensors MM256, MM263, and MM264, the gas indicator to be predicted, all-indicators, specified indicators (gas concentration, temperature, and humidity indicators, i.e., MM, T, and H), dynamically optimized indicators are respectively established in the following four models: Bi-LSTM single indicator model, Bi-LSTMs all-indicators combined model, Bi-LSTMs specified indicators combined model, and Bi-LSTMs dynamic optimized indicators combined model. The changes of Loss with Epoch during the training of different models are shown in Figure 12.
As can be seen from Figure 12, with the increase in Epoch, the Loss fluctuates and decreases during the model training. When Epoch = 10, the Loss gradually tends to be flat, and the model convergence performance is better. The convergence of the single indicator model is slightly poor, with large fluctuations. The all-indicators model has inferior convergence performance than the specified indicators model and the optimized indicators model due to the introduction of more redundant indicators. The specified indicators model is better than the optimized indicators model and converges slightly faster for MM256, but when Epoch = 30, the convergence of the two models is consistent. The optimized indicators model converges significantly better than the specified indicators model for MM263 and MM264. Indicators dynamic optimization is helpful to improve the convergence performance of the model, which verifies the necessity of indicators dynamic optimization to determine the input indicators of the model. Meanwhile, the changes of Loss with Epoch show that the Bi-LSTMs dynamic optimized combined indicators model is better for gas concentration prediction, and the prediction results of which can reflect the gas abnormal or normal emission in the future.

4.3. Prediction Results

Using the trained Bi-LSTMs dynamic optimized combined indicators model, 325 groups of test set data are used for model testing, and the prediction results are shown in Figure 13. It can be seen from Figure 13 that the fluctuation trend of the predicted value is basically consistent with the real value, and the proposed model has good prediction performance for gas concentration. In the prediction of gas abnormal emission, there is no false alarm for the three key sensors, but there are a few missing alarms.

4.4. Model Comparison

In order to comprehensively evaluate the prediction performance of the model, the predicted and real value of gas concentration, gas emission prediction, and actual gas emission situation of different models are compared. According to Equations (15)–(17), the R2, MAE, and RMSE are applied to evaluate the model’s gas concentration prediction performance. According to the gas concentration prediction value, the threshold value (1.0%), and Equations (18)–(20), the FAR, MAR, and R are applied to evaluate the model’s gas emission prediction performance. The prediction results of different indicators models are shown in Figure 14 the gray area in Figure 14 indicates that the prediction of gas abnormal or normal emission events is accurate. The evaluation of prediction results is specifically shown in Table 2 and Table 3, and R2′, MAE′, RMSE′, FAR′, MAR′, and R′ represent the average values of the corresponding evaluation indexes.
Compare and analyze the prediction results in Figure 14, Table 2 and Table 3. As shown in Figure 14a, Bi-LSTM single indicator model has poor prediction performance for gas concentration, the prediction performance of MM256 and MM264 is slightly better than that of MM263. The average R2 of predicted value and real value is 0.667, and the average MAE and RMSE are 0.156 and 0.143, respectively. There is a large difference between the predicted value and the real value. The prediction performance of a single indicator on gas abnormal or normal emission is also poor, with a high missing alarm rate, among as high as 100.0% for MM263, and the average missing alarm rate is 67.2%, which is quite serious. At the same time, there are also a few false alarms, with an average false alarm rate of 1.2% and an average prediction efficiency R of 31.6%.
As shown in Figure 14b, the prediction performance of the Bi-LSTMs all indicators combined model for gas concentration is slightly better than that of the single indicator model, but the results are also unsatisfactory. The average R2 is 0.675, the average MAE and RMSE are 0.095 and 0.138, respectively, and the error between the predicted and real value is reduced compared with the single indicator model. However, the prediction performance of all indicators on gas abnormal or normal emission is inferior to that of a single indicator, and the missing alarm rate of MM263 is still 100.0%, and the average missing alarm rate is 70.8%, which is still serious instead of decreasing compared to the single indicator model. The false alarm rate is reduced, with an average of 0.4%, and the average prediction efficiency R of 28.7%, which is lower than that of the single indicator.
As shown in Figure 14c, the prediction performance of the Bi-LSTMs specified indicators combined model for gas concentration is similar to that of single indicator and all-indicators models. The average R2 is 0.687 and the average MAE and RMSE are 0.100 and 0.138, respectively. However, the prediction performance of gas abnormal or normal emission is better than that of single indicator and all-indicators models. The missing alarm rate of MM263 is reduced to 68.2%, and the average missing alarm rate is 58.8%, the missing alarm is still serious. Compared with all indicators, the average false alarm rate is still 0.4%. The average prediction efficiency R is 40.7%, which is higher than that of single indicator and all indicators models.
As shown in Figure 14d, the Bi-LSTMs dynamic optimized indicators combined model has a better prediction performance for gas concentration, R2 of the three sensors is all above 0.90, with an average value of 0.965, average MAE and RMSE is 0.039 and 0.046 respectively. Compared with the single indicator, all indicators, and specified indicators models, the average R2 is significantly improved. The prediction performance for gas abnormal or normal emission is also better than that of other models. It significantly reduces the missing alarm and eliminates false alarms, with a false alarm rate of 0.0% and an average missing alarm rate of 20.1%. The average prediction efficiency R is 79.9%, which is 152.8%, 178.4%, and 96.3% higher than that of applying the single indicator, all-indicators, or specified indicators models, respectively.

5. Discussion

Mining working faces are important production places of a coal mine, and safe and efficient coal mining is the key to coal mine production and operation. Gas concentration usually increases abnormally before the occurrence of gas disasters such as gas explosions and coal and gas outbursts. In order to prevent gas accidents and ensure safe production during coal mining, the Coal Mine Safety Regulations of China stipulate that the gas concentration in the coal mining working face shall not exceed 1%; otherwise, an alarm will be given. The setting of this threshold value has effectively reduced gas accidents. However, with the gradual increase in the mine depth, the ground stress, and gas content increase. When the gas concentration exceeds the threshold value, coal mining processes such as coal cutters will be interrupted, and the operators will also fall into panic, which seriously affects the safety production efficiency of the coal mine and the psychological health of the operators. Therefore, the dynamic and real-time prediction of the gas concentration can improve the working efficiency of the coal cutters and effectively ensure the safe and continuous production of the mining working face.
In the process of coal mining, gas flows into the roadway from the fallen coal, coal wall, and mined-out area, and the gas continuously diffuses in the roadway. In this process, the gas emission from the coal wall and adjacent coal seams is relatively stable, while the gas from the fallen coal and mined-out area is easy to accumulate, which is greatly influenced by the degree of coal fragmentation, coal cutter speed, wind speed, etc. The influences of this information are largely reflected in other types of sensors, such as gas, temperature, humidity, and air pressure sensors. However, not every sensor information is useful for gas concentration prediction, and the correlation between other sensors and gas concentration changes with time. Therefore, this paper proposes a dynamic optimization method for prediction indicators. This method not only considers the spatiotemporal relationship of indicators but also the change of correlation between other indicators and gas concentration with time. On this basis, establishes a combined prediction model of gas concentration based on indicators dynamic optimization and Bi-LSTMs, which further improves the accuracy of gas concentration prediction and its abnormal emission alarming. The prediction model is based on multivariate real-time monitoring data, which can realize real-time and dynamic prediction of gas concentration. However, the prediction model is more suitable for the prediction of gas or other indicators monitored by multiple sensors compared with a few sensors, which is also the limitation of the model. Admittedly, the alarming effect of the proposed prediction model is significantly better than that of the single indicator model, all-indicators combined model, and specified indicators combined model.
Analyzing the above reasons, the single indicator model only uses its own historical data for prediction, with less input information, which can only reflect the gas time features, but not the spatial topology features. Compared with the single indicator model, all-indicators combined model introduces more information for the model, but this information is not all useful for prediction, and some information may even have a negative effect, causing dimensional disasters. The specified indicators combined model considers the spatial topology features, and the prediction indicators are relatively better, but the quality of these indicators changes with time and is not optimal at the prediction moment. Therefore, compared with single indicator and all-indicators, the performance of the model is only slightly improved. The dynamic optimized indicators combined model optimizes strongly correlated indicators, and different moments correspond to different indicators, which ensures the optimal sample indicators at the prediction moment. Compared with a single indicator, all-indicators, and specified indicators, the model introduces the most useful information for prediction at the prediction moment. Therefore, indicators dynamic optimization significantly improves the performance of the model.
The prediction model established in this paper sets many adjustable parameters, including time window length l, lag step h, prediction length T, prediction step p, etc. In this paper, set l = 1 d, lag step h = 5, T = 30 min, and p = 30 min, calculate the features matrix of optimized indicators with the length of 5 based on the multivariate monitoring data with 1 d before the current moment, apply optimized indicators to predict the gas abnormal or normal emission within 30 min after the current moment and the prediction is made once every 30 min. Of course, these parameters can be adjusted according to the actual demand of the coal mine, but this may affect the prediction performance of the model. The influence of these parameters on model prediction performance will be also the focus of future research. In addition, this study is a “Point prediction” of gas concentration for three key gas sensors. We plan to expand the gas concentration prediction to “Two-dimension prediction”, or even “Three-dimension prediction”, so as to visualize the gas concentration prediction area in the mining face, which is also the key direction of future research. The research in this paper lays the foundation for future research.

6. Conclusions

Accurate prediction of gas concentration in the mining face and early detection of gas abnormal emissions are the keys to preventing gas accidents. In order to avoid the contingency of a single indicator early alarming method and capture the precursor information of gas disaster more accurately, this paper proposes a new combined prediction model of gas concentration based on indicators dynamic optimization and Bi-LSTMs, which is proved to be effective. The main conclusions are as follows:
(1)
The Spearman’s rank correlation between gas concentration and other indicators is analyzed. The gas concentration, temperature, and humidity indicators have a strong correlation with the gas concentration to be predicted, and Spearman’s rank correlation coefficient can reach 0.92 at most. For a certain indicator, its correlation with gas concentration fluctuates with time, with a maximum range of 1.53. The prediction effect of a single indicator on gas concentration is not stable, it is necessary to consider that the correlation changes with time at the prediction moment to dynamically optimize indicators;
(2)
A dynamic optimization method for prediction indicators based on Spearman’s rank correlation coefficient is proposed. This method not only considers the spatiotemporal relationship of indicators, but also the change of correlation between other indicators and gas concentration with time, reducing the training cost and model complexity, and improving the model convergence speed and training efficiency;
(3)
A new combined prediction model of gas concentration based on indicators dynamic optimization and Bi-LSTMs is established, which can predict the gas concentration for the next 30 min. Using the proposed model, the average R2 of predicted value and real value is 0.965, and the average prediction efficiency R for gas abnormal or normal emission is 79.9%. The proposal of this new prediction model significantly reduces the missing alarm and eliminates false alarms. Compared with other models, the gas concentration prediction is more accurate, and the missing alarm of gas abnormal emission is greatly alleviated, which greatly improves the alarming accuracy;
(4)
The research of this paper is a “Point prediction” of gas concentration for three key gas sensors. “Two-dimension prediction” or “Three-dimension prediction” for gas concentration could realize to visualize the gas concentration prediction area in the mining face, and it is also the key direction of future research. In addition, the influence of adjustable parameters of the model on prediction performance will be also the focus of future research. The research in this paper lays the foundation for future research.

Author Contributions

Conceptualization, Y.P., D.S. and H.W.; methodology, Y.P. and H.W.; software, Y.P.; validation, Y.P., D.S., H.W. and Q.L.; formal analysis, Y.P.; investigation, Y.P., D.S. and L.Q.; resources, H.W.; data curation, Y.P. and H.W.; writing—original draft preparation, Y.P.; writing—review and editing, D.S., L.Q., H.W., X.H. and Q.L.; visualization, Y.P.; supervision, D.S. and X.H.; project administration, D.S. and X.H.; funding acquisition, D.S. and L.Q. All authors have read and agreed to the published version of the manuscript.

Funding

This work was supported by the National Natural Science Foundation of China (Grant No. 52174162, 52004016).

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Data Availability Statement

The processed data cannot be shared at this time as the data also forms part of an ongoing study.

Acknowledgments

The research data in this paper are obtained from IJCRS’15 Data Challenge: Mining Data from Coal Mines, organized at Knowledge Pit web platform. Thanks to Knowledge Pit web platform.

Conflicts of Interest

The authors declare no competing financial interest or personal relationship that could have appeared to influence the work in this paper.

References

  1. Wang, E.; Zhang, G.; Zhang, C.; Li, Z. Research progress and prospect on theory and technology for coal and gas outburst control and prevention in China. J. China Coal Soc. 2022, 47, 297–322. [Google Scholar] [CrossRef]
  2. Yu, M.; Yang, X.; Zheng, K.; Luan, P. Progress and development of coal mine gas explosion suppression and disaster reduction technology in China. J. China Coal Soc. 2020, 45, 168–188. [Google Scholar] [CrossRef]
  3. Li, B.; Wang, E.; Shang, Z.; Li, Z.; Li, B.; Liu, X.; Wang, H.; Niu, Y.; Wu, Q.; Song, Y. Deep learning approach to coal and gas outburst recognition employing modified AE and EMR signal from empirical mode decomposition and time-frequency analysis. J. Nat. Gas Sci. Eng. 2021, 90, 103942. [Google Scholar] [CrossRef]
  4. Zhang, G.; Wang, E.; Ou, J.; Li, Z. Regional prediction of coal and gas outburst under uncertain conditions based on the spatial distribution of risk index. Nat. Resour. Res. 2022, 31, 3319–3339. [Google Scholar] [CrossRef]
  5. Chen, L.; Wang, E.; Feng, J.; Kong, X.; Li, X.; Zhang, Z. A dynamic gas emission prediction model at the heading face and its engineering application. J. Nat. Gas Sci. Eng. 2016, 30, 228–236. [Google Scholar] [CrossRef]
  6. Song, D.Z.; He, X.Q.; Qiu, L.M.; Zhao, Y.J.; Cheng, X.H.; Wang, A.H. Study on real time dynamic monitoring and early warning technology of regional and local outburst danger. Coal Sci. Technol. 2021, 49, 110–119. [Google Scholar] [CrossRef]
  7. Wang, G.F.; Ren, H.W.; Zhao, G.R.; Gong, S.X.; Du, Y.B.; Xue, Z.X.; Pang, Y.H.; Zhang, X. Digital model and giant system coupling technology system of smart coal mine. J. China Coal Soc. 2022, 47, 61–74. [Google Scholar] [CrossRef]
  8. Chen, H.Z.; Wang, Y. Design and realization of integrated emergency response system in safety monitoring system of coal mine based on GIS. Min. Res. Dev. 2019, 39, 121–125. [Google Scholar] [CrossRef]
  9. Ge, S.; Hao, S.; Zhang, S.; Zhang, X.F.; Zhang, L.; Wang, S.B.; Wang, Z.B.; Bao, J.S.; Yang, X.L.; Yang, J.J. Status of intelligent coal mining technology and potential key technologies in China. Coal Sci. Technol. 2020, 48, 28–46. [Google Scholar] [CrossRef]
  10. Lai, X.W.; Xia, Y.N.; Zheng, W.B.; Cui, J.F.; Wu, Y.Q.; Shi, Y.X. Improved grey prediction of gas concentration sequence based on integrated learning. J. Saf. Sci. Technol. 2021, 17, 16–21. [Google Scholar] [CrossRef]
  11. Dai, W.; Ji, C.P.; Wang, Y.J. The prediction of multi-scale entropy on gas concentration time series based on variational mode decomposition. J. Nanoelectron. Optoelectron. 2019, 14, 1313–1325. [Google Scholar] [CrossRef]
  12. Wang, P.; Wu, Y.P.; Wang, S.L.; Song, C.; Wu, X.M. Study on Lagrange-ARIMA real-time prediction model of mine gas concentration. Coal Sci. Technol. 2019, 47, 141–146. [Google Scholar] [CrossRef]
  13. Fu, H.; Liu, Y.Z.; Xu, N.; Zhang, J.N. Research on gas concentration prediction based on multi-sensor-deep long short-term memory network fusion. Chin. J. Sens. Actuators 2021, 34, 784–790. [Google Scholar] [CrossRef]
  14. Zhang, Z.; Zhu, Q.; Li, Q.; Liu, Y.; Zhang, E.H.; Zhao, Q.M.; Qing, X.F. Prediction of mine gas concentration in heading face based on Keras long short time memory network. Saf. Environ. Eng. 2021, 28, 61–67+78. [Google Scholar] [CrossRef]
  15. Zhang, T.; Song, S.; Li, S.; Ma, L.; Pan, S.; Han, L. Research on gas concentration prediction models based on LSTM multidimensional time series. Energies 2019, 12, 161. [Google Scholar] [CrossRef] [Green Version]
  16. Lyu, P.; Chen, N.; Mao, S.; Li, M. LSTM based encoder-decoder for short-term predictions of gas concentration using multi-sensor fusion. Process Saf. Environ. Prot. 2020, 137, 93–105. [Google Scholar] [CrossRef]
  17. Li, S.G.; Ma, L.; Pan, S.B.; Shi, X. Research on prediction model of gas concentration based on RNN in coal mining face. Coal Sci. Technol. 2020, 48, 33–38. [Google Scholar] [CrossRef]
  18. Song, S.; Li, S.; Zhang, T.; Ma, L.; Pan, S.; Gao, L. Research on a multi-parameter fusion prediction model of pressure relief gas concentration based on RNN. Energies 2021, 14, 1384. [Google Scholar] [CrossRef]
  19. Zhang, Y.; Guo, H.; Lu, Z.; Zhan, L.; Hung, P.C. Distributed gas concentration prediction with intelligent edge devices in coal mine. Eng. Appl. Artif. Intell. 2020, 92, 103643. [Google Scholar] [CrossRef]
  20. Wu, J.; Liu, G.; Wang, J.; Zuo, Y.; Bu, H.; Lin, H. Data Intelligence: Trends and Challenges. Syst. Eng.-Theory Pract. 2020, 40, 2116–2149. [Google Scholar] [CrossRef]
  21. Cheng, L.; Li, L.; Li, S.; Ran, S.; Zhang, Z.; Zhang, Y. Prediction of gas concentration evolution with evolutionary attention-based temporal graph convolutional network. Expert Syst. Appl. 2022, 200, 116944. [Google Scholar] [CrossRef]
  22. Wang, X.; Xu, N.; Meng, X.; Chang, H. Prediction of gas concentration based on LSTM-LightGBM variable weight combination model. Energies 2022, 15, 827. [Google Scholar] [CrossRef]
  23. Wu, H.; Shi, S.; Lu, Y.; Liu, Y.; Huang, W. Top corner gas concentration prediction using t-distributed Stochastic Neighbor Embedding and Support Vector Regression algorithms. Concurr. Comput.-Pract. Exp. 2020, 32, e5705. [Google Scholar] [CrossRef]
  24. Meng, X.R.; Chang, H.Q.; Wang, X.Q. Methane concentration prediction method based on deep learning and classical time series analysis. Energies 2022, 15, 2262. [Google Scholar] [CrossRef]
  25. Dey, P.; Saurabh, K.; Kumar, C.; Pandit, D.; Chaulya, S.K.; Ray, S.K.; Prasad, G.M.; Mandal, S.K. t-SNE and variational auto-encoder with a bi-LSTM neural network-based model for prediction of gas concentration in a sealed-off area of underground coal mines. Soft Comput. 2022, 25, 14183–14207. [Google Scholar] [CrossRef]
  26. Xu, N.; Wang, X.; Meng, X.; Chang, H. Gas concentration prediction based on IWOA-LSTM-CEEMDAN residual correction model. Sensors 2022, 22, 4412. [Google Scholar] [CrossRef] [PubMed]
  27. Janusz, A.; Sikora, M.; Wróbel, L.; Stawicki, S.; Grzegorowski, M.; Wojtas, P.; Ślęzak, D. Mining Data from Coal Mines: IJCRS’15 Data Challenge. Rough Sets Fuzzy Sets Data Min. Granul. Comput. 2015, 9437, 429–438. [Google Scholar] [CrossRef]
  28. Zagorecki, A. Prediction of methane outbreaks in coal mines from multivariate time series using random forest. Rough Sets Fuzzy Sets Data Min. Granul. Comput. 2015, 9437, 494–500. [Google Scholar] [CrossRef]
  29. Ślęzak, D.; Grzegorowski, M.; Janusz, A.; Kozielski, M.; Nguyen, S.H.; Sikora, M.; Stawicki, S.; Wróbel, Ł. A framework for learning and embedding multi-sensor forecasting models into a decision support system: A case study of methane concentration in coal mines. Inf. Sci. 2018, 451, 112–133. [Google Scholar] [CrossRef]
  30. Janusz, A.; Grzegorowski, M.; Michalak, M.; Wróbel, Ł.; Sikora, M.; Ślęzak, D. Predicting seismic events in coal mines based on underground sensor measurements. Eng. Appl. Artif. Intell. 2017, 64, 83–94. [Google Scholar] [CrossRef]
  31. Li, H.M.; Cao, Y.C.; Su, L.M. Pythagorean fuzzy multi-criteria decision-making approach based on Spearman rank correlation coefficient. Soft Comput. 2022, 26, 3001–3012. [Google Scholar] [CrossRef]
  32. Wang, D.L.; Su, A.X.; Liu, W.P. Trend analysis of vegetation cover changes based on Spearman Rank correlation coefficient. J. Appl. Sci. 2019, 37, 519–528. [Google Scholar] [CrossRef]
  33. Greff, K.; Srivastava, R.K.; Koutník, J.; Steunebrink, B.R.; Schmidhuber, J. LSTM: A Search Space Odyssey. IEEE Trans. Neural Netw. Learn. Syst. 2017, 28, 2222–2232. [Google Scholar] [CrossRef] [Green Version]
  34. Ma, Y.; Peng, H.; Khan, T.; Cambria, E.; Hussain, A. Sentic LSTM: A Hybrid Network for Targeted Aspect-Based Sentiment Analysis. Cogn. Comput. 2018, 10, 639–650. [Google Scholar] [CrossRef]
  35. Yin, J.; Deng, Z.; Ines, A.V.; Wu, J.; Rasu, E. Forecast of short-term daily reference evapotranspiration under limited meteorological variables using a hybrid bi-directional long short-term memory model (Bi-LSTM). Agric. Water Manag. 2020, 242, 106386. [Google Scholar] [CrossRef]
  36. Liang, T.; Zhao, Q.; Lv, Q.; Sun, H. A novel wind speed prediction strategy based on Bi-LSTM, MOOFADA and transfer learning for centralized control centers. Energy 2021, 230, 120904. [Google Scholar] [CrossRef]
  37. Zhou, Z.; Zhang, J.J.; Ding, H.H.; Li, F. Settlement prediction model of shield tunnel under-crossing existing tunnel based on GA-Bi-LSTM. Chin. J. Rock Mech. Eng. 2022, 1–11. [Google Scholar] [CrossRef]
  38. Guo, H.P.; Liu, B.; Xiao, Y. Short-term railway freight volume prediction based on Bi-LSTM network. Railw. Freight Transp. 2022, 40, 52–58. [Google Scholar] [CrossRef]
  39. Qiu, L.M.; Peng, Y.J.; Song, D.Z. Risk prediction of coal and gas outburst based on abnormal gas concentration in blasting driving face. Geofluids 2022, 2022, 3917846. [Google Scholar] [CrossRef]
  40. Ouyang, T.H.; Zha, X.M.; Qin, L. A combined multivariate model for wind power prediction. Energy Convers. Manag. 2017, 144, 361–373. [Google Scholar] [CrossRef]
  41. Wang, H.; Li, Z.; Song, D.; He, X.; Sobolev, A.; Khan, M. An intelligent rockburst prediction model based on scorecard methodology. Minerals 2021, 11, 1294. [Google Scholar] [CrossRef]
  42. Tian, X.H.; Li, Z.L.; Song, D.Z.; He, X.Q.; Jiao, B.; Cao, A.Y.; Ma, Y.K. Study on microseismic precursors and early warning methods of rockbursts in a working face. Chin. J. Rock Mech. Eng. 2020, 39, 2471–2482. [Google Scholar] [CrossRef]
  43. Xia, Y.X.; Kang, L.J.; Qi, Q.X.; Mao, D.B.; Ren, Y.; Lan, H.; Pan, J.F. Five indexes of microseismic and their application in rock burst forecastion. J. China Coal Soc. 2010, 35, 2011–2016. [Google Scholar] [CrossRef]
Figure 1. The layout of various sensors. The red arrows show the directions of air flow, and the blue arrows show the directions of gas flow.
Figure 1. The layout of various sensors. The red arrows show the directions of air flow, and the blue arrows show the directions of gas flow.
Sensors 23 02883 g001
Figure 2. The basic structure of LSTM.
Figure 2. The basic structure of LSTM.
Sensors 23 02883 g002
Figure 3. The structure of Bi-LSTM.
Figure 3. The structure of Bi-LSTM.
Sensors 23 02883 g003
Figure 4. Combined prediction method of gas concentration based on indicators dynamic optimization and Bi-LSTMs.
Figure 4. Combined prediction method of gas concentration based on indicators dynamic optimization and Bi-LSTMs.
Sensors 23 02883 g004
Figure 5. The frequency of each indicator in the top three indicators.
Figure 5. The frequency of each indicator in the top three indicators.
Sensors 23 02883 g005
Figure 6. The change of SRC with time.
Figure 6. The change of SRC with time.
Sensors 23 02883 g006
Figure 7. A dynamic optimization method for prediction indicators based on SRC.
Figure 7. A dynamic optimization method for prediction indicators based on SRC.
Sensors 23 02883 g007
Figure 8. Time series monitoring data of different sensors.
Figure 8. Time series monitoring data of different sensors.
Sensors 23 02883 g008
Figure 9. Bi-LSTM single indicator model network structure.
Figure 9. Bi-LSTM single indicator model network structure.
Sensors 23 02883 g009
Figure 10. LSTM-optimized indicators combined model network structure.
Figure 10. LSTM-optimized indicators combined model network structure.
Sensors 23 02883 g010
Figure 11. Time series data processing.
Figure 11. Time series data processing.
Sensors 23 02883 g011
Figure 12. The changes of Loss with Epoch.
Figure 12. The changes of Loss with Epoch.
Sensors 23 02883 g012
Figure 13. Prediction results of Bi-LSTMs dynamic optimized combined indicators model.
Figure 13. Prediction results of Bi-LSTMs dynamic optimized combined indicators model.
Sensors 23 02883 g013
Figure 14. Prediction results of different indicators models. (a) Bi-LSTM single indicator model. (b) Bi-LSTMs all-indicators combined model. (c) Bi-LSTMs specified indicators combined model. (d) Bi-LSTMs dynamic optimized indicators combined model.
Figure 14. Prediction results of different indicators models. (a) Bi-LSTM single indicator model. (b) Bi-LSTMs all-indicators combined model. (c) Bi-LSTMs specified indicators combined model. (d) Bi-LSTMs dynamic optimized indicators combined model.
Sensors 23 02883 g014
Table 1. The top three indicators of the strongest correlation with MM256.
Table 1. The top three indicators of the strongest correlation with MM256.
DayThe Top 1 Indicator and SRCThe Top 2 Indicator and SRCThe Top 3 Indicator and SRC
day 22MM2630.79TP17110.76TC8620.73
day 23MM2620.51MM2520.45MM2110.34
day 24MM2640.81RH17120.77MM2630.71
day 25MM2630.66MM2640.59MM2110.58
day 26RH17120.71MM2640.51AMP10.5
day 27RH17120.78MM2640.66MM2630.65
day 28MM2630.71AMP10.49RH17120.48
Table 2. The evaluation of gas concentration prediction results by different indicators models.
Table 2. The evaluation of gas concentration prediction results by different indicators models.
Indicator(s)SensorR2MAERMSER2MAERMSE
Single IndicatorMM2560.7110.0980.1290.6670.1560.143
MM2630.5710.1100.151
MM2640.7180.1090.148
All IndicatorsMM2560.7590.0850.1180.6750.0950.138
MM2630.4880.1140.166
MM2640.7780.0870.131
Specified IndicatorsMM2560.7430.0860.1220.6870.1000.138
MM2630.5790.1130.150
MM2640.7390.1020.142
Optimized IndicatorsMM2560.9800.0270.0340.9650.0390.046
MM2630.9740.0270.037
MM2640.9400.0620.068
Table 3. The evaluation of gas abnormal or normal emission prediction results by different indicators models.
Table 3. The evaluation of gas abnormal or normal emission prediction results by different indicators models.
Indicator(s)SensorTPFNTNFPFARMARRFARMARR
Single IndicatorMM25671530210.3%68.2%31.5%1.2%67.2%31.6%
MM2630831700.0%100.0%0.0%
MM264168291103.3%33.3%63.3%
All IndicatorsMM256111130120.7%50.0%49.0%0.4%70.8%28.7%
MM2630831700.0%100.0%0.0%
MM26491529920.7%62.5%36.8%
Specified IndicatorsMM25671530210.3%68.2%31.5%0.4%58.8%40.7%
MM2633531610.3%62.5%37.2%
MM264131129920.7%45.8%53.5%
Optimized IndicatorsMM25617530300.0%22.7%77.3%0.0%20.1%79.9%
MM2636231700.0%25.0%75.0%
MM26421330100.0%12.5%87.5%
Disclaimer/Publisher’s Note: The statements, opinions and data contained in all publications are solely those of the individual author(s) and contributor(s) and not of MDPI and/or the editor(s). MDPI and/or the editor(s) disclaim responsibility for any injury to people or property resulting from any ideas, methods, instructions or products referred to in the content.

Share and Cite

MDPI and ACS Style

Peng, Y.; Song, D.; Qiu, L.; Wang, H.; He, X.; Liu, Q. Combined Prediction Model of Gas Concentration Based on Indicators Dynamic Optimization and Bi-LSTMs. Sensors 2023, 23, 2883. https://doi.org/10.3390/s23062883

AMA Style

Peng Y, Song D, Qiu L, Wang H, He X, Liu Q. Combined Prediction Model of Gas Concentration Based on Indicators Dynamic Optimization and Bi-LSTMs. Sensors. 2023; 23(6):2883. https://doi.org/10.3390/s23062883

Chicago/Turabian Style

Peng, Yujie, Dazhao Song, Liming Qiu, Honglei Wang, Xueqiu He, and Qiang Liu. 2023. "Combined Prediction Model of Gas Concentration Based on Indicators Dynamic Optimization and Bi-LSTMs" Sensors 23, no. 6: 2883. https://doi.org/10.3390/s23062883

APA Style

Peng, Y., Song, D., Qiu, L., Wang, H., He, X., & Liu, Q. (2023). Combined Prediction Model of Gas Concentration Based on Indicators Dynamic Optimization and Bi-LSTMs. Sensors, 23(6), 2883. https://doi.org/10.3390/s23062883

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop