Abstract
Deep mining in Ordos Basin faces the threat of bed-separation water inrush (BWI), so it is necessary to carry out BWI risk assessment before mining. However, the existing BWI risk assessment methods fail to provide a universal risk classification standard. In this paper, taking three coal mines in Ordos Basin as research cases, the water resistance of fractured rock mass was analyzed, and a BWI coefficient (BWIC) method for BWI risk assessment was established. Firstly, by comparing the weight of each BWI-related factor, the formula for calculating the thickness of equivalent water-resisting layer of fractured rock mass was derived. The weight of each BWI-related factor was obtained by entropy weight method. Secondly, based on the stress arch theory, the overall BWIC in bed separation zone corresponding to different excavation lengths was obtained. BWIC was expressed by the ratio of the Cretaceous water pressure to the total thickness of the unbroken water-resisting layer and the equivalent water-resisting layer of the fractured rock mass. Finally, by comparing the value range of the overall BWIC in the bed separation zone corresponding to the excavation length with and without BWI, the standard of BWI risk classification was determined. The application results verified the scientificity of BWIC method and the universality of the BWI risk classification standard proposed in this paper.
1. Introduction
Bed separation is a layered cavity formed between adjacent strata due to uneven settlement of strata in overburden during underground coal mining [1–5]. When the water in the surrounding aquifer accumulates into the bed separation, bed-separation water is formed. Once the bed-separation water breaks through the water-resisting layer below it and flows into the mining site, it will cause bed-separation water inrush (BWI) [6]. BWI has the characteristics of large instantaneous water volume and strong destructiveness. Coal is still the most important fossil energy in China, and Ordos Basin has become China’s main coal producing area [7, 8]. In underground coal mining, because China has the most complex geological conditions and the largest coal mining scale in the world, almost all accidents caused by BWI occur in China. For example, in 2016, a BWI occurred at Zhaojin coal mine in China, resulting in 11 deaths. Therefore, for those coal mines facing BWI threat in Ordos Basin, BWI risk assessment has important reference value for the formulation of safe mining scheme.
At present, mine engineers still have different views on the formation mechanism of BWI. Some views are that the power to induce BWI comes from the impact of the fracture of the rock stratum constituting the bed separation on the bed-separation water [9–11], while some views are that the bed-separation water pressure is the power to drive the bed-separation water to break through the aquiclude [12, 13]. In fact, in the overburden, the inrush of bed-separation water in the fracture zone should be caused by the fracture of rock stratum, but the inrush of bed-separation water above the fracture zone should be driven by bed-separation water pressure. Although BWI is not a new type of roof water hazards in coal mining in recent years, but due to the lack of attention, there are few research results on risk assessment of BWI. Multifactor-weighted superposition method is a mature risk assessment method, and it is often used in mine roof/floor water inrush risk assessment [14, 15]. However, in the multifactor-weighted superposition method, the result of data standardization processing will be affected by the value range of the data set, which will lead to the BWI risk may be wrongly evaluated and the boundary values between different risk levels changing with the change of evaluation area, for example, assuming that the values of one of the BWI-related factors in mining areas A and B are data set and data set , respectively. The maximum and minimum values in data set A are and , respectively, and the maximum and minimum values in data set B are and , respectively. The value of the BWI-related factor at coordinate point X in mining area A is , and the value of the BWI-related factor at coordinate point Y in mining area B is , and . If and cannot exist at the same time, the standardized values of and will not be equal, and the BWI risk at coordinate point X and coordinate point Y will be evaluated as unequal, which is not in line with the reality. The water inrush coefficient method, which was improved from the empirical formula suitable for floor water inrush risk assessment and suitable for roof water inrush risk assessment, has been used in BWI risk assessment [16]. However, the evaluation results were not convincing because the water-resisting capacity of the fractured Jurassic strata in Ordos Basin due to the fracture self-healing effect was not considered.
The above analysis shows that there is still a lack of practical methods suitable for BWI risk assessment in Ordos Basin. Moreover, the BWI risk corresponding to different excavation lengths cannot be obtained by using the existing BWI risk assessment methods, which is not conducive to guiding relevant disaster prevention and control. Therefore, the purpose of this paper is to analyze the role of BWI-related factors in BWI and establish a BWI risk assessment method that can provide water inrush risk corresponding to different excavation lengths and universal risk classification standard.
2. Engineering Geological Background of Coal Mining in Ordos Basin
Ordos Basin is the main coal producing area in China, and Jurassic coal seams are the main excavation objects. In the middle of the basin, a thick layer of Cretaceous strata is deposited above the Jurassic Anding Formation (Figure 1). The Cretaceous strata are mainly composed of thick sandstone with a small amount of sandy mudstone intercalation, which has good integrity (Figures 2 and 3(a)) and is not easy to deform, while the Jurassic Anding Formation is mainly composed of mudstone and sandy mudstone and has weathering zone inside, which is easy to deform (Figures 3(b) and 4). In the subsidence caused by underground mining, due to the asynchronous deformation between sandstone and sandy mudstone in Cretaceous and between Cretaceous and Jurassic Anding Formation, layered cavities (i.e., bed separations) appear in Cretaceous and the interface between Cretaceous and Jurassic Anding Formation (Figure 5). When the water in the aquifer around a bed separation converges into the bed separation, bed-separation water is formed, which can lead to the decline of the water level in the surrounding aquifer (Figure 6). When the bed-separation water suddenly burst into the mining area, BWI is formed. Due to the deformation of sandy mudstone intercalation in Cretaceous is clamped by its adjacent thick sandstone layers, bed separations in Cretaceous have smaller scale. On the contrary, because the overall antideformation ability of Cretaceous is greater than that of Jurassic Anding Formation, the deformation difference between Cretaceous and Jurassic Anding Formation is larger; bed separations at the interface of Cretaceous and Jurassic Anding Formation have larger scale. Therefore, the bed-separation water at the interface between Cretaceous and Jurassic Anding Formation is the focus of this study.
(a)
(b)
3. Methods and Materials
3.1. BWI-Related Factors and Their Quantification
As shown in Figure 5, BWI is caused by the breakthrough of the Jurassic protective layer by the Cretaceous bed-separation water above it. Therefore, the research objects directly involved in BWI process include the Cretaceous bed-separation water and the Jurassic protective layer. In the Jurassic protective layer, the upper part is the unbroken Jurassic strata (UJS), and the lower part is the fractured Jurassic strata (FJS). The water-resisting capacity of FJS is due to the closure of fractures in FJS under the action of water flow. Furthermore, the strata in FJS can be divided into fractured mudstone and fractured sandstone. The role and quantification of each BWI-related factor in BWI are described below. (1)Bed-separation water pressure ()
The bed-separation water pressure () provides the inducing power for the occurrence of BWI. Since the bed-separation water comes from the aquifer around the bed separation, when the bed separation is filled with water, will gradually increase to be consistent with the water pressure in the aquifer around the bed separation. Since the bed separations at the interface between Cretaceous and Jurassic strata are the focus of this study, in this study, the Cretaceous can be expressed by the water pressure in the aquifer at the bottom of Cretaceous. (2)UJS thickness ()
UJS is the key factor affecting BWI because it can prevent the release of bed-separation water. Because the integrity of UJS is not damaged, it has greater strength and better water-resisting capacity than FJS. The greater the thickness of UJS () is, the stronger its water-resisting capacity is. can be calculated as follows: where is the thickness of Jurassic strata between coal seam and Cretaceous (that is, the thickness of Jurassic protective layer); is the thickness of FJS; is the mining thickness; and is the ratio of and , which can be obtained through field monitoring. (3)Closure potential of mudstone fractures in FJS ()
Due to the low strength (with an average saturated uniaxial compressive strength of 4.75 MPa) and low softening coefficient (with an average softening coefficient of 0.34), the Jurassic mudstone can be easily argillated and deformed after encountering water (Figure 7), resulting in the self-healing effect of fractures in Jurassic mudstone under the action of water flow.
In FJS, the better the closure effect of mudstone fractures, the stronger the water-resisting capacity of mudstone, and the smaller the probability of BWI. Here, the closure potential of mudstone fractures in FJS () was used to evaluate the closure effect of mudstone fractures. Because is positively correlated with the volume of mudstone and the volume of mudstone is positively correlated with the mudstone thickness, is positively correlated with the mudstone thickness. On the contrary, because is negatively correlated with the opening degree of mudstone fractures, and the opening degree of mudstone fractures is positively correlated with the height of available subsidence space of mudstone before mudstone fracture, is negatively correlated with the height of available subsidence space of mudstone before mudstone fracture. Therefore, can be quantified by the following formula: where , , and refer to the closure potential, thickness, and height of available settlement space of the th mudstone layer in FJS, respectively; refers to the average expansion coefficient of the fractured rock mass, which can be obtained through field monitoring; and refers to the total thickness of all strata between the th mudstone layer in FJS and the coal seam. (4)Closure potential of sandstone fractures in FJS ()
Due to the low cementation degree, the Jurassic sandstone is easy to disintegrate after encountering water (Figure 8), and the fractures in the Jurassic sandstone are easy to be filled and closed by the disintegrated rock debris.
In FJS, the better the closure effect of sandstone fractures, the stronger the water-resisting capacity of sandstone, and the smaller the probability of BWI. Here, the closure potential of sandstone fractures in FJS () was used to evaluate the closure effect of sandstone fractures. Because is positively correlated with the sandstone fracture length and the sandstone fracture length is positively correlated with the sandstone thickness, is positively correlated with the sandstone thickness. On the contrary, because is negatively correlated with the opening degree of sandstone fractures, and the opening degree of sandstone fractures is positively correlated with the height of available subsidence space of sandstone before sandstone fracture, is negatively correlated with the height of available subsidence space of sandstone before sandstone fracture. Therefore, can be quantified by the following formula: where , , and refer to the closure potential, thickness, and height of available settlement space of the th sandstone layer in FJS, respectively; and refers to the total thickness of all strata between the th sandstone layer in FJS and the coal seam.
3.2. Bed-Separation Water Inrush Coefficient (BWIC) Method
3.2.1. Obtaining the BWI-Related Factors at Each Coordinate Point
According to the strata information revealed by each borehole, the values of BWI-related factors in each borehole are calculated. Then, through interpolation calculation, the values of each BWI-related factor at any coordinate point in the area to be evaluated can be obtained. Interpolation calculation can be realized by using a drawing software named Surfer.
3.2.2. Identifying the Weight of Each BWI-Related Factor
UJS with a certain thickness and the same water-resisting capacity as FJS was called the equivalent strata of FJS (E-FJS). In order to calculate the thickness of E-FJS, it is necessary to determine the weight of each BWI-related factor in BWI.
In this paper, the mining areas 106A and 108 in Shilawusu coal mine are taken as examples to explain how to obtain the weight of each BWI-related factor. Mining in mining area 106A has been completed, mining in mining area 108 is still in progress, and 1000 m has been excavated. When the excavation length was 554 m, a BWI occurred in mining area 106A, but no BWI occurred in mining area 108 during mining. Because the total amount of water in a bed separation is limited, the water inflow in a BWI has the characteristics of rapid increase and rapid decrease (Figure 9). However, the instantaneous water inflow of BWI is very large, so BWI is very destructive.
Even though mining areas 106A and 108 have the same mining width, similar mining depth, and similar mining thickness, there are differences in BWI-related factors between the two mining areas, which eventually lead to differences in BWI risk between the two mining areas. Therefore, by comparing the differences of BWI-related factors between mining area 106A and mining area 108, the weight of each BWI-related factor in the formation of BWI can be obtained under the premise of excluding the influence of mining width, mining depth, and mining thickness on BWI. The principle of obtaining the weight of each BWI-related factor is consistent with the basic principle of entropy weight method. Therefore, in this paper, entropy weight method was used to obtain the weight of each BWI-related factor.
The data used for weight calculation obtained from sampling points and BWI-related factors is written as:
After standardization, the following results can be obtained: where is the standardization value of the th sampling point for the th BWI-related factor. For the BWI-related factors positively correlated with BWI, is calculated as follows: whereas for the BWI-related factors negatively correlated with BWI, is calculated as follows:
The information entropy of the th BWI-related factor () is defined as: where is calculated as follows:
Because when infinitely approaches 0, approaches 0; here, when , the value of was set to 0. The entropy weight of the th BWI-related factor () can be calculated by the following formula:
In this paper, there are four BWI-related factors, so , and , , , and refer to the weight of , , , and , respectively.
3.2.3. Calculating the E-FJS Thickness ()
After calculating the weight of each BWI-related factor in BWI, the water-resisting capacity index () of Jurassic protective layer at the th sampling point can be calculated by the following formula:
Equation (14) can be transformed into: where , , and are constants.
Assuming that the values of , , and at the coordinate point A are , , and , respectively, and the values of , , and at the coordinate point B are , , and , respectively, then according to Equation (15), the water-resisting capacity index of Jurassic protective layer at the coordinate points A and B can be shown as follows: where and are the water-resisting capacity index of Jurassic protective layer at the coordinate points A and B, respectively. Then, the restriction conditions (18) and (19) were set as follows:
By synthesizing Equations (16) and (17) and restriction condition (18), the following equation can be obtained: where is a constant and its value is the ratio of to . From Equation (20), it can be concluded that:
The physical meaning of Equation (21) can be expressed as when is dimensionless and the unit of is m, the water-resisting capacity of UJS with is equal to that of one unit of . UJS with a certain thickness and the same water-resisting capacity as mudstone in FJS was called equivalent strata of mudstone in FJS (E-M-FJS). Therefore, the E-M-FJS thickness () can be obtained as:
By synthesizing Equations (16) and (17) and restriction condition (19), the following equation can be obtained: where is a constant and its value is the ratio of to . From equation (23), it can be concluded that:
The physical meaning of Equation (24) can be expressed as when is dimensionless and the unit of is m, the water-resisting capacity of UJS with is equal to that of one unit of . UJS with a certain thickness and the same water-resisting capacity as sandstone in FJS was called equivalent strata of sandstone in FJS (E-S-FJS). Therefore, the E-S-FJS thickness () can be obtained as:
The relationship among , , and satisfies the following equation:
According to Equations (22), (25), and (26), the formula for calculating can be obtained as follows:
3.2.4. Calculating the BWIC at Each Coordinate Point
The formula for calculating the BWIC at each coordinate point can be obtained as follows: where , ,, and , respectively, refer to the BWIC, , , and at the coordinate point . According to Equation (27), Equation (28) can be transformed into: where and , respectively, refer to the and at the coordinate point .
3.2.5. Obtaining the Overall BWIC Corresponding to Different Excavation Lengths
Since any excavation length corresponds to a stress balance arch (Figure 10), the average BWIC in the bed separation zone in the stress balance arch corresponding to each excavation length was used as the overall BWIC corresponding to each excavation length. Therefore, the overall BWI risk when the excavation length is can be expressed by the following formula: where and z, respectively, refer to the overall BWIC and the number of selected coordinate points within the distribution range of bed separation zone when the excavation length is . In this paper, for the convenience of calculation, refers to the number of grid points after each bed separation zone distribution range was divided into squares with side length of 10 m.
With the increase of the excavation length, the size of the stress balance arch expands continuously. When the excavation length reaches a specific value, the size of the stress balance arch develops to the maximum, which is called the limit arch. The maximum height of limit arch is about 0.7 times of mining depth. Then, as the excavation continues, the stress balance arch exists in the form of moving arch, which has the same size as the limit arch [17, 18]. In order to observe the shape of the stress balance arch, some mining cases in Cuimu coal mine, Buerdong coal mine, and Daliuta coal mine in Ordos Basin were studied. The mining process of mining area 21301 in Cuimu coal mine was simulated by using the discrete element software UDEC. The results show that the deformed and broken rock strata are located in a triangular area as a whole, the two sides of the triangle are rock fracture lines, and the bottom of the triangle is the distribution range of goaf. In the triangular area, the deformation and fracture degree of the rock stratum below the bed separation is large, while the deformation of the rock stratum above the bed separation is small (Figure 11(a)). Similar material physical simulation experiments of mining cases in Buerdong coal mine and Daliuta coal mine also show that the rock strata that can deform and fracture are located in the triangular area surrounded by the two rock fracture lines and the distribution area of goaf (Figures 11(b) and 11(c)). Because the deformation and failure of overburden occur in the stress balance arch, the area surrounded by the triangle composed of rock fracture lines and goaf can be equivalent to the area in the stress balance arch, that is, the stress balance arch can be simplified into a triangle with a base angle of about 60° [19, 20].
(a)
(b)
(c)
3.2.6. Grading BWI Risk
In actual mining, the excavation length with or without BWIs is called danger excavation length () and safety excavation length (), respectively. The corresponding to and are and , respectively. greater than indicates a high risk of BWI, while smaller than indicates a low risk of BWI. Therefore, the risk level with greater than was identified as danger, the risk level with smaller than was identified as safety, and the risk level with greater than less than was identified as transition. Further division of the transition can be achieved by means of the average value of and . Here, the average value of and was marked as . The risk level with greater than less than was identified as relatively safe, and risk level with greater than less than was identified as relatively dangerous. In addition, when the excavation length is less than , the BWI risk was determined as safe because the stress balance arch has not expanded to the Cretaceous bed separation zone (Figure 10). The principle of BWI risk classification is shown in Figure 12.
4. Results and Discussion
Taking mining areas 106A and 108 in Shilawusu coal mine as examples, the establishment process of BWIC method was introduced.
4.1. Obtaining the BWI-Related Factors at Each Coordinate Point
The distribution of related factors in mining areas 106A and 108 in Shilawusu coal mine is shown in Figure 13.
(a)
(b)
(c)
(d)
4.2. Identifying the Weight of Each BWI-Related Factor
According to the analysis in Section 3.2, the entropy weight method can be used to determine the weight of BWI-related factors, and the sampling points in mining areas 108 and 106A should be located in the range of the bed separation zone corresponding to the excavation length of 554 m. Based on the triangle stress arch theory, it can be determined that when the excavation length is 554 m, the sampling points in mining area 108 are A1-A9, and the sampling points in mining area 106A are B1-B9 (Figure 13). The values of BWI-related factors at each sampling point are shown in Table 1.
By analyzing the data in Table 1 according to Equations (9)–(13), the weight of each BWI-related factor can be obtained (Table 2).
Because bed-separation water pressure is the power source of BWI, has a large weight. As the integrity of UJS was not damaged, UJS plays an important role in preventing water inrush, so the weight of should be the largest. Because the integrity of the rock strata in FJS has been destroyed, the water-resisting capacity of FJS is still weaker than that of UJS, so the weight of and is less than that of . In addition, because the fracture closure degree of mudstone is relatively better than that of sandstone, the water-resisting capacity of mudstone in FJS is stronger than that of sandstone, so the weight of is greater than that of .
The above analysis results are obtained by reasonable inference based on the actual engineering geological conditions, which are consistent with the results in Table 2, indicating that the calculation results of the weights of BWI-related factors are in line with the reality.
4.3. Calculating the E-FJS Thickness ()
According to the weights of BWI-related factors in Table 2, Equation (15) can be expressed as:
It can be seen from Equation (31) that , , and are , , and , respectively. Then, and can also be calculated, which are 2.039 and 1.301, respectively, and Equations (22) and (25) can be expressed as follows:
Taking mining area 106A in Shilawusu coal mine as an example, when , , and are 10 m, 23, and 1.032, respectively, the of mudstone with a thickness of 1 m at different locations in FJS was calculated by using Equations (4) and (32) (Table 3), and the of sandstone with a thickness of 1 m at different locations in FJS was calculated by using Equations (6) and (33) (Table 4).
Although the rock in FJS still have water-resisting capacity due to the closure of internal fractures, the restored water-resisting capacity is still weaker than that of the unfractured rock because the integrity of the fractured rock has been destroyed. This conclusion can also be obtained by analyzing the data in Table 3 and Table 4. It can be seen from Table 3 and Table 4 that the water resistance capacity of mudstone and sandstone in FJS is weaker than that of intact rock. Due to the larger opening and weaker closure of the fractures, the water-resisting capacity of the rock at the bottom of FJS is less than that of the rock at the top of FJS. In addition, in FJS, due to the self-healing effect of fractures in mudstone, the fractures in mudstone are easier to close than those in sandstone, so that the water-resisting capacity of mudstone at the same location is stronger than that of sandstone.
The above analysis shows that the method of obtaining determined in this paper is practical, and the formula for calculating is as follows:
4.4. Calculating the BWIC at Each Coordinate Point
Since and are 2.039 and 1.301, respectively, Equation (29) can be transformed into:
According to formula (35), the BWIC at each coordinate point in mining areas 108 and 106A can be obtained (Figure 14).
4.5. Obtaining the Overall BWIC Corresponding to Different Excavation Lengths
Based on the triangle stress arch theory, the overall BWIC corresponding to different excavation lengths can be obtained by using Equations (30) and (35) (Figure 15).
4.6. Grading BWI Risk
Because BWI occurred in mining area 106A when , but no BWI occurred in the completed 1000 m mining process in mining area 108, the corresponding to the excavation length of 554 m in mining area 106A was regarded as , and the maximum value of in mining area 108 was used as . According to Figure 15, , , so . Therefore, Figure 12 can be improved to Figure 16.
5. Application and Validation
5.1. Shilawusu Coal Mine
According to formula (35), the BWIC at each coordinate point in mining areas 201 and 201A in Shilawusu coal mine can be obtained (Figure 17). Based on the triangle stress arch theory, the overall BWIC corresponding to different excavation lengths can be obtained by using Equations (30) and (35) (Figure 18). It can be seen from Figure 18 that the BWI risk in mining areas 201 and 201A during mining belongs to the “safety zone.” In fact, the mining in mining areas 201 and 201A has been completed, and no BWI occurred during mining, which verifies the accuracy of the evaluation results shown in Figure 18.
(a)
(b)
5.2. Yingpanhao Coal Mine
According to formula (35), the BWIC at each coordinate point in mining areas 2202 and 2201 in Yingpanhao coal mine can be obtained (Figure 19). Based on the triangle stress arch theory, the overall BWIC corresponding to different excavation lengths can be obtained by using Equations (30) and (35) (Figure 20). It can be seen from Figure 20 that the BWI risk in mining areas 2202 and 2201 during mining belongs to the “relatively safe zone.” In fact, mining has been completed in mining area 2201, and 1000 m excavation has been completed in mining area 2202. No BWI occurred during mining in these two mining areas, which verifies the accuracy of the assessment results shown in Figure 20.
(a)
(b)
5.3. Cuimu Coal Mine
According to formula (35), the BWIC at each coordinate point in mining areas 21301 and 21302 in Cuimu coal mine can be obtained (Figure 21). Based on the triangle stress arch theory, the overall BWIC corresponding to different excavation lengths can be obtained by using Equations (30) and (35) (Figure 22). It can be seen from Figure 22 that the BWI risk in mining areas 21301 and 21302 during mining belongs to the “danger zone.” In fact, the mining in mining areas 21301 and 21302 has been completed, and BWIs occurred during mining (Figure 23), which verifies the accuracy of the evaluation results shown in Figure 22.
(a)
(b)
(a)
(b)
6. Conclusions
The objects directly involved in BWI process include the Cretaceous bed-separation water and the Jurassic protective layer. The bed-separation water pressure () provides the inducing power for the occurrence of BWI and can be expressed by the water pressure in the aquifer at the bottom of Cretaceous. In the Jurassic protective layer, the upper part is the unbroken Jurassic strata (UJS), and the lower part is the fractured Jurassic strata (FJS). The water-resisting capacity of FJS is due to the closure of fractures in FJS under the action of water flow.
The greater the UJS thickness () is, the stronger its water-resisting capacity is. In FJS, the closure potential of mudstone fractures () was used to evaluate the closure effect of mudstone fractures, and the closure potential of sandstone fractures () was used to evaluate the closure effect of sandstone fractures. is positively correlated with the mudstone thickness and negatively correlated with the height of available subsidence space of mudstone before mudstone fracture. is positively correlated with the sandstone thickness, and negatively correlated with the height of available subsidence space of sandstone before sandstone fracture. UJS with a certain thickness and the same water-resisting capacity as mudstone in FJS was called equivalent strata of mudstone in FJS (E-M-FJS). UJS with a certain thickness and the same water-resisting capacity as sandstone in FJS was called equivalent strata of sandstone in FJS (E-S-FJS). By comparing the weight of each BWI-related factor, the formulas for calculating the E-M-FJS thickness () and the E-S-FJS thickness () were derived. And the corresponding to one unit of is 2.039 m, while the corresponding to one unit of is 1.301 m. The weight of each BWI-related factor was obtained by entropy weight method. And the weights of , , , and are 0.354, 0.390, 0.158, and 0.098, respectively.
In bed-separation water inrush (BWI) risk assessment, the existing multifactor-weighted superposition method cannot provide a universal risk classification standard, and the existing water inrush coefficient method ignores the water separation capacity of FJS, so there is still a lack of scientific BWI risk assessment method. In this paper, three coal mines in Ordos Basin were taken as research cases, and a bed-separation water inrush coefficient (BWIC) method for BWI risk assessment was established. The value of BWIC is the ratio of to the sum of , and . Based on the stress arch theory, the overall BWIC in bed separation zone corresponding to different excavation lengths () was obtained. By comparing the value range of with and without BWI, the standard of BWI risk classification was determined. In the BWI risk classification discrimination diagram established based on BWIC method, the area with was divided into safety zone, the area with was divided into relatively safe zone, the area with was divided into relatively dangerous zone, and the area with was divided into danger zone. The application results verified the scientificity of BWIC method and the universality of the BWI risk classification standard proposed in this paper.
Data Availability
The authors confirm that the data supporting the findings of this study are available within the article.
Conflicts of Interest
The authors declare that they have no conflicts of interest.
Acknowledgments
The authors are grateful to the many people who assisted with this study. The study was jointly supported by the Scientific Research Start-Up Fund for High-Level Introduced Talents of Anhui University of Science and Technology (grant no. 2022yjrc21), the Key Projects of the Natural Science Foundation of China (grant no. 41931284), the China Postdoctoral Science Foundation (grant no. 2019M652015), and the Natural Science Foundation of Jiangsu Province (grant no. BK20190646).