Next Article in Journal
Emulation Strategies and Economic Dispatch for Inverter-Based Renewable Generation under VSG Control Participating in Multiple Temporal Frequency Control
Next Article in Special Issue
Simulation of Mining-Induced Ground Damage Using Orthogonal Experiments to Determine Key Parameters of Super-Large Coalface: A Case Study in Shendong Coalfield in China
Previous Article in Journal
Numerical Investigation of Erosion Wear in the Hydraulic Amplifier of the Deflector Jet Servo Valve
Previous Article in Special Issue
Numerical Parametric Study of Countermeasures to Alleviate the Tunnel Excavation Effects on an Existing Tunnel in a Shallow-Buried Environment near a Slope
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Application of Artificial Neural Networks in Assessing Mining Subsidence Risk

1
Department of Earth Resources and Environmental Engineering, Hanyang University, Seoul 04763, Korea
2
AICONT, Seoul 08378, Korea
*
Author to whom correspondence should be addressed.
Appl. Sci. 2020, 10(4), 1302; https://doi.org/10.3390/app10041302
Submission received: 6 December 2019 / Revised: 30 January 2020 / Accepted: 11 February 2020 / Published: 14 February 2020
(This article belongs to the Special Issue Land Subsidence: Monitoring, Prediction and Modeling)

Abstract

:
Subsidence at abandoned mines sometimes causes destruction of local areas and casualties. This paper proposes a mine subsidence risk index and establishes a subsidence risk grade based on two separate analyses of A and B to predict the occurrence of subsidence at an abandoned mine. For the analyses, 227 locations were ultimately selected at 15 abandoned coal mines and 22 abandoned mines of other types (i.e., gold, silver, and metal mines). Analysis A predicts whether subsidence is likely using an artificial neural network. Analysis B assesses a mine subsidence risk index that indicates the extent of risk of subsidence. Results of both analyses are utilized to assign a subsidence risk grade to each ground location investigated. To check the model’s reliability, a new dataset of 22 locations was selected from five other abandoned mines; the subsidence risk grade results were compared with those of the actual ground conditions. The resulting correct prediction percentage for 13 subsidence locations of the abandoned mines was 83–86%. To improve reliability of the subsidence risk, much more subsidence data with greater variations in ground conditions is required, and various types of analyses by numerical and empirical approaches, etc. need to be combined.

Graphical Abstract

1. Introduction

Various types of subsidence can occur at abandoned mines, especially when proper post-mining measures, such as filling the mine cavern with infilling material, are not implemented to prevent future collapse. Subsidence sometimes leads to significant consequences, including devastating ground settlement resulting in damage to the local area and also human casualties. A significant obstacle in handling subsidence at abandoned mines is predicting when and where the ground movement will occur, which significantly hinders development in the area. Accurate predictions of dynamic surface movement and deformation are exceptionally difficult because mining activity is a space–time process and dynamic surface movement is related to a wide range of factors [1].
Accordingly, many attempts have been made to more precisely predict subsidence and validate the stability of mining and land areas by numerical analyses [1,2,3,4,5], theoretical and experimental solutions [6,7,8,9,10], and/or statistical approaches [11,12], but the results have been limited because abandoned mines are generally too large and the goafs are too complex for effective analysis using the aforementioned approaches.
Against this background, the use of neural networks in fields of engineering has become prevalent in the last few years [13,14,15,16]. An artificial neural network (ANN) is a computational mechanism capable of acquiring, representing, and computing a mapping from one multivariate space of information to another, given a set of data representing that mapping [17].
ANN provides tools for optimising operations, equipment selection, and problems involving large amounts of information that humans cannot easily assimilate in the process of decision-making. Hence, the neural network can serve as a tool for determining the relative importance of the factors influencing the stability of underground objects [18].
Since the 1990s, when a number of ANN-based systems were introduced in the mining industry [19], ANN has significantly improved the current approach towards predicting abandoned mine subsidence risk [20,21,22,23]. In particular, Oh et al. [24] evaluated the performance of predictive Bayesian, functional, and meta-ensemble machine learning models in generating land subsidence susceptibility maps. Kim et al. [25] also constructed a hazard map for possible ground subsidence around abandoned underground coal mines in Korea using an ANN together with a geographic information system (GIS). Zhao and Chen [26] used an ANN to predict ground subsidence at a metal mine because ground subsidence is influenced by many factors which are fuzzy and nonlinear, and the ANN has a great function to be able to handle such nonlinear problems, without knowing specific mine conditions.
Meanwhile, Mine Reclamation Corporation (MIRECO) who is responsible for the safety and management of abandoned mines and implements mine reclamation projects in Korea has made several attempts to reliably predict subsidence by numerical, statistical, probabilistic, and/or model approaches [27,28,29,30,31] prior to this study. However, the prediction results were not up to par with MIRECO’s standards because they were not appropriate for the assessment of many as well as extensive abandoned mines.
Finally, MIRECO decided to adopt ANN as well to predict subsidence considering the merits of ANN mentioned above. Hence, the main objective of this study is to predict the possibility of subsidence on any ground under which a goaf exists by subsidence risk grade that represents the likelihood of subsidence using ANN. This study can be used as an initial method screening the specific locations that are at highest risk of subsidence and thus require immediate detailed investigation, from an extensive mine area for subsequent detailed analyses by the other approaches.
Meanwhile, from a civil engineering point of view, mining induced ground subsidence is described by parameters measuring ground deformation i.e., radius of curvature, tilt, and vertical as well as horizontal ground deformation [32,33]. The subsidence in this study signifies the big downward movement of ground’s surface enough for one to be able to easily recognize, and Figure 1 shows several examples of the subsidence.

2. Methodology

The overall study process is summarized in Figure 2 and is classified into four categories; collection and preprocessing of subsidence data, construction of ANN model and training, establishment of subsidence risk grade, and test of the model with new data.

2.1. Basic Principle of ANN

The capability of an ANN model to optimize, determine, and predict a given situation depends heavily on the reliability of the database it draws from as well as the amount of data it has processed. This is because an ANN system derives its output from patterns recognized during the pre-learning stage so that a more reliable outcome can be acquired from a model that has learned from a more extensive variety of situations.
Figure 3 shows a model of the artificial neuron or processing element (PE) forming the basis of the ANN structure. This layered structure is the most common structure in ANNs and is usually called the fully connected feedforward or acyclic network. The starting point of the ANN structure is a layer of input units that allows entrance of information into the network. Normalization is the process of equalizing the signal range of different inputs. Normalization ensures that changes in the signals of different inputs have the same effect on the network’s behavior regardless of their magnitude [19].
Learning from examples is the main operation of any ANN. In this case, ‘learning’ means the ability of an ANN to improve its performance through an interactive process of adjusting its free parameters. The adjustment of an ANN’s free parameters is stimulated by a set of examples presented to the network during the application of a set of well-defined rules, called a learning algorithm, for improving its performance. There are many different learning algorithms for ANNs, each with different ways of adjusting the synaptic weights of PEs and different ways of formalizing the measurement of the ANN’s performance [19].
In general, there are three major learning algorithms: supervised learning, unsupervised learning, and reinforcement learning [34]. In supervised learning, the network is trained by providing it with input and matching output patterns. These input output pairs can be provided by an external teacher or by the system that contains the neural network (self-supervised). This study adopted supervised learning. Furthermore, there are two different approaches: error back-propagation algorithms, a gradient method, and a genetic algorithm, which is a stochastic search method [35]. The back-propagation training algorithm is the most frequently used neural network method [36,37,38]. It is trained using a set of examples of associated input and output values.
After the ANN reaches the required performance by learning examples, it can be used to predict the outcome of situations given specific input parameters. It merely behaves as a series of functions that produce an output for a particular input.

2.2. ANN Software for this Study

The software used for this study was the Generalized Data Analyzer and Predictor (GDAP, Ver. 2.0, 2005) developed by the University of Wales Swansea in Wales, which has been used for prediction, supervised learning, and design of structure [39]. Figure 4 shows a typical example of GDAP training, which consists of the learning rate, architecture of the hidden layer, training algorithm, normalization scheme, etc.

3. Information on/about Abandoned Mines and Datasets

The subsidence data for this study were collected from 467 locations (one location means one subsidence point), among which subsidence had occurred at 388 locations, and there were 79 locations where subsidence had not occurred (Table 1). However, more than a few data points did not have any information on the goafs—e.g., the depth from ground to goaf, goaf height/width, and rock mass quality (RMR)—which has an essential impact on subsidence, or the information was very unclear. Such data were excluded from the analyses.
Eventually, 227 locations from 37 abandoned mines were selected for the study. Among these, 118 locations were from 15 abandoned coal mines and 109 locations were from other types of mines (i.e., gold, silver, metal). The numbers of locations that did and did not display subsidence were 166 and 61, respectively.
The locations of the selected abandoned mines are shown in Figure 5.
Furthermore, 207 of the total 227 locations were used for machine learning, and the other 20 were used for validation of the prediction model during the machine learning process. The information on the mines and subsidence locations was collected from various mine reports published by the Mine Reclamation Corporation in Korea. The reports present the results of extensive investigations on subsidence conducted at various abandoned mines across the country.
Meanwhile, many previous studies have identified important factors that contribute to ground subsidence around mines, including the depth and height of the mined cavities, excavation method, degree of inclination of the excavation, scope of mining, structural geology, and flow of groundwater [40,41]. Salmi et al. [3] found that gradual deterioration due to weathering plays a significant role in the propagation of the cavity and the formation of sinkholes above abandoned mine workings. Tzampoglou and Loupasakis [4] have analyzed the causes of 4 km long land subsidence around the mine causing damage to villages of Greece using the PLAXIS three-dimensional finite element code, and concluded that the factors which affect land subsidence (from the most dangerous to the least dangerous) are the groundwater drawdown, the geotechnical behavior of the formation, the tectonic structure, and the offset of the faults. Manekar and Chaudhari [2] mentioned that subsidence is a time-dependent process and controlled by various factors including methods of mining, depth of extraction, thickness of deposit, and topography, as well as the in-situ properties of the rock mass. Yu et al. [9] concluded that factors of the relationship between surface dynamic subsidence and overburden separated strata include the requirement of abscission layer development (stress, stiffness, mining condition, and deflection) and surface subsidence conditions (mining thickness, mining depth, mining overburden rock properties, speed). Wang et al. [10] developed a new time function model to predict the dynamic subsidence that occurs in the grout-injected overburden of a coal mine during overburden grout injection mining. This function includes the mining parameters (the width of the working face and the advancement length of the filled area), a geological parameter (the angle of full subsidence), grout injection parameters, and other relevant coefficients. Li et al. [1] mentioned that factors related to subsidence are mining and geological conditions, mining methods, working face advance speeds, and mining period, amongst others.
Based on the literature study regarding the factors affecting subsidence, the following 11 potential influence factors for subsidence were identified: ① elapsed time in years from the closing of the mine until an investigation was performed; ② RMR; ③ groundwater level; ④ existence of an outcrop at the investigated location; ⑤ types of loads near the investigated location (bridge, railway, house, etc.); ⑥ depth from ground to the nearest goaf; ⑦ goaf height; ⑧ horizontal distance from the investigated location to the nearest goaf; ⑨ uniaxial compressive strength (UCS) of the rock mass; ⑩ rainfall of the local area; and ⑪ occurrence or lack of occurrence of subsidence.
As the second step of data pre-processing, correlation analyses between subsidence occurrence and the potential influential factors were conducted and the multicollinearity among the factors was assessed. Multicollinearity is when two or more of the independent variables being fed into the model are highly correlated. In regression models, this may cause regression coefficients to change erratically in response to small changes in the model or the data because of confounding between the influences of highly correlated factors. Thus, if the highly correlated input variables are identified and eliminated in advance, then the network can determine the most appropriate weights more efficiently [42].
The results showed that the correlation between subsidence occurrence and UCS, or rainfall was very low and that the multicollinearity between RMR and UCS, as well as between groundwater level and rainfall, was high (Figure 6 and Figure 7). These mean that the effects of UCS and rainfall on subsidence are comparatively low, and RMR and groundwater level are closely related to UCS and rainfall respectively. Finally, UCS and rainfall were excluded from the factors for the analyses, and the other nine influence factors were selected, as depicted in Figure 8.
Figure 9 presents the value ranges of the selected data for each factor. The time elapsed from the closing of the mine until an investigation was performed ranged between 7 and 23 years, with an average of 10 years, whereas the time elapsed for the other types of mines ranged between 0 and 71 years, with an average of 25 years. The average depth from ground level to the nearest goaf at a subsidence location was 36 m for coal mines and 41 m for other types of mines. Furthermore, the average horizontal distance was approximately 7 m and 5 m, respectively. The number of investigated locations near a railway or road was 15 and 14, corresponding to about 12% and 13% of the total investigated locations (118 and 108 locations), respectively. The number of locations located on an outcrop was 20 and 25, which corresponds to about 24% and 18% of the total investigated locations, respectively.

4. Analysis and Discussion

The analysis of this study consisted of two stages: analysis A and analysis B. Analysis A sought to predict whether or not subsidence would occur using an ANN. Analysis B established a mine subsidence risk (MSR) index using an interaction matrix to indicate the extent of risk for subsidence at any given location. The results of both analyses were integrated into a single classification system referred to as the subsidence risk grade, which categorizes the location based on risk level.

4.1. Analysis A

4.1.1. Normalization of the Data

Normalization is a technique often applied as a part of data preparation for machine learning. The goal of normalization is to ensure that every data point has the same scale so that each factor is rated as equally important. Hence, it is only required when features have different ranges. In this study, normalization was used to compare the contribution of each identified influence factor towards the risk of subsidence and was performed by scaling the values of each factor from 0 to 10 points.

4.1.2. Structure of the Model

For optimum learning and prediction of an ANN, the structure of the ANN model and initial conditions in the ANN program need to be established first. Thus, a series of analyses were performed for various cases, varying the number of hidden layers and nodes to determine the optimum structure of the ANN for this study. For this purpose, the same dataset was applied for each case, and the case with minimized error % was adopted as the optimum structure.
Finally, it was discovered that the optimum structure occurred when the number of nodes for the input, first hidden, second hidden, and output layer was 8, 10, 10, and 1 for coal mines and 8, 9, 9, and 1 for other types of mines, respectively, as depicted in Figure 10.

4.1.3. Setting of the Training Parameters

In the learning and prediction stage of the ANN development, the learning rate in the GDAP program was set to 0.1, and the initial weights were randomly selected between 0.1 and 0.9. As a training algorithm for this analysis, a resilient propagation algorithm that is a variant of the standard back-propagation neural network was used because it leads to much faster training than others and minimizes the error between the predicted output values and the calculated output values. The algorithm propagates the error backwards and iteratively adjusts the weights. The maximum number of iterations and maximum system errors to determine when training should be terminated were set to 20,000 and 0.000001, respectively, applying ‘early stopping’ to avoid overfitting on the training dataset.

4.1.4. Results of Analysis A

After training the model, the accuracy of the model was measured separately for the training datasets (total 207), and validation datasets (total 20) to check for overfitting of the model. The results showed that the accuracies for coal mines for training datasets and validation datasets were 97% and 83%, respectively, and were 97% and 75% for other type mines, respectively, as presented in Table 2. Logically, because the training was performed using the training datasets, their accuracy was higher than the accuracy for the validation datasets. It is expected that, once additional training datasets are added, in the future, the accuracy for validation will also improve, reducing the overfitting on the training dataset.

4.2. Analysis B

For analysis B, the MSR index was assessed using the eight factors, which were divided into two groups: a direct influence factor group and an indirect influence factor group. The direct influence group consisted of three factors, depth to a goaf, goaf height, and horizontal distance to the nearest goaf, which directly influence subsidence. The indirect influence group included the other five factors: elapsed time, RMR, load type, outcrop, and groundwater level.

4.2.1. Direct Influence Factor Group

Jung et al. [43] suggested a method to evaluate abandoned mine subsidence risk using only the depth to a gangway, which can be collected easily from existing materials. Many other studies also placed much importance on the depth to a gangway or goaf as a major factor inducing subsidence [30,44,45,46,47].
Accordingly, it can be stated that the depth to a goaf is a vital factor to the occurrence of subsidence and that the magnitude of subsidence is greatly affected by goaf height and the horizontal distance to a goaf. Thus, the depth to goaf was classified into 11 groups by using the depth range to allocate the points between 10 and 100, and this was adjusted by goaf height and distance (Table 3). Finally, the subsidence risk for direct influence factors was calculated by Equation (1).
F D = F a · F b · F c

4.2.2. Indirect Influence Factor Group

An interaction matrix is a tool used to represent correlations among multiple components simultaneously. The given components are aligned horizontally in the matrix so that their correlations are expressed in parallel links rather than in one-to-one relations. The approach is based on the concept that every factor of an actual location is related to every other factor, albeit to varying degrees [48].
Hudson [49] positioned subject components as diagonal entries and represented their interactions in the clockwise direction. As shown in Figure 11, the major components A and B are placed as diagonal entries, and the components are linked in the clockwise direction. Accordingly, the components to the right of A and above B represent the sensitivity of A toward B. Conversely, the components to the left of B and below A represent the sensitivity of B toward A. The significance of each component of the matrix can be expressed with partial derivatives. Although interactions can be represented with the matrix in a simple manner, the critical question is how the sensitivity of each component should be determined for best use of the concept for analysis and application of field data [49,50]. Although only two factors were used to explain the concept in Figure 11, the interaction matrix can be constructed for multiple factors based on the identical concept. Although the matrix is mathematically capable of expressing the correlations among all the factors in diagonal entries, during the initial stage, Hudson [49] constructed each component of the matrix, indicating sensitivity by independently considering the linear relationship between linked components (A and B in Figure 11) and intuitively assigning values to the matrix.
Based on the concept that two or more components on the matrix diagonal are related to each other with varying degrees of influence, Jiao and Hudson [51] proposed an algorithm for constructing an interaction matrix for which all of the component values of the matrix are fully coupled. However, each component of the interaction matrix is still determined based on the intuition of experts, making its objective assessment difficult. Yang and Zhang [50] proposed constructing an interaction matrix by using the ANN. To do so, they introduced the concept of the relative strength effect (RES) and an ANN-based calculation algorithm to assess the linkage strength among the components. For this study, a RES value between 0 and 1 was used to represent the mutual relationship between components.
Figure 12 shows the architecture of the ANN specially organized for sensitivity analysis, which consists of the same number of nodes for one input layer, one output layer, and two hidden layers.
Once training is complete, a partial derivative of an input node can be calculated for each output node using the ANN sensitivity analytical system [39]. All five components of the interaction matrix can then be assigned with the calculated partial derivatives. Figure 13 shows an interaction matrix constructed using this process. Every value on the diagonal calculated by the trained ANN is 1.0, which is natural because, in theory, a partial derivative of an identical component should always be 1.0.
According to the basic concept of the applied interaction matrix, the row of the subject factor represents the causes it has on other factors, and the sum of the rows of the matrix indicates its causes on the entire system. Similarly, the column of the factor represents the influence of other factors, corresponding to the effects. The sum of the columns, therefore, represents the effects on the entire system [49,52]. Based on these concepts, a representative value can be assessed to express the degree of influence of each selected factor on the location. The representative value is calculated as follows.
First, the rows and columns of the interaction matrix are added so that each factor represents the cause and effect of the influence on the entire system. The degree of influence of each factor (i) on the entire system as a cause is denoted by Cpi, and as an effect by Epi. In Figure 13, Cpi is indicated on the right of each row heading and Epi is indicated below each column heading. The two categories can be expressed as Equation (2)
C p i = n = 1 n m a x ( I m n ) ,   m = i E p i = m = 1 m m a x ( I m n ) ,   n = i
where Imn represents each component of the m × n interaction matrix, nmax represents the total number of columns, and mmax represents the total number of rows.
Consequently, the interaction matrix must be able to comprehensively express the links among the influence factors of the given location related to subsidence. The influence weights assigned to each factor are expressed as a percentage of the sum of the system’s cause (Cpi) and effect (Epi) according to the following equation. The total sums of causes and effects are expressed as Equation (3)
W i = ( C p i + E p i ) (   C p i +   P p i ) × 100 N g r a d e
where Ngrade is the number of grades assigned to the factor i.
Figure 14 provides the influence weights for the five influence factors calculated from Equation 3. The results show that the most sensitive factor affecting the subsidence for both types of mines was RMR, while the load type and groundwater level were least sensitive for subsidence for coal mines and other types of mines, respectively.
Accordingly, the subsidence risk for indirect influence factors can be defined as
Subsidence   risk   for   indirect   factors   ( F I ) =     W i · V i
where Vi denotes the normalized value for each factor i at the location in question, as indicated in Section 4.1.1. The final possible FI ranges from 0 to 100.

4.2.3. Mine Subsidence Risk Index

The mine subsidence risk (MSR) index is defined as in Equation 5 by combining the two subsidence risks by the direct influence group and the indirect influence group.
MSR index = WD·FD + WI·FI
Here, WD and WI are the weights of subsidence risk for the direct influence factor group and the indirect group, respectively. The two weights can be determined by a series of sensitivity analyses on subsidence for the direct group and the indirect group.
Finally, the MSR index is defined as in Equations (6) and (7), and the final MSR index for all of the subsidence locations is summarized as Figure 15. A higher MSR index implies a higher possibility of subsidence.
MSR index for coal mines = 0.5·FD + 0.5·FI
MSR index for other types of mines = 0.65·FD + 0.35·FI

4.3. Mine Subsidence Grade

Analysis A is used to predict whether subsidence can occur for the given factors, but it always includes the possibility of incorrect prediction as well. Analysis B is necessary to back up the results of analysis A by assessing the MSR index to indicate the extent of risk for subsidence at the location in question.
Finally, in this study, the grade of MSR was classified into the three categories, level 0 (L0, low risk), level 1 (L1, moderate risk), and level 2 (L2, high risk), by combining the results of analyses A and B, as indicated in Figure 16. In this diagram, the MSR index for deciding L0 or L1 is determined considering the average values of the MSR index presented in Figure 15.
The subsidence risk grades determined by the diagram for the training data and the validation data are presented in Figure 17. For the training data of coal mines, only one of the 82 locations where subsidence actually occurred was predicted to be safe, which means that the correct prediction percentage for the training data is 98.8% (=81/82 subsidence locations). For the other types of mines, the correct prediction percentage was 98.5% (=66/67 subsidence locations). For the validation data, the correct prediction percentage was 90% (=9/10 subsidence locations) and 100% (=7/7 subsidence locations) for the coal mines and the other types of mines, respectively.
The overall procedure for deciding subsidence risk grade is summarized in Figure 18.

5. Testing the Model Using New Data

The applicability of the subsidence risk grade was tested with newly collected 22 subsidence data from five abandoned mines. The data consist of 12 locations (6 subsidence and 6 no subsidence locations) at 3 coal mines and 10 locations (7 subsidence and 3 no subsidence locations) at 2 other types of mines. The information of data is summarized in Table 4.
The subsidence risk grade analyses were carried out through the procedure illustrated in Figure 19. The results showed that five of six coal mine locations where subsidence had occurred were predicted to be L2 and that six of the seven locations for other types of mines were predicted as L2. For both types of mines, no location was predicted to be L1. Finally, it was found that the correct prediction percentage was 83% (5/6 subsidence locations) for coal mines and 86% (6/7 subsidence locations) for the other types of mines.
It is thought that the main reason for the lower correct prediction percentage for the new data than the data for training was that the number of newly investigated locations was much smaller than the that for the training data.
Finally, the results imply that the subsidence risk grade model could generally be used to assess the subsidence risk for an abandoned mine. Nevertheless, much more subsidence data with greater variations in ground conditions will be essential and various types of analyses by numerical and empirical approaches, etc. need to be combined to improve the correct prediction percentage of subsidence.
Meanwhile, this analysis process can be used as an initial method screening the specific locations that are at high risk of subsidence, for subsequent detailed analyses by other approaches, for example, the amplitude or rate of the subsidence at high risk can be calculated by empirical methods, methods employing influence functions, and/or methods employing theoretical models [32].

6. Conclusions

This paper establishes a mine subsidence risk index and a subsidence risk grade based on two separate analyses using ANN to predict whether subsidence can occur at a ground location of an abandoned mine. The primary conclusions of the study are summarized below:
(1) From a total of 467 investigated locations in this study, 227 locations from 37 abandoned mines were selected for analysis. Among these, 118 locations were from 15 abandoned coal mines, and 109 locations were from other types of mines (i.e., gold, silver, metal). The numbers of locations where subsidence had and had not occurred were 166 and 61, respectively.
(2) The correlation and multicollinearity analyses among 10 factors related to mine subsidence were performed to improve the prediction accuracy. The results showed that the correlation between subsidence occurrence and UCS, or rainfall was very low and that the multicollinearity between RMR and UCS as well as between groundwater level and rainfall was high. This means that the effects of UCS and rainfall on subsidence are comparatively low, and RMR and groundwater level are closely related to UCS and rainfall respectively. Finally, eight factors excluding UCS and rainfall (i.e., elapsed time, depth to goaf, goaf height, distance to the nearest goaf, RMR (rock mass rating), groundwater level, load types on the investigated location, and the existence of an outcrop on the investigated location) were selected.
(3) The analysis consisted of two separate stages, analysis A and analysis B. Analysis A was conducted to predict whether subsidence could occur using an ANN model of machine learning. Analysis B involved assessing the mine subsidence risk index to indicate the extent of risk for subsidence at the location. After training the model in analysis A, the results showed that the correct prediction percentages for the validation datasets of coal mines and other types of mines were 83% and 75%, respectively.
In analysis B, the eight factors were divided into two groups, a direct influence factor group (depth to a goaf, goaf height, and horizontal distance to the nearest goaf) and an indirect influence factor group (elapsed time, RMR, load type, outcrop, and groundwater level), to assess the mine subsidence risk (MSR) index.
The final determination of subsidence possibility is made using a subsidence risk grade that is established based on the two analyses. The analysis results showed that the correct prediction percentage was 90% (=9/10 subsidence locations) for the 10 validation data of coal mines and 100% (=7/7 subsidence locations) for the 7 validation data of other types of mines.
(4) To verify the prediction model, new datasets for 22 locations—including 9 locations with no subsidence—were selected from 5 abandoned mines and investigated. The correct prediction percentage for coal mines was 83% (=5/6 subsidence locations), and the correct prediction percentage for other types of mines was 86% (=6/7 subsidence locations).
Although this model performed reasonably well in the assessment of potential subsidence at various locations, the sample size was relatively small; thus, additional analyses should be performed at other mine locations and in a larger variety of ground conditions to more clearly verify the model.
(5) It is thought that more reliable results can be obtained when various types of analyses by numerical, statistical, and theoretical approaches are combined. This study can be used as an initial method for screening the specific locations that are at high risk of subsidence and thus require immediate additional investigation, from an extensive mine area for subsequent detailed analyses by the other approaches.

Author Contributions

Conceptualization, Y.K.; Validation, S.S.L.; Formal analysis, Y.K.; Data curation, Y.K.; Writing—original draft preparation, Y.K.; Writing—review, S.S.L.; Supervision, S.S.L.; Project administration, S.S.L.; Funding acquisition, Y.K. All authors have read and agreed to the published version of the manuscript.

Funding

This research was funded by the Mine Reclamation Corporation (MIRECO) in Korea.

Acknowledgments

This study was performed as a part of a research development project and supported by the research fund of Hanyang University (HY-201700000002431). The authors greatly appreciate the financial support.

Conflicts of Interest

The authors declare no conflict of interest. The funders had no role in the design of the study; the collection, analyses, or interpretation of data; the writing of the manuscript; or the decision to publish the results.

References

  1. Li, H.; Zha, J.; Guo, G. A new dynamic prediction method for surface subsidence based on numerical model parameter sensitivity. J. Clean. Prod. 2019, 233, 1418–1424. [Google Scholar] [CrossRef]
  2. Manekar, G.G.; Shome, D.; Chaudhari, M.P. Prediction of subsidence parameters & 3-D analysis at Balaghat Underground Manganese Mine of MOIL Limited, India. Procedia Eng. 2017, 191, 1075–1086. [Google Scholar]
  3. Salmi, E.F.; Karaku, M.; Nazem, M. Assessing the effects of rock mass gradual deterioration on the long-term stability of abandoned mine workings. Tunn. Undergr. Space Technol. 2019, 88, 169–185. [Google Scholar] [CrossRef]
  4. Tzampoglou, P.; Loupasakis, C. Numerical simulation of the factors causing land subsidence due to overexploitation of the aquifer in the Amyntaio open coal mine, Greece. HydroResearch 2019, 1, 8–24. [Google Scholar] [CrossRef]
  5. Zingano, A.; Weiss, A. Subsidence over room and pillar retreat mining in a low coal seam. Int. J. Min. Sci. Technol. 2019, 29, 51–57. [Google Scholar] [CrossRef]
  6. Piggott, R.J.; Eynon, P. Ground movements arising from the presence of shallow abandoned mine workings. In Proceedings of the Conference on Large Ground Movements and Structures, UWIST; Geddes, J.D., Ed.; Pentech Press: London, UK, 1977; pp. 749–780. [Google Scholar]
  7. Brady, B.; Brown, E. Rock Mechanics for Underground Mining; Chapman & Hall: London, UK, 1993. [Google Scholar]
  8. Yue, J.; Cheng, W.; Fan, L. The Study on mathematical model of urban land subsidence based on statistical analysis. In Proceedings of the International Conference on Management and Service Science, Wuhan, China, 20–22 September 2009. [Google Scholar]
  9. Yu, G.; Mi, W.; Wang, D.; Gao, L.; Lu, S.; Li, G. Research on the relationship between the surface dynamic subsidence and overburden separated strata of coal mine and its model. Procedia Eng. 2017, 191, 196–205. [Google Scholar] [CrossRef]
  10. Wang, B.; Xu, J.; Xuan, D. Time function model of dynamic surface subsidence assessment of grout-injected overburden of a coal mine. Int. J. Rock Mech. Min. Sci. 2018, 104, 1–8. [Google Scholar] [CrossRef]
  11. Hejmanowski, R.; Malinowska, A.A. Evaluation of reliability of subsidence prediction based on spatial statistical analysis. Int. J. Rock Mech. Min. Sci. 2009, 46, 432–438. [Google Scholar] [CrossRef]
  12. Goel, S.C.; Page, C.H. An empirical method for predicting the probability of chimney cave occurence over a mining area. Int. J. Rock Mech. Min. Sci. 1982, 19, 325–337. [Google Scholar] [CrossRef]
  13. Cha, Y.; Choi, W.; Buyukozturk, O. Deep learning-based crack damage detection using convolutional neural networks. Comput. Aided Civ. Infrastruct. Eng. 2017, 32, 361–378. [Google Scholar] [CrossRef]
  14. Waziri, B.S.; Bala, K.; Bustani, S.A. Artificial neural networks in construction engineering and management. Int. J. Arch. Eng. Constr. 2017, 6, 50–60. [Google Scholar] [CrossRef]
  15. Juszczyk, M.; Le’sniak, A. Modelling construction site cost index based on neural network ensembles. Symmetry 2019, 11, 411. [Google Scholar] [CrossRef] [Green Version]
  16. Singh, V.; Bano, S.; Yadav, A.K.; Ahmad, S. Feasibility of artificial neural network in civil engineering. Int. J. Trend Sci. Res. Dev. 2019, 3, 724–728. [Google Scholar] [CrossRef]
  17. Garrett, J. Where and why artificial neural networks are applicable in civil engineering. Comput. Aided Civ. Infrastruct. Eng. 1994, 8, 129–130. [Google Scholar] [CrossRef]
  18. Yang, Y.; Zhang, Q. A hierarchical analysis for rock engineering using artificial neural networks. Rock Mech. Rock Eng. 1997, 30, 207–222. [Google Scholar] [CrossRef]
  19. Kapageridis, K. Artificial neural network technology in mining and environmental applications. Mine Plan. Equip. Sel. 2002, 172–179. [Google Scholar]
  20. Rafie, M.; Samimi Namin, F. Prediction of subsidence risk by FMEA using artificial neural network and fuzzy inference system. Int. J. Min. Sci. Technol. 2015, 25, 655–663. [Google Scholar] [CrossRef]
  21. Pishro, M.; Khosravi, S.; Tehrani, S.M.; Mousavi, S.R. Modeling and zoning of land subsidence in the southwest of Tehran using artificial neural networks. Int. J. Hum. Cap. Urban Manag. 2016, 1, 159–168. [Google Scholar]
  22. Pham, B.T.; Tien Bui, D.; Prakash, I.; Dholakia, M.B. Hybrid integration of multilayer perceptron neural networks and machine learning ensembles for landslide susceptibility assessment at Himalayan area (India) using GIS. Catena 2017, 149, 52–63. [Google Scholar] [CrossRef]
  23. Tien Bui, D.; Shahabi, H.; Shirzadi, A.; Chapi, K.; Pradhan, B.; Chen, W.; Khosravi, K.; Panahi, M.; Bin Ahmad, B.; Saro, L. Land subsidence susceptibility mapping in south Korea using machine learning algorithms. Sensors 2018, 18, 2464. [Google Scholar] [CrossRef] [Green Version]
  24. Oh, H.; Syifa, M.; Lee, C.; Lee, S. Land subsidence susceptibility mapping using bayesian, functional, and Meta-Ensemble Machine Learning Models. Appl. Sci. 2019, 9, 1248. [Google Scholar] [CrossRef] [Green Version]
  25. Kim, K.; Lee, S.; Oh, H. Prediction of ground subsidence in Samcheok City, Korea using artificial neural networks and GIS. Environ. Geol. 2009, 58, 61–70. [Google Scholar] [CrossRef]
  26. Zhao, K.; Chen, S. Study on artificial neural network method for ground subsidence prediction of metal mine. Procedia Earth Planet. Sci. 2011, 2, 177–182. [Google Scholar] [CrossRef] [Green Version]
  27. Suh, J.; Choi, Y.; Park, H.; Kwon, H.; Yoon, S.; Go, W. Application of frequency ratio and analytic hierarchy process to subsidence hazard assessment around abandoned coal mine. J. Korean Soc. Geosystem Eng. 2010, 47, 690–704. [Google Scholar]
  28. MIRECO (Mine Reclamation Corporation). Development of a Korean Software and Safety Assessment Technique of Subsidence in Mine Area; MIRECO: Wonju, Korea, 2008. [Google Scholar]
  29. MIRECO (Mine Reclamation Corporation). Development of an Assessment and Prediction Model of Land Subsidence; MIRECO: Wonju, Korea, 2010. [Google Scholar]
  30. MIRECO (Mine Reclamation Corporation). Development of Integrated Hazard Assessment System of Surface Subsidence in Abandoned Mine Area; MIRECO: Wonju, Korea, 2011. [Google Scholar]
  31. MIRECO (Mine Reclamation Corporation). Development of a Subsidence Hazard Assessment Technique in Mine Area; MIRECO: Wonju, Korea, 2014. [Google Scholar]
  32. Kratzsch, H. Mining Subsidence Engineering; Springer: Berlin/Heidelberg, Germany, 1983. [Google Scholar]
  33. Li, W.; Wen, L.; Liu, X. Ground movements caused by deep underground mining in Guan-Zhuang iron mine, Luzhong, China. Int. J. Appl. Earth Obs. Geoinf. 2010, 12, 175–182. [Google Scholar] [CrossRef]
  34. Raja, S.K.; Srinivasa, K.V.; Phanindra, K.B. Training of artificial neural networks in data mining. Int. J. Innov. Technol. Explor. Eng. IJITEE 2013, 3, 214–217. [Google Scholar]
  35. Goldberg, D.E. Genetic Algorithms in Search, Optimization, and Machine Learning; Addison-Wesley: Reading, MA, USA, 1989; p. 412. [Google Scholar]
  36. Lee, S.; Ryu, J.; Won, J.; Park, H. Determination and application of the weights for landslide susceptibility mapping using an artificial neural network. Eng. Geol. 2004, 71, 289–302. [Google Scholar] [CrossRef]
  37. Tunusluoglu, M.; Gokceoglu, C.; Sonmez, H.; Nefeslioglu, H. An artificial neural network application to produce debris source areas of Barla, Besparmak, and Kapi Mountains (NW Taurids, Turkey). Nat. Hazard Earth Syst. Sci. 2007, 7, 557–570. [Google Scholar] [CrossRef] [Green Version]
  38. Nefeslioglu, H.; Gokceoglu, C.; Sonmez, H. An assessment on the use of logistic regression and artificial neural networks with different sampling strategies for the preparation of landslide susceptibility maps. Eng. Geol. 2008, 97, 171–191. [Google Scholar] [CrossRef]
  39. Shin, H. Neural Network Based Constitutive Models for Finite Element Analysis. Ph.D. Thesis, University of Wales Swansea, Swansea, UK, 2001. [Google Scholar]
  40. Coal Industry Promotion Board (CIPB). A Study on the Mechanism of Subsidence over Abandoned Mine Area and the Construction Method of Subsidence Prevention; Coal Industry Promotion Board: Seoul, Korea, 1997; Volume 97, pp. 1–67. [Google Scholar]
  41. Waltham, A.C. Ground Subsidence; Blackie & Son Ltd.: New York, NY, USA, 1989; pp. 49–97. [Google Scholar]
  42. Matignon, R. Neural Network Modeling Using SAS Enterprise Miner; AuthorHouse: Boomington, IN, USA, 2005; p. 150. [Google Scholar]
  43. Jung, Y.; Song, W.; Kang, S. Development of subsidence hazard estimation method based on the depth of gangway. Tunn. Undergr. Space 2008, 18, 272–279. [Google Scholar]
  44. Sahua, P.; Lokhande, R.D. An Investigation of sinkhole subsidence and its preventive measures in underground coal mining. Procedia Earth Planet. Sci. 2015, 11, 63–75. [Google Scholar] [CrossRef] [Green Version]
  45. Sinclair Knight Merz Pty Ltd. Subsidence from Coal Mining Activities; Department of the Environment: Canberra, Australia, 2014.
  46. Ahn, S.; Kim, K. Study on the principle factors related to ground subsidence at abandoned underground coal mine area using probability and sensitivity analysis. In Proceedings of the KSRS Conference, Daejeon, Korea, 29–30 March 2007; pp. 296–300. [Google Scholar]
  47. Choi, J.; Kim, K.; Lee, S.; Kim, I.; Won, J. Prediction of ground subsidence hazard area using GIS and probability model near abandoned underground coal mine. Econ. Environ. Geol. 2007, 40, 295–306. [Google Scholar]
  48. Shin, H.; Kwon, Y.; Jung, Y.; Bae, G.; Kim, Y. Methodology for quantitative hazard assessment for tunnel collapses based on case histories in Korea. Int. J. Rock Mech. Min. Sci. 2009, 46, 1072–1087. [Google Scholar] [CrossRef]
  49. Hudson, J.A. Rock Engineering Systems: Theory and Practice; Ellise Horwood: Chichester, UK, 1992. [Google Scholar]
  50. Yang, Y.; Zhang, Q. The application of neural networks to rock engineering systems (RES). Int. J. Rock Mech. Min. Sci. 1998, 35, 727–745. [Google Scholar] [CrossRef]
  51. Jiao, Y.; Hudson, J.A. The fully-coupled model for rock engineering systems. Int. J. Rock Mech. Min. Sci. Geomech. Abstr. 1995, 32, 491–512. [Google Scholar] [CrossRef]
  52. Cancelli, A.; Crosta, G. Hazard and Risk Assessment in Rockfall Prone Areas. In Proceedings of Risk and Reliability in Ground Engineering; Institute of Civil Engineering: London, UK, 1993; pp. 177–190. [Google Scholar]
Figure 1. Typical examples of subsidence in this study.
Figure 1. Typical examples of subsidence in this study.
Applsci 10 01302 g001
Figure 2. Study process.
Figure 2. Study process.
Applsci 10 01302 g002
Figure 3. Artificial neural network structure.
Figure 3. Artificial neural network structure.
Applsci 10 01302 g003
Figure 4. An example of GDAP training.
Figure 4. An example of GDAP training.
Applsci 10 01302 g004
Figure 5. Locations of selected and investigated abandoned mines.
Figure 5. Locations of selected and investigated abandoned mines.
Applsci 10 01302 g005
Figure 6. Correlation between influence factors and subsidence.
Figure 6. Correlation between influence factors and subsidence.
Applsci 10 01302 g006
Figure 7. Correlation matrix for checking multicollinearity.
Figure 7. Correlation matrix for checking multicollinearity.
Applsci 10 01302 g007
Figure 8. Selected influence factors.
Figure 8. Selected influence factors.
Applsci 10 01302 g008
Figure 9. Ranges of investigated data for each factor.
Figure 9. Ranges of investigated data for each factor.
Applsci 10 01302 g009
Figure 10. Analyzed optimum structure.
Figure 10. Analyzed optimum structure.
Applsci 10 01302 g010
Figure 11. Principle of interaction matrix (modified after [49]).
Figure 11. Principle of interaction matrix (modified after [49]).
Applsci 10 01302 g011
Figure 12. ANN structure for the analysis of indirect influence factors.
Figure 12. ANN structure for the analysis of indirect influence factors.
Applsci 10 01302 g012
Figure 13. Interaction matrix for (a) coal mines and (b) other types of mines.
Figure 13. Interaction matrix for (a) coal mines and (b) other types of mines.
Applsci 10 01302 g013
Figure 14. Influence weights for (a) coal mines and (b) other types of mines.
Figure 14. Influence weights for (a) coal mines and (b) other types of mines.
Applsci 10 01302 g014
Figure 15. Range of MSR index: (a) coal mines (b) other types of mines.
Figure 15. Range of MSR index: (a) coal mines (b) other types of mines.
Applsci 10 01302 g015
Figure 16. Subsidence risk grade diagram.
Figure 16. Subsidence risk grade diagram.
Applsci 10 01302 g016
Figure 17. Subsidence risk grade for coal mines and other types of mines: (a) training data and (b) validation data.
Figure 17. Subsidence risk grade for coal mines and other types of mines: (a) training data and (b) validation data.
Applsci 10 01302 g017
Figure 18. Summary of the subsidence risk grade analysis procedure.
Figure 18. Summary of the subsidence risk grade analysis procedure.
Applsci 10 01302 g018
Figure 19. Subsidence risk grade for the new data.
Figure 19. Subsidence risk grade for the new data.
Applsci 10 01302 g019
Table 1. Number of data investigated.
Table 1. Number of data investigated.
Total (37)At Coal Mines (15)At Other Mines (22)
SumSubsidenceSumSubsidenceSumSubsidence
YesNoYesNoYesNo
Total collected data467388793002574316713136
Selected dataSum2271666111892261097435
Training2071495810682241016734
Validation2017312102871
Table 2. Accuracy of analysis A by the data type.
Table 2. Accuracy of analysis A by the data type.
Data TypeAccuracy
(Number of Data)
Coal MinesOther Type Mines
Training data97%97%
(106)(101)
Validation data83%75%
(12)(8)
Table 3. Points for direct influence factors with their ranges.
Table 3. Points for direct influence factors with their ranges.
Depth (m, Fa)PointsHeight (m, Fb)Correction FactorDistance (m, Fc)Correction Factor
CoalOtherCoalOtherCoal and Other
< 1000< 100010< 1000< 10001.00< 10000.1
<=60<=3010<=27<=450.95<=200.1
<=54<=2720<=24<=400.90<=180.2
<=48<=2430<=21<=350.85<=160.3
<=42<=2140<=18<=300.80<=140.4
<=36<=1850<=15<=250.75<=120.5
<=30<=1560<=12<=200.70<=100.6
<=24<=1270<=9<=150.65<=80.7
<=18<=980<=6<=100.60<=60.8
<=12<=690<=3<=50.55<=40.9
<=6<=3100=0<=00.00<=21.0
Table 4. New subsidence data.
Table 4. New subsidence data.
TypeMine
Name
No.SubsidenceElapsed Time 1 (Year)RMRGroundwater Level (m)OutcropLoad Type 2Depth 3 (m)Height 4 (m)Distance 5 (m)
CoalCM11Yes15 25 18 Yes 0 3 16 0
2Yes 15 20 13 Yes1 2 20 0
3Yes 19 16 10 No0 6 13 0
4No19 26 25 No 0 7 13 3
5No19 32 29 No 0 15 8 8
CM21Yes26 22 27 Yes3 14 2 0
2Yes 26 21 31 No 3 30 2 42
3No 26 17 20 No 0 30 2 126
4No 26 20 25 No 0 35 1 126
CM31Yes24 33 20 Yes3 19 13 9
2 No24 26 22 No3 32 16 44
3No24 26 26 No 3 38 11 55
Other typesOM11Yes 39 53 4 No3 3 2 14
2Yes 39 53 2 No3 1 2 5
3No 39 37 4 No 0 3 2 22
4No39 37 5 No0 10 2 32
OM21Yes 31 64 12 No 0 33 5 5
2Yes 31 30 1 No 4 1 16 1
3Yes 31 40 4 No 3 4 11 2
4Yes 31 50 7 No2 7 6 3
5Yes 31 60 10 No1 10 1 4
6No31 42 7 No 0 27 5 7
1 Elapsed time: The time elapsed since the corresponding abandoned mine was closed. 2 Load type: 0 (none), 1 (house), 2 (concrete structure), 3 (railway/road), 4 (bridge). 3 Depth: The depth to goaf crown. 4 Height: The height of the goaf beneath the corresponding subsidence. 5 Distance: The distance from the nearest underground goaf to the ground location to be investigated.

Share and Cite

MDPI and ACS Style

Kim, Y.; Lee, S.S. Application of Artificial Neural Networks in Assessing Mining Subsidence Risk. Appl. Sci. 2020, 10, 1302. https://doi.org/10.3390/app10041302

AMA Style

Kim Y, Lee SS. Application of Artificial Neural Networks in Assessing Mining Subsidence Risk. Applied Sciences. 2020; 10(4):1302. https://doi.org/10.3390/app10041302

Chicago/Turabian Style

Kim, Yangkyun, and Sean S. Lee. 2020. "Application of Artificial Neural Networks in Assessing Mining Subsidence Risk" Applied Sciences 10, no. 4: 1302. https://doi.org/10.3390/app10041302

APA Style

Kim, Y., & Lee, S. S. (2020). Application of Artificial Neural Networks in Assessing Mining Subsidence Risk. Applied Sciences, 10(4), 1302. https://doi.org/10.3390/app10041302

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