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

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.


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][2][3][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][6][7][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. erefore, 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][10][11][12][13][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][16][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. e 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.

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.

Evaluation Index of Key Factors of Surface Subsidence in Paste Backfilling
Mining. e 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.

Fuzzy Matter Element.
e concept of matter element is introduced by extenics, R � (thing, feature, quantity) � (N, c, v), in which thing N, feature c, and thing v form a triple. Matter element correctly reflects the relationship between quality and quantity and is the logical cell of extenics [20][21][22]. If N is characterized by c 1 , c 2 , . . ., c n and the corresponding fuzzy values v 1 , v 2 , . . ., v n , then R is called n-dimensional fuzzy matter element, R � (N, c, v). If the ndimensional matter elements of m objects are combined, the n-dimensional compound matter element R mn of m objects is formed.
If the value of R mn 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 R mn is the n-dimensional compound fuzzy matter element of m objects, M i is the i-th object (i � 1, 2, . . ., m), C i is the j-th feature (j � 1, 2, . . . , n), and V ij is the membership degree corresponding to the j-th feature of the i-th object.

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 e 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. en, 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. e 10-point system is a common scoring method and means. erefore, 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.
, the bigger, the better, , the smaller, the better, where μ ij is the fuzzy membership degree, and max μ ij and min μ ij are the maximum and minimum values of the characteristic values of each evaluation index in each object. e comprehensive membership degree of the first level qualitative index is calculated as follows: where ω ij 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 e determination method is shown in Table 2.

Association Analysis
(1) Determine the matter element to be evaluated (a) Matter element to be evaluated where v i (i � 1, 2, . . ., n) is the quantity value of matter element p with respect to the i-th evaluation index c i . (b) Furthermore, the classical domain and node domain of the criterion layer can be obtained Order   where R 0j is the matter element and N 0j is the j-th evaluation category; is the value range of the i-th evaluation index c i in the j-th evaluation category, that is, the data range of the index c i classic domain, and there is where P is the whole of the evaluation categories; V Pi (i � 1, 2, . . ., n) is the value range of P with respect to the i-th evaluation index c i , that is, the nodal region of P.
(2) e solution of correlation degree is the lateral distance between point v i and interval v 0ji with V 0ji as the optimal value. e lateral distance is divided into left distance and right distance, which are calculated according to equations (8) and (9), respectively; ρ(v i , V 0ji ) and ρ(v i , V Pi ) are the lateral distances between point v i and interval V 0ji and V Pi , which are calculated according to equations (10) and (11) respectively. 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.

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.

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.

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.

Advances in Civil Engineering
(3) Determination of evaluation index weight for key factors of surface subsidence in paste backfilling mining e 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 1 , u 2 , . . ., u n , and let the relative weights of each element to criterion C be ω 1 , ω 2 , . . ., ω n . e 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 u i and u j is assigned to qualitatively reflect which is more important. e meaning of proportional scale is shown in Table 3.
rough pairwise comparison, a paired comparison judgment matrix A � (a ij ) n×n is obtained, where a ij is the proportional scale of the importance of elements u i and u j relative to C. e judgment matrix A has the following properties: a ij > 0, a ji � 1/a ij , a ii � 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 [ω 1 , ω 2 , . . ., ω 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). e evaluation criteria are shown in Table 5. erefore, if then the matter element p belongs to grade j 0 .

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 e weight of u i is greater than u j 5 e weight of u i is larger than that of u j 7 e weight of u i is stronger than that of u j 9 e weight of u i is significantly higher than that of u j 2, 4, 6, 8 e median of 1, 3, 5, 7, and 9 Reciprocal If the ratio of u i and u j is a ij , the ratio of u j and u i is 1/a ij 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. rough 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. e 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."

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. erefore, 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.

Control of Roof Subsidence after Mining and before Paste Backfilling.
e 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.

Increase of Backfilling Rate.
e 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. e lower the plane backfilling rate, especially the honeycomb backfilling, the more serious the surface subsidence. e 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.

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. e 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. erefore, it is necessary to study the material ratio, pumpability, filling equipment, and conveying distance.

Optimization of Backfilling and Production Process.
Reasonable backfilling step and roof control distance are very important to the stability of the roof. e minimum roof control distance depends on the mining ventilation conditions, and the maximum roof control distance depends on the backfilling step. e 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 Advances in Civil Engineering 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. erefore, for the paste backfilling with 6 m large mining height and slicing mining, two-miningand 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.

Conclusions
e 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. e 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. e 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." e 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 Advances in Civil Engineering 7 the roof and floor damage scope and gas emission in goaf and improve the safety guarantee of coal mining. e 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
e 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. e right to using data belongs to the authors before the article gets published, but the data can be referenced after publishing,.

Conflicts of Interest
e authors declare no conflicts of interest.