Next Article in Journal
Using the Life Cycle Approach for Multiobjective Optimization in the Context of the Green Supply Chain: A Case Study of Brazilian Coffee
Next Article in Special Issue
A Numerical Method for Evaluating the Collapse of High-Steep Scarp Slopes Based on the Bonded Block Model–Discrete Fracture Network Model
Previous Article in Journal
Unveiling Job Satisfaction of Teachers through a Blend of Methodologies
Previous Article in Special Issue
Experimental Study on Mechanical Properties and Acoustic Emission Characteristics of Dry and Water-Saturated Soft Rocks under Different Dynamic Loadings
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Prediction of the Water Inrush Risk from an Overlying Separation Layer in the Thick Overburden of a Thick Coal Seam

1
College of Earth Science and Engineering, Shandong University of Science and Technology, Qingdao 266590, China
2
College of Energy and Mining Engineering, Shandong University of Science and Technology, Qingdao 266590, China
3
Yankuang Energy (Ordos) Co., Ltd., Ordos 017010, China
4
Shandong Energy Group Luxi Mining Co., Ltd., Heze 274700, China
*
Author to whom correspondence should be addressed.
Sustainability 2023, 15(18), 13988; https://doi.org/10.3390/su151813988
Submission received: 2 August 2023 / Revised: 13 September 2023 / Accepted: 18 September 2023 / Published: 20 September 2023
(This article belongs to the Special Issue Deep Mining Engineering in Sustainability)

Abstract

:
With the expansion of coal mining westward in China, water inrush from seam roofs has become a prominent safety problem during mining. The roof rock of the coal seam in the Shilawusu coal mine has the characteristics of a double-layer structure, and the overlying separation space formed in the mining process of the coal seam poses a risk of water inrush. To ensure the safety of coal mine production, considering the geological and hydrogeological data of the mining area, the core recovery rate, lithologic assemblage index, key aquifer thickness, hydrostatic head and lithologic structure index of the Zhidan Formation are selected as evaluation indexes. The index weights are calculated based on the attribute hierarchical model and coefficient of variation methods, and subjective and objective preference coefficients are introduced to determine the ranking of comprehensive indexes. The catastrophe progression method is improved, and a zoning prediction model for water inrush risk is established by the improved catastrophe progression method. The results show that only a tiny part of the mining area is in danger, and most areas are in the safe and transition zones. The model realizes the prediction of the risk of water inrush from the overlying separation layer in the study area and provides a theoretical basis for the prevention and control of water inrush from the overlying separation layer in coal mining.

1. Introduction

The Shaanxi and Mongolia region in Western China is an important coal base in China and plays a pivotal role in energy cooperation along the “Belt and Road Initiative” route and in the construction of national energy transportation channels. However, in 2006, secondary water from bed separation occurred in the 745 workfaces of Haizi coal mine, with an inflow of 3887 m3/h, resulting in the death of five workers [1]. On 25 April 2016, the ZF202 working face of Zhaojin coal mine in Tongchuan caused many deaths due to the increasing volume of space between layers of overlying strata and water that had accumulated in the space [2]. When the mining face in the Hongliu coal mine advanced different distances from the initial cut, the maximum water inflow increased to 3000 m3/h [3]. Therefore, the threat of water inrush from an overlying separation layer has been given increasingly more attention in mining engineering. It is significant to accurately predict the risk of water inrush from the overlying separation layer for the safety of coal mine production. Due to the excellent occurrence state of the coal seams in the northern Ordos Basin, the coal seam roof can be described as a double-layer structure with a coal-measure composite aquifer and Cretaceous Zhidan Formation aquifer, which indicates that the coal mining in this region will face the threat of water inrush from separated layers. As a newly identified type of roof water hazard, water inrush from layer separation can influence the average production of coal mines and the safety of personnel underground. Therefore, it is significant to accurately predict the risk of water inrush from an overlying separation layer for coal mine safety and production.
At present, many experts and scholars focus on the location of the separation layer water damage and the mechanism of water inrush from the separation layer and have achieved substantial results. Qiao et al. [4] summarized the main theories to explain and judge the development of overburden separation in China, including key layer theory, arch beam balance theory, “upper four belts” theory and elastic thin plate theory. Gui et al. [5] analyzed the key influencing factors of the formation of a separation layer and the occurrence of water disasters and provided a reference for further research on water inrush and disaster prevention in coal mines. Yan et al. [6] proposed a new method to predict the location of roof strata separation and derived the calculation formula of bedding separation in roof strata by combining the theoretical analysis and test results of continuous and discontinuous beams and layered strata loads. He et al. [7] proposed an improved stepwise comparative combination method based on composite beam theory to determine the location of a separation layer. According to key strata theory, Hu et al. [8] determined the position of separated beds in overlying rocks and introduced the trapezoidal platform model of the cracks in overlying strata. Ji et al. [9] studied the mechanism of four large water inrush events in the mining process of the Hongliu coal mine 1121 workface in the Ningdong mining area. Wu et al. [10] used the rock plate method to study the deformation and internal force of coal seam roofs and proposed an identification method for isolated water damage. Gui et al. [11] studied the restriction effect of geological and mining factors such as roof rock lithology, coal seam dip angle, separation layer water source, roof management mode, working face size and mining schedule on the formation of water damage to the separation layer. Numerical simulation and physical simulation also play an important role in the study of separation layers. Li et al. [3] studied the evolution of water inrush through numerical simulation and underground full-space transient electromagnetic detection and concluded that water inrush from a separation layer is due to unique geological conditions, lithology and mining conditions. Gui et al. [12] determined the location of the separation layers and the location of the most dangerous separation layer of the overlying strata through numerical simulation and similar material simulation. Wang et al. [13] developed a new type of similar material to simulate the separation of rock layers and observe the water flow in the waterproof layer. Wang et al. [14] experimentally revealed the evolution of layer separation, the process of water accumulation due to layer separation and the damage characteristics of water inrush due to layer separation. By means of numerical simulation and analysis, Li et al. [15] studied the influencing factors and mechanism of water inrush from the separation layer and proposed effective measures for the prevention and control of separation layer water inrush. In addition, the risk assessment system of water inrush of coal seams tends to be mature in China [16,17,18,19,20,21,22].
The above work has carried out a series of theoretical studies on the formation of the conventional separation layer and the mechanism of water inrush from the separation layer. However, coal seam occurrence conditions are relatively good in China’s northern mining area of the Ordos Basin. The coal seam roof can be divided into a double-layer structure, which includes a coal-measure composite aquifer and the Cretaceous Zhidan Formation aquifer. There are few studies on the prediction of water inrush of high separation layer water, which seriously threatens the safety of mine mining. Therefore, this paper establishes a more comprehensive and objective risk prediction model that combines subjective and objective weights with the catastrophe progression method. The risk of water inrush from the separation layer in the 221 mining area of Shilawusu coal mine is predicted, and the evaluation results have important guiding significance for the mining area in the north of Ordos Basin.

2. Materials and Methods

2.1. Geological Setting

The Dongsheng coal field, Shilawusu coal mine, is located in Ordos city, and the administrative division is under the jurisdiction of the Ejin Horo Banner of Ordos city. The length of the coal mine is 7.35 km from north to south and 9.40 km from east to west on average, covering an area of 70.644 km2 (Figure 1). The coal seam is adjacent to the Hulusu well field to the west and the No. 2 exploration area to the south. Strata in the area include, from old to young, the upper Triassic Yanchang Formation (T3y), middle–lower Jurassic Yan’an Formation (J1–2y), middle Jurassic Zhiluo Formation (J2z), middle Jurassic Anding Formation (J2a), lower Cretaceous Zhidan Formation (K1zh), upper Quaternary Pleurian Malan Formation (Q3m), Quaternary Holocene Aeolian Layer (Q4eol) and alluvial sand (Q4al+pl); see Table 1 for a brief list of stratum characteristics in the mining area.
This paper selects the 221 mining area, which is located in the middle of the Shilawusu coal mine, as the study area. The 221 mining area is approximately 7.4 km from north to south and 5.1 km from east to west, covering an area of 37.74 km2. The main coal-bearing strata are the Yan’an Formation (J1–2y) of the middle Lower Jurassic, with no surface outcrop and stable coal seam deposition. The main coal seam is 2-2 coal seam, and the coal thickness is approximately 9.18 m. The base of the coal seam is the Yanchang Formation of the upper Triassic system.
In the study area, the middle and upper parts of the Yan’an Formation, Zhiluo Formation and Anding Formation are interbedded with sandy mudstone and sandstone (59.3% sandstone on average). Coal measures are dominated by aquifers. The rock beam combination is mainly composed of coarse sandstone, medium sandstone and fine sandstone, and the rocking beam overlying rock is generally a secondary-grained rock layer. The rock beam combination is relatively complex, and the mechanical properties of rock strata in this group are more obvious than those in the Zhidan Formation. The overall water richness is weak.
Here, the Cretaceous system is mainly the lower Cretaceous Zhidan Formation strata. The lithology is dominated by fine sandstone, medium sandstone and coarse sandstone (90% sandstone), and sandy mudstone and mudstone are rarely developed (Table 2). The rock beam is mainly composed of thick, gravelly strata and the combined structure of the rock beam is more straightforward compared to the coal-measure composite aquifer. The brittleness characteristic of mechanical properties is relatively apparent, and the overall water richness is high.
Based on the above differences above the roof of the coal seam in the study area, the whole stratum in the study area is divided into a double-layer structure (Figure 2), that is, from the bottom to top of the 2–2 coal roof of the main coal seam, two combined water-bearing rock groups can be distinguished: the first layer of the coal-measure composite aquifer and the second layer of the Cretaceous Zhidan Formation aquifer.

2.2. Evaluation Model and Weight Determination Method

As the workface advances, different mechanical properties, different degrees of softness and hardness of the rock layer and the interlayer cohesion will cause the lower rock layer under the action of self–weight on the upper rock layer to produce tensile force. When the tensile force reaches the limit of interlayer tensile strength, the upper and lower layers of the rock layer is separated, resulting in the occurrence of varying degrees of subsidence and their independent bending and deformation [23]. When the tensile force reaches the ultimate tensile strength between layers, the upper and lower rock layers are separated, different degrees of settlement occur, and they independently bend. When the deformation reaches a certain degree, a separation space is formed between the soft and hard rock layers. Water inrush is likely to cause a disaster when the following four conditions are met: a closed and sustainable separation space, a water replenishment channel, a stable water source and a water inrush channel.
Under the condition of double-layer overlying rock, the mining of deep and thick coal seams results in the development of large–scale layer separation between the double-layer overlying coal measures and the Cretaceous series on the basis of the traditional “upper three zones”. With the continuous mining of the lower coal seam, the range of the overlying separation continues to expand. However, due to the weak cementation and brittle deformation characteristics of the upper overburden of the separated roof, new fracture movement occurs within a certain range of the overburden, namely the secondary movement (Figure 3). The intensity and scale of the movement of the extremely thick sandstone layer in the overlying Zhidan Formation are both large due to secondary movement. The deformation is no longer the fracture and bending of the rock layer in a single layer or small range, but is instead the fracture of a group of strata. The bottom rock layer of the Zhidan Formation is broken, the upper rock layer is curved, and a separation space is generated. When the aquiclude at the lower part of the separation layer breaks under the double action of mine pressure and water pressure, the water in the separation layer will flow into the working face and cause a water inrush accident (Figure 4).
The possibility of water inrush is closely related to the key aquiclude and the properties of the underlying barrier layer for the overlying separation layer containing confined water. Therefore, in view of the double-layer structure characteristics of overburden strata of the 2-2 coal seam in the 221 mining area, this paper selects five indexes, namely, core recovery rate of the key aquiclude, lithologic assemblage index of the key aquiclude, thickness of the key aquiclude, static pressure head and lithologic structure index of the Zhidan Formation, as the basis for studying water inrush from an overlying separation layer of a coal seam roof.

2.2.1. Core Recovery Rate of the Key Aquiclude (C1)

The core recovery rate refers to the ratio of the length of the core recovered to the length drilled in the drilling process, which can reflect the fragmentation degree of the stratum to a certain extent [24]. The larger the core recovery rate is, the more complete the rock mass, the less developed the cracks and the stronger the water insulation performance. In contrast, the water insulation performance is weaker, the water permeability and conductivity are stronger and the corresponding water insulation performance is worse for smaller core recovery rates. Therefore, the core recovery rate of the key aquiclude is negatively correlated with the risk of water damage from upper layers (Figure 5).

2.2.2. Lithologic Assemblage Index of Key Aquiclude (C2)

A large number of studies have shown that the unit water inflow in lithologic assemblages is generally larger if the proportion of sandstone is larger, while the unit water inflow in lithologic assemblages is smaller if the proportion of sandstone is smaller. It is generally believed that as the proportion of sandstone decreases, the porosity will gradually decrease, so the capacity of water storage and water conduction will gradually decrease, and then the aquiclude performance will gradually strengthen. Therefore, different structural coefficients can be assigned to different rock types (Table 3). The larger the structural coefficient is, the larger the particle size of the rock, and thus the poorer the barrier performance of the rock. Hou et al. [25] proposed the lithologic assemblage of key aquifers as an evaluation index of waterproofing. Therefore, in this paper, the thicknesses and structural coefficients of different lithologies in the key aquiclude are used to construct the lithologic assemblage index z, reflecting the key aquiclude to evaluate the key aquiclude performance (Equation (1)).
Table 3. Structural coefficient table.
Table 3. Structural coefficient table.
Rock TypesCoarse SandstoneMedium SandstoneFine SandstoneSiltstoneSand MudstoneMudstone
Structural coefficient10.80.60.40.30.2
z = h i s i h i i = 1 , 2 , 6
where hi represents the thickness of coarse sandstone, medium sandstone, fine sandstone, siltstone, sandy mudstone and mudstone for i = 1, 2, 3, 4, 5 and 6, respectively, and Si is the corresponding structural coefficient.
The larger the z value is, the larger the proportion of sandstone in the key waterproof layer, the better the water storage and water conductivity, and the lower the aquiclude performance; in contrast, the larger the z value is, the stronger the water resistance. Therefore, the lithologic assemblage index of the key aquiclude is positively correlated with the risk of water damage from the upper layers.

2.2.3. Key Aquiclude Thickness (C3)

The thickness of the key aquiclude is an important factor affecting water inrush from overlying separation layers. The greater the thickness of the key aquiclude, the stronger the ability to resist water pressure in overlying separation layers, the lower the possibility of connecting a water channel below the separation layers to the lower key separation layer and the lower the risk of water inrush in the overlying separation layers. Therefore, the thickness of the key aquiclude is negatively correlated with the risk of water damage from upper layers.

2.2.4. Lithologic Structure Index of the Zhidan Formation (C4)

The lithologic structure index is used to multiply the thickness of the aquifer in the Zhidan Formation, including the thicknesses of the medium sandstone and fine sandstone by an equivalent coefficient, and then convert it into the thickness of coarse sandstone. Considering the number of alternating sandstone–mudstone layers, the more mudstone–mudstone (siltstone) interlayer combinations there are, the weaker the water-carrying performance, so the number of sandstone–mudstone interlayers is negatively correlated with S when other conditions remain unchanged.
The water-richness of the Zhidan Formation sandstone aquifer, as the main aquifer of the overlying separation layer cavity, will affect the outburst risk of overlying separation layer water. The greater the water-richness, the greater the water volume and the greater the outburst risk of overlying separation layer water. In contrast, the lower the water-richness, the lower the water outburst risk of overlying separation layer water. Therefore, the lithologic structure index of the Zhidan Formation is positively correlated with the risk of water damage in the upper strata.
Based on the above analysis, the expression of the lithologic structure index was constructed as follows:
S = ( a + b × i + c × j ) n
where a, b and c are the thicknesses of the coarse sandstone, medium sandstone and fine sandstone, respectively; m, and i and j are the equivalent coefficients of the medium sandstone and fine sandstone, respectively; and n is the number of sandstone–mudstone (siltstone) interbeds.

2.2.5. Hydrostatic Head (C5)

There are many unexploited areas in the mining area, so the highest water level of the sandstone aquifer of the Zhidan Formation in the mining area is used as its head elevation. The higher the water head in the overlying separation area, the greater the water pressure, and the more likely the water is to surge into the lower channel fracture zone, resulting in the risk of water damage; otherwise, the risk of water inrush is smaller. Therefore, the hydrostatic water head of the sandstone aquifer in the Zhidan Formation is positively correlated with the risk of water damage to upper layers. The elevation of the sandstone aquifer floor in the Zhidan Formation was calculated using the borehole data in the mining area, and then the static pressure head at the overlying separation space was calculated.

2.2.6. Data Statistics

According to risk factor analysis, 47 boreholes in the 221 mining area were selected for statistical analysis, and the original data are shown in Table 4.

2.3. Attribute Hierarchical Model (AHM)

As a subjective weighting method, the AHM is an improvement of the AHP. There is no need to calculate the feature vector and perform a consistency test, thus avoiding the shortcomings of the AHP, and so, an AHM is a more practical subjective analysis method. The decision-making steps using AHM are divided into three steps: (1) building a recursive hierarchy, (2) constructing a judgment matrix and calculating the relative weights, and (3) calculating the synthetic weights of the scenarios concerning the system’s objectives to make a decision.
The hierarchical model of water inrush risk is constructed from an overlying separation layer, as shown in Figure 6. The flow chart of the evaluation process is shown in Figure 7.
The Saaty scaling method was adopted to obtain the n-order discriminant matrix R = (rij)n×n via the expert scoring method, in which Rij represents the importance of factor i compared with factor j. The discriminant matrix R = (rij)n×n has the following properties:
R i j > 0 R i j = 0 R i j = 1 k i j
where ij, 1 ≤ in and ijn.
The attribute discrimination matrix was constructed, which can be expressed as follows:
k i j = 2 m 2 m + 1 R i j = m ,   i j 1 2 m + 1 R i j = 1 m ,   i j     0 . 5   R i j = 1 ,   i j   0   R i j = 1 ,   i j
where m is a positive integer not less than 2.
The relative attribute weight of each indicator was constructed, which can be expressed as follows:
ω A H M = 2 n ( n 1 ) j = 1 n k i j
where i = 1, 2, …, n, n is the number of evaluation indicators.

2.4. Coefficient of Variation Method (CVM)

As an objective evaluation method, the CVM is usually based on the standard deviation of each index’s data and the average quotient to obtain a method to measure the degree of variation among the indicators. When the value of a certain index is more different in terms of the overall index data, it means that the index plays a greater role in the evaluation of the overall goal and can be given a greater weight. In contrast, the weight is smaller. Equations (6) and (7) were used to carry out dimensionless processing on the positive correlation index and negative correlation index of the original matrix, respectively, and then the mean Equation (8), standard deviation Equation (9) and coefficient of variation Equation (10) were calculated successively to obtain the weight of each index.
The dimensionless data were processed via the following steps.
For the positive correlation index, kij can be calculated as follows:
k i j = x i j min j ( x i j ) max j ( x i j ) min j ( x i j )
For the negative correlation index, kij can be calculated as follows:
k i j = max j ( x i j ) x i j max j ( x i j ) min j ( x i j )
The average value of each index can be expressed as follows:
χ j ¯ = 1 n i = 1 m χ i j
The standard deviation of each index can be expressed as follows:
D = 1 n i = 1 m ( χ i j χ ¯ ) 2
The coefficient of variation of each index can be expressed as follows:
δ j = D χ j ¯
The weight of each index can be expressed as follows:
ω j = δ j j = 1 n δ j
Among them, j = 1 n ω j = 1 .

2.5. Catastrophe Progression Method (CPM)

The CPM method is an important application of catastrophe theory to evaluation systems. It is an evaluation method with fuzzy function characteristics derived from catastrophe theory and calculus [26]. Zhang et al. [27] predicted the risk of coal and gas explosions in coal mines using the catastrophe progression method. Compared with the commonly used multi-index comprehensive evaluation methods at the present stage, the calculation complexity of the CPM is relatively small, and the operation is simpler. It considers the relative importance of the indicators. It uses the normalized formula to calculate the indicators, solving the problems that the indicators cannot be quantified and the complex calculation procedures. Therefore, it is generally applicable to the study of evaluation index systems.

2.5.1. Determining the Catastrophe Type of Each System of Indicators

In the CPM, according to the number of variables controlled by basic factors, there are four commonly used catastrophe types, including folding catastrophes, pointy catastrophes, dovetail catastrophes and butterfly catastrophes, as shown in Table 5 below.
According to the catastrophe theory, through simultaneous equations f(x)′ and f(x)″, the normalized formula for each catastrophe type was derived, as shown in Table 6 below.

2.5.2. Nondimensionalization of the Raw Data

In the application of the CPM, the value range of control variables should be 0–1. Therefore, dimensionless processing should be carried out on the original data to nondimensionalize the index.

2.5.3. Evaluation with the Normalized Formula

The catastrophe progression of each level should follow the complementary and noncomplementary principles, that is, when there is a correlation between the subordinate indicators, the complementary principle is adopted, and the catastrophe progression of the upper level is the normalized average of the subordinate indicators; otherwise, the noncomplementary principle is adopted, and the catastrophe progression of the upper level is the minimum value of the normalized value of the subordinate indicators. The Pearson correlation coefficient method was used to determine the correlation between each index, and the following Equation (12) was adopted:
r = i = 1 m ( x i x ¯ ) ( y i y ¯ ) i = 1 m ( x i x ¯ ) 2 i = 1 m ( y i y ¯ ) 2
where m represents the number of samples; xi′ and yi′ represent standardized index values; x ¯ and y ¯ represent the mean value of indicators after standardization; and r is the correlation coefficient. When r is greater than or equal to 0.8, it indicates that the two indicators are highly correlated; 0.8 > r ≥ 0.5 indicates a moderate correlation; when 0.5 > r ≥ 0.3, the correlation is low; and r < 0.3 indicates little correlation.

3. Results

3.1. Weight Determination via the AHM

The discriminant matrix R is obtained according to the expert scoring method, and the results are shown in Table 7 and Table 8.
According to Equation (4) above, the attribute discrimination matrix K was calculated using scale Rij, and the results are shown in Table 9 and Table 10.
According to Equation (5), the relative attribute weight ωAHM of the factors influencing the risk of water inrush in the upper layer was calculated, and the results are shown in Table 11.

3.2. Coefficient of Variation Method to Determine the Weight

According to Equation (6) through to Equation (11), the original data were calculated, and the weight values of each index of the risk of water inrush in the upper layer were obtained, as shown in Table 12.

3.3. Determine the Weight Ranking of Comprehensive Indicators

The following Equation (13) was used to combine the subjective weight value calculated via AHM and the objective weight value calculated via the CVM to obtain the comprehensive weight ωi of each influencing factor:
ω i = α ω A H M + ( 1 α ) ω j
where α is the preference coefficient, which represents the degree of preference of the analyzer and judge of the influencing factors in terms of the proportion of subjective weight in the comprehensive weight. In this paper, the value is 0.7 after considering the practical significance of the evaluation system of the risk of water inrush in the upper layer. The comprehensive weight of each factor was obtained, as shown in Table 13 below.
The ranking results of comprehensive indexes were as follows: the thickness of the aquiclude C3, the hydrostatic head of the Zhidan Formation C5, the lithologic assemblage index of the aquiclude C2, the lithologic structure index of the Zhidan Formation C4 and the core recovery rate from the aquiclude C1.

3.4. Establish an Improved Catastrophe Progression Method Evaluation Model

The core recovery rate C1, lithologic assemblage index C2 and thickness of the aquiclude C3 can be used to composite a dovetail catastrophe. The lithologic index of the Zhidan Formation C4 can be used to composite a folding catastrophe. The hydrostatic head of the Zhidan Formation C5 can be used to composite a folding catastrophe. The water-retaining capacity of key layer B1, water-retaining capacity of the Zhidan Formation B2 and water inrush pressure B3 can be used to composite the dovetail catastrophe.
Equations (6) and (7) above were used to process the original data and the standardized data were obtained, as shown in Table 14.
The correlation coefficient matrix of indicators was obtained by calculating the original data in Table 14 according to Formula (11), as shown in Table 15 below.
As seen from the data in Table 15, the relationship among indicators is either a weak correlation or noncorrelation, which meets the principle of non-complementarity. Therefore, the principle of taking the minimum value is adopted. The results of the water inrush risk from an overlying separation layer calculated using borehole data are shown in Table 16.
The natural breakpoint method is a method that classifies the dataset according to the discontinuous places in the dataset. According to the natural breakpoint method, the risk level of water inrush in the overburden is divided into three levels, namely, 0–0.7 is the safe zone, 0.7–0.8 is the transition zone and 0.8–0.9 is the danger zone; thus, a zoning diagram of water inrush risk in the overburden is established. The areas with a high risk of water inrush are mainly concentrated in the southern part of the mining area, while the risk in the northern area is relatively small (Figure 8).

4. Conclusions

(1)
The core recovery rate, lithologic assemblage index, key aquiclude thickness, hydrostatic head and lithologic structure index of the Zhidan Formation are selected as the basis for studying of the water inrush from an overlying separation layer. A risk prediction model based on combination weighting and an improved catastrophe progression method is constructed to predict the risk of water inrush from the overlying separation layer of the double-layer structure in the 221 mining area of Shilawusu coal mine.
(2)
This paper applies the subjective and objective weighting method to improve the evaluation index locally. Then, the typical sample of high water inrush risk levels is established. To a certain extent, it reduces the disadvantage of difficult weight allocation caused by complex water inrush disaster factors, enriches the evaluation and research system of water inrush risk prediction and has great significance for avoiding the threat of water inrush from the separation layer.
(3)
Due to the complexity of the induced mechanism and disaster factors of overlying separation layer water inrush, the quantification and weight determination of indicators still need to be discussed and studied. In addition, the indicators selected in this paper and the risk prediction of overlying separation layer water inrush are all based on the Shilawusu coal mine, and more robust data need to be collected for the whole mining area and the critical parameters of the mathematical model for further discussion.

Author Contributions

All authors contributed to the study conception and design. Writing—review and editing, D.X.; Investigation, C.H.; Writing—original draft, Z.D.; Formal analysis, J.H. and J.Y.; Visualization, J.W. All authors have read and agreed to the published version of the manuscript.

Funding

This research was financially supported by the National Natural Science Foundation of China (Grant No. 41702305).

Conflicts of Interest

On behalf of all authors, the corresponding author states that there are no conflict of interest.

References

  1. Wang, M.-J.; Yu, D.-H. Simulation of water hazards caused by burst of cells formed by overburden stratum separation. Chin. J. Geotech. Eng. 2010, 32, 231–236. (In Chinese) [Google Scholar]
  2. He, J.-H.; Li, W.-P.; Qiao, W. A rock mechanics calculation model for identifying bed separation position and analyzing overburden breakage in mining. Arab. J. Geosci. 2020, 13, 920. [Google Scholar] [CrossRef]
  3. Li, H.; Chen, Q.; Shu, Z.; Li, L.; Zhang, Y. On prevention and mechanism of bed separation water inrush for thick coal seams: A case study in China. Environ. Earth Sci. 2018, 77, 759. [Google Scholar] [CrossRef]
  4. Qiao, W.; Wang, Z.W.; Li, W.P.; Lv, Y.G.; Li, L.G.; Huang, Y.; He, J.H.; Li, X.Q.; Zhao, S.L.; Liu, M.N. Formation mechanism, disaster-causing mechanism and prevention technology of roof bed separation water disaster in coal mines. J. China Coal Soc. 2021, 46, 507–522. (In Chinese) [Google Scholar]
  5. Gui, H.; Lin, M.; Song, X. Features of separation water hazard in China coalmines. Water Pract. Technol. 2017, 12, 146–155. [Google Scholar] [CrossRef]
  6. Yan, H.; He, F.; Yang, T.; Li, L.; Zhang, S.; Zhang, J. The mechanism of bedding separation in roof strata overlying a roadway within a thick coal seam: A case study from the Pingshuo Coalfield, China. Eng. Fail. Anal. 2016, 62, 75–92. [Google Scholar] [CrossRef]
  7. He, J.H.; Li, W.P.; Liu, Y.; Yang, Z.; Liu, S.L.; Li, L.F. An improved method for determining the position of overlying separated strata in mining. J. Eng. Fail. Anal. 2018, 83, 17–29. [Google Scholar] [CrossRef]
  8. Hu, R.; Wu, J.; Zhai, X.; Liu, W. Recognition and Prevention of Bed Separation Water: Based on Trapezoid Platform Model. Geofluids 2021, 2021, 6626268. [Google Scholar] [CrossRef]
  9. Ji, Y.; Cao, H.; Zhao, B. Mechanism and Control of Water Inrush from Separated Roof Layers in the Jurassic Coalfields. Mine Water Environ. 2021, 40, 357–365. [Google Scholar] [CrossRef]
  10. Wu, L.; Bai, H.; Ma, D. Prediction and Prevention of Water Inrush Hazards from Bed Separation Space. Mine Water Environ. 2021, 40, 657–670. [Google Scholar] [CrossRef]
  11. Gui, H.; Tong, S.; Qiu, W.; Lin, M. Research on preventive technologies for bed-separation water hazard in China coal mines. Appl. Water Sci. 2018, 8, 7. [Google Scholar] [CrossRef]
  12. Gui, H.; Lin, M.; Song, X. Identification and Application of Roof Bed Separation (Water) in Coal Mines. Mine Water Environ. 2018, 37, 376–384. [Google Scholar] [CrossRef]
  13. Wang, Z.; Zhang, Q.; Shao, J.; Zhang, W.; Wu, X.; Zhu, X. New type of similar material for simulating the processes of water inrush from roof bed separation. J. ACS Omega 2020, 5, 30405–30415. [Google Scholar] [CrossRef]
  14. Wang, Z.; Zhang, Q.; Zhang, W. A novel collaborative study of abnormal roof water inrush in coal seam mining based on strata separation and wing crack initiation. Eng. Fail. Anal. 2022, 142, 106762. [Google Scholar] [CrossRef]
  15. Li, L.; Shu, Z.-Y.; Feng, Y.-F. Analysis and prevention on mechanism of water inrush from bed separation water of overburden by fully-mechanized top coal caving mining in ultra thick seam. Coal Sci. Technol. 2018, 46, 175–182. (In Chinese) [Google Scholar] [CrossRef]
  16. Xie, D.; Han, J.; Zhang, H.; Wang, K.; Du, Z.; Miao, T. Risk Assessment of Water Inrush from Coal Seam Roof Based on Combination Weighting-Set Pair Analysis. Sustainability 2022, 14, 11978. [Google Scholar] [CrossRef]
  17. Yuan, S.; Sui, W.; Han, G.; Duan, W. An optimized combination of mine water control, treatment, utilization, and reinjection for environmentally sustainable mining: A case study. J. Mine Water Environ. 2022, 41, 828–839. [Google Scholar] [CrossRef]
  18. Zhu, T.; Li, W.; Chen, W. Risk assessment of Cretaceous water inrush in the Ordos Basin based on the FAHP-EM. Water Policy 2021, 23, 1249–1265. [Google Scholar] [CrossRef]
  19. Li, B.; Wu, Q.; Meng-yu, C. Risk analysis model of water inrush through the seam floor based on set pair analysis. J. Mine Water Environ. 2018, 37, 281–287. [Google Scholar]
  20. Yuan, S.; Sun, B.; Han, G.; Duan, W.; Wang, Z. Application and Prospect of Curtain Grouting Technology in Mine Water Safety Management in China: A Review. Water 2022, 14, 4093. [Google Scholar] [CrossRef]
  21. Fan, G.; Zhang, S.; Zhang, D.; Zhang, C.; Chen, M.; Li, Q. An Index of Aquiclude Destabilization for Mining-Induced Roof Water Inrush Forecasting: A Case Study. Water 2019, 11, 2170. [Google Scholar] [CrossRef]
  22. Li, L.; Xie, D.; Wei, J.; Yin, H.; Li, G.; Man, X.; Zhang, W. Analysis and control of water inrush under high-pressure and complex karstic water-filling conditions. Environ. Earth Sci. 2020, 79, 493. [Google Scholar] [CrossRef]
  23. Zhang, P.; Zhu, H.; Wu, Y.; Duan, Z.; Niu, H.; Li, F. State-of-the-art of mechanism of water inrush from bed separation and key technology of prevention and pre-control in China. J. Eng. Geol. 2021, 29, 1057–1070. (In Chinese) [Google Scholar]
  24. Li, Z.; Wang, Y.; Liu, S.-Q. Analysis of influence factors on rock quality designation and core recovery percentage. J. Min. Sci. Technol. 2017, 2, 511–518. (In Chinese) [Google Scholar]
  25. Hou, E.K.; Tong, R.J.; Wang, S.J.; Feng, J.; Chen, T. Prediction method for the water enrichment of weathered bedrock based on Fisher model in Northern Shaaxi Jurassic coalfield. J. China Coal Soc. 2016, 41, 2312–2318. (In Chinese) [Google Scholar]
  26. Tian, L.-G.; Hu, Q.-Y. Study on the Influencing Factors of Ecological-Health in Small Watershed by Catastrophe Progression Method. Yellow River 2020, 42, 71–75+80. (In Chinese) [Google Scholar]
  27. Zhang, T.J.; Ren, S.X.; Li, S.G.; Zhang, T.C.; Xu, H.J. Application of the catastrophe progression method in predicting coal and gas outburst. Min. Sci. Technol. 2009, 19, 430–434. [Google Scholar] [CrossRef]
Figure 1. Location map of the study area.
Figure 1. Location map of the study area.
Sustainability 15 13988 g001
Figure 2. Schematic diagram of double–layer structure.
Figure 2. Schematic diagram of double–layer structure.
Sustainability 15 13988 g002
Figure 3. Schematic diagram of secondary movement of double-layer structure overlying strata.
Figure 3. Schematic diagram of secondary movement of double-layer structure overlying strata.
Sustainability 15 13988 g003
Figure 4. The formation process of the separation layer water hazard: (a) Formation of an overlying separation layer water; (b) Outburst of the separated water body.
Figure 4. The formation process of the separation layer water hazard: (a) Formation of an overlying separation layer water; (b) Outburst of the separated water body.
Sustainability 15 13988 g004
Figure 5. Thematic maps of each main controlling factor. (a) Core recovery rate in key aquiclude. (b) Lithologic assemblage index of key aquiclude. (c) Key aquiclude thickness. (d) Lithologic structure index of the Zhidan Formation. (e) Hydrostatic head.
Figure 5. Thematic maps of each main controlling factor. (a) Core recovery rate in key aquiclude. (b) Lithologic assemblage index of key aquiclude. (c) Key aquiclude thickness. (d) Lithologic structure index of the Zhidan Formation. (e) Hydrostatic head.
Sustainability 15 13988 g005
Figure 6. A hierarchical model for the risk of water inrush from an overlying separation layer.
Figure 6. A hierarchical model for the risk of water inrush from an overlying separation layer.
Sustainability 15 13988 g006
Figure 7. Flow chart of the prediction of water inrush risk.
Figure 7. Flow chart of the prediction of water inrush risk.
Sustainability 15 13988 g007
Figure 8. Prediction map of the overlying water inrush hazard zones.
Figure 8. Prediction map of the overlying water inrush hazard zones.
Sustainability 15 13988 g008
Table 1. Summary table of regional stratigraphic characteristics.
Table 1. Summary table of regional stratigraphic characteristics.
EraSystemSeriesFormationThickness (m)Lithologic Description
CenozoicQuaternaryHoloceneAlluvial sand
(Q4al+pl)
<200Light yellow–brown, yellow medium fine sand and silty sand. Lacustrine sedimentary layer, alluvial–diluvium layer and eolian layer.
Aeolian layer (Q4eol)Brown gravel, grayish yellow sand and silty sand. The thickness of sand layer in the western desert is 0~180 m.
Upper PleistoceneMalan Formation
(Q3m)
0~40Yellowish sandy loess with calcareous nodules, columnar joints, unconformable above all older strata.
MesozoicCretaceousLowerZhidan Formation
(K1zh)
40~230The upper part is light gray, gray–purple, gray–yellow, yellow, purple and red mudstone, siltstone, fine sandstone, sand conglomerate mudstone sand mudstone interbedded with thin layers of argillaceous limestone. Cross-bedding is more developed. Large cross-bedding and channel migration are common at the bottom. Unconformable contact with the underlying floor.
30~280The lower part is light gray, gray–green, brown, red, gray, and purple mudstone, siltstone, sandy mudstone and various grained sandstones and conglomerates, with thin layers of calcareous fine sandstone. Argillaceous cement, relatively loose, with oblique bedding development. Large cross-bedding is common at the bottom. Unconformable contact with the underlying floor.
JurassicMiddleAnding Formation
(J2a)
10~151Light gray, gray–green, yellow, purple and brown mudstone, sandy mudstone, medium sandstone. Calcareous nodules and argillaceous masses with parallel bedding and cross bedding.
Zhiluo Formation
(J2z)
10~400Gray, gray–yellow, gray–green and purplish red mudstone, sandy mudstone, fine sandstone, medium sandstone and coarse sandstone, with cross-bedding and wavy bedding. The lower part sandwiches a thin coal seam and oil shale, containing the 1 coal group. Coal bearing layers 1~3. Pseudoconformity contact with the underlying floor.
LowerYan’an Formation (J1–2y)78~458The sandstones of various gray–gray to grayish white grade are interbedded with dark gray and grayish black sandy mudstone and mudstone. Including coal seams with industrial mining value. Contains the 2, 3, 4, 5, 6, 7 coal groups, 27 layers of coal; among them, the main recoverable coal seams are 2-2 middle, 3-1, 4-1, 4-2 middle, 5-1, 6-1 middle, 6-2. Conformable contact with the underlying floor
TriassicUpperYanchang Formation (T3y)35~312Gray–green, yellow, purple, gray–black massive coarse and medium sandstone, local containing fine gravel, mudstone, siltstone and coal line. Pseudoconformity contact with the underlying floor.
Table 2. Lithologic statistics.
Table 2. Lithologic statistics.
StratumJ1–2yJ2zJ2aK1zh
LithologyAverage Thickness (m)Proportion (%)Average Thickness (m)Proportion (%)Average Thickness (m)Proportion (%)Average Thickness (m)Proportion (%)
Coarse Sandstone3.5879.071111.98726.868
Medium Sandstone10.652114.671719.511286.5125
Fine Sandstone13.312615.91931.4519142.9242
Siltstone6.591310.671222.751450.5915
Mudstone16.923335.064179.734834.2310
Table 4. Statistical table of risk assessment indexes.
Table 4. Statistical table of risk assessment indexes.
BoreholeC1 (%)C2C3 (m)C4 (m)C5 (m)
S0158.190.43191.9318.42288.85
S0373.500.34155.5143.91328.60
S0462.500.36166.4643.14306.08
S0563.200.30128.2671.37344.27
S0767.570.46142.90208.82322.35
S0967.670.54166.33192.63304.54
S1078.440.48169.33226.25302.72
S1374.000.45216.74148.59250.50
S1569.200.4563.54240.85311.26
S1670.800.3674.53104.20300.87
S1863.200.3257.4563.94299.67
N1876.430.45103.8034.03366.79
N4267.000.42148.3518.66345.27
N1975.500.39133.6426.23347.56
N4678.710.45137.05121.11335.24
N4779.500.50159.67212.97304.54
N5164.750.5879.5229.18286.04
N5689.000.4052.15133.88335.96
N5877.170.42139.53233.49324.82
k374.120.46195.89230.67306.67
k784.750.56111.4028.12346.17
k874.710.48176.5024.47322.89
k1289.270.4398.6262.42366.32
k1378.800.47183.7850.69309.37
k1472.640.48201.5155.87263.78
k1578.190.53179.41215.75280.90
k1677.310.46168.23103.43301.40
k2188.130.46190.28125.00280.69
k2264.130.47182.1755.87294.83
k2388.430.64206.75255.27283.76
k2868.000.54161.8489.74306.97
k2986.930.42196.31108.47281.51
k3080.500.45185.9347.37297.37
k3186.850.56186.37117.57289.28
k3982.500.38166.2299.85296.05
k4072.640.48149.14109.50326.47
k4667.000.6637.90109.36310.54
k4761.670.3935.5950.62297.55
k4874.500.51196.9038.70269.62
k5458.670.6955.5317.04290.48
k5555.200.5761.8214.80286.78
k6167.750.5773.5664.70286.87
k6259.430.58106.5043.41278.41
k6363.570.46108.1619.93280.27
k7172.570.55106.6958.03266.25
k7281.000.53118.9167.70256.61
k7570.250.4192.9021.03287.87
Table 5. The four catastrophe models.
Table 5. The four catastrophe models.
Type of CatastropheControl VariableState VariablePotential Function
Folding11 f ( x ) = x 3 + μ x
Pointy21 f ( x ) = x 4 + μ x 2 + ν x
Dovetail31 f ( x ) = x 5 + μ x 3 + ν x 2 + ω x
Butterfly41 f ( x ) = x 6 + μ x 4 + ν x 3 + ω x 2 + t x
Table 6. Normalized value of catastrophe progression method.
Table 6. Normalized value of catastrophe progression method.
Type of CatastropheControl VariableState VariableNormalized Formula
Folding11 x μ = μ
Pointy21 x u = u , x v = v 3
Dovetail31 x u = u , x v = v 3 , x w = w 4
Butterfly41 x u = u , x v = v 3 , x w = w 4 , x t = t 5
Table 7. Discriminant matrix R (A–Bi, i = 1, 2, 3).
Table 7. Discriminant matrix R (A–Bi, i = 1, 2, 3).
Evaluation IndexB1B2B3
B1153
B21/511/3
B31/331
Table 8. Discriminant matrix R (B1–Ci, i = 1, 2, 3).
Table 8. Discriminant matrix R (B1–Ci, i = 1, 2, 3).
Evaluation IndexC1C2C3
C1134
C21/313
C31/41/31
Table 9. Attribute discriminant matrix K (A–Bi, i = 1, 2, 3).
Table 9. Attribute discriminant matrix K (A–Bi, i = 1, 2, 3).
Evaluation IndexB1B2B3
B1010/116/7
B21/1101/7
B31/76/70
Table 10. Attribute discriminant matrix K (B1–Ci, i = 1, 2, 3).
Table 10. Attribute discriminant matrix K (B1–Ci, i = 1, 2, 3).
Evaluation IndexC1C2C3
C106/78/9
C21/706/7
C31/91/70
Table 11. AHM analytic hierarchy process to calculate index weights.
Table 11. AHM analytic hierarchy process to calculate index weights.
Influencing FactorsC1C2C3C4C5
Attribute weight ωAHM0.04980.19620.34270.080.33
Table 12. Coefficient of variation method to calculate index weight.
Table 12. Coefficient of variation method to calculate index weight.
Influencing FactorsC1C2C3C4C5
Index weight ωj0.18000.15590.20670.29300.1644
Table 13. Comprehensive weight value of influencing factors.
Table 13. Comprehensive weight value of influencing factors.
Influencing FactorsC1C2C3C4C5
Index weight ωi0.088860.184110.30190.14390.28032
Table 14. Standardization of indicator data.
Table 14. Standardization of indicator data.
BoreholeC1C2C3C4C5
S010.91230.33160.13690.01510.3298
S030.46290.11860.33800.12110.6716
S040.78580.15410.27750.11790.4779
S050.76520.00000.48840.23530.8064
S070.63690.40990.40760.80690.6179
S090.63410.60680.27830.73950.4647
S100.31800.46610.26170.87930.4490
S130.44820.37860.00000.55640.0000
S150.58910.38330.84570.94010.5225
S160.54220.15750.78500.37180.4331
S180.76520.05860.87930.20440.4228
N180.37700.39010.62340.08001.0000
N420.65370.30450.37750.01610.8150
N190.40420.24600.45870.04760.8347
N460.30990.38220.43990.44210.7287
N470.28680.52130.31500.82410.4647
N510.71970.71870.75750.05980.3056
N560.00800.26210.90850.49520.7349
N580.35530.30240.42620.90940.6391
k30.44480.41370.11510.89770.4830
k70.13270.68030.58150.05540.8227
k80.42730.45540.22210.04020.6225
k120.00000.33660.65200.19810.9960
k130.30740.43210.18200.14930.5062
k140.48810.46030.08410.17080.1142
k150.32530.58330.20610.83570.2614
k160.35120.41790.26780.36860.4377
k210.03370.41930.14610.45830.2596
k220.73810.42830.19080.17080.3812
k230.02480.87140.05511.00000.2860
k280.62430.62600.30310.31170.4856
k290.06880.32120.11280.38950.2667
k300.25750.39500.17010.13540.4031
k310.07110.67030.16770.42740.3335
k390.19880.20550.27890.35370.3917
k400.48810.45990.37310.39380.6533
k460.65370.91690.98720.39320.5163
k470.81020.24341.00000.14900.4046
k480.43360.54380.10950.09940.1644
k540.89831.00000.88990.00930.3438
k551.00000.68930.85520.00000.3120
k610.63170.70340.79040.20750.3128
k620.87590.72570.60860.11900.2400
k630.75430.40330.59940.02130.2560
k710.49020.65590.60750.17980.1354
k720.24280.59930.54010.22000.0525
k750.55830.29240.68360.02590.3214
Table 15. Correlation coefficients between indicators.
Table 15. Correlation coefficients between indicators.
C1C2C3C4C5
C11.0000
C20.03081.0000
C30.37430.08381.0000
C4−0.36540.0705−0.25591.0000
C5−0.1528−0.34340.2174−0.03691.0000
Table 16. Numerical calculation results of water inrush hazard risk from an overlying separation layer.
Table 16. Numerical calculation results of water inrush hazard risk from an overlying separation layer.
BoreholeABoreholeABoreholeABoreholeA
S010.5919N420.5967k140.5384k460.8898
S030.7009N190.6833k150.6737k470.8081
S040.7258N460.6512k160.7193k480.5752
S050N470.7491k210.6182k540.5575
S070.799N510.7032k220.6609k550.5523
S090.7263N560.5469k230.4845k610.8215
S100.7152N580.6733k280.7419k620.7663
S130k30.5824k290.5795k630.6182
S150.8522k70.6965k300.6421k710.7166
S160.7348k80.6692k310.6398k720.6119
S180.6231k120k390.7267k750.6334
N180.7292k130.6531k400.7815
Disclaimer/Publisher’s Note: The statements, opinions and data contained in all publications are solely those of the individual author(s) and contributor(s) and not of MDPI and/or the editor(s). MDPI and/or the editor(s) disclaim responsibility for any injury to people or property resulting from any ideas, methods, instructions or products referred to in the content.

Share and Cite

MDPI and ACS Style

Xie, D.; Du, Z.; Han, C.; Han, J.; Wei, J.; Yan, J. Prediction of the Water Inrush Risk from an Overlying Separation Layer in the Thick Overburden of a Thick Coal Seam. Sustainability 2023, 15, 13988. https://doi.org/10.3390/su151813988

AMA Style

Xie D, Du Z, Han C, Han J, Wei J, Yan J. Prediction of the Water Inrush Risk from an Overlying Separation Layer in the Thick Overburden of a Thick Coal Seam. Sustainability. 2023; 15(18):13988. https://doi.org/10.3390/su151813988

Chicago/Turabian Style

Xie, Daolei, Zhongwen Du, Chenghao Han, Jie Han, Jiuchuan Wei, and Jiulei Yan. 2023. "Prediction of the Water Inrush Risk from an Overlying Separation Layer in the Thick Overburden of a Thick Coal Seam" Sustainability 15, no. 18: 13988. https://doi.org/10.3390/su151813988

APA Style

Xie, D., Du, Z., Han, C., Han, J., Wei, J., & Yan, J. (2023). Prediction of the Water Inrush Risk from an Overlying Separation Layer in the Thick Overburden of a Thick Coal Seam. Sustainability, 15(18), 13988. https://doi.org/10.3390/su151813988

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop