Water Richness Zoning and Evaluation of the Coal Seam Roof Aquifer Based on AHP and Multisource Geological Information Fusion

Roof water disaster is one of the most serious disasters in the process of coal mine safety mining in China. The thickness of modern comprehensive mechanized coal mining is large, which has a great impact on the roof aquifer, and the threat degree of water disaster is also increasing. Therefore, the evaluation of the water richness of the roof aquifer is an important work for coal mine water prevention and control. In order to systematically evaluate the water-rich property of an aquifer in a coal seam roof, this paper couples multisource geological information based on an analytic hierarchy process (AHP). The characteristics of the roof aquifer, hydrogeology, and structure are determined to be the main control factors of K 2 limestone water richness evaluation. Under the main control factors, seven independent water richness evaluation indexes, including aquifer depth, aquifer thickness, borehole water level, borehole water consumption, faults, folds, and collapse columns, are divided, and a normalized water richness coupling evaluation model is constructed, which realizes the zoning and evaluation of a water-rich aquifer in a coal seam roof.


Introduction
Mine water disaster is a major threat to coal mine safety production, among which water inrush from the roof aquifer is one of the main types of coal mine water disaster. With the popularization of modern comprehensive mechanized coal mining, in order to make full use of coal resources to the maximum extent, the mining thickness of coal seam increases, which leads to the increase in the influence range of the caving zone and water-conducting fracture zone of the coal seam roof, the damage to the roof aquifer increases, and the threat degree of water disaster increases accordingly. The water richness of the aquifer is one of the important influencing factors of mine roof water disaster intensity [1][2][3][4]. Before coal mining, it can systematically and accurately evaluate the water richness of the aquifer and predict the water inrush threat of the aquifer during mining, which will accurately guide the coal mine water prevention and control work and reduce the probability of roof water inrush and water damage.
Generally, the strength of water richness of the aquifer can be visually characterized by unit water inflow [5], but there are great differences in coal seam occurrence conditions and hydrogeological conditions in different mining areas, so it is obviously not comprehensive to rely solely on unit water inflow. The distribution law of water richness of the aquifer in the coal mine roof is controlled by the comprehensive influence of multiple geoscience factors, including the original stratum background factors, complex hydrogeology, tectonic geology, and other control mechanisms as well as linear and nonlinear characteristics. The evaluation of water richness of the aquifer is a systematic evaluation of multisource geological information factors and artificial mining effects [6][7][8]. Based on this, many scholars use various methods to evaluate the water richness of the aquifer. Wu et al. [9][10][11] put forward a GIS-based information fusion aquifer water richness evaluation method, which was named "water-rich index method", and realized the multiangle and all-round evaluation of water richness of the aquifer. Lv et al. [12] propose the "multitype four-double method" evaluation and prediction technology, which improved the accuracy and sensitivity of evaluation and overcome the deficiency of the empirical formula. Bu et al. [13] completed the water richness zoning of the roof by using an information fitting method and made fuzzy prediction on this basis. Qiu et al. [14] put forward a comprehensive grey correlation method, fuzzy Delphi analytic hierarchy process, and geophysical exploration-combined evaluation and prediction method of Ordovician limestone water richness. Han et al. [15] introduced a set pair analysis-variable fuzzy set coupling evaluation method to predict and evaluate the water richness of the coal seam roof aquifer and constructed the aquifer water richness evaluation model of the Zhiluo Formation aquifer based on the set pair analysis-variable fuzzy set coupling evaluation method. Hou et al. [16] aiming at the problem of water richness prediction of weathered bedrock established a Fisher discriminant analysis model of water richness of the aquifer type of the weathered bedrock aquifer by using the measured unit water inflow and index data of the borehole. Wei et al. [17] combined the improved entropy weight method and analytic hierarchy process by introducing the distance function, obtained the comprehensive weight of the main control factors of water inrush, and constructed the calculation model of the water inrush risk index. Xue [18] put forward the water richness index method based on fractal theory to evaluate aquifer water richness and gave the evaluation process. Compared with the traditional water richness evaluation method, this method reduces the number of structural description elements by 2~3 items and solves the problem of reasonable and accurate evaluation and zoning of aquifer water richness under the condition of low hydrogeological exploration. Guo et al. [19] divided the water richness grade of this mining area into three grades by the FAHP method and ArcGIS as the platform, combined with the GRA grey correlation degree, and determined the threshold value of water richness classification objectively.
In this paper, through comprehensive analysis of multisource geological information such as aquifer thickness, borehole water level, borehole water consumption, faults, folds, and collapse columns of K 2 limestone in the direct roof aquifer of the no. 15 coal in the lower group of the coal seam in a mine field, AHP is used to carry out multifactor coupling analysis, and a normalized evaluation model of water richness of the aquifer is established, finally realizing the zoning evaluation of water richness of the aquifer in the roof of the mine field.

Overview of the Study Area
The study area is located in the deep circulating runoff area of karst water in the Niangziguan spring area and is divided into six aquifer rock groups from bottom to top, which are the limestone karst fissure aquifer group of the middle Ordovician Shangmajiagou Formation, limestone karst fissure aquifer group of the middle Ordovician Fengfeng Formation, karst and sandstone fissure aquifer group of the upper Carboniferous Taiyuan Formation, sandstone fissure aquifer group of the lower Permian Shanxi Formation, sandstone fissure aquifer group of the Permian, and pore aquifer group of the Quaternary glutenite layer. Among them, the K 2 limestone aquifer in the karst and sandstone fissure aquifer group of the Taiyuan Formation of the Upper Carboniferous System is the direct roof aquifer for mining no. 15 coal in the future mining field, which poses a great flood threat to the mining of no. 15 coal. According to statistics, K 2 limestone water gushed during the construction of an inlet and return air shaft from 2009 to 2014, and the water gushed from different locations ranging from 5 to 102.30 m 3 /h. From 2015 to 2020, 10 K 2 limestone water exploration holes were constructed, of which 3 holes had water gushing about 1 m 3 /h and the other water exploration holes had no water outflow. According to the actual exposure, it shows that the K 2 limestone aquifer of the Taiyuan Formation in the mining field is characterized by uneven water richness and abnormal local water richness.

Analysis of the Main Controlling Factors of
Water Richness in the K 2 Limestone Karst Aquifer For the K 2 limestone karst aquifer, the water richness of the limestone aquifer is largely affected by the development of karst space in limestone. Therefore, the factors that affect the karst development degree are also the factors that determine the water richness of the karst aquifer. In this paper, based on the factors affecting karst development and combined with the specific conditions of the study area, the aquifer characteristics, hydrogeological characteristics, and structural characteristics are selected as the main control factors to study the water richness of the karst aquifer.
3.1. Aquifer Characteristics. K 2 limestone water mainly occurs in karst space, so the properties of K 2 limestone and the development characteristics of karst fissures will greatly affect the aquifer characteristics. First of all, its nature is the thickness of the limestone aquifer. Under normal circumstances, a thicker aquifer has the potential to store more groundwater, and under the condition of satisfying karst development, large dissolution caves will be formed. Another characteristic of the limestone aquifer is the buried depth of limestone. The greater the buried depth, the slower the groundwater circulation, the smaller the dissolution power provided by groundwater movement, and the worse the dissolution degree. When the buried depth of K 2 limestone is greater, its water storage space is compressed under the action of the self-weight stress of overlying strata, which reduces the water richness degree of limestone.

Hydrogeological
Characteristics. Hydrogeological characteristics are an important factor affecting limestone water richness. Hydrogeological parameters, such as permeability coefficient, directly reflect K 2 limestone water abundance, while groundwater level reflects runoff conditions in the mining field. Another measurable indicator is water consumption of drilling holes, which can well reflect the location and degree of karst fracture development. The greater water consumption of drilling holes indicates fracture development in this horizon. Under the condition of karst development, 2 Geofluids there are abundant dissolution fractures in this layer, which makes the water richness degree of this layer good.
3.3. Structure Characteristics. Faults, folds, karst collapse columns, and other structures are mainly developed in the mine field. The structural development degree greatly affects the stability of K 2 limestone development, and it is easy to form a flow channel for groundwater in limestone. Under the condition of satisfying karst development, a larger karst water storage space will be formed. Karst development is the result of multiple factors such as lithology, structure, hydrogeological conditions, and landform, but in fact, topography and hydrogeological characteristics in karst areas are controlled by structure, lithology is the material basis, hydrology and landform are important conditions for karst formation, and structure is the main controlling factor for karst development, which is described as follows: (1) Faults provide the migration channel of groundwater, and the rocks in the fault zone are generally broken and have good water conductivity under the condition of loose cementation, and the development degree of faults is positively correlated with the water abundance of K 2 limestone (2) The rock mass in the collapse column is broken, the stress is released, and the in situ stress redistribution causes the local stress concentration around the collapse column, which leads to the development of fracture zones with a certain width around and above the collapse column. The thickness of the effective water-resisting layer is reduced, which generally leads to multiple aquifers. In addition, the crushing strength of the column is reduced, and it is easy to shear under limestone water pressure and water infiltration. The pores and fissures formed by column breakage become the channels of groundwater seepage. On the whole, collapse column development and karst development promote each other (3) The influence of the development degree of the fold structure on the karst development degree is mainly reflected in the fact that a large number of tensile joints are associated with the fold crankshaft and the turning end, which develop in the roof and floor strata of the coal seam, forming a good water storage space for the K 2 limestone aquifer. And some fissures are also channels for groundwater storage and migration, which will form a large karst space under suitable conditions in the long process of geological evolution. In general, the fold structure promotes karst development

Establishment of the Water Richness Evaluation Model
Aquifer characteristics, hydrogeological characteristics, and structural characteristics are determined as the three main control factors to control K 2 limestone water richness evalu-ation, and seven independent factors such as aquifer thickness, buried depth, borehole water level, borehole water consumption, faults, collapse columns, and folds are established under the three main control factors. Geoinformation thematic maps of independent factors are established, respectively, and geoinformation data are analyzed by spatial fusion based on the analytic hierarchy process [10].

Quantitative Expression Method of Each
Factor. By analyzing drilling data of K 2 limestone, the buried depth value of the K 2 limestone roof, thickness value of K 2 limestone, and initial water level value of the aquifer are counted. According to a simple hydrological chart, the maximum water consumption of the drilling hole in K 2 limestone is counted and the degree of structural development is quantitatively evaluated by fractal dimension theory [20].
In this study, the study area is divided into 350 square units of 500 m × 500 m (among which 117 grids have structures). Every 500 m long grid is subdivided into 250 m, 125 m, and 62.5 m square grids, and the number of grids passing through each level of the construction trace is counted one by one, so as to obtain the Nðr i Þ values when R = 250 m, 125 m, and 62.5 m, which are put into the lg NðrÞ − lg r coordinate system to generate fitting straight lines, and their slopes can be calculated by the least square method.
After calculating one by one, assign the absolute value of fractal dimension D of all grids to the center point of the corresponding grid unit and then interpolate one by one to make the contour map of fractal dimension in the whole region.

Determination of Weights among Factors.
In this study, an analytic hierarchy process was used to determine the weight value of each independent factor in the K 2 limestone aquifer water richness evaluation model. Based on the basic principle of the analytic hierarchy process, firstly, the hierarchical structure analysis model is established, and the judgment matrix is constructed. Then, the importance of each influencing factor is scored by the expert scoring system, and then, the consistency test is carried out, and finally the weights among each influencing factor are obtained [21][22][23].
According to the analysis of multivariate information reflecting aquifer water abundance, it can be divided into three levels. The ultimate goal of this problem is to evaluate the water abundance of the K 2 limestone aquifer, which is taken as the target layer (level A) of the model. Aquifer characteristics, hydrogeological characteristics, and structural characteristics reflect the water richness degree of the aquifer, but its influence mode needs to be reflected by specific factors related to it, which is the intermediate link to solve the problem. As the criterion layer (level B) of the model, each specific multivariate geoscience information constitutes the decisionmaking level (level C) of this model, and the required solution can be finally achieved by making decisions on the 3 Geofluids problems at this level. See Figure 1 for the analysis structure model.

Construct the Judgment Matrix.
According to the analysis of the factors affecting water inrush in the Xinyuan Mine, using the method of "collecting experts to score," the opinions of field experts and researchers in universities and scientific research units are collected and consulted and their opinions are widely listened to. According to their personal experience and experience in field production practice and scientific research and specific methods to deal with problems, many factors affecting water inrush are scored. The scoring standard is to list the proposed factors affecting water inrush according to the scale method of 1-9 created by Saaty (the scale and meaning of the judgment matrix are shown in Table 1) [24] and invite domain experts to evaluate the relative importance of each factor according to the role of each factor in numerous water inrush accidents and their personal experience in production practice and scientific research. The quantitative score of each factor is given. According to the final cumulative score, the total scores of each factor are compared, and the expert evaluation set of each influencing factor is formed, thus constructing the judgment matrix of AHP evaluation of roof and floor water inrush risk.

Hierarchical Single Sorting and Consistency Check.
According to the judgment matrix, the weights of single ordering of each layer are calculated as shown in Tables 2-5.
The above judgment matrix CR is less than 0.1, which has satisfactory consistency and can pass the consistency test.

Hierarchical Total Sorting and Consistency
Test. The general ranking of local academic information in aquifer water richness analysis is shown in  Figure 1: Hierarchical analysis structure model for water richness evaluation of the K 2 limestone aquifer. Compared with the two factors, the former is slightly more important than the latter 5 Compared with the two factors, the former is obviously more important than the latter 7 Compared with two factors, the former is strongly more important than the latter 9 Compared with two factors, the former is extremely more important than the latter 2, 4, 6, 8 The intermediate value of the above adjacent judgments

Reciprocal
If the ratio of importance of factor i to factor j is a ij , then the ratio of importance of factor j to factor i is a ji = 1/a ij   index Ci is relative to total target A, and WA/Ci is the weight of each index Ci to total target A. Calculate and get the random consistency ratio of the total sorting of layer C: Therefore, CR < 0:1 has satisfactory consistency, and W A/Ci can be used as the final decision-making basis, so as to determine the weight value superimposed by multiple information in the study of water richness zoning of water-filled aquifers ( Table 7).

Establishment of the Water Richness Evaluation Model.
The water richness index is the sum of the superposition effects of various independent information in a certain area or a certain position, and it is a normalized measure of the water richness degree of the aquifer [10,11].
In fact, the evaluation model of K 2 limestone water richness of the aquifer is to establish a mathematical model to express the comprehensive influence degree of various factors. The calculated value of the model can reflect the water richness of the aquifer in a certain geographical position. The establishment of the preliminary model must be based on the analysis of geological conditions, water richness factors, and the contribution mechanism of each factor to water inrush.
The initial model of water richness index CI is introduced to evaluate the water richness of the aquifer. It can be expressed by the following model formula: where CI is the water richness index, W k is the weight of geoscience information, f k ðx, yÞ is the normalized single information influence value function, x, y is the geographical coordinate, and n is the number of multivariate information. f k ðx, yÞ is the normalized value of the kth geoscience information value in the evaluation of the water-rich aquifer of the coal seam roof. The weight values of each influencing factor are W 1 = 0:111, W 2 = 0:055, W 3 = 0:064, W 4 = 0:256, W 5 = 0:282, W 6 = 0:108, and W 7 = 0:124, respectively. Therefore, it can be concluded that the water richness evaluation model of the coal seam roof aquifer is

Evaluation of Water Richness Zoning of the K 2 Aquifer
The determined structural fractal dimension data coordinates are cast into the normalized thematic map of aquifer characteristics and the normalized thematic map of hydrogeological characteristics, and the subitems of aquifer characteristic and hydrogeological characteristic evaluation are read out, respectively, and recorded in the structural fractal dimension value list. The normalized fractal dimension value constructed at that same coordinate point and the normalized values of aquifer characteristics and hydrogeological characteristics are multiplied by corresponding weight values and then added to obtain the comprehensive evaluation index table of water richness. Import the final data list into the analysis and drawing software, and get the comprehensive evaluation map of water richness.

Thematic Map of Aquifer Characteristics
(1) Thematic Map of K 2 Aquifer Thickness. On the basis of fully collecting and sorting out the previous drilling data, the thickness and buried depth of limestone above K 2 or K 2 are counted, respectively, and the thematic map of K 2 limestone aquifer thickness is drawn by using the powerful interpolation function of Surfer, which is a Kriging interpolation method (shown in Figure 2). It can be seen from Figure 2 that the thickness of the K 2 limestone aquifer changes greatly and the thickness of the K 2 limestone aquifer generally becomes thinner from northwest to southeast, and K 2 limestone is missing locally. The thickness of K 2 limestone is mostly distributed in 1-3 m, and the thickness of limestone in some areas is more than 3 m, among which the thickest limestone area is 12 m.
(2) Thematic Map of K 2 Limestone Burial Depth. According to the data of K 2 limestone exposed by drilling statistics, a thematic map of K 2 limestone floor buried depth is drawn,   Figure 3. It can be seen from the thematic map of limestone floor buried depth that K 2 limestone is buried deeper in the north than in the south. The maximum buried depth data is 890.31 m. The thickness of limestone decreases from north to south, and the buried depth is mainly between 300 m and 500 m.

Thematic Map of Hydrogeological Characteristics
(1) K 2 Limestone Water Level Thematic Map. According to drilling data, analyze the corresponding K 2 limestone water level value and draw the contour map of K 2 limestone initial water level, as shown in Figure 4. It can be seen from the figure that the initial water level of K 2 limestone is high in the middle, low in the northwest and southeast, and the water level in the southwest of the mine field is the highest.
(2) Thematic Map of Water Consumption of the K 2 Limestone Borehole. For the statistical data of coal seam drainage holes, mine hydrogeology holes, underground exploration holes, etc., study the water consumption data of K 2 limestone drilling holes and draw a thematic map of K 2 limestone drilling water consumption (as shown in Figure 5). The water con-sumption in the central part of the mining field is large, and the karst fissure development degree is high and the groundwater mobility is good. On the whole, the karst fissure development in this area is good, and the groundwater conduction is smooth.

Structural Thematic Map
(1) Thematic Map of Fault Fractal Dimension Value. The fractal dimension value of fault structure in the mining field is shown in Figure 6. According to the thematic map analysis of fault fractal dimension value, it can be seen that the fault fractal dimension value decreases from northeast to southwest. The fault fractal dimension value in the whole mining field is 1.35 at the maximum (38417648.7005, 4196650.361) and 0 at the minimum, so the fault does not develop in this area. Generally speaking, the fault development degree is high in the east and northeast of the mine field, but undeveloped in the west.
(2) Thematic Map of Fractal Dimension Value of the Collapse Column. The fractal dimension of the collapse column in the mining field is shown in Figure 7. By analyzing the thematic     12 Geofluids map of fractal dimension value of the collapse column, it can be seen that the maximum fractal dimension value of the collapse column is 1.585 and the minimum fractal dimension value is 0. On the whole, the development degree of the collapse column is high in the middle of the mining field and low in the northeast.
(3) Thematic Map of Fold Fractal Dimension Value. The fractal dimension of fold in the mining field is shown in Figure 8. By analyzing the thematic map of fold fractal dimension, it can be seen that the fold development degree in the north of the mine field is higher than that in the south.

Establishment of the Normalized Thematic Map.
The water richness index model of the coal seam roof aquifer is actually a mathematical model to express the comprehensive influence of various geological information. The calculation results of this model can reflect the water richness degree of the water-filled aquifer in a certain area. The establishment of the initial model must be based on various geological information to reflect the water richness of the aquifer. Because a single thematic map contains only one geoscience information, it cannot meet the requirements of multiinformation comprehensive processing through digital models. Therefore, an initial model must be built first, which can basically reflect the comprehensive effect of local learning information, and then, through repeated parameter adjust-ment and fitting operation, it can gradually approach the target, and finally a model that can reflect the actual situation of water richness of the aquifer in the coal seam roof is established. Therefore, in order to fuse geoscience information data, it is necessary to establish a normalized theme map of each factor.

Data Normalization.
The evaluation adopts linear function transformation, and the expression is as follows: 5.2.2. Single-Factor Normalized Thematic Map. After the single-factor data is normalized, each single-factor attribute database is established. Using the normalized data at the same borehole coordinate position, each single-factor normalized thematic map is made as shown in Figure 9.

Comprehensive Evaluation
Chart of Water Richness. By constructing grid and drilling distribution coordinate position as shown in Figure 10, the encrypted normalized thematic map isoline is projected to read out the normalized thematic map values of each factor at the same position. According to the water richness evaluation model established above, the comprehensive evaluation of water richness is shown in Figure 11. It can be seen from the figure that the  13 Geofluids water richness property in the northwest of the mining field is poor, while the water richness property in the middle and northeast is good, especially in the middle.

Water Richness Zoning
Map. Through the mathematical statistical analysis of the spatial data of the influence degree of water richness, combined with the water inflow data related to K 2 limestone drilling in the past, water richness indexes of 0.12, 0.2, and 0.3 are the basis of division. According to the water richness index, the whole area is divided into four areas with different relative water richness degrees, and the water richness areas are in order from large to small: weak water richness areas, medium water richness areas, strong water richness areas, and extremely strong water richness areas. See Figure 12 for the water richness zoning map of K 2 limestone. It can be seen from the figure that the water  14 Geofluids richness property in the middle, southeast, and central part of the mine field is strong~extremely strong, which is the focus of mine water disaster prevention and control work.

Conclusions
(1) Based on the factors affecting karst development and combined with the specific conditions of the study area, three main control factors are determined, including aquifer characteristics, hydrogeological characteristics, and structural characteristics, and seven water richness evaluation indexes are divided, namely, K 2 limestone buried depth, K 2 limestone aquifer thickness, borehole water level, borehole water consumption, faults, folds, and collapse columns (2) The weight of seven independent influencing factors is assigned by AHP, and the water richness evaluation model of the normalized K 2 limestone aquifer is constructed (3) The geoinformation thematic map of each independent factor is established, and on this basis, the com-prehensive evaluation map of water richness is established (4) The water richness evaluation results show that the water richness property in northwest of the mine field is poor, while the water richness property in the middle and northeast of the mine field is better, especially in the middle part. It provides scientific guidance for mine water disaster prevention and control

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

Conflicts of Interest
The author declares that there is no conflict of interest regarding the publication of this paper.