Next Article in Journal
Effects of Flooding Duration and Growing Stage on Soybean Growth Based on a Multi-Year Experiment
Next Article in Special Issue
Application of a Combined FEM/DEM Approach for Teaching a Deep Rock Mass Mechanics Course
Previous Article in Journal
Internal Control, Environmental Uncertainty and Total Factor Productivity of Firms—Evidence from Chinese Capital Market
Previous Article in Special Issue
Mechanical Properties and Energy Characteristics of Flawed Samples with Two Non-Parallel Flaws under Uniaxial Compression
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Study on Height Prediction of Water Flowing Fractured Zone in Deep Mines Based on Weka Platform

1
Department of Mining Engineering, Lyuliang University, Lvliang 033001, China
2
State Key Laboratory of Mining Response and Disaster Prevention and Control in Deep Coal Mines, Anhui University of Science and Technology, Huainan 232001, China
3
State Key Laboratory Cultivation Base for Gas Geology and Gas Control, Henan Polytechnic University, Jiaozuo 454003, China
4
Lvliang Engineering Research Center of Intelligent Coal Mine, Lvliang 033001, China
5
School of Mechanicas Engineering, Shandong University of Science and Technology, Qingdao 266590, China
6
School of Architecture & Civil Engineering, Liaocheng University, Liaocheng 252000, China
*
Author to whom correspondence should be addressed.
Sustainability 2023, 15(1), 737; https://doi.org/10.3390/su15010737
Submission received: 8 December 2022 / Revised: 22 December 2022 / Accepted: 28 December 2022 / Published: 31 December 2022
(This article belongs to the Special Issue Challenges and Strategies for Sustainable Development in Deep Mines)

Abstract

:
Accurately predicting the height of water flowing fractured zone is of great significance to coal mine safety mining. In recent years, most mines in China have entered deep mining. Aiming at the problem that it is difficult to accurately predict the height of water flowing fractured zone under the condition of large mining depth, the mining depth, height mining, inclined length of working face and coefficient of hard rock lithology ratio are selected as the main influencing factors of the height of water flowing fractured zone. The relationship between various factors and the height of water flowing fractured zone is analyzed by SPSS software. Based on the data mining tool Weka platform, Bayesian classifier, artificial neural network and support vector machine model are used to mine and analyze the measured data of water flowing fractured zone, and the detailed accuracy, confusion matrix and node error rate are compared. The results show that, the accuracy rate of instance classification of the three models is greater than 60%. The accuracy of the artificial neural network model is the highest and the node error rate is the lowest. In general, the training effect of the artificial neural network model is the best. By predicting engineering examples, the prediction accuracy of the model reaches 80%, and a good prediction effect is obtained. The height prediction system of water flowing fractured zone is developed based on VB language, which can provide a reference for the prediction of the height failure grade of water flowing fractured zone.

1. Introduction

After the mining of the coal seam, under the action of the mine pressure, the overlying roof strata will move, which will lead to the fissure and fracture of the rock strata. According to the theory of “upper three zones” [1], the water flowing fractured zone refers to the sum of the caving zone and the fracture zone. The fracture channel is easily formed in the water flowing fractured zone. Once these fractures are penetrated, the water-conducting channel will be formed, which will lead to the increase of mine water inflow, which will lead to water inrush and other mine flood accidents, which will seriously threaten the safety of coal mining. Therefore, it is very important for coal mine safety production to accurately predict the development height of water flowing fractured zone.
Domestic and foreign research focus on the development height of water flowing fractured zone after coal mining. In order to ensure the safe mining under the aquifer, foreign countries have carried out long-term research on this. In the 1970s, Britain formulated regulations on coal mining under water bodies; in 1973, Russian researchers also proposed the calculation method of the height of the water flowing fractured zone and formulated the safety regulations for mining under water bodies; Japan has more than ten mines threatened by water damage, developed a special waterproof measure, for the composition and thickness of alluvium developed a safety procedure [2]. At present, the empirical formula given in ‘Regulations of Buildings, Water Bodies, Railways, Main Roadway Coal Pillar Setting and Coal Mining’ [3] is widely used in the calculation formula of water flowing fractured zone in China. These empirical formulas only consider the influence of mining thickness on the height of water flowing fractured zone, and are obtained by regression statistics. For example, the calculation formula of the height Hf of water flowing fractured zone in hard rock is as follows:
H f = 100 M 1.2 M + 2.0 ± 8.9
or
H f = 30 M + 10
In the formula: M is mining thickness.
With the continuous exploitation of coal resources, many mines in China have entered the deep mining stage of coal resources [4,5,6]. For example, the maximum mining depth of Suncun Coal Mine in Xinwen Coalfield has even reached 1350 m [7]. A large number of studies have shown [8,9,10,11,12] that the factors affecting the height of water flowing fractured zone are buried depth of coal seam, inclined length of working face, coefficient of hard rock lithology ratio and mining advance speed. Therefore, it is obviously unreasonable to use the empirical formula to calculate the height of the water flowing fractured zone.
V. Palchik [13,14,15,16] used the method of borehole detection to study the crack development law of overlying strata after coal mining in Donetsk coalfield, Ukraine, and divided the overlying strata into three zones, namely caving zone, fractured zone and continuous deformation zone. Domestic scholars mainly use theoretical analysis [17,18,19,20,21], numerical simulation [22,23,24] and similar material simulation [25,26,27,28], field measurement [29,30] and other methods to study the height of water flowing fractured zone. The calculation method in the ‘three under’ regulation is simple, but only one influencing factor of coal seam mining thickness is considered. Due to the large difference of mine occurrence conditions, it is only suitable for preliminary estimation. The theoretical calculation and numerical simulation method are better than the calculation method in the procedure, but there are shortcomings such as single calculation and simple model. The accuracy of similar material simulation test and field measurement is high, but the workload is heavy, the operation is complex and the cost is high. The prediction method proposed in this paper has comprehensive considerations, simple operation and high prediction accuracy. In this paper, on the basis of previous studies, collected 43 groups of coal seam mining depth greater than 400m of water flowing fractured zone development height of the measured data, using SPSS software to analyze the various factors and water flowing fractured zone relationship. Based on the data mining tool Weka platform, Bayesian classifier, artificial neural network and support vector machine model are used to mine and analyze the measured data of water flowing fractured zone. After comparing and analyzing the three models, the optimal model is obtained and the engineering example is predicted.

2. Modeling

In order to realize the accurate prediction of the height of the water flowing fractured zne [31,32,33], the prediction process is constructed, as shown in Figure 1.

2.1. Selection of Raw Data

By collecting and collating the actual data of the observation of the water flowing fractured zone in China, four influencing factors of mining depth, coefficient of hard rock lithology ratio, height mining and inclined length of working face are selected. Finally, 43 [34] sets of original data of water flowing fractured zone are selected. The first 33 sets of data are used as training samples. The specific data are shown in Table 1, and the last 10 sets of data are used as prediction samples. It can be seen from Table 1 that the mining depth of coal seam is mainly between 400~700 m, the thickness of coal seam is mainly medium thick coal seam and thick coal seam, and the inclined length of working face is between 110~230 m.

2.2. Correlation Analysis

In order to study the development height of water flowing fractured zone under large buried depth (mining depth >400 m), mining depth (X1), coefficient of hard rock lithology ratio (X2), height mining (X3) and inclined length of working face (X4) are selected as the main influencing factors.
(1)
Mining depth
According to the theory of mine pressure control, in a certain range, the greater the mining depth, the greater the mine pressure, mine pressure is proportional to the size and depth of coal mining.
(2)
Coefficient of hard rock lithology ratio
The coefficient of hard rock lithology ratio refers to the ratio of hard rock to statistical height above the roof of coal seam. The hard rock participating in the statistics refers to sandstone, mixed rock and igneous rock. The specific calculation formula is as follows:
b = h ( 15 ~ 20 ) M
In the formula: M is height mining; h is the cumulative thickness of hard rock strata within the height range of the estimated water flowing fractured zone.
(3)
Mining thickness
When the working face advances, periodic pressure will be generated, resulting in roof caving. With the increase of coal seam mining thickness, the plastic zone of overlying strata becomes larger, resulting in the height of caving zone is also larger.
(4)
Inclined length of working face
Before the coal seam is fully mined, the development height of the water flowing fractured zone gradually increases with the mining of the working face; when the coal seam is fully mined, the influence of the inclined length of the working face on the development of the high belt is not obvious.
In order to determine the relationship between the influencing factors and the height of the water flowing fractured zone, each scatter plot was established for research [35,36,37,38,39]. It can be seen from Figure 2 that there is a certain linear relationship between the mining depth (Figure 2a), the coefficient of hard rock lithology ratio (Figure 2b), the height mining (Figure 2c) and the inclined length of working face (Figure 2d) and the height of the water flowing fractured zone (y), as shown below:
y1 = 0.083X1 + 9.009
y2 = 12.103X2 + 44.711
y3 = 2.667X3 + 38.359
y4 = 0.142X4 + 26.698

2.3. Normalization

In order to better retain valid data, it is necessary to reduce the dimension and noise of the raw data, that is, normalization and discretization. The purpose of normalization is to concentrate the values between 0 and 1, and the specific results are shown in Table 2. The normalized calculation formula is as follows:
G i j = X i j min ( X j ) max ( X j ) min ( X j )
In the formula: X i j is the sample before normalization, G i j is the normalized sample, min ( X j ) is the minimum value in the original sample, max ( X j ) is the maximum value in the original sample.
Table 2. Normalization results.
Table 2. Normalization results.
NOX1X2X3X4Y
10.0340.0530.0560.38735.400
20.3110.5530.3750.41554.790
30.2520.6320.3750.15157.450
40.0770.2760.2220.03845.100
50.6770.6051.0000.98176.370
60.0630.6180.1670.87752.010
70.8920.2370.1670.66042.990
80.2610.3030.2920.87749.050
90.6000.7890.2570.15160.140
100.5590.5260.5560.66065.250
110.0340.0390.0560.38735.200
121.0000.5390.0420.60444.540
130.0000.0660.0000.03822.000
140.9490.1840.7921.00070.300
150.1080.6180.7220.49147.550
160.1120.3950.1670.27438.410
170.2970.4080.4170.55743.430
180.1410.4080.2220.03828.630
190.1230.0000.2220.03886.400
200.2170.7500.8060.00022.610
210.0000.0390.0280.18957.490
220.2660.7630.2570.15155.000
230.4080.3950.3750.29286.800
240.1700.8820.8610.50951.400
250.3510.5530.4170.32145.000
260.3150.6180.3060.17945.000
270.0610.1180.1670.27430.290
280.4090.9080.1600.85054.500
290.1130.5390.2220.18945.100
300.1530.0260.3060.74538.810
310.3510.5530.4170.32154.000
320.0580.1450.1670.27432.830
330.5321.0000.0830.60455.320

2.4. Discretization

The discretization is divided into supervised and unsupervised discretization of numerical attributes, which is used to discretize some numerical attributes in the data set to the classification attributes. The ‘mining depth’, ‘coefficient of hard rock lithology ratio’, ‘height mining’ and ‘inclined length of working face’ are equidistantly divided into 3 sections. Similarly, the height of water flowing fractured zone is also divided into 3 sections, 0~40 m is denoted by ‘1’ (water flowing fractured zone height grade ‘low’), 40~60 m is denoted by ‘2’ (water flowing fractured zone height grade ‘medium’), >60 m is denoted by ‘3’ (water flowing fractured zone height grade ‘high’). The calculation formula of discretization is as follows:
L i j = { Q = max ( G j ) min ( G j ) 3 0 , min ( G j ) < G i j < min ( G j ) + Q 1 , min ( G j ) + Q < G i j < min ( G j ) + 2 Q 2 , min ( G j ) + 2 Q < G i j < max ( G j )
In the formula: L i j is the discretized sample, max ( G j ) is the maximum value of the normalized sample data, min ( G j ) is the minimum value of the normalized sample data, Q is the step size.
The repeated data from the discretization results are as follows: the first group and the 11th group of sample data are repeated, the second group and the 17th group of sample data are repeated, the fourth group and the 21st group of sample data are repeated, the 13th group, the 27th group and the 32nd group of sample data are repeated, the 16th group and the 18th group of sample data are repeated, the 25th group and the 31st group of sample data are repeated, and the 26th group and the 29th group of sample data are repeated. The first group, the second group, the fourth group, the thirteenth group, the sixteenth group, the twenty-fifth group, the twenty-sixth group and the twenty-seventh group were removed, and the remaining 25 groups of data were used as training samples, as shown in Table 3.

3. Comparative Analysis of Model Prediction

This paper mainly from the confusion matrix, node error rate, detailed accuracy of the three aspects of comparative analysis.

3.1. Confusion Matrix

The confusion matrix is a special matrix used to show the performance of the algorithm. The larger the diagonal value of the confusion matrix, the more examples of classification. The confusion matrix results of the three models are shown in Table 4. Visualize Classifier errors, in which the instances of correct classification are represented by crosses, and the instances of wrong classification are represented by blocks. Blue indicates a low forecast grade, red indicates a medium forecast grade, green indicates a high forecast grade, The error scatter plots of the three models are shown in Figure 3.
The training samples are a total of 25 sets of data. It can be seen from Table 4 that there are 5 data with ‘low’ height grade of water flowing fractured zone in Naive Bayes model, 3 of which are predicted as ‘medium’; There are 14 data with a height grade of ‘medium’ in the water flowing fractured zone, of which 2 are predicted to be ‘low’ and 1 is predicted to be ‘high’; There are six data with a height grade of ‘high’ in the water flowing fractured zone, of which one is predicted to be ‘low’ and two are predicted to be ‘medium’. There are 16 correct classification examples and 9 wrong classification examples in Naive Bayes model. The correct rate is 64% and the error rate is 36%.
It can be seen from Table 4 that there are 5 data with ‘low’ height grade of water flowing fractured zone in artificial neural network model, of which 1 is predicted as ‘medium’; there are 14 data with a height grade of ‘medium’ in the water flowing fractured zone, of which 2 are predicted to be ‘low’. There are six data with a height grade of ‘high’ in the water flowing fractured zone, of which one is predicted to be ‘low’ and one is predicted to be ‘medium’. There are 20 correct classification examples and 5 wrong classification examples in the artificial neural network model. The correct rate is 80% and the error rate is 20%.
It can be seen from Table 4 that there are 5 data with ‘low’ height grade of water flowing fractured zone in support vector machine model, all of which are predicted as ‘medium’; there are 14 data with a height grade of ‘medium’ in the water flowing fractured zone, and one is predicted to be ‘high’. There are 6 data with a height grade of ‘high’ in the water flowing fractured zone, and 2 are predicted to be ‘medium’. There are 17 correct classification examples and 8 wrong classification examples in the support vector machine model. The correct rate is 68% and the error rate is 32%.
The accuracy of instance classification is: artificial neural network > support vector machine > Naive Bayes.

3.2. Node Error Rate

The node error rate is mainly reflected in the mean absolute error, root mean square error, absolute relative error and root relative square error. The node error rates of the three models are shown in Figure 4.
It can be seen from Figure 4 that the Naive Bayesian model is slightly larger than the support vector machine in the mean absolute error and the absolute relative error. In terms of root mean square error and root relative square error, support vector machine is slightly larger than Naive Bayesian; the value of the node error rate of the two models is not much different. However, it is clear from the diagram that the value of the artificial neural network model is the lowest in the node error rate, and the training effect is the best.

3.3. Detailed Accuracy

The detailed accuracy is mainly reflected in TP Rate (true positive ratio), FP Rate (false positive ratio), Precision, Recall (recall ratio), F-Measure (harmonic average of precision and recall rates), MCC, Kappa statistics and characteristic curve area. The detailed accuracy of the three models is shown in Table 5.
It can be seen from the detailed accuracy that in the TP Rate (TP Rate represents the proportion of the predicted positive class. The higher the value, the higher the accuracy of the positive class prediction.), artificial neural network > support vector machine > naive Bayes; the FP Rate represents the proportion of the negative class contained in the predicted positive class. The smaller the value of the FP Rate is, the better the effect is. The FP Rate of the artificial neural network model is the smallest, and the training effect is the best. The Precision, F-Measure, and MCC values are artificial neural network > naive Bayes > support vector machine. In terms of Recall and Kappa statistics, artificial neural network > support vector machine > naive bayes. In general, the artificial neural network model is optimal.
ROC Area (receiver operating characteristic curve area): Display the ROC area, the decimal range of [0, 1]. The ROC area is generally greater than 0.5, and the closer to 1, the better the classification effect of the model. When the value is between 0.5 and 0.7, the accuracy is low. When the value is between 0.7 and 0.9, it shows a certain accuracy. When the value is greater than 0.9, it shows a higher accuracy. The ROC values of the three models are shown in Table 6, and the ROC curves of the three models are shown in Figure 5.
It can be seen from Table 6 that the average ROC value of artificial neural network model is 0.954, with high accuracy. The ROC values of Naive Bayes model and support vector machine model are between 0.7 and 0.9, with certain accuracy. Overall, Artificial Neural Network > Naive Bayes > Support Vector Machine.

3.4. Prediction Using Artificial Neural Network Models

Through the above comparative analysis, it is concluded that the effect of the artificial neural network model is the best. The artificial neural network model is used to predict the measured data of 10 groups of water flowing fractured zone height. The prediction results are shown in Table 7.
It can be seen from Table 7 that there are 10 groups of prediction data, only 2 groups of data prediction errors. The 7th group of samples to be tested is predicted to be ‘low’ (the actual damage level is ‘high’), and the 8th group of samples to be tested is predicted to be ‘medium’ (the actual damage level is ‘low’). The correct rate of prediction reaches 80%, and good prediction results are obtained.

4. Water Flowing Fractured Zone Height Prediction System Based on VB Language

Based on VB language, the height prediction system of water flowing fractured zone is developed. The software design flow chart is shown in Figure 6.
Some interfaces are shown in Figure 7. Mainly realized the following functions:
(1)
User management module includes user’s department, login account, login password, user name, user role, mobile phone number and notes. User management interface is designed to give managers and ordinary users different permissions, user managers can also be real-time input measured data.
(2)
The factors affecting the height of water flowing fractured zone can be divided into three categories: geological factors, time factors and mining factors. In order to comprehensively collect the measured data of the height of water flowing fractured zone under various influencing factors, the measured data module of the height of water flowing fractured zone is designed, including the name of the mining area, the number of working face and the number of holes, the mining depth, the mining thickness, the dip angle of coal seam, the coefficient of hard rock lithology ratio, the inclined length of working face, the advance speed, the lithology combination structure of overburden rock, the compressive strength of roof, the hardness of coal seam and the measured value.
(3)
Prediction module mainly includes four methods, namely the neural network model, the empirical formula in the ‘three under’ procedure, theoretical calculation formula and multiple linear regression formula.

5. Conclusions

(1)
Based on the Weka platform, 33 groups of measured data of water flowing fractured zone height are used as training samples, and Naive Bayes, artificial neural network and support vector machine models are used for training respectively. The prediction results are compared and analyzed by confusion matrix, node error rate and detailed accuracy. It is concluded that the artificial neural network model is the best.
(2)
Based on artificial neural network, taking mining depth, coefficient of hard rock lithology ratio, height mining and inclined length of working face as the main factors, the accuracy of the height prediction model of water flowing fractured zone is 80%, which has achieved good prediction effect and can provide reference for the prediction of height failure grade of water flowing fractured zone.
(3)
Based on VB computer language, the height prediction system of water flowing fractured zone is developed. The system is simple to operate, practical and convenient, and realizes the height prediction of water flowing fractured zone under different methods.

Author Contributions

All the authors contributed to publishing this paper. L.B. and C.W. conceived the main idea of the paper; C.L. contributed to the theoretical analysis; F.M. analyzed the data; C.W. and M.Z. wrote the paper; M.F. and B.Z. modified the figures and proofread the revised version. All authors have read and agreed to the published version of the manuscript.

Funding

This paper was supported by the High-level Talent Scientific Research Launch Fund of Anhui University of Technology (No. 2021yjrc17), Shandong Provincial Natural Science Foundation (No. ZR2022ME165), Basic Research Program of Shanxi Province (20210302124633); Lvliang platform base construction project (2021GCZX-1-46); Key Laboratory of Gas Geology and Gas Control in Henan Province-Open Fund Project of National Key Laboratory Cultivation Base (WS2021B07); Open Fund Project of Key Laboratory of Safety and High-efficiency Coal Mining, Ministry of Education (Anhui University of Science and Technology) (JYBSYS2021101).

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Data Availability Statement

All the data in this paper are available from the corresponding author upon request.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Bai, L.; Liao, C.; Wang, C.; Zhang, M.; Meng, F.; Fan, M.; Zhang, B. Overburden Migration Law of Inclined Coal Seam and Determination of “Upper Three Zones” Height in Tengda Coal Mine. J. Min. Strat. Control Eng. 2022, 4, 70–79. [Google Scholar]
  2. Shengyou, L. Study on the Overburden Movement Principle and Integrated Water-Sand Prevention and Control Technology of Fully-mechanized Top Coal Caving under Soft Overburden and Multi-Aquifers; China University of Mining & Technology: Beijing, China, 2017. [Google Scholar]
  3. State Bureau of Coal Industry, Regulations of Buildings, Water Bodies, Railways, Main Roadway Coal Pillar Setting and Coal Mining; China Coal Industry Publishing House: Beijing, China, 2000.
  4. Hu, Y.; Li, W.; Wang, Q.; Chen, X.; Zheng, G. Evaluation Method of Water Hazard Control Effect of Coal Seam Floor in Deep Mining: Sequence Verification Evaluation. Geofluids 2022, 2022, 6728045. [Google Scholar] [CrossRef]
  5. Zhang, G.; Wang, Z.; Guo, G.; Wei, W.; Wang, F.; Zhong, L.; Gong, Y. Study on Regional Strata Movement during Deep Mining of Erdos Coal Field and Its Control. Int. J. Environ. Res. Public Health 2022, 19, 14902. [Google Scholar] [CrossRef] [PubMed]
  6. Pang, D.; Niu, X.; He, K.; Li, C.; Pang, D.; Hu, T.; Chen, Z.; Luo, X. Study on the Deformation Mechanism of the Bottom Plate along the Empty Lane of Deep Mining and the Control Technology of the Bottom Drum. Geofluids 2022, 2022, 3429063. [Google Scholar] [CrossRef]
  7. Shi, L.Q.; Xin, H.Q.; Zhai, P.H.; Li, S.C.; Liu, T.B.; Yan, Y.; Wei, W.X. Calculating the Height of Water Flowing Fracture Zone in Deep Mining. J. China Univ. Min. Technol. 2012, 41, 37–41. [Google Scholar]
  8. Xu, X.; Wang, X.; Li, Y. Prediction of the Height of Water Flowing Fractured Zone Based on the MPSO-BP Neural Network Model. Math. Probl. Eng. 2022, 2022, 2133695. [Google Scholar] [CrossRef]
  9. Yin, H.; Dong, F.; Zhang, Y.; Cheng, W.; Zhai, P.; Ren, X.; Liu, Z.; Zhai, Y.; Li, X. Height Prediction and 3D Visualization of Mining-Induced Water-Conducting Fracture Zone in Western Ordos Basin Based on a Multi-Factor Regression Analysis. Energies 2022, 15, 3850. [Google Scholar] [CrossRef]
  10. Hou, E.; Wen, Q.; Ye, Z.; Chen, W.; Wei, J. Height prediction of water-flowing fracture zone with a genetic-algorithm support-vector-machine method. Int. J. Coal Sci. Technol. 2020, 7, 740–751. [Google Scholar] [CrossRef]
  11. Dai, S.; Han, B.; Liu, S.; Li, N.; Geng, F.; Hou, X. Neural network–based prediction methods for height of water-flowing fractured zone caused by underground coal mining. Arab J Geosci 2020, 13, 495. [Google Scholar] [CrossRef]
  12. Liu, Y.; Yuan, S.; Yang, B.; Liu, J.; Ye, Z. Predicting the height of the water-conducting fractured zone using multiple regression analysis and GIS. Environ. Earth Sci. 2019, 78, 422. [Google Scholar] [CrossRef]
  13. Palchik, V. Localization of mining-induced horizontal fractures along rock layer interfaces in overburden: Field measurements and prediction. Environ. Geol. 2005, 48, 68–80. [Google Scholar] [CrossRef]
  14. Palchik, V. Formation of fractured zones in overburden due to longwall mining. Environ. Geol. 2003, 44, 28–38. [Google Scholar] [CrossRef]
  15. Palchik, V. Influence of physical characteristics of weak rock mass on height of caved zone over abandoned subsurface coal mines. Environ. Geol. 2002, 42, 92–101. [Google Scholar] [CrossRef]
  16. Palchik, V. Prediction of hollows in abandoned underground workings at shallow depth. Geotech. Geol. Eng. 2000, 18, 39–51. [Google Scholar] [CrossRef]
  17. Tan, Y.; Cheng, H.; Lv, W.; Yan, W.; Guo, W.; Zhang, Y.; Qi, T.; Yin, D.; Wei, S.; Ren, J.; et al. Calculation of the Height of the Water-Conducting Fracture Zone Based on the Analysis of Critical Fracturing of Overlying Strata. Sustainability 2022, 14, 5221. [Google Scholar] [CrossRef]
  18. Chengcheng, C. Investigation on the Height Calculation of Water Flowing Fracture Zone during Coal Mining under Ordovician Limestone Nappe Aquifer. Adv. Civ. Eng. 2022, 2022, 6997343. [Google Scholar]
  19. Dong, Y.; Huang, Y.; Du, J.; Zhao, F. Study on Overburden Stability and Development Height of Water Flowing Fractured Zone in Roadway Mining with Cemented Backfill. Shock Vib. 2021, 2021, 6661168. [Google Scholar] [CrossRef]
  20. Guo, W.; Zhao, G.; Lou, G.; Wang, S. A New Method of Predicting the Height of the Fractured Water-Conducting Zone Due to High-Intensity Longwall Coal Mining in China. Rock Mech. Rock Eng. 2019, 52, 2789–2802. [Google Scholar] [CrossRef]
  21. Zhang, Y.; Cao, S.; Wan, T.; Wang, J. Field Measurement and Mechanical Analysis of Height of the Water Flowing Fracture Zone in Short-Wall Block Backfill Mining beneath the Aquifer: A Case Study in China. Geofluids 2018, 2018, 7873682. [Google Scholar] [CrossRef] [Green Version]
  22. Fan, H.; Wang, L.; Lu, Y.; Li, Z.; Li, W.; Wang, K. Height of water-conducting fractured zone in a coal seam overlain by thin bedrock and thick clay layer: A case study from the Sanyuan coal mine in North China. Environ. Earth. Sci. 2020, 79, 125. [Google Scholar] [CrossRef]
  23. Zhu, T.; Li, W.; Wang, Q.; Hu, Y.; Fan, K.; Du, J. Study on the Height of the Mining-Induced Water-Conducting Fracture Zone Under the Q 2l Loess Cover of the Jurassic Coal Seam in Northern Shaanxi, China. Mine Water Environ 2020, 39, 57–67. [Google Scholar] [CrossRef]
  24. Zhang, W.; Li, B.; Zhang, G.; Li, Z. Investigation of Water-Flow Fracture Zone Height in Fully Mechanized Cave Mining Beneath Thick Alluvium. Geotech. Geol. Eng. 2017, 35, 1745–1753. [Google Scholar] [CrossRef]
  25. He, C.; Lu, W.; Zha, W.; Wang, F. A geomechanical method for predicting the height of a water-flowing fractured zone in a layered overburden of longwall coal mining. Int. J. Rock Mech. Min. Sci. 2021, 143, 104798. [Google Scholar] [CrossRef]
  26. Hu, T.; Hou, G.; Bu, S.; Zhu, Z.; Wang, Y.; Hu, Z.; Li, Z. A Novel Approach for Predicting the Height of Water-Conducting Fracture Zone under the High Overburden Caving Strength Based on Optimized Processes. Processes 2020, 8, 950. [Google Scholar] [CrossRef]
  27. Ren, Z.; Wang, N. The Overburden Strata Caving Characteristics and Height Determination of Water Conducting Fracture Zone in Fully Mechanized Caving Mining of Extra Thick Coal Seam. Geotech. Geol. Eng. Int. J. 2020, 38, 329–341. [Google Scholar] [CrossRef]
  28. Wang, F.; Xu, J.; Chen, S.; Ren, M. Method to Predict the Height of the Water Conducting Fractured Zone Based on Bearing Structures in the Overlying Strata. Mine Water Environ. 2019, 38, 767–779. [Google Scholar] [CrossRef]
  29. Gao, W.; Li, Y.; He, Q. Determination of Fractured Water-Conducting Zone Height Based on Microseismic Monitoring: A Case Study in Weiqiang Coalmine, Shaanxi, China. Sustainability 2022, 14, 8385. [Google Scholar] [CrossRef]
  30. Feng, J.; Wang, S.; Hou, E.; Ding, X.; Duan, H. Determining the Height of Water-Flowing Fractured Zone in Bedrock-Soil Layer in a Jurassic Coalfield in Northern Shaanxi, China. Adv. Civ. Eng. 2021, 2021, 9718802. [Google Scholar] [CrossRef]
  31. Wang, S.; Xia, Q.; Xu, F. Investigation of collector mixtures on the flotation dynamics of low-rank coal. Fuel 2022, 327, 125171. [Google Scholar] [CrossRef]
  32. Yao, D.; Jiang, N.; Wang, X.; Jia, X.; Lv, K. Mechanical behaviour and failure characteristics of rocks with composite defects of different angle fissures around hole. Bull. Eng. Geol. Environ. 2022, 81, 290. [Google Scholar] [CrossRef]
  33. Ji, S.; Wang, Z.; Karlovšek, J. Analytical study of subcritical crack growth under mode I loading to estimate the roof durability in underground excavation. Int. J. Min. Sci. Technol. 2022, 32, 375–385. [Google Scholar] [CrossRef]
  34. Shi, L.Q.; Wu, H.B.; Li, Y.L.; Lv, W.K. Optimization Model of PCA-GA-Elman for Development Height Prediction of Water-Conducting Fissure Zone. J. Henan Polytech. Univ. (Nat. Sci.) 2021, 40, 10–18. [Google Scholar]
  35. Ji, S.; Karlovšek, J. Calibration and uniqueness analysis of microparameters for DEM cohesive granular material. Int. J. Min. Sci. Technol. 2022, 32, 121–136. [Google Scholar] [CrossRef]
  36. Pan, H.; Jiang, N.; Gao, Z.; Liang, X.; Yin, D. Simulation study on the mechanical properties and failure characteristics of rocks with double holes and fractures. Geomech. Eng. 2022, 30, 93–105. [Google Scholar]
  37. Cui, W.; Chen, J. Insight into mineral flotation fundamentals through the DFT method. Int. J. Min. Sci. Technol. 2021, 31, 983–994. [Google Scholar] [CrossRef]
  38. Feng, F.; Chen, S.; Wang, Y.; Huang, W.; Han, Z. Cracking mechanism and strength criteria evaluation of granite affected by intermediate principal stresses subjected to unloading stress state. Int. J. Rock Mech. Min. 2021, 143, 104783. [Google Scholar] [CrossRef]
  39. Chen, S.J.; Feng, F.; Wang, Y.J.; Li, D.Y.; Huang, W.P.; Zhao, X.D.; Jiang, N. Tunnel failure in hard rock with multiple weak planes due to excavation unloading of in-situ stress. J. Cent. South Univ. 2020, 27, 2864–2882. [Google Scholar] [CrossRef]
Figure 1. Prediction procedure.
Figure 1. Prediction procedure.
Sustainability 15 00737 g001
Figure 2. Relationship between influencing factors and height of water flowing fractured zone.
Figure 2. Relationship between influencing factors and height of water flowing fractured zone.
Sustainability 15 00737 g002
Figure 3. Error scatter plot.
Figure 3. Error scatter plot.
Sustainability 15 00737 g003
Figure 4. Node error rate.
Figure 4. Node error rate.
Sustainability 15 00737 g004
Figure 5. ROC.
Figure 5. ROC.
Sustainability 15 00737 g005aSustainability 15 00737 g005b
Figure 6. Process design of the software.
Figure 6. Process design of the software.
Sustainability 15 00737 g006
Figure 7. System interface.
Figure 7. System interface.
Sustainability 15 00737 g007aSustainability 15 00737 g007b
Table 1. Raw data.
Table 1. Raw data.
Serial NumberMining Depth s/mCoefficient of Hard Rock Lithology Ratio bHeight Mining M/mInclined Length of Working Face l/mHeight of Water Flowing Fractured Zone H/m
1412.40.092.215735.4
24890.474.516054.79
3472.50.534.513257.45
4424.420.263.412045.1
55900.51922076.37
6420.50.52320952.01
7649.10.23318642.99
8475.20.283.920949.05
9568.60.653.6513260.14
10557.250.455.818665.25
11412.550.082.215735.2
126790.462.118044.54
13403.20.11.812022
146650.197.522270.3
154330.52716847.55
16434.10.35314538.41
174850.364.817543.43
18441.970.363.412028.63
19437.170.053.412086.4
204630.627.611622.61
21403.10.08213657.49
22476.40.633.6513255
23515.70.354.514786.8
244500.72817051.4
25499.90.474.815045
264900.52413545
27420.060.14314530.29
285160.742.95206.154.5
29434.40.463.413645.1
30445.40.07419538.81
31499.920.474.815054
32419.030.16314532.83
335500.812.418055.32
Table 3. Training samples.
Table 3. Training samples.
NOX1X2X3X4Y
112212
232333
312132
431122
511132
623113
722223
811121
932122
1031333
1112322
1212222
1312111
1411113
1513311
1611112
1713112
1822213
1913322
2023132
2112112
2211131
2322212
2411111
2523122
Table 4. Confusion matrix.
Table 4. Confusion matrix.
Naive BayesArtificial Neural NetworkSupport Vector Machine
abcClassified AsabcClassified AsabcClassified As
230a = 1410a = 1050a = 1
2111b = 22120b = 20131b = 2
123c = 3114c = 3024c = 3
Table 5. Detailed accuracy.
Table 5. Detailed accuracy.
ModelTP RateFP RatePrecisionRecallF-MeasureMCCKappa
Naive Bayes0.640.2970.6450.640.6350.3670.3608
Artificial neural network0.80.1320.8340.80.8050.680.6622
Support vector machine0.680.3690.5560.680.6030.360.3651
Table 6. ROC value of the models.
Table 6. ROC value of the models.
Naive BayesArtificial Neural NetworkSupport Vector MachineClass
0.860.950.751
0.8150.9450.6462
0.8290.9780.8513
Average weight0.8270.9540.716
Table 7. Comparison of prediction results.
Table 7. Comparison of prediction results.
Serial NumberMining Depth/mCoefficient of Hard Rock Lithology RatioHeight Mining/mInclined Length of Working Face/mHeight of Water Flowing Fractured Zone/mHeight Grade of Water Flowing Fractured ZoneComparison
ActualPrediction
14200.713.77056.8mediummediumcorrect
2478.30.543.8520952.15mediummediumcorrect
3568.40.852.94180.457mediummediumcorrect
4453.60.16419544.96mediummediumcorrect
5412.50.242.213635.2lowlowcorrect
6411.70.32.213635.21lowlowcorrect
74750.376.117064.6highlowerror
8453.50.421.618030.3lowmediumerror
9818.50.457.5523074.57highhighcorrect
10427.30.434.612056.6mediummediumcorrect
Disclaimer/Publisher’s Note: The statements, opinions and data contained in all publications are solely those of the individual author(s) and contributor(s) and not of MDPI and/or the editor(s). MDPI and/or the editor(s) disclaim responsibility for any injury to people or property resulting from any ideas, methods, instructions or products referred to in the content.

Share and Cite

MDPI and ACS Style

Bai, L.; Liao, C.; Wang, C.; Zhang, M.; Meng, F.; Fan, M.; Zhang, B. Study on Height Prediction of Water Flowing Fractured Zone in Deep Mines Based on Weka Platform. Sustainability 2023, 15, 737. https://doi.org/10.3390/su15010737

AMA Style

Bai L, Liao C, Wang C, Zhang M, Meng F, Fan M, Zhang B. Study on Height Prediction of Water Flowing Fractured Zone in Deep Mines Based on Weka Platform. Sustainability. 2023; 15(1):737. https://doi.org/10.3390/su15010737

Chicago/Turabian Style

Bai, Liyang, Changlong Liao, Changxiang Wang, Meng Zhang, Fanbao Meng, Mingjin Fan, and Baoliang Zhang. 2023. "Study on Height Prediction of Water Flowing Fractured Zone in Deep Mines Based on Weka Platform" Sustainability 15, no. 1: 737. https://doi.org/10.3390/su15010737

APA Style

Bai, L., Liao, C., Wang, C., Zhang, M., Meng, F., Fan, M., & Zhang, B. (2023). Study on Height Prediction of Water Flowing Fractured Zone in Deep Mines Based on Weka Platform. Sustainability, 15(1), 737. https://doi.org/10.3390/su15010737

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop