Abstract
Water inrush is a serious geological disaster in deep coal mining and rock engineering in China. The occurrence of water inrush is affected by the interaction of many controlling factors in underground reservoir areas. In order to improve the accuracy of water inrush risk evaluation, the changes in index weights due to the change in influencing factors and the uncertain factors should be both considered. In this study, an evaluation model that integrates variable weight and unascertained measure (VWUM) theories was established. First, a modified credible degree recognition was established to make the evaluation results represent both the evaluation grade and the difference in the degree to which the indexes belong to the same grade. Then, the VWUM model was established by coupling the variable weight, unascertained measure, modified credible degree recognition, and geographic information system (GIS). Last, the VWUM model is applied to a case study to evaluate water inrush risk. Results indicates that more detailed risk grades and their spatial distribution zones can be obtained by the VWUM model compared with those obtained by the water inrush coefficient. The VWUM model can provide a more accurate evaluation result.
1. Introduction
The extraction of energy sources such as coal, oil, and natural gas hydrate is a global concern [1, 2]. China is the largest coal producer and consumer [3]. The complex geological conditions in coal mines can cause many disasters, such as roof accidents, gas accidents, coal dust accidents, water inrush, and fire. Mine water inrush is a major threat to mine production safety [4]. The impacts of high temperature, high groundwater pressure, high ground stress, and complex hydro chemical erosion have a great influence on rock damage [5–7], which increase the possibility of water inrush. The risk of water inrush needs to be evaluated to reduce the water inrush accidents. The risk of water inrush can be evaluated by studying the relationship between the water pressure and the properties of the water-resisting strata. Before the middle of the 20th century, Slesarev’s formula was adopted to predict the risk of water inrush from coal seam floor. This period is called Slesarev’s formula period [8]. Z. T. Bieniawski and C. F. Santos studied the mechanism of the coal seam floor, and further studied the water-resisting capacity of the coal seam floor by using the modified Hoek-Brown strength criterion combined with the critical energy release point [9]. In 1964, researchers proposed the water inrush coefficient (WIC) based on the statistical analysis of long-term water inrush data. The water inrush coefficient denotes the ratio of the water pressure to the thickness of the water-resisting strata [10]. This period is called the water inrush coefficient period. At the end of the 20th century, based on the rapid development of hydraulics, some researchers focused on the relationship between stress and seepage and tried to establish a relevant risk evaluation system through a finite element calculation program [11]. With the development of modern computer information technology, various scientific systems are integrated to form a new water inrush prediction method, which can be called the pan-decision analysis theory of water inrush forecast [12]. Pan-decision-making method refers to the fields of engineering and natural science, includes information theory, artificial neural network, vulnerability index method, expert system, and so on. It has high accuracy in the prediction and evaluation of water inrush. For example, Wu Qiang proposed the vulnerability index method. This method creatively coupled geographic information system (GIS) with the artificial neural network, the weight of evidence, the logistic regression, or the analytic hierarchy process to evaluate the risk of floor water inrush [13].
Index weights reflect the relative strength or importance of related factors [14]. Constant weights are usually adopted in the pan-decision-making method. However, index weights will change with the change of index values. Constant weights cannot reflect the change of index values. In order to take index weights to respond to the change in the values of controlling factors of water inrush, variable weight theory has been used in the above-mentioned evaluation methods. For example, Qiang Wu et al. used variable weight model to partition and evaluate the risk of floor water inrush [14]. Wenping Li et al. improved the water inrush vulnerability evaluation model of confined aquifers based on the water inrush coefficient method [10].
The prediction of geological disasters and underground engineering accidents are often vague and uncertain because the influencing factors are often difficult to measure accurately and the interaction of the influencing factors is complex. The evaluation of geological disasters and underground engineering accidents usually considers the uncertain factors. Unascertained measure theory can effectively and quantitatively analyze various uncertain factors and has been widely used in many fields [15–17]. It was first proposed by Guangyuan Wang [18]. Kaidi Liu et al. [19] established unascertained mathematical model to quantitatively describe the unascertained state or the unascertained size for unascertained measure theory. Water inrush is affected by many uncertain factors. The unascertained measure theory can be adopted to analyze uncertainty in the evaluation of water inrush risk. In order to consider the variable weights and the uncertain factors at the same time, the models integrating variable weights theory and unascertained measure theory have been proposed and applied to the evaluation of water inrush risk [20]. The water inrush risk grades can be determined by a credible degree recognition of unascertained measure theory. Some indexes belonging to the same risk grades may have a different degree to which these indexes belong to this grade. However, the risk grades determined by the credible degree recognition cannot reflect this difference. Moreover, this grade cannot be processed as zoning risk maps by GIS because it is an integer. The integer only shows the evaluation grade and cannot show the difference in the same evaluation grade. If decimals are added, the difference can be described. Meanwhile, they can be processed as zoning risk maps by GIS to predict the risk grades of all zones. Therefore, the credible degree recognition should be modified to make the evaluation results represent both the evaluation grade and the difference in degree and be coupled with GIS.
In summary, the evaluation of water inrush risk needs to consider the changes in index weights due to the change in the values of controlling factors and the uncertainty in controlling factors. Because of this, an evaluation model that integrates variable weight and unascertained measure (VWUM) theories was established. In which, the credible degree recognition is modified by defining a new parameter to make the evaluation results represent both the evaluation grade and the difference in the degree to which the indexes belong to the same grade. The parameter can be processed as zoning risk maps by GIS. Then, the VWUM model is applied to a case study to evaluate water inrush risk.
2. Methods
2.1. Unascertained Measure Theory
There are m evaluation objects and n evaluation indexes, which can be expressed as an evaluation matrix . xij denotes the measured value of the i-th evaluation object with respect to the j-th evaluation index. The evaluation grade space is expressed as C = {C1, C1, …, Cp}, where Ck (k = 1, 2, …, p) is the k-th evaluation grade and Ck > Ck+1.
2.1.1. Single-Index Unascertained Measure
The degree to which the measured value xij belongs to the k-th evaluation grade Ck is called the unascertained measure, denoted as μ. Then μ should satisfy the following conditions [15]:
Equations (2)–(4) are called nonnegative boundedness, normalization, and additivity, respectively. The value of μ is determined by the values of evaluation indexes x and threshold of the grades a. The function to describe the relationship of μ, x, and a is called the single-index unascertained measure function. The construction methods of a single-index unascertained measure function mainly include linear, exponential, parabolic, and sinusoidal methods [15]. The linear unascertained measure function is adopted in this study, and its calculation expression is as follows [15]:
The evaluation objects, indexes, and grades form the single-index unascertained measure matrix:
2.1.2. Multi-Index Comprehensive Unascertained Measure
The degree to which the evaluation object belongs to the k-th evaluation grade is determined by the single-index unascertained measure and the index weight , and is calculated by [15]
The vector {μi1, μi1,…, μip,} is called the multi-index comprehensive unascertained measure vector. A credible degree recognition is introduced to determine the evaluation grade of the object, which is expressed as [15]where λ is the credible degree, and is usually set as 0.6 or 0.7. k0 denotes the evaluation object belonging to the evaluation grades Ck0, and it is the minimum of k that satisfies .
2.2. Modified Credible Degree Recognition
The evaluation grade Ck0 can be obtained by the credible degree recognition. However, the evaluation grade Ck0 calculated by the credible degree recognition corresponds to different values of cumulative sum term (i.e., ). The cumulative sum term indicates the degree to which the indexes belong to the same grade. For example, there are two evaluation objects X1 and X2, and their multi-index comprehensive unascertained measure vectors are μ1k = {0.3, 0.35, 0.2, 0.15} and μ2k = {0.3, 0.45, 0.15, 0.1}, respectively. If λ is set as 0.6, k0 is 2 and is 0.65 and 0.75, respectively. The degree to which the X1 belong to C2 is smaller than the degree to which the X2 belongs to C2 since 0.65 < 0.75. However, k0 cannot represent the difference in cumulative sum term of evaluation objects that belong to the same evaluation grade. That is, it cannot represent the difference in the degree to which the indexes belong to same grade. Therefore, the credible degree recognition needs to be modified to make the evaluation results represent both the evaluation grade and the difference in cumulative sum term. Because of this, the credible degree recognition is modified by the flowing steps in this study (see Figure 1). First, the value of cumulative sum term and k are calculated by equation (7). Second, a parameter, named the degree parameter, is introduced to represent the difference in cumulative sum term of evaluation objects that belong to the same evaluation grade, defined aswhere D is the degree parameter. Note that of (7) should be changed to to satisfy D > 0. Last, a parameter, named the grade and degree parameter is introduced to represent both the evaluation grade and the difference in cumulative sum term, defined aswhere G is the grade and degree parameter. The evaluation grade and the difference in cumulative sum term can be obtained by the grade and degree parameter. For example, the grade and degree parameters of the above-mentioned two objects X1 and X2 are 1.125 and 1.375, respectively. The evaluation grades of X1 and X2 are both C2. The degree parameters of X1 and X2 are 0.125 and 0.375, which indicates that the degree to which the X1 belongs to C2 is smaller than the degree to which the X2 belongs to C2.

2.3. Determination of Index Variable Weight
The constant weights of major controlling factors should be determined before the determination of index variable weights. The constant weights can be calculated by the entropy weight method [15, 21, 22]:where Hj is the index entropy value, which denotes the chaos and uncertainty of the evaluation objects. It can be calculated bywhere . In order to avoid , can be expressed as [23]
The variable weights of major controlling factors are calculated by [14]where is the state variable weight vector, which can be calculated bywhere β and u are called punishment level and incentive level, respectively, which values are set as 0.7 and 0.5, respectively [24]. Then, the variable weights can be obtained and expressed as a matrix:
2.4. The VWUM Model for the Evaluation of Water Inrush Risk
Through the above analysis, we can see that: the variable weight can respond to the change of index values; the modified credible degree recognition can represent both the evaluation grade and the difference in cumulative sum term; and the grade and degree parameter G can be analyzed by GIS to obtain the risk grade zoning maps because G contains integers and decimals. Combing the variable weight theory, modified credible degree recognition, and unascertained measure theory, the VWUM model is established in this study. The steps of using the VWUM model and GIS to evaluate the water inrush are as follows: (1) the variable weights are determined; (2) a standard is determined to divide index values into different risk grades; (3) the single-index unascertained measure are calculated by equation (4); (4) the multi-index comprehensive unascertained measure are calculated by putting the variable weights and single-index unascertained measure into equation (6); (5) the grade and degree parameter are determined by the modified credible degree recognition; and (6) The grade and degree parameters show the risk grades of discrete points. The risk grades of all other zones are obtained by the interpolation analysis GIS based on the risk grades of discrete points. Then, the water inrush risk grade zoning map and evaluation results are obtained, which can show both the evaluation grade and the difference in degree in all zones.
3. Case Study
3.1. Study Area
Yangcheng coal mine is located about 17 km northwest of Jining City, Shandong Province. The mining area is located in the monsoon climate zone with four distinct seasons. The average annual precipitation is about 600 mm. The precipitation is concentrated from June to September. The faults are relatively developed and cover the whole mining area. There is no collapse column in the mining area. The target coal seam for the risk evaluation of water inrush is the No. 3 coal seam. Ordovician aquifers are a major threat to the mining safety of the No. 3 coal seam.
3.2. Evaluation Factors (Indexes)
Five major controlling factors are used in this study to evaluate the floor water inrush risk according to geological and hydrogeological conditions of the Yangcheng coal mine.
3.2.1. Coal Seam Mining Depth
Mine pressure increases with the increase of mining depth. The floor is damaged under the impact of rock pressure, and thus, water inrush channels are formed. Therefore, the increase of mining depth will increase the risk of water inrush, which is a main controlling factor for evaluating the risk of water inrush.
3.2.2. Dip Angle of Coal Seam
Dip angles of coal seam cause a water pressure difference between the upper and lower roadways of the working face. The water pressure difference and mine pressure have an impact on mining-induced damage of coal seam floor, and thus changes the risk of water inrush.
3.2.3. Thickness of Water-resisting Strata
Water-resisting strata can block groundwater from entering the working face and prevent the occurrence of water inrush. Therefore, the increase in the thickness of the water-resisting strata reduces the risk of water inrush.
3.2.4. Water Pressure in Confined Aquifers
Water pressure in confined aquifers is the prerequisite for water inrush. The risk of water inrush increases with increasing water pressure.
3.2.5. The Fault Intensity Index
Faults can connect aquifers and coal seams, and increase the possibility of a water inrush. The distribution and amounts of faults can be denoted by the fault intensity index [25]:where FII is the fault intensity index, H is the fault throw, L is the corresponding strike length, S is the area of the grid cell; and n is the number of faults encountered in the grid cell. The risk of water inrush increases with increasing FII.
3.3. Geographic Data of Evaluation Factors (Indexes)
The borehole data is shown in Table 1. The geographic data of each major controlling factor was first input into GIS and then was output as a zoning map. The zoning maps of major controlling factors are shown in Figure 2.

(a)

(b)

(c)

(d)

(e)
3.4. The Evaluation of Water Inrush Risk Using VWUM Model
The original data of boreholes should be normalized before the evaluation of water inrush risk using the VWUM model. The original data of boreholes are normalized bywhere is the data after normalization, is the original data, and are the maximum and minimum of each index values, respectively. The normalization results are shown in Table 2.
3.4.1. The Variable Weights
The constant weights are determined by the entropy weight method. can be calculated by putting the normalization data of boreholes to equation (12). can be calculated by putting to equation (11), then, the constant weights can be obtained by equation (10) (see Table 3).
The state variable weight vectors are calculated by putting the punishment level, the incentive level, and the normalization of original data into equation (14). Then, the variable weights can be obtained by putting the state variable weights and constant weights to equation (13) (see Table 4).
3.4.2. The Index Risk Grade Standards
The index risk is divided into four grades. The index risk grade standards are determined by the K-means clustering method. For example, the threshold values of each group for mining depth, obtained by the K-means clustering method are 449.96, 511.39, 701.7, 733.78, 929.44, and 955.74, respectively. Then, the threshold values of index risk grade standards are 480.675, 717.54, and 942.59, respectively. The index risk grade standards are shown in Table 5.
3.4.3. The Single-Index Unascertained Measure
The single-index unascertained measure of evaluation objects can be obtained by putting the original data of boreholes and index risk grade standard to equation (4). The calculation process of the single-index unascertained measure matrix is as follows by taking borehole number 1st as an example. 750 belongs to Index grades 3. Putting the threshold values of the index risk grade standards for mining depth into equation (4) can obtain the single-index unascertained vector for mining depth, i.e., (0, 0.3467, 0.6533, 0). Then, the single-index unascertained measure matrix for all major controlling factors can be obtained as follows:
3.4.4. The Multi-Index Comprehensive Unascertained Measure
The multi-index comprehensive unascertained measure can be obtained by putting single-index unascertained measure and variable weights to (6) as shown in Table 6.
4. Results and Discussion
The grade and degree parameter can be determined by the modified credible degree recognition (see Table 7).
The grade and degree parameters are put into GIS and processed as a risk grade zoning map (see Figure 3(a)). We can see that water inrush risk is divided into four levels. Risk grades 0-1, 1-2, 2-3, and 3-4 denote extremely safe zone, safe zone, dangerous zone, and extremely dangerous zone, respectively.

(a)

(b)
The water inrush coefficient (WIC) is an important value and can be used to predict the occurrence of water inrush. It was proposed based on the statistical analysis of long-term water inrush data and was defined in the Detailed Regulations for Coal Mine Water Prevention and Control as an empirical formula [10]:where is the WIC (MPa/m), is the water pressure (MPa), and is the thickness of the water-resisting strata (m). The critical water inrush coefficient is 0.1 MPa/m, i.e., water inrush will occur when the coefficient is greater than 0.1 MPa/m, and will not occur when it is smaller than or equal to 0.1 MPa/m. The zoning maps of water inrush risk obtained by the WIS method is shown in Figure 3(b). We can see that the water inrush risk is only divided into two levels, i.e., safety and danger. The results obtained by the WIC method cannot show the degree of safety and danger.
Comparing the evaluation results obtained by the VWUM model and the WIC method, it can be seen that there is the same part in the dangerous zone of the two results, which can validate the evaluation results obtained by the VWUM model. In addition, the more detailed risk grades and their distribution zones can be obtained by the VWUM model compared with those obtained by the WIC method. The more dangerous zones are found by the VWUM model. This is because the VWUM model considers five controlling factors, while the WIC method only considers two factors, i.e., the thickness of water-resisting strata and water pressure. The evaluation results obtained by the VWUM model can show more detailed risk grades and the difference in the degree. Therefore, the VWUM model can improve the accuracy of water inrush risk evaluation.
There might be a better method to modify the credible degree recognition to represent both the evaluation grade and the difference in degree. In addition, the constant weights only consider the objective weights. In some cases, the subjective weights should be also considered.
5. Conclusions
The VWUM model was established to evaluate water inrush risk caused by deep coal mining and rock engineering by coupling the variable weight, unascertained measure, modified credible degree recognition, and GIS in this study. The VWUM model can consider the changes in index weights due to the change in controlling factors and the uncertain factors. The main conclusions are as follows:(1)Compared with the evaluation results obtained by the credible degree recognition, the grade and degree parameter obtained by the modified credible degree recognition can represent both the evaluation grade and the difference in the degree to which the indexes belong to the same grade(2)Compared with the traditional method, WIC, the VWUM model can consider more controlling factors, and the modified credible degree recognition can make the VWUM model be coupled with GIS(3)Compared with the evaluation results obtained by WIC, the water inrush risk obtained by the VWUM model can be divided into more levels
The VWUM model can provide more comprehensive and accurate evaluation results, and thereby help reduce water inrush and increase mining safety.
Data Availability
The data used to support the findings of this study are included within the article.
Conflicts of Interest
The authors declare that there are no conflicts of interest regarding the publication of this paper.
Acknowledgments
This research was funded by the Shandong Higher Education Teaching Reform Research Project of China (2015M138), the National Natural Science Foundation of China (41807211), and Ministry of Education, Higher Education Department, Industry, Education and Research Collaborative Education Project of China (201802197019).