^{1}

^{*}

^{1}

The traditional research considers that the effective range of forced local ventilation in single head tunnel is the horizontal distance from the jet exit section to the rotary section, but there is a low wind speed area at the end of the range which can not completely and effectively release toxic and harmful gases. In this area, due to the small jet energy, there is no way to effectively remove the toxic and harmful gases produced in the driving process, which seriously affects the safety of the mine . Therefore, the reasonable definition of the effective range of the forced local ventilation of the single head tunnel can not only improve the theory of the local ventilation of the mine, but also prevent the accumulation of gas and ensure the life safety of the employees. In this paper, the definition of effective range of local ventilation in single head heading face is put forward, and the rationality of the definition of effective range is verified by numerical simulation. On the basis of reasonable definition of effective range, the relationship between wind speed and effective range under different ventilation modes is studied by numerical simulation. The results show that when the wind speed at the outlet of the air duct is less than 20 m/s, the effective range changes rapidly with the outlet wind speed, and after more than 20 m/s, the effect of wind speed on the effective range is gradually weakened. Under the same exit wind speed, the effective range of parallel and coplanar air inlet is much greater than that of single air inlet.

With the development of China’s social economy, the demand for coal resources is also increasing. In order to meet the needs of development, China’s coal mining enterprises are also constantly exploring potential and improving production capacity. Coal mine gas disaster seriously restricts the production capacity of China’s coal mines. Only by controlling gas accidents and stabilizing the overall situation of safety production, can we ensure the healthy development of China’s coal industry. In the process of mining, the gas concentration in the accumulation area is high and widely distributed, which is easy to cause gas explosion and gas combustion accidents, and it is a major hazard source of the mine [

Since the 1980s, scholars at home and abroad have carried out a lot of experimental research on the air flow state of driving face. Internationally, Nakayama and other scholars in Japan have obtained the distribution of air flow in the driving tunnel by using the method of laboratory measurement, and Ichinose Masatomo and other scholars have analyzed the gas distribution in the driving tunnel by using the method of field measurement in Japan [

In China, Wang Haiqiao et al. based on the hydrodynamics and jet theory, deeply analyzed the confined attached jet ventilation, elaborated the generation process and air flow structure of the driving ventilation jet, and provided the calculation formula of the average wind speed, maximum return speed and jet action distance of the air flow in the return area of the working face space flow field [

To sum up, previous studies on the effective range of the pressure type local ventilation in the single head tunnel have found that when the wall jet in the single head tunnel develops to a certain extent, part of the air in the jet begins to return. This section is called the rotary section. The jet flow of the rotary section is the largest. After the rotary section, the flow gradually decreases. When the jet reaches the maximum range, the air flow is full the wind speed is zero [

The minimum jet energy in the effective range area must be able to effectively eliminate the toxic and harmful gases produced in the process of driving, so as to ensure the basic requirements of safety and health of the operators. Therefore, on the basis of the previous definition of range, the requirements of the coal mine safety regulations on the lower limit of wind speed are added, that is, the minimum wind speed of the coal roadway in the driving face is 0.25 m/s, and the minimum wind speed of the rock roadway is 0.15 m/s. The low wind speed area that does not meet the requirements is removed, and the remaining range area length is defined as the effective range. In short, the effective range is the horizontal distance between the intersection point of the isoline of the lowest wind speed and the center line of the tunnel and the jet outlet.

In order to prove the scientificity of the definition of effective range, the effective range and gas distribution of single air duct forced local ventilation in single head roadway are simulated by using fluid dynamics fluent simulation software.

1) Simulation analysis model of forced local ventilation in single head roadway

A physical model of single air duct forced ventilation in a single head roadway is established. The air duct is hung at the middle waist of the side wall of the roadway. The cross section size of the roadway is 3.3 m × 5 m, the length is 200 m, the diameter of the air duct is 1 m, and the length is 150 m.

2) Mathematical model

Considering that the air flow in driving tunnel is generally turbulent, the mass conservation equation, momentum conservation equation, standard k − ε turbulence model equation and component transport equation are established, which are suitable for turbulence.

3) Setting of boundary conditions

The exit of the air duct is defined as vehicle inlet; the exit of the roadway is defined as outflow; the boundary condition around the roadway is set as wall. The fluid area is the space occupied by the tunnel to remove the air duct, and UDF files are compiled for the gas distribution during the coal tunnel driving. Assuming that the gas emission from the driving face, the roof and the side wall of the roadway is uniform at a distance of 50 m from the excavation head, the calculation formula of the gas emission source term is as follows:

S i = C Q ρ φ 60 V

Formula 1 Calculation formula of gas emission source term

S_{i}: gas gushing source term, kg·m^{−3}·s^{−1};

C: gas concentration, %;

Q: air volume of tunnel, m^{3}/min;

ρ: gas density, 0.719 kg/m^{3};

V: volume of gas gushing source, m^{3};

φ: proportional coefficient of gas gushing source term.

According to the actual measurement of a mine, it can be concluded that q = 1200 m^{3}/min, C = 0.83%, ρ = 0.719 kg/m^{3}, φ = 1, v = 58.536 m^{3}, in which V is calculated by the sum of the volume with the thickness of 0.1 M of the head coal wall, the roof of the tunnel and the side wall of the tunnel. Substituting formula, S_{i} = 0.002038 kg·M^{−3}·S^{−1} is obtained.

4) Mesh division

According to the characteristics of the model, tetrahedron is used to divide the mesh, and the number of mesh is about 120,000.

5) Numerical simulation

The unsteady state is used for simulation, and the time step is set as 20 s. Since the model is a three-dimensional model, in order to study the gas migration in the tunnel, a horizontal section is made from the center of the air duct. The contour distribution of the wind speed above 0.25 m/s is shown in

6) Result analysis

From

can be verified that the definition of effective range is in accordance with the actual law.

According to the three different air duct arrangements of single air duct, double air ducts in common plane and two air ducts in parallel with each other with a distance of 5 m, the relationship between the wind speed and the effective range of the driving face under different air duct arrangements is analyzed by using fluent fluid simulation software.

1) Establishment of physical model

The air duct of the forced ventilation is generally arranged on the side wall of the single head roadway, and the outlet of the air duct is generally located in the middle or on the side wall of the roadway, ensuring a certain distance from the head-on. In order to analyze the relationship between wind speed and effective range, a simplified model of wind speed and effective range is established (

2) Mathematical model

Arrangement of air duct | Size of tunnel (width × height × length) | Fan drum length and diameter | Model section diagram |
---|---|---|---|

Single fan drum wall-attached jet ventilation | 5m × 3.3 m × 200 m | length 100 m diameter 1 m | |

Parallel double fan drums wall-attached jet ventilation (The difference between the two air ducts is 5 m) | 5 m × 3.3 m × 200 m | length 100 m and 95 m, diameter 1 m | |

Double fan drums coplanar wall-attached jet ventilation (The two air ducts are the same length) | 5 m × 3.3 m × 200 m | length 100 m and 100 m, diameter 1 m |

In order to simplify the study, the following assumptions are made: the fluid in the tunnel can be regarded as incompressible fluid, ignoring the heat dissipation caused by the work of fluid viscosity, assuming the wall is adiabatic and isothermal ventilation; assuming the flow field has a high turbulence Reynolds number, the turbulence viscosity of the fluid is isotropic, and the turbulence viscosity coefficient can be treated as scalar; the flow is steady-state turbulence, meeting the Boussinesq false Set up. Therefore, the mathematical models of wind speed field simulation are mass conservation equation, momentum conservation equation and standard k − ε turbulence model equation.

3) Boundary conditions

The entrance boundary of the model is defined as the exit of the forced draft fan, the exit boundary is defined as the roadway section at the back of the fan, and the exit velocity of the fan is v. Therefore, the boundary conditions of the model can be described as follows:

Inlet boundary: V x = 0 , V y = 0 , V z = v , turbulent kinetic energy K = 3 2 ( v I ) 2 , turbulent kinetic energy dissipation rate ε = C u 3 4 k 3 2 l , in the formula, I = 0.16 ( Re D H ) − 1 8 , l = 0.07 L , L is equivalent diameter, C u = 0.084 .

Outlet boundary: ∂ V i ∂ X i = 0 , P = P o u t relative pressure P = 0 .

Wall boundary: All of the wall is non-slip boundary, V i = 0 ; Wall is adiabatic; Pressure gradient which is vertical to wall surface is zero.

4) Grid division

According to the characteristics of the model, the Boolean operation is carried out, and the fluid region is obtained by subtracting the model of the air duct from the tunnel model. Under the Cartesian coordinate system, the calculation grid of the fluid region is divided, and the solution region is divided into tetrahedral grid. The length of the tunnel (i.e. the length of the calculation area) is l, the width is B, the height is h, and the perimeter of the air duct is c. the grid density of the calculation area is L/0.5, the width direction is B/0.1, the height direction is H/0.1, and the perimeter position is C/0.1. The finite element grid is generated by mapping grid. The tetrahedron mesh is used to divide the single air duct model, and the mesh number is about 30,000; the tetrahedron and hexahedron hybrid mesh is used to divide the parallel 5 m model of the double air ducts, and the mesh number is about 50,000; the tetrahedron and hexahedron hybrid mesh is used to divide the coplanar model of the double air ducts, and the mesh number is about 50,000.

5) Analysis of simulation results

According to the above-mentioned local ventilation model of single head tunnel, when the wind speed at the outlet of the air duct is 4, 6, 8, 10, 12, 14, 16, 18, 20, 30 and 40 m/s respectively, the relationship between the jet wind speed at the outlet of the air duct and the effective range under different air duct layout is studied. The simulation results are shown in

According to the effective range obtained from

Weed speed m/s | Single fan drum effective range (m) | Parallel double fan drums effective range (m) | Double fan drums coplanar effective range (m) | |||
---|---|---|---|---|---|---|

More than 0.15 m/s | More than 0.25 m/s | More than 0.15 m/s | More than 0.25 m/s | More than 0.15 m/s | More than 0.25 m/s | |

4 | 18.1 | 12.4 | 19.3 | 19.16 | 17.7 | 13.54 |

6 | 23.5 | 14.4 | 22.96 | 19.21 | 27.24 | 17.11 |

8 | 25.9 | 19.3 | 29.2 | 20.31 | 33.2 | 20.79 |

10 | 27.9 | 20.7 | 35.15 | 24.57 | 34.24 | 27.26 |

12 | 30 | 26.5 | 39.25 | 29.09 | 42.9 | 33 |

14 | 32 | 28.6 | 42.14 | 35.01 | 43.76 | 34.01 |

16 | 31.9 | 28 | 43.7 | 38.34 | 44.37 | 34.19 |

18 | 35.1 | 29.7 | 45.24 | 39.4 | 45.14 | 37.67 |

20 | 35.4 | 28.3 | 47.67 | 41 | 47.46 | 43.05 |

30 | 36.1 | 33.4 | 49.87 | 43.65 | 55.2 | 47.15 |

40 | 36.4 | 35.4 | 49.6 | 43.04 | 59.34 | 48.64 |

Fan drum arrangement mode | The lowest weed speed isoline | Fitting curves | Correlation coefficient R 2 |
---|---|---|---|

Single fan drum | 0.15 m/s | y = 0.0009 x 3 − 0.0811 x 2 + 2.5332 x + 10.004 | 0.9869 |

0.25 m/s | y = 0.0011 x 3 − 0.0983 x 2 + 2.9786 x + 1.1508 | 0.9693 | |

Parallel double fan drums | 0.15 m/s | y = 0.0014 x 3 − 0.135 x 2 + 4.3948 x + 2.7193 | 0.9943 |

0.25 m/s | y = − 0.0002 x 3 − 0.0241 x 2 + 2.2148 x + 7.5535 | 0.953 | |

Double fan drums coplanar | 0.15 m/s | y = 0.0018 x 3 − 0.1504 x 2 + 4.543 x + 3.8194 | 0.9767 |

0.25 m/s | y = 0.0005 x 3 − 0.0732 x 2 + 3.266 x + 0.8003 | 0.9837 |

The above fitting curves are obtained by cubic polynomial fitting, and the correlation coefficient R^{2} is greater than 0.95, so the fitting has a high correlation. From the above fitting curve, it can be seen that when the wind speed at the outlet of the air duct reaches 20 m/s, the effective range gradually keeps stable, that is, the influence of the wind speed on the effective range gradually weakens.

The effective range of forced ventilation is redefined by the contour of the minimum allowable wind speed in the horizontal section of the tunnel center, and the numerical simulation method of the relationship between the effective range of forced local ventilation and the jet wind speed at the outlet of the single air duct, the double air duct and the coplanar and parallel air duct are put forward, and the corresponding relationship model is established.

The rationality and scientificity of the definition of effective range are verified by numerical simulation. It is found that the gas concentration isoline in front of the effective range position is very dense, and the average gas concentration rises sharply; while in the back, the gas concentration isoline is sparse, and the average gas concentration drops sharply. It shows that the redefined effective range can effectively eliminate toxic and harmful gases, and can meet the requirements of safety production.

When the wind speed at the outlet of the air duct is less than 20 m/s, the effective range changes faster with the outlet wind speed. When the wind speed exceeds 20 m/s, the effective range changes less with the speed, that is, the effect of wind speed on the effective range decreases gradually. Under the same exit wind speed, the effective range of parallel and coplanar air inlet is much greater than that of single air inlet.

The authors declare no conflicts of interest regarding the publication of this paper.

Hao, Y.W. and Wang, Y. (2020) Study on the Effective Range of Local Ventilation in the Driving Face of Single Head Roadway. World Journal of Engineering and Technology, 8, 50-59. https://doi.org/10.4236/wjet.2020.81006