Investigation on the dynamic evolution and evaluation of floor microseismic responses during extra-thick coal seam mining
Haorui Wang1*, Shangxian Yin2, Huiqing Lian2, Min Cao2*, Enke Hou3, Xu Wang1, Guofu Xing4
1https://ror.org/046fkpt18grid.440720.50000 0004 1759 0801College of Geology and Environment, Xi’an University of Science and Technology, Xi’an, 710054 China; https://ror.org/0096c7651grid.443279.f0000 0004 0632 3206Hebei State Key Laboratory of Mine Disaster Prevention, North China Institute of Science and Technology, Beijing, 101601 China
2https://ror.org/0096c7651grid.443279.f0000 0004 0632 3206Hebei State Key Laboratory of Mine Disaster Prevention, North China Institute of Science and Technology, Beijing, 101601 China
3https://ror.org/046fkpt18grid.440720.50000 0004 1759 0801College of Geology and Environment, Xi’an University of Science and Technology, Xi’an, 710054 China
4Inner Mongolia Zhongtai Suancigou Mining Co., Ltd, Ordos, 017000 China
* Corresponding author.
Trans. Intell. Comput. Syst. 2025, 3(2), 2775; https://doi.org/10.58304/24kr5770
Abstract
The mining of extremely thick coal seams causes strong disturbances to the floor strata. When a strong aquifer exists beneath the coal seam and the effective aquiclude is weak, the risk of floor water inrush increases significantly. In recent years, microseismic (MS) monitoring technology has been widely applied, and using MS monitoring to assess floor water inrush risk has proven to be a reliable and effective approach. To address the high-risk problem of floor water inrush during the extraction of extremely thick coal seams, this study employed MS monitoring to investigate the evolution characteristics of floor MS activity and its associated water inrush risk. The MS response characteristics induced by thick seam mining were analyzed, and a floor water inrush risk evaluation model was established based on the Analytic Hierarchy Process–Entropy Weight Method (AHP–EWM). The evaluation results were consistent with field monitoring data. The results indicate that most floor MS events in the 61607 working face of the Longwanggou Coal Mine were concentrated in the Taiyuan Formation. The MS monitoring indicators, such as the extent of the advanced influence zone, floor failure depth, and total daily released energy, showed strong correlations with observed data including water inflow and mine pressure. By summarizing the key controlling factors of floor water-inrush risk, this study departs from traditional static geological parameter indicators and instead focuses on dynamic indicators derived from the activity, spatiotemporal characteristics, and intensity characteristics of floor MS events. Integrating multiple influential factors of floor water inrush, a quantitative evaluation system for floor water-inrush hazard in extremely thick coal seam mining was established using the AHP–EWM combined subjective–objective weighting method, which shows good consistency with field monitoring results. This study provides theoretical and data support for analyzing the dynamic process and influencing factors of floor damage and fracture evolution under thick seam mining. It offers a new approach for predicting and forecasting water inrush risk, enriches the application of MS monitoring in mine floor water hazard prevention, and holds significant practical importance for promoting safe coal mine production.
Introduction
The coal resources in the Jin–Shan–Meng mining area of China are characterized by large reserves and thick coal seams1. With the continuous progress of coal mining technology, a high-intensity and high-efficiency mining mode has been gradually established. Compared with conventional thin and medium-thick coal seams, the mining of extra-thick coal seams is associated with higher strata pressure, stronger ground stress effects, and greater disturbance to the floor, often resulting in extensive and deep floor fracturing2–6. Moreover, the Carboniferous–Permian coal seam floors in China generally overlie Ordovician limestone aquifers with high water pressure. The fractures in the floor can easily connect with the aquifers and trigger water inrush disasters, posing serious threats to mine safety. Floor water hazards caused by high-pressure confined aquifers have become one of the main factors restricting the safe and efficient extraction of coal resources in China7,8. Therefore, it is necessary to investigate the dynamic failure evolution of the floor during extra-thick coal seam mining.
Many scholars have conducted studies on extra-thick coal seam mining and achieved valuable results. Zheng et al.9 analyzed the water inrush mechanism of the Huangyuchuan Mine during top-coal caving and clarified the coupling mechanism between Ordovician limestone permeable channels and floor failure, proposing a directional long-borehole exploration technique to prevent Ordovician limestone water hazards during mining. Shao10 systematically analyzed the roof and floor structure and failure characteristics of fully mechanized top-coal caving faces and established an elastic mechanical model of the floor. The study showed that the maximum failure depth of the floor is positively correlated with the face width and surrounding rock stress and negatively correlated with rock compressive strength. Li et al.11 based on the gas–water two-phase flow similarity theory in fractured rock masses, established a CFDEM numerical model to simulate the fault-induced floor water inrush and roof–floor coupling failure in extra-thick coal seam mining. Their results indicated that the root cause of water inrush lies in the instability of key high-level strata, subsidence of masonry beams, and secondary fracturing of the floor. Wang et al.12–14 used theoretical analysis, numerical simulation, and field measurements to study the stress and fracture evolution of the floor during fully mechanized top-coal caving. Zhang et al.15–18 applied distributed optical fiber sensing and cross-hole resistivity CT testing to obtain in situ data on strain and geoelectric field responses during floor failure in extra-thick coal seam mining. Guo19 analyzed the three-dimensional dynamic evolution of floor disturbance using borehole stress–strain monitoring. Cao et al.20 derived a mechanical relationship between floor failure depth and coal seam thickness and established a floor water inrush risk evaluation model using interval grey optimal clustering analysis.
In addition to floor failure research, Zhu et al.21 studied the mechanism of floor heave in fully mechanized top-coal caving faces and proposed an underground regional hydraulic fracturing technique for controlling floor heave. Zhang et al.22,23 investigated the characteristics of roof caving structure and dynamic loading mechanisms, analyzing the strong dynamic pressure phenomena in working faces of extra-thick seams. Kong et al.24,25 examined the evolution of fractures and seepage characteristics during mining and developed a predictive model for the height of the water-conducting fracture zone. Yu et al.26 explored the multi-field coupling mechanism of disaster formation and proposed a collaborative control technology.
Most existing studies have been based on theoretical analysis, numerical simulation, and field measurements, which have limited continuity and fail to fully represent the dynamic evolution process of floor failure. As a real-time, dynamic, and efficient geophysical monitoring technique, MS observation has been increasingly applied in floor safety studies27. Yin et al.28 used machine learning and deep learning models to analyze MS events and successfully predicted water inrush incidents. Xin et al.29 applied cumulative and kernel density indicators of MS events to interpret fault activation characteristics in the floor. Jia et al.30,31 established a water inrush risk evaluation method based on MS events, while Gai et al.32,33 introduced weighted variation coefficients and floor safety indices to improve the risk evaluation model and zoning accuracy. Xu et al.34,35 proposed an all-time–space early warning method for floor water inrush using MS monitoring, achieving accurate prediction at different mining stages. Zuo et al.36–38 combined theoretical analysis, MS monitoring, and numerical simulation to identify and predict potential water inrush channels based on floor rock fracture characteristics.
Although these studies have explored various aspects of floor water hazard prevention using MS technology and developed related theories and prediction methods, the understanding of floor failure characteristics and evaluation systems under special mining conditions remains incomplete. This is especially true for extra-thick coal seam mining in western China, where the large mining height produces strong and complex disturbances to the floor. Most hazard evaluations are based on static geological parameters such as aquiclude thickness, aquifer water pressure, and geological structures. The core logic relies on assessing the potential for water inrush from inherent geological conditions, but dynamic indicators are generally lacking. Therefore, taking the Longwanggou Coal Mine in the Zhungeer Coalfield as the research area, this study applies MS monitoring technology to analyze the dynamic evolution of floor failure. Based on the AHP-EWM method and combined with field water inrush data, a MS dynamic indicator system was established to evaluate floor water inrush risk, enabling real-time monitoring and quantitative assessment of the hazard. The research results can further enhance the application of MS technology in floor water hazard prevention and provide theoretical support for similar mining operations, offering practical significance for improving coal mine safety.
Study area overview
Fig. 1 Working Face 61607, Longwanggou Coal Mine, Inner Mongolia, China.
Fig. 2 Comprehensive stratigraphic column of roof and floor at working face 61607.
The distance between the floor of the working face and the top boundary of the Ordovician limestone ranges from 33 to 53.3 m, with an average of 44.57 m. The Bauxite-bearing mudstone of the Carboniferous Benxi Formation, with an average thickness of about 8.22 m, serves as an effective aquiclude between the coal seam floor and the Ordovician aquifer.
Research methods
Microseismic monitoring system
Fig. 3 Schematic diagram of the MS monitoring system.
The KJ551 MS system employs high-sensitivity sensors (100 V/m/s) with a wide frequency range (0.1–600 Hz). For signal transmission, an advanced combination of cable and optical fiber technology is used. Cable transmission can cover distances of at least 10 km, while optical fiber transmission can reach up to 60 km, meeting the requirements of large-scale mining operations.
Based on the actual conditions of the 61607 working face, two MS substations were deployed, each installed in a well-conditioned chamber. The boreholes were drilled to a vertical depth of approximately 2.4 m, and the geophones were anchored along the full length using φ22 mm cement bolts with resin adhesive at the bottom. A cement pad was placed at the borehole mouth, leaving the bolt exposed by 0.05–0.1 m. The geophones were connected to the ends of the bolts using universal joints. Except for the first sensor, which was located 50 m from the cut-through, the spacing between adjacent sensors was approximately 100 m. After installation, each sensor underwent on-site hammer testing to ensure signal quality.
Once the substations and sensors were installed, calibration blasting tests were conducted in the two main entries of the working face to ensure precise sensor positioning. The elastic wave velocity of the surrounding rock was accurately determined based on the blasting signals received by the sensors, allowing the calculation of the actual average wave velocity of the rock mass. Using this wave velocity, the source location of the test blasts was analyzed to validate the system’s positioning accuracy.
The obtained velocity model and the actual sensor positions were imported into the processing software. During actual mining operations, the time differences of arrival at each MS sensor were used to calculate the source location and energy of events. This enabled determination of the spatial–temporal characteristics and energy levels of each MS events, providing the basis for corresponding early-warning measures.
Analytic hierarchy process–entropy weight method (AHP–EWM)
Analytic hierarchy process
Constructing a hierarchical structure model: The decision problem is decomposed into a hierarchical structure of goal layer–criteria layer–indicator layer.
Building the judgment matrix: For elements at the same hierarchical level, a relative importance scale is used. Experts are invited to score the importance of each element, forming a judgment matrix A=(aij) n×n, where n is the number of elements at that level. The matrix satisfies reciprocity and reflexivity.
- Weight calculation and consistency test: The maximum eigenvalue λmax and its corresponding eigenvector ω of matrix A are determined. The weights ωAHP,j for each element are calculated accordingly. The consistency index (CI) and consistency ratio (CR) are computed to ensure CR<0.1, avoiding logical contradictions in the judgments. The formulas are:12
Entropy weight method
Data normalization: Since the indicators differ in dimension and magnitude, the raw data are normalized.
Calculation of information entropy: For the i-th sample of the j-th indicator, the probability pij is calculated. The entropy value Ej of each indicator is then obtained using the information entropy formula.
- Determination of indicator weights: Smaller entropy Ej indicates larger data dispersion and more effective information. The weight of the j-th indicator is calculated as:345
Vulnerability index method
Results
Overview of microseismic monitoring
Table 1 Frequency statistics of MS events at working face 61607 (Apr 1–Aug 31, 2023).
| Layer | Frequency | Percentage/% | Cumulative percentage/% | |
|---|---|---|---|---|
| Roof | Above the roof | 280 | 4.33 | 4.33 |
| Floor | Carboniferous Taiyuan Formation | 6008 | 92.92 | 97.25 |
| Carboniferous Benxi Formation | 112 | 1.73 | 98.98 | |
| Ordovician | 48 | 1.02 | 100 |
Before mining, hydraulic pre-fracturing was applied to the roof, causing partial roof failure. The compromised roof integrity interrupted the stress transmission paths, resulting in a relatively small number of roof MS events, accounting for only 4.33% of the total events. Therefore, this study focuses on floor failure, and roof-related events are not further discussed.
Statistical analysis shows that MS events in the floor accounted for 92.92% of all recorded events within the study interval. This indicates that the closer the rock mass is to the coal seam floor, the more significant the disturbance under mining-induced stress, making it more susceptible to failure and stress-concentration fracturing.
The Carboniferous Benxi Formation and Ordovician Formation accounted for 1.73% and 1.02%, respectively. Although floor failure is mainly concentrated in the Taiyuan Formation, in areas with concealed structures, local floor disturbance can extend deeper, breaking through the Benxi Formation aquiclude and reaching the Ordovician aquifer. According to MS monitoring, the deepest disturbed point in the floor was 83.55 m below the coal seam floor, reaching the Ordovician limestone strata.
Fig. 4. Frequency statistics of MS events at working face 61607 (Apr 1–Aug 31, 2023).
Spatiotemporal response characteristics of floor microseismic events
Floor MS events represent energy-releasing fractures in the floor rock and can directly reflect the degree of floor failure across different strata in both space and time. In the 61607 working face, the floor aquiclude is relatively thin. The propagation of primary fractures or the development of concealed structures can easily connect to the floor aquifer, increasing water inflow or even causing water inrush. To better understand the impact of mining on floor disturbance, the spatiotemporal distribution characteristics of floor MS events were analyzed from the perspectives of time, space, and energy.
Temporal distribution characteristics
The Carboniferous Benxi Formation lies below the floor failure zone and above the Ordovician limestone aquifer. This aquiclude mainly consists of bauxite-bearing mudstone and mudstone, with the bauxite-bearing mudstone layers being thick, lithologically stable, and possessing high uniaxial compressive strength, making it an effective floor aquiclude. If this layer is compromised, the integrity of the floor aquiclude may be lost, forming a water-conducting channel through which pressurized water can rise along fractures, potentially causing water inrush. Therefore, the Benxi Formation can be treated as a monitoring alert layer in MS studies.
Fig. 5 Histogram of MS events frequency in different floor strata.
Comparison with real-time monitoring of water inflow and Ordovician limestone water levels (Figure 6) shows two significant water level drops during the study period. The first occurred from April 30 to May 10, without a corresponding increase in inflow. MS events frequencies in the Benxi alert layer and Ordovician strata remained low, suggesting that floor failure was confined to the Taiyuan Formation and had not yet connected to the Ordovician aquifer. The second drop occurred from June 30 to July 20, during which water inflow increased significantly between June 30 and July 31. From July 15, MS events frequency in the Benxi alert layer and Ordovician strata rose sharply, indicating that floor failure had developed into the Benxi Formation and the top of the Ordovician aquifer. Consequently, water inflow increased markedly, demonstrating a strong hydraulic connection between floor fractures and the aquifer. However, no water inrush occurred, suggesting that although mining-induced fractures communicated with the aquifer, some fractures were reclosed under high-stress conditions due to the large mining height, limiting water-conducting capacity.
Planar distribution characteristics
The planar distribution of floor MS events within the 61607 working face study area is shown in Figure 7. Events are primarily concentrated 640 m ahead of the mining face, with an outward extension of approximately 200 m along both roadway sides.
Examining the planar distribution of MS events by stratigraphic layer, the Taiyuan Formation exhibits the widest distribution. Events cluster near the SF10 fault and adjacent 61605 working face goaf, with concentration along the roadway sides significantly higher than in the center. This indicates that the post-mining surrounding rock experiences substantial vertical stress, leading to severe plastic deformation and release of potential energy. The clustering of MS events near the SF10 fault suggests that geologically weak zones are more susceptible to disturbance and failure.
MS events in the Benxi Formation and Taiyuan Formation are mainly concentrated in the central part of the working face. The distribution of MS events in the Benxi Formation is broader than that in the Ordovician strata. During the early stage of mining, MS events clustering in both layers was low, gradually increasing toward the mid-section around mid-July. Combined with water inflow and Ordovician water level data (Figure 6 and Figure 7), the increase in MS events frequency in the Benxi alert layer corresponds with a drop in Ordovician water level and an increase in water inflow, lasting until the end of July. After the mining face passes through this region, the Ordovician water level rises while water inflow decreases.
Fig. 6 Variation trend of water inflow and ordovician limestone water level in the working face.
Fig. 7 Plan distribution map of floor MS events.
The MS events clusters in the Benxi and Ordovician strata form irregular elliptical patterns. This is attributed to roof pressure relief after mining, which causes the main stress to concentrate in the sidewalls. Consequently, lateral shear stress is imposed on the floor beneath the goaf, resulting in floor compression and expansion toward the goaf. This lateral stress concentration in the central part of the working face leads to floor disturbance and fracturing.
Fig. 8 Schematic diagram of the planar influence range of MS events.
Fig. 9 influence distance before and after the mining line.
From April 1 to August 31, the advance influence distance continuously expanded with mining progress. During the early mining stage (April), the advance influence distance was relatively short, mostly within 200 m, indicating that early mining-induced stress had limited effect on floor disturbance. In the mid-mining stage (May–June), the advance influence distance extended to 300–400 m, and during the late mining stage (July–August), it reached a maximum of approximately 500 m. This indicates that the large mining height and strong disturbance of thick coal seam mining result in an advance influence range far exceeding that of conventional seams, a characteristic that becomes more pronounced as mining progresses.
In contrast, the rear influence distance remained smaller and relatively stable, generally within 200 m, suggesting that after stress release in the goaf, the continued disturbance on the floor is weak. Therefore, MS events are primarily dominated by stress concentration in the advance zone.
Vertical distribution characteristics
Fig. 10 Longitudinal section distribution of MS events along the roadway.
The Carboniferous Taiyuan Formation (colored in blue) exhibits the highest number and widest distribution of MS events, representing the primary layer of floor MS activity. The large mining height and associated stress from thick coal seam extraction directly act on this layer, causing frequent rock fracturing and the strongest MS responses.
In the Carboniferous Benxi Formation (yellow), MS events are significantly reduced and distributed sparsely beneath the Taiyuan Formation. The vertical transmission of mining-induced stress to this layer is weaker than in the Taiyuan Formation, resulting in a lower frequency of rock fracturing.
The Ordovician Formation (red) exhibits the fewest MS events, which are sporadic and scattered. The high rock strength and greater distance from the coal seam floor, coupled with the influence of mining-induced stress and floor water pressure, make this layer less susceptible to large-scale fracturing.
Fig. 11 Statistical chart of floor MS events development depth.
MS events are highly concentrated in the shallow region beneath the coal seam floor, with the 0–28 m depth interval exhibiting peak frequency, accounting for approximately 96% of total events. This interval is directly below the coal seam floor and is subjected to mining-induced stress, resulting in rock failure primarily through tensile fracturing and shear slip, which represents the core zone of mining energy release and the strongest MS responses.
With increasing depth, MS events frequency decays exponentially. The 26–28 m interval accounts for about 0.59%, and the 28–30 m interval accounts for only 0.31%, indicating very few MS events beyond 28 m. Therefore, the concentrated floor failure zone based on MS monitoring is considered to extend up to 28 m below the coal seam floor.
The MS monitoring analysis results show a maximum deviation of less than 5% from the theoretical empirical calculation, indicating that the floor failure depth derived from MS monitoring is reasonable and reliable.
Energy distribution characteristics
Fig. 12 Daily total energy of floor MS events and variation of maximum hydraulic support convergence.
The daily cumulative MS energy exhibits an overall “low–high–low–high” trend. During the early and mid-mining stages (April–June), daily cumulative energy remained at a relatively low level with small fluctuations. During this period, the working face was in the “initial stress adjustment” stage, and mining-induced disturbances to the floor were relatively mild. Rock fracturing was limited in scale and frequency, the maximum subsidence of hydraulic supports was stable, mining-induced stress on the floor was low, and water inflow remained stable without significant increase.
In the late mining stage (July–August), daily cumulative energy experienced a sudden surge around mid-July, coinciding with a corresponding increase in maximum hydraulic support subsidence. This period was accompanied by strong manifestation of mining-induced stress, a drop in limestone water level, and an increase in water inflow, indicating that as mining progressed, stress accumulation and release in the floor entered a high-disturbance stage. Large-scale, high-energy rock fracturing occurred, and MS responses intensified significantly. In late August, as the working face entered the “secondary square” stage, both daily cumulative energy and mining-induced stress increased again.
Fig. 13 Energy distribution of ms events in the benxi formation and deeper strata: (a) June; (b) July; (c) August.
The temporal evolution of floor MS energy exhibits a characteristic pattern of “early low-energy stability, late high-energy burst”, which is closely coupled with mining-induced stress manifestation and water inflow changes. High-energy MS events are the key factor driving the extension of floor failure to greater depths and provide a direct basis for establishing energy-threshold-based early warning for floor water inrush risk.
Floor water inrush hazard evaluation based on AHP-EWM
Indicators for floor water inrush hazard
Based on the spatiotemporal characteristics of floor MS events analyzed in Section "Spatiotemporal response characteristics of floor microseismic events", the MS-based indicators directly reflecting the floor water inrush hazard can be classified into three categories: activity indicators, spatiotemporal indicators, and intensity indicators.
(1) Activity Indicators. As discussed in section "Temporal distribution characteristics", the activity of floor MS events reflects the degree of rock fracturing and floor disturbance. The selected activity indicators include: Daily frequency of floor MS events and daily frequency of MS events in the Benxi Formation and deeper strata.
(2) Spatiotemporal Indicators. The formation of water-conducting channels is a necessary but not sufficient condition for water inrush. The spatiotemporal distribution of MS events can reflect the development and connectivity of these channels. Therefore, the selected spatiotemporal indicators include: Maximum daily development depth of floor MS events and advance influence distance of MS events.
The Benxi Formation serves as a key floor aquiclude, and when the maximum daily development depth of floor MS events exceeds the Benxi Formation, it indicates an increased risk of water inflow, potentially leading to enhanced water discharge or even floor water inrush. The dynamic trend of this indicator can be directly linked to the evolution of mining-induced water-conducting fracture channels.
Fig. 14 Point-line chart of the daily maximum development depth of floor MS events.
(3) Intensity Indicators. As analyzed in Section "Energy Distribution Characteristics", the energy of floor MS events reflects high-energy, large-scale fracturing and its impact on mining-induced stress manifestation, illustrating that MS events energy is closely associated with mining stress and the development of water-conducting channels. Therefore, the daily cumulative energy of floor MS events is selected as the intensity indicator.
Based on the above analysis, five MS events indicators were identified to evaluate the floor water inrush hazard, reflecting the activity, spatiotemporal, and intensity characteristics of floor MS events: daily frequency of floor MS events, daily frequency of MS events in the Benxi Formation and deeper strata, maximum daily development depth of floor MS events, advance influence distance of MS events, and daily cumulative energy of floor MS events.
Fig. 15 The evaluation indicator system.
Weight calculation of floor water inrush hazard indicators
Table 2 Judgment matrix for first-level indicators.
| Indicators | B1 | B2 | B3 |
|---|---|---|---|
| B1 | 1 | 1.15 | 3.90 |
| B2 | 1/1.15 | 1 | 3.38 |
| B3 | 1/3.90 | 1/3.38 | 1 |
Table 3 Judgment matrix for activity indicators.
| Indicators | C1 | C2 |
|---|---|---|
| C1 | 1 | 1/5.58 |
| C2 | 5.58 | 1 |
Table 4 Judgment matrix for spatiotemporal indicators.
| Indicators | C3 | C4 |
|---|---|---|
| C3 | 1 | 1.55 |
| C4 | 1/1.55 | 1 |
Table 5 EWM calculation results.
| Indicators | Information Entropy | Information Utility Value | ωEWM |
|---|---|---|---|
| Daily frequency of floor MS events | 0.9126 | 0.0874 | 0.1409 |
| Daily frequency of MS events in the Benxi Formation and deeper strata | 0.7944 | 0.2056 | 0.3315 |
| Maximum depth of floor MS activity | 0.8464 | 0.1536 | 0.2476 |
| Advance influence distance of MS events | 0.9538 | 0.0462 | 0.0745 |
| Daily cumulative energy of floor MS events | 0.8725 | 0.1275 | 0.2056 |
Table 6 Table of weight calculation results.
| Indicators | ωAHP | ωEWM | ωAHP×ωEWM | ω |
|---|---|---|---|---|
| Daily frequency of floor MS events | 0.0716 | 0.1409 | 0.010145 | 0.049015 |
| Daily frequency of MS events in the Benxi Formation and deeper strata | 0.3998 | 0.3315 | 0.107738 | 0.520538 |
| Maximum depth of floor MS activity | 0.2478 | 0.2476 | 0.023691 | 0.114464 |
| Advance influence distance of MS events | 0.1597 | 0.0745 | 0.040111 | 0.193799 |
| Daily cumulative energy of floor MS events | 0.1210 | 0.2056 | 0.025289 | 0.122184 |
According to the composite weights, the ranking of the indicators from highest to lowest importance is as follows: daily frequency of MS events in the Benxi Formation and deeper strata > advance influence distance of MS events > daily cumulative energy of floor MS events > maximum development depth of floor MS events > daily frequency of floor MS events.
Comprehensive evaluation of floor water inrush hazard
The daily comprehensive hazard index was calculated and compared with the observed water inflow data of the working face, as shown in Figure 16. The comprehensive index was classified into risk levels using equal interval division: low risk (0–0.3), medium risk (0.3–0.6) and high risk (0.6–0.9).
Fig. 16 Comprehensive evaluation of floor water inrush risk.
Fig. 18 Linear regression fitting.
Fig. 17 Verification of the comprehensive evaluation of floor water inrush risk (February–June 2024).
Discussion and outlook
Stress transfer characteristics of floor fracturing.
(2) Dynamic evaluation using MS monitoring
(3) Limitations and future research directions
The current MS monitoring dataset covers a limited temporal span, excluding the full mining cycle (e.g., initial and periodic stress manifestation stages), which may introduce bias in weight calculations. Moreover, the selected evaluation indicators do not consider multi-factor coupling effects, such as rock mechanical properties or grouting reinforcement coverage, potentially causing deviations between predicted and actual water inrush risk under complex geological conditions.
Future research could focus on:
Extending monitoring duration and expanding the indicator system: Cover the full mining cycle and incorporate factors such as floor rock integrity coefficients and grouting reinforcement effectiveness, constructing a “dynamic MS – static geological – engineering intervention” multidimensional evaluation system to improve comprehensiveness.
Mechanism study with numerical simulation: Develop a floor stress-fracture numerical model for thick coal seam mining, and compare the simulated fracture depth and energy distribution with MS monitoring data to reveal the “mining-induced stress – rock fracturing – hydraulic lift” coupled disaster thresholds, providing a theoretical basis for weight optimization.
Development of intelligent early-warning systems: Based on the AHP-EWM evaluation model and real-time MS monitoring data, develop an integrated system of data acquisition → weight calculation → risk classification → early-warning issuance, enabling dynamic and automated water inrush risk prediction, thereby supporting safe and efficient mining of thick coal seams.
Conclusions
Based on continuous MS monitoring of the 61607 working face in the Longwanggou Coal Mine, this study systematically characterizes the spatiotemporal response of floor MS activity induced by extremely thick coal-seam mining. The results show that MS events are mainly concentrated within the Taiyuan Formation and exhibit coordinated increases in frequency, depth, spatial extent, and energy release as the working face advances. In combination with observed inflow rates and mine pressure data, MS indicators such as the extent of the advance-affected zone, floor damage depth, and daily cumulative energy show clear responses, revealing the full spatiotemporal evolution of floor failure.
The characteristics of the recorded MS events indicate that floor microseismicity induced by extremely thick coal-seam mining exhibits the features of “deep occurrence–wide spatial distribution–high energy release.” Deep MS events within the Benxi Formation and underlying strata are strongly associated with aquifer disturbance, structural reactivation, and groundwater movement. These events serve as critical indicators for identifying deep water-related hazards and significantly enhance the accuracy of floor water-inrush early warning.
By summarizing the key controlling factors of floor water-inrush risk, a dynamic MS indicator system—incorporating event activity, spatiotemporal distribution, and intensity characteristics—was established. Using the AHP–EWM combined weighting method, a quantitative evaluation model for floor water-inrush risk was constructed. The evaluation results show strong consistency with field monitoring data, confirming the applicability and effectiveness of the model under extremely thick coal-seam mining conditions.
Publisher’s note
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
Acknowledgements
The author gratefully acknowledges the funding and assistance from proof-readers and editors.
Author contributions
Haorui Wang conceived the research and wrote the initial draft of the paper, Shangxian Yin and Huiqing Lian provided guidance for the on-site experiments, Min Cao was responsible for data processing and chart drawing, Enke Hou and Xu Wang proofread and translated the paper, and Guofu Xing provided the experimental site. All authors reviewed, revised, and approved the final version of the manuscript.
Funding
National Key Research and Development Program of China (No. 2022YFC3005905-1); Project under the Central-Guided Fund for Local Science and Technology Development (No. 246Z7604G); Hebei Natural Science Foundation Project (No. D2025508011); The Fundamental Research Funds for the Central Universities (No. 3142025007).
Data availability
The authors declare that all data presented in this manuscript are available to researchers for noncommercial use, provided participant confidentiality is maintained. Data requests can be directed to the corresponding author.
Declarations
Competing interests
The authors declare no competing interests.
Informed consent
All experts involved in the experimental part of this study were fully informed of the research objectives, procedures, and intended use of the collected data. Their participation was entirely voluntary, and each provided explicit informed consent before taking part in the study.
References
- 1. TianGShangXLiJWangYA strategic study of coal exploitation in Shanxi, Shannxi and Inner MongoliaCoal Geol. China200820115
- 2. ZhangZet al.Overburden destruction structure model and mechanism of mine pressure in Caojiatan coal mine with 10 m super large mining heightJ. China Coal Soc.20255019511964
- 3. PanCXiaBZuoYYuBOuCMechanism and control technology of strong ground pressure behaviour induced by high-position hard roofs in extra-thick coal seam miningInt. J. Min. Sci. Technol.20223249951110.1016/j.ijmst.2022.01.006
- 4. FengDet al.Research on overburden structural characteristics and support adaptability in cooperative mining of sectional coal pillar and bottom coal seamSci. Rep2024141145810.1038/s41598-024-62375-738769414PMC11106068
- 5. LinHet al.Study on the energy evolution mechanism and fractal characteristics of coal failure under dynamic loadingACS Omega.202510547105471910.1021/acsomega.5c0816241280835PMC12631684
- 6. CaoZet al.Disaster-causing mechanism of spalling rock burst based on folding catastrophe model in coal mineRock Mech. Rock Eng.202558710.1007/s00603-025-04497-6
- 7. DongSet al.Theoretical framework and key technology of advance regional control of water inrush in coal seam floorCoal Geol. Explor.202351185195
- 8. Yin, S. et al. Mechanism and full-time-space prevention and control technology of water inrush from Ordovician and thin limestone in deep mines. J. China Coal Soc.45 (2020).
- 9. ZhengGXuXHeYRenCLiuBMechanism of Ordovician limestone water inrush from the floor of caving mining face in thick coal seam and water hazard detection technology of nearly horizontal boreholeCoal Geol. Explor.201947713
- 10. Shao, Q. H. Study on roof and floor failure law of fully mechanized caving mining under the extra-thick coal seam—A case study of Longwanggou Coal Mine. Master’s thesis Anhui Univ. Sci. Technol. (2022) (in Chinese)
- 11. LiHTangSBKangZQYangDMaLQCFDEM simulation on the linkage effect between floor fault failure and roof collapse in ultra-thick coal seamJ. China Coal Soc.20244926152629
- 12. Wang, S. Study on the evolution law of mining-induced fractures in the floor of extra-thick coal seam fully mechanized top-coal caving mining. Master’s thesis Henan Polytech. Univ. (2024) (in Chinese)
- 13. LiAMaQLianYMaLMuQChenJNumerical simulation and experimental study on floor failure mechanism of typical working face in thick coal seam in Chenghe mining area of WeibeiChina. Environ. Earth Sci.20207911810.1007/s12665-020-8839-2
- 14. HeTLiGDSunCLuoFLiXQFloor failure characteristics of thick coal seam mining above confined aquiferMin. Metall. Explor.20223915531562
- 15. ZhangPOuYSunBLiuCA case study of floor failure characteristics under fully mechanised caving mining conditions in extra-thick coal seamsJ. Geophys. Eng.20201781382610.1093/jge/gxaa031
- 16. LiuCet al.Comprehensive research on the failure evolution of the floor in upper mining of deep and thick coal seamJ. Appl. Geophys.202220610477410.1016/j.jappgeo.2022.104774
- 17. ZhouWZhangPWuRHuXDynamic monitoring the deformation and failure of extra-thick coal seam floor in deep miningJ. Appl. Geophys.201916313213810.1016/j.jappgeo.2019.02.007
- 18. ZhangPLiuCOuYSunBXuSLiSComprehensive testing research on floor damage characteristics of mining extra-thick seam in Jungar CoalfieldCoal Geol. Explor.20214930
- 19. GuoGFloor failure law of extra-thick coal seam in fully mechanized caving miningCoal Geol. Explor.20225013
- 20. Cao, M. et al. Failure mechanism and risk evaluation of water inrush in floor of extra-thick coal seam. Water17 (2025).
- 21. ZhuZWangPLiRQinHShiQChenKMechanism and control technology of floor heave in fully mechanized top-coal caving faces of ultra-thick coal seamsChina Saf. Sci. J.20253573
- 22. ZhangZHuangZLiuXLiuQZhaoTLinXCharacteristics and mechanism of strong dynamic load mine pressure in large mining height full-mechanized caving face with multiple thick and hard roofJ. China Coal Soc.20255033693387
- 23. LiJZhaoBResearch on strong mine pressure behavior and control technology in an extra-thick coal seam working faceCoal Eng.2025577581
- 24. KongXet al.Fracture-seepage characteristic zones of the overburden of extra-thick coal seams: evolutionary patterns and their engineering applicationsCoal Geol. Explor.202553102113
- 25. MengHLChengXGQiaoWResearch on the evolution and prediction of heights of water-conducting fracture zones in overlying rocks during layered mining of extremely thick coal seamsInd. Mine Autom.2024506775
- 26. YuBet al.Disaster mechanism and collaborative control technology of multi-field coupling in extra-thick coal seam with hard roofCoal Sci. Technol.202553116
- 27. LiNWangEYGeMCMicroseismic monitoring technique and its applications at coal mines: present status and future prospectsJ. China Coal Soc.2017428396
- 28. YinHet al.A deep learning-based data-driven approach for predicting mining water inrush from coal seam floor using microseismic monitoring dataIEEE Trans. Geosci. Remote Sens.20236111510.1109/TGRS.2023.3331904
- 29. XinCJiangFZhaiCChenYAnalysis of coal floor fault activation inducing water inrush using microseismic monitoring—a case study in Zhaogu No. 1 coal mine of Henan Province, ChinaSustainability202315736110.3390/su15097361
- 30. JiaJZhouJNiuQDynamic evolution characteristics and risk assessment of floor microseismic in close repeated miningCoal Sci. Technol.2025534657
- 31. HuangLet al.Research on the development law of pre-mining microseisms and risk assessment of floor water inrush: a case study of the Wutongzhuang coal mine in ChinaSustainability202214977410.3390/su14159774
- 32. Gai, Q., Gao, Y., Zhang, X. & He, M. A new method for evaluating floor spatial failure characteristics and water inrush risk based on microseismic monitoring. Rock Mech. Rock Eng.57 (2024).
- 33. GaiQGaoYHuangLShenXLiYMicroseismic response difference and failure analysis of roof and floor strata under dynamic load impactEng. Fail. Anal.202314310687410.1016/j.engfailanal.2022.106874
- 34. XuYHuangLFull-time and space early-warning method for floor water inrush in working face based on microseismic monitoringCoal Sci. Technol.202351369382
- 35. Huang, L. Study on full-time and space prediction and evaluation for floor water inrush at the working face based on microseismic monitoring. Doctoral dissertation China Univ. Mining Technol. (Beijing) (2023) (in Chinese)
- 36. ZuoJWuGDuJLeiBLiYRock strata failure behavior of deep Ordovician limestone aquifer and multi-level control technology of water inrush based on microseismic monitoring and numerical methodsRock Mech. Rock Eng.2022554591461410.1007/s00603-022-02891-y
- 37. MaKSunXYTangCAYuanFZWangSJChenTFloor water inrush analysis based on mechanical failure characters and microseismic monitoringTunn. Undergr. Space Technol.202110810369810.1016/j.tust.2020.103698
- 38. LiuRZhiGYangSXuXInvestigation on mining-induced floor water inrush from column and its control based on microseismic monitoringAdv. Civ. Eng.202320233754079
- 39. NieBSYuanSFWangLKZhangJLSuXQRoadway roof stability analysis based on microseismic monitoring technologyCoal Mine Saf.201647178180
- 40. Su, Z. Z. Research on the monitoring and early-warning system of mine water in Jining Mine. Master’s thesis North China Inst. Sci. Technol. (2021) (in Chinese)
- 41. BaiJJLüYLComparative analysis of two kinds of microseismic monitoring systems under collaborative verification modeShandong Coal Sci. Technol.202011187189
- 42. YangQJWangDPHeQFResearch and application of microseismic monitoring technology for large dip angle coal seamEnergy Conserv.202211160164
- 43. ChangJEJiangTLResearch on the weight of coefficient through analytic hierarchy processJ. Wuhan Univ. Technol.200729153166
- 44. DengXLiJMZengHJChenJYZhaoJFAnalysis and application of weight calculation methods in analytic hierarchy processMath. Pract. Theory20124293100
- 45. IshizakaALabibAReview of the main developments in the analytic hierarchy processExpert Syst. Appl.2011381433614345
- 46. YangLLiuCCSongLEvaluation of coal mine emergency rescue capability based on entropy weight methodChina Soft Sci.201311185192
- 47. ZhuYTianDYanFEffectiveness of entropy weight method in decision-makingMath. Probl. Eng.202020203564835
- 48. LiJLZhangHYWangXYFengYLLiuZHanLApplication and suggestion on vulnerable index method of coal seam floor water burst evaluationJ. China Coal Soc.201439725730
- 49. WuQZhangSA new practical methodology of the coal floor water bursting evaluating II—the vulnerable index methodJ. China Coal Soc.20073211211126
- 50. WuQLiuYLiuDZhouWPrediction of floor water inrush: the application of GIS-based AHP vulnerable index method to Donghuantuo coal mineChina. Rock Mech. Rock Eng.20114459160010.1007/s00603-011-0146-5
- 51. QiangWLiYLiuYUsing the vulnerable index method to assess the likelihood of a water inrush through the floor of a multi-seam coal mine in ChinaMine Water Environ.201130546010.1007/s10230-010-0125-1
- 52. ZengYWuQLiuSZhaiYZhangWLiuYVulnerability assessment of water bursting from Ordovician limestone into coal mines of ChinaEnviron. Earth Sci.201675143110.1007/s12665-016-6239-4
- 53. MengHet al.Mechanism and criterion of water burst in thick plate zone with pressure miningChina Coal Soc.202550S1250262
Copyright and License
© 2025