Risk Assessment of Water Inrush of a Coal Seam Floor Based on the Combined Empowerment Method

: With the exploitation of the lower coal seams of the Taiyuan Formation, the Ordovician limestone water inrush in the ﬂoor became more serious. This paper considers the 162 and 163 mining areas of the Jiangzhuang Coal Mine, in Shandong Province, China. A comprehensive analysis of the geological and hydrogeological conditions of the mining area revealed the following: water pressure and water richness provide the water source and power for the ﬂoor water inrush; the thickness of the effective aquifer and the ratio of brittle rock can restrain ﬂoor water inrush; fault structures provide water inrush channels; and mining damage is an artiﬁcial interference and increases the probability of water inrush. Therefore, six factors: the water pressure of the Ordovician limestone aquifer, water abundance of the Ordovician limestone aquifer, equivalent thickness of effective aquifuge, brittle rock ratio, fracture structure, and mining destruction were selected as the inﬂuencing factors on water penetration of the bottom plate, and drawing software was used to establish a mining area map of related factors. The improved fuzzy hierarchical analysis method is more suitable for analyzing multi-objective decisions than the traditional hierarchical analysis method, but the weighting of results is inﬂuenced by expert experience. The entropy weight method is data-driven, and the empowerment results are objective. The improved fuzzy analytic hierarchy process and entropy weight method were coupled together, to determine the weight of each factor. The new method is not only data driven, but also takes empirical experience into consideration, making the empowerment results more reasonable. An evaluation of coal ﬂoor water inrush was established using MapGIS10.6, which is a general tool-type geographic information system software developed by the China University of Geosciences; and the risk of Ordovician limestone water inrush in the ﬂoor of the study area was classiﬁed into four levels: dangerous, relatively dangerous, relatively safe, and safe. The whole evaluation process is simple, but the evaluation results have practical importance and are very efﬁcient, providing theoretical support for coal mine water prevention and control engineering.


Introduction
Coal is the most important energy source in China for the medium and long term. With the rapid development of the economy, the demand for coal has also greatly increased, resulting in accelerated coal mining. With the gradual depletion of the resource reserves of shallow coal seams with relatively good mining conditions, the proportion of deep coal seam mining has increased. The old mining area has developed into deep mining. Mining the lower coal seams of the Taiyuan Formation in the coal-bearing area of North it is not overly sophisticated and complicated, so it is easy for people to master and use. The shortcoming of the fuzzy evaluation method is that it applies subjective judgment in the evaluation process. At the same time, it is still necessary to check the consistency of the judgment matrix in the calculation process. Some scholars [26] have improved the fuzzy analytic hierarchy process, changing the reciprocal judgment matrix into the fuzzy consistency judgment matrix, and combining the sum row normalization method or square root method with the eigenvector method. The improved fuzzy analytic hierarchy process not only solves the consistency problem of the judgment matrix, but also improves the convergence speed and accuracy of evaluation, so as to obtain a ranking vector consistent with reality. Nevertheless, because the determination of the weight of each factor is subjective, the improved fuzzy analytic hierarchy process is an evaluation method based on subjective information. The entropy weight method is a decision method that objectively quantifies and synthesizes multi-factor information, assigns values to each factor, and determines the degree of influence of each factor on the target. With its strong objectivity and relatively simple calculation, many scholars [27,28] have adopted this method for evaluations.
In this paper, the two most common methods (analytic hierarchy process and entropy weight method) are used to improve the final coupling, to establish a more functional and accurate prediction model. Finally, the improved fuzzy analytic hierarchy process (IFAHP) and entropy weight method (EW) are coupled to evaluate and predict the risk of Ordovician limestone water inrush, in the lower coal seam of the Taiyuan Formation in the Jiangzhuang Coal Mine of Tengzhou Mining Area, in Shandong Province, China. The new method provides a reliable theoretical basis for the prevention and control of water disasters in mining areas.

General Situation of the Mine
Jiangzhuang Coal Mine is located in the Tengnan mining area of the Tengxian coalfield of southwestern Shandong Province, China. The Jiangzhuang minefield is affected by the Gaomiao fault, Liuxianzhuang fault, and Yinjiawa fault. The overall structural line is distributed in a NNE direction, forming a wide and gentle fold area with a developed fault structure, and main the characteristics of a graben and barrier. There are many faults in the field. Coal bearing strata change along the strike and tendency, and magmatic intrusion occurs in some sections. The Jiangzhuang Mine is in the middle of the Tengnan Coalfield.
The coalfield strata comprise North China type sedimentary rock, with coal seams in the Carboniferous-Permian formations. The strata's ages, from old to new, are Cambrian, Ordovician, Carboniferous, Permian, Jurassic, and Quaternary. The Cambrian strata are underlain by Quaternary sediments, except for sporadic outcrops outside the coalfield.
With the depletion of the upper coal resources in the Jiangzhuang Coal Mine, the mine will gradually start the mining of lower coal seams. The Taiyuan Group #16 coal will be the main mining object in the future. The 162 and 163 mining areas of the #16 coal seam are the initially developed and mined areas of the lower coal group in the mine. The burial depth of the 162 and 163 mining areas increases gradually from west to east, and the threat of floor water inrush increases gradually. In particular, the potentiometric pressure of the underlying Ordovician limestone confined aquifer increases, and the water inrush risk from the Ordovician limestone aquifer increases. Therefore, the study of the water inrush risk of the floor of the #16 coal seam has importance for the future prevention and control of water in the mine. The 162 and 163 mining areas are the research object of this paper, and their location is shown in Figure 1.

Analysis of Influencing Factors on Floor Water Inrush
The proper selection of main control factors is conducive to establishing a hierarchical structure model, which plays a key role in the water inrush risk assessment of a coal seam floor. The water of a coal bed is affected by various factors, such as the geological structure, hydrogeology, and mining conditions [29]. As shown in Figure 2, the existence of a confined aquifer under a coal seam is the material basis of a water outburst, hydraulic pressure and mine pressure are the source of forces, and the water isolation formation is the inhibitory condition. The inhibitory capacity depends on the thickness, strength, and the combination of the water separator [30]. A fault is controlled when the water pressure, mine pressure, and stability of the isolation strata are in relative equilibrium. Through the study of the existing hydrological, geological, and production experience data of the Jiangzhuang Coal Mine and the analysis of water inrush factors, it was determined that the Ordovician limestone aquifer underlying the coal in the mine field is the main aquifer threatening the mining of the lower coal group. Coal seam #16 is seriously threatened by the Ordovician limestone water, because the seam is close to the Ordovician limestone, and the potentiometric pressure in the Ordovician limestone is high. Therefore, the water pressure of the Ordovician limestone aquifer, water richness of the Ordovician limestone aquifer, equivalent thickness of the effective water barrier, brittle rock ratio, fracture structure, and depth of failure of the bottom plate were selected as the main control factors for the water inrush prediction of the coal seam floor.

Analysis of Influencing Factors on Floor Water Inrush
The proper selection of main control factors is conducive to establishing a hierarchical structure model, which plays a key role in the water inrush risk assessment of a coal seam floor. The water of a coal bed is affected by various factors, such as the geological structure, hydrogeology, and mining conditions [29]. As shown in Figure 2, the existence of a confined aquifer under a coal seam is the material basis of a water outburst, hydraulic pressure and mine pressure are the source of forces, and the water isolation formation is the inhibitory condition. The inhibitory capacity depends on the thickness, strength, and the combination of the water separator [30]. A fault is controlled when the water pressure, mine pressure, and stability of the isolation strata are in relative equilibrium. Through the study of the existing hydrological, geological, and production experience data of the Jiangzhuang Coal Mine and the analysis of water inrush factors, it was determined that the Ordovician limestone aquifer underlying the coal in the mine field is the main aquifer threatening the mining of the lower coal group. Coal seam #16 is seriously threatened by the Ordovician limestone water, because the seam is close to the Ordovician limestone, and the potentiometric pressure in the Ordovician limestone is high. Therefore, the water pressure of the Ordovician limestone aquifer, water richness of the Ordovician limestone aquifer, equivalent thickness of the effective water barrier, brittle rock ratio, fracture structure, and depth of failure of the bottom plate were selected as the main control factors for the water inrush prediction of the coal seam floor.

Water Pressure of the Ordovician Limestone Aquifer
The water pressure of the Ordovician limestone aquifer is the source of the power of the floor water inrush. The greater the water pressure, the greater the possibility of water inrush. According to the mine's 2006-1 hydrological drilling, the Ordovician limestone is 57.0 m above the #16 coal floor. The Ordovician limestone water head is −130.0 m below mean sea level, resulting in a water pressure of 4.564 MPa. The Ordovician limestone water pressure of the #16 coal seam floor in the 162 and 163 mining areas is between 2.5 and 4.7 MPa ( Figure 3A). It can be seen from the figure that the water pressure in the northeastern mining area is the largest.

Water Abundance of the Ordovician Limestone Aquifer
The water abundance of the Ordovician limestone aquifer is the material basis for the floor water inrush. Its water abundance and supply conditions determine the water quantity of the floor water inrush and whether the water inrush point is sustained. Its water-richness is affected by its lithology, structure, hydrodynamic circulation conditions, burial depth, and other factors. The water-richness of the Ordovician limestone can be divided into three levels: weak, medium, and rich ( Figure 3B). The water-richness of the Ordovician limestone aquifer in the mining area is weak to medium. It can be seen from the figure that the water abundance in the mining area increases from southwest to northeast.

Equivalent Thickness of the Effective Aquifuge
The aquifuge has an inhibitory effect on the water inrush of the coal seam floor. According to the "lower three zones" (Figure 4) of water inrush from the coal seam floor [31], the effective water-resisting layer plays a important role in water blocking, and the thickness of the effective water-resisting layer is equal to the total thickness of the water-resisting layer, minus the depth of the rock pressure failure zone and the height of the water pressure lifting, which is expressed as Equation (1). The equivalent thickness of the effective waterresisting layer between coal seam #16 and the Ordovician limestone gradually becomes thinner from west to east, between 34.62 m and 37.32 m ( Figure 3C).
the floor water inrush. The greater the water pressure, the greater the possibility of water inrush. According to the mine's 2006-1 hydrological drilling, the Ordovician limestone is 57.0 m above the #16 coal floor. The Ordovician limestone water head is −130.0 m below mean sea level, resulting in a water pressure of 4.564 MPa. The Ordovician limestone water pressure of the #16 coal seam floor in the 162 and 163 mining areas is between 2.5 and 4.7 MPa ( Figure 3A). It can be seen from the figure that the water pressure in the northeastern mining area is the largest.

Water Abundance of the Ordovician Limestone Aquifer
The water abundance of the Ordovician limestone aquifer is the material basis for the floor water inrush. Its water abundance and supply conditions determine the water quantity of the floor water inrush and whether the water inrush point is sustained. Its waterrichness is affected by its lithology, structure, hydrodynamic circulation conditions, burial depth, and other factors. The water-richness of the Ordovician limestone can be divided into three levels: weak, medium, and rich ( Figure 3B). The water-richness of the Ordovician limestone aquifer in the mining area is weak to medium. It can be seen from the figure that the water abundance in the mining area increases from southwest to northeast. According to the empirical formula of the depth of the bottom plate, and considering the mining depth, the inclination angle and slope length of the working face can be calculated according to Equation (2).

Equivalent Thickness of the Effective Aquifuge
In this formula, H 2 -the thickness of the effective water-resisting layer, m; H-the total thickness of the water-resisting layer, m; H 1 -the depth of the rock pressure failure

Brittle Rock Ratio
Brittle rock plays a key role in preventing floor water inrush. Due to the heterogeneity of the rock strata, the ratio of brittle rock in different regions is different, and therefore the water insulation ability is also very different. Figure 3D shows the brittle rock ratio isoline in the 162 and 163 mining areas. The brittle rock ratio gradually increases from southeast to northwest.

Fracture Structure
The fracture structure not only destroys the rock integrity, but also increases water conductivity. The greater the degree of fracture structure development, the higher the risk of water inrush. Therefore, the fracture structure is one of the important factors affecting the floor water inrush. The fractal dimension method [32] is used to synthesize many faultrelated influencing factors, and the fault fractal dimension value is used to represent the influence of fault-related factors, so as to improve the accuracy of evaluation. First, the unit grids are established according to the size ), when the fault scale index is established. Then, each tectonic unit grid is divided into three levels, according to the continuation, and the number of grids containing faults in each cell is counted step by step, when the edge length of each grid is different. Through calculation, four groups of lgr and lgN (r) values (r is the cell edge length and N is the number of lattices with faults) were obtained. The point (lgr, lgN (r)) was plotted in a two-dimensional coordinate system, and a double logarithmic linear regression analysis was carried out. The slope k of the linear function is the fractal value, and a fault fractal dimension contour map of the mining area was drawn ( Figure 3E).

Mining Damage
Mining damage refers to the stope space after coal seam mining, which causes the original stress change of the surrounding rock. When the stress borne by the surrounding rock exceeds its ultimate strength, displacement, cracking, fracture, and breakage will occur. When overburden deformation or displacement changes the water conductivity and water isolation performance of the rock strata, this results in mine damage effects. According to the degree and form of damage, mining damage causes caving and cracking. The damage to a coal seam floor can be expressed by the depth of the floor damage zone. Unlike the roof mining cracks, the floor damage depth of the floor mining cracks makes the effective protective thickness of the water-resisting layer smaller, resulting in a greater risk of water inrush from the floor of coal seam 16, which, in turn, has a great influence on the water filling of coal seam 16. The greater the floor damage depth, the higher the risk of water inrush from coal seam 16. The lower the floor failure depth, the lower the

Brittle Rock Ratio
Brittle rock plays a key role in preventing floor water inrush. Due to the heterogeneity of the rock strata, the ratio of brittle rock in different regions is different, and therefore the water insulation ability is also very different. Figure 3D shows the brittle rock ratio isoline in the 162 and 163 mining areas. The brittle rock ratio gradually increases from southeast to northwest.

Fracture Structure
The fracture structure not only destroys the rock integrity, but also increases water conductivity. The greater the degree of fracture structure development, the higher the risk of water inrush. Therefore, the fracture structure is one of the important factors affecting the floor water inrush. The fractal dimension method [32] is used to synthesize many fault-related influencing factors, and the fault fractal dimension value is used to represent the influence of fault-related factors, so as to improve the accuracy of evaluation. First, the unit grids are established according to the size l 0 × l 0 (l 0 = 200 m), when the fault scale index is established. Then, each tectonic unit grid is divided into three levels, according to the continuation, and the number of grids containing faults in each cell is counted step by step, when the edge length of each grid is different. Through calculation, four groups of lgr and lgN (r) values (r is the cell edge length and N is the number of lattices with faults) were obtained. The point (lgr, lgN (r)) was plotted in a two-dimensional coordinate system, and a double logarithmic linear regression analysis was carried out. The slope k of the linear function is the fractal value, and a fault fractal dimension contour map of the mining area was drawn ( Figure 3E).

Mining Damage
Mining damage refers to the stope space after coal seam mining, which causes the original stress change of the surrounding rock. When the stress borne by the surrounding rock exceeds its ultimate strength, displacement, cracking, fracture, and breakage will occur. When overburden deformation or displacement changes the water conductivity and water isolation performance of the rock strata, this results in mine damage effects. According to the degree and form of damage, mining damage causes caving and cracking. The damage to a coal seam floor can be expressed by the depth of the floor damage zone. Unlike the roof mining cracks, the floor damage depth of the floor mining cracks makes the effective protective thickness of the water-resisting layer smaller, resulting in a greater risk of water inrush from the floor of coal seam 16, which, in turn, has a great influence on the water filling of coal seam 16. The greater the floor damage depth, the higher the risk of water inrush from coal seam 16. The lower the floor failure depth, the lower the risk of Ordovician limestone water inrush from coal seam 16. The northeast of the mining area is seriously damaged ( Figure 3F).

Steps of IFAHP
IFAHP not only conforms to human mental logic, but also has a simple and accurate form. It is more suitable for multi-objective decision-making analyses than traditional AHP [33,34]. The improved fuzzy analytic hierarchy process not only solves the consistency problem of the judgment matrix, but also solves the accuracy problem and convergence speed problem of the solution; so as to obtain the sorting vector with a higher calculation accuracy and greater consistency with the actual situation.
(1) Appropriate selection of the evaluation index system is a key step, which directly affects the evaluation results. Many scholars [35][36][37]  risk of Ordovician limestone water inrush from coal seam 16. The northeast of the mining area is seriously damaged ( Figure 3F).

Steps of IFAHP
IFAHP not only conforms to human mental logic, but also has a simple and accurate form. It is more suitable for multi-objective decision-making analyses than traditional AHP [33,34]. The improved fuzzy analytic hierarchy process not only solves the consistency problem of the judgment matrix, but also solves the accuracy problem and convergence speed problem of the solution; so as to obtain the sorting vector with a higher calculation accuracy and greater consistency with the actual situation.
(1) Appropriate selection of the evaluation index system is a key step, which directly affects the evaluation results. Many scholars [35][36][37]   (2) A complementary fuzzy judgment matrix (priority judgment matrix) was established using the three-scale method with element values of 0, 0.5, and 1. In order to accurately describe the relative importance of any two factors for a certain criterion, this paper used the 0.1~0.9 scale method to determine the element value (Table 1), and establish a priority judgment matrix F = f ij n×n .  (4) Use and line normalization method , obtaining an Order Vector.
(5) Using conversion formula e ij = r ij r ji , transform the complementary judgment matrix R = r ij n×n into a reciprocal matrix E = e ij n×n .
(6) Taking the sorting vector W (0) as the initial value V 0 of the eigenvalue method, the sorting vector W (k) with higher accuracy is further obtained, namely: 1 Taking V 0 = V 0 (v 01 , v 02 , · · ·v 0n ) T as the initial value of iteration, the eigenvector V k+1 is obtained using the iteration formula V k+1 = EV k , and the infinite norm V k+1 ∞ of V k+1 is obtained.
T is normalized after V k+1 , and the obtained vector W (k) = V i+1 is the scheme sorting vector, and the iteration ends. 3 Otherwise, take as the new initial value, reiterate.
(7) Use the method of "soliciting expert ratings" [6][7][8] in Figure 5, to score the factors affecting water inrush, and combine with Table 1 to obtain the corresponding 16 coal floor water inrush index structure system, with three levels of the priority judgment matrix, respectively: The priority judgment matrix of A-B is The priority judgment matrix of B 2 -C is 0.5 0.6 0.4 0.5 Through a MATLAB calculation, the weight of each factor was obtained using IFAHP, see Table 2. The results calculated using IFAHP and AHP are shown in Table 2. The weight ranking of each index obtained by the two methods is different. IFAHP increased the weight of aquifer influencing factors, and the weight of fracture structure and mining failure influencing factors was the largest in AHP. In reality, the water pressure and water yield of an aquifer often play a key role in the factors causing water inrush from a coal seam floor. Compared with the weight distribution method of IFAHP, it is more reasonable and highlights the contribution of the main factors.

Determination of Weight by Entropy Weight Method
The entropy weight method is a decision-making method used to objectively quantify and synthesize multi-factor information, assign values to each factor, and determine the degree of influence of each factor on the target. In information theory, entropy is a measure of uncertainty. The greater the amount of information, the smaller the uncertainty and the smaller the entropy. The smaller the amount of information, the greater the uncertainty and the greater the entropy [42]. The smaller the entropy value of the evaluation index, the greater the information provided, and the greater the role played in the comprehensive judgment of the evaluation results. In recent years, the entropy weight method has been widely used in aquifer evaluation, ecological vulnerability, and economic development investigations [43][44][45].
Taking the standardized matrix R as the research object, with one evaluation object and b evaluation indexes, there is a matrix R = (X ij ) a×b , e i as the entropy value of the ith factor in i influencing factors, which is expressed in Equation (3): In this expression, P ij is the weight of the jth evaluation index value under the ith index, which is expressed as Equation (4): On this basis, all the influencing factors are weighted, and a more objective comparison result can be obtained. The entropy weight of the ith index can be expressed as Equation (5): Through this calculation, the weight of each factor can be obtained using the entropy weight method (Table 3).

Determination of Comprehensive Weight using IFAHP-EW Coupling
The subjective weight W α obtained using the improved fuzzy analytic hierarchy process and the objective weight W β obtained with the entropy weight method are combined to obtain the comprehensive weight W ( Table 4), which is expressed as in Equation (6): In order to eliminate the influence of different dimensions on the superposition calculation, the quantitative data of related factors should be standardized. In this paper, the range standardization Equation (7) is used to standardize the original data.
where Z i is the original data of a factor, and max(x i ) (or min(x i ))is the maximum (or minimum) value in the original data for the same factor.

Establishing the Prediction Model
The linear weighting method is used to evaluate each factor, and the results are obtained using Equation (8): where ξ is the water inrush index; C 1 , C 2 , C 3 , C 4 , C 5 , C 6 are the standardized data of various factors; and ω is the standardized weight in different ranges.

Results
The assessment of water inrush risk consisted of the calculation of the mathematical model, obtainment of the water inrush index, determination of the partition threshold of the water inrush index, and the evaluation of the water inrush risk of the Ordovician limestone floor. The normalized data of the #16 coal seam were substituted into the prediction model to obtain the water inrush index. A total of 49,490 data points were extracted from the contour map of the water inrush index at an interval of 50 m (Table 5). Through an Excel analysis, the distribution of the water inrush index was determined, and a frequency histogram of the water inrush index was established ( Figure 6). According to the distribution law of the frequency histogram, it suddenly decreased at 0.19, 0.31, and 0.56; gradually decreased after 0.46, and the changes of other parts were relatively gentle. The partition threshold of the Ordovician limestone water inrush index of the 16 coal seam floor was determined using the natural breakpoint method, which were 0.19, 0.31, and 0.56. According to this threshold, the mine field was divided into four areas: safe area (<0.19), safer area (0.19-0.31), more dangerous area (0.31-0.56), and dangerous area (>0.56). Figure 7 shows the water inrush risk zoning map. Table 5. Partial water inrush index data (standardization).  At the same time, the thickness of the waterproof layer in this area is the smallest, which further increases the probability of water inrush. The risk of water inrush is the greatest in the northeast and east of the 162 and 163 mining areas. It is suggested that the workers in this mining area need to increase their awareness and be effective in preventing water inrush.   The Chinese geodetic coordinate system 2000 was used in both Figures 7 and 8. It can be seen from Figure 7 that the safety zone in the mining area is relatively large, mainly distributed in the west and southwest, and accounting for about 45% of the mining area. The relatively safe areas are mainly distributed in the central and southeast, accounting for about 20% of the mining area. The dangerous area is relatively large, mainly distributed in the east and northeast, and accounting for about 30% of the mining area. The dangerous area is the smallest, and is distributed in the northeast, accounting for about 5% of the mining area. From the evaluation model, it can be seen that the northeast and the east of the mining area have a greater possibility of water inrush. By analyzing the factors of water inrush from the coal floor, it can be seen that the northeast and east of the mining area have the largest water abundance and water pressure, which provide the source power for the water inrush from the coal floor. Moreover, the maximum mining damage leads to an increase of cracks in the coal floor, which provides a channel for water inrush. At the same time, the thickness of the waterproof layer in this area is the smallest, which further increases the probability of water inrush. The risk of water inrush is the greatest in the northeast and east of the 162 and 163 mining areas. It is suggested that the workers in this mining area need to increase their awareness and be effective in preventing water inrush.

Discussion
The model needed to be tested. In this paper, the evaluation model established based on the water inrush coefficient method [46] (Figure 8) and the actual working situation of the mine were used to test the evaluation results. Chinese scholars first put forward the concept of the water inrush coefficient method during the Jiaozuo Mining Area hydrogeology Conference in 1964, and its expression is in Equation (9): In the formula: Ts-water inrush coefficient, MPa/m; P-aquifer water pressure, MPa; M-thickness of the aquifer from mining coal seam to the floor, M.
It can be seen from the figure that the evaluation results of the model based on the combination weighting method and the results of the model based on the water inrush coefficient method are basically similar. By comparison, the water inrush coefficient method takes too few factors into account, and the range of zoning results obtained is fuzzy, so it cannot accurately reflect the actual situation of the mining area. At the same time, the 16301, 16303, and 16306 working faces in the mining area were mined safely, and no water inrush occurred during the mining period. These working faces are also located in the safety zone of the evaluation results, which further shows that the evaluation results are consistent with the actual situation. Multi-factor evaluation of mining areas, under the coupling of subjective and objective weights, can more specifically reflect the actual situation of a mining area, and the evaluation results are more accurate and scientific, which

Discussion
The model needed to be tested. In this paper, the evaluation model established based on the water inrush coefficient method [46] (Figure 8) and the actual working situation of the mine were used to test the evaluation results. Chinese scholars first put forward the concept of the water inrush coefficient method during the Jiaozuo Mining Area hydrogeology Conference in 1964, and its expression is in Equation (9): In the formula: T s -water inrush coefficient, MPa/m; P-aquifer water pressure, MPa; M-thickness of the aquifer from mining coal seam to the floor, M.
It can be seen from the figure that the evaluation results of the model based on the combination weighting method and the results of the model based on the water inrush coefficient method are basically similar. By comparison, the water inrush coefficient method takes too few factors into account, and the range of zoning results obtained is fuzzy, so it cannot accurately reflect the actual situation of the mining area. At the same time, the 16301, 16303, and 16306 working faces in the mining area were mined safely, and no water inrush occurred during the mining period. These working faces are also located in the safety zone of the evaluation results, which further shows that the evaluation results are consistent with the actual situation. Multi-factor evaluation of mining areas, under the coupling of subjective and objective weights, can more specifically reflect the actual situation of a mining area, and the evaluation results are more accurate and scientific, which gives them a guiding role in the risk evaluation of coal seam water inrush.
At present, as mining has not commenced in the northeast of the mining area, it is impossible to verify the accuracy of the evaluation results in the northeast. If this can be accurately verified in the future, it can be used in future mine safety procedures.

Conclusions
(1) Based on an analysis of the water inrush factors of the #16 coal seam floor in the Jiangzhuang Coal Mine, six main control factors affecting the water inrush of the Ordovician limestone confined water in the floor were determined, including the water pressure of the Ordovician limestone aquifer, water abundance of the Ordovician limestone aquifer, equivalent thickness of effective aquifuge, brittle rock ratio, fragile structure, and depth of failure of the bottom plate. These factors interact with each other, to jointly affect and control the dynamic process of the floor water inrush. (2) IFAHP was introduced to weight the main control factors of water inrush, which not only avoided the cumbersome consistency test of traditional AHP, but also highlighted the advantageous factors in determining the weight and determined the weight of each index more effectively, which provides a new method for the evaluation and prediction of water inrush from a floor. (3) In the calculation process of AHP, the determination of weight depends on expert experience, and the results are easily affected by the subjective factors of experts or evaluators. The entropy weight method determines the weight based on historical water inrush data. The calculation results are relatively objective and can avoid the influence of an evaluator's subjective factors on the weight. However, the weight of the entropy weight method only represents the relative importance of each influencing factor, rather than the actual importance. (4) The IFAHP-EW method combines the entropy weight method and the fuzzy analytic hierarchy process (FAHP), utilizing the respective advantages of the FAHP and the entropy weight methods in determining the weight. This not only uses the actual water inrush data as the basis for calculation, but also reduces the influence of human factors. It effectively combines the subjective and objective, making the evaluation results closer to the actual situation. At the same time, a water inrush risk zoning model was established in combination with MapGIS10.6 software, to more directly and clearly represent the water inrush risk in the mining area. (5) The risk of the Ordovician limestone water inrush in the study area was divided into four grades: dangerous, relatively dangerous, relatively safe, and safe. The evaluation model predicted that the northeast and east of the mining area have the greatest probability of water inrush, and water prevention and control measures should be taken in future coal mining. It was verified that the prediction was in line with the actual situation. Compared with the prediction results of the water inrush coefficient method, it was more detailed and specific, and can be used to guide the prevention and control of mine water. Informed Consent Statement: Not applicable.

Data Availability Statement:
No new data were created or analyzed in this study. Data sharing is not applicable to this article.