Characteristics of Airﬂow Reversal of Excavation Roadway after a Coal and Gas Outburst Accident

: Determining the inﬂuence scope of the airﬂow disorder is an important problem after coal and gas outburst accidents in ventilation systems. This paper puts forward the indexes of airﬂow disorder, including the length of the excavation roadway, the outburst pressure, the pressure difference of the air door, and the air quantity of the auxiliary fan. Using the orthogonal table of L 9 (3 4 ) and numerical simulation method, the characteristics of airﬂow reversal are studied, and the outburst airﬂow reversal degree is calculated should the ventilation facility fail. Furthermore, on the basis of fuzzy comprehensive optimization theory, the comprehensive evaluation model of the airﬂow disorder is established. The results show that the length of the excavation roadway is the most important factor affecting the stability of the ventilation system, followed by the outburst pressure, pressure difference of the air door, and air quantity of the auxiliary fan. The inﬂuence of a gas outburst accident on the return air system is greater than that on the inlet air system, and a larger air velocity has a greater impact on the ventilation system, especially the air inlet part. Moreover, the airﬂow reversal degree of the ventilation system increases with the increase of the outburst pressure or decreases with the length of the excavation roadway. This paper provides a basis for the prevention of gas outburst accidents. On the basis of the fuzzy comprehensive optimization theory, the evaluation model of the airﬂow disorder is established. The results show that the order of effect factors for the stability of the mine ventilation system are the length of excavation roadway > outburst pressure > pressure difference of air door > air quantity of auxiliary fan. The most inﬂuential result is that the length of the excavation roadway is 200 m, the outburst pressure is 0.3 MPa, the pressure difference of the air door is 1500 Pa, and the air quantity of the auxiliary fan is 600 m 3 /min. The airﬂow reversal degree of the ventilation system increases with the increase of the outburst pressure or decreases with the length of the excavation roadway.


Introduction
A coal and gas outburst is an accident that is caused by a large amount of coal and gas being ejected into the underground roadway in a single moment [1]. In an outburst accident, the high-energy shock wave destroys the roadway facilities instantly and changes the roadway resistance and the structure of the ventilation network [2]. After the dynamic effect of the shock wave disappears, the gas continues to flow and diffuse unsteadily [3], generating ventilation pressure and additional force [4], causing airflow disorder [5], casualties, or secondary gas explosion accidents.
A coal and gas outburst is studied from three aspects: shock wave [6][7][8], outburst gas [9][10][11][12], and outburst coal rock [13,14]. The theoretical analysis [15] and numerical simulation are popular methods [16,17] in a single roadway. After the outburst power disappears in the upwind roadway, Feng et al. [18] concluded, the influencing factors of the airflow reversal in the parallel branch include the pressure of the main fan, the height difference of roadway, the initial air velocity, length, and the sectional area of roadway. Yu et al. [19] pointed out that the airflow is affected by the seam inclination in the inlet roadway. However, the underground mine ventilation network is formed by connected roadways, and the study is often not limited to a single roadway. Zhou et al. [20] established 45 • and 135 • crossover tunnels to simulate the flow and attenuation rules of shock waves and gas in a gas outburst accident. The influence of air pressure increases with the increase of the angle between the excavation roadway and the adjacent excavation roadway, but the intersection of the roadway has no effect on the airflow reversal. They also found that the outburst pressure could not only lead to the airflow reversal in the downward ventilation roadway, but also may lead to the airflow reversal in parallel branches of the upward ventilation roadway [21,22]. The phenomenon of gas flow restriction and retention appears in the roadways with the air door and windshield, so draining the gas immediately is the key to prevent secondary disasters in the underground roadways.
In recent years, scholars have conducted a series of studies on the disaster process of the whole mine ventilation system during the outburst period. Due to the complexity of the ventilation system after an gas outburst accident, the tunneling face is usually seen as a unidirectional burst source, and a three-dimensional simulation method is adopted to simulate the flow process of air and gas in the ventilation system during the disaster period [23]. The result shows that the disturbance of the ventilation system is affected by the gas force, natural wind pressure, and ventilation power. Then, the characteristics of bidirectional gas outburst in the coal mining face are studied [24]. Different from the unidirectional outburst mode, the impact strength of this type is reduced. Ventilation facilities are the weak point of the ventilation system. The research mainly focuses on the characteristics of air door failure [25,26], the reasonable position, the strength of the air door [27], as well as the automatic air door [28,29]. There is still a lack of research on the influencing factors of airflow disorders in ventilation systems when ventilation facilities fail. There are many factors that affect the ventilation system of a coal and gas outburst [30][31][32], it is particularly important to establish the index system of the airflow disorder and determine relatively important influencing factors when ventilation facilities fail, which could provide theoretical support for disaster prevention.
The whole paper is structured as follows: Section 2 introduces the ventilation system. In Section 3, the influencing factors of the airflow disorder induced by gas outburst accident are put forward. Using the orthogonal experiment analysis and numerical simulation method, the sensitivity factors of the airflow disorder when the ventilation facility fails are analyzed in Section 4. Section 5 conducts a comprehensive evaluation of the airflow disorder on the basis of fuzzy comprehensive optimization theory.

Ventilation System and Its Simplification
Most of the coal and gas outburst accidents occur in the heading face [24]. The ventilation system is formed by the excavation roadway, heading face, auxiliary fan, and air ducts. Figure 1 illustrates a forced ventilation system, which is the most commonly used ventilation form. It has the advantages of a small air leakage and large air supply. The fresh air is transported from the intake airway to the heading face through the auxiliary fan and air duct, and the polluted air enters the return airway through the excavation roadway. The air door obstructs fresh air and polluted air, and the air duct needs to pass through the air door. The blue arrow indicates the airflow direction of the fresh air, and the red one indicates the polluted air.
In general operation conditions, the order of the air gauge pressure is P 20 > P 21 > P 22 , P 10 > P 11 > P 12 > P 13 , P 11 > P 31 > P 21 .
According to graph theory, the air duct can be considered as a separate roadway for the air supply on the heading face, so it is an independent branch. There are nine branches and nine nodes in a general ventilation system. The ventilation system is simplified, as shown in Figure 2.
In a gas outburst accident, the air pressure of the roadway is redistributed as the gas concentration increases rapidly in the ventilation system, and the shock wave usually destroys the ventilation facilities. As a result, if the air door is invalid in a tunneling ventilation system, the outburst gas will flow into the intake airway. Then, the airflow will be reversed if the air pressure at the upwind side is lower than the outburst pressure.  In general operation conditions, the order of the air gauge pressure is P20 > P21 > P22 P10 > P11 > P12 > P13, P11 > P31 > P21.
According to graph theory, the air duct can be considered as a separate roadway for the air supply on the heading face, so it is an independent branch. There are nine branches and nine nodes in a general ventilation system. The ventilation system is simplified, as shown in Figure 2.  In a gas outburst accident, the air pressure of the roadway is redistributed as the gas concentration increases rapidly in the ventilation system, and the shock wave usually destroys the ventilation facilities. As a result, if the air door is invalid in a tunneling ventilation system, the outburst gas will flow into the intake airway. Then, the airflow wil be reversed if the air pressure at the upwind side is lower than the outburst pressure.
The outburst airflow reversal degree is defined as the influence of outburst gas on the airflow disorder in an airway, which is the ratio of the air quantity variation after the gas outburst accident and the original value. A positive reversal degree means that the airflow direction of the outburst gas is the same as the original airflow, and a negative reversal degree means that the airflow is reversed.

Indexes of Airflow Disorder
There are 3 factors that affect the safety and stability of the mine ventilation system as shown in Figure 3.    In general operation conditions, the order of the air gauge pressure is P20 > P21 > P22, P10 > P11 > P12 > P13, P11 > P31 > P21.
According to graph theory, the air duct can be considered as a separate roadway for the air supply on the heading face, so it is an independent branch. There are nine branches and nine nodes in a general ventilation system. The ventilation system is simplified, as shown in Figure 2.  In a gas outburst accident, the air pressure of the roadway is redistributed as the gas concentration increases rapidly in the ventilation system, and the shock wave usually destroys the ventilation facilities. As a result, if the air door is invalid in a tunneling ventilation system, the outburst gas will flow into the intake airway. Then, the airflow will be reversed if the air pressure at the upwind side is lower than the outburst pressure.
The outburst airflow reversal degree is defined as the influence of outburst gas on the airflow disorder in an airway, which is the ratio of the air quantity variation after the gas outburst accident and the original value. A positive reversal degree means that the airflow direction of the outburst gas is the same as the original airflow, and a negative reversal degree means that the airflow is reversed.

Indexes of Airflow Disorder
There are 3 factors that affect the safety and stability of the mine ventilation system, as shown in Figure 3. The outburst airflow reversal degree is defined as the influence of outburst gas on the airflow disorder in an airway, which is the ratio of the air quantity variation after the gas outburst accident and the original value. A positive reversal degree means that the airflow direction of the outburst gas is the same as the original airflow, and a negative reversal degree means that the airflow is reversed.

Indexes of Airflow Disorder
There are 3 factors that affect the safety and stability of the mine ventilation system, as shown in Figure 3.

Parameters of Outburst Accident
The outburst intensity indicates the amount of coal and rock thrown into the roadway and the amount of gas emissions in a coal and gas outburst accident. Because it is difficult to calculate the amount of gas emission, the outburst intensity is represented by the

Parameters of Outburst Accident
The outburst intensity indicates the amount of coal and rock thrown into the roadway and the amount of gas emissions in a coal and gas outburst accident. Because it is difficult to calculate the amount of gas emission, the outburst intensity is represented by the amount of outburst coal and rock. Usually, an outburst accident is divided into a small outburst accident (<100 t), medium outburst accident (<500 t), large outburst accident (<1000 t), and extra-large outburst accident (≥1000 t).
The outburst position refers to the place where the outburst accident occurs. If the air pressure of the outburst position is low, the outburst gas can be faster exhausted to the ground, as the outburst position is closer to the air shaft. On the contrary, it has a great impact on the ventilation system, especially the air intake roadway. Therefore, we can reduce the impact of the outburst accident on the ventilation system by reducing the air resistance of the return airway or increasing the intake airway.

Physical Conditions of Roadway
The characteristics of the roadway are determined by the physical conditions, including the section and length of roadway. Specifically, the outburst gas is easy to exhaust from the large section or the shorter length of roadway.

Safety and Stability of Mine Ventilation System
The high-pressure outburst gas suddenly flows into the roadway space in a gas outburst accident, which has a great impact on the mine ventilation system. It is easy to cause a large range of airflow disorders in the mine ventilation system after the shock wave destroys the air doors.
Among the above factors, the outburst pressure, length of roadway, pressure difference of the air door, and supply air quantity of the roadway are chosen to analyze the laws of the airflow disorder and the degree of influence in a mine ventilation system.

Orthogonal Design
The orthogonal experiment is used for the analysis of the mutual influence of multiple factors [33]. Using the orthogonal table, an appropriate amount of representative points is selected from a large number of test points. Through the overall design, comprehensive comparison, and statistical analysis, the balanced sampling within the range of factor change is carried out. Then, the analysis of the representative results is achieved by a small number of tests.
In the orthogonal design, the table of L 9 (3 4 ) is adopted to study the effect of 4 factors, i.e., A is the length of roadway, B is the outburst pressure, C is the pressure difference of air door, and D is the air quantity of the auxiliary fan on the airflow disorder in the mine ventilation system. The specific parameters of each test are illustrated in Table 1.  Table 1, the level of each factor is in the common value range, such that the length of roadway is 200 m, 1000 m, and 2000 m, the outburst pressure is 0.1 MPa, 0.2 MPa, and 0.3 MPa respectively, the pressure difference of the air door is 500 Pa, 1000 Pa, and 1500 Pa, the air quantity of the auxiliary fan is 240 m 3 /min, 600 m 3 /min, and 1200 m 3 /min. A total of nine combined schemes are designed.

Physical Model and Meshing
According to the orthogonal tests, it is necessary to establish three physical models with different roadway lengths. The spatial relationship model and grid meshing of the ventilation system are shown in Figure 4. The parameters of each test model are shown in Table 2.  Table 1, the level of each factor is in the common value range, such tha length of roadway is 200 m, 1000 m, and 2000 m, the outburst pressure is 0.1 MP MPa, and 0.3 MPa respectively, the pressure difference of the air door is 500 Pa, 100 and 1500 Pa, the air quantity of the auxiliary fan is 240 m 3 /min, 600 m 3 /min, and m 3 /min. A total of nine combined schemes are designed.

Physical Model and Meshing
According to the orthogonal tests, it is necessary to establish three physical m with different roadway lengths. The spatial relationship model and grid meshing o ventilation system are shown in Figure 4. The parameters of each test model are show Table 2.    There are three basic equations of the fluid numerical simulation, which includes the mass conservation equation, the momentum conservation equation, and the energy conservation equation. After the gas outburst accident in the excavation roadway, the airflow should also be satisfied by gas migration. Furthermore, a chemical component exchange exists in this numerical simulation, i.e., CH 4 , O 2 , and N 2 , so the component transport equation should be introduced.
A turbulent flow is a common flow phenomenon in nature, and the fluid is always in the state of turbulence in most engineering problems. The realizable k-ε model is adopted because its good performance for solving the adverse pressure gradient and vortex problems.

Boundary Conditions
The boundary conditions of the model are shown in Table 3. From Table 3, the physical model area is the Fluid Type, and the pressure of operation condition is 101.325 kPa. The inlet boundary condition includes Inlet 1 of the intake airway, Inlet 2 of the return airway, outburst position, and auxiliary fan. The inlet boundary type of auxiliary fan inlet is the Velocity-inlet, and the other is the Pressure-inlet. The species of the outburst position is 100% CH 4 , and the other is 20% oxygen and 80% nitrogen. The outlet boundary includes Outlet 1 of the intake airway and Outlet 2 of the return airway, the outlet boundary type is also the Pressure-outlet. Furthermore, the species are 20% oxygen and 80% nitrogen either.
The gauge pressure of each boundary condition is calculated as follows: Assuming that both the sectional area of the return airway and excavation roadway is 25 m 2 , the friction resistance coefficient is 100 × 10 −4 , the inlet air quantity of the intake airway and return airway Q 10 = Q 20 is 75 m 3 /s, and the air quantity of the auxiliary fan Q fan is 4 m 3 /s, 10 m 3 /s, and 20 m 3 /s, respectively. The outlet of the return airway is Q 21 = Q 20 + Q fan , Q 12 = Q 10 − Q fan . The relationship between the air resistance and the friction resistance coefficient of the roadway is where R represents the air resistance, kg/m 7 ; L represents the length of the roadway, m; U represents the perimeter of the roadway section, m; α represents coefficient of frictional resistance, kg/m 3 ; S represents the sectional area of the roadway, m 2 . From Equation (1), the air resistance value of the return airway is R 10 = R 12 = R 20 = R 21 = 0.0032 kg/m 7 , the excavation roadway from the air duct outlet to the return airway is R 31 = 0.0256 kg/m 7 . The air pressure H 10 = 16.13 Pa, H 12 = H 20 = 18 Pa, H 21 = 19.97 Pa and the air pressure of the excavation roadway are calculated by the air quantity of the auxiliary fan and the length of the excavation roadway.
Furthermore, the initial gauge total pressure at Inlet 1 is the sum of the pressure difference of the air door and the air pressure of the upper wind side of the intake airway H 10 , and the initial pressure of Outlet 1 is the sum of the pressure difference of air door and H 12 . Suppose the operation pressure is the standard atmospheric pressure (101.325 kPa), the gauge total pressure at Outlet 2 of the return airway is 0 Pa. Then, the air pressure of each inlet and outlet are calculated as shown in Table 3.

Controlling Parameters of Fluid Dynamics
(1) Solver In the Pressure-based type, the absolute velocity formulation and Pressure-velocity coupling algorithm are usually adopted to solve the problem. Furthermore, using the SIMPLE scheme to calculate the mathematical model, it computes the mass conservation and obtains the pressure field by the mutual correction of pressure and velocity.
(2) Convergence accuracy A reasonable accuracy is an important parameter to ensure the convergence of the model. The convergence residuals are 10 −6 in this model, and the number of iterations is 500, so the calculation will be finished when the residuals of each variable are less than Energies 2021, 14, 3645 7 of 13 10 −6 or the iterations reach to 500 steps. Furthermore, the convergence of the calculation results can be dynamically monitored by checking the iterative residual of each variable.

Results of Numerical Simulation
Using Fluent numerical simulation software, the flow field distributions of the outburst gas of the above nine tests are simulated. Figure 5 illustrates the path-line near the air door of each test on the horizontal plane 1 m away from the floor, which is colored by the velocity magnitude. SIMPLE scheme to calculate the mathematical model, it computes the mass conservation and obtains the pressure field by the mutual correction of pressure and velocity.
(2) Convergence accuracy A reasonable accuracy is an important parameter to ensure the convergence of the model. The convergence residuals are 10 −6 in this model, and the number of iterations is 500, so the calculation will be finished when the residuals of each variable are less than 10 −6 or the iterations reach to 500 steps. Furthermore, the convergence of the calculation results can be dynamically monitored by checking the iterative residual of each variable.

Results of Numerical Simulation
Using Fluent numerical simulation software, the flow field distributions of the outburst gas of the above nine tests are simulated. Figure 5 illustrates the path-line near the air door of each test on the horizontal plane 1 m away from the floor, which is colored by the velocity magnitude. From Figure 5, the air velocity of L3 is the highest near the air door after the gas burst accident, followed by L2 and L1, and L7 is the smallest. The excavation roadway length of L1, L2, and L3 is 200 m, which has the greatest air velocity magnitude. The larger air From Figure 5, the air velocity of L 3 is the highest near the air door after the gas burst accident, followed by L 2 and L 1 , and L 7 is the smallest. The excavation roadway length of L 1 , L 2 , and L 3 is 200 m, which has the greatest air velocity magnitude. The larger air velocity has a greater impact on the ventilation system, especially in the air inlet part. Figures 6 and 7 show the air velocity distributions of the two inlet boundaries. velocity has a greater impact on the ventilation system, especially in the air inlet part. Figures 6 and 7 show the air velocity distributions of the two inlet boundaries.   As can be seen from the contour map in Figures 6 and 7, all inlet boundaries are reversed in each test, and the air velocity of Inlet 2 is larger than Inlet 1. Therefore, the influence of the gas outburst accident on the return air system is greater than that on the inlet air system. Extracting the maximum air velocity of the section, the airflow reversal degree is calculated, as shown in Table 4. Meanwhile, the variation curve of air velocity and airflow reversal degree is shown in Figure 8.  As can be seen from the contour map in Figures 6 and 7, all inlet boundaries are reversed in each test, and the air velocity of Inlet 2 is larger than Inlet 1. Therefore, the influence of the gas outburst accident on the return air system is greater than that on the inlet air system. Extracting the maximum air velocity of the section, the airflow reversal degree is calculated, as shown in Table 4. Meanwhile, the variation curve of air velocity and airflow reversal degree is shown in Figure 8.  From Table 4 and Figure 8, the airflow reversal degree of Inlet 2 is larger than Inlet 1, that is, because the air pressure distribution of the ventilation system can restrain the outburst gas, and the air pressure variation of the intake airway is smaller than the return airway.  Table 4 and Figure 8, the airflow reversal degree of Inlet 2 is larger than Inlet 1, that is, because the air pressure distribution of the ventilation system can restrain the outburst gas, and the air pressure variation of the intake airway is smaller than the return airway.
The numerical simulation results show that the airflow reversal degree of 200 m roadway is greater than that of 1000 m roadway and 2000 m roadway. Therefore, the airflow reversal degree decreases with the increase of the roadway length. Furthermore, the outburst air flow pressure of L 3 , L 6 , and L 9 is 0.3 MPa, and the absolute value of their airflow reversal degree is greater than the outburst pressure of 0.1 MPa and 0.2 MPa in the same length of the excavation roadway. Therefore, the airflow reversal degree increases with the increase of the outburst pressure and the outburst energy.

Fuzzy Comprehensive Optimization Theory
According to the main factors that affect the safety and stability of the ventilation system as well as the basic data of orthogonal experiments, the air quantity of the auxiliary fan, the air pressure difference of the air door, and the air flow reversal degree of Inlet 1 and Inlet 2 are used to establish the evaluation model of the airflow disorder on the basis of variable fuzzy theory.
Assuming F = {F ij }, i = 1,. . . . , m, j = 1,. . . . , n indicates that the parameter of the m test model corresponds to the set of evaluation factors of n. Here, m = 4. The method of relative membership is adopted to deal with the initial data. There are three situations, (a) The larger the better where Fr ij is the relative membership of the set F ij and its range is (c) The value is equal to 1 The fuzzy partition matrix is defined as follows, where u 1i is subordinate to the superior index and u 2i is the inferior index, The objective function is the minimal sum of the weight difference squared of the m models.
where ω j is the weight of each factor; gr j represents the standard superior membership, which is the optimal value of evaluation factor j; br j represents the standard inferior subordinate, which is the worst value of evaluation factor j. Supposing = 0, the optimal fuzzy partition matrix is calculated. Then, the influence of different factors is evaluated by u 1i .

Factors
The pressure difference of the air door has a strong inhibition of the outburst energy, which belongs to situation (a). The greater the pressure difference of the air door, the greater the pressure energy of the air inlet roadway. The air quantity of the auxiliary fan belongs to situation (a). Because it is opposite to the direction of the outburst gas, the resistance to the migration of the outburst gas increases with the larger air quantity of the auxiliary fan. The airflow reversal degree of Inlet 1 and Inlet 2 also belongs to situation (a), the larger the better.
The above factors consist of the F matrix.

Weight
Since the airflow reversal degree is more important than the pressure difference of air door and the air quantity of the auxiliary fan, the weight of the airflow reversal degree is 0.4 and 0.4, respectively, at Inlet 1 and Inlet 2. The weight of the pressure difference of the air door is equal to the air quantity of the auxiliary fan, their value is 0.1, respectively. Then the weight vector ω j = (0.1, 0.1, 0.4, 0.4) T .

Results and Discussion
The fuzzy comprehensive optimization theory is used to evaluate the airflow disorder degree of the mine ventilation network after a gas outburst accident. From Tables 1 and 4  From Figure 9, the evaluation result of L3 is the largest, illustrating that this parameter condition has the strongest influence on the ventilation system, and L7 is the smallest. In L3, the length of excavation roadway is 200 m, the outburst pressure is 0.3 MPa, the pressure difference of air door is 1500 Pa, and the air quantity of the auxiliary fan is 600 m 3 /min. In L7, the length of the excavation roadway is 2000 m, the outburst pressure is 0.1 MPa, the pressure difference of air door is 1500 Pa, and the air quantity of the auxiliary fan is 1200 m 3 /min.
According to the extreme difference of the range analysis method, the results of the orthogonal experiment are determined. The extreme difference is the difference between the maximum and minimum value of the airflow disorder of each factor. With the increase of the influence range, the disturbance degree of the factor increases. The result is shown From Figure 9, the evaluation result of L 3 is the largest, illustrating that this parameter condition has the strongest influence on the ventilation system, and L 7 is the smallest. In L 3 , the length of excavation roadway is 200 m, the outburst pressure is 0.3 MPa, the pressure difference of air door is 1500 Pa, and the air quantity of the auxiliary fan is 600 m 3 /min. In L 7 , the length of the excavation roadway is 2000 m, the outburst pressure is 0.1 MPa, the pressure difference of air door is 1500 Pa, and the air quantity of the auxiliary fan is 1200 m 3 /min.
According to the extreme difference of the range analysis method, the results of the orthogonal experiment are determined. The extreme difference is the difference between the maximum and minimum value of the airflow disorder of each factor. With the increase of the influence range, the disturbance degree of the factor increases. The result is shown in Table 5 and Figure 10. K i represents the sum of levels of factor i; k i represents the average value of the levels of factor; R represents the polar difference.  From Figure 10, A > B > D > C, the length of the excavation roadway and the outburst pressure are relatively important factors affecting the ventilation system. The airflow reversal degree of the ventilation system increases with the increase of the outburst pressure or the decrease of the length of the excavation roadway.
The energy is continuously attenuated while the outburst gas migrates in the roadway. Adopting the method of reducing the air resistance of the airway or changing the air pressure distribution of the main fans, the initial air pressure in both the roadway and its parallel roadway is easily affected and will be increased in a gas outburst accident. Meanwhile, it will contribute to increase the total air quantity and improve the kinetic energy of the airflow under the condition of the constant capacity of the main fans, so it is harder to reverse the airflow.

Conclusions
This paper analyzes the reversal principle of airflow induced by a coal and gas outburst in the excavation roadway and proposes the indexes of the airflow disorder. Through an orthogonal experiment, it is found that all the inlet boundaries are reversed in each test, and the air velocity of Inlet 2 is larger than Inlet 1, indicating that the influence of the gas outburst accident on the return air system is greater than that on the inlet air system. Moreover, the larger air velocity has a greater impact on the ventilation system, especially the air inlet part.
On the basis of the fuzzy comprehensive optimization theory, the evaluation model of the airflow disorder is established. The results show that the order of effect factors for the stability of the mine ventilation system are the length of excavation roadway > From Figure 10, A > B > D > C, the length of the excavation roadway and the outburst pressure are relatively important factors affecting the ventilation system. The airflow reversal degree of the ventilation system increases with the increase of the outburst pressure or the decrease of the length of the excavation roadway.
The energy is continuously attenuated while the outburst gas migrates in the roadway. Adopting the method of reducing the air resistance of the airway or changing the air pressure distribution of the main fans, the initial air pressure in both the roadway and its parallel roadway is easily affected and will be increased in a gas outburst accident. Meanwhile, it will contribute to increase the total air quantity and improve the kinetic energy of the airflow under the condition of the constant capacity of the main fans, so it is harder to reverse the airflow.

Conclusions
This paper analyzes the reversal principle of airflow induced by a coal and gas outburst in the excavation roadway and proposes the indexes of the airflow disorder. Through an orthogonal experiment, it is found that all the inlet boundaries are reversed in each test, and the air velocity of Inlet 2 is larger than Inlet 1, indicating that the influence of the gas outburst accident on the return air system is greater than that on the inlet air system. Moreover, the larger air velocity has a greater impact on the ventilation system, especially the air inlet part.
On the basis of the fuzzy comprehensive optimization theory, the evaluation model of the airflow disorder is established. The results show that the order of effect factors for the stability of the mine ventilation system are the length of excavation roadway > outburst pressure > pressure difference of air door > air quantity of auxiliary fan. The most influential result is that the length of the excavation roadway is 200 m, the outburst pressure is 0.3 MPa, the pressure difference of the air door is 1500 Pa, and the air quantity of the auxiliary fan is 600 m 3 /min. The airflow reversal degree of the ventilation system increases with the increase of the outburst pressure or decreases with the length of the excavation roadway.