# Prediction of Mine Dust Concentration Based on Grey Markov Model

Authors: Zhou Xu; Guo Liwen; Zhang Jiuling; Qin Sijia; Zhu Yi
Journal: Shock and Vibration
 (2021)

---

## Abstract

Accurate quantitative analysis and prediction of dust concentration in mines play a vital role in avoiding pneumoconiosis to a certain extent, improving industrial production efficiency, and protecting the ecological environment. The research has far-reaching significance for the prediction of dust concentration in mines in the future. Aiming at the shortcomings of the grey GM (1, 1) model in forecasting the data sequence with large random fluctuation, a grey Markov chain forecasting model is established. Firstly, considering the timeliness of monitoring data, the new dust concentration data is supplemented by using the method of cubic spline interpolation in the original data sequence. Therefore, the GM (1, 1) model is established by the method of metabolism. Then, the GM (1, 1) model is optimized by the theory of the Markov chain model. According to the relative error range generated during the prediction, the state interval is divided. Subsequently, the corresponding state probability transition matrix is constructed to obtain the grey Markov prediction model. The model was applied to the prediction of mine dust concentration and compared with the prediction results of the BP neural network model, grey prediction model, and ARIMA (1, 2, 1) model. The results showed that the prediction accuracy of the grey Markov model was significantly improved compared with other traditional prediction models. Therefore, the rationality and accuracy of this model in the prediction of mine dust concentration were verified.


---

## Body

## 1. Introduction

The continuous development of industrial technology and high technology has made the global economy develop rapidly and brought great convenience to our lives. However, in industrial production, with the improvement of coal mine tunneling and mining mechanization level, the dust production of mine working face is increasing [1]. In mine accidents, the occurrence of most accidents is related to the dust concentration in the mine. Mine dust will not only damage the mechanical life and endanger the safety and health of miners but also cause explosions when the dust reaches a certain concentration and even cause more serious underground disaster accidents [2]. Therefore, mine dust prevention work is becoming more and more important; the correct prediction of mine dust concentration is of great significance to the safe production under the mine.Domestic and foreign scholars have conducted a lot of research on the prediction of mine dust concentration, and neural network, grey prediction model, and time series prediction model are widely used prediction methods. Wang [3] analyzed the influencing factors of dust concentration in the fully mechanized working face of coal roadway and used LMBP neural network to predict the dust concentration. Zhao and Ma [4] and others used the BP neural network optimized by particle swarm optimization algorithm to predict the dust concentration. Li [5] quantitatively analyzed the dust in the construction stage and used BP neural network to predict the concentration of PM10 so as to provide targeted protection measures for practitioners and provide effective protection for the health of construction personnel. Wang [6] et al. established a grey-generalized regression neural network combination model to study the accuracy of the model in predicting pneumoconiosis. However, the neural network model needs a large number of original sample data for long-term training and optimization, and the results are satisfactory, and the prediction ability for fluctuation data is generally weak.Li [7] and Zhao [8] established a grey prediction model to predict and comprehensively evaluate the dust concentration and the operation status of petrochemical equipment. Wang [9] used the grey prediction model GM (1, 1) to predict and optimize PID optimization, which improved the accuracy of prediction data and the accuracy of shape control. Many experts and scholars have also continuously explored and improved the grey theory [10, 11]. However, the grey GM (1, 1) model is based on the exponential curve, which has poor fitting accuracy for large random fluctuation prediction objects. Through accumulation of the original data, its volatility is weakened. There is a certain deviation between the prediction results and the actual value.Wang et al. [12] used the research method of time series to construct the prediction model of mine dust. Although the traditional time series model has a good fitting effect on the data samples with time variables as independent variables, it can only achieve short-term prediction, and the error increases rapidly when predicting.Although the above model shows high fitting and accuracy in the field of traditional fitting prediction, due to the forgetting mechanism in the model, there will be some errors in the prediction of this kind of data with obvious time series characteristics. Liu et al. [13] constructed a grey prediction model based on the direct reporting data of the new coronavirus epidemic in Hubei Province from March 1 to March 21, 2020, and used the Markov model to correct the prediction results to predict the daily discharge number in Hubei Province from March 22 to March 24. It is concluded that the prediction accuracy of the grey Markov model is higher than that of the single grey prediction model. Shi [14, 15] established a grey Markov model to explore its application in the field of occupational disease prediction. Lu et al. [16], Wang and Peng [17], Weng et al. [18], and other people apply the grey Markov model to the prediction of subway passenger flow, production safety accidents, and surface subsidence. It is concluded that the grey Markov model is better than the prediction accuracy of the traditional grey GM (1, 1) model and more in line with the actual situation.At the same time, the grey Markov prediction model can analyze the data with a small amount of data and high prediction accuracy and can handle the data with large random fluctuation [19, 20], but it has not been applied in the field of mine prediction.In order to solve these problems, based on the dust concentration of a mine, this paper constructs the grey Markov model and the neural network model and compares them with the ARIMA model and the grey prediction model. The results show that the grey Markov model is superior in the prediction of dust concentration.The rest of this article is organized as follows. Section2 introduces the basic theory, including the neural network model, grey model, and Markov model. Section 3 uses ARIMA, neural network model, grey prediction model, and grey Markov model to predict and compare mine dust concentration. Section 4 solves the measurement performance evaluation index and tests the residuals of the four types of models. Finally, Section 5, draws conclusions, finds the most suitable evaluation model, and makes suggestions for future work.


## 2. Theoretical Basis of the Model

### 2.1. Theoretical Basis of the Neural Network

The learning and training process of the neural network method includes two stages: forward propagation and backward propagation. In the forward propagation process, the sample data are transmitted from the input layer to the output layer through the transfer function of the hidden layer [4]. If the output layer does not get the required output, it enters the backward propagation process and returns the error signal along the original forward propagation path. The mean square error and gradient descent method are used to modify the network connection weight and adjust the mean square error between the actual output of the network and the guided learning signal [21, 22]. This process has to be repeated until the specified error requirements or the maximum number of training are reached [23]. Figure 1 shows the concept map of neural networks.Figure 1
Neural network diagram.


### 2.2. GreyGM1,1 Prediction Model

Grey system theory is a kind of control theory proposed by Professor Deng [24]. It takes the grey system as the research object and makes scientific quantitative predictions for the prediction and control of the development of various grey systems. It is widely used in agriculture, industry, meteorology, and other fields [25, 26]. The concrete steps of the grey GM1,1 prediction model are as follows:


#### 2.2.1. Data Processing

For a given raw data sequence, theGM1,1 model of a first-order differential equation is generated by one accumulation, and the original sequence is processed as follows:(1)X1=x11,x12,…x1n,whereX1 is 1-AGO sequence of X0.


#### 2.2.2. Construction of Differential Equations and Calculation of Parameters

Accumulate the time series of samples and construct the grey differential equation as follows:(2)dx1dt+ax1=b,Z1=z11,z12,…z1n,wherea is the development coefficient of the model, reflecting the trend of x1 and x0, and u is the grey action of the model, reflecting the change relationship between data.Ifα^=a,bT is a parameter column, then(3)Y=x02x03⋯x0nB=−z121−z131⋯⋯−z1n1.


#### 2.2.3. Model Construction and Solution

The least-squares estimation parameter sequence of a grey differential equationx0k+az1k=b satisfies(4)α^=BTB−1BTY.The time response sequence ofGM1,1 grey differential equation is x^1k+1 for x10=x01; then,(5)x^1k+1=x0−bae−ak+ba.The predicted datax^1k+1 is reduced to the predicted value as follows:(6)x^0k+1=x^1k+1−x^1k.


### 2.3. Grey Markov Model

The basic Markov model is(7)xk+1=xk⋅R1,where xk represents the state of time t=k, R1 represents a one-step transition probability matrix, and xk+1 represents the state of time t=k+1. The grey Markov chain model is based on the prediction results of the grey prediction model to predict the correct stationary process of the overall trend, and then the Markov model is used to optimize and correct the prediction results. [27] The flow chart of the grey Markov model is shown in Figure 2.Figure 2
Grey Markov model process.The details of the Grey Markov model are introduced briefly as follows:Step1: the result state division of GM (1, 1) model predictionThe fitting data obtained by the model are compared with the actual sample data, and their relative values are used as the correction valueλ (k) for interval division. Each interval represents a state [28].(8)λk=x^0kx0k,Ei∈Qi1,Qi2,i=1,2,…,k,whereQi1 represents the lower limit of the relative value in the state interval and Qi2 represents the upper limit of the relative value in the state interval.Step2
: construction of the state probability matrixAccording to the interval state of the prediction results, the one-step transition probability is determined, and the state matrix of the research system is constructed [27–30]. The state transition probability matrix of the above model is as follows:(9)pk=p11kp21k⋯p1nkp21kp22k⋯p2nk⋯⋯⋯⋯pn1kpn2k⋯pnnk.Step3
: model correctionThe residuals between the measured value and the simulation value are calculated, and the average residuals of all states in each state interval are taken to correct the predicted value of the grey Markov model.


## 2.1. Theoretical Basis of the Neural Network

The learning and training process of the neural network method includes two stages: forward propagation and backward propagation. In the forward propagation process, the sample data are transmitted from the input layer to the output layer through the transfer function of the hidden layer [4]. If the output layer does not get the required output, it enters the backward propagation process and returns the error signal along the original forward propagation path. The mean square error and gradient descent method are used to modify the network connection weight and adjust the mean square error between the actual output of the network and the guided learning signal [21, 22]. This process has to be repeated until the specified error requirements or the maximum number of training are reached [23]. Figure 1 shows the concept map of neural networks.Figure 1
Neural network diagram.


## 2.2. GreyGM1,1 Prediction Model

Grey system theory is a kind of control theory proposed by Professor Deng [24]. It takes the grey system as the research object and makes scientific quantitative predictions for the prediction and control of the development of various grey systems. It is widely used in agriculture, industry, meteorology, and other fields [25, 26]. The concrete steps of the grey GM1,1 prediction model are as follows:


### 2.2.1. Data Processing

For a given raw data sequence, theGM1,1 model of a first-order differential equation is generated by one accumulation, and the original sequence is processed as follows:(1)X1=x11,x12,…x1n,whereX1 is 1-AGO sequence of X0.


### 2.2.2. Construction of Differential Equations and Calculation of Parameters

Accumulate the time series of samples and construct the grey differential equation as follows:(2)dx1dt+ax1=b,Z1=z11,z12,…z1n,wherea is the development coefficient of the model, reflecting the trend of x1 and x0, and u is the grey action of the model, reflecting the change relationship between data.Ifα^=a,bT is a parameter column, then(3)Y=x02x03⋯x0nB=−z121−z131⋯⋯−z1n1.


### 2.2.3. Model Construction and Solution

The least-squares estimation parameter sequence of a grey differential equationx0k+az1k=b satisfies(4)α^=BTB−1BTY.The time response sequence ofGM1,1 grey differential equation is x^1k+1 for x10=x01; then,(5)x^1k+1=x0−bae−ak+ba.The predicted datax^1k+1 is reduced to the predicted value as follows:(6)x^0k+1=x^1k+1−x^1k.


## 2.2.1. Data Processing

For a given raw data sequence, theGM1,1 model of a first-order differential equation is generated by one accumulation, and the original sequence is processed as follows:(1)X1=x11,x12,…x1n,whereX1 is 1-AGO sequence of X0.


## 2.2.2. Construction of Differential Equations and Calculation of Parameters

Accumulate the time series of samples and construct the grey differential equation as follows:(2)dx1dt+ax1=b,Z1=z11,z12,…z1n,wherea is the development coefficient of the model, reflecting the trend of x1 and x0, and u is the grey action of the model, reflecting the change relationship between data.Ifα^=a,bT is a parameter column, then(3)Y=x02x03⋯x0nB=−z121−z131⋯⋯−z1n1.


## 2.2.3. Model Construction and Solution

The least-squares estimation parameter sequence of a grey differential equationx0k+az1k=b satisfies(4)α^=BTB−1BTY.The time response sequence ofGM1,1 grey differential equation is x^1k+1 for x10=x01; then,(5)x^1k+1=x0−bae−ak+ba.The predicted datax^1k+1 is reduced to the predicted value as follows:(6)x^0k+1=x^1k+1−x^1k.


## 2.3. Grey Markov Model

The basic Markov model is(7)xk+1=xk⋅R1,where xk represents the state of time t=k, R1 represents a one-step transition probability matrix, and xk+1 represents the state of time t=k+1. The grey Markov chain model is based on the prediction results of the grey prediction model to predict the correct stationary process of the overall trend, and then the Markov model is used to optimize and correct the prediction results. [27] The flow chart of the grey Markov model is shown in Figure 2.Figure 2
Grey Markov model process.The details of the Grey Markov model are introduced briefly as follows:Step1: the result state division of GM (1, 1) model predictionThe fitting data obtained by the model are compared with the actual sample data, and their relative values are used as the correction valueλ (k) for interval division. Each interval represents a state [28].(8)λk=x^0kx0k,Ei∈Qi1,Qi2,i=1,2,…,k,whereQi1 represents the lower limit of the relative value in the state interval and Qi2 represents the upper limit of the relative value in the state interval.Step2
: construction of the state probability matrixAccording to the interval state of the prediction results, the one-step transition probability is determined, and the state matrix of the research system is constructed [27–30]. The state transition probability matrix of the above model is as follows:(9)pk=p11kp21k⋯p1nkp21kp22k⋯p2nk⋯⋯⋯⋯pn1kpn2k⋯pnnk.Step3
: model correctionThe residuals between the measured value and the simulation value are calculated, and the average residuals of all states in each state interval are taken to correct the predicted value of the grey Markov model.


## 3. Model Establishment and Solution

### 3.1. Prediction of Mine Dust Concentration Based on the ARIMA Model

Taking the dust concentration of 5424 working face in a mine as the original data source [12], the cubic spline interpolation method was used to interpolate the data collected between 9:31 and 10:17, and the model was established, so that the data were sufficient to support the analysis. Since the original sequence is a nonstationary sequence [29], the second-order difference method is selected to smooth the original data, and the time sequence diagram is drawn. It is found that the difference data fluctuate up and down at a certain value, so the sequence is considered to be stable. Take 2 in ARIMA model d. The sequence diagram after the difference is shown in Figure 3.Figure 3
Sequence diagram after difference.Establish the original time series model as follows:(10)1−∑i=1pαiLi1−Ldyt=α0+1+∑i=1qβiLiεt.L is a delay operator of y, satisfying(11)Liyt=yt−i.Calculated values ofp,d,q are 1, 2, and 1, respectively.Formula (10) can be reduced to(12)yt=α0+1+β1L1εt1−L21−α2L2.The fitting degree of the ARIMA model is given by SPSS. The results are shown in Table1.Table 1
ARIMA model-fitting degree.
Fitting statisticsAverageStandard errorMinMaxSmooth R0.444.0.4440.444R−0.238.-0. 238-0.238RMSE14.326.14.32614.326MAPE34.244.34.24434.244MAE9.792.9.7929.792Normalizing BIC5.886.5.8865.886Mine dust concentration is predicted by ARIMA (1, 2, 1) model. The forecast results of ARIMA are shown in Figure4.Figure 4
Forecasting results of the ARIMA model.


### 3.2. Prediction of Dust Concentration in Mines Based on the Neural Network Model

From the change trend of dust concentration reflected in Figure3, the time factor occupies a great weight in the influencing factors of dust concentration change. The relationship between the concentration change trend and the time variable is not a simple linear relationship, so the neural network model is used to predict this variable.


#### 3.2.1. Determination of Hidden Layers

The hidden layer is responsible for automatically learning the input feature information, and the multilayered hidden layer is a complex network to solve complex problems. With the increase of network layers, the accuracy and speed of calculation will also be improved [30, 31]. However, the complexity of the model also increases, and the risk of overfitting increases. In order to improve the fitting degree of the model prediction results, considering that the actual data sample size is small and it is not easy to overfit, this paper selects a ten-layer network to establish the model.


#### 3.2.2. Selection of Activation Functions

The essence of continuous learning of neural networks is to optimize its model parameters in receiving the feedback of continuous training errors. Among them, the error of the current layer is closely related to the selection of the activation function [32, 33]. Since the derivative values of the sigmoid function and tanh function are not greater than 1, the gradient will inevitably disappear in the iterative process, leading to the premature termination of training. The ReLU function can well avoid this phenomenon, so the ReLU function is selected as the activation function of the model in this paper.


#### 3.2.3. Selection of Training Algorithm

In the process of neural network training, Levenberg–Marquardt, Bayesian regularization, and quantitative conjugate gradient training algorithms are generally used for training. The sample data used in this paper is less, not easy to appear overfitting phenomenon, and the training time is short, in order to improve the fitting degree; this paper uses the training method of regularization. The results of neural network model prediction are shown in Table2.Table 2
Prediction value of BP neural network.
NumberConcentrationPre1Pre2Pre3119.610016.578919.655222.7850215.577516.594315.431223.1299315.690016.707215.610023.7446418.742517.501418.456924.7542523.530022.307823.316626.2917629.318637.127430.041928.7765……………4233.700933.270248.401744.60744331.370031.942158.313456.71074430.919731.873765.161361.74354531.370031.864763.240959.24204630.840330.815048.760049.87084727.450027.314630.790135.2170The weights and thresholds of neural network initialization are random, resulting in different results each time, so experiments are repeated. The prediction results were obtained by BP neural network, and the optimal model was selected in the three experiments.BP neural network model prediction results of the first group of data relative error is 0.1418; mean square relative error is 0.0353; the second group of prediction results relative error is 0.0874; mean square relative error is 0.0269; the third group of prediction results relative error is 0.1566; and mean square relative error is 0.0410. It can be seen that the second group has the best prediction results.


### 3.3. Prediction of Dust Concentration Based on the Grey Markov Model

#### 3.3.1. Dust Concentration Model Based on MetabolicGM1,1

Taking the dust concentration of 5424 working face in a mine as the original data source [12], the cubic spline interpolation method was used to interpolate the data collected between 9:31 and 10:17, and the metabolic model was established. According to the comparison between the predicted results of the model and the actual values, the model prediction results were finally obtained. Build the original sequence according to the mine dust concentration measured every two minutes and accumulate it once to get X1. Matrix B is calculated from formula (3):(13)B=−9.80501−17.59381⋮⋮−701.8331Y=15.5815.69⋮27.45.Calculate available by formulas (3)–(5) using MATLAB a=0.0013172 and grey action b=28.5269. Five groups of development coefficient and grey action will appear in the fourth phase of backward prediction.Bring thea,b parameters into formula (6):(14)x^1k+1=x0−bae−ak+ba=−21637.6175e−0.0013172k+21657.2275.The simulation value of mine dust concentration generated by the GM (1, 1) model is simulated by MATLAB, and the results are shown in Table3.Table 3
Prediction of mine dust concentration by the GM (1, 1) model.
NumberConcentrationGM (1, 1)Relative errorRelative valueStata119.6119.6101E2215.577528.57160.8341580.545209E1315.6928.60920.8234030.548425E1418.742528.6470.5284510.654257E1523.5328.68470.2190690.820298E1629.318628.72250.02031.0208E2737.260028.76040.22811.2955E2848.067928.79830.40091.6691E3958.820028.83620.50982.0398E41065.424828.87420.55872.2659E4………………3743.140029.91960.08251.0899E23841.504129.95900.30651.4419E33935.290029.99850.27821.3854E24034.256530.03800.14991.1764E24135.290030.07760.12311.1404E24233.700930.11730.1063351.118988E24331.3730.1570.0386681.040223E24430.919730.19670.0233831.023943E24531.3730.23650.0361331.037488E24630.840330.27640.0182851.018625E24727.4530.31630.1044190.905454E1From the GM (1, 1) model, it can be predicted that the dust concentrations in the following four periods are 30.3562 g/cm3, 30.3962 g/cm3, 30.4363 g/cm3, and 30.4764 g/cm3, respectively. The total relative error of the grey prediction model is 0.2849, and the total mean square error is 0.1431.


#### 3.3.2. Numerical Simulation of the Grey Markov Model

Division of Results forGM1,1Model Prediction. Based on the mine dust concentration index from 9:31 to 10:17, the relative value between the actual concentration and the predicted value of the grey prediction model is calculated. The upper and lower limits of the relative value are divided into intervals, and the interval distance is equalized. Among them, the distribution area of relative values is relatively concentrated in [0.52 and 2.28]. The above data are divided into four states by selecting 0.52, 0.96, 1.40, 1.84, and 2.28 as the critical values, as shown in Table 4.Table 4
State partition table.
StateInterval rangeNumberE1[0.52, 0.96]23E2[0.96, 1.4]18E3[1.4, 1.84]3E4[1.84, 2.28]3Construction of State Probability Matrix. The state divided by Table 2 can be represented as follows:(15)M1=23,M2=18,M3=3,M4=3.In the calculation of the transition probability matrixP, it is impossible to determine the next state of the last state turning, so it not included in the calculation. One-step transition frequency matrix and one-step transition probability matrix are as follows:(16)fij4×4=184005112002010012,P1=911211005181118190023013001323.P1, P2, P3, and P4 matrices are obtained by using a four-step transfer probability matrix to calculate 10:19 dust concentration data. Select the nearest four time points: 10:17, 10:16, 10:15, and 10:14; according to the distance from the time point in turn to take the number of transfer steps in the transfer matrix corresponding to the number of transfer steps in the initial state of the corresponding row vector to form a new probability matrix and the column vector summation of the probability matrix to determine the future steering state of the system, the largest sum is the state of E2. The results of mine dust concentration prediction at 10:19 are shown in Table 5.Table 5
10:19 Prediction of mine dust concentration.
TimeInitial stateTransfer stepsE1E2E3E410:17E210.27780.61110.1111010:16E220.39700.49800.06790.037010:15E230.46310.42180.06770.047310:14E240.49610.38710.06260.0541Total1.6341.9180.30930.1384The probability of state E2 is the highest in the total row in Table3, so it can be predicted that dust concentration at 10:19 is in state E2. The predicted value of the grey Markov chain is as follows:(17)x^t=0.96+1.4×30.3562÷2=35.820316.The prediction value and relative error of dust concentration under the ARIMA model and grey Markov chain model are listed in Table6.Table 6
Comparison between ARIMA model and grey Markov model.
NumberConcentrationARIMARelative errorMean square relative errorGrey MarkovRelative errorMean square relative error119.610023.13980.18000.0324215.577521.14300.35730.1276315.690015.59000.00640.000021.17080.34930.1220418.742513.57000.27600.076221.19880.13110.0172523.530023.72000.00810.000121.22670.09790.0096629.318624.75000.15580.024333.89260.15600.0243737.260034.81000.06580.004333.93730.08920.0080848.067939.88000.17030.029046.65320.02940.0009958.820048.03000.18340.033759.40260.00990.00011065.424861.64000.05780.003359.48090.09090.00831162.750053.33000.15010.022559.55930.05080.00261248.325464.36000.33180.110146.89960.02950.00091331.370049.32000.57220.327434.20650.09040.00821421.672334.43000.58870.346521.47980.00890.00011519.610034.48000.75830.575021.50810.09680.00941622.337721.98000.01600.000321.53640.03590.00131723.530021.51000.08580.007421.56490.08350.00701818.886722.93000.21410.045821.59330.14330.02051915.690021.39000.36330.132021.62180.37810.14292020.755615.87000.23540.055421.65030.04310.00192127.450020.60000.24950.062321.67870.21020.04422228.001126.18000.06500.004221.70740.22480.05052323.530021.24000.09730.009521.73590.07620.00582417.730022.37000.26170.068521.76460.22760.05182515.690022.95000.46270.214121.79330.38900.15132620.699016.26000.21450.046021.82200.05430.00292727.450021.37000.22150.049121.85080.20400.04162829.874226.25000.12130.014734.88910.16790.02822927.450023.28000.15190.023121.90840.20190.04083022.299326.14000.17220.029721.93730.01620.00033121.570027.43000.27170.073821.96620.01840.00033229.703522.15000.25430.064735.07340.18080.03273335.290025.45000.27880.077735.11960.00480.00003428.731833.53000.16700.027935.16590.22390.05013523.530031.47000.33740.113922.08230.06150.00383632.565724.44000.24950.062335.25860.08270.00683743.140032.05000.25710.066148.46980.12350.01533841.504141.59000.00210.000035.35160.14820.02203935.290035.94000.01840.000335.39820.00310.00004034.256535.05000.02320.000535.44480.03470.00124135.290036.08000.02240.000535.49160.00570.00004233.700935.48000.05280.002835.53840.05450.00304331.370033.68000.07360.005435.58530.13440.01814430.919731.70000.02520.000635.63210.15240.02324531.370031.35000.00060.000035.67910.13740.01894630.840331.39000.01780.000335.72620.15840.02514727.450029.80000.08560.007322.43410.18270.0334The average interval data are used to replace the Markov prediction interval as its predicted value. At the same time, ARIMA (1, 2, 1) model was used to predict the dust concentration of mine, and the prediction results were compared with those of the grey Markov model. The results show that the overall relative error of the ARIMA model is 0.1865, the total mean square relative error is 0.0348, the overall relative error of the grey Markov model is 0.1262, and the total mean square relative error is 0.0259. The comparison results of prediction accuracy of the three models show that the prediction effect of the GM (1, 1) model is the most unsatisfactory, the prediction effect of the ARIMA model is in the middle, and the prediction effect of the grey Markov model is the best.It can be seen from Figure5 that the predicted values of dust concentration during the period from 9:31 to 10:17 obtained by the four models are compared with the actual values during this period. It can be obviously seen that the values predicted by the grey Markov model are closer to the actual values of comprehensive evaluation. Therefore, it is concluded that the grey Markov model can well fit the development trend and volatility characteristics of the sequence.Figure 5
Comparison of measured and predicted dust concentration.


## 3.1. Prediction of Mine Dust Concentration Based on the ARIMA Model

Taking the dust concentration of 5424 working face in a mine as the original data source [12], the cubic spline interpolation method was used to interpolate the data collected between 9:31 and 10:17, and the model was established, so that the data were sufficient to support the analysis. Since the original sequence is a nonstationary sequence [29], the second-order difference method is selected to smooth the original data, and the time sequence diagram is drawn. It is found that the difference data fluctuate up and down at a certain value, so the sequence is considered to be stable. Take 2 in ARIMA model d. The sequence diagram after the difference is shown in Figure 3.Figure 3
Sequence diagram after difference.Establish the original time series model as follows:(10)1−∑i=1pαiLi1−Ldyt=α0+1+∑i=1qβiLiεt.L is a delay operator of y, satisfying(11)Liyt=yt−i.Calculated values ofp,d,q are 1, 2, and 1, respectively.Formula (10) can be reduced to(12)yt=α0+1+β1L1εt1−L21−α2L2.The fitting degree of the ARIMA model is given by SPSS. The results are shown in Table1.Table 1
ARIMA model-fitting degree.
Fitting statisticsAverageStandard errorMinMaxSmooth R0.444.0.4440.444R−0.238.-0. 238-0.238RMSE14.326.14.32614.326MAPE34.244.34.24434.244MAE9.792.9.7929.792Normalizing BIC5.886.5.8865.886Mine dust concentration is predicted by ARIMA (1, 2, 1) model. The forecast results of ARIMA are shown in Figure4.Figure 4
Forecasting results of the ARIMA model.


## 3.2. Prediction of Dust Concentration in Mines Based on the Neural Network Model

From the change trend of dust concentration reflected in Figure3, the time factor occupies a great weight in the influencing factors of dust concentration change. The relationship between the concentration change trend and the time variable is not a simple linear relationship, so the neural network model is used to predict this variable.


### 3.2.1. Determination of Hidden Layers

The hidden layer is responsible for automatically learning the input feature information, and the multilayered hidden layer is a complex network to solve complex problems. With the increase of network layers, the accuracy and speed of calculation will also be improved [30, 31]. However, the complexity of the model also increases, and the risk of overfitting increases. In order to improve the fitting degree of the model prediction results, considering that the actual data sample size is small and it is not easy to overfit, this paper selects a ten-layer network to establish the model.


### 3.2.2. Selection of Activation Functions

The essence of continuous learning of neural networks is to optimize its model parameters in receiving the feedback of continuous training errors. Among them, the error of the current layer is closely related to the selection of the activation function [32, 33]. Since the derivative values of the sigmoid function and tanh function are not greater than 1, the gradient will inevitably disappear in the iterative process, leading to the premature termination of training. The ReLU function can well avoid this phenomenon, so the ReLU function is selected as the activation function of the model in this paper.


### 3.2.3. Selection of Training Algorithm

In the process of neural network training, Levenberg–Marquardt, Bayesian regularization, and quantitative conjugate gradient training algorithms are generally used for training. The sample data used in this paper is less, not easy to appear overfitting phenomenon, and the training time is short, in order to improve the fitting degree; this paper uses the training method of regularization. The results of neural network model prediction are shown in Table2.Table 2
Prediction value of BP neural network.
NumberConcentrationPre1Pre2Pre3119.610016.578919.655222.7850215.577516.594315.431223.1299315.690016.707215.610023.7446418.742517.501418.456924.7542523.530022.307823.316626.2917629.318637.127430.041928.7765……………4233.700933.270248.401744.60744331.370031.942158.313456.71074430.919731.873765.161361.74354531.370031.864763.240959.24204630.840330.815048.760049.87084727.450027.314630.790135.2170The weights and thresholds of neural network initialization are random, resulting in different results each time, so experiments are repeated. The prediction results were obtained by BP neural network, and the optimal model was selected in the three experiments.BP neural network model prediction results of the first group of data relative error is 0.1418; mean square relative error is 0.0353; the second group of prediction results relative error is 0.0874; mean square relative error is 0.0269; the third group of prediction results relative error is 0.1566; and mean square relative error is 0.0410. It can be seen that the second group has the best prediction results.


## 3.2.1. Determination of Hidden Layers

The hidden layer is responsible for automatically learning the input feature information, and the multilayered hidden layer is a complex network to solve complex problems. With the increase of network layers, the accuracy and speed of calculation will also be improved [30, 31]. However, the complexity of the model also increases, and the risk of overfitting increases. In order to improve the fitting degree of the model prediction results, considering that the actual data sample size is small and it is not easy to overfit, this paper selects a ten-layer network to establish the model.


## 3.2.2. Selection of Activation Functions

The essence of continuous learning of neural networks is to optimize its model parameters in receiving the feedback of continuous training errors. Among them, the error of the current layer is closely related to the selection of the activation function [32, 33]. Since the derivative values of the sigmoid function and tanh function are not greater than 1, the gradient will inevitably disappear in the iterative process, leading to the premature termination of training. The ReLU function can well avoid this phenomenon, so the ReLU function is selected as the activation function of the model in this paper.


## 3.2.3. Selection of Training Algorithm

In the process of neural network training, Levenberg–Marquardt, Bayesian regularization, and quantitative conjugate gradient training algorithms are generally used for training. The sample data used in this paper is less, not easy to appear overfitting phenomenon, and the training time is short, in order to improve the fitting degree; this paper uses the training method of regularization. The results of neural network model prediction are shown in Table2.Table 2
Prediction value of BP neural network.
NumberConcentrationPre1Pre2Pre3119.610016.578919.655222.7850215.577516.594315.431223.1299315.690016.707215.610023.7446418.742517.501418.456924.7542523.530022.307823.316626.2917629.318637.127430.041928.7765……………4233.700933.270248.401744.60744331.370031.942158.313456.71074430.919731.873765.161361.74354531.370031.864763.240959.24204630.840330.815048.760049.87084727.450027.314630.790135.2170The weights and thresholds of neural network initialization are random, resulting in different results each time, so experiments are repeated. The prediction results were obtained by BP neural network, and the optimal model was selected in the three experiments.BP neural network model prediction results of the first group of data relative error is 0.1418; mean square relative error is 0.0353; the second group of prediction results relative error is 0.0874; mean square relative error is 0.0269; the third group of prediction results relative error is 0.1566; and mean square relative error is 0.0410. It can be seen that the second group has the best prediction results.


## 3.3. Prediction of Dust Concentration Based on the Grey Markov Model

### 3.3.1. Dust Concentration Model Based on MetabolicGM1,1

Taking the dust concentration of 5424 working face in a mine as the original data source [12], the cubic spline interpolation method was used to interpolate the data collected between 9:31 and 10:17, and the metabolic model was established. According to the comparison between the predicted results of the model and the actual values, the model prediction results were finally obtained. Build the original sequence according to the mine dust concentration measured every two minutes and accumulate it once to get X1. Matrix B is calculated from formula (3):(13)B=−9.80501−17.59381⋮⋮−701.8331Y=15.5815.69⋮27.45.Calculate available by formulas (3)–(5) using MATLAB a=0.0013172 and grey action b=28.5269. Five groups of development coefficient and grey action will appear in the fourth phase of backward prediction.Bring thea,b parameters into formula (6):(14)x^1k+1=x0−bae−ak+ba=−21637.6175e−0.0013172k+21657.2275.The simulation value of mine dust concentration generated by the GM (1, 1) model is simulated by MATLAB, and the results are shown in Table3.Table 3
Prediction of mine dust concentration by the GM (1, 1) model.
NumberConcentrationGM (1, 1)Relative errorRelative valueStata119.6119.6101E2215.577528.57160.8341580.545209E1315.6928.60920.8234030.548425E1418.742528.6470.5284510.654257E1523.5328.68470.2190690.820298E1629.318628.72250.02031.0208E2737.260028.76040.22811.2955E2848.067928.79830.40091.6691E3958.820028.83620.50982.0398E41065.424828.87420.55872.2659E4………………3743.140029.91960.08251.0899E23841.504129.95900.30651.4419E33935.290029.99850.27821.3854E24034.256530.03800.14991.1764E24135.290030.07760.12311.1404E24233.700930.11730.1063351.118988E24331.3730.1570.0386681.040223E24430.919730.19670.0233831.023943E24531.3730.23650.0361331.037488E24630.840330.27640.0182851.018625E24727.4530.31630.1044190.905454E1From the GM (1, 1) model, it can be predicted that the dust concentrations in the following four periods are 30.3562 g/cm3, 30.3962 g/cm3, 30.4363 g/cm3, and 30.4764 g/cm3, respectively. The total relative error of the grey prediction model is 0.2849, and the total mean square error is 0.1431.


### 3.3.2. Numerical Simulation of the Grey Markov Model

Division of Results forGM1,1Model Prediction. Based on the mine dust concentration index from 9:31 to 10:17, the relative value between the actual concentration and the predicted value of the grey prediction model is calculated. The upper and lower limits of the relative value are divided into intervals, and the interval distance is equalized. Among them, the distribution area of relative values is relatively concentrated in [0.52 and 2.28]. The above data are divided into four states by selecting 0.52, 0.96, 1.40, 1.84, and 2.28 as the critical values, as shown in Table 4.Table 4
State partition table.
StateInterval rangeNumberE1[0.52, 0.96]23E2[0.96, 1.4]18E3[1.4, 1.84]3E4[1.84, 2.28]3Construction of State Probability Matrix. The state divided by Table 2 can be represented as follows:(15)M1=23,M2=18,M3=3,M4=3.In the calculation of the transition probability matrixP, it is impossible to determine the next state of the last state turning, so it not included in the calculation. One-step transition frequency matrix and one-step transition probability matrix are as follows:(16)fij4×4=184005112002010012,P1=911211005181118190023013001323.P1, P2, P3, and P4 matrices are obtained by using a four-step transfer probability matrix to calculate 10:19 dust concentration data. Select the nearest four time points: 10:17, 10:16, 10:15, and 10:14; according to the distance from the time point in turn to take the number of transfer steps in the transfer matrix corresponding to the number of transfer steps in the initial state of the corresponding row vector to form a new probability matrix and the column vector summation of the probability matrix to determine the future steering state of the system, the largest sum is the state of E2. The results of mine dust concentration prediction at 10:19 are shown in Table 5.Table 5
10:19 Prediction of mine dust concentration.
TimeInitial stateTransfer stepsE1E2E3E410:17E210.27780.61110.1111010:16E220.39700.49800.06790.037010:15E230.46310.42180.06770.047310:14E240.49610.38710.06260.0541Total1.6341.9180.30930.1384The probability of state E2 is the highest in the total row in Table3, so it can be predicted that dust concentration at 10:19 is in state E2. The predicted value of the grey Markov chain is as follows:(17)x^t=0.96+1.4×30.3562÷2=35.820316.The prediction value and relative error of dust concentration under the ARIMA model and grey Markov chain model are listed in Table6.Table 6
Comparison between ARIMA model and grey Markov model.
NumberConcentrationARIMARelative errorMean square relative errorGrey MarkovRelative errorMean square relative error119.610023.13980.18000.0324215.577521.14300.35730.1276315.690015.59000.00640.000021.17080.34930.1220418.742513.57000.27600.076221.19880.13110.0172523.530023.72000.00810.000121.22670.09790.0096629.318624.75000.15580.024333.89260.15600.0243737.260034.81000.06580.004333.93730.08920.0080848.067939.88000.17030.029046.65320.02940.0009958.820048.03000.18340.033759.40260.00990.00011065.424861.64000.05780.003359.48090.09090.00831162.750053.33000.15010.022559.55930.05080.00261248.325464.36000.33180.110146.89960.02950.00091331.370049.32000.57220.327434.20650.09040.00821421.672334.43000.58870.346521.47980.00890.00011519.610034.48000.75830.575021.50810.09680.00941622.337721.98000.01600.000321.53640.03590.00131723.530021.51000.08580.007421.56490.08350.00701818.886722.93000.21410.045821.59330.14330.02051915.690021.39000.36330.132021.62180.37810.14292020.755615.87000.23540.055421.65030.04310.00192127.450020.60000.24950.062321.67870.21020.04422228.001126.18000.06500.004221.70740.22480.05052323.530021.24000.09730.009521.73590.07620.00582417.730022.37000.26170.068521.76460.22760.05182515.690022.95000.46270.214121.79330.38900.15132620.699016.26000.21450.046021.82200.05430.00292727.450021.37000.22150.049121.85080.20400.04162829.874226.25000.12130.014734.88910.16790.02822927.450023.28000.15190.023121.90840.20190.04083022.299326.14000.17220.029721.93730.01620.00033121.570027.43000.27170.073821.96620.01840.00033229.703522.15000.25430.064735.07340.18080.03273335.290025.45000.27880.077735.11960.00480.00003428.731833.53000.16700.027935.16590.22390.05013523.530031.47000.33740.113922.08230.06150.00383632.565724.44000.24950.062335.25860.08270.00683743.140032.05000.25710.066148.46980.12350.01533841.504141.59000.00210.000035.35160.14820.02203935.290035.94000.01840.000335.39820.00310.00004034.256535.05000.02320.000535.44480.03470.00124135.290036.08000.02240.000535.49160.00570.00004233.700935.48000.05280.002835.53840.05450.00304331.370033.68000.07360.005435.58530.13440.01814430.919731.70000.02520.000635.63210.15240.02324531.370031.35000.00060.000035.67910.13740.01894630.840331.39000.01780.000335.72620.15840.02514727.450029.80000.08560.007322.43410.18270.0334The average interval data are used to replace the Markov prediction interval as its predicted value. At the same time, ARIMA (1, 2, 1) model was used to predict the dust concentration of mine, and the prediction results were compared with those of the grey Markov model. The results show that the overall relative error of the ARIMA model is 0.1865, the total mean square relative error is 0.0348, the overall relative error of the grey Markov model is 0.1262, and the total mean square relative error is 0.0259. The comparison results of prediction accuracy of the three models show that the prediction effect of the GM (1, 1) model is the most unsatisfactory, the prediction effect of the ARIMA model is in the middle, and the prediction effect of the grey Markov model is the best.It can be seen from Figure5 that the predicted values of dust concentration during the period from 9:31 to 10:17 obtained by the four models are compared with the actual values during this period. It can be obviously seen that the values predicted by the grey Markov model are closer to the actual values of comprehensive evaluation. Therefore, it is concluded that the grey Markov model can well fit the development trend and volatility characteristics of the sequence.Figure 5
Comparison of measured and predicted dust concentration.


## 3.3.1. Dust Concentration Model Based on MetabolicGM1,1

Taking the dust concentration of 5424 working face in a mine as the original data source [12], the cubic spline interpolation method was used to interpolate the data collected between 9:31 and 10:17, and the metabolic model was established. According to the comparison between the predicted results of the model and the actual values, the model prediction results were finally obtained. Build the original sequence according to the mine dust concentration measured every two minutes and accumulate it once to get X1. Matrix B is calculated from formula (3):(13)B=−9.80501−17.59381⋮⋮−701.8331Y=15.5815.69⋮27.45.Calculate available by formulas (3)–(5) using MATLAB a=0.0013172 and grey action b=28.5269. Five groups of development coefficient and grey action will appear in the fourth phase of backward prediction.Bring thea,b parameters into formula (6):(14)x^1k+1=x0−bae−ak+ba=−21637.6175e−0.0013172k+21657.2275.The simulation value of mine dust concentration generated by the GM (1, 1) model is simulated by MATLAB, and the results are shown in Table3.Table 3
Prediction of mine dust concentration by the GM (1, 1) model.
NumberConcentrationGM (1, 1)Relative errorRelative valueStata119.6119.6101E2215.577528.57160.8341580.545209E1315.6928.60920.8234030.548425E1418.742528.6470.5284510.654257E1523.5328.68470.2190690.820298E1629.318628.72250.02031.0208E2737.260028.76040.22811.2955E2848.067928.79830.40091.6691E3958.820028.83620.50982.0398E41065.424828.87420.55872.2659E4………………3743.140029.91960.08251.0899E23841.504129.95900.30651.4419E33935.290029.99850.27821.3854E24034.256530.03800.14991.1764E24135.290030.07760.12311.1404E24233.700930.11730.1063351.118988E24331.3730.1570.0386681.040223E24430.919730.19670.0233831.023943E24531.3730.23650.0361331.037488E24630.840330.27640.0182851.018625E24727.4530.31630.1044190.905454E1From the GM (1, 1) model, it can be predicted that the dust concentrations in the following four periods are 30.3562 g/cm3, 30.3962 g/cm3, 30.4363 g/cm3, and 30.4764 g/cm3, respectively. The total relative error of the grey prediction model is 0.2849, and the total mean square error is 0.1431.


## 3.3.2. Numerical Simulation of the Grey Markov Model

Division of Results forGM1,1Model Prediction. Based on the mine dust concentration index from 9:31 to 10:17, the relative value between the actual concentration and the predicted value of the grey prediction model is calculated. The upper and lower limits of the relative value are divided into intervals, and the interval distance is equalized. Among them, the distribution area of relative values is relatively concentrated in [0.52 and 2.28]. The above data are divided into four states by selecting 0.52, 0.96, 1.40, 1.84, and 2.28 as the critical values, as shown in Table 4.Table 4
State partition table.
StateInterval rangeNumberE1[0.52, 0.96]23E2[0.96, 1.4]18E3[1.4, 1.84]3E4[1.84, 2.28]3Construction of State Probability Matrix. The state divided by Table 2 can be represented as follows:(15)M1=23,M2=18,M3=3,M4=3.In the calculation of the transition probability matrixP, it is impossible to determine the next state of the last state turning, so it not included in the calculation. One-step transition frequency matrix and one-step transition probability matrix are as follows:(16)fij4×4=184005112002010012,P1=911211005181118190023013001323.P1, P2, P3, and P4 matrices are obtained by using a four-step transfer probability matrix to calculate 10:19 dust concentration data. Select the nearest four time points: 10:17, 10:16, 10:15, and 10:14; according to the distance from the time point in turn to take the number of transfer steps in the transfer matrix corresponding to the number of transfer steps in the initial state of the corresponding row vector to form a new probability matrix and the column vector summation of the probability matrix to determine the future steering state of the system, the largest sum is the state of E2. The results of mine dust concentration prediction at 10:19 are shown in Table 5.Table 5
10:19 Prediction of mine dust concentration.
TimeInitial stateTransfer stepsE1E2E3E410:17E210.27780.61110.1111010:16E220.39700.49800.06790.037010:15E230.46310.42180.06770.047310:14E240.49610.38710.06260.0541Total1.6341.9180.30930.1384The probability of state E2 is the highest in the total row in Table3, so it can be predicted that dust concentration at 10:19 is in state E2. The predicted value of the grey Markov chain is as follows:(17)x^t=0.96+1.4×30.3562÷2=35.820316.The prediction value and relative error of dust concentration under the ARIMA model and grey Markov chain model are listed in Table6.Table 6
Comparison between ARIMA model and grey Markov model.
NumberConcentrationARIMARelative errorMean square relative errorGrey MarkovRelative errorMean square relative error119.610023.13980.18000.0324215.577521.14300.35730.1276315.690015.59000.00640.000021.17080.34930.1220418.742513.57000.27600.076221.19880.13110.0172523.530023.72000.00810.000121.22670.09790.0096629.318624.75000.15580.024333.89260.15600.0243737.260034.81000.06580.004333.93730.08920.0080848.067939.88000.17030.029046.65320.02940.0009958.820048.03000.18340.033759.40260.00990.00011065.424861.64000.05780.003359.48090.09090.00831162.750053.33000.15010.022559.55930.05080.00261248.325464.36000.33180.110146.89960.02950.00091331.370049.32000.57220.327434.20650.09040.00821421.672334.43000.58870.346521.47980.00890.00011519.610034.48000.75830.575021.50810.09680.00941622.337721.98000.01600.000321.53640.03590.00131723.530021.51000.08580.007421.56490.08350.00701818.886722.93000.21410.045821.59330.14330.02051915.690021.39000.36330.132021.62180.37810.14292020.755615.87000.23540.055421.65030.04310.00192127.450020.60000.24950.062321.67870.21020.04422228.001126.18000.06500.004221.70740.22480.05052323.530021.24000.09730.009521.73590.07620.00582417.730022.37000.26170.068521.76460.22760.05182515.690022.95000.46270.214121.79330.38900.15132620.699016.26000.21450.046021.82200.05430.00292727.450021.37000.22150.049121.85080.20400.04162829.874226.25000.12130.014734.88910.16790.02822927.450023.28000.15190.023121.90840.20190.04083022.299326.14000.17220.029721.93730.01620.00033121.570027.43000.27170.073821.96620.01840.00033229.703522.15000.25430.064735.07340.18080.03273335.290025.45000.27880.077735.11960.00480.00003428.731833.53000.16700.027935.16590.22390.05013523.530031.47000.33740.113922.08230.06150.00383632.565724.44000.24950.062335.25860.08270.00683743.140032.05000.25710.066148.46980.12350.01533841.504141.59000.00210.000035.35160.14820.02203935.290035.94000.01840.000335.39820.00310.00004034.256535.05000.02320.000535.44480.03470.00124135.290036.08000.02240.000535.49160.00570.00004233.700935.48000.05280.002835.53840.05450.00304331.370033.68000.07360.005435.58530.13440.01814430.919731.70000.02520.000635.63210.15240.02324531.370031.35000.00060.000035.67910.13740.01894630.840331.39000.01780.000335.72620.15840.02514727.450029.80000.08560.007322.43410.18270.0334The average interval data are used to replace the Markov prediction interval as its predicted value. At the same time, ARIMA (1, 2, 1) model was used to predict the dust concentration of mine, and the prediction results were compared with those of the grey Markov model. The results show that the overall relative error of the ARIMA model is 0.1865, the total mean square relative error is 0.0348, the overall relative error of the grey Markov model is 0.1262, and the total mean square relative error is 0.0259. The comparison results of prediction accuracy of the three models show that the prediction effect of the GM (1, 1) model is the most unsatisfactory, the prediction effect of the ARIMA model is in the middle, and the prediction effect of the grey Markov model is the best.It can be seen from Figure5 that the predicted values of dust concentration during the period from 9:31 to 10:17 obtained by the four models are compared with the actual values during this period. It can be obviously seen that the values predicted by the grey Markov model are closer to the actual values of comprehensive evaluation. Therefore, it is concluded that the grey Markov model can well fit the development trend and volatility characteristics of the sequence.Figure 5
Comparison of measured and predicted dust concentration.


## 4. Residual Test

### 4.1. Evaluation Indexes of Forecasting Performance

To compare the forecast obtained using the proposed model with those obtained using other models, the mean absolute percentage error (MAPE), the root mean squared error (RMSE), and the mean absolute error (MAE) are used to evaluate forecasting accuracy.(18)MAPE=1N∑i=1Nyi−fiyix100%,RMSE=∑i=1Nyi−fi2N,MAE=∑i=1Nyi−fi2N,where N is the total number of forecasts, yi is the actual value at forecasting point i, and fi is the forecasted value at forecasting point i. A comparison of prediction accuracy of the above models is shown in Table 7.Table 7
Comparison of prediction accuracy of the above models.
ModelMAEMAPERMSEGM (1, 1)8.120.2811.58Grey Markov3.260.133.87ARIMA5.190.196.82Neural network2.350.094.31The comparison of the prediction accuracy of the above four models shows that the grey Markov model and the neural network model have better prediction accuracy, and their prediction accuracy is much higher than that of the other two models.


### 4.2. Residual Verification

The posterior error ratioC and small error probability P test were used to test the pros and cons of the prediction model and the accuracy level. The accuracy of the results is calculated by the above two prediction models. The accuracy classification table is shown in Table 8:Table 8
Classification table of prediction accuracy.
Probability of small error (P)Mean variance ratio (C)Prediction accuracy>0.95<0.35Excellent>0.80<0.50Qualification>0.70<0.65Reluctance≤0.70≥0.65NonconformityThe grey Markov model considers the dynamic transfer influence components; the dynamic transfer influence components reflect the influence degree of random factors and volatility factors, which has the advantage of predicting large random volatility problems. The variance ratio ofGM1,1 model is C=0.9928, and that of small probability error is P=0.5532. The variance ratio of Grey Markov model is C=0.3256, and that of small probability error is P=1. The accuracy of four prediction models is compared as shown in Table 9.Table 9
Comparison of prediction accuracy between four prediction models.
Forecasting modelProbability of small errorMean variance ratioPrediction accuracy(P)(C)ARIMA0.72340.5879ReluctanceNeural network0.91490.3699QualificationGM (1, 1)0.55320.9928NonconformityGrey Markov model10.3256ExcellentIt can be seen from Table9 that the prediction accuracy of the grey Markov model is higher than that of GM1,1 model. The prediction of mine dust concentration by grey Markov model is more reliable and accurate and can accurately describe the state of mine dust concentration.


## 4.1. Evaluation Indexes of Forecasting Performance

To compare the forecast obtained using the proposed model with those obtained using other models, the mean absolute percentage error (MAPE), the root mean squared error (RMSE), and the mean absolute error (MAE) are used to evaluate forecasting accuracy.(18)MAPE=1N∑i=1Nyi−fiyix100%,RMSE=∑i=1Nyi−fi2N,MAE=∑i=1Nyi−fi2N,where N is the total number of forecasts, yi is the actual value at forecasting point i, and fi is the forecasted value at forecasting point i. A comparison of prediction accuracy of the above models is shown in Table 7.Table 7
Comparison of prediction accuracy of the above models.
ModelMAEMAPERMSEGM (1, 1)8.120.2811.58Grey Markov3.260.133.87ARIMA5.190.196.82Neural network2.350.094.31The comparison of the prediction accuracy of the above four models shows that the grey Markov model and the neural network model have better prediction accuracy, and their prediction accuracy is much higher than that of the other two models.


## 4.2. Residual Verification

The posterior error ratioC and small error probability P test were used to test the pros and cons of the prediction model and the accuracy level. The accuracy of the results is calculated by the above two prediction models. The accuracy classification table is shown in Table 8:Table 8
Classification table of prediction accuracy.
Probability of small error (P)Mean variance ratio (C)Prediction accuracy>0.95<0.35Excellent>0.80<0.50Qualification>0.70<0.65Reluctance≤0.70≥0.65NonconformityThe grey Markov model considers the dynamic transfer influence components; the dynamic transfer influence components reflect the influence degree of random factors and volatility factors, which has the advantage of predicting large random volatility problems. The variance ratio ofGM1,1 model is C=0.9928, and that of small probability error is P=0.5532. The variance ratio of Grey Markov model is C=0.3256, and that of small probability error is P=1. The accuracy of four prediction models is compared as shown in Table 9.Table 9
Comparison of prediction accuracy between four prediction models.
Forecasting modelProbability of small errorMean variance ratioPrediction accuracy(P)(C)ARIMA0.72340.5879ReluctanceNeural network0.91490.3699QualificationGM (1, 1)0.55320.9928NonconformityGrey Markov model10.3256ExcellentIt can be seen from Table9 that the prediction accuracy of the grey Markov model is higher than that of GM1,1 model. The prediction of mine dust concentration by grey Markov model is more reliable and accurate and can accurately describe the state of mine dust concentration.


## 5. Conclusion

The exploration of dust concentration changing with time in mines has been one of the research hotspots in the field of mine safety. Dust is extremely harmful to human health, industrial production, and living environment. It is of great significance to accurately predict the dust concentration by building a model. The following conclusions can be drawn based on the study of the dust concentration in mines:(1)
The Markov chain theory is used to optimize the GM (1, 1) model. According to the range of relative error generated by GM (1, 1) model prediction, the state interval is reasonably divided, and the corresponding state transition probability matrix is determined. The GM (1, 1) model based on Markov optimization is established.(2)
The model combines the advantages of the grey prediction model and the Markov model to improve the fitting accuracy of prediction objects with large random fluctuation. The experimental results show that the grey Markov model has better prediction accuracy than grey prediction, ARIMA, and neural network models.(3)
In the analysis of the prediction results of mine dust concentration, it can be concluded that the grey Markov model is suitable for short-term prediction with less sample data. The experimental results show that this model is more accurate than the traditional model. It has great application space in mine dust concentration prediction.



---

*Source: 5859249-2021-12-23.xml*