Distribution Law and Prediction Model of Dust Concentration under Airflow Adjustment in Fully Mechanized Heading Face

The current situation of extensive ventilation management for the fully mechanized heading face cannot dynamically control air velocity and make reasonable dust migration distribution, resulting in serious disaster of dust and environmental pollution in the roadway. In this paper, the fluid mechanics, finite element numerical simulation, and underground measurement were combined to analyze the distribution of dust concentration under the variable airflow state at the duct outlet to obtain the massive correlation data of duct outlet parameters and dust concentration. For the pedestrian breathing-height in the backflow side and driver location, the double-objective BP prediction model for dust concentration under airflow adjustment was established, and the performance parameters and prediction accuracy of the BP prediction model were analyzed by using the relative error and fitting analysis. In Ningtiaota coal mine, located in Northern Shaanxi province of China, the self-developed control device is installed on the duct outlet with underground test and application verification to refine the model. The results indicated the dynamic control of airflow optimizes dust concentration distribution. The dust concentration at the pedestrian breathing-height in the backflow side and driver location was significantly decreased after the installation of adjustment device. Dust concentration at the pedestrian breathing-height and driver location was decreased by 31% and 34%, respectively, compared with the results before adjustment, which achieved the safe, environment-friendly, and energy-saving ventilation and the dust removal function in fully mechanized heading face.


Introduction
With the continuous breakthroughs of cross section size, driving speed, and the length of coal roadway, the traditionally extensive control and mode of ventilation become inadequate to meet the requirements with higher security level, efficient and environmental protection [1][2][3]. In the process of excavation, the airflow cannot be dynamically regulated according to actual ventilation demand, resulting in the unreasonable airflow distribution; besides, dust concentration reaches 1500 mg/m 3 , which is far exceeding the standard required by The Coal Mine Safety Regulations of 4 mg/m 3 [4,5], causing the dust explosion accidents and pollution risks. If the air supply demand is configured according to the initial stage of excavation, the wind energy loss will increase as the roadway extended under the specific ventilation scheme and air volume, and the air velocity often fails to meet the requirements of the regulations, which cause the problem of serious dust accumulation. If it is configured according to the maximum demand, the air velocity will be too large in the initial stage of excavation and exceed the specified range, which causes discomfort feelings of workers, secondary blowing dust, and serious pollution at breathing-height area [6][7][8][9]. At present, many researchers have shown that the duct position, diameter, and distance from the duct outlet to heading face have a great influence on the distribution of dust concentration [10][11][12][13][14][15][16]. In the research, the authors found that the change of airflow state at the duct outlet has a significant effect on controlling dust migration distribution and creatively put forward a new method to optimize the distribution of airflow, gas, and dust fields by dynamic adjustment of airflow [ [17][18][19][20]. The new control device for airflow adjustment at the duct outlet is developed [21], and the installation position is shown in Figure 1. According to the theory of fluid mechanics, numerical simulation and underground measurement method, analyzing the influence of airflow adjustment on dust concentration distribution, and the correlation mechanism among airflow adjustment parameters, air velocity and dust concentration are established. Based on the large amount of sample data obtained, the BP neural network algorithm is used for data mining [22,23]; then, the dust safety hazard and pollution comprehensive prediction neural network model and optimal adjustment scheme are established. In the coal-bearing coal mine of the Shaanxi Coal Group, the self-developed control device is used for the underground installation test and application. Hence, the distribution of dust concentration under airflow adjustment and the accuracy of neural network prediction model are verified, which provides a new theoretical basis for the safe, environment-friendly, and efficient ventilation in fully mechanized heading face.

Theoretical Analysis
(1) Airflow Field. The airflow motion in the roadway complies with the law of the physical conservation. It is assumed that there is no heat source in the roadway, and the thermal radiation from the wall is not considered. The airflow motion meets the law of mass conservation and momentum conservation, which are embodied as the commonly used continuity equation and the N-S equation in the fluid mechanics. The airflow is treated as the incompressible continuous phase; therefore the continuity equation can be expressed as For a control volume, the fluid element motion obeys Newton's second law; therefore the N-S equation of the incompressible viscous fluid can be described by Among them, Considering that the airflow is incompressible ideal fluid, ∇ • = 0, = 0. Then, (2) can be simplified as follows: where ( / 3 ) is the air density, ( ) is the time, ( / ) is the air velocity, ( / ) is the velocity of the surrounding fluid, ( ) is the hydrostatic pressure, and ( • ) is the dynamic viscosity of air.
(2) Dust Field. In dust field, the dust emission source is set at the heading face, and its injection type is dynamic according to the specific cutting mode and sequence of the roadheader. The dust particles are considered as discrete phase, and air flows continuously. The Euler-Lagrange method is used to simulate the distribution of airflow and dust migration under the dynamic change of the dust emission source. The force among dust particles is neglected, but dust particles themselves should be loaded with gravity, buoyancy, drag force, etc. Therefore, gas-solid two-phase flow mathematical model is used to analyze the dust migration law. The unit mass dust motion equation in the Lagrangian coordinate system can be expressed as Among them, ( − ) is the unit mass drag force, and has the following form, which is shown as where ( ) is the drag force, ( / ) is the duct particle velocity, ( / ) is the air velocity, ( / 3 ) is the dust particle density, ( / 3 ) is the air density, ( ) is the Saffman force, ( • ) is the dynamic viscosity of air, is the drag coefficient, is the particle Reynolds number, and ( ) is the dust particle diameter.

Finite Element Analysis Method and Verification.
According to the actual condition of the heading face in Ningtiaota coal mine, the finite element model under the dynamic change of dust emission source is established. The geometric model of the roadway is a cuboid of 6.1m × 4.0m × 40m, the diameter of the duct is 1m, and the distance from the central axis of the duct to the top and side walls is 0.95 m and 0.7 m. In terms of the roadheader, it is simplified into two parts: the machine part (cuboid) and the cutting part (cylinder). The origin of the coordinate is located on the near-duct side of the boundary between the head section and the bottom plate. The X-axis represents the width of the roadway, the Yaxis represents the height of the roadway, and the negative direction of the Z-axis is the digging direction. The main area of the pedestrian movement in the roadway is represented by two lines, where the red line (X=3.05m, Y=2m, and Z=0∼ 40m) indicates the variation of driver location along the roadway and the blue line (X=5m, Y=1.5m, and Z=0∼40m) indicates the variation of pedestrians location in the backflow side. The finite element model of heading face is shown in Figure 2.
The geometric model is imported into Fluent to set boundary conditions. The end of the outlet is set as speed inlet, the tail section of the roadway is set as outlet with standard atmospheric pressure, and the total air volume is 360 m 3 /min, assuming that the walls surface are not slipped, and the airflow is incompressible. The dust emission source is set at the heading face and its injection type is dynamic. The other numerical simulation parameters of airflow and dust field are set as shown in Tables 1 and 2.
In order to verify the accuracy and feasibility of the numerical simulation scheme, the numerical simulation  results of airflow field and dust field at different distances from the duct outlet to heading face are compared with the data measured in the roadway. Air velocity and dust concentration at the distance of 8m are selected as the verification objects, and the measuring points at the pedestrian breathing-height in the backflow side (X=5m, Y=1.5m) are arranged in the sections of 5m, 7.5m, 10m, 12.5m, 15m, 17.5m, 20m, 22.5m, 25m, and 30m from the heading face. The simulated data and measured data are shown in Tables 3 and  4, and the simulated data is compared with the measured data (see Figure 3). It can be seen from Table 3 and Figure 3(a) that the simulated values are similar to the underground measured   values. Although there is a certain error because the actual ventilation is affected by complicated factors such as the roadheader, the error is basically controlled at about 10%, which proves that the numerical simulation of airflow field is feasible. From Table 4 and Figure 3(b), it can be seen that the numerical simulated results are consistent with the measured results, and the error is basically controlled at about 8%, which verifies the feasibility and accuracy of the numerical simulation in the dust field.

Influence of Airflow Adjustment on Air Velocity and
Dust Concentration

Simulation Schemes of Airflow Field and Dust Field.
Aiming at the problems of airflow and dust field under the original ventilation condition, the influence of airflow adjustment parameters on the distribution of airflow and dust migration is observed. There are dust pollution and safety hazards at the driver location and backflow side; we found that, at distance of 5∼10m from the duct outlet to heading face, the outlet parameters on caliber of 0.7∼1.2m, horizontal right deviation of 0 ∘ ∼25 ∘ , and vertical upward deflection of 0 ∘ ∼6 ∘ , there is a great influence on the distribution of airflow and dust migration in the roadway. In order to analyze influence of airflow adjustment on air velocity and dust concentration, the simulation schemes are designed according to the orthogonal experiment (see Table 5).

Influence of Airflow Adjustment on Air
Velocity. The airflow field under the different adjustment schemes is simulated and calculated. Here list the nearest distance of 5m and the farthest distance of 10m from the duct outlet to heading face, the air velocity along the pedestrian breathing-height in the backflow side (X=5m, Y=1.5m) and the driver's breathingheight (X=3.05m, Y=2m) are shown in Figures 4 and 5.
It can be seen from Figure 4 that, at the nearest distance of 5m, the air velocity in some areas has exceeded the most comfortable range of 0.5∼1.0m/s, causing discomfort feelings  Air velocity (m/s)  for workers. Increasing the caliber can reduce the air velocity and expand the range of airflow radiation. However, if the caliber is too large, the air velocity will be too low to effectively discharge the dust. Therefore, the caliber should be within a reasonable range. From Figure 5, it can be seen that, at the farthest distance of 10m, the air velocity is too low, and the air velocity in some areas is even lower than the minimum limit of 0.25m/s specified by The Coal Mine Safety Regulations [24]. By reducing the caliber, the effective range is increased, and the air velocity in the backflow side is increased. But the boundary layer of the jet is too narrow when the caliber is too small, causing dust accumulated. The vertical upward deflection of the duct outlet can effectively reduce dust concentration at the air duct side and the backflow side, and the horizontal right deviation can reduce dust concentration at the backflow side. However, there is a safety hazard of dust accumulation in the air duct side when the horizontal right angle is too large.

Influence of Airflow Adjustment on Dust Concentration.
The dust field under different adjustment schemes is simulated and calculated. Here list the nearest distance of 5m and the farthest distance of 10m from the duct outlet to heading face, dust concentration distribution at pedestrian breathingheight area (Y=1.5m) and driver's breathing-height area (Y=2m) are shown in Figures 6 and 7, and dust concentration along the pedestrian breathing-height (X=5m, Y=1.5m) and   the driver's breathing-height (X=3.05m, Y=2m) are shown in Figures 8 and 9. It can be seen from Figures 6 and 8 that, at the nearest distance of 5m, the secondary blowing dust easily occurs due to the high air velocity at the outlet. At the caliber of 1.1m, horizontal right deviation of 10 ∘ , and vertical upward deflection of 6 ∘ , dust concentration is obviously decreased. Significantly, dust concentration at the pedestrian breathingheight in the backflow side is decreased by 30% compared with the original field, and dust concentration at the driver location is decreased by 35% compared with the original one. From Figures 7 and 9, it can be seen that, at the farthest distance of 10m, the outlet parameters on caliber of 0.9m, horizontal right deviation of 15 ∘ , and vertical upward deflection of 0 ∘ , there is a great influence on the overall dust concentration. Respectively, dust concentration at the pedestrian breathing-height is reduced by 38%, and dust concentration at the driver location is reduced by 27%, compared with the results in the original dust field.
In order to meet the more stringent requirements in safety, high efficiency, and environmental protection, it is necessary to analyze dust concentration distribution under the change of the airflow state, so as to dynamically control the airflow for optimizing the distribution of airflow and dust migration in the heading face, thereby improving the underground working environment. Besides, predicting dust concentration in the roadway is also significant, and thus grasping dust concentration at the key position. Through the analysis of the distribution law of dust concentration under airflow adjustment, the factors that have a great influence on dust concentration are determined, specifically the distance from the duct outlet to heading face, the duct outlet caliber, the horizontal right deviation, and the vertical upward deflection, which lays a foundation for the establishment of the dust concentration prediction model under the control of airflow.

Prediction Model of Dust Concentration under
Airflow Adjustment

Design of BP Prediction Model
Structure. The BP neural network generally consists of three-layer structure, namely, the input layer, the hidden layer, and the output layer. Each neuron is connected to all neurons in the next layer, and there is no connection between neurons in the same layer [25]. The factors affecting dust concentration of the heading face are extremely complicated, including the coal seam and the mining technology. Besides, there is a complex nonlinear relationship between dust concentration and various influencing factors. The analytic hierarchy process (AHP) method is used to analyze the key factors affecting dust concentration in the heading face. The hierarchy general ranking is shown in Table 6.
From Table 6, it can be found that the distance from the duct outlet to heading face, the duct outlet caliber, the horizontal right deviation, the vertical upward deflection, the air velocity at the duct outlet, and dust production play a main role in the dust concentration distribution of roadway; then, we take them as input layer parameters of BP neural network prediction model.  the greatest impact on the workers which is the breathingheight area, the dust safety hazard and pollution risks comprehensive prediction neural network model under airflow adjustment are established. Thus, the number of nodes in the output layer is 2, including both dust concentration at the driver location and dust concentration at the pedestrian breathing-height in the backflow side. The determination of the number of nodes in hidden layer is a difficult problem, and the method affects the prediction results and the complexity of the neural network. Equation (6) is commonly used to calculate the neurons number of the hidden layers.
where is the number of nodes in hidden layer, is the number of nodes in input layer, is the number of nodes in output layer, and is the constant between 1 and 10.
Using empirical equation and MATLAB software program, a network with variable number of hidden layer nodes is designed, and the optimal number of nodes is determined by error comparison. The results of MATLAB operation are shown in Table 7.
It can be seen from Table 7 that, under the premise that the number of nodes in hidden layer is 5, the network error is 1.81e-2, which is the minimum. Thus, the number of nodes in the hidden layer is determined as 5, and the topology of the network model is 6-5-2. The specific structure is shown in Figure 10      concentration at the driver location is the value of the point (X=3.05m, Y=2m, and Z=7.5m). According to the air velocity between 0.25 and 4 m/s specified by The Coal Mine Safety Regulations [24], sample data preprocessing is performed to eliminate unreasonable adjustment schemes. Therefore, the input and output data are obtained (see Tables 8 and 9). The maximum value is 586.11, and the minimum value is 0.0076.
In order to facilitate the analysis and the network convergence, the data is normalized by (7). The normalized data Mathematical Problems in Engineering 9  are used as training and testing data for the BP prediction model.
where is the input or output data, is the maximum value, and is the minimum value.

Determining Functions of Neural Network.
Combined with the sample data in Tables 8 and 9, the transfer, training, and learning functions are used to train the double-objective BP prediction model. According to the mean square error and iteration times, the most suitable transfer function, training function, and learning function are selected. The comparison results of different functions are shown in Table 10.
Under the premise that the mean square error is the minimum, the function with the smallest iteration times is selected. For the double-objective prediction model of dust concentration under airflow adjustment, the transfer function is logsig, the training function is trainlm, and the learning function is learngd.

Double-Objective Prediction Model of Dust Concentration under Airflow
Adjustment. The number of input nodes at double-objective prediction model is 6, including the distance from the duct outlet to heading face, the duct outlet caliber, the horizontal right deviation, the vertical upward deflection, the air velocity at the duct outlet, and dust production. The number of nodes in hidden layer is 5, and the number of output nodes is 2, involving dust concentration at the driver location and dust concentration at the pedestrian breathing-height in the backflow side. The transfer function is logsig, the training function is trainlm, and the learning function is learngd. According to the designed double-objective prediction model, the top 74 data is trained by using the neural network toolbox. The max training times for the network is set as 1000 and the learning rate is 0.5. The training results are shown in Figure 11.
It can be seen from Figure 11(a) that the network training has reached the accuracy requirement after 13 iterations, and the number of verification steps is 6. From Figure 11(b), it can be seen that the network model needs 13 training steps to achieve convergence, and the mean square error is 0.015532. The prediction model of dust concentration under airflow adjustment is trained by using BP neural network, the mean square error is minimal, and the network tends to be stable. From Figure 11(c), it can be seen that the network training is ended when the gradient approaches 0.01; the value of mu determines whether the learning is done according to the Newton method or the gradient method; val fail is the maximum number of verification failures. When the error is tested at 6 times, the error does not decrease or rise, indicating that the training set error is no longer decreased, and the network stops training.
The weights between nodes in different layers are shown in Tables 11 and 12, and thresholds of hidden layer nodes are shown in Table 13.
In addition, the thresholds of the output layer nodes are 0.82971 and 1.2618.
The latter 10 data in Tables 8 and 9 is used as test sample to examine the effect of the dust concentration prediction model under the airflow adjustment. The comparison between test values and predicted values is shown in Table 14. It can be seen that the maximum relative error is 9.11%, and the minimum value is 0.74%. The prediction model has high precision and meets the requirements; therefore, it can accurately predict dust concentration at the pedestrian breathing-height and driver location under airflow adjustment.

Analysis on Practical Application.
Taking Ningtiaota coal mine in Northern Shaanxi province of China as the research object, the analysis of dust concentration prediction and optimal adjustment scheme are carried out. At the nearest distance of 5m and the farthest distance of 10m from the duct outlet to heading face, the adjustment schemes are established, and the BP prediction model is used to predict dust concentration under different adjustment schemes. By analyzing and comparing the minimum value, the optimal adjustment scheme of the best dust removal effect is obtained.
At the nearest distance of 5m from the duct outlet to heading face, the effective range of the jet region is large, resulting in high air velocity at the end of the jet, and air velocity at the outlet is reduced by increasing the caliber of duct outlet. Therefore, the duct outlet caliber is 1.1m or 1.2m, the horizontal right deviation is 5∼25 ∘ , and the vertical upward deflection is 2∼6 ∘ . The airflow adjustment schemes at the nearest distance of 5m are shown in Table 15.
According to the practical investigation in Ningtiaota coal mine, the dust production is 0.0096 kg/s, air velocity at the outlet is 5.94 m/s when the caliber is 1.1m, and air velocity at the outlet is 5.27 m/s when the caliber is 1.2m. The established BP prediction model is used to predict dust concentration of 30 schemes; thus, the dust concentration at the pedestrian breathing-height and driver location is obtained (see Table 16).
According to scheme 9 from Table 16, the optimal adjustment scheme at the nearest distance of 5m is obtained. Dust concentration at the driver location is 325.36 mg/m 3 after the adjustment of scheme 9, and dust concentration at the pedestrian breathing-height is 304.32 mg/m 3 , which all reach the minimum relatively. And the corresponding distance is 5m from the duct outlet to heading face, the caliber is 1.1m, the horizontal right deviation is 15 ∘ , and the vertical upward deflection is 6 ∘ . The dust concentration is compared with that before adjustment (see Table 17).
At the farthest distance of 10m from the duct outlet to heading face, the established neural network model is used to predict dust concentration under different schemes, and the optimal adjustment scheme is obtained. After adjustment, dust concentration at the driver location is 376.09 mg/m 3 , which is decreased by 35%, and dust concentration at the pedestrian breathing-height is 323.49 mg/m 3 , which is decreased by 33%. The corresponding distance is 10m from the duct outlet to heading face, the caliber is 0.8m, the horizontal right deviation is 15 ∘ , and the vertical upward deflection is 2 ∘ .

Underground Test Verification.
We have developed the duct outlet control device and carried out the underground actual tests in Ningtiaota coal mine last year. The main hazard of dust in the heading face is reflected in the impact on the health of underground workers. Considering the location that has the greatest impact on workers, which is breathingheight area, dust concentration at the driver location and pedestrian breathing-height in the backflow side is selected as the test objects. The underground measuring points are shown in Figure 12.
(1) Underground Test at the Nearest Distance of 5m. According to the arrangement of underground measuring points in Figure 12, the original dust concentration measured data is obtained at the nearest distance of 5m (see Table 18). Taking the optimal adjustment scheme obtained by the prediction model as the underground test scheme, the measured data is shown in Table 19. Then, dust concentration before and after adjustment is compared (see Figure 13).
From Tables 18 and 19 and Figure 13, it can be seen that the dust removal effect is obvious after adjustment. Dust concentration at the driver location is decreased by 34%, and dust concentration at the pedestrian breathing-height in the backflow side is decreased by 29% compared with that before adjustment. (2) Underground Test at the Farthest Distance of 10m. According to the underground measuring points in Figure 12, the original dust concentration measured data are obtained at the farthest distance of 10m (see Table 20). Taking the optimal adjustment scheme obtained by the prediction model as the underground test scheme, the measured data are shown in Table 21, and dust concentration before and after adjustment is compared (see Figure 14).

Mathematical Problems in Engineering
It can be seen from Tables 20 and 21 and Figure 14 that dust concentration at the driver location is decreased by 33%, and dust concentration at the pedestrian breathing-height in the backflow side is decreased by 31%.

Conclusions
Based on the airflow and dust field numerical simulation analysis and the large amount of sample data under the change of airflow adjustment parameters, the doubleobjective BP prediction model of dust concentration is established for the safety hazard in the backflow side and the pollution risks at driver location. The distribution law and    Figure 12: Arrangement of underground measuring points. prediction model of dust concentration under airflow adjustment are verified in S1212 of Ningtiaota coal mine through the self-developed control device. The specific conclusions are as follows.   (1) The problems of airflow and dust field under the original ventilation mode and the influence of airflow adjustment parameters on the distribution of airflow and dust migration are analyzed. We come to the conclusion of the distribution law of air velocity and dust concentration under various adjustment schemes. At the nearest distance of 5m, the caliber of duct outlet is 1.1m or 1.2m, the horizontal right deviation is 5∼15 ∘ , the vertical upward deflection is 2∼6 ∘ , and the airflow and dust field are distributed reasonably. At the farthest distance of 10m, the caliber of duct outlet is 0.8m or 0.9m, the horizontal right deviation is 5∼15 ∘ , and the vertical upward deflection is 0∼4 ∘ , the minimum air velocity is higher than the minimum limit of 0.25 m/s specified by The Coal Mine Safety Regulations, and the dust field is better distributed.
(2) The input layer parameters of the double-objective prediction model are determined as the distance from the duct outlet to heading face, the duct outlet caliber, the horizontal right deviation, the vertical upward deflection, the air velocity at the duct outlet, and dust production. The optimal number of hidden layer nodes is 5. And the output layer includes dust concentration at the driver location and the pedestrian breathing-height in the backflow side. The      neural network prediction model is established and tested by function with the smallest mean square error. The maximum predicted error is 9.11%, which proves that it can accurately predict dust concentration at the driver location and the pedestrian breathing-height. And the practical application of the double-objective prediction model is carried out.
(3) Taking Ningtiaota coal mine in Northern Shaanxi province of China as the research object, the airflow adjustment schemes are designed, and the established prediction model is used to predict dust concentration under various schemes; then, the optimal adjustment scheme is obtained. At the nearest distance of 5m, the caliber of 1.1m, the horizontal right deviation of 15 ∘ , and the vertical upward deflection of 6 ∘ , dust concentration at the driver location is decreased from 541.1 mg/m 3 to 357.4 mg/m 3 , which is decreased by 34%, and dust concentration at the pedestrian breathing-height in the backflow side decreased from 486.2 mg/m 3 to 347.3 mg/m 3 , which is decreased by 29%. At the farthest distance of 10m, the caliber of 0.8m, the horizontal right deviation of 15 ∘ , and the vertical upward deflection of 2 ∘ , dust concentration at the driver location is decreased to 397.8 mg/m 3 , and the dust removal effect reached 33%. And dust concentration at the pedestrian breathing-height in the backflow side is decreased to 346.2 mg/m 3 ; the dust removal effect reached 31% compared with that before adjustment.

Data Availability
The data used to support the findings of this study are available from the corresponding author upon request.