Abstract

Damage to mine shafts located below thick loose aquifers, caused by coal seam mining, can seriously affect the operational safety of mines. In view of such a problem, the comprehensive weight method and the fuzzy matter-element analysis method were used to analyze and evaluate the stability of mine shafts below thick, loose aquifers. Based on data relating to shaft structure parameters collected from fifteen damaged shafts in coal mines in East China, a comprehensive analysis was performed. The results showed that (1) the two indexes of surface subsidence velocity and water level drop value at the bottom of a loose aquifer have a large influence on the stability of shafts below thick loose aquifers, and (2) the predicted results of thirteen sets of samples were the same as the results measured during actual production; the accuracy of the model was 86.67%. The comprehensive weight method and the fuzzy matter-element analysis model both show good reliability for evaluating the stability of mine shafts below thick loose aquifers and can provide a scientific reference for the analysis of the stability of mine shafts in areas with thick loose strata, as well as for the design of comprehensive control plans that could subsequently be implemented.

1. Introduction

Indirect mining damage to vertical mine shafts located below thick loose aquifers is a newly identified potential geological disaster facing mines. Underground mines located below such aquifers are widely distributed in Huainan, Huaibei, Xuzhou, Jining, and other regions of China. The geological conditions around coal seams beneath thick loose aquifers are complex. The thickness of thick loose strata exerts great stress on any shaft located below them. At the same time, the bottom stratum of thick loose strata is a layered structure, with the subaquifer and the subaquiclude alternately connected, through which water can easily flow following mining disturbance, posing a hazard to shafts [14]. During the mining process in the mining face, the overburden layer of the coal seam gradually bends and sinks, generating mining fissures. It is easy to run through the bottom loose aquifer at the upper part of the bedrock, causing water and sand to gradually penetrate downward along these mining fissures, resulting in consolidation of the loose aquifer due to water losses [5]. In the long-term, this water loss process leads to consolidation and compression occurring in the aquifer due to it being affected by the load of the upper loose strata, resulting in structural damage to shaft walls [69]. A schematic diagram showing structural damage to a shaft following disturbance by mining is shown in Figure 1.

Following instability deformation of loose sections of a vertical shaft, an affected mine will generally be shut down for maintenance. In severe cases, a large amount of water and sand can burst through the shaft wall and damage the shaft, putting the lives of underground workers in serious danger and causing huge economic losses for the mine. According to statistics, since the 1980s more than twenty mine shafts in Huaibei, Huainan, Xuzhou, Jining, and other regions have suffered various degrees of damage, which have seriously affected the safe production of coal in each mine (see Table 1).

There is a need, therefore, for an in-depth analysis of nonmining failures of mine shafts. It is important to explore the various factors influencing this and their degrees of correlation. The accurate evaluation of the stability of mining shafts under thick loose aquifers will be useful for the prediction and prevention of disasters and accidents caused by nonmining damage in shafts located below thick loose aquifers, and will be of great practical significance to ensure the safety and reliability of mining processes [1216].

To date, most studies of mine shafts in China and other countries have focused on their failure mechanisms and ways to address these failures, which usually involves the use of solid filling technology to control the deformation of shafts [17, 18]. There have been relatively few studies that have evaluated the stability of mining shafts below thick loose aquifers. Liu et al. [19] constructed a system for predicting nonmining fractures in coal mine shafts, based on the concept of an artificial neural network. Kong [20], using theoretical analysis, numerical simulation, and a field measurement method, performed an in-depth analysis of the disaster-causing factors in frozen vertical shafts below deep, thick overburden, and determined the spatiotemporal dynamic processes involved in shaft instability. Huang and Guo [21] established a system for monitoring the stability of vertical shafts by means of observation, through the arrangement of a number of on-site observation points. Liu et al. [22] analyzed a number of factors affecting the stability of shafts and built a vertical shaft stability evaluation model in combination with fuzzy clustering analysis theory. Xu et al. [23] established shaft safety evaluation models based on the minimum distance discrimination method to predict the ranking of the stability of mine shafts after three years, using as examples the main shaft, auxiliary shaft, east air shaft, and west air shaft of Xinglongzhuang Coal Mine. Lyu and Wang [24] established an automatic, multilayer, deformation monitoring system for the complex geological effects resulting from the weakly cemented formation of a shaft wall, to evaluate changes in the stability of vertical shaft walls. These scholars studied and analyzed mine shafts using various theoretical methods and technical means, and their work has greatly improved the prediction and control of disasters and accidents caused by nonmining failures of shafts below thick loose strata. However, in view of the complex geological conditions that influence mining areas under thick loose aquifers, some problems remain around shaft stability evaluation and prediction that urgently need to be solved. With respect to numerical simulations, previously, the simulation process has been somewhat simplified, and it is difficult to represent the true geological changes. The main factors affecting the deformation and failure of shafts and their ranking have yet to be determined theoretically. Existing shaft failure evaluation and prediction studies have not comprised any in-depth, comprehensive consideration of positive and negative factors that can influence shaft stability. The importance of human subjective and objective factors requires further exploration, and the study of methods for evaluating predictive models needs to be improved.

To address some of the problems that exist with the study results described above, this paper compiles various research from the perspective of both subjective and objective analysis, with reference to multiple examples of field data, and combines them with the input of many experts’ years of work experience. These various data were integrated using the combination weight method, and the comprehensive weight value [2527] of each factor causing the deformation of mine shafts below thick loose aquifers was calculated. In addition, the degrees of correlation between the influencing factors were determined using fuzzy matter-element analysis theory. Using this approach enabled the effective evaluation of the stability of mine shafts below thick loose aquifers.

2. Evaluation Methods

2.1. Comprehensive Weight Method
2.1.1. Entropy Weight Method

The entropy weight method [2831] is an objective analysis method for determining the weight of evaluation indexes. The basic principle is based on the degree of variation between evaluation indexes. If the degree of variation is larger and the information entropy is smaller, the weight of the corresponding index is larger. This method is based on actual geological mining parameters and is conducive to producing more accurate results. If the number of evaluation samples is “” and the number of evaluation indexes is “,” the evaluation matrix can be established as follows: (1)Standardization of indexesIn a multiindex evaluation system, it is not conducive to obtaining accurate results, due to the influence of factors such as the units and the dimensions of each index. Therefore, in order to ensure the reliability of the evaluation results, all evaluation results should undergo a process of standardization.Standardization of positive indexes is as follows: Standardization of negative indexes is as follows: (2)Weight assignment of each index in the matrix where refers to the weight of No. “” evaluation index of No. “” evaluation sample in the matrix.(3)Calculate the entropy value of each index where is a constant value, .(4)Solve the entropy weight of each index where , . From the above equation, it can be determined that the larger the entropy value and the smaller the entropy weight, the lower the degree of importance of the evaluation index.

2.1.2. Analytical Hierarchy Process (AHP) Method

The analytical hierarchy process (AHP) is a subjective analysis method, first proposed by Saaty in the early 1970s, to quantify qualitative factors. It is described as a way to convert intangible factors into numerical values, following which a series of pairwise comparisons can be carried out based on the degree of importance of each influencing factor to calculate the weight of each factor. This method solves the problem of evaluating factors that are difficult to quantify or prioritize [3234]; the basic steps are as follows: (1)Establishment of a hierarchical structure system, which can generally be divided into three grades, namely, the target grade, the intermediate classification grade, and the factor layer, as shown in Figure 2(2)Construction of a judgment matrix. When the weight analysis is carried out for each factor of each grade, two factors are selected to compare their relative scales and determine their grades according to their degrees of importance. For details, please refer to Table 2(3)Calculation of the weight of each factor. First, the degrees of importance of each grade’s influencing factors are quantified then normalized. Second, fixed values of the judgment matrix are obtained, followed by an analysis of the proportion contributed by the degrees of importance of each grade’s influencing factors, to obtain each factor’s weight value .(4)Consistency test. AHP requires that the judgment matrix has a general level of consistency, which can ensure the calculated results are basically reasonable. The calculation formula is as follows: where is the consistency ratio, is the maximum eigenvalue of the consistency matrix, is the number of factors compared in pairs, and is the random consistency index; the value of which is determined according to Table 3.

2.1.3. Combined Weighting Method Based on Game Theory

Subjective and objective weight information in relation to the stability of shafts in areas with thick loose strata were comprehensively analyzed using game theory, so that the calculated results of the comprehensive weighting of the evaluation indexes were comprehensive and scientific. A single weight vector can be expressed as , and the formula of the combined weight method is as follows:

Based on the concept of game theory, the coefficient was optimized to minimize the variance between the integrated weight vector and the single weight vector , thereby obtaining an ideal weight vector.

The coefficient obtained was then normalized:

The integrated weight vector is as follows:

2.2. Fuzzy Matter-Element Analysis Method

The concept of Extenics first appeared in a paper published by the Chinese scholar Professor Cai Wen, entitled “Extension set and non-compatible problems”. Its theoretical background is in matter-element theory and extension set theory. Through the qualitative analysis of a matter-element, combined with the quantitative calculation of a correlation function, it transforms various indexes that influence the object under evaluation into a compatible problem; thus, a comprehensive evaluation model [3541] that is more in line with objective reality is obtained. The calculation process of the fuzzy matter-element analysis method is as follows: (1)Determination of the classical fieldA matter-element analysis model generally consists of three parts: matter, character, and value. If the number of evaluation indexes is “,” it is an -dimensional matter-element model. where is the classical field matter element, is the th stability evaluation grade of the shaft, is the th index of the evaluation index system, is the value range of the th evaluation index in the th grade, and and , respectively, represent the lower and upper limits of the value.(2)Determination of the segment field where is the classical field matter element, is all grades of the shaft stability evaluation, is an extended quantization range of index , and the variation range is from to .(3)Determination of the matters to be evaluatedThe stability evaluation grade of the mine shaft to be evaluated by the matter element is as follows: where is the object element to be evaluated, is the th eigenvalue of the th grade, and is the actual value of the stability evaluation index of the mine shaft to be evaluated.(4)Calculation of the correlation functionWhen constructing the correlation function, the concepts of the extension distance and position are introduced to describe the positional relationship between the point and the interval, so that it can express the degrees of certain attributes of the object. The expression is as follows: where is the distance from to the classical field interval , is the distance from to the segment field interval , and is the correlation index function.(5)Calculation of the comprehensive correlation degree where is the weight coefficient of the shaft stability evaluation index and is the comprehensive correlation degree of the object to be tested at the th grade.(6)Determination of the category and grade variable eigenvalues of the object to be evaluated: if , then is evaluated to belong to the category denoted as follows: Then, the grade variable eigenvalue of is as follows:

3. Model Preparation

3.1. Evaluation Index System and Determination of the Evaluation Grade

Since mine shafts located below thick loose aquifers are subject to the influence of complex geological mining factors, which may lead to instability and failure, it is undoubtedly very important to establish a standard evaluation index system when studying the degree of stability of shafts. Based on previous study results and field monitoring data of shaft deformations from recent years, it was concluded following comprehensive analysis that the factors affecting shaft instability can be divided into geological factors (), shaft factors (), and human factors (). Geological factors can be further divided into surface subsidence velocity (1), cumulative surface subsidence (2), loose strata thickness (3), and water level drop value at the bottom of a loose aquifer (4). Shaft factors can be further divided into shaft wall thickness (1), construction method and type of shaft wall (2), and the inner-outer diameter ratio of the shaft (3). Human factors include the quality of the shaft wall construction (1). (1)Surface subsidence velocity: the surface subsidence velocity reflects the magnitude of the subsidence of the stratum per unit time, which can directly affect the stability of a shaft structure. Under normal circumstances, the greater the surface subsidence velocity, the greater the magnitude of subsidence of the formation, and the greater the stress on the shaft; this is more likely to damage the shaft structure or cause the shaft to be deflected, increasing the likelihood of shaft instability [42, 43](2)Cumulative surface subsidence: this can reflect the historical total amount influenced by underground mining and the cumulative ground subsidence (including the amount of cumulative compression of the loose strata). If other factors are kept unchanged, the greater the cumulative surface subsidence, the greater the additional stress applied to a shaft, which is likely to damage the fixed structure of the shaft outer wall, resulting in lower stability of the shaft(3)Loose strata thickness: loose strata is a type of unsaturated soil composed of three phases, solid, liquid, and gas. Its mechanical properties are extremely complex, and the skeletal structure composed of solids plays a supporting role [44]. Generally speaking, if the loose strata is thicker, the additional stress action (compression) it applies to the shaft is more obvious, and the more likely damage will occur in the shaft at the bottom of the loose section, thus reducing the stability of the shaft(4)The water level drop value at the bottom of a loose aquifer: the magnitude of any drop in water level at the bottom of a loose aquifer plays a very important role in shaft stability. The change in water content of loose soil has a great influence on that soil’s mechanical properties. With a decrease in water level, the sand in an aquifer becomes compacted and consolidated, leading to additional stress being applied to a shaft; this can cause serious damage to the outer wall structure of a shaft at the aquifer, such as generating cracks on the shaft wall, and sudden water burst accidents are more likely to occur in this case [45, 46](5)Shaft wall thickness: the shaft wall thickness can effectively bear the asymmetric additional stress generated by the vertical displacement of loose strata, which can protect the shaft. In general, as shaft thickness increases, the resistance of the shaft structure becomes stronger and the stability of the shaft is increased [47](6)Construction method and shaft wall types: methods of shaft construction include the freezing method and the drilling method, the latter being more conducive to maintaining shaft stability. Shaft wall types include single-layer shaft walls, double-layer shaft walls, prefabricated shaft walls, and composite shaft walls. The drilling method is used for prefabricated shaft walls, while the freezing method can be adopted for all four shaft wall types. In terms of waterproofing and the ability to resist deformation, the composite shaft wall is the best, followed by the prefabricated shaft wall, then the double-layer shaft wall. In this paper, a comprehensive qualitative analysis of the factors influencing shaft stability was carried out, as shown in Table 4Among them, it is 0.7 for the single-layer shaft wall and the freezing method, 0.5 for the double-layer shaft wall and the freezing method, 0.3 for the prefabricated shaft wall and the drilling method, and 0.1 for the composite shaft wall and the freezing method.(7)Inner-outer diameter ratio of the shaft: the larger the outer diameter of a shaft, the greater the contact area with the stratum, and the greater the additional stress applied. The smaller the inner diameter of a shaft, the greater its ability to resist deformation, and the more stable it is. So, it is more convincing to use the inner-outer diameter ratio of the shaft as an influencing factor. If the inner-outer diameter ratio of a shaft is larger, the shaft’s stability is lower, as is its ability to resist deformation and its waterproof effectiveness(8)Shaft wall construction quality: the construction quality of a shaft wall is an important factor affecting the quality of a shaft and represents a human factor. In order to highlight the importance of construction quality on shaft stability, a qualitative analysis was performed, combined with the situations of previous field examples, as shown in Table 5

There have been few studies on the degree of deformation and damage to shafts, so until now, there has been no clear qualitative division between the levels of damage. Based on many years’ monitoring data of the postmining deformation of shafts below thick loose strata, together with the results of relevant studies, we categorized the stability of shafts below thick loose aquifers into five grades, following a comprehensive consideration of variations in a number of factors. These grades are outlined as follows. Grade I refers to new shafts that have just been put into production and have very low deformation values; the stability of such shafts is defined as being safe. Grade II shafts are defined as having small deformations due to influencing factors and undamaged shaft wall surfaces; the stability of such shafts is also defined as being safe. Grade III shafts are those shafts in which deflection or instability has occurred but which have been repaired in a timely manner; they belong to the critical safety grade. Grade IV includes shafts in which instability has occurred and the shaft wall has been affected to a certain extent, but the fracture height and depth are less than 1 m and 100 mm, respectively; these shafts belong to the dangerous grade. Grade V shafts are shafts in which instability has occurred, large pieces of concrete have fallen from the shaft walls, and the depth of the damage is more than 100 mm; these shafts are classified as extremely dangerous.

At the same time, with reference to relevant scholars’ analyses of shaft deformation and shaft stability evaluation, and geological parameter monitoring results from occasions when shaft instability has occurred, we performed a standard range classification of the parameter values according to the abovementioned shaft stability grade classification (Table 6). This is more conducive to the evaluation of shaft failure in areas with thick loose strata, and provides safe and reliable guidance for the prediction and prevention of shaft failure when mining coal seams below thick loose strata.

3.2. Model Research Process

Figure 3 displays the present research process. First, the corresponding weight vectors of various evaluation indexes were calculated using the entropy weight method and AHP. Next, coefficient optimization was conducted on the weight vectors by means of the combination weight method and comprehensive weight vectors under both subjective and objective conditions were obtained. Finally, the comprehensive weight vectors were substituted into fuzzy matter-element model to calculate the evaluation grades of different shafts.

3.3. Selection of the Example Data

Mines where shaft rupture accidents have occurred in East China are generally covered by thick loose strata or very thick loose strata. In view of the large differences between the mechanisms and phenomena of mine shaft failures under thick loose strata compared with those under very thick loose strata, this paper focuses on analyzing the influence of thick loose strata on the stability of mine shafts. This will help to ensure that the evaluation model we have established is more representative and appropriate (the influence of very thick loose strata on the stability of mine shafts below them will be investigated in future work). Taking the four mining regions of Huaibei, Huainan, Xuzhou, and Zhangzhou as examples, the parameters of factors influencing typical mine shafts were selected for comprehensive analysis (please see Table 7 for details).

4. Model Calculation Results

4.1. The Process of Calculation Using the Comprehensive Weight Method
4.1.1. Evaluation Indexes Affecting Shaft Stability when Calculated by the Entropy Weight Method

Many different shaft instability phenomena occur in mines in East China; therefore, we collected geological data with distinctive characteristics from each mining area (Table 7) and calculated the factors influencing shafts using the entropy weight method to determine the major factors that affected shaft stability.

First, the indicated value of each influencing factor was determined. The indicated value of the positive index was 1 and that of the negative index was 0. The indicated value of each influencing factor is shown in Table 8.

The parameters of each index were standardized with reference to the values indicated above. Equations (8), (9), (10), (11) and (12) were then coded using MATLAB software. The original geological parameters collected were substituted, and the indicated values were used to calculate the results. Finally, entropy values and entropy weights for each index were obtained, as shown in Table 9.

It can be seen from Table 9 that the water level drop value at the bottom of a loose aquifer has the largest weight value, followed by surface subsidence velocity, loose strata thickness, shaft wall construction quality, shaft wall thickness, construction method, and shaft type, with the cumulative surface subsidence and the inner-outer diameter ratio of the shaft weighted least. The total weighting of the first six influencing factors is 0.8364, which means they play a decisive role in the stability of shafts below thick loose aquifers. Therefore, an in-depth study of these leading factors that influence shaft stability will help prevent and control shaft failure when mining coal seams located below thick loose aquifers.

4.1.2. Evaluation Indexes Affecting Shaft Stability when Calculated by the AHP Method

The AHP method was used to qualitatively analyze evaluation indexes affecting shaft stability from a subjective perspective. By combining these findings with the empirical analyses of factors affecting the shaft stability reported by many experts and scholars, the degree of importance of each index was calculated and ranked.

First, based on the opinions of various experts, the results of their analyses of each influencing factor were summarized, and a judgment matrix was constructed according to the grades (please see Tables 1013 for details).

4.1.3. Evaluation Indexes Affecting Shaft Stability when Calculated by the Combined Weight Method

Based on the concept and theoretical method of game theory, the results calculated using the subjective weight method and the objective weight method were substituted into Equations (9), (10), (11) and (12), and the weight coefficients and were determined. Finally, the comprehensive weight value, , was calculated (please see Figure 3 for the results).

It can be seen from Figure 4 that the weights of the influencing factors were 0.1826 for the water level drop at the bottom of a loose aquifer, 0.1576 for the surface subsidence velocity, 0.1273 for the loose strata thickness, 0.1249 for the shaft wall construction quality, 0.1185 for the shaft wall thickness, and 0.1181 for the construction method and shaft wall type. The total weighting of these six influencing factors was 0.829, so they can be used to effectively analyze and evaluate the stability of shafts below thick loose aquifers. The comprehensive weight method can reasonably summarize the results of the statistical analysis carried out using the entropy weight method and the AHP method. This means it can more effectively interpret the actual situation of shaft stability under both subjective and objective conditions, and help to reliably rank the degrees of importance of each influencing factor. Thus, the leading factors affecting shaft deformation were explored, providing a reference for reducing shaft deformation and controlling it within a safe range.

4.2. Fuzzy Matter-Element Analysis Process Calculation

On the basis of the established evaluation index system and evaluation grades, and after the weight of each influencing factor was determined by the comprehensive weight method, further comprehensive evaluation and analysis of the stability of shafts below thick loose aquifers was carried out. The specific calculation procedures were as follows: (1)According to the evaluation index standard range for the stability of shafts below thick loose strata (Table 6), the classical field and the segment field were constructed as shown in Table 14(2)Calculation of the degree of correlation. In this paper, we chose the auxiliary shaft of Linhuan Coal Mine as an example. With reference to Table 14 and combined with Equations (16) and (17), five grades of correlation degree of the eight evaluation indexes were calculated. Then, the actual category of each index was determined, as shown in Table 15(3)Calculation of the comprehensive degree of correlation and determination of the stability grade that each shaft belongs to the following:

From the results of the calculation shown in Table 16, it can be seen that by using a comprehensive analysis based on the comprehensive weight method and the matter-element analysis model, the actual situations in fifteen mine shafts could be compared, and the evaluated safety grades of thirteen of these shafts were correctly judged, giving a judgment accuracy of 86.67%. The comprehensive analysis results of 15 mine shafts using the entropy weight method and fuzzy matter-element analysis model were compared with the actual conditions. In terms of the shaft safety and stability evaluation grade, 12 shafts were accurately judged, with a judgement accuracy of 80%. The comprehensive analysis results of 15 mine shafts using the AHP and the fuzzy matter-element analysis model were also compared with the actual results. The evaluation grades based on the shaft safety and stability of 12 shafts were accurately identified, with a judgment accuracy of 80%. The combination of the comprehensive weighting method and the fuzzy matter-element analysis model was superior to the other two evaluation models in terms of identification accuracy, i.e., it can truly reflect the stability of a vertical shaft below a thick loose aquifer.

For shafts where a Class IV or V failure occurred, the surface subsidence velocity and the water level drop value at the bottom of a loose aquifer were larger than those of other shafts: the surface subsidence velocity was generally greater than 36 mm/a, and the water level drop value at the bottom of a loose aquifer was greater than 45 m. Some shafts were heavily influenced by the shaft wall construction quality, which led to misjudgment.

Above all, the shaft safety and stability evaluation model established by using the comprehensive weight method and the fuzzy matter-element analysis method showed good reliability, and can be used to predict the deformation of mine shafts in regions of East China where there are overlying thick loose strata. This will be beneficial for the judgment of stability and failure prevention in mine shafts in these regions.

The misjudgment of the auxiliary shaft of Linhuan Coal Mine and the main shaft of Dongtan Coal Mine in this analysis may be due to the following reasons: (1) since there are relatively few statistical cases of mine shaft accidents, the matter-element analysis model established still needs to be improved; (2) the evaluation index standard range (Table 6), established with reference to the results of previous studies and field experience, is not sufficiently accurate and needs further refinement.

5. Conclusions and Prospects

5.1. Main Conclusions

In this paper, the comprehensive weight method and the fuzzy matter-element analysis method were used to analyze the main factors influencing the stability of shafts below thick loose aquifers, and an analysis model for determining a shaft stability evaluation grade was established. This method provided a novel way to quantitatively evaluate the stability of shafts below the thick loose strata of the Huaihe River and Jining regions, and has important practical significance and safety benefits. The main conclusions are as follows: (1)The combined weight method based on game theory was used to rationally distribute the weights of the evaluation indexes obtained using the entropy weight method and the AHP method. This avoided the objective influence of a single geological condition and reduced the influence of subjective factors, thus enabling a more scientific and accurate conclusion to be drawn(2)The ranking of various influencing factors obtained by the combined weight method was as follows: water level drop value at the bottom of a loose aquifer, surface subsidence velocity, loose strata thickness, shaft wall construction quality, shaft wall thickness, construction method and shaft wall type, cumulative surface subsidence, and the inner-outer diameter ratio of a shaft(3)The matter-element analysis model constructed in this paper showed good reliability and can provide scientific reference values for the prediction and analysis of the stability of mine shafts below thick loose strata in the Huaihe River and Jining regions

5.2. Prospects

(1)The matter-element analysis model was based on existing sample data, so the shaft stability evaluation index standard may not reflect the same factors that influence single large deformations of shafts in some mining areas, which may have an influence on the final evaluation results(2)The thick loose strata in the Huaihe River and Jining regions are composed of Quaternary strata, whereas the thick loose strata in Juye Coal Mine in the Heze region is composed of Quaternary and Neogene strata and has a thickness of about 550 to 850 m. Shaft failures in different regions have their own particular characteristics. In the future, studies into such failures in different regions should be carried out to expand the compatibility and accuracy of the matter-element analysis model described in this paper(3)In view of the complexity of evaluating the stability of shafts below thick loose aquifers, more example field data should be collected to further study the actions and mechanisms of the factors influencing the stability of shafts below thick loose aquifers, thus improving the accuracy of the model’s predictions and judgments

Data Availability

All data, models, or codes generated or used during the study are available from the corresponding author by request.

Conflicts of Interest

The authors declare no conflicts of interest.

Acknowledgments

The authors of the paper would like to extend grateful thanks to the National Natural Science Foundation of China (51774199) for its support, the Haizi Coal Mine of Huaibei Mining Group, the Linhuan Coal Mine, Zhuxianzhuang Coal Mine, and Zhangshuanglou Coal Mine of Xuzhou Mining Group, the Pansan Coal Mine and Panyi Coal Mine of Huainan Mining Group, and the Xinglongzhuang Coal Mine, Nantun Coal Mine, and Dongtan Coal Mine of Yankuang Group Co., Ltd. for providing relevant shaft failure data, and also my senior university fellow Zhu Xianxiang for his suggestions on how to improve this paper.