Advances in Civil Engineering
Volume 2021 (2021), Article ID 6634331, 8 pages
https://doi.org/10.1155/2021/6634331

Research Article

Research on Key Factors Influencing Surface Subsidence of Paste Backfilling Mining in Thick Coal Seam of Deep Mine

Meng Zhang,¹ Hui He,¹ Xin Jin,¹ Yonglong Qu,¹ and Hongjun Guo²

¹School of Civil & Architecture Engineering, Xi’an Technological University, Xi’an 710021, China
²Jiangsu Vocational Institute of Architectural Technology, Xuzhou 221116, China

Correspondence should be addressed to Meng Zhang; zhangmeng@xatu.edu.cn

Received 25 December 2020; Accepted 26 May 2021; Published 7 June 2021

Academic Editor: Hualei Zhang

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

Abstract.

The research on the key factors of surface subsidence in paste backfilling mining of thick coal seam in deep mine is a complex system engineering, which involves mining, backfilling, support, subsidence, safety, and other aspects. At present, there is no systematic research on the key factors of surface subsidence in paste backfilling mining of thick coal seam more than 6 meters in a deep mine. In this paper, field research, laboratory experiments, theoretical analysis, and other research methods are used to carry out the research about 3# coal seam under buildings in the Lu’an area. The main conclusions are as follows: through the construction of the fuzzy extension model of surface subsidence in paste backfilling mining, five types of surface subsidence are obtained, including overburden structure, roof subsidence before backfilling, nonbackfilling account, the strength of filling body, and backfilling technology. It is the key factor to control the surface subsidence; the optimization measures are given to provide a reference for the reasonable design of paste backfilling mining in other working faces. The research conclusion has a certain reference value for solving a series of problems such as coal resources under buildings, solid waste disposal, controlling strata and surface movement, improving the recovery rate of mine resources, and extending the service life of mine.

1. Introduction

Paste backfilling mining in coal mines can effectively solve the serious waste of resources, low recovery rate, and ecological and environmental damage caused by high-intensity mining, which is an important technical means to realize the sustainable development of coal mine production [1–4]. In 1979, paste backfilling was first carried out in Germany for the purpose of waste treatment. Due to fewer resources, the technology has not been further developed [5–8]. In 2006, it was the first industrial test of paste backfilling mining in China. At present, the technology has been applied to nearly 20 pairs of mines in 6 mining areas in China, realizing green mining. Practice has proved that the effect of paste backfilling is good, but the technology still has some prominent problems, such as low production of filling face, not suitable for deep mine thick seam paste backfilling mining, high filling cost, and so on, which seriously restricts the wide application of the technology. Therefore, it is urgent to study the paste backfilling technology of thick coal seam in a deep mine.

Mining scholars and industry workers have carried out a lot of research on strata control and overlying strata migration law of filling mining in coal seam more than 6 meters in deep mine [9–14]. Mining technology under buildings has been developed for many years at home and abroad and has made a lot of systematic achievements in rock stratum control and surface subsidence of paste backfilling mining [15–17], but the research on surface subsidence of paste backfilling mining with large mining height above 6 meters is relatively few. As for the surface subsidence problem of paste backfilling, the traditional methods are mostly based on the equivalent mining height and probability integral method to predict, without considering the limit conditions such as too low paste backfilling rate, too large undercutting amount, and insufficient support strength, which may also lead to roof fracture. In this case, there is no relevant special theoretical research on how to predict surface subsidence. In addition, there are few research studies on the combination of rock control and filling mining technology and method; that is, the current theoretical research lags behind the field practice of filling technology.

At present, with the deepening of the research on the influencing factors of paste backfill mining, many evaluation methods have been put forward, such as the analytic hierarchy process and the neural network method. However, it is difficult to effectively solve the problem of the combination of qualitative and quantitative indicators by sticking to a single evaluation method or model, and the evaluation results are often uncertain and fuzzy. The trend and direction of filling mining theory research and development should be based on the different production practices in the mining area. To solve the above problems, these methods are improved. Usually, different evaluation methods are combined and combined with their advantages to build a combination evaluation model.

Extenics is a new discipline that aims to study the possibility of business expansion and the rules and methods of innovation by using formal models. It was first proposed by Chinese scholar Cai et al. [18, 19]. It can unify quantitative analysis and qualitative analysis, which is suitable for mathematical measurement under the condition of limited information, and can reasonably reflect the development of evaluation objects. It is feasible to evaluate the influence factors of coal seam filling mining. In this paper, the fuzzy extension model is used to study the key factors of surface subsidence in deep mine thick seam paste backfilling mining, and the improvement measures are put forward.

2. Fuzzy Extension Analysis on Key Factors of Surface Subsidence in Paste Backfilling Mining

From the surface subsidence coefficient, we can see that the factors affecting the filling mining are affected by the mining technical conditions, roof and floor displacement, backfilling body parameters, and mining and backfilling technology, including both qualitative and quantitative indicators, and have variability, randomness, and fuzziness. In this section, combined with the characteristics of surface subsidence factors, a fuzzy extension analysis model will be established to find out the key factors and evaluate the stability and reliability of paste backfilling mining to control surface subsidence.

2.1. Evaluation Index of Key Factors of Surface Subsidence in Paste Backfilling Mining

The matter element system diagram of the key factors of surface subsidence in paste filling mining, namely, the evaluation index system, is shown in Figure 1.

Figure 1: Evaluation index of key factors of surface subsidence in paste backfilling mining.

2.2. Fuzzy Extension Evaluation Model

2.2.1. Fuzzy Matter Element

The concept of matter element is introduced by extenics, R = (thing, feature, quantity) = (N, c, ), in which thing N, feature c, and thing form a triple. Matter element correctly reflects the relationship between quality and quantity and is the logical cell of extenics [20–22]. If N is characterized by , , …, and the corresponding fuzzy values , , …, , then R is called n-dimensional fuzzy matter element, . If the n-dimensional matter elements of m objects are combined, the n-dimensional compound matter element Rmn of m objects is formed.

If the value of Rmn is rewritten as the value of fuzzy matter element, it is called the n-dimensional compound fuzzy matter element of m objects, recorded as

where Rmn is the n-dimensional compound fuzzy matter element of m objects, Mi is the i-th object (i = 1, 2, …, m), Ci is the j-th feature (j = 1, 2, … , n), and Vij is the membership degree corresponding to the j-th feature of the i-th object.

2.2.2. Fuzzy Membership Degree

To effectively determine the stability level of the key factors causing surface subsidence, its evaluation criteria must be established. It can be seen from Figure 1 that there are quantitative and qualitative indexes in the evaluation indexes of key factors of surface subsidence in paste filling mining, while the fuzzy membership function has great advantages in quantitative processing of qualitative indexes.

(1)Fuzzy Membership Degree of Qualitative Index The fuzzy membership function is determined by expert scoring to obtain the membership degree of the secondary qualitative indicators. Further, by virtue of the weight between the indicators, the comprehensive membership degree of the primary qualitative indicators is obtained according to the weighted average calculation results. People are subjective, so are experts. Then, in the actual scoring process, we usually take multiple (>5) experts to evaluate the results and then screen them. After removing the abnormal values, we take the average value as the result, to be consistent with the situation on-site as much as possible. The 10-point system is a common scoring method and means. Therefore, the specific scoring principles are shown in Table 1. As shown in the following formula, the fuzzy membership degree is calculated according to the characteristics of each index.

 where is the fuzzy membership degree, and and are the maximum and minimum values of the characteristic values of each evaluation index in each object. The comprehensive membership degree of the first level qualitative index is calculated as follows:

 where is the weight of the j-th secondary qualitative index of the i-th object and K is the number of secondary qualitative indexes.(2)Fuzzy Membership Degree of Quantitative Index The determination method is shown in Table 2.

Table 1: Point system.

Influence degree
Slight
Lighter
Secondary
Little serious
Serious

Score
0∼2
2∼4
4∼6
6∼8
8∼10

Table 2: Fuzzy membership degree of quantitative indicators.

No
Quantitative index
Fuzzy membership degree

1
Mining depth
When the mining depth is greater than 800, the membership degree is 0, less than 200 is 1, and the middle is 0.5.

2
Coal seam
When the dip angle of the coal seam is greater than 54°, the membership degree is 0; when the dip angle of the coal seam is 0°, the membership degree is 1.

3
Overburden structure
When the uniaxial compressive strength (unit: MPa) ∈ (40, 80) (hard overburden), the membership degree is 0.3; when ∈ (20, 40) (medium hard overburden), the membership degree is 0.5; when ∈ (10, 20) (weak overburden), the membership degree is 0.7; when < 10 (very weak overburden), the membership degree is 0.9.

4
Subsidence before filling
When the roof has no subsidence, the membership degree is taken as 1; when the roof is larger than 100 mm, it is taken as 0.

5
Nonbackfilling account
When there is no subsidence, the membership degree is taken as 1; when it is greater than 60 mm, it is taken as 0

6
Floating coal
When there is no compression, the membership degree is 1; when the compression is greater than 50 mm, the membership degree is 0.

7
Rock compression
When there is no compression, the membership degree is 1; when the compression is greater than 50 mm, the membership degree is 0.

8
Volume shrinkage
When there is no shrinkage, the membership degree is 1; when it is greater than 30 mm, it is 0.

9
Volume compression
When there is no compression, the membership degree is 1; when the compression is greater than 50 mm, the membership degree is 0.

10
Strength of filling body
When the strength is greater than 5 MPa, the membership degree is 0; when the strength is less than 1 MPa, the membership degree is 1.

11
Slicing height
When the slice mining height is greater than 5 m, the membership degree is 0; when the slice mining height is less than 1 m, the membership degree is 1.

Table 3: The meaning of proportion quotient.

Proportional scale
Meaning

1
The weights of ui and uj are the same

3
The weight of ui is greater than uj

5
The weight of ui is larger than that of uj

7
The weight of ui is stronger than that of uj

9
The weight of ui is significantly higher than that of uj

2, 4, 6, 8
The median of 1, 3, 5, 7, and 9

Reciprocal
If the ratio of ui and uj is aij, the ratio of uj and ui is 1/aij

Table 4: RI value table.

n
RI

1
0

2
0

3
0.52

4
0.89

5
1.12

6
1.26

7
1.36

8
1.41

9
1.46

10
1.49

11
1.52

12
1.54

13
1.56

14
1.58

15
1.59

Table 5: Evaluation guidelines.

Stability grades
Membership degree
Meaning of each grade

I
0.8∼1.0
Good stability

II
0.6∼0.8
Fairly good stability

III
0.4∼0.6
General stability

IV
0.2∼0.4
Relatively poor stability

V
0∼0.2
Poor stability

2.2.3. Association Analysis

(1)Determine the matter element to be evaluated

(a)Matter element to be evaluated

 where (i = 1, 2, …, n) is the quantity value of matter element p with respect to the i-th evaluation index ci.(b)Furthermore, the classical domain and node domain of the criterion layer can be obtained Order

 where R0j is the matter element and N0j is the j-th evaluation category; ci (i = 1, 2, …, n) is the i-th evaluation index; (j = 1, 2, …, m) is the value range of the i-th evaluation index ci in the j-th evaluation category, that is, the data range of the index ci classic domain, and there is . Order

 where P is the whole of the evaluation categories; (i = 1, 2, …, n) is the value range of P with respect to the i-th evaluation index ci, that is, the nodal region of P.

(2)The solution of correlation degree

 where is the lateral distance between point and interval with as the optimal value. The lateral distance is divided into left distance and right distance, which are calculated according to equations (8) and (9), respectively; and are the lateral distances between point and interval and , which are calculated according to equations (10) and (11) respectively.

(3)Determination of evaluation index weight for key factors of surface subsidence in paste backfilling mining The analytic hierarchy process is selected here to determine the weight. Let the element C of the upper layer be the criterion, let the elements of the next layer be u₁, u₂, …, un, and let the relative weights of each element to criterion C be ω₁, ω₂, …, ωn. The main process of determining the weight of the evaluation index based on the analytic hierarchy process (AHP) is as follows:

(a)Construction of pairwise comparison matrix According to the scale of 1–9, the importance degree of elements ui and uj is assigned to qualitatively reflect which is more important. The meaning of proportional scale is shown in Table 3. Through pairwise comparison, a paired comparison judgment matrix A = (aij)n×n is obtained, where aij is the proportional scale of the importance of elements ui and uj relative to C. The judgment matrix A has the following properties: aij > 0, aji = 1/aij, aii = 1.(b)Weight calculation by eigenvalue method Solve the eigenvalue of judgment matrix A; that is, the following formula is satisfied:

 where λmax is the largest eigenvalue of A, and W is the corresponding eigenvector. After normalization, W can be used as the weight vector [ω₁, ω₂, …, ωn]T.(c)Consistency test(i)    Calculate the consistency index (CI)

     where n is the order of matrix A, that is, the number of evaluation indexes.(ii)    Find the corresponding average random consistency index RI     Table 4 gives the average random consistency index of order 1–15.(iii)    Calculate the consistency ratio CR

     When CR is less than 0.1, the consistency of the judgment matrix is acceptable; otherwise, repeat steps (1)∼(3) after making an appropriate correction to the judgment matrix until the consistency test is met.

(4)Calculation of correlation degree Relevance degree K is a measure to describe the relevance between things. Calculate the correlation degree of matter element p to be evaluated with respect to grade j.

(5)Determine the evaluation level According to the expert experience and corresponding standards and based on the evaluation index value principle and value range, the stability grade of surface subsidence in paste backfilling mining is divided into five grades, namely, good stability (grade I), fairly good stability (grade II), general stability (grade III), relatively poor stability (grade IV), and poor stability (grade V). The evaluation criteria are shown in Table 5. Therefore, if

 then the matter element p belongs to grade j₀.

3. Study on the Influence of the Key Factors of Surface Subsidence in Paste Backfilling Mining

3.1. Analysis on Key Factors of Surface Subsidence in Paste Backfilling Mining

According to the above established fuzzy extension analysis model and evaluation index, we calculate the fuzzy membership degree of the index, establish the fuzzy matter element to be evaluated, determine the classical domain and node domain, determine the weight, calculate the correlation degree, and comprehensively evaluate the key factors and reliability of the surface subsidence of paste backfilling mining. Limited to the length of the article, the calculation process directly gives the weight value and the final overall evaluation; the results are shown in Table 6.

Table 6: Weight of evaluation indicators.

Criterion layer
Weight coefficient
Quantitative/primary qualitative indicators
Weight coefficient
Secondary qualitative indicators
Weight coefficient

Mining geological conditions
0.15
Mining depth
0.20
Original rock stress
0.54

Tectonic stress
0.46

Coal seam
0.21
Coal seam dip angle
0.48

Coal seam thickness
0.52

Surface water drainage
0.23

Overburden structure
0.36
Strength of roadway
0.31

Key layer location
0.34

Topsoil thickness
0.35

Approach of roof and floor
0.28
Subsidence before filling
0.36
Roof subsidence
0.66

Floor heave
0.27

Stress expansion
0.07

Nonbackfilling account
0.47

Floating coal
0.12

Rock compression
0.05

Filling body
0.32
Volume shrinkage
0.17
Bleeding rate
0.45

Water content
0.55

Volume compression
0.35
Bulk concentration
0.53

Volumetric strain
0.47

Strength of filling body
0.48
Poisson’s ratio
0.49

Elastic modulus
0.51

Mining and filling technology
0.17
Mining technology
0.33
Mining speed
0.36

Support mode
0.38

Coal pillar strength
0.26

Filling technology
0.47
Filling step
0.40

Unsupported roof span
0.42

Filling sequence
0.18

Slicing height
0.20

Overburden retained voids
0.08
Mining method
0.38
Strip mining
0.67

Longwall mining
0.33

Stratigraphic structure
0.34

Surface structure
0.28

Through the calculation of single index and multi-index correlation degree, the correlation degree of mining geological conditions, roof and floor approach, mining and backfilling technology, and overburden retained gap of the criterion layer and the target layer is determined. The results show that the stability and reliability of surface subsidence of paste backfilling mining are “good,” and the key factors are the overburden structure, the roof subsidence before backfilling, nonbackfilling account, the strength of the filling body, and the backfilling technology. At the same time, it is also the main reason that the control stability of surface subsidence is “good.”

3.2. Technical Path and Improvement Measures

To sum up, combined with theoretical research and production practice, the technical path of controlling surface subsidence by backfilling mining is mainly to control the four aspects of mining geological conditions, roof and floor approaching amount, filling body, and overburden retained gap, especially the five key factors of overburden structure, roof subsidence before backfilling, nonbackfilling account, filling body strength, and backfilling technology. But for example, in mining geological conditions, such as overburden structure, it is a given factor for a certain mine. Although mining depth, coal seam thickness, key layer position, and so on all have greater weight, which has a great impact on surface subsidence, they are all internal factors that cannot be changed. Therefore, this paper proposes to control the surface subsidence of paste backfilling mining from the following technical approaches and external key factors affecting the movement of overlying strata.

3.2.1. Control of Roof Subsidence after Mining and before Paste Backfilling

The roof subsidence before backfilling after mining is composed of three parts: roof subsidence, floor heave, and stress expansion of roof and floor rock strata, which is mainly influenced by time interval, roof pressure, unsupported roof span, and working face length. At present, in the paste backfilling mining, we should try to shorten the time interval after mining and before backfilling, shorten unsupported roof span, improve the active support force of the support, parallel mining and backfilling operation as far as possible, and avoid designing long working face mining to reduce and control the roof subsidence before backfilling after mining.

3.2.2. Increase of Backfilling Rate

The backfilling rate of goaf is divided into plane backfilling rate and top connection rate of the filling body. To increase the backfilling rate is to increase the backfilling plane to make it fully filled and to increase the height of contact with the roof to make it fully contacted. The lower the plane backfilling rate, especially the honeycomb backfilling, the more serious the surface subsidence. The reason for the low amount of contact between the backfill and roof is mostly due to the limitation of low backfilling technology and insufficient backfilling materials.

Increasing the backfilling rate is the most effective measure to reduce the surface subsidence, but it will increase the cost of backfilling mining. In the actual production, improving the backfilling technology, strengthening the monitoring, and even secondary grouting can avoid honeycomb backfilling, make the filling material fill the whole roof gap, reduce the potential roof subsidence space, and make the filling body fully supported to achieve the purpose of reducing subsidence control.

For the compensation of nonbackfilling account, the sealing of filling space in goaf should be strengthened in time to improve the liquid level of filling slurry. In addition, the support strength should be ensured to prevent the phenomenon of sinking net pocket and false top roof connection. For the new backfilling mining face, advance planning, reasonable design of roadway location, support mode, and coal pillar size, the working face implement inclined mining and inclined backfilling.

3.2.3. Improvement of the Physical and Mechanical Properties of Backfilling Slurry and Filling Body

In paste backfilling mining, the slurry is required to be free from precipitation, segregation, and bleeding. To improve the physical and mechanical properties of backfilling slurry and filling body, it is necessary to reduce the shrinkage and compression of the filling body. The methods to reduce the shrinkage and compression of the filling body are to reduce the bleeding rate and volume strain of the filling slurry and increase the volume concentration of the filling slurry, which is also the form of improving the physical and mechanical properties of the filling body. Therefore, it is necessary to study the material ratio, pumpability, filling equipment, and conveying distance.

3.2.4. Optimization of Backfilling and Production Process

Reasonable backfilling step and roof control distance are very important to the stability of the roof. The minimum roof control distance depends on the mining ventilation conditions, and the maximum roof control distance depends on the backfilling step. The longer the roof control time is, the longer the roof hanging time is, and the greater the roof subsidence is. Under the condition of three-mining- and one-backfilling-cycle mode, the roof of the area to be filled has become a completely tensile state, and the risk of roof damage increases. Therefore, for the paste backfilling with 6 m large mining height and slicing mining, two-mining- and one-backfilling-cycle mode is reasonable.

In addition, there are some subkey factors to control the surface subsidence of paste backfilling mining. For example, the height of slicing mining is optimized. When slicing mining, even if the mining height is 4 m, the maximum roof approaching amount is about 300 mm, which is also within the controllable range. In the later stage, the mining height can be adjusted appropriately according to the site conditions to increase the production of the working face. For example, strip mining instead of longwall mining can reduce the surface subsidence rate by 7.5∼66.7%, which makes more void volume stay in the overlying strata and controls the surface subsidence obviously. However, the compression of the roof and floor rock and floor floating coal can be ignored.

4. Conclusions

The research on the key factors affecting the surface subsidence of thick seam paste filling mining in deep mines is a complicated system engineering, which is a new exploration to liberate the coal resources under buildings combined with the mining conditions of the Lu’an mining area and has very important strategic significance. The main conclusions are as follows:

(1)A fuzzy extension model for the stability of surface subsidence caused by paste backfilling mining is established, and many factors of surface subsidence caused by paste backfilling mining are analyzed. The limitation of the traditional model which makes evaluation only from a qualitative or quantitative point of view is overcome, and the reliability grade of controlling surface subsidence caused by backfilling mining is “good.” The key factors to control the surface subsidence are the overburden structure, the amount of roof subsidence before backfilling, nonbackfilling account, the strength of the filling body, and the backfilling technology. In addition, optimization suggestions in design and mining are given.(2)It is preliminarily estimated that paste backfilling can make the recovery rate of 3# coal close to 80% in Gaohe coal mine, increase the recovery rate close to 40%, liberate more than 210 million tons of coal resources under buildings, effectively extend the service life of the mine, and increase employment opportunities. Paste backfilling can effectively reduce the roof and floor damage scope and gas emission in goaf and improve the safety guarantee of coal mining. The main raw material of backfilling mining is coal gangue, which reduces the occupied land area of gangue while solving the problem of gangue discharge. Backfilling mining can protect groundwater resources and promote the ecological environment construction of mining areas.

Data Availability

The data used to support the findings of this study are included within the article. All the data are obtained through experiments and tested by our research team in Gaohe Mine and laboratory and are true and effective. The right to using data belongs to the authors before the article gets published, but the data can be referenced after publishing,.

Conflicts of Interest

The authors declare no conflicts of interest.

Acknowledgments

The authors would like to express their sincere thanks to all those who have helped writing this paper. They have done a lot of preliminary work in the related projects of the paper. Here the authors express their sincere thanks. Without their help, it would be much harder for the authors to finish this study. This paper was supported by Jiangsu University Natural Science Research Project (20KJB560032) and Jiangsu Construction System Science and Technology Project (Guidance) (2020ZD30, 2019ZD080, and 2019ZD070)

References

T. Liu, Backfilling Mining Technology and Its Application, Metallurgical Industry Press, Beijing, China, 2001.

A. Zhou, Cemented Backfilling of Mine Waste, Metallurgical Industry Press, Beijing, China, 2007.

M. Zhang, Theoretical Study on Overburden Control of High Efficiency Paste Backfilling Working Face, China University of Mining and Technology, Xuzhou, China, 2019.

X. Huang, “Leading energy industry reform with green development,” People’s Forum, vol. 31, pp. 26-27, 2017.

R. Cowling, “Twenty-five years of mine filling-developments and directions,” in Proceedings of the 6th International Symposium on Mining with Backfill, Brislane, Australia, Apri1 1998.

M. Cullen, Geotechnical Studies of Retreat Pillar Coal Mining at Mining at Shallow Depth, McGill University, Montreal, Canada, 2002.

K. B. Singh and T. N. Singh, “Ground movements over long wall workings in the Kamptee coalfield, India,” Engineering Geology, vol. 50, no. 1, pp. 125–139, 1998.

X. L. Yao, B. N. Whittaker, and D. J. Reddish, “Influence of overburden mass behavioral properties on subsidence limit characteristics,” Mining Science and Technology, vol. 13, no. 2, pp. 167–173, 1991.

M. Qian and P. Shi, Ground Pressure and Strata Control, China University of Mining and Technology Press, Xuzhou, China, 2003.

M. Qian, “Structural model of overlying strata and its application,” Journal of China University of Mining and Technology, vol. 2, no. 2, pp. 6–16, 1982.

A. Philip, Engineering Design of Paste Backfill Systems, Queen’s University Kingston, Kingston, Canada, 2003.

X.-x. Miao, J.-x. Zhang, and M.-m. Feng, “Waste-filling in fully-mechanized coal mining and its application,” Journal of China University of Mining and Technology, vol. 18, no. 4, pp. 479–482, 2008.

Z. Song and J. Jiang, “Research focus and direction of coal mine strata control,” Journal of Rock Mechanics and Engineering, vol. 15, no. 2, pp. 33–39, 1996.

S. Widisinghe and N. Sivakugan, “Vertical stress isobars for trenches and mine stopes containing granular backfills,” International Journal of Geomechanics, vol. 14, no. 2, pp. 313–318, 2014.

G. He and L. Yang, Mining Subsidence, China University of Mining and Technology Press, Xuzhou, China, 1994.

L. Wu and J. Wang, Theory and Practice of Strip Mining under Buildings, China University of Mining and Technology Press, Xuzhou, China, 1994.

K. Deng, “Application of rock mass structure effect in mining subsidence,” China University of Mining and Technology, Xuzhou, China, 1993, Doctoral dissertation.

W. Cai, C. Yang, and B. He, “Some problems in extenics research,” Journal of Guangdong University of technology, vol. 18, no. 1, pp. 1–5, 2001.

Y. Wang, “Extension fuzzy evaluation method for rock mass quality of underground mining engineering,” Journal of Rock Mechanics and Engineering, vol. 21, no. 1, pp. 18–22, 2002.

W. Cai, “The extension set and incompatibility problems,” Journal of Scientific Exploration, vol. 1, no. 1, pp. 81–93, 1983.

M. Kamber and J. Han, Data Mining Concepts and Techniques, Morgan Kaufmann Publishers, San Francisco, CA, USA, 2001.

G. Stoilos, G. Stamou, J. Z. Pan, V. Tzouvaras, and I. Horrocks, “Reasoning with very expressive fuzzy description logics,” Journal of Artificial Intelligence Research, vol. 30, no. 8, pp. 273–320, 2007.