Disastrous Mechanism of Water Burst by Karst Roof Channel in Rocky Desertification Mining Area in Southwest China

School of Geoscience and Technology, Southwest Petroleum University, Chengdu, 610500 Sichuan, China School of Energy Science and Engineering, Henan Polytechnic University, Jiaozuo, 454000 Henan, China Collaborative Innovation Center of Coal Work Safety and Clean High Efficiency Utilization, Jiaozuo, 454000 Henan, China Institute of Mining Engineering, Guizhou Institute of Technology, Guiyang, 550003 Guizhou, China Xintian Coal Mine, Yonggui Energy Development Co., Ltd., Qianxi, 551500 Guizhou, China


Introduction
As one of the most important coal bases in Southwest China, Guizhou is known as the "Southwest Coal Sea" [1,2]. At the same time, it is also the largest and most concentrated area of karst distribution in China. The exposed area of limestone reaches 73% of the total area of the province [3,4]. The underground karst fissures are highly developed in this area. Most of the atmospheric rainfall are not stored on the surface, which flows into the underground karst aquifer [5].
When mining coal seams in karst mining areas, roof karst water seriously threatens the safe and efficient mining of coal in mines [6,7]. Therefore, it is of great significance to study the evolution characteristics of water gushing channels in karst roof of mining areas, which is beneficial for guiding coal mining in rocky desertification mining area in Southwest China.
Many scholars at home and abroad have done lots of research work on the development height and characteristics of the water-conducting fracture zone in the karst coal roof.
In terms of the development characteristics and laws of the water-conducting fractured zone, researchers represented by Liu Tianquan [8] have summarized and obtained the empirical formula for the development height of the waterconducting fractured zone through years of field research and analysis; Qiao et al. [9] systematically summarized the research progress of aerosol water from four aspects, namely, the formation mechanism of roof aerosol water, disastercausing mechanism, water disaster prediction and early warning, and key prevention and control technologies; Yang and Xu [10] comprehensively employed the methods of theoretical analysis, similar material tests, and numerical simulation to obtain the evolution law of the water-conducting fracture zone in a large mining height face; Lai et al. [11] used physical similar material simulation experiments, combined with the total station and borehole peep monitoring, 3DEC, and SPSS statistical analysis software, and obtained the migration law of overlying strata in coal seam mining, the development and evolution of fractures, and the distribution characteristics of water-conducting fracture zones; Zou et al. [12] used FLAC3D software to calculate and analyze the plastic failure zone, displacement field, and stress field distribution characteristics of the surrounding rock above the stope before and after the fully mechanized caving work passes the fault and obtained the formation mechanism of the water channel; Wang et al. [13] used theoretical analysis and field detection to obtain the characteristics of the development height of roof water-conducting fractures under the influence of the key layer structure with the mining thickness, and the development height is affected by both the mining thickness and the key layer structure; Zhang et al. [14][15][16] used the theory of elastic foundation beams to establish a mechanical analysis model for the height of the overburden water-conducting fissure zone of the blockfilled stope; Zhu et al. [17] built a coal seam mining model in the karst cave area and obtained the development characteristics of roof mining cracks during coal seam mining; Wang et al. [18] used similarity simulation and theoretical analysis to propose a composite mechanism model of "elastic thin plate" and "parallel pressure arch" for the migration of overlying strata in high-strength mining under threedimensional spatial conditions; Liu et al. [19] used the borehole television system and borehole simple hydrological observation method, combined with similarity simulation and numerical simulation, and obtained the development characteristics of the overlying water-conducting fissure zone in fully mechanized caving mining in deep and extrathick coal seams. In terms of the height of the waterconducting fissure zone, Y. P. Zhang et al. [20] used a combination of field measurement, numerical simulation, and similarity simulation to obtain the overburden failure height of the deep thick coal seam in the west of Mongolia with the large mining height; Guo et al. [21] used the method of onsite ground drilling flushing fluid leakage and theoretical analysis methods to obtain the height of the watertransmitting fracture zone in top coal mining under soft and hard alternate overburden conditions; Yang et al. [22] comprehensively used downhole borehole water injection loss observation, borehole television, and numerical simula-tion technology to obtain the development height of waterconducting fracture zone in fully mechanized caving mining under thick loose layer and weak overburden. In terms of water flow height prediction, Shi et al. [23,24] combined principal component analysis (PCA), genetic algorithm (GA), and optimized Elman neural network to establish PCA-GA-Elman for the height prediction of water flow fracture zone development. Based on neural network algorithms, Z. H. Li et al. [25] selected mining thickness, mining depth, working face inclination length, coal seam inclination, and overlying rock structure characteristics as the main influencing factors for the height of the waterconducting fissure zone. Based on the particle swarm (POS)-support vector regression (SVR) research method, Xue et al. [26] constructed the Ordos Basin Jurassic coal field water-transmitting fractured zone height prediction model.
The above research results have an important guiding significance for coal mining under water bodies, but there are few studies on the development height and development rules of water-conducting fissures in karst mining areas [27,28]. This paper employs a combination of theoretical analysis, similarity simulation, numerical simulation, and field microseismic monitoring to analyze the source and volume of water inrush from coal seam mining in karst areas. Through similarity simulation and numerical simulation, the process of connection between the water conduction fissures and the original karst fissures is proved, and the water inrush model for coal seam mining in karst areas is proposed, and the coal seam mining in karst mining areas is mastered. Therefore, the evolution characteristics of roof water gushing channels provide an important reference for safe and efficient coal mining in rocky desertification mining area in Southwest China.

Mining and Hydrogeological Conditions
2.1. Water Source of Mine Water Filling. Xintian Mine is located in Qianxi County, Bijie City, Guizhou Province. This area is a typical karst mining area. Mine water is mainly filled with two major water sources, namely, atmospheric precipitation and underground karst water. Atmospheric precipitation is the main source of replenishment for surface water and underground karst water, which restricts the dynamic changes in the flow of surface rivers and mines. According to the meteorological data provided by the Qianxi County Meteorological Bureau, the annual rainfall in the mining area is about 940~1090 mm, and the rainy season is from May to September, accounting for about 80% of the annual rainfall. Under normal circumstances, atmospheric precipitation mainly replenishes the Yulongshan limestone aquifer through shallow weathering fissures, structural fissures, and sinkholes. Continuous rainfall increases the water supply to the underground Yulongshan limestone aquifer, so atmospheric precipitation is the direct water source for filling water in the shallow limestone aquifer of the mine. When the mining fissure is connected to the aquifer, the aquifer becomes a direct source of water for the mine.

Geofluids
The water richness of the aquifer in the Yulongshan section of the Yelang Formation in the mining area is generally not strong and extremely uneven. However, affected by the structure of the fault and lateral fracture zone, the karst is strongly developed, and the underground karst space is connected with the surface karst, especially during rainfall. Atmospheric precipitation and surface water leak along surface creeks and sinkholes to replenish the ground, resulting in a large amount of groundwater enrichment in the limestone aquifer in the Yulongshan section, which becomes an enrichment zone for groundwater. The lower part contains the limestone of the Changxing Formation with an average thickness of 35 m. Because there is a water barrier in the middle, and the stratum is not exposed, deep buried, and poor replenishment conditions, the limestone of the Changxing Formation is a weak aquifer.

Mine Water Filling
Channel. The water-filled channel in the karst mining area is composed of the original karst fissures and mining fissures, which is different from the water-filled channels in ordinary mining areas only by the mining fissures. The original fissures include weathered zones, faults, collapse pits, sinkholes, and underground karst fissures. Mining fissures are the fractures formed by the collapse of the roof of underground coal mining, forming the water-filling channel of the karst mining area. When the mining fissures and karst fissures are connected, karst water continuously flows into the well, seriously affecting the safe and efficient mining of the working face.

Similarity Simulation of Evolution
Characteristics of Water Gushing Channel 3.1. Model Design and Establishment. According to the mining geological conditions of the 1402 working face, the buried depth of working face is about 340 m, and the mining thickness of the coal seam is only 3 m. In consideration, if the similarity ratio is small, the geometric similarity ratio is selected as 1 : 100; and the length, width, and height of the selected similarity simulation test bench are2500 mm × 200 mm × 1300 mm, and the simulated rock layer height is 120 m. According to similarity simulation principles, the top loading of the model fails to simulate the weight of the rock formation. According to the origi-nal model, the loading value is The load value q m on the model is where q p is the prototype unsimulated rock formation pressure, KPa; H is the mining depth, m;  Figure 2. This experiment uses sand as aggregate, calcium carbonate and gypsum as cementing materials, and borax as retarder. According to the calculation method of the simulated strength value of similarity materials, a reasonable ratio of similarity materials in each layer is obtained. According to the cross-sectional area of the model frame, the thickness of the rock (coal), and the geometric similarity ratio, the weight of the similarity material of each rock (coal) layer is calculated (considering the richness factor of 1.2), and the proportion number and parameters are shown in Table 1.       Figure 4. The overburden fissure field is divided into fissure opening area and fissure closure area. The opening angle of the overburden rock near the working face and the open cut is relatively large, and the fissures are more developed. The water conductivity is strong, and the upper fissures in the middle of the gob are closed due to compaction, and the water conductivity is poor. After the coal seam is fully mined, the mining fissures have developed to the top of the model and are connected to the Yulongshan limestone aquifer. The roof water enters the gob through mining fissures. Because it is inclined mining, the karst water in the gob flows to the working face, which is consistent with the water gushing phenomenon that occurs in the actual coal seam mining process.

Field Measurement of Breaking Height in
Overlying Rock

Layout of Microseismic Monitoring Points.
According to the general principle that all measuring points form a spatial body, candidate points for the station layout of the microseismic monitoring system of Xintian Mine have been

Geofluids
selected. The schematic diagram of the designed layout of measuring points is shown in Figure 5.
The optimal plan for the location of the measuring point determined by analysis and measurement calculation is shown in Table 2.

Analysis of Microseismic Monitoring
Results. The distribution of microseismic events was monitored from November 11, 2019, to January 5, 2020. The microseismic events mostly occurred in front of the work and were mostly biased towards the 1402 belt lane, and the energy range was in the range of 0~1000 J. There were only four microseismic events in the mined-out area behind the working face, with energy three times greater than 1000 J. During the monitoring period, microseismic events mostly occur on the leading working face. Large energy events exceeding 100 J are mainly concentrated in the middle of the working face and lagging behind the working face and occur in the upper part of the gob behind the working face. According to the microseismic event profile, the maximum height of the microseismic event is located at the Yulongshan limestone 125 m away from the 4# coal seam. Therefore, the height of the water gushing channel develops to the Yulongshan limestone, and the water gushing channel is connected to the Yulongshan limestone cave.

Analysis of Mine Water Inflow Process
After the working face is fully recovered, there is continuous water gushing to the working face behind the working face. Temporary water pumps are installed on the working face to pump the water to the storage tank of the stop line. After the precipitation, it is discharged into the underground silo through the water pump, shown in Figure 6.
5.1. Analysis of Mine Water Inflow. By analyzing the average monthly water inflow from April 2015 to October 2019 and the atmospheric rainfall data in the area, the relationship between atmospheric rainfall and mine water inflow is obtained, shown in Figure 7. The change in mine water inflow is closely related to atmospheric rainfall, indicating that there is a water inflow channel between the mine and the ground, resulting in an increase in mine water inflow after atmospheric rainfall.

Mine Water Gushing
Process. Based on the hydrogeological conditions of the mining area, a conceptual model of karst roof gushing water is proposed, shown in Figure 8. There are many bead-shaped sinkholes on the surface. The limestone karst fissures in the upper part of the Yulongshan section are relatively developed, and the karst fissures in the lower part are poorly developed. Gas and rain tend to enter underground karst caves and karst fissures through sinkholes and surface karst cracks, becoming a potential threat to coal mining. Xintian Coal Mine mainly mines the 4# and 9# coal seams. With the increase of the mining space    When the distance between the coal seam and the aquifer reaches the critical value of water gushing, the original karst fissures and mining fissures are connected to form a water gushing channel between the surface-karst caveworking face or gob. When 4# coal seam is mined in Xintian Mine, the working face is threatened by roof karst water, and the mining work of 9# coal seam is also threatened by roof karst water.

Conclusion
(1) Atmospheric precipitation in karst areas is the key supply water source for the Yulongshan limestone aquifer, which enters the mine through the original karst fissure and mining-induced fissure. Besides, the Yulongshan limestone aquifer is the direct source of water gushing in the mine (2) There are ultrahigh-conductivity fracture zones in coal seam mining in karst mining areas. When the working face advances 135 m, the fissure develops to the Yulongshan limestone floor. With the continuous advancement of the working face, the height of the water-conducting fissure zone continues to develop upwards. The final monitoring results show that the development height of water-conducting fissure zone is 125 m (3) When the 4# coal seam is fully mined, the mining fissure is connected to the karst fissure, and a water gushing channel is formed between the surface, the karst cave, and the working face (or gob). On this basis, a conceptual model of karst roof gushing in coal mining in rocky desertification mining area is proposed

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.