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

Authors: Xiaoyan Gong; Congcong Jia; Kang Sun; Jian Cui; Kefan Lei; Yuxuan Xue; He Xue
Journal: Mathematical Problems in Engineering
 (2019)

---

## Abstract

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.


---

## Body

## 1. 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–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/m3, which is far exceeding the standard required by The Coal Mine Safety Regulations of 4 mg/m3 [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–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–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–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.Figure 1
Schematic diagram of airflow structure and installation position of the adjustment device. (a) Schematic diagram of airflow structure. (b) Installation position of the adjustment device.
(a) (b)


## 2. Analysis Methods of Airflow and Dust Field

### 2.1. 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 (1)∂ρ∂t+∂∂xiρui=0For 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(2)ρdudt=ρf-∇p+μ∇2uAmong them,∇2=∂2/∂x12+∂2/∂x22+∂2/∂x32. Considering that the airflow is incompressible ideal fluid, ∇∙u=0, μ=0. Then, (2) can be simplified as follows: (3)∂u∂t+u∙∇=f-1ρ∇p+v∇2uwhere ρkg/m3 is the air density, ts is the time, um/s is the air velocity, uim/s is the velocity of the surrounding fluid, pPa is the hydrostatic pressure, and μPa∙s 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(4)dupdt=FDu-up+gρp-ρρp+FsAmong them,FDu-up is the unit mass drag force, and FD has the following form, which is shown as (5)FD=18μρpdp2CDRep24where FDN is the drag force, upm/s is the duct particle velocity, um/s is the air velocity, ρpkg/m3 is the dust particle density, ρkg/m3 is the air density, FsN is the Saffman force, μPa∙s is the dynamic viscosity of air, CD is the drag coefficient, Rep is the particle Reynolds number, and dpm is the dust particle diameter.


### 2.2. 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 Y-axis 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 Figure2.Figure 2
The finite element model of fully mechanized heading face. (a) Geometric model. (b) Meshing model.
(a) (b)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 m3/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.Table 1
Setting of airflow field parameters.
Condition Define Inlet Boundary Type Velocity-inlet Viscous Model Realizable k-ε Inlet Velocity Magnitude (m/s) 7.64 Turbulence Intensity 3.08% Hydraulic Diameter (m) 1 Outlet Boundary Type Outflow Gravitational Acceleration (m/s2) -9.8Table 2
Setting of dust field parameters.
Condition Define Turbulent Dispersion Discrete Random Walk Model Diameter Distribution Rosin-Rammler Min. Diameter (m) 1e-06 Max. Diameter (m) 1e-04 Mid. Diameter (m) 1.2e-05 Number of Diameters 10 Spread Parameter 1.77 Total Flow Rate (kg/s) 0.0076 Time Scale Constant 0.15In 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 Tables3 and 4, and the simulated data is compared with the measured data (see Figure 3).Table 3
Air velocity comparison of measured data and simulated data.
Distance (m) 5 7.5 10 12.5 15 17.5 20 22.5 25 30 Measured data (m/s) 1.52 — 0.18 0.21 0.20 0.17 0.12 — — 0.14 Simulated data (m/s) 1.41 0.03 0.17 0.19 0.19 0.16 0.11 0.06 0.07 0.13 Error rate (%) 7.2 — 5.6 9.5 5.0 5.9 8.3 — — 7.1Table 4
Dust concentration comparison of measured data and simulated data.
Distance (m) 5 7.5 10 12.5 15 17.5 20 22.5 25 30 Measured data (mg/m3) 249 230 246 225 327 258 271 272 304 232 Simulated data (mg/m3) 238 218 233 208 314 244 257 254 285 217 Error rate (%) 4.4 5.2 5.3 7.6 4.0 5.4 5.2 6.6 6.3 6.5 Note: “—” means that the air velocity is too small to be tested.Figure 3
Comparison of measured data and simulated data. (a) Air velocity along the pedestrian breathing-height. (b) Dust concentration along the pedestrian breathing-height.
(a) (b)It can be seen from Table3 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.


## 2.1. 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 (1)∂ρ∂t+∂∂xiρui=0For 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(2)ρdudt=ρf-∇p+μ∇2uAmong them,∇2=∂2/∂x12+∂2/∂x22+∂2/∂x32. Considering that the airflow is incompressible ideal fluid, ∇∙u=0, μ=0. Then, (2) can be simplified as follows: (3)∂u∂t+u∙∇=f-1ρ∇p+v∇2uwhere ρkg/m3 is the air density, ts is the time, um/s is the air velocity, uim/s is the velocity of the surrounding fluid, pPa is the hydrostatic pressure, and μPa∙s 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(4)dupdt=FDu-up+gρp-ρρp+FsAmong them,FDu-up is the unit mass drag force, and FD has the following form, which is shown as (5)FD=18μρpdp2CDRep24where FDN is the drag force, upm/s is the duct particle velocity, um/s is the air velocity, ρpkg/m3 is the dust particle density, ρkg/m3 is the air density, FsN is the Saffman force, μPa∙s is the dynamic viscosity of air, CD is the drag coefficient, Rep is the particle Reynolds number, and dpm is the dust particle diameter.


## 2.2. 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 Y-axis 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 Figure2.Figure 2
The finite element model of fully mechanized heading face. (a) Geometric model. (b) Meshing model.
(a) (b)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 m3/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.Table 1
Setting of airflow field parameters.
Condition Define Inlet Boundary Type Velocity-inlet Viscous Model Realizable k-ε Inlet Velocity Magnitude (m/s) 7.64 Turbulence Intensity 3.08% Hydraulic Diameter (m) 1 Outlet Boundary Type Outflow Gravitational Acceleration (m/s2) -9.8Table 2
Setting of dust field parameters.
Condition Define Turbulent Dispersion Discrete Random Walk Model Diameter Distribution Rosin-Rammler Min. Diameter (m) 1e-06 Max. Diameter (m) 1e-04 Mid. Diameter (m) 1.2e-05 Number of Diameters 10 Spread Parameter 1.77 Total Flow Rate (kg/s) 0.0076 Time Scale Constant 0.15In 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 Tables3 and 4, and the simulated data is compared with the measured data (see Figure 3).Table 3
Air velocity comparison of measured data and simulated data.
Distance (m) 5 7.5 10 12.5 15 17.5 20 22.5 25 30 Measured data (m/s) 1.52 — 0.18 0.21 0.20 0.17 0.12 — — 0.14 Simulated data (m/s) 1.41 0.03 0.17 0.19 0.19 0.16 0.11 0.06 0.07 0.13 Error rate (%) 7.2 — 5.6 9.5 5.0 5.9 8.3 — — 7.1Table 4
Dust concentration comparison of measured data and simulated data.
Distance (m) 5 7.5 10 12.5 15 17.5 20 22.5 25 30 Measured data (mg/m3) 249 230 246 225 327 258 271 272 304 232 Simulated data (mg/m3) 238 218 233 208 314 244 257 254 285 217 Error rate (%) 4.4 5.2 5.3 7.6 4.0 5.4 5.2 6.6 6.3 6.5 Note: “—” means that the air velocity is too small to be tested.Figure 3
Comparison of measured data and simulated data. (a) Air velocity along the pedestrian breathing-height. (b) Dust concentration along the pedestrian breathing-height.
(a) (b)It can be seen from Table3 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.


## 3. Results

### 3.1. Influence of Airflow Adjustment on Air Velocity and Dust Concentration

#### 3.1.1. 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).Table 5
The simulation schemes under airflow adjustment.
Scheme L (m) d (m) a (°) b (°) 1 5 1.1 5 4 2 5 1.1 10 6 3 5 1.1 15 2 4 5 1.2 5 6 5 5 1.2 10 2 6 5 1.2 15 4 7 6 1.0 5 4 ⋮ ⋮ ⋮ ⋮ ⋮ 34 10 0.9 5 2 35 10 0.9 10 4 36 10 0.9 15 0 where L (m) represents the distance from the duct outlet to heading face, d (m) represents the duct outlet caliber, a (°) represents the horizontal right deviation of duct outlet, b (°) represents the vertical upward deflection of duct outlet.


#### 3.1.2. 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 breathing-height (X=3.05m, Y=2m) are shown in Figures4 and 5.Figure 4
Comparison of air velocity at the nearest distance of 5m. (a) Air velocity along the pedestrian breathing-height. (b) Air velocity along the driver’s breathing-height.
(a) (b)Figure 5
Comparison of air velocity at the farthest distance of 10m. (a) Air velocity along the pedestrian breathing-height. (b) Air velocity along the driver’s breathing-height.
(a) (b)It can be seen from Figure4 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 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.


#### 3.1.3. 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 breathing-height area (Y=1.5m) and driver’s breathing-height area (Y=2m) are shown in Figures6 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.Figure 6
Contours of dust concentration under different adjustment schemes at the nearest distance of 5m.Figure 7
Contours of dust concentration under different adjustment schemes at the farthest distance of 10m.Figure 8
Comparison of dust concentration at the nearest distance of 5m. (a) Dust concentration along the pedestrian breathing-height. (b) Dust concentration along the driver’s breathing-height.
(a) (b)Figure 9
Comparison of dust concentration at the farthest distance of 10m. (a) Dust concentration along the pedestrian breathing-height. (b) Dust concentration along the driver’s breathing-height.
(a) (b)It can be seen from Figures6 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 breathing-height 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.


### 3.2. Prediction Model of Dust Concentration under Airflow Adjustment

#### 3.2.1. 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 Table6.Table 6
The hierarchy general ranking.
Rule hierarchy Weights Factor hierarchy Weights Total sequencing weight Adjustment parameters 0.6250 Distance from the duct outlet to heading face 0.4673 0.292 Caliber of the duct outlet 0.2772 0.173 Horizontal right deviation 0.1601 0.100 Vertical upward deflection 0.0954 0.060 Original gas and dust parameters 0.3285 Gas emission quantity 0.5715 0.186 Dust production 0.2857 0.094 Particle size 0.1428 0.047 Ventilation parameters 0.1365 Air velocity at the duct outlet 0.75 0.102 Ventilation volume 0.25 0.034From Table6, 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 main hazard caused by dust is impact on the health of the underground workers. Considering the location that has the greatest impact on the workers which is the breathing-height 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. (6)m=n+l+awhere m is the number of nodes in hidden layer, n is the number of nodes in input layer, l is the number of nodes in output layer, and a 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 Table7.Table 7
Error under the different number of nodes in hidden layer.
The number of nodes Error Iteration times The number of nodes Error Iteration times 4 5.76e-2 7 9 6.20e-2 4 5 1.81e-2 16 10 1.01e-1 6 6 1.14e-1 8 11 8.85e-2 18 7 6.60e-2 6 12 5.58e-2 14 8 6.13e-2 18 13 3.11e-2 12It can be seen from Table7 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.Figure 10
The network model structure of dust concentration prediction.


#### 3.2.2. Data Sample and Preprocessing

Considering the simulation results of the airflow and dust field under the comprehensive changes of influencing factors, dust concentration at the pedestrian breathing-height in the backflow side is the average at three points of Z=7m, Z=10m, and Z=15m on the blue line (X=5m, Y=1.5m, and Z=0~40m), and dust 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.Table 8
The input layer data.
Number Distance (m) Caliber (m) Horizontal right deviation (°) Vertical upward deflection (°) Air velocity (m/s) Dust production (kg/s) 1 5 1.1 5 2 5.43 0.0076 2 5 1.1 5 6 5.68 0.0076 3 5 1.1 15 2 5.17 0.0076 4 5 1.1 15 6 5.31 0.0076 5 5 1.2 5 6 4.82 0.0076 6 5 1.2 15 2 5.27 0.0076 7 5 1.2 25 2 5.19 0.0076 8 5 1.2 5 2 4.98 0.0076 9 6 1 5 6 5.61 0.0076 ⋮ ⋮ ⋮ ⋮ ⋮ ⋮ ⋮ 76 9 0.9 25 6 6.51 0.0116 77 10 0.7 5 2 6.81 0.0116 78 10 0.7 15 6 6.94 0.0116 79 10 0.7 25 6 6.87 0.0116 80 10 0.7 15 2 7.04 0.0116 81 10 0.8 5 2 6.51 0.0116 82 10 0.8 15 6 6.43 0.0116 83 10 0.8 25 6 6.72 0.0116 84 10 0.8 15 2 6.86 0.0116Table 9
The output layer data.
Number Dust concentration at the Dust concentration at the pedestrian driver location (mg/m3) breathing-height (mg/m3) 1 293.40 280.45 2 307.82 294.34 3 274.41 286.35 4 220.64 270.35 5 326.64 243.83 6 283.89 304.80 7 331.24 318.07 8 322.52 237.90 9 319.91 221.07 ⋮ ⋮ ⋮ 76 508.42 532.20 77 462.64 528.20 78 461.91 491.39 79 403.29 487.09 80 447.94 507.06 81 445.01 516.14 82 454.15 478.11 83 499.62 483.24 84 452.98 507.23In order to facilitate the analysis and the network convergence, the data is normalized by (7). The normalized data are used as training and testing data for the BP prediction model. (7)xi∧=xi-xminxmax-xminwhere xi is the input or output data, xmax is the maximum value, and xmin is the minimum value.


#### 3.2.3. Determining Functions of Neural Network

Combined with the sample data in Tables8 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.Table 10
The comparison results of different functions.
Function type Function name Iteration times Mean square error Transfer function logsig 9 2.52e-2 tansig 11 4.65e-2 purelin 9 6.68e-1 Training function traingdm 1000 3.17e-1 trainglm 11 5.57e-2 traingdx 137 8.52e-1 Learning function learngdm 139 7.41e-2 learngd 25 3.46e-2Under 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.


#### 3.2.4. 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 Figure11.Figure 11
Training results of dust concentration prediction model. (a) Training interface. (b) Mean squared error. (c) Training state.
(a) (b) (c)It can be seen from Figure11(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 Tables11 and 12, and thresholds of hidden layer nodes are shown in Table 13.Table 11
The weights from input layer nodes to hidden layer nodes.
Input layer nodes X1 X2 X3 X4 X5 X6 Hidden layer nodes 2.4314 -0.0622 -1.2612 1.0935 -2.5342 1.5057 1.3383 -1.0115 2.0125 -0.7830 -0.7608 0.6870 -0.5912 -1.0331 -0.2196 0.1168 -1.4287 -1.9597 2.8285 -4.8427 0.0719 -0.5322 -0.0306 -1.8395 0.6135 -1.1564 0.4116 0.2652 1.8586 -0.3303Table 12
The weights from hidden layer nodes to output layer nodes.
Hidden layer nodes 1 2 3 4 5 Output layer nodes 2.2022 -0.4141 -2.1345 -0.7098 0.9605 -0.3397 -0.3385 -1.9866 -0.4598 -1.0131Table 13
Thresholds of hidden layer nodes.
Hidden layer nodes 1 2 3 4 5 Thresholds -4.2854 -3.2851 0.4871 0.3751 -4.8327In addition, the thresholds of the output layer nodes are 0.82971 and 1.2618.The latter 10 data in Tables8 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.Table 14
The comparison between test value and predicted value.
Test value at the driver location (mg/m3) Test value at the pedestrian breathing-height (mg/m3) Predicted value at the driver location (mg/m3) Predicted value at the pedestrian breathing-height (mg/m3) Absolute error Relative error (%) 455.35 476.21 480.74 508.49 25.39 32.28 5.58 6.78 508.42 532.20 512.19 483.70 3.77 48.50 0.74 9.11 462.64 528.20 483.55 533.60 20.91 5.40 4.52 1.02 461.91 491.39 449.80 515.04 12.11 23.65 2.62 4.81 403.29 487.09 415.54 481.65 12.25 5.44 3.04 1.12 447.94 507.06 452.53 523.94 4.59 16.88 1.02 3.33 445.01 516.14 468.57 521.58 23.56 5.44 5.29 1.05 454.15 478.11 467.89 493.32 13.74 15.21 3.03 3.18 499.62 483.24 475.09 490.14 24.53 6.90 4.91 1.43 452.98 507.23 465.85 545.50 12.87 38.27 2.84 7.54


## 3.1. Influence of Airflow Adjustment on Air Velocity and Dust Concentration

### 3.1.1. 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).Table 5
The simulation schemes under airflow adjustment.
Scheme L (m) d (m) a (°) b (°) 1 5 1.1 5 4 2 5 1.1 10 6 3 5 1.1 15 2 4 5 1.2 5 6 5 5 1.2 10 2 6 5 1.2 15 4 7 6 1.0 5 4 ⋮ ⋮ ⋮ ⋮ ⋮ 34 10 0.9 5 2 35 10 0.9 10 4 36 10 0.9 15 0 where L (m) represents the distance from the duct outlet to heading face, d (m) represents the duct outlet caliber, a (°) represents the horizontal right deviation of duct outlet, b (°) represents the vertical upward deflection of duct outlet.


### 3.1.2. 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 breathing-height (X=3.05m, Y=2m) are shown in Figures4 and 5.Figure 4
Comparison of air velocity at the nearest distance of 5m. (a) Air velocity along the pedestrian breathing-height. (b) Air velocity along the driver’s breathing-height.
(a) (b)Figure 5
Comparison of air velocity at the farthest distance of 10m. (a) Air velocity along the pedestrian breathing-height. (b) Air velocity along the driver’s breathing-height.
(a) (b)It can be seen from Figure4 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 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.


### 3.1.3. 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 breathing-height area (Y=1.5m) and driver’s breathing-height area (Y=2m) are shown in Figures6 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.Figure 6
Contours of dust concentration under different adjustment schemes at the nearest distance of 5m.Figure 7
Contours of dust concentration under different adjustment schemes at the farthest distance of 10m.Figure 8
Comparison of dust concentration at the nearest distance of 5m. (a) Dust concentration along the pedestrian breathing-height. (b) Dust concentration along the driver’s breathing-height.
(a) (b)Figure 9
Comparison of dust concentration at the farthest distance of 10m. (a) Dust concentration along the pedestrian breathing-height. (b) Dust concentration along the driver’s breathing-height.
(a) (b)It can be seen from Figures6 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 breathing-height 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.


## 3.1.1. 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).Table 5
The simulation schemes under airflow adjustment.
Scheme L (m) d (m) a (°) b (°) 1 5 1.1 5 4 2 5 1.1 10 6 3 5 1.1 15 2 4 5 1.2 5 6 5 5 1.2 10 2 6 5 1.2 15 4 7 6 1.0 5 4 ⋮ ⋮ ⋮ ⋮ ⋮ 34 10 0.9 5 2 35 10 0.9 10 4 36 10 0.9 15 0 where L (m) represents the distance from the duct outlet to heading face, d (m) represents the duct outlet caliber, a (°) represents the horizontal right deviation of duct outlet, b (°) represents the vertical upward deflection of duct outlet.


## 3.1.2. 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 breathing-height (X=3.05m, Y=2m) are shown in Figures4 and 5.Figure 4
Comparison of air velocity at the nearest distance of 5m. (a) Air velocity along the pedestrian breathing-height. (b) Air velocity along the driver’s breathing-height.
(a) (b)Figure 5
Comparison of air velocity at the farthest distance of 10m. (a) Air velocity along the pedestrian breathing-height. (b) Air velocity along the driver’s breathing-height.
(a) (b)It can be seen from Figure4 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 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.


## 3.1.3. 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 breathing-height area (Y=1.5m) and driver’s breathing-height area (Y=2m) are shown in Figures6 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.Figure 6
Contours of dust concentration under different adjustment schemes at the nearest distance of 5m.Figure 7
Contours of dust concentration under different adjustment schemes at the farthest distance of 10m.Figure 8
Comparison of dust concentration at the nearest distance of 5m. (a) Dust concentration along the pedestrian breathing-height. (b) Dust concentration along the driver’s breathing-height.
(a) (b)Figure 9
Comparison of dust concentration at the farthest distance of 10m. (a) Dust concentration along the pedestrian breathing-height. (b) Dust concentration along the driver’s breathing-height.
(a) (b)It can be seen from Figures6 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 breathing-height 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.


## 3.2. Prediction Model of Dust Concentration under Airflow Adjustment

### 3.2.1. 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 Table6.Table 6
The hierarchy general ranking.
Rule hierarchy Weights Factor hierarchy Weights Total sequencing weight Adjustment parameters 0.6250 Distance from the duct outlet to heading face 0.4673 0.292 Caliber of the duct outlet 0.2772 0.173 Horizontal right deviation 0.1601 0.100 Vertical upward deflection 0.0954 0.060 Original gas and dust parameters 0.3285 Gas emission quantity 0.5715 0.186 Dust production 0.2857 0.094 Particle size 0.1428 0.047 Ventilation parameters 0.1365 Air velocity at the duct outlet 0.75 0.102 Ventilation volume 0.25 0.034From Table6, 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 main hazard caused by dust is impact on the health of the underground workers. Considering the location that has the greatest impact on the workers which is the breathing-height 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. (6)m=n+l+awhere m is the number of nodes in hidden layer, n is the number of nodes in input layer, l is the number of nodes in output layer, and a 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 Table7.Table 7
Error under the different number of nodes in hidden layer.
The number of nodes Error Iteration times The number of nodes Error Iteration times 4 5.76e-2 7 9 6.20e-2 4 5 1.81e-2 16 10 1.01e-1 6 6 1.14e-1 8 11 8.85e-2 18 7 6.60e-2 6 12 5.58e-2 14 8 6.13e-2 18 13 3.11e-2 12It can be seen from Table7 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.Figure 10
The network model structure of dust concentration prediction.


### 3.2.2. Data Sample and Preprocessing

Considering the simulation results of the airflow and dust field under the comprehensive changes of influencing factors, dust concentration at the pedestrian breathing-height in the backflow side is the average at three points of Z=7m, Z=10m, and Z=15m on the blue line (X=5m, Y=1.5m, and Z=0~40m), and dust 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.Table 8
The input layer data.
Number Distance (m) Caliber (m) Horizontal right deviation (°) Vertical upward deflection (°) Air velocity (m/s) Dust production (kg/s) 1 5 1.1 5 2 5.43 0.0076 2 5 1.1 5 6 5.68 0.0076 3 5 1.1 15 2 5.17 0.0076 4 5 1.1 15 6 5.31 0.0076 5 5 1.2 5 6 4.82 0.0076 6 5 1.2 15 2 5.27 0.0076 7 5 1.2 25 2 5.19 0.0076 8 5 1.2 5 2 4.98 0.0076 9 6 1 5 6 5.61 0.0076 ⋮ ⋮ ⋮ ⋮ ⋮ ⋮ ⋮ 76 9 0.9 25 6 6.51 0.0116 77 10 0.7 5 2 6.81 0.0116 78 10 0.7 15 6 6.94 0.0116 79 10 0.7 25 6 6.87 0.0116 80 10 0.7 15 2 7.04 0.0116 81 10 0.8 5 2 6.51 0.0116 82 10 0.8 15 6 6.43 0.0116 83 10 0.8 25 6 6.72 0.0116 84 10 0.8 15 2 6.86 0.0116Table 9
The output layer data.
Number Dust concentration at the Dust concentration at the pedestrian driver location (mg/m3) breathing-height (mg/m3) 1 293.40 280.45 2 307.82 294.34 3 274.41 286.35 4 220.64 270.35 5 326.64 243.83 6 283.89 304.80 7 331.24 318.07 8 322.52 237.90 9 319.91 221.07 ⋮ ⋮ ⋮ 76 508.42 532.20 77 462.64 528.20 78 461.91 491.39 79 403.29 487.09 80 447.94 507.06 81 445.01 516.14 82 454.15 478.11 83 499.62 483.24 84 452.98 507.23In order to facilitate the analysis and the network convergence, the data is normalized by (7). The normalized data are used as training and testing data for the BP prediction model. (7)xi∧=xi-xminxmax-xminwhere xi is the input or output data, xmax is the maximum value, and xmin is the minimum value.


### 3.2.3. Determining Functions of Neural Network

Combined with the sample data in Tables8 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.Table 10
The comparison results of different functions.
Function type Function name Iteration times Mean square error Transfer function logsig 9 2.52e-2 tansig 11 4.65e-2 purelin 9 6.68e-1 Training function traingdm 1000 3.17e-1 trainglm 11 5.57e-2 traingdx 137 8.52e-1 Learning function learngdm 139 7.41e-2 learngd 25 3.46e-2Under 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.


### 3.2.4. 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 Figure11.Figure 11
Training results of dust concentration prediction model. (a) Training interface. (b) Mean squared error. (c) Training state.
(a) (b) (c)It can be seen from Figure11(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 Tables11 and 12, and thresholds of hidden layer nodes are shown in Table 13.Table 11
The weights from input layer nodes to hidden layer nodes.
Input layer nodes X1 X2 X3 X4 X5 X6 Hidden layer nodes 2.4314 -0.0622 -1.2612 1.0935 -2.5342 1.5057 1.3383 -1.0115 2.0125 -0.7830 -0.7608 0.6870 -0.5912 -1.0331 -0.2196 0.1168 -1.4287 -1.9597 2.8285 -4.8427 0.0719 -0.5322 -0.0306 -1.8395 0.6135 -1.1564 0.4116 0.2652 1.8586 -0.3303Table 12
The weights from hidden layer nodes to output layer nodes.
Hidden layer nodes 1 2 3 4 5 Output layer nodes 2.2022 -0.4141 -2.1345 -0.7098 0.9605 -0.3397 -0.3385 -1.9866 -0.4598 -1.0131Table 13
Thresholds of hidden layer nodes.
Hidden layer nodes 1 2 3 4 5 Thresholds -4.2854 -3.2851 0.4871 0.3751 -4.8327In addition, the thresholds of the output layer nodes are 0.82971 and 1.2618.The latter 10 data in Tables8 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.Table 14
The comparison between test value and predicted value.
Test value at the driver location (mg/m3) Test value at the pedestrian breathing-height (mg/m3) Predicted value at the driver location (mg/m3) Predicted value at the pedestrian breathing-height (mg/m3) Absolute error Relative error (%) 455.35 476.21 480.74 508.49 25.39 32.28 5.58 6.78 508.42 532.20 512.19 483.70 3.77 48.50 0.74 9.11 462.64 528.20 483.55 533.60 20.91 5.40 4.52 1.02 461.91 491.39 449.80 515.04 12.11 23.65 2.62 4.81 403.29 487.09 415.54 481.65 12.25 5.44 3.04 1.12 447.94 507.06 452.53 523.94 4.59 16.88 1.02 3.33 445.01 516.14 468.57 521.58 23.56 5.44 5.29 1.05 454.15 478.11 467.89 493.32 13.74 15.21 3.03 3.18 499.62 483.24 475.09 490.14 24.53 6.90 4.91 1.43 452.98 507.23 465.85 545.50 12.87 38.27 2.84 7.54


## 3.2.1. 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 Table6.Table 6
The hierarchy general ranking.
Rule hierarchy Weights Factor hierarchy Weights Total sequencing weight Adjustment parameters 0.6250 Distance from the duct outlet to heading face 0.4673 0.292 Caliber of the duct outlet 0.2772 0.173 Horizontal right deviation 0.1601 0.100 Vertical upward deflection 0.0954 0.060 Original gas and dust parameters 0.3285 Gas emission quantity 0.5715 0.186 Dust production 0.2857 0.094 Particle size 0.1428 0.047 Ventilation parameters 0.1365 Air velocity at the duct outlet 0.75 0.102 Ventilation volume 0.25 0.034From Table6, 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 main hazard caused by dust is impact on the health of the underground workers. Considering the location that has the greatest impact on the workers which is the breathing-height 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. (6)m=n+l+awhere m is the number of nodes in hidden layer, n is the number of nodes in input layer, l is the number of nodes in output layer, and a 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 Table7.Table 7
Error under the different number of nodes in hidden layer.
The number of nodes Error Iteration times The number of nodes Error Iteration times 4 5.76e-2 7 9 6.20e-2 4 5 1.81e-2 16 10 1.01e-1 6 6 1.14e-1 8 11 8.85e-2 18 7 6.60e-2 6 12 5.58e-2 14 8 6.13e-2 18 13 3.11e-2 12It can be seen from Table7 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.Figure 10
The network model structure of dust concentration prediction.


## 3.2.2. Data Sample and Preprocessing

Considering the simulation results of the airflow and dust field under the comprehensive changes of influencing factors, dust concentration at the pedestrian breathing-height in the backflow side is the average at three points of Z=7m, Z=10m, and Z=15m on the blue line (X=5m, Y=1.5m, and Z=0~40m), and dust 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.Table 8
The input layer data.
Number Distance (m) Caliber (m) Horizontal right deviation (°) Vertical upward deflection (°) Air velocity (m/s) Dust production (kg/s) 1 5 1.1 5 2 5.43 0.0076 2 5 1.1 5 6 5.68 0.0076 3 5 1.1 15 2 5.17 0.0076 4 5 1.1 15 6 5.31 0.0076 5 5 1.2 5 6 4.82 0.0076 6 5 1.2 15 2 5.27 0.0076 7 5 1.2 25 2 5.19 0.0076 8 5 1.2 5 2 4.98 0.0076 9 6 1 5 6 5.61 0.0076 ⋮ ⋮ ⋮ ⋮ ⋮ ⋮ ⋮ 76 9 0.9 25 6 6.51 0.0116 77 10 0.7 5 2 6.81 0.0116 78 10 0.7 15 6 6.94 0.0116 79 10 0.7 25 6 6.87 0.0116 80 10 0.7 15 2 7.04 0.0116 81 10 0.8 5 2 6.51 0.0116 82 10 0.8 15 6 6.43 0.0116 83 10 0.8 25 6 6.72 0.0116 84 10 0.8 15 2 6.86 0.0116Table 9
The output layer data.
Number Dust concentration at the Dust concentration at the pedestrian driver location (mg/m3) breathing-height (mg/m3) 1 293.40 280.45 2 307.82 294.34 3 274.41 286.35 4 220.64 270.35 5 326.64 243.83 6 283.89 304.80 7 331.24 318.07 8 322.52 237.90 9 319.91 221.07 ⋮ ⋮ ⋮ 76 508.42 532.20 77 462.64 528.20 78 461.91 491.39 79 403.29 487.09 80 447.94 507.06 81 445.01 516.14 82 454.15 478.11 83 499.62 483.24 84 452.98 507.23In order to facilitate the analysis and the network convergence, the data is normalized by (7). The normalized data are used as training and testing data for the BP prediction model. (7)xi∧=xi-xminxmax-xminwhere xi is the input or output data, xmax is the maximum value, and xmin is the minimum value.


## 3.2.3. Determining Functions of Neural Network

Combined with the sample data in Tables8 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.Table 10
The comparison results of different functions.
Function type Function name Iteration times Mean square error Transfer function logsig 9 2.52e-2 tansig 11 4.65e-2 purelin 9 6.68e-1 Training function traingdm 1000 3.17e-1 trainglm 11 5.57e-2 traingdx 137 8.52e-1 Learning function learngdm 139 7.41e-2 learngd 25 3.46e-2Under 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.


## 3.2.4. 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 Figure11.Figure 11
Training results of dust concentration prediction model. (a) Training interface. (b) Mean squared error. (c) Training state.
(a) (b) (c)It can be seen from Figure11(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 Tables11 and 12, and thresholds of hidden layer nodes are shown in Table 13.Table 11
The weights from input layer nodes to hidden layer nodes.
Input layer nodes X1 X2 X3 X4 X5 X6 Hidden layer nodes 2.4314 -0.0622 -1.2612 1.0935 -2.5342 1.5057 1.3383 -1.0115 2.0125 -0.7830 -0.7608 0.6870 -0.5912 -1.0331 -0.2196 0.1168 -1.4287 -1.9597 2.8285 -4.8427 0.0719 -0.5322 -0.0306 -1.8395 0.6135 -1.1564 0.4116 0.2652 1.8586 -0.3303Table 12
The weights from hidden layer nodes to output layer nodes.
Hidden layer nodes 1 2 3 4 5 Output layer nodes 2.2022 -0.4141 -2.1345 -0.7098 0.9605 -0.3397 -0.3385 -1.9866 -0.4598 -1.0131Table 13
Thresholds of hidden layer nodes.
Hidden layer nodes 1 2 3 4 5 Thresholds -4.2854 -3.2851 0.4871 0.3751 -4.8327In addition, the thresholds of the output layer nodes are 0.82971 and 1.2618.The latter 10 data in Tables8 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.Table 14
The comparison between test value and predicted value.
Test value at the driver location (mg/m3) Test value at the pedestrian breathing-height (mg/m3) Predicted value at the driver location (mg/m3) Predicted value at the pedestrian breathing-height (mg/m3) Absolute error Relative error (%) 455.35 476.21 480.74 508.49 25.39 32.28 5.58 6.78 508.42 532.20 512.19 483.70 3.77 48.50 0.74 9.11 462.64 528.20 483.55 533.60 20.91 5.40 4.52 1.02 461.91 491.39 449.80 515.04 12.11 23.65 2.62 4.81 403.29 487.09 415.54 481.65 12.25 5.44 3.04 1.12 447.94 507.06 452.53 523.94 4.59 16.88 1.02 3.33 445.01 516.14 468.57 521.58 23.56 5.44 5.29 1.05 454.15 478.11 467.89 493.32 13.74 15.21 3.03 3.18 499.62 483.24 475.09 490.14 24.53 6.90 4.91 1.43 452.98 507.23 465.85 545.50 12.87 38.27 2.84 7.54


## 4. Discussion

### 4.1. 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 Table15.Table 15
The adjustment schemes at the nearest distance of 5m.
Scheme Caliber (m) Horizontal right deviation (°) Vertical upward deflection (°) 1 1.1 5 2 2 1.1 5 4 3 1.1 5 6 4 1.1 10 2 ⋮ ⋮ ⋮ ⋮ 28 1.2 25 2 29 1.2 25 4 30 1.2 25 6According 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 Table16).Table 16
Dust concentration under different adjustment schemes at the nearest distance of 5m.
Scheme Dust concentration at the Dust concentration at the pedestrian driver location (mg/m3) breathing-height (mg/m3) 1 329.07 342.38 2 347.82 363.36 3 411.41 429.07 4 375.64 363.06 ⋮ ⋮ ⋮ 28 331.76 429.29 29 334.04 409.71 30 345.38 457.73According to scheme 9 from Table16, the optimal adjustment scheme at the nearest distance of 5m is obtained. Dust concentration at the driver location is 325.36 mg/m3 after the adjustment of scheme 9, and dust concentration at the pedestrian breathing-height is 304.32 mg/m3, 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).Table 17
Analysis of adjustment effect at the nearest distance of 5m.
Predictive object Dust concentration at the Dust concentration at the pedestrian driver location (mg/m3) breathing-height (mg/m3) Before adjustment 517 447 After adjustment 325.36 304.32 Reduced percentage 37 32At 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/m3, which is decreased by 35%, and dust concentration at the pedestrian breathing-height is 323.49 mg/m3, 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°.


### 4.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 breathing-height 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 Figure12.Figure 12
Arrangement of underground measuring points.(  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).Table 18
Test results at the nearest distance of 5m before adjustment.
Measuring object Dust concentration (mg/m3) Test section  1# 2# 3# 4# 5# 6# 7# 8# Measuring points 1 486.2 431.6 364.4 321.3 284.6 256.2 239.4 205.8 2 541.1 471.5 369.6 310.8 233.1 225.8 233.1 215.3Table 19
Test results at the nearest distance of 5m after adjustment.
Measuring object Dust concentration (mg/m3) Test section  1# 2# 3# 4# 5# 6# 7# 8# Measuring points 1 347.3 309.8 300.3 239.4 223.7 187.9 183.8 176.4 2 357.4 322.4 189.5 203.7 177.3 145.5 126.8 115.5Figure 13
Comparison of dust concentration before and after adjustment at the nearest distance of 5m. (a) Dust concentration along the pedestrian breathing-height. (b) Dust concentration along the driver’s breathing-height.
(a) (b)From Tables18 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).Table 20
Test results at the farthest distance of 10m before adjustment.
Measuring object Dust concentration (mg/m3) Test section  1# 2# 3# 4# 5# 6# 7# 8# Measuring points 1 504.2 443.1 384.3 686.7 592.2 501.9 478.8 510.3 2 590.6 436.8 396.9 499.8 558.6 661.5 575.4 537.6Table 21
Test results at the farthest distance of 10m after adjustment.
Measuring object Dust concentration (mg/m3) Test section 1# 2# 3# 4# 5# 6# 7# 8# Measuring points 1 346.2 291.9 323.4 354.1 256.2 278.2 255.7 265.5 2 397.8 322.1 352.8 304.5 266.7 277.2 220.5 222.6Figure 14
Comparison of dust concentration before and after adjustment at the farthest distance of 10m. (a) Dust concentration along the pedestrian breathing-height. (b) Dust concentration along the driver’s breathing-height.
(a) (b)It can be seen from Tables20 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%.


## 4.1. 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 Table15.Table 15
The adjustment schemes at the nearest distance of 5m.
Scheme Caliber (m) Horizontal right deviation (°) Vertical upward deflection (°) 1 1.1 5 2 2 1.1 5 4 3 1.1 5 6 4 1.1 10 2 ⋮ ⋮ ⋮ ⋮ 28 1.2 25 2 29 1.2 25 4 30 1.2 25 6According 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 Table16).Table 16
Dust concentration under different adjustment schemes at the nearest distance of 5m.
Scheme Dust concentration at the Dust concentration at the pedestrian driver location (mg/m3) breathing-height (mg/m3) 1 329.07 342.38 2 347.82 363.36 3 411.41 429.07 4 375.64 363.06 ⋮ ⋮ ⋮ 28 331.76 429.29 29 334.04 409.71 30 345.38 457.73According to scheme 9 from Table16, the optimal adjustment scheme at the nearest distance of 5m is obtained. Dust concentration at the driver location is 325.36 mg/m3 after the adjustment of scheme 9, and dust concentration at the pedestrian breathing-height is 304.32 mg/m3, 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).Table 17
Analysis of adjustment effect at the nearest distance of 5m.
Predictive object Dust concentration at the Dust concentration at the pedestrian driver location (mg/m3) breathing-height (mg/m3) Before adjustment 517 447 After adjustment 325.36 304.32 Reduced percentage 37 32At 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/m3, which is decreased by 35%, and dust concentration at the pedestrian breathing-height is 323.49 mg/m3, 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°.


## 4.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 breathing-height 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 Figure12.Figure 12
Arrangement of underground measuring points.(  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).Table 18
Test results at the nearest distance of 5m before adjustment.
Measuring object Dust concentration (mg/m3) Test section  1# 2# 3# 4# 5# 6# 7# 8# Measuring points 1 486.2 431.6 364.4 321.3 284.6 256.2 239.4 205.8 2 541.1 471.5 369.6 310.8 233.1 225.8 233.1 215.3Table 19
Test results at the nearest distance of 5m after adjustment.
Measuring object Dust concentration (mg/m3) Test section  1# 2# 3# 4# 5# 6# 7# 8# Measuring points 1 347.3 309.8 300.3 239.4 223.7 187.9 183.8 176.4 2 357.4 322.4 189.5 203.7 177.3 145.5 126.8 115.5Figure 13
Comparison of dust concentration before and after adjustment at the nearest distance of 5m. (a) Dust concentration along the pedestrian breathing-height. (b) Dust concentration along the driver’s breathing-height.
(a) (b)From Tables18 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).Table 20
Test results at the farthest distance of 10m before adjustment.
Measuring object Dust concentration (mg/m3) Test section  1# 2# 3# 4# 5# 6# 7# 8# Measuring points 1 504.2 443.1 384.3 686.7 592.2 501.9 478.8 510.3 2 590.6 436.8 396.9 499.8 558.6 661.5 575.4 537.6Table 21
Test results at the farthest distance of 10m after adjustment.
Measuring object Dust concentration (mg/m3) Test section 1# 2# 3# 4# 5# 6# 7# 8# Measuring points 1 346.2 291.9 323.4 354.1 256.2 278.2 255.7 265.5 2 397.8 322.1 352.8 304.5 266.7 277.2 220.5 222.6Figure 14
Comparison of dust concentration before and after adjustment at the farthest distance of 10m. (a) Dust concentration along the pedestrian breathing-height. (b) Dust concentration along the driver’s breathing-height.
(a) (b)It can be seen from Tables20 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%.


## 5. 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 double-objective 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 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/m3 to 357.4 mg/m3, which is decreased by 34%, and dust concentration at the pedestrian breathing-height in the backflow side decreased from 486.2 mg/m3 to 347.3 mg/m3, 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/m3, 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/m3; the dust removal effect reached 31% compared with that before adjustment.



---

*Source: 6210704-2019-07-21.xml*