Table of Contents Author Guidelines Submit a Manuscript
Volume 2019, Article ID 5485731, 9 pages
Research Article

Research on Water-Filled Source Identification Technology of Coal Seam Floor Based on Multiple Index Factors

1Institute of Resources & Environment, Henan Polytechnic University, Jiaozuo 454000, China
2Collaborative Innovation Center of Coalbed Methane and Shale Gas for Central Plains Economic Region, Jiaozuo, Henan Province 454000, China
3Institute of Resources & Environment, North China University of Water Resources and Electric Power, Zhengzhou, Henan 450000, China
4Energy and Chemical Industry Group of China Pingmei Shenma, Pingdingshan 467000, China
5Shandong Weifang Zongheng Building Materials Company Limited, Weifang, Shandong 262404, China

Correspondence should be addressed to Guang Yang; moc.qq@293067916 and Qi Wang; nc.ude.uwcn@iqgnaw

Received 26 October 2018; Revised 3 January 2019; Accepted 14 February 2019; Published 31 March 2019

Academic Editor: Carlos R. De Souza Filho

Copyright © 2019 Xinyi Wang et al. This is an open access article distributed under the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.


According to the No. 13 Mine of Pingdingshan Coal Co. Ltd., the fault dimension of the study area is calculated based on fractal theory. The four impact factors such as water inflow, water pressure, water inrush coefficient, and fault fractal dimension of 21 water boreholes are used as evaluation indices, and a mathematical model for identifying the water-filled source of the coal seam floor is established by coupling an entropy weight method and fuzzy variable set theory. The model is used to identify the water-filled source of 7 boreholes, which provide a reliable reference for the identification of the water-filled source. According to the calculation result of the entropy weight method, the water volume per unit of the borehole and the groundwater pressure have a significant impact on the water source identification, which accounts for 89.93% of the weight value. In the fuzzy variable set model, when the distance parameter is 1 and the optimization criterion parameter is 2, the accuracy of the water sample identification of the water source category to be identified is 85.71%, which is at a high recognition level. The more typical the impact factors selected and the more samples, the accuracy of water-filled source identification is much higher.

1. Introduction

The Pingdingshan coalfield is located in a North China typical coalfield area. At present, the No. 2-1 coal seams of the Permian Shanxi formation are studied. The coal distribution has 5–7 layers of Carboniferous thin limestone aquifer and a thick layer of Cambrian limestone aquifer. Due to the influence of geological structures and coal mining disturbance, faults and fissures are developed in thin layer of limestone aquifers and aquicludes, which lead to a close hydraulic connection between the thin layer of limestones and the thick layer of limestone aquifers, so the coal seam mining process is threatened by the high-pressure aquifer of the bottom limestone. In order to reduce the influence of water inrush on safe underground production, the key to controlling mine water disasters is to identify the water source quickly and accurately. At present, the main methods for identifying water-filled sources include geological analysis, hydrodynamic analysis, hydrochemical analysis, water temperature analysis, and geophysical prospecting [1].

Liu et al. considered the influence of water temperature and water level combined with a QLT mathematical model to identify the water-filled source in the Panxie mining area [2]. Gui and Lu analyzed the main water inrush and the hydraulic connection between aquifers in the Wanbei mining area by using the radioactive isotope tritium as the discriminant index [3]. Yuan and Gui established the ground temperature equation according to the geothermal characteristics when judging the water source of the Renlou coal mine and achieved good results based on the temperature analysis method [4]. Chen established a three-dimensional (3D) geological model based on the geological structure of the Sunan mining area and identified the water-filled source by simulating the groundwater flow [5]. Wang et al. used six conventional ions as the discriminating factor for the Jiaozuo mining area. Based on the distance discriminant analysis and grey system correlation degree method, the discriminant model of water-filled source identification was established, and the application test was carried out in the Xin’an mining area under similar hydrogeological conditions [6]. Wang et al. analyzed the 6 conventional ions in the groundwater of Pingdingshan coalfield and refined the hydrogeological unit. Based on this, the water source discrimination model was established [7]. Gao carried out hydrogeological pumping and water injection tests in Qiganlou Iron Mine, analyzed its hydrogeological parameters, and predicted the water-filled source of the mine [8]. Wang et al. quantitatively analyzed the structural faults in the Luan mining area and, combined with the improved analytic hierarchy process (AHP), evaluated the dangerous degree of water inrush in the Luan mining area [9]. With the continuous development and improvement of science and technology, a large amount of high-precision equipment has been applied to water source identification. Based on geographic information system technology and drilling depth water temperature fitting results, Ma et al. constructed a model to identify water-filled sources [10]. Zhang used laser-induced fluorescence technology as a means for identifying the water-filled source based on changes in water fluorescence caused by different parameters, such as water temperature and the flow rate of the water source [11]. These achievements provide a reference for future generations to study the water-filled source of mines.

Previous studies have used modern mathematical methods to establish identification models and identify water-filled sources using typical water chemical components, temperature, water level, etc. These methods are simple and fast. However, due to the difficulty of water sample collecting and the external interference of water sample test results, the established identification model is inconsistent with the reality. When geological structure and hydrogeological parameters are used as indicators, they are not quantitatively analyzed. They have space limitations in applicability and cannot be used as a general method to solve problems in other places, which in turn makes it difficult to accurately judge individual water-filled sources.

In this paper, the Carboniferous thin layer and the Cambrian thick limestone aquifer in the coal seam floor of Pingmei No. 13 Mine are used as the object, and the fractal theory is used to quantitatively evaluate the fault complexity of the mining area. Taking groundwater pressure, water inrush coefficient, water inflow per borehole, and fractal dimension of faults as index factors, the weight of each index is calculated based on the entropy weight method, and then the fuzzy variable set theory is used to obtain the comprehensive relative membership degree of the water sample to be discriminated. Finally, the water samples are classified according to the principle of maximum membership degree to identify the water-filled source. The key factors of water inrush are fully considered when the model is built, and the external disturbance is removed to the greatest extent, so as to accurately identify the water-filled source of the coal seam floor.

2. Mathematical Method

2.1. Entropy Weight Method

Determining the weight of each index is the key to calculating the relative membership degree of each water sample. The accuracy of the identification results is directly related to the selection of appropriate methods. The methods for determining weights include AHP, fuzzy inverse equation method, and entropy weight method [12]. Among them, the weight value of the entropy weight method is objective and reliable and can quantify the effective information of each index and provide a new idea for the comprehensive evaluation of each discriminant index. Therefore, this paper uses the entropy weight method to calculate the weight of each index. The calculation steps are as follows [13].

There are kinds of water sources to judge, and the identification indices are . Based on the average value of each identification index , the evaluation matrix is set up. where the entropy value is as follows: When , then .

In the formula,

The weight value is as follows:

2.2. Fuzzy Variable Set Theory

Fuzzy variable set theory is established based on fuzzy set theory. Using the relative membership function, the theory constructs a relative difference function with relativity and dynamic variability, and the concept and model of fuzzy variable set are depicted. This theory expands the static concept and definition of the fuzzy set membership function, overcomes the defect of the static and unique membership function of the fuzzy set, and is widely used in the field of hydrological resources and other engineering examples [1417].

2.2.1. Determination of Relative Membership Degree

Starting with a fuzzy concept (phenomenon or thing) on the domain, let be any element on and satisfy . Let denote the relative membership degree of the attraction property of the fuzzy concept as , and the relative membership degree of the repellent property is . If , it means that the two properties reach dynamic equilibrium; if , it means that equals and attracts the dominant property; if , the opposite is true.

is the relative difference degree of to on the continuous number axis, in the formula .

This is the relative difference function of to.

In addition, the definition of the redundant set of fuzzy variable sets is as follows:

Therefore, the relative membership degree is as follows:

2.2.2. Construction of Relative Differential Function Model

Let be the interval range of on the number axis (Figure 1), i.e., the attraction domain of the fuzzy set ; is an interval within a finite range of ; and are the interval ranges of on the axis , the repulsive domain of the fuzzy set .

Figure 1: Position graph of .

is located in the middle of the interval, when the position of . When any point is on the left of , the relative difference function model is as follows:

When any point is on the right of , the relative difference function model is as follows:

In formulas (8) and (9), is generally used to denote when the model function is a linear function. Formulas (8) and (9) should satisfy (1) when or , then ; (2) when , then ; (3) when or , then ; or (4) when , , .

2.2.3. Comprehensive Relative Membership Calculation

Compared to the standard samples of the level, the comprehensive relative membership degree of the identified sample is . The weight value of each evaluation factor determined by the entropy weight method is . According to equations (4)–(7), the relative membership degree of each standard sample to the evaluation factor can be calculated, and then the comprehensive relative membership degree of each sample to the different standard samples can be obtained according to model (10) [18, 19].

In the above formula, and , where is the number of evaluation factors, is the evaluation factor weight, is the distance parameter, and is the optimization criterion parameter.

2.3. Fractal Theory

Research results at home and abroad show that the distribution of geological faults in the deep part of the strata is a complex and self-similar fractal system [20]. Therefore, the fractal dimension that can quantitatively describe the irregularity of the fractal structure is used as an index for evaluating fault complexity [21, 22].

Fractal theory includes a variety of research methods, in which fault analysis mainly uses the network coverage method to study the fractal structure of geometric objects in a region. Based on this theory, the study is divided into several square lattices with side length . Since there are different levels of cracks and faults inside the fractal, some lattices are empty, and some lattices cover part of the fractal, counting the number of points or lines entering the grid. According to the ratio of 1/2, is reduced, the number of corresponding lattice numbers is counted, and so on. The obtained data are plotted in a double logarithmic coordinate graph , and then the line is fit based on the least squares method. The slope of the fitted line is the fractal dimension . The fractal dimension can be directly related to the complexity of the fault, i.e., the larger the , the more complex the fault.

3. Construction of Index Factors

3.1. Mining Situation

Pingdingshan Coal Co. Ltd. 13 Mine (hereinafter referred to as No. 13 Mine) is located to the northeast of Pingdingshan City, Henan Province. The mine is about 15 km long from East to West and about 4 km wide from North to South, with an area of 54 km2. At present, the No. 2-1 coal seam of the Shanxi formation from the Permian is mainly mined, with an annual output of 2.1 million tons. The northwest boundary of the mine starts from the normal fault of Xingguo Temple, the southeast stops at the Goulifeng normal fault, the northeast is bounded by the normal fault of No. 1 Xiangjia, and the southwest is bounded by the contour line of the -800 m floor of the No. 2-1 coal seam. The minefield is divided into four mining areas, Ji-1, Ji-2, Ji-3, and Ji-4. At present, the Ji-1 and Ji-2 mining areas have been excavated, and the Ji-3 and Ji-4 mining areas are in the process of mining production. The main water hazards in mine excavation come from the Carboniferous limestone aquifer and the Cambrian limestone aquifer under the No. 2-1 coal seam. The thickness of the aquifers is 50-75 m and 90-130 m, respectively, and the groundwater temperature is between 28 and 50°C.

3.2. Complexity of Fault

The study area is divided into several rectangular grids with side length  m according to the longitude and latitude, and they are, respectively, labeled (Figure 2). The similarity ratios , , , and are, respectively, taken (the square grid is further refined into 1, 4, 16, and 64 square grids), and the number of square grids , including the fault traces under different scales, is counted. The fractal dimension of each block is obtained according to the method described above (Table 1).

Figure 2: Grid partition diagram for fractal evaluation of the faults in the No. 13 Mine.
Table 1: Fractal dimension and the correlation coefficient of each block.
3.3. Water Inflow and Water Pressure

The water-filled source is identified based on the geological analysis method and the hydrodynamic analysis method. The data from the 21 boreholes of the No. 13 Mine Cambrian limestone aquifer and Carboniferous limestone aquifer are taken as standard water sample data, and the discriminant indices are water inflow, water pressure, water inrush coefficient, and fault fractal dimension (Table 2). The water-filled Carboniferous aquifer and the Cambrian limestone aquifer are represented by I and II, respectively. The distribution characteristics of the water pressure and water inflow are shown in Figure 3.

Table 2: Discriminant index of water inrush water source.
Figure 3: The contour map of water inflow and water pressure.
3.4. Water Inrush Coefficient

According to the water inrush coefficient of each borehole, the contour map of the water inrush coefficient for the No. 13 Mine water inrush point can be drawn using Surfer software (Figure 4). The water inrush coefficients of the Ji-1 and Ji-3 mining areas in the study area are higher, the water inrush coefficients of the Ji-4 and Ji-2 mining areas are lower, and the overall trend from the northwest to the southeast is gradually increasing.

Figure 4: The contour map of water inrush coefficient.

4. Establishment of Evaluation Model

4.1. Entropy and Weight Calculation

For the 21 water boreholes of the coal seams in the No. 13 Mine, the entropy and weight values of each discriminant are calculated according to equations (2) and (3) (Table 3). According to the calculation results in Table 3, the weights of water inflow and water pressure are 0.6367 and 0.2626, respectively, and the weight ratio of the two indexes is 89.93%, which is much larger than the other two indices, indicating that the water inflow and water pressure have a greater impact on the water-filled source identification in the study area.

Table 3: Calculation results of entropy and weight values.
4.2. Calculation of Comprehensive Relative Membership Degree

A relative difference function model was established based on the 21 water samples of the two types of water sources, and the data points with large deviations in each index data were removed before analysis. The midpoint of the model is the average value of each discriminant index of the water source, i.e., , and the matrix can be obtained. The values of and are calculated based on the mean-standard deviation method. If the standard deviation is , then . When , then make , and matrices and of and can be calculated accordingly.

According to the calculated matrices , , and , the relative position of the index value of the sample in Figure 1 and can be discriminated, and the relative difference degree can be calculated according to a suitable algorithm in formulas (8) and (9) of the index value. Then the relative membership degree can be obtained by substituting the transformed formula (7). Finally, according to formula (10), the comprehensive relative membership degree of each sample for different levels of the standard samples can be obtained, and the water-filled source can be identified according to the principle of maximum membership degree.

5. Discrimination of Water-Filled Source

5.1. Water Source to Be Identified

The identification data of the water burst holes for the 7 water sources to be identified in the No. 13 Mine are shown in Table 4.

Table 4: Parameters of water sample to be identified.

According to the calculation steps listed in Section 4.2, water sample Z1 is taken as an example to explain the process of solving the comprehensive relative membership degree.

As shown in Table 4, the discrimination parameter of water sample is , and is compared with , , and of the standard water sample matrix in formula (10).

According to the position of each parameter value of in the established relative difference function model, the relative difference degree is calculated by using formulas (8) and (9). Since , , and , i.e., . According to this, , and then according to formula (7), the relative membership degree can be calculated, so that we can get the relative membership degree of at all levels. The following matrix is the relative membership degree matrix of water sample Z1.

For formula (10), generally is 1 (Haiming distance) and is 2, and the comprehensive membership degree of water sample Z1 for different aquifers can be calculated [18].

From the above discussion, we find that the comprehensive relative membership degree of water sample Z1 to class I water and class II water is 0.149 and 0.249, respectively. According to the principle of maximum membership degree, water sample Z1 belongs to class II water, which is consistent with the actual situation.

5.2. Discriminant Result Analysis

Referring to the steps in Section 5.1, the water-filled sources of remaining water samples to be identified are analyzed. The comprehensive relative membership degree matrix is as follows. The recognition results are shown in Table 5. According to the results of the table, except for water sample Z6, the water-filled source of other water samples was accurately identified, and the accuracy rate was 85.71%, which belongs to a higher level.

Table 5: Identification results of water inrush.

The weight values of each index calculated by the entropy weight method remain unchanged. Formula (10) uses four different parameter transformation combinations, , ; , ; , ; and , , to calculate the comprehensive relative membership degree of the water sources to be identified in the No. 13 Mine. The results are shown in Table 6.

Table 6: Recognition results under different parameters.

According to Table 6, after changing the parameters of the fuzzy variable set model, the comprehensive relative membership of each water sample for aquifers of different levels changes, but it is basically stable in a small range, indicating that the calculation method used in this paper is more reliable. When and , the recognition accuracy of model parameters is higher than that of other parameters, indicating that the entropy weight-fuzzy variable set method is suitable for water-filled source identification under this parameter combination.

5.3. Discussion

The authors believe that the accuracy of the water-filled source is determined by the representativeness and reliability of the model.

The representativeness of the established water-filled source identification model is related to the typical degree of the selected impact factor. If the impact factors are more typical, the recognition results are more consistent with the actual situation. Therefore, it is necessary to systematically analyze the hydrogeological conditions and mining conditions in the study area and select key identification factors to ensure the representativeness of the model.

The reliability of the model is closely related to the number of samples of the impact factors. The more the number of samples, the more obvious the marker characteristics indicating the change of the impact factor, the easier it is to quantify, and the more reliable it is to identify. Therefore, it is necessary to fully collect and comprehensively monitor the information of various water-filled sources, establish a sound data management system, and lay a foundation for establishing a reliable identification model.

6. Conclusion

(1)Four factors, such as water inflow, water pressure, water inrush coefficient, and the fractal dimension of faults, are taken as the evaluation indices to judge the water-filled source, and the weights of each index are calculated by the entropy weight method. The weights of water inflow and water pressure are 63.67% and 26.26%, respectively, which are much higher than the other two indicators, indicating that they have a greater impact on the identification of the water-filled sources(2)Based on the drilling data of the 21 target strata for the No. 13 Mine reaching the thin layer of the coal seam and the thick layer of limestone, respectively, the fuzzy variable set theory was used to construct a model to identify the water-filled source and applied to the 7 water sources to be identified. The accuracy rate is 85.71%, which provides a strong support for the reliable identification of water-filled sources(3)For the No. 13 Mine, the accuracy of the entropy weight-fuzzy variable set is higher when the distance parameter is 1 and the optimization criterion parameter is 2 in the water-filled source model of the coal seam floor

Data Availability

Data for underlying the findings of the study may therefore be accessed directly from the author or Professor Wang or from the No. 13 Mine of Pingdingshan Coal Co. Ltd.


This study was undertaken as part of a master’s study by the first author.

Conflicts of Interest

The authors declare that they have no conflicts of interest.


This work was supported by the National Natural Science Foundation of China (Grant 41672240, 41802186), the Innovation Scientists and Technicians Troop Construction Projects of Henan Province (Grant CXTD2016053), the Fundamental Research Funds for the Universities of Henan Province (Grant NSFRF1611), and the Henan Province Technological Innovation Team of Colleges and Universities (Grant 15IRTSTHN027). The authors sincerely thank Mr. Wei Zhao and Mr. Pinghua Huang (Henan Polytechnic University) and Mr. Guosheng Chen (Energy and Chemical Industry Group of China Pingmei Shenma) for their valuable assistance.


  1. X. Y. Wang, W. Zhao, X. M. Liu et al., “Identification of water inrush source from coalfield based on entropy weight-fuzzy variable set theory,” Journal of China Coal Society, vol. 42, no. 9, pp. 2433–2439, 2017. View at Google Scholar
  2. W. M. Liu, H. R. Gui, X. F. Sun, J. Luo, and H. Z. Wang, “Discrimination of mine water-filled source in Panxi mining area by QLT method,” China Coal, vol. 27, no. 5, pp. 31–34, 2001. View at Google Scholar
  3. H. R. Gui and Z. W. Lu, “Study on hydrogeological feature of the main pouring water aquifers within the mining area in northern Anhui,” Journal of China Coal Society, vol. 29, no. 3, pp. 323–327, 2004. View at Google Scholar
  4. W. H. Yuan and H. R. Gui, “The characteristics of geothermal temperature and its application in distinguishing the source of water in Ren Lou mine,” Journal of ANHUI University of Science and Technology (Natural Science), vol. 25, no. 4, pp. 9–11, 2005. View at Google Scholar
  5. S. Chen, Evolution and Recognition of Groundwater Based on 3D Geological Structure Model in Sunan Mining Area, China University of Mining and Technology, 2014.
  6. X. Y. Wang, T. Xu, and D. Huang, “Application of distance discriminance in identifying water inrush resource in similar coalmine,” Journal of China Coal Society, vol. 8, pp. 1354–1358, 2011. View at Google Scholar
  7. X. Wang, H. Ji, Q. Wang et al., “Divisions based on groundwater chemical characteristics and discrimination of water inrush sources in the Pingdingshan coalfield,” Environmental Earth Sciences, vol. 75, no. 10, p. 872, 2016. View at Publisher · View at Google Scholar · View at Scopus
  8. S. Gao, Water Source Identification and Prediction of Mine Water Inflow for Xiaoxian Qiganlou Iron Ore, Hefei University of Technology, 2017.
  9. X. Wang, T. Wang, Q. Wang, X. Liu, R. Li, and B. J. Liu, “Evaluation of floor water inrush based on fractal theory and an improved analytic hierarchy process,” Mine Water and the Environment, vol. 36, no. 1, pp. 87–95, 2017. View at Publisher · View at Google Scholar · View at Scopus
  10. L. Ma, J. Z. Qian, and W. D. Zhao, “An approach for quickly identifying water-inrush source of mine based on GIS and groundwater chemistry and temperature,” Coal Geology & Exploration, vol. 2, pp. 49–53, 2014. View at Google Scholar
  11. K. Y. Zhang, Development of coal mine water inrush source identification system based on laser induced fluorescence, AnHui University of Science and Technology, 2016.
  12. T. F. Hu and B. Z. Zhu, “Fuzzy comprehensive evaluation method for complex slope stability based on entropy weight method and analytic hierarchy process,” Railway Engineering, vol. 12, pp. 69–73, 2013. View at Google Scholar
  13. W. D. Gao, “Application of entropy weight fuzzy comprehensive evaluation method in mine water-filled source discrimination,” Mining Safety & Environmental Protection, vol. 39, no. 2, pp. 22–24, 2012. View at Google Scholar
  14. C. Shouyu and F. Guangtao, “A drastic-based fuzzy pattern recognition methodology for groundwater vulnerability evaluation,” Hydrological Sciences Journal, vol. 48, no. 2, pp. 211–220, 2003. View at Publisher · View at Google Scholar · View at Scopus
  15. S. Chen and Z. Hou, “Multicriterion decision making for flood control operations: theory and applications,” Journal of the American Water Resources Association, vol. 40, no. 1, pp. 67–76, 2004. View at Publisher · View at Google Scholar · View at Scopus
  16. S. Y. Chen and X. J. Han, “Engineering method of variable fuzzy set for assessment of surrounding rock stability,” Chinese Journal of Rock Mechanics and Engineering, vol. 25, no. 9, pp. 1857–1861, 2006. View at Google Scholar
  17. S. Y. Chen, “Theory and model of engineering variable fuzzy set—mathematical basis for fuzzy hydrology and water resources,” Journal of Dalian University of Technology, vol. 45, no. 2, pp. 308–312, 2005. View at Google Scholar
  18. J. S. Chen, H. X. He, and T. Wang, “Dam leakage detection based on entropy weight-variable fuzzy set model,” Journal of Hohai University (Natural Sciences), vol. 44, no. 4, pp. 358–363, 2016. View at Google Scholar
  19. S. Y. Chen and M. Li, “Assessment model of water resources reproducible ability based on variable fuzzy set theory,” Journal of Hydraulic Engineering, vol. 37, no. 4, pp. 431–435, 2006. View at Google Scholar
  20. Q. Q. Kang, L. W. Chen, and W. B. Li, “Application of fractal theory in evaluating the complexity of structure in Jinhuashan coal mine,” Shaanxi Coal, vol. 31, no. 1, pp. 65-66, 2012. View at Google Scholar
  21. R. Li, Q. Wang, X. Wang, X. Liu, J. Li, and Y. Zhang, “Relationship analysis of the degree of fault complexity and the water irruption rate, based on fractal theory,” Mine Water and the Environment, vol. 36, no. 1, pp. 18–23, 2017. View at Publisher · View at Google Scholar · View at Scopus
  22. L. X. Fang, S. Y. Lin, F. Zeng, C. H. Wang, and Y. C. Xie, “Research on an alarm method for initial loose parts in a nuclear reactor based on fractal theory,” Nuclear Technology, vol. 195, no. 1, pp. 71–78, 2016. View at Publisher · View at Google Scholar · View at Scopus