Abstract

The existence of large-scale weak zone will have a great adverse impact on coal mining in high confined aquifer. Taking the Wutongzhuang Coal Mine which is threatened by high-pressure water as an example, this paper studies the difference between the microseismic events before and after mining and analyzes the influence of the large size weak zone on the coal mining on the confined aquifer. The research results show that the microseismic characteristics of the large soft weak belt are small number of events, the spatial distribution of events is concentrated, and the energy level is large. The amplitude of microseismic events is higher, and the proportion of large events in microseismic events is larger than that of small events; the characteristics of microseismic events caused by mining face mining are that the number of events is more, the distribution of events is loose, the distribution of roof and floor is more, the energy level is less, the amplitude is smaller, and the proportion of small and medium events in microseismic events is larger than that of large events. Due to the joint influence of large-scale weak zone of floor and mining, the floor in the middle area of working face is affected by mining, the number of microseismic events in each aquifer increases suddenly, the karst fissures between the aquifers are further developed, and there is a trend of transfixion. Therefore, measures such as floor grouting should be taken to reinforce the large-scale weak zone before mining.

1. Introduction

Mine water disaster is one of the five major disasters in coal mines in China [1]. Coal mine water inrush accidents account for a large proportion of coal mine disasters, causing heavy losses to the country and people’s lives and property [2]. According to statistics, about 60% of coal mines in China are affected by confined water to varying degrees, and the area and the severity of water disaster rank first among major coal-producing countries in the world [3]. Take Wutongzhuang Coal Mine as an example. It is typical coal mining on confined water in the “three underground and one waterborne circumstances” special mining [4]. In recent years, with the increase of mining depth, the Wutongzhuang Coal Mine has been more and more threatened by high-pressure o.l. water [5, 6]. Even if water disaster prevention measures are taken, water inrush accidents still occur [7]. Therefore, studying the mining of coal seam on confined water and avoiding the occurrence of water inrush from the floor are an urgent problem to be solved [8].

China is a country that started to study coal mining on high confined water relatively early in the world and has achieved a lot of results. Xie et al. [9] started from the mechanism of hydraulic fracturing and analyzed the permeability changes to infer the risk of water inrush from the floor and introduced the concept of damage. Xu et al. [10] studied the mining failure of coal mining floor under high confined water and the formation and evolution of water inrush channels through on-site monitoring and numerical simulation. Li and Xu [11] established a physical model to prevent water inrush disasters from the collapse column using the comprehensive methods of theoretical analysis, similar simulation, and digital simulation. Zhao et al. [12] put forward the key technology of advanced treatment of o.l. water damage in the ground area, and it has been widely used. Dong et al. [13] and Wang et al. [14] put forward the idea of using a weathering zone at the top of the Ordovician limestone which has better water blocking performance as a water-repellent layer, and they conducted a study on its availability and established a criterion for the thickness of the rock section on top of Ordovician limestone. Xu et al. [15] used borehole ultrasonic detection, core indoor ultrasonic transmission, and point load test methods to study the damage and destruction degree of the floor rock before mining working face in the upper and lower coal seams of the confined water, and they obtained the rock elastic modulus. The density, lithology, and completeness of rock mass show a nonlinear positive correlation. Ma [16] combined method of transient electromagnetic method and radio wave perspective method to comprehensively detect the internal water-conducting structure of the working face under the effect of high confined water and achieved good results. Wei et al. [17] and others used numerical simulation and similar simulation methods to analyze the degeneration characteristics of the roadway under the effect of high confined water. Chen et al. [18] used the knowledge of mechanics combined with the “drilling double-ended side leakage plugging device” to study the evolution of cracks in the mining rock mass under high confined water. Xu and Gao [19] established a mechanical model to theoretically analyze the mechanism of water inrush from a working face cross-fault under high confined water conditions. Wang et al. [20] analyzed the failure law of complex rock during UCA mining and proposed a method for predicting the danger zone of support fractured water inrush. Yu et al. [21] used the fractal dimension method to partition the structural complexity of the study area and evaluated the risk of water inrush from the deep coal seam floor by mining above the confined aquifer in the mining area. Duan and Zhao [22], by analyzing some water inrush cases, put forward four main factors that affect water inrush from coal seam floor and established a mechanical model for the homogeneously broken and weak rock layer of the floor water-insulating layer. Li [23] introduced the development history of various geophysical prospecting technologies such as underground electrical methods, transient electromagnetics, electrical perspective, radio wave tunnel perspective, and channel wave seismic in China and discussed the development direction of the discipline.

The microseismic monitoring method has been widely used in coal mines. In recent years, many scholars have made a series of achievements in the field of microseismic events. Ge and Han [24] proposed an imaging method for reverse-time migration of the structure below the microseismic zone. Lou et al. [25] used microseismic monitoring methods to analyze and study the stress conditions in the goaf and the rearrangement of the fracture zone. Wang et al. [26] proposed a source-independent objective function based on convolving reference traces, which avoids the cycle jump phenomenon caused by unknown sources and improves the accuracy of estimating the location of the seismic source. Hao et al. [27, 28] studied the fracture characteristics of coal under dynamic load through uniaxial compression tests and studied the influence of bedding cleats on it, and they found in subsequent studies that microseismic events occur when the coal bedding dip is 45° at most. Hao et al. [29] used a quantitative analysis method to construct three physical models and combined numerical simulation to analyze the influence of the main principal stress directions on the long axis of the tunnel on the stability of the underground tunnel. Zheng et al. [30], based on the limit equilibrium method and genetic algorithm, proposed a new method for stability assessment of ABRS considering seismic inertial forces. The effects of seismic and mechanical properties of rock strata on the failure mechanism of ABRS are discussed. Zhang et al. [31] proposed a new method based on VMD to distinguish between coal and rock fracture and blasting vibration microseismic signals by processing the signal. Chen et al. [32] evaluated the effect of hydraulic fracturing by analyzing the number of microseismic events and energy changes caused by high-pressure hydraulic fracturing. The Brune model and grid search method were used to calculate the source parameters, and the relationship between M-0 and other source parameters was analyzed. The relationship between M-w (moment magnitude) and M-L (local magnitude) of microseisms induced by hydraulic fracturing is obtained. Amad et al. [33], through the joint inversion of source location and distance tensor, reconstructed the space and role of (micro) seismic events and constructed a full-waveform technology. Estay et al. [34] obtained the b value of seismic activity related to mining by studying the time decay law of seismic activity after rockburst and large-magnitude events occurred in the aftershock sequence of the mine. Liu et al. [35] proposed a microseismicity-based method for the dynamic estimation of the potential rockburst scale during tunnel excavation.

In summary, it can be seen that there are a large number of results in exploration, mechanical models, and similar simulations in high-pressure water mining. However, there are few studies on the microseismic characteristics of excessively large-sized weak zones before mining. Therefore, this paper uses the data collected by the microseismic monitoring system to study the difference between the large-scale weak zone near the F702-11 fault in the front floor before and after the mining of Wutongzhuang Coal Mine 182703 working face, and it analyzes the impact of large-size weak zones on coal mining above confined aquifers and provides technical guidance for coal mine safety production.

2. Project Profile

Wutongzhuang Coal Mine is located in Handan City, Hebei Province. It is one of the main mines of Jizhong Energy Fengfeng Group, with a design production capacity of 3.0 Mt/a. The hydrogeological conditions are extremely complex, and the ground temperature of the mine affected by the geological structure is higher than that of other mines in the group [36, 37].

The 182703 working face is located in the seventh mining area, the east is the designed 182704 working face, the west is the 182702 mined area, the south is the F2-2 fault, and the north is the south belt transportation lane with no other excavation condition. The No. 2 coal seam is mined at the working face, the coal bulk density is 1.35 t/m3, the mining area is 145071 m2, the industrial reserves are 0.722 Mt, and the recoverable reserves are 0.686 Mt. The strike length of the working face is 683 m, with an average of 683 m; the inclination length is 94 m~252 m, with an average of 212 m; the inclination of the working face is 1°∼36°, with an average of 25°; the maximum coal seam thickness is 4.2 m, the minimum is 3.0 m, and the average is 3.69 m. The buried depth of the working face is about 660 meters, and the direct top lithology is sandy shale, which is brittle and compact.

3. Distribution of Large-Size Weak Zones in the Working Face

According to the data provided by the Wutongzhuang Coal Mine Geology Department, in the 182703 working face, there are many weak zones dominated by faults, and the fault occurrence is shown in Table 1. Figure 1 shows the distribution of abnormal weak bands on the working face measured by the transient electromagnetic method. Figure 2 shows the cross-sectional view of the weak fault zone in the working face.

4. Comparative Analysis of Microseismic Events before and after Mining

4.1. Construction of the Microseismic Monitoring System
4.1.1. Layout of the Microseismic Monitoring System

There are 11 geophones arranged in the upper and lower tunnels and 2 underground monitoring substations for a total of 15 channels on the 182703 working surface of Wutongzhuang Mine. Among them, 3# and 7# are three-component detectors, the others are single-component detectors, 1#∼5# detectors are located in the lower groove, and 6#∼11# detectors are located in the upper groove. There are 24-bit digital microseismic acquisition substations, respectively, in the 182703 transportation roadway and the 182703 return airway, with a sampling frequency of 5 kHz. The monitoring layout is shown in Figure 3. In order to avoid the influence of the loose circle of the roadway, the geophones are buried at a depth of 6 m and the spacing is about 100 m.

4.1.2. Transmission and Processing of Microseismic Monitoring Data

The microseismic monitoring data is continuously collected in real time, and the underground optical fiber is directly transmitted to the surface switch of Xingdong Coal Mine and transmitted to the microseismic monitoring and monitoring center of the Academy of Coal Science through the Jizhong Energy LAN (Figure 4).

The 182703 working face monitoring system was set up in April and began to collect data. The working face started mining on May 17th. We compared the number of events, energy, amplitude, and spatial distribution in the microseismic monitoring data collected from April 1st to April 10th and May 18th to May 27th to analyze the impact of mining work on microseisms.

4.2. Comparative Analysis of the Number of Microseismic Events before and after Mining

Figure 5 shows the statistics of the number of events before and after mining. From April 1st to April 10th, before mining, there were a total of 340 events collected. The number of events on April 10 was the least, with 15 events, and the number of events on April 3 was the most, with 55 events, with an average of 34 events per day. From May 18 to May 27 after the mining, there were a total of 820 microseismic events collected, of which the number of events on May 19 was the least, with 47 events, and the number of events on May 24 was the most, with 132 events, with an average of 82 events per day. After querying the data of the mine, it was concluded that the increase in microseisms on May 24 was caused by the roof first weighting.

From the comparison of the data in Figure 5, in terms of quantity, it can be seen that the number of microseismic events affected by mining activities has increased significantly in May.

4.3. Comparative Analysis of Microseismic Energy before and after Mining

Figure 6 shows the comparison of microseismic energy before and after mining. Figure 6 shows that, from April 1st to April 10th, before mining, the minimum energy was 4691 J on April 8, the maximum energy was 4614966 J on April 5, and the average daily energy was 1002882 J. After mining from May 18 to May 27, the lowest energy was 2801 J on May 21, the highest energy was 658864 J on May 27, and the average daily energy was 119664 J. It can be seen that the energy in April is an order of magnitude higher than that in May, so it can be concluded that before the working face was mined in April, the cause of the microseismic event was due to the existence of a weak zone in front of the working face, which was affected by the geological tectonic movement of the coal mine, causing a large microseismic event, so the energy is huge. The microseismic event in May is caused by mining activities, although the number is large, but the energy is much smaller than the energy generated by the weak zone.

Figure 7 is a graph of cumulative energy and cumulative events. It can be seen from Figure 7 that the cumulative number of microseismic events is higher after mining, but in terms of cumulative energy of microseisms, the large-size weak zone was affected by the geological tectonic movement resulting in the dislocation of the large-size weak zone before mining, which caused a large microseismic event, resulting in a surge in energy, which was an order of magnitude higher in energy data than after mining.

It can be seen from the comparison diagrams of Figures 57 that the microseismic events monitored in April were caused by the large-scale weak zone in front of the working face that was affected by the movement of geological tectonic movement and caused the large microseismic event. Its manifestation is large energy and small number of events. The microseismic events monitored in May are microruptures caused by the mining activities. The manifestation is low energy and large numbers of microseismic events.

4.4. Comparison of the Amplitude of Microseismic Events before and after Mining

Figure 8 is a comparison diagram of the amplitude of microseismic events before and after mining. It can be seen that, from April 1st to April 10th, before mining, there was a microseismic event with a minimum amplitude of 0.8597 on April 1st and a maximum microseismic amplitude of 4.2951 on April 7th with an average amplitude of 2.0751, and the microseismic amplitude distribution was more even. From May 18th to May 27th, after mining, there was a microseismic event with a minimum amplitude of 0.0934 on May 23rd and a microseismic event with a maximum amplitude of 4.8481 on May 18th, with an average amplitude of 1.3903, and there were more small amplitude events. Therefore, it can be seen from the comparison diagram of the amplitude of the microseismic event that before the working face was mined in April, the large-scale weak zone was affected by the geological tectonic movement and caused a large microseismic event. Although the number of events is small, the energy and amplitude are relatively large. In May, the working face was affected by the continued mining activities to cause small ruptures. There were many microseismic events, but most of them were small events with small amplitude and low energy.

4.5. Analysis of the Spatial Distribution of Microseismic Events before and after Mining
4.5.1. Distribution Characteristics of Microseismic Events in Each Aquifer

Figures 9 and 10 are, respectively, the spatial distribution of microseismic events. The sphere in the figure represents the microseismic event, the size of the sphere represents the amplitude of the microseismic event, and the black curved surface is the No. 2 coal seam. By comparing the spatial distribution map of microseismic events before and after mining, it can be seen that the microseismic events before mining are mainly floor events, and the distribution range of the events is concentrated between −400 and −650 meters. Although the microseismic events caused by mining activities are still mainly floor events, the roof events have obviously increased. From the distribution range, it can be seen that the distribution range of microseismic events after mining is mainly concentrated between −400 and −800 meters, but there are still a small number of microseismic events extending to the top −200 meters and the bottom −1000 meters. Therefore, the distribution of microseismic events caused by mining is wider, with a large number of microseismic events on the floor and the roof, but the magnitude of the microseismic events is smaller than that of the large-scale weak band movement.

The microseismic event distribution diagram (Figures 1114) was established through the microseismic data collected by the microseismic monitoring system. From Figure 11, it can be seen that the microseismic events generated by tectonic activities in April were distributed in the middle of the area on the plane, and the events that occurred in the Ordovician limestone were mainly concentrated in the area of the F702-11 fault and extending to F2-2 fault area, the distribution range of events that occurred in the Ordovician limestone is relatively concentrated. It can be seen from Figure 12 that the number of events that occurred in the Ordovician limestone after mining in May increased significantly compared to Figure 11, and they were mainly distributed near the F702-11 fault area. However, the number of events that occurred in the Ordovician limestone was more dispersed than that in May, from the shallow part of 2# coal to the Yeqing layer mainly concentrated near the transport roadway and the return air trough and extending to the open-off cut. The distribution of events is relatively concentrated. The analysis reason is that, from the shallow part of 2# coal to Yeqing layer mainly affected by mining pressure, Ordovician limestone interval is mainly affected by the o.l. water pressure.

Figures 13 and 14 are cross-sectional views perpendicular to the transport roadway and the return air trough. From Figures 13 and 14 it can be seen that microseismic events occur in all aquifers, but the proportion of events that occurred in the Ordovician limestone in the distribution of microseismic events in May has increased significantly compared to the proportion of events that occurred in the Ordovician limestone in the distribution of microseismic events in April, and the proportion of occurrences of microseismic events in the floor of Shanfuqing-the floor of Daqing has decreased. At the same time, combining the microseismic events of each aquifer in Figure 15, it can be seen that the large-scale weak zone was distorted by the geological tectonic movement in April, which mainly occurred in the floor of Yeqing-the floor of Fuqing area, which caused a significant increase in the proportion of microseismic events in this area, and after the mining in May, since no tectonic activity has occurred so far, it is only affected by mining activities, so the distribution of microseismic events in each aquifer is relatively even.

4.5.2. The Relationship between the Variation of Water Pressure in Different Aquifers and Microearthquakes

The water pressure of aquifers with different depths is different, which leads to different microseismic monitoring data. According to the water pressure data collected by the mine hydrological monitoring unit, there is a clear difference in water pressure between the Yeqing limestone aquifer and the mountain blue limestone aquifer. The water pressure of the Yeqing limestone aquifer ranges from 1.6 to 7.1 MPa. The water pressure of the rock aquifer ranges from 1.3 to 8.8 MPa. According to the microseismic profile, it can be seen that the number of microseismic events in the mountain bladder limestone aquifer is more than that of the Yeqing limestone aquifer, which can indicate that there is a certain connection between water pressure and microseismic events, and the increase in water pressure will lead to the increase in the number of microseismic events.

4.6. Microseismic Characteristics When Passing through the Grouting Abnormal Area

There are two abnormal grouting areas on the working face. When the working face passes through the abnormal grouting area, the microseismic monitoring data has changed significantly compared with the normal grouting area. When the working face passes the 1# grouting anomaly zone, the accumulated energy of microseismic events from June 26th to July 10th, 2019, is 16575760 J, of which the minimum energy of the microseismic event on July 9th is 7322 J and the maximum energy of July 5th is 13242289 J. The average is 1105051 J. When the working face passes the 2# grouting abnormal area, the accumulated energy of microseismic events from September 6th to September 20th, 2019, is 29697036 J, of which the lowest energy on September 6th is 1906 J and the highest energy on September 15th is 7154687 J, with an average of 1979802 J. Therefore, it can be seen that when the working face gradually approaches and passes through the grouting abnormal area, the number of microseismic events will increase significantly. However, from the perspective of the energy of microseismic events, most of the increased microseismic events (when the working face approaches the grouting abnormal area) are mining activities. The roof collapse event caused by the incident is manifested as a large number of small energy events; when passing through the grouting abnormal area, the surrounding rock deep in the floor ruptures, resulting in microseismic events, manifested in the form of large energy and large number of events.

5. Conclusion

Through the quantitative processing of the data collected by the microseismic monitoring system, a spatial distribution diagram of microseismic events, a comparison diagram of the number of events, and a comparison diagram of amplitudes have been established, combined with the actual measurement before and after mining in the Wutongzhuang Coal Mine. The following can be concluded:(1)There is a large-scale weak zone in the floor of Yeqing-the floor of Fuqing area in front of the working face, and the cause of its dislocation is geological tectonic movement.(2)When the large-scale weak zone under the floor of the working face is misaligned, the microseismic feature is that the number of events is small, the event spatial distribution is concentrated, the energy magnitude is high, the amplitude is high, and the proportion of large events in microseismic events is greater than that of small ones. The characteristics of microseismic events caused by working face stopping are the following: the number of events is large, the distribution of events is relatively loose, the roof and floor are more distributed, but the energy level is small, the amplitude is also small, and the proportion of small and medium events in microseismic events is greater than that of large ones.

It can be seen that there is a weak zone dominated by faults under the floor of Wutongzhuang Coal Mine 182703 working face. When it is misaligned by the influence of geological tectonic movement, it will trigger a large microseismic event, which is different from the microseismic event caused by mining activities. The microseismic event caused by the weak zone dislocation has a large energy level and is more dangerous to the mine. It will aggravate the changes in the stability of the floor aquifer during the mining and increase the possibility of water inrush from the coal floor. Therefore, measures such as floor grouting should be taken in time to enhance coal mine safety.

At the same time, it should also be noted that there is no tectonic activity caused by mining at the current working face excavation. Therefore, this article only compares and analyzes the microseisms caused by tectonic activity and those caused by mining. If there is tectonic activity caused by mining, it is inferred that the slight earthquake may be more severe.

Data Availability

The data used to support the findings of this study are available from the corresponding author upon request.

Conflicts of Interest

The authors declare that they have no conflicts of interest.

Acknowledgments

This work was supported by the National Natural Science Foundation of China (51804309), Yue Qi Young Scholar Project (2019QN02) and Distinguished Scholar Project (2017JCB02) from China University of Mining and Technology (Beijing), and Joint Fund of State Key Laboratory of Coal Resources and Safe Mining and Beijing Outstanding Young Scientist Program (Grants nos. SKLCRSM20LH01 and BJJWZYJH01201911413037).