Abstract
Setting up a waterproof coal pillar is an important measure to prevent water inrush from the Weibei mining through fault floor. Based on the plastic slip line field theory, a mechanical model of floor water inrush induced by confined water in the through fault zone was established. The mechanical expressions of confined water pressure and the width of the waterproof coal pillar under the state of limit equilibrium were derived. Combining the laws of floor deformation, failure and fault activation under two kinds of coal pillar width, the safety width of the waterproof coal pillar was determined. Furthermore, the safety threshold is better than the empirical value mentioned in the “coal mine safety regulations.” Following this, grouting transformation was carried out on the K2 sand layer of the cut roadway floor. This provided a theoretical basis and engineering practice for water disaster prevention and the control of the structural floor under similar conditions in the Weibei mining area for future benefit.
1. Introduction
The North China-type coalfield is rich in coal resources and accounts for about 90% of the total coal [1, 2]. However, the complex hydrogeological conditions of minefields have caused the frequent occurrence of limestone karst water hazards [2, 3]. With the large-scale development of coal resources in the recent decades, the shallow North China-type coalfield has gradually entered the stage of deep mining and lower group coal mining [4, 5]. Furthermore, the coal seam floor is increasingly threatened by Ordovician ash highly confined water [6–8].
In some old mining areas, coal reserves threatened by karst water disasters account for 49%-87.6% of total mine reserves [9]. The Carboniferous-Permian coalfield in Weibei, known as the “Weibei black belt,” is distributed in Tongchuan, Pubai, Chenghe, and Hancheng mining areas. Faults on the floor generally exist in the geological structure of the Weibei coalfield. Once mining induces fault activation, Ordovician confined water can easily rise along the fault zone and connected fractures. This results in water inrush from the floor.
Currently, remarkable theory and practice bases have been achieved by research carried out on the failure law, water inrush mechanism, and water inrush prediction [10–18]. Yin et al. [19] developed a numerical model to predict the time and the longwall locations of flood occurrences. Furthermore, Odintsev and Miletenko [20] studied the mine water inrush caused by spontaneous hydrofracturing of surrounding rocks and gave the influencing factors of spontaneous hydraulic fracturing. Additionally, Wang et al. [21] established an evaluation model for predicting water inrush in the Lu-an mining area based on fractal theory and improved the analytic hierarchy progress. This model provided a basis for water inrush research from the floor strata through faults. However, the study on the mechanism of water inrush and seepage from the through fault floor in the Weibei mining area has not yet taken a corresponding karst confined water prevention technology system and supporting measures. In the case of lower group coal resources through fault floors being blindly exploited, it is very easy for floor water inrush to proceed, endangering the safety of the construction of the mine and personnel safety.
The primary floor water inrush prevention and control measures include grouting transformation, water exploration, drainage, and setting of waterproof coal pillars [22–27]. Setting waterproof coal pillar measures is often adopted to prevent water inrush from the Weibei coal mine floor. However, the setting width of the waterproof coal pillar obtained from the empirical formula in “coal mine safety regulations” failed to consider the fault structure with complex conditions. Therefore, scientifically setting up a reasonable width of a waterproof coal pillar is one of the urgent problems to be solved through fault coal seam mining.
This paper established a mechanical model of floor water inrush induced by confined water in the through fault zone. Furthermore, this paper gave an analytical solution of the limit width of the waterproof coal pillar. By combining analysis of the floor’s deformation, failure, and the law of fault activation under two kinds of coal pillar width with engineering examples and numerical simulation, it is proven that the mechanical model is reasonable. Moreover, the grouting of floor K2 sandstone can also be used as a reference for the treatment of other similar geological structures.
2. Theoretical Analysis
2.1. Analysis of the Mechanical Model of the Plastic Slip Line
Figure 1 shows the abutment pressure distribution curve of coal in front of the working face. According to the failure degree, the coal body of the working face can be divided into the inelastic zone, the elastic zone, and the original rock stress zone under abutment pressure with the coal seam mining. Among them, the width of the inelastic zone is [28] where is the buried depth of the coal seam (m), is the mining height of the coal seam (m), is the peak coefficient of abutment pressure, is the average bulk density of the coal seam overlying rock (N/m3), is the deformation angle of the coal seam (degree), is the cohesive force on the contact surface between the coal seam and the roof and floor (MPa), is the frictional force between the interface between the coal seam and the roof, , is the ultimate compressive strength of the coal under uniaxial compression (MPa), and is the strength of the coal under uniaxial compression residual strength (MPa).

Figure 2 shows the mechanical model of the plastic slip line for mining failure of floor rock mass caused by abutment pressure in front of the coal mining face [29], where is the coal mining face position, is the goaf boundary, is the horizontal distance from the peak abutment pressure to the working face, and represents the maximum depth of mining failure of the floor. With the advancement of the working face, the maximum failure depth of floor mining is

The horizontal distance between the maximum failure depth of the floor rock mass and the front end of the work is
The maximum failure length of floor rock mass in the goaf along the horizontal plane is where is the width of the inelastic zone of the coal seam (m) and is the weighted average internal friction angle of the floor rock mass.
2.2. Mechanical Analysis of the Activated Fault Floor Structure
As the dip angle of the through fault is generally large in the Weibei mining area, this paper only discussed the fault with a large dip angle. For the convenience of derivation, the shaping slip line of the base plate is simplified to a straight line, as shown in Figure 3. Then, the angle between shaping slip line and coal seam can be expressed as

The mechanical model in Figure 4 is taken from the triangular region in Figure 3. By analyzing the force of the triangular HGM rock block above fault, the upper rock mass pressure and confined water pressure of the rock block are regarded as a uniformly distributed load, as shown in Figure 4. Among them, represents the confined water pressure, represents the uniformly distributed load of the upper rock mass acting on the rock block, and represent the shear stress and normal stress provided by the lateral rock mass, respectively. The length of each segment can be expressed as

By decomposing the forces in the diagram, the simplified force expression in the -axis direction in the limit equilibrium state is where is the weighted average cohesion of floor rock mass.
The mechanical model in Figure 5 is taken from the rectangular area in Figure 3. The force of rectangular MFGH rock above the triangle is analyzed, as shown in Figure 5. By taking any microelement in the -axis direction, the expression of force acting on the -axis in the equilibrium state is where is the stress of the element in the -axis direction (MPa).

Since the failure of the floor rock mass should conform to the principle of shear failure, the M-C criterion is selected as the yield criterion, and the limit ultimate equilibrium condition of the rock mass when it is about to fail is
Let ; then,
Substituting formulas (9) and (10) into formula (12), then solving the differential equation, and simplifying, one can obtain
When , take the unit thickness as . The following formula can be obtained: where is the average bulk density of the floor strata (kN/m3).
The and formula (14) are substituted into formula (15) and simplified as
Substituting equation (15) into the differential equation (13), we can get
By substituting into equation (16) and taking unit thickness as , the simplified expression is
Substituting equation (16) into equation (17), one can obtain
In this limit state, can be expressed as where is the width of the coal pillar corresponding to the direction of the working face (m).
By substituting equations (8) and (19) into equation (18) and simplifying it, the expression of the confined water pressure of the through fault floor and the width of the corresponding coal pillar under the limit equilibrium state can be obtained: because this equation is transcendental and cannot be simplified to the form manually. However, it can be given a value by the method of a single variable solution by software. Therefore, the width of waterproof coal pillar in the through fault floor is
3. Engineering Practice and Application
3.1. Working Face Overview
The Dongjiahe coal mine is located 3.5 kilometers southwest of Chengcheng County, Weinan City, Shaanxi Province. This mine has a high terrain in the north and low terrain in the south (Figure 6). It belongs to the Weibei Loess Plateau in geomorphology. The main coal seam of the mine is the no. 5 coal, which has a coal thickness of 2.9 m~4.7 m; the dip angle of the coal seam is 3°~15°. The 23503 working face is located in the second level and third mining area of the Dongjiahe coal mine, with a strike length of 1190 m and a dip length of 180 m. The terrain is lower than the static water level of Ordovician ash +370 m. This is classified as mining under pressure.

The 23503 working face is near the eastern boundary of the mining area. It is only about 150 m away from the normal fault of Dongjiahe, which is a big fault outside the minefield. According to drilling, the floor fault structure of the working face is well developed, and the hydrogeological conditions are complex. The distance between the 5# coal seam floor and the top face of the Ordovician limestone is about 62.82 m (CH187)~31.47 m (DB6). There is a normal fault structure on the right boundary of the working face. This fault structure has a fault dip angle of 40°, and the fault crosses the K2 limestone. If mining starts at the working face, mining disturbances would affect the normal fault and result in fracture conduction. This will then cause floor water inrush. Therefore, in the mine, we decided to reduce the disturbance to the fault by setting up a waterproof coal pillar to avoid safety accidents.
3.2. Calculation of Setting Width of the Fault Coal Pillar
3.2.1. Empirical Formula of “Coal Mine Safety Regulations”
In the “provisions on prevention and control of water in coal mine,” when considering the pressure of fault water along the direction of the coal seam, the width of the waterproof coal pillar can be calculated by referring to the following empirical formula. This calculation applies in the case whereby coal seam is located above the aquifer and the fault is water-conducting. where is the width of the coal pillar (m), is the mining height (m), is the hydrostatic pressure (MPa), is the tensile strength of coal (MPa), and is the safety factor; take 2-5.
Substitute the geological parameters of the 23503 working face in the Dongjiahe coal mine into the empirical formula (22), where the basic parameters are as follows: is 3.8 m, is 5, is 0.25 MPa, is 1.4 MPa, and the minimum width of the waterproof coal pillar is 25 m.
3.2.2. Theoretical Solution of the Setting Width of the Waterproof Coal Pillar
Based on the theoretical derivation in the previous section, the mechanical expressions of the confined water pressure with fault floor and the width of the retained waterproof coal pillar were obtained. By substituting the basic parameters of the 23503 working face in Table 1 into formulas (20) and (21), the minimum width theoretical solution of the waterproof coal pillar is 31.92 m.
In summary, the empirical and theoretical values of the waterproof coal pillar width in the 23503 working face were 25.13 m and 31.92 m, respectively. The empirical value was smaller than the theoretical value. Considering that some folds and small fault structures may exist in the actual mining of the working face, there will be unsafe factors in the empirical value of the regulations. Hence, in the mine, we initially proposed to set up a waterproof coal pillar of 28 m. To make a comprehensive study and judgment of the above results, a numerical simulation analysis was carried out on the lower cut mining at the distance of 28 m from the fault.
4. Numerical Simulation
4.1. Numerical Simulation Study on 28 m Width of the Waterproof Coal Pillar
4.1.1. Model Building
Taking the 23503 working face of the Dongjiahe coal mine as a research object and taking the simulated mining depth as 450 m, the mining height as 4.0 m, and the near-horizontal coal seam mining, the size of the model is . According to the actual lithology of typical drilled columnar strata in the 23503 working face, it was divided into 23 layers of model materials. The fault cut through the coal seam and floor aquifer and entered the Ordovician limestone aquifer. Table 2 shows some rock physical and mechanical parameters of the 23503 working face.
4.1.2. Boundary Conditions and Excavation Steps
The uniformly distributed load was applied at the top of the model to replace unsimulated overburden pressure. Water pressure of 1.4 MPa was applied at the bottom to simulate the effect of Ordovician ash confined water on the floor. The original rock stress was a static stress field, and the rock stratum was a continuous medium. As shown in Figure 7, the fault of the model was normal with a distance of 28 m from the right fault (along the coal seam) and 334 m from the left boundary. By mining from right to left, there were a total of 10 steps. For the first and second steps of mining, each step was 32 m. For the third step to the tenth step, each step was 23 m, for a total of 248 m.

4.1.3. Analysis of Simulation Results of the Plastic Zone
When the working face was mined up to 32 m (step 1-15-1), plastic failure on both sides of the goaf was almost symmetrically distributed (Figure 8). The direct roof and floor mainly experienced tensile failure, and the coal and rock masses on both sides of the goaf mainly experienced shear failure. At this time, the floor fault was almost unaffected by mining disturbance, and only local shear failure occurred inside the fault. As the working face continued advancing to 64 m (step 2-15-1), the plastic failure range of the roof and floor continued expanding and deepening. The plastic failure area was formed in the fault zone due to mining failure of the floor. Subsequently, in the process of advancing to 110 m, the plastic failure range and depth of roof and floor of goaf continued increasing, and the fault continued activating. When advancing to 110 m (step 4-11-1), the cracks between the floor mining failure zone and the pressurized water uplift zone penetrated, resulting in water inrush from the floor.

4.1.4. Analysis of Simulation Results of the Stress Field
The picture shows the variation law of the vertical stress field (Figure 9). With the continuous advancement of the working face, the stress variation range of the cutting side was larger than that of the mining face, resulting in stress concentration on the cutting side. In the process of advancing to 110 m (step 4-15-1), it was evident that the stress on the left side of the fault continued to decrease due to failure of floor mining. Plastic failure occurred when the stress range exceeded the strength of the rock mass. This provided a water inrush channel for the confined water lift.

4.1.5. Analysis of Seepage Field Simulation Results
By combining Figures 10 and 11, it can be seen that the pore water pressure in the fault increases slightly with the advancement of the working face to 64 m (step 2-15-1). Subsequently, when the working face advanced to 110 m (step 4-5-1), the fault zone was activated by disturbance and damage along the fault and confined water pressure continued to increase and rose along the fault zone. When the simulation reached step 4-10-1, the development degree of the floor crack and the seepage vector near the mining failure area of the fault zone continued to increase. Until simulation to step 4-11-1, the cracks between the floor mining failure zone and the confined water guide zone developed and led to water inrush channels. The seepage vector in the water inrush channel was huge, and the water pressure increased sharply from the previous 0.23 MPa to 0.97 MPa. Then, in the process of simulation to step 4-14-1, water pressure decreased gradually and stabilized gradually. The change of water pressure in the fault zone is shown in Figure 11.


4.1.6. Comprehensive Comparison of Simulation and Theoretical Results
It can be seen from the above simulation results that when the 28 m waterproof coal pillar is left, the fault confined water guide zone and floor mining failure zone are connected when the working face advances to 110 m. The rapid increase in seepage pressure leads to the inflow of confined water into the goaf, resulting in water inrush. Therefore, the calculation result of the coal pillar empirical formula for the 23503 working face is smaller than the actual situation. However, the theoretical calculation result is 32 m. Moreover, the normal fault of Dongjiahe is on the right side of the 23503 working face. The secondary fractures and structures of the floor are complicated. According to the theoretical calculation value of waterproof coal pillar with 2-3 times safety factor, the cut hole is arranged at 80 m from the fault.
4.2. Numerical Simulation Study of Waterproof Coal Pillar Width 80 m (Working Condition 2)
To verify whether it is reasonable to arrange the cut hole at 80 m away from the fault, a numerical simulation analysis was carried out in the case of setting up an 80 m waterproof coal pillar at the 23503 working face. Model establishment and excavation were similar to the 28 m modeling process of working condition 1. In this simulation, the cut hole is 80 m away from the right fault and 282 m from the left boundary. 10 steps are taken from the right to the left. For the first and second steps, each mining was 32 m. From the third step to the tenth step, each mining was 23 m and the cumulative mining was 248 m.
4.2.1. Analysis on the Law of the Floor Plastic Zone
As shown in Figure 12, when the working face advanced to 64 m (step 2-15-1), local shear failure occurred at the fault floor under the influence of mining disturbance. Fractures occurred at the interface between the Ordovician ash roof and the fault. When the working face advanced to 110 m (step 4-15-1), the plastic zone of the roof and floor of the mining face expanded. The floor fracture depth on the cut hole side was greater than the floor fracture depth before the working face. At this time, the fault zone was continuously activated upward along the fault. When the working face continued advancing to 133 m (step 5-16-1), the working face reached the limit mining width, and the plastic failure depth of the roof and the floor reached the maximum. At this time, the fault zone was still activated upward along the fault but did not break through the fault zone and hence did not develop outward. When the working face entered the full mining stage until mining was completed, the floor failure depth was less than the maximum failure depth in the nonfull mining stage. Fault activation and confined water uplift only occurred in the fault zone and did not break through the fault zone to develop into the floor mining fracture area.

4.2.2. Analysis on Variation Law of the Floor Seepage Field
As shown in Figure 13, when the working face advanced to 64 m (step 2-15-1), coal seam mining barely affected the change of water pressure in the fault, and high confined water in the floor did not rise along the fault zone. When the working face advanced to 110 m (step 4-15-1), water pressure in the fault changes indicated that the fracture occurred in the fault zone due to the disturbance of floor mining. The confined water was uplifted along the fault zone, and the fault was activated into a water-conducting fault. As the working face continued advancing to 15 6 m (step 6-15-1), Ordovician ash confined water continued to rise upward along the fault zone, and the water pressure in the fault zone reached the maximum value of 0.213 MPa. At this time, the confined water guide zone did not break through the fault zone and therefore did not develop outward. Until the mining of the 248 m (step 10-15-1) working face was completed, the fault activation along the fault zone was higher, and the fracture in the fault zone was more fully developed. However, the confined water does not break through the scope of the fault zone and always moved within the fault zone.

4.2.3. Simulation Result Analysis
From the above simulation results, it can be seen that when an 80 m waterproof coal pillar was left, the floor Ordovician confined water continuously rose along the fault zone. The fault did not extend to the floor mining damage area. The floor mining failure zone was far away from the confined water uplift zone, and no floor water inrush had occurred.
5. Engineering Practice
5.1. Hole Cutting Construction and Groove Wave Detection
Based on the actual situation of the site and through comprehensive research, analysis, and expert meeting to exchange analysis, it was decided that according to the theoretical calculation value of the waterproof coal pillar with 2-3 times safety factor, the cut hole should be arranged at the distance of 80 m from the fault. During the open-cut roadway excavation, three fault structures were revealed, namely, Q1, Q2, and Q3, respectively. Because the Q1 fault was very small and the drop was less than 0.3 m, it did not affect the later excavation. On the other hand, the drops of the roadway where the Q2 and Q3 faults were located were relatively large, and the possible existence of other structures could not be ruled out. Therefore, the groove wave detection was carried out on the working face, and the groove wave observation results are shown in Figure 14.

The detection results showed that the abnormal area is 30 m from 96 m to 129 m from the track roadway to the transport roadway, 81 m from 1071 m to 1132 m and from 1160 m to 1180 m in the transportation roadway, and 102.3 m from 158 m to 255 m and 1060.6 m to 1066 m in the track roadway. According to detection results, it was estimated that the fault at the cut hole extends to the west for a range of 20 m.
The hydrogeological conditions of the 23503 working face are complex, and the driving areas are all located below the Ordovician ash static water level. These areas belong to the pressure-bearing area. According to drilling data, there are different water gushing degrees in the drilling site. Judging from the amount of water gushing and its static pressure, the gushing water recharge had no hydraulic relationship with the underlying Ordovician limestone solution fissure water. It should be fracture water in the lower aquifer group (K2 quartz sandstone) of coal measure strata. However, it cannot be ruled out that the impact of mining causes Ordovician water to flow into the K2 sandstone layer. This forms a large amount of indirect replenishment to the water inrush source. To not affect the safe mining of the working face, the mine carried on the grouting transformation to the floor K2 sandstone.
5.2. Implementation Scheme and Application Effect of Grouting Reinforcement in the Working Face
There are two roadways in the 23503 working face, which are the transport roadway and track roadway. The length of the two grooves in the transport lane and track lane was 1200 meters, and the cut was 180 meters. There are 13 drilling sites and 39 holes in the working face. The specific gravity of the grout was between 1.13 and 1.30 g/cm2, and the total grouting volume was 7970.7 cubic meters (Figure 15). The grouting reinforcement project of the working face starts from 28 November 2018 to 13 October 2019. After completing grouting, the whole working face was explored by the electrical method to meet water prevention and control requirements. Because the Q2 and Q3 faults exposed by the open-cut roadway have not developed to the Ordovician ash roof interface and are coupled with the previous floor grouting transformation, the water inflow of the open-cut roadway was reduced to 20 m3/h. It met the requirements of the detailed rules of water prevention and control and the existing drainage capacity of 200 m3/h of the mine. Currently, the 23503 working face has mined 700 m. From the observation data of the underground water inrush point, it can be seen that the average water inflow was maintained at 40-60 m3/h in the process of mining. The highest water inflow in the whole working face was over 105 m3/h. However, it still meets the existing drainage capacity of the mine and does not cause a water inrush accident.

6. Conclusions
Based on the plastic slip line field theory, this paper puts forward the analytical solution of the limit width of the waterproof coal pillar. Moreover, this paper analyzed the activation law of floor fault in the 23503 working face of the Dongjiahe coal mine. It also verified the theory by numerical simulation and engineering practice. The conclusions are as follows: (1)The numerical models under two different conditions of setting coal pillar width were established using the finite difference software. The laws of deformation and failure of the through fault floor and fault activation were obtained. It was determined that the safety width of the waterproof coal pillar is 80 m. This provided more reasonable data than the empirical value mentioned in the “coal mine safety regulations” for the location of new cut holes(2)When the waterproof coal pillar was set at 28 m, the existence of the through fault hindered the transmission of mining stress in the floor. This made the stress concentration in the floor rock mass between the goaf and the fault to increase, thereby deepening the floor failure depth near the fault. With the advancement of the working face, the fault zone continued to activate upward, and the seepage vector in the zone continued to grow and gradually extended to the failure depth of the floor(3)In the practical application of the research results, three fault structures were exposed during the process of open-hole excavation. Combining these results with grooving wave detection results in the 23503 working face, grouting transformation was carried out on the working face floor. This measure reduced the risk of floor water inrush and had a good effect of preventing water inrush. It also provided engineering references and practical experiences for water disaster prevention and for the control of structural floors under similar conditions in the Weibei mining area for future benefits
Data Availability
Some or all data, models, or codes generated or used during the study are available from the corresponding author by request.
Conflicts of Interest
The authors declare that there is no conflict of interest regarding the publication of this paper.
Acknowledgments
The authors are grateful for the financial assistance provided by the Open Projects of Key Laboratory of Coal Resources Exploration and Comprehensive Utilization, Ministry of Land and Resources (KF2020-5); the National Natural Science Foundation of China (Nos. 41402265 and 51874229); and the Natural Science Foundation of Shaanxi Province (2020JZ-52).