Abstract

The raising and flowing of groundwater caused by coal mining threaten the stability of mining faces, which cause casualties and machine damage accidents. Among the above accidents, the water inrush disaster caused by the water-rich water-conducting fault zone is the largest. Considering the complexity of geological structure and the suddenness of water inrush, reserving a reasonable thickness of waterproof coal pillars in front of the fault tectonic belt can effectively predict and control the occurrence of water inrush. The excellent adaptability of the numerical model to the geological conditions makes it an effective research method for simulating waterproof coal pillars. Based on the analysis of the background of on-site mining, this paper proposes a three-zone waterproof coal pillar calculation theory and establishes a numerical model for comparative analysis. The comparison results show that (1) the elastic-plastic theory and fracture theory can be used to calculate the thickness of the disturbed zone and the water-resisting zone, and the thickness of the fractured zone is positively correlated with the accuracy of the existing detection technology and equipment. (2) For the numerical model results, the increase of tangential stress is positively correlated with the distance of coal seam mining and the thickness of fault; the large plastic zone of the fault causes a higher increase in pore pressure, which ultimately increases the risk of water inrush. (3) The two results are in good agreement. The theoretical results have a safety margin, indicating that the three-zone theory is reasonable, which are used to guide the actual mining of the project to ensure the smooth passage of the project through the fault area.

1. Introduction

When underground coal seam mining encounters faults with complex structure, low strength, and instability under external interference, it is easy to cause collapse, roof caving, and water seepage accidents [15]. In particular, when mining on a confined aquifer, high water pressure and stress redistribution activate existing faults and further improve the permeability of faults; thus, the water inrush channel is formed by a groundwater conduction fault, in which a significant amount of confined water can suddenly invade the mining face [68]. The risk and casualty of water inrush disasters have increased due to the existence of faults [9, 10], and more than 10000 deaths have occurred due to water inrush accidents worldwide [11, 12]. Therefore, it is of great theoretical significance to study the fault activation and water conduction mechanism to prevent the occurrence of water inrush accidents.

It is widely recognised that both water pressure and mining activities have considerable influences on fault activation [1315]. The stress redistribution resulting from the underground mining process exerts an enormous influence on the stress state of nearby rock masses, which could lead to a decrease in the principal stress on regional faults [16, 17]. Meanwhile, the original crack propagation of the rock mass in a fault further reduces the stress intensity of the rock mass and becomes a new path for hydraulic migration [1820]. Furthermore, the development of the failure zone formed by coal seam mining breaks over the hydraulic crack zone of the fault, which indicates that a water inrush channel is formed [2123]. In addition, mining height also has an important influence on water inrush in underground mines where faults are present. Therefore, many numerical analyses of mining models, including fault-water geological structures, have been employed by scholars to solve the problem, for instance, finite element and discrete element methods [8, 2429]. Among these studies, the conclusions about fault water inrush obtained by the FLAC3D software are important for understanding the mechanical properties of faults [12, 3033]. Furthermore, few theories consider the water pressure, as well as the waterproof coal pillar in the process of coal seam excavation. In order to better understand the problem of fault water inrush encountered in underground mining, it is necessary to incorporate the mechanism of fault action into the simulation of coal mining.

Based on the field monitoring data and geological conditions, this paper aims to quantitatively evaluate the influence of fault size and water content on the width of waterproof coal pillar in coal seam mining activities by combining numerical model analysis with theory. The width acquired within the theoretical approach is then quantitatively compared to that of the FLAC3D model and other approaches estimated from mining engineering regulations, which help to gain useful insight into the mechanical properties of the faults and to provide knowledge about the numerical simulation of waterproof coal pillar retaining.

2. Engineering Background

The Wugou Mine is located 50 km southwest of Huaibei city, Anhui Province, China, covering a minable area of approximately 15 km2. There are multiple coal seams in the mining area, and the thickness of the #10 coal seam is 3.95 m with a burial depth of 450 m. A large compressive normal fault F14 tends to the NW, with an inclination angle of 70°, a drop between 10 and 380 m, and an extension length greater than 4300 m. The roof of the #10 coal seam is a medium and silty sand composite stratum, the floor is sandstone and mudstone, and the limestone of the Taiyuan Formation is an aquifer. The Wugou coal mine is mainly affected by the fourth aquifer of the Cenozoic, and the water pressure at the bottom of the #10 coal seam is 2.6 MPa. During the mining of the #1011 working face, the water inflow in the goaf is expected to be 35 m3/h. Thus, it is necessary to evaluate the degree of confined water inrush near-fault F14 and consider a reasonable coal pillar height. The geological structure map and geological profile of the Wugou mine field are shown in Figure 1, and Table 1 shows the physical and mechanical properties of roof and floor strata in coal seams.

3. Theoretical Content Results

3.1. Model of Three-Zone Calculation Theory

The strata between the mining face and the fault are divided into three parts: the influence area of construction disturbance, the main water-resisting area, and the water pressure fracture area. The mechanical model diagram is shown in Figure 2.

3.2. Construction Disturbance Influence Area

The failure zone of the rock mass at the edge of the excavation face under the plane stress state is larger than that of the plane strain state, and the failure domain under the plane strain state considering stress relaxation is close to that of the plane stress state. Therefore, in practical applications, the range of the failure zone of the surrounding rock edge of the mining face is calculated according to the plane strain of the stress relaxation. The mining section is circular, and the rock mass is assumed to be isotropic, continuous, and homogeneous. Based on the two-way compression plate of elastic theory, the plane strain state at a point on the centre line of the mining face is analysed, as is shown in Figure 3.

3.2.1. Stress State Analysis

According to Chelsea’s answer in elasticity theory, we get the following equation:where and are the radial and circumferential normal stresses, respectively; is the shearing stress; and are the maximum and minimum principal stresses, respectively; a is the excavation radius; and r and θ are the distance and included angle from point B to the centre, respectively. Set the lateral pressure coefficient of the surrounding rock as follows:where each parameter has the same meaning as in equation (1).

At this point, the elastic stresses near the circular excavation face are shown in Figure 4.

As shown in Figure 4, the principal stress difference around the excavation face is the largest, so is the shear stress. The deformation and failure begin to develop gradually from the periphery of the excavation face to the depth of the surrounding rock. The thickness of the plastic zone is used as the thickness of the rock layer in the fracture zone.

3.2.2. Stress Distribution of Inelastic Deformation Zone

Based on the Mohr–Coulomb failure criterion, we get the following equation:where c and ϕ are the cohesion and friction angle of the rock stratum, is the support reaction force, and the other parameters have the same meaning.

Plastic Zone of Circular Mining Face. H. Kastner proposed that the plastic failure zone of roadway excavation is a circular plastic zone around the roadway, and its calculation formula is as follows:where the meaning of each parameter in the formula is the same as in the above equations.

3.3. Main Water-Resisting Zone

Affected by mining activities, the stress release of the main water-resisting strata reduces the ability of the rock to resist the splitting of the surrounding fissure water, resulting in the growth and even penetration of the original cracks in the strata, and the splitting failure causes water inrush.

3.3.1. Initiation and Propagation of Fracture

When brittle failure occurs in closed cracks, the stress concentration at the crack tip leads to crack propagation and failure, which is shown as follows:where is the tensile stress and the other parameters have the same meaning.

The compressive stress required for fracture closure is very small, which can be negligible, and the equation can be simplified as follows:where f is the coefficient of friction and the other parameters have the same meaning.

3.3.2. Thickness of the Hydraulic Fracturing Zone

When the mining is advancing to a certain distance, the critical water pressure of the hydraulic fracturing of strata is reached, and fracture failure of the face is about to occur.

According to the D-B model of the plastic zone, the size of the plastic zone is determined by the stress state at the crack tip. When the stress intensity factor at the crack tip is zero, the size of the plastic zone is obtained as follows:where is the radius of the plastic zone, is the initial stress, and a is the excavation radius.

Under the small-scale yield, we get the following equation:

Expand the formula in a Taylor series as follows:where the parameters have the same meaning as in equation (7).

The first item was choosen only as follows:where K is the stress intensity factor and the other parameters have the same meaning as in equation (7).

The same application of the above formula was performed under large yield conditions. The stress state of the surrounding rock crack during hydraulic fracturing is as follows:where and are the normal stress and shear stress of rock under water pressure, respectively; and are the maximum and minimum principal stresses, respectively; β is the friction angle of rock under water; and is the critical water pressure under tension or compression-shear stress, which is expressed as follows:where the parameters have the same meaning as before. It can be obtained as follows:where the parameters have the same meaning as in equation (12).

3.4. Rock Thickness in Fracture Zone

Due to the complex nature of rock fractures around the fault, drilling or geophysical prospecting methods are used to obtain a more detailed evaluation of groundwater erosion of rock strata during the excavation of rock strata. The hydraulic characteristics of rock strata in front of mining are obtained by integrating the advanced prediction method of the geological radar. The detection results are used as the thickness value of the rock strata in the fracture area.

3.4.1. The Water-Rich Exploration of the F14 Fault

The transient electromagnetic method (TEM) is used to detect water abundance in the F14 region, and the water abundance in the detection area is shown in Figure 5. The purple virtual frame area is a water-rich area, and the length of water-rich area between fault and coal seam is 2.6 m.

3.4.2. Fault Information Revealed by Gallery and Ground Drilling

The profile results of fault drilling are shown in Figures 6 and 7 and Table 2. From Table 2, it can be seen that the length of water-rich fracture zone is 1.3 m.

3.5. Comparison of Theoretical and Experimental Results

The waterproof rock layer reserved thickness can be expressed as follows:

The values of each parameter are shown in Table 3.

The reasonable coal pillar thickness obtained from the failure theory of the three zones is 30.6 m. The width of the fault waterproof coal pillar obtained by the theoretical method is more in line with actual mining than with similar test results. However, more verification methods are needed, such as numerical simulations, theoretical comparisons, and example verification, to further determine the adaptability of the three zones theory.

4. Numerical Model Validation Method

4.1. Numerical Modelling

The FLAC-3D finite difference software was used to study the characteristics of water inrush and the reasonable safety width of the coal pillar. This simulation adopts the Mohr–Coulomb criterion of an ideal elastic-plastic rock mass. As shown in Figure 8, the size of the established 3D numerical model is 800 m (length) × 300 m (width) × 144 m (height), and the fault dip angle is 70°, which was used to simulate the length of 150 m fully mechanized 1011 mining face in Wugou Mine. The simulated actual formation conditions are shown in Table 1.

The boundary conditions of the model are the horizontal constraint around the model, the fixed constraint at the bottom, and the free boundary of the upper surface.

In this simulation, the overburden load is optimised to a uniform load, and the calculation parameters are selected in Table 4.

4.2. Results of Numerical Simulation
4.2.1. Variation Law of the Overburden Stress Field

Figure 9 shows the change in shear stress and vertical stress of the overlying rock with fault widths of 5 m, 8 m, and 12 m when the mining advances 90 m.

Figure 10 shows the change in shear stress and vertical stress of the overlying rock when the mining advances 210 m, 150 m, and 30 m.

As the advancing distance increases, the vertical stress and shear stress increase slowly. Under the same advancing distance, the shear stress of the overlying strata has a positive correlation with the increase in fault width, while the growth in vertical stress is the opposite. For example, when coal seam mining reaches a fault depth of 90 cm, the vertical stress caused by the fault fracture zone with widths of 5 m, 8 m, and 12 m is 39.791 MPa, 39.342 MPa, and 38.852 MPa, and the maximum shear stress is 8.603 MPa, 8.618 MPa, and 8.632 MPa. The changes in shear stress and vertical stress with the mining advance distance under three fault widths are shown in Figure 11.

The time and stress curves of the other two faults (5 m and 8 m) are shown in Figure 12. The increase of shear stress and normal stress is proportional to the distance between the excavation face and the fault. When the excavation distance reaches T0, the stress growth curve changes from gentle to sharp.

4.2.2. Variation of the Plastic Zone

Figure 13 shows the plastic strain range zone of the 5 m fault with the advance of the mining distance. It is shown that with the development of mining, the plastic zone is expanding accordingly. The roof and floor of the mined-out area are tensile-shear stress zones. At this distance, the roof is shear stress, and the floor is a tensile stress zone. When the distance between the mining face and fault is 30 m, the shear stress in the fault is obvious.

Figure 14 shows the plastic zones of 8 m and 12 m faults at 30 m distances under mining faces and faults. It is shown that with the same mining distance, there is much similarity in the range of the plastic strain zones in the goaf of the models with 8 m and 12 m faults, and the shear stress zone of the fault first appears in the lower plate position, which is shown in the red line section. The range of the plastic zone decreases and shear stress increases with increasing fault width, and at the same mining level, the footwall of the fault is in the state of compression and shear.

Figure 15 shows the displacement curves of the upper and lower plates at the same horizontal position under the width of the three fault fracture zones during the mining advance. It can be seen from Figure 15 that the vertical displacement increases step by step when mining advances to a certain position. Furthermore, the displacement of the upper plate is greater than that of the lower plate, where differences are expressed as , , and . The vertical displacement of the upper and lower plates increases with increasing fault thickness, so is the difference between the upper and lower plates, that is,  <  < . Moreover, for the upper side, the displacements of the 5 m, 8 m, and 12 m faults are 3.9 mm, 4.3 mm, and 4.8 mm, respectively.

4.2.3. Pore Water Pressure

Figure 16 shows the pore pressure of the 5 m fault under different propulsion distances. It is seen from Figure 16 that with the advance of the mining face to the fault, the pore pressure value increases continuously, which increases from 0.8 MPa to 1.7 MPa, and the maximum value appears at both ends of the goaf. When the mining advances to 30 m from the fault, the pore pressure in the fault increases to 0.4 MPa, which means water pressure is high, and seepage is considered to have occurred.

Figure 17 shows the pore pressure change cloud chart of the 8 m and 12 m faults for a 30 m distance between the mining face and fault. When the mining advances to 30 meters from the fault, with the increase in fault width, the water pressure in the fault is increasing, and the greater the growth range of water pressure is, the greater the possibility of water inrush. The fault becomes an important water inrush channel, which leads to underground aquifers and the mining face.

Figure 18 shows the pore pressure and flow velocity of the fault under different mining advances. It can be seen from Figure 18 that the pore water pressure increases nonlinearly with mining distance and fault width, where peak value is 0.3 MPa when the mining distance is 30 m from the fault. Taking the inflection point of three pore pressure curves as the reference points , , and for determining the width of waterproof coal pillar can prevent further increase of pore pressure in fault F14.

4.2.4. Fault Seepage Characteristics

When the width of the fault is 12 m, the flow vector along the mining advancing fault is shown in Figure 19. The flow vector nephogram of the 5 m and 8 m fault fracture zones when the mining distance is 30 m is shown in Figure 20.

In Figure 19, when the distances between the working face and the fault are 210 m, 150 m, 90 m, and 30 m, the flow vector is 7.207 L, 9.637 L, 13.84 L, and 20.58 L, respectively. With the continuous advancement of the working face, the width of the corresponding coal pillar decreases. The plastic damage degree of the fault and the surrounding rock of the fault caused by mining are increasing, and the water pressure of the aquifer is also increasing. Finally, with the gradual decrease in the coal pillar width, the flow rate of water is gradually accelerated. Figure 20 shows that when the working face is mined to 30 m away from the fault, the flow rate of water increases due to the influence of the width of the fault fracture zone. This phenomenon proves that there is a positive correlation between the risk of water inrush and the width of the fault fracture zone; that is, the risk of water inrush increases with the increase in the width of the fault fracture zone, and with the increase in the width of the fault fracture zone, the width of the safety pillar required for mining also increases accordingly.

5. Discussion

The recommended height of waterproof coal pillar in the #10 coal seam of Wugou coal mine obtained from three-zone theory and numerical model is shown in Table 5.

For the thickness of waterproof coal pillar, the theoretical calculation results are greater than the numerical model results, and the width of waterproof coal pillar increases with the increase of fault width. Therefore, for engineering applications, the theoretical calculation results of the three zones are relatively safe, and the theoretical results of the three zones can be applied to engineering applications.

6. Conclusions

During coal seam mining, groundwater enters the excavation face through faults, and the water inrush disaster caused by the tunnel face is the main problem faced by mining near confined water. The engineering practice shows that the reasonable thickness of waterproof coal pillar in front of the tunnel face is an effective method to solve this problem. Based on the mining background of Wugou coal mine, this paper establishes the mechanical model of waterproof coal pillar, analyses the mechanical properties and water conductivity of faults under different excavation distances and fault thicknesses, and obtains the calculation model of three-zone waterproof thickness. Finally, the rationality of the calculation results is verified by the numerical model results. The main conclusions are as follows:(1)Comparing the three-zone theory and numerical model, it is concluded that the width of waterproof coal pillar is positively correlated with the influence depth of fault and water pressure. When mining advances, the tangential stress and normal stress increase. In front of the working face, the normal stress is tensile stress, and the compressive stress is at the roof and floor. However, the peak value of tangential stress mainly appears at both ends of the excavation face, and it shows compression-shear failure at the initial section and tensile-shear failure at the mining end.(2)The pore water pressure and flow rate increase nonlinearly with the increasing mining distance and fault thickness, with peak values of 0.8 MPa and 0.35 L/s for a fault width of 12 m. The peak value of the upper and lower plates of the fault increases sharply, which means that fault activation occurs and a water rush channel is formed.(3)The comparison of theoretical and numerical simulations indicates that when the mining distance approaches the fault, a reasonable waterproof coal pillar can effectively control the occurrence of water inrush, and the numerical simulation results of coal pillar width are larger than theoretical values in three regions. In this study, the recommended thickness of the waterproof coal pillar in the No. 1011 working face of the Wugou coal mine is 80.2 m.(4)The three-zone theoretical model in this paper provides a reasonable solution for the study of waterproof coal pillars, which is helpful to understand the hydraulic change process of the rock strata on the working face during near-fault mining, and can be used as an auxiliary measure to prevent and control water disasters on the excavation face. However, further research is needed on the water inrush mechanism of the fault tectonic zone and the countermeasures for the prevention and control of water on the working face.

Data Availability

The datasets generated and analysed during the current study are not publicly available due to continuous studying but are available from the corresponding author on reasonable request.

Conflicts of Interest

The authors declare that they have no conflicts of interest.

Authors’ Contributions

Y.Q.W. wrote the main manuscript text, B.S.Y. guide the theme and theory of the manuscript, and J.Z. prepared Figures 814. All authors reviewed the manuscript.

Acknowledgments

This research was supported by the Qingdao Huanghai University Doctoral Research Foundation. The authors thank AJE (https://www.aje.com/) for English language editing and review services.