Abstract

In U-shaped ventilation working face, different tunnel section shapes are one of the important factors affecting the propagation of gas explosion shock wave. In order to study the propagation law of gas explosion shock wave in working face, the numerical simulation study was carried out by using Fluent simulation software combined with the actual situation of gas explosion in #415 working face of Chenjiashan Coal Mine in Shaanxi Province. By constructing a three-dimensional mathematical and physical model, a simulation study of the upper-corner gas explosion was carried out. The results are described as follows. (1) After the gas explosion shock wave propagates 40 m, the overpressure peak equidistant difference tends to be stable and attenuates and propagates in the form of a single shock wave. The study determines that the effective length of the U-shaped ventilation inlet/return tunnel is 40 m. (2) When the tunnel section is trapezoidal, the initial overpressure of the gas explosion shock wave propagating to the inlet/return airway is the highest, followed by rectangular and semicircular arches, but the internal overpressure attenuation trend of different cross-sectional shapes is the same. (3) The gas explosion shock wave propagates radially along the working face section during the working face propagation. The farther away the location is from the upper corner of the tunnel during a gas explosion with different cross-sectional shapes, the closer the cutoff overpressure peak is. The attenuation trend of overpressure with the propagation distance conforms to the power function law. The research results provide an important theoretical direction for the numerical simulation of gas explosions in coal mining faces.

1. Introduction

With the increase in coal mining depth, underground gas emission and gas accumulation are becoming more and more serious. The upper corner of a coal mining face is the hazardous area of gas accumulation in coal mine and is the hot spot for gas explosion accidents in coal mines. After the gas explosion accident in upper corner, the huge gas explosion shock wave will spread along the inlet and outlet tunnel, causing serious damage to underground facilities and miners.

The propagation process of explosion shock wave is a process of constant conversion between energy and the external environment. In this process, as the energy of explosion shock wave is gradually converted into gas energy, the propagation length of explosion shock wave is limited. Therefore, studying the effective length can determine the energy dissipation capacity of explosion shock wave propagation. In addition, the cross-section shape is also one of the factors affecting the propagation of the blast wave. Different cross-section shapes have an impact on the propagation speed of the blast wave and the external energy conversion capacity.

Scholars have conducted a significant amount of research on the propagation law of gas explosion shock waves. Both et al. adjusted and optimized the grid model parameters of turbulence and combustion by simulating gas explosions and significantly improved the prediction of the development model [1]. Jing et al. deduced the relationship between the peak value of overpressure and the gas velocity and propagation distance in the process of a gas explosion [2]. Some scholars studied the influence of various conditions (initial temperature, ignition conditions, gas concentration, and barrier type) on gas explosion parameters [312]. Some scholars obtained the mechanism of suppressing a gas explosion by ultrasonic water mist and a vacuum chamber through experiments and numerical simulations [13, 14].

Some scholars studied the distribution and propagation law of gas explosion shock waves in different media of pipelines [1518]. Some scholars studied the influences of pipe length, pipe diameter, and space distance on gas explosion overpressure [1921]. Li and Hao simulated the ignition and propagation characteristics of gas explosions by introducing basic chemical reactions [22]. Some scholars studied the propagation law of gas explosion shock waves in roadways with different turning angles and on bifurcated roadways [2327]. Some scholars simulated and studied the propagation law of gas explosion shock waves in the underground parallel roadway network [28, 29]. Some scholars simulated and analyzed the propagation law of gas explosion shock waves in coal mining face roadways [30, 31].

Many scholars have focused on the pipeline shock wave propagation law and its influencing factors in gas explosion studies, but the actual tunnel effective length and cross-section shape impacts of the gas explosion shock wave law have not been reported. The explosion energy generated by gas accumulation under different geological conditions and different gas pressures is different, so the effective length of blast wave propagation is also different. Studying the effective length of blast wave propagation is beneficial to determine the energy exchange ability of the energy generated by explosion with the outside world. Therefore, the author uses Fluent simulation software, combined with the actual situation of a gas explosion in the upper corner of the #415 working face in Chenjiashan coal mine, Shaanxi Province. To determine the effective length of the air inlet/return roadway, a numerical simulation of an upper-corner gas explosion of different cross-sectional shapes of a U-type ventilation roadway is carried out. The research results provide theoretical guidance for the prevention and control of gas explosion disasters in coal mining faces.

2. Construction and Verification of Mathematical Model

2.1. Basic Assumptions and Boundary Conditions of Numerical Simulation

Gas explosion is a very complex fluid elastic-plastic process, accompanied by chemical reaction, turbulence changes, and other phenomena, so it is difficult to simulate the law of mine gas explosion after considering all conditions. In order to simulate the law of mine gas explosion as much as possible and ensure the reliability of numerical simulation of gas explosion, some basic assumptions must be made:(1)Both premixed gas and combustion products satisfy the ideal gas state equation(2)The specific heat capacity of the mixed gas follows the mixing rule, and the specific heat capacity of each component is a function of temperature(3)The wall of the physical model is a rigid adiabatic wall, and there is no relative displacement(4)The gas explosion reaction is one-way and irreversible(5)In the physical model, the gas filling area is a normal uniform mixture of gas and air, and it is an ideal state

The model parameters and boundary conditions of gas explosion numerical simulation are as follows. The boundary conditions: the pipe boundary is set as the adiabatic wall with the temperature of 300 K, and the outlet above the vertical pipe is set as the pressure outlet.

Basic model parameter: the initial conditions of the burned zone are T = 1600 K, H2O volume fraction is 0.118, CO2 volume fraction is 0.145, and initial pressure is 101325 Pa; the initial conditions of unburned zone are as follows: T = 300 K, CH4 volume fraction is 0.053, O2 volume fraction is 0.21, and H2O and CO2 volume fraction is 0; the initial conditions of the air zone are as follows: the volume fraction of CH4, H2O, and CO2 is 0 and the volume fraction of O2 is 0.233.

2.2. Construction of Gas Explosion Mathematical Model

A gas explosion is a combustion reaction process, and the reaction is rapid. Assuming that the mixture of methane and air is an ideal gas, the hydrodynamic mathematical model of methane gas is established without considering the unstable factors such as viscosity, heat conduction, diffusion, and detonation; the conservation equation of the process is described as follows.

Equation of state is as follows:

Continuity equation is as follows:

Energy equation is as follows:

Momentum equation is as follows:where P is the pressure; t is the time; x, y, and z are rectangular coordinate system parameters; u, , and are the velocities in three coordinate directions; is the fluid density; T is the temperature; R is the gas constant; e is the specific energy; and (γ is the gas index).

2.3. Turbulence Model and Combustion Model

Gas explosion will lead to turbulent combustion, so the turbulence model should be selected in numerical simulation, and the correlation terms and turbulence viscosity coefficients can be calculated by constructing differential equations. According to the number of differential equations required for solution, the turbulence model can be divided into zero equation model, single equation model, and two-side equation model [32]. The large eddy numerical model can better simulate the required results by referring to the literature and several times of practical calculation and debugging [3335]. Therefore, the LES turbulence model is used in this paper to describe the characteristics of turbulent flow field in combustion process.

Chemical reaction rates can be calculated using four models: laminar finite rate model, laminar finite rate/vortex dissipation model, vortex diffusion model, and vortex diffusion conceptual model.

The laminar flow finite velocity model ignores the influence of turbulent pulsation and uses Arrhenius to calculate the combustion rate. The condition is laminar flame, which is not suitable for turbulent combustion. The laminar finite velocity/eddy current dissipation model is characterized by the fact that the reaction rate is smaller than Arrhenius and eddy current dissipation equation, which is suitable for single-stage and two-stage reaction mechanism but not for multistage reaction mechanism. The reaction rate of the eddy current dissipation model is controlled by turbulent mixing for non-premixed flames. Combining the advantages and disadvantages of the four models, the vortex dissipation vortex diffusion model suitable for the LES turbulence model is adopted in this paper.

2.4. Verification of Calculation Model

According to previous gas explosion experiments in pipelines [36], a full-scale simulation comparison of 1 : 1 was carried out using Fluent simulation software to verify the reliability of the large eddy simulation turbulence equation and eddy-dissipation vortex diffusion model. In the test, the width and height of the straight pipe are both 80 mm, the length is 19.2 m, and the length of the vertical pipe is 5 m. According to the setting of the experimental monitoring points, measuring point 1 is 19 m away from the leftmost end and is located on the central line of the horizontal pipeline, and measuring point 2 is located 0.5 m away from the central line of the horizontal pipeline. A section diagram of the physical model and the locations of the measuring points are shown in Figure 1. The left end of the pipeline is closed and filled with gas, and the pressure outlet is above the vertical pipeline.

The initial conditions and boundary conditions are consistent with those in the gas explosion experiment [37]; specifically, the gas concentration in the gas filling area is 0.053; O2 concentration is 0.212; CO2 concentration in the air area is 0.1456; and H2O concentration is 0.11925. The pipeline boundary is set as an adiabatic wall, the temperature is 300 K, and the outlet above the vertical pipe is set as the pressure outlet.

The selected calculation model is as follows: the fluid flow is unsteady turbulent flow, the turbulence model adopts the large eddy simulation turbulence equation, the wall adopts a standard wall function, the chemical reaction of methane combustion adopts an eddy-dissipation vortex diffusion model that is suitable for turbulent combustion, and a simple algorithm is used for an iterative solution. The simple algorithm is a widely used numerical method to solve flow field in computational fluid mechanics. It can not only solve nonlinear problems but also complete the coupling of velocity and pressure. Therefore, the simple algorithm is used for iterative solution.

The experimental and simulation results are shown in Figure 2.

As can be seen from Figure 2, in measuring point 1, when the inflatable length is, respectively, 4 m, 5.5 m, and 7 m, the errors between the simulation results and the mean values of the four groups of experimental data are 2.9%, 7.5%, and 10.9%, respectively. In point 2, when the inflating length is 4 m, 5.5 m, and 7 m, respectively, the errors between the simulation results and the average values of the four groups of experimental data are 1.9%, 9.5%, and 5.4%, respectively. It can be concluded that the maximum error between the simulated results and the average values of the four groups of experimental data is 10.9% and the minimum error is only 1.9%. When comparing the numerical simulation data with the experimental data, due to the limitation of the numerical simulation conditions, there will be some error with the actual experimental data. Combined with Figure 2, when the inflatable lengths of pipes are 4 m, 5.5 m, and 7 m, respectively, the overpressure peak errors obtained by numerical simulation and experiment are all less than 15%, which are within the allowable range of international engineering error standards. The reliability of the LES + ED model used in gas explosion simulation is verified experimentally.

3. Determination of Effective Length of Intake/Return Air Roadway

The intake/return air roadway of the #415 working face in Chenjiashan coal mine of Shaanxi Province is arranged along the coal seam floor. The length of the r inlet roadway is 1192 m, and the length of the return air roadway is 1170 m. The cross-section of roadway is a semicircular arch with a parallel arrangement, and the inclined length of the working face is 175 m. In the numerical simulation study, the reasonable length of the intake/return air roadway not only makes the calculation results more accurate but also significantly reduces the workload of numerical calculation. Therefore, the author carried out a simulation study on the gas explosion with U-type ventilation to determine the reasonable effective length of the inlet/outlet roadway.

3.1. Establishing the 3D Model and Meshing

According to the specific parameters of the #415 working face in Chenjiashan coal mine of Shaanxi Province and the actual situation of the explosion, geometry and mesh modules in workbench are used to establish a three-dimensional physical model and grid division. As shown in Figure 3, the length of the inlet/return air roadways is 100 m, and the mesh size is 0.5 m. The mesh is refined by local mesh densification technology in the gas filling area and ignition point. As shown in Figure 3, the length of the return air roadway is 100 m, and 9 monitoring points are arranged in the return air roadway, each with a distance of 10 m.

3.2. Initial Conditions and Boundary Conditions

The explosion site of the November 28th major gas explosion accident in Chenjiashan coal mine was the upper corner of the #415 working face. The inlet and return air outlets are set as pressure outlets, and the others are standard wall surfaces. According to the accident analysis report of the government department of Chenjiashan coal mine gas explosion accident, the mass fractions of CH4 and O2 in the gas filling area of the upper corner are 0.095 and 0.212, respectively, and the burnt area is a ball with a radius of 0.5 m in the center of the upper corner. In this area, the mass fraction of CO2 is 0.1456, the mass fraction of H2O is 0.11925, the standard atmospheric pressure is 101325 Pa, and the temperature is 1600 K. A simple algorithm is used for an iterative solution; nine monitoring points are arranged in the return air roadway, each of which is 10 m apart.

3.3. Analysis of Numerical Simulation Results

The whole explosion can be regarded as a process that starts from the ignition zone, detonates the gas/air mixture in turn, and propagates from the combustion wave and shock wave to the unburned area. This can be divided into three stages: ignition stage, compression combustion stage, and single shock wave propagation stage [37]. The propagation cloud diagrams of blast wave with different time scales shown in Figure 4 are given based on simple algorithms and simulation results, respectively, describing the propagation characteristics of blast wave at different moments. The propagation nephogram of the explosion shock wave and the pressure curve of each measuring point are shown in Figures 4 and 5.

The overpressure peak value, peak time, and variation of overpressure peak value between adjacent measuring points and equidistant attenuation coefficient are listed in Table 1.

According to the parameter data in Table 1, the variation of the overpressure peak value and arrival time with distance is shown in Figure 6. The variation of the overpressure peak’s 10 m isometric difference between adjacent measuring points with distance is shown in Figure 7. With an increase in distance, the attenuation coefficient of the overpressure peak and data fitting prediction of each measuring point is shown in Figure 8.(1)Figure 6 shows that the peak value of overpressure varies greatly within 40 m from the explosion area. The overpressure peak value of each measuring point basically attenuates linearly beyond 40 m from the explosion area. The time to reach the peak value increases linearly with an increase in distance.(2)Figure 7 shows that the difference in the overpressure peak value of each measuring point at measuring points #1–4 is relatively large, and the difference in the overpressure peak value between adjacent measuring points at measuring points #4–9 basically fluctuates up and down with a very small amplitude. That is, starting from measuring point #4 (40 m away from the explosion area), the maximum overpressure attenuation is approximately 0.02 MPa for every 10 m propagation of the gas explosion shock wave.(3)Figure 8 shows that the attenuation coefficient of the overpressure peak value of measuring point #4 (40 m away from the explosion area) to measuring point #9 (90 m away from the explosion area) presents an obvious linear increasing relationship, and its functional relationship is y = 0.0066x + 0.1244.

According to the results in Figures 68, the gas explosion shock wave starts to propagate stably in the form of a single shock wave from a distance of 40 m from the explosion area. Therefore, the length of intake/return air roadway for the gas explosion simulation with U-type ventilation is determined to be 40 m. When the length of U-shaped gas explosion simulation inlet and return air tunnel is determined to be 40 m, we can provide a certain reference for the specific location of explosion-proof door according to the effective length and choose different support modes to ensure the stability of roadway structure. In addition, it not only ensures the accuracy of the calculation results but also greatly reduces the workload of numerical simulation.

4. Study of Influence of Different Roadway Cross-Section Shapes on Shock Wave Law of Gas Explosion

According to different mine types, mining methods, geological conditions, hydrological conditions, and other factors, the roadway adopts different sectional shapes. The most commonly used cross-sectional shapes of roadways in China are rectangular, trapezoidal, and semicircular arch. Therefore, in the study of U-shaped ventilation gas explosion with different cross-section shapes, rectangle, trapezoid, and semicircle arch are adopted, respectively, as shown in Figure 9.

4.1. Physical Model and Initial/Boundary Conditions

Three-dimensional physical models of the intake/return air roadway and working face with rectangular, trapezoidal, and semicircular arch cross-sections are established as shown in Figure 10. The vertical height of the roadway is 0.3 m higher than that of the working face when the cross-sectional shape is a semicircular arch, and the roadway height is the same as that of the working face when the cross-sectional shape is rectangular or trapezoidal. The length of the inlet/return air roadway is 40 m. The initial and boundary conditions of the numerical simulation are the same as those in Section 3.2.

4.2. Analysis of Numerical Simulation Results

The numerical simulation of different cross-section shapes of roadway is carried out, the upper-corner gas is detonated by high temperature ignition point, the rectangular section is taken as an example, and the pressure nephogram of gas explosion shock wave is obtained at different times as shown in Figure 11. The pressure cloud image of gas explosion shock wave at typical time was selected to better explain the propagation characteristics of explosion shock wave in the upper corner of the U-shaped ventilation model.

According to the above pressure nephogram of gas explosion shock wave at different times, it can be seen that after the gas explosion occurs in the upper corner of U-type ventilation mode, the shock wave first propagates through the return air roadway and coal mining face, and the shock wave pressure gradually attenuates in the process of propagation. When the shock wave propagates to the corner where the working face is connected with the air inlet roadway, it can be seen that the shock wave is reflected and superimposed, which causes the overpressure of the shock wave to rise in a very short time.

The data of monitoring points are sorted out to get the change of overpressure at each point of 40 m inlet/return air roadway with different cross-section shapes, as shown in Figures 12 and 13.

Figures 12 and 13 show that when the upper-corner gas explosion occurs, the initial overpressure of the inlet/return air roadways with different cross-sectional shapes is different. When the cross-sectional shape is a trapezoid, the initial overpressure of the explosion shock wave propagating to the inlet/return air roadway is the largest, followed by the rectangular shape, and the minimum occurs with the semicircular arch. When the explosion shock wave propagates in roadways of different cross-sectional shapes, it is consistent with the initial overpressure at all positions in the roadway. The results show that the maximum overpressure occurs with the trapezoidal shape, the second is rectangular, and the smallest is the semicircular arch. However, the attenuation trend of overpressure propagation is the same for all three cross-sectional shapes.

In the three-dimensional physical models of the above three roadway cross-sectional shapes, the size of the coal mining face is always consistent with the actual situation. Under the conditions of different cross-sectional shapes of the inlet/return air roadway, the overpressure changes at various working face position, which is shown in Figure 14.

Figure 14 shows that under the condition of different cross-sectional shapes of air inlet/return roadways, the initial overpressure of the working face cut is largest when the cross-sectional shape of the roadway is rectangular, followed by the semicircular arch and the trapezoid. With the propagation of the shock wave along the working face cut, the farther away from the upper corner of different cross-sectional shapes, the closer the overpressure at each point of the cut. The pressure attenuation trend is in line with the power function law. The overpressure attenuation model of the upper-corner gas explosion working face is obtained as follows: .

5. Conclusions

(1)Through an analysis of the attenuation coefficient and overpressure difference of the adjacent measuring points, it is concluded that the blast wave of a gas explosion propagates stably in the form of a single shock wave at 40 m away from the explosion source, and the effective length of the inlet/return air roadway in the upper corner with U-type ventilation is determined to be 40 m.(2)The simulation results of a gas explosion in the upper corner of three different cross-sectional shapes were analyzed. When the cross-sectional shape is a trapezoid, the initial overpressure of the explosion shock wave propagating to the inlet/return air roadway is the largest, followed by the rectangle and semicircle arch, but the attenuation trend of overpressure propagation in the three different cross-sectional shapes is the same.(3)When the gas explosion shock wave propagates along the cut hole of the working face, the farther the distance from the upper-corner angle, the closer the overpressure at each point of the cut hole. The decay trend of the overpressure conforms to the law of power function. The attenuation model is .

Data Availability

The data used to support the findings of this study are included within the article.

Conflicts of Interest

The authors declare that they have no conflicts of interest.

Acknowledgments

This work was supported by National Natural Science Foundation of China (52074106 and 51734007). This project was also supported by Zhongyuan Postdoctoral Innovative Talent (ZYQR201810171), China Postdoctoral Science Foundation (2019M652536), and Henan Postdoctoral Foundation (001801016) and funded by the Doctoral Fund of Henan Polytechnic University (B2018-59).