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

Authors: Meng Zhang; Hui He; Xin Jin; Yonglong Qu; Hongjun Guo
Journal: Advances in Civil Engineering
 (2021)

---

## 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.


---

## Body

## 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 Figure1.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, 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–22]. If N is characterized by c1, c2, …, cn and the corresponding fuzzy values v1, v2, …, vn, then R is called n-dimensional fuzzy matter element, R=N,c,v. 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 ofRmn is rewritten as the value of fuzzy matter element, it is called the n-dimensional compound fuzzy matter element of m objects, recorded as(1)Rmn=M1M2⋅⋅⋅Mmc1v11v21⋅⋅⋅vm1c2v12v22⋅⋅⋅vm2⋅⋅⋅⋅⋅⋅⋅⋅⋅⋅⋅⋅⋅⋅⋅cnv1nv2n⋅⋅⋅vmn,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 Figure1 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 IndexThe 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 Table1.As shown in the following formula, the fuzzy membership degree is calculated according to the characteristics of each index.(2)μij=μijmaxμij,the bigger, the better,μij=minμijμij,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.The comprehensive membership degree of the first level qualitative index is calculated as follows:(3)μ=∑i=1kωij⋅μij,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 IndexThe determination method is shown in Table2.Table 1
Point system.
Influence degreeSlightLighterSecondaryLittle seriousSeriousScore0∼22∼44∼66∼88∼10Table 2
Fuzzy membership degree of quantitative indicators.
NoQuantitative indexFuzzy membership degree1Mining depthWhen the mining depth is greater than 800, the membership degree is 0, less than 200 is 1, and the middle is 0.5.2Coal seamWhen 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.3Overburden structureWhen 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.4Subsidence before fillingWhen 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.5Nonbackfilling accountWhen there is no subsidence, the membership degree is taken as 1; when it is greater than 60 mm, it is taken as 06Floating coalWhen there is no compression, the membership degree is 1; when the compression is greater than 50 mm, the membership degree is 0.7Rock compressionWhen there is no compression, the membership degree is 1; when the compression is greater than 50 mm, the membership degree is 0.8Volume shrinkageWhen there is no shrinkage, the membership degree is 1; when it is greater than 30 mm, it is 0.9Volume compressionWhen there is no compression, the membership degree is 1; when the compression is greater than 50 mm, the membership degree is 0.10Strength of filling bodyWhen 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.11Slicing heightWhen 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 scaleMeaning1The weights ofui and uj are the same3The weight ofui is greater than uj5The weight ofui is larger than that of uj7The weight ofui is stronger than that of uj9The weight ofui is significantly higher than that of uj2, 4, 6, 8The median of 1, 3, 5, 7, and 9ReciprocalIf the ratio ofui and uj is aij, the ratio of uj and ui is 1/aijTable 4
RI value table.
nRI102030.5240.8951.1261.2671.3681.4191.46101.49111.52121.54131.56141.58151.59Table 5
Evaluation guidelines.
Stability gradesMembership degreeMeaning of each gradeI0.8∼1.0Good stabilityII0.6∼0.8Fairly good stabilityIII0.4∼0.6General stabilityIV0.2∼0.4Relatively poor stabilityV0∼0.2Poor stability


#### 2.2.3. Association Analysis

(1)
Determine the matter element to be evaluated(a)
Matter element to be evaluated(4)R=pc1v1c2v2⋯⋯cnvn,wherevi (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 obtainedOrder(5)R0j=N0j,C,V0j=N0jc1V0j1c2V0j2⋯⋯cnV0jn=N0jc1a0j1,b0j1c2a0j2,b0j2⋯⋯cna0jn,b0jn,whereR0j is the matter element and N0j is the j-th evaluation category; ci (i = 1, 2, …, n) is the i-th evaluation index; V0ji (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 V0ji=a0ji,b0ji.Order(6)RP=P,C,VP=Pc1VP1c2VP1⋯⋯cnVPn=Pc1aP1,bP1c2aP2,bP2⋯⋯cnaPn,bPn,whereP is the whole of the evaluation categories; VPi (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(7)Kjvi=−ρvi,v0ji,V0jib0ji−a0ji,vi∈V0ji,ρvi,v0ji,V0jiρvi,VPi−ρvi,V0ji,vi∉V0ji,whereρvi,v0ji,V0ji is the lateral distance between point vi and interval v0ji with V0ji 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; ρvi,V0ji and ρvi,VPi are the lateral distances between point vi and interval V0ji and VPi, which are calculated according to equations (10) and (11) respectively.(8)ρvi,v0ji,V0ji=a0ji−vi,vi≤a0ji,vi−b0ji,vi≥v0ji,b0ji−v0ji⋅vi−a0jia0ji−v0ji,vi∈a0ji,v0ji,(9)ρvi,v0ji,V0ji=a0ji−vi,vi≤v0ji,vi−b0ji,vi≥b0ji,a0ji−v0ji⋅b0ji−vib0ji−v0ji,vi∈b0ji,v0ji,(10)ρvi,V0ji=vi−a0ji+b0ji2−b0ji−a0ji2,(11)ρvi,VPi=vi−aPi+bPi2−bPi−aPi2.(3)
Determination of evaluation index weight for key factors of surface subsidence in paste backfilling miningThe analytic hierarchy process is selected here to determine the weight.Let the elementC of the upper layer be the criterion, let the elements of the next layer be u1, u2, …, un, and let the relative weights of each element to criterion C be ω1, ω2, …, ω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 matrixAccording to the scale of 1–9, the importance degree of elementsui 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 matrixA = (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 methodSolve the eigenvalue of judgment matrixA; that is, the following formula is satisfied:(12)AW=λmaxW,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)(13)CI=λmax−nn−1,wheren is the order of matrix A, that is, the number of evaluation indexes.(ii)
Find the corresponding average random consistency index RITable4 gives the average random consistency index of order 1–15.(iii)
Calculate the consistency ratio CR(14)CR=CIRI.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 degreeRelevance degreeK 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.(15)Kjp=∑i=1nωiKjvi.(5)
Determine the evaluation levelAccording 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 Table5.Therefore, if(16)Kj0p=maxj0∈1,2,…,mKjp,then the matter elementp belongs to grade j0.


## 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 Figure1.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, 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–22]. If N is characterized by c1, c2, …, cn and the corresponding fuzzy values v1, v2, …, vn, then R is called n-dimensional fuzzy matter element, R=N,c,v. 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 ofRmn is rewritten as the value of fuzzy matter element, it is called the n-dimensional compound fuzzy matter element of m objects, recorded as(1)Rmn=M1M2⋅⋅⋅Mmc1v11v21⋅⋅⋅vm1c2v12v22⋅⋅⋅vm2⋅⋅⋅⋅⋅⋅⋅⋅⋅⋅⋅⋅⋅⋅⋅cnv1nv2n⋅⋅⋅vmn,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 Figure1 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 IndexThe 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 Table1.As shown in the following formula, the fuzzy membership degree is calculated according to the characteristics of each index.(2)μij=μijmaxμij,the bigger, the better,μij=minμijμij,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.The comprehensive membership degree of the first level qualitative index is calculated as follows:(3)μ=∑i=1kωij⋅μij,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 IndexThe determination method is shown in Table2.Table 1
Point system.
Influence degreeSlightLighterSecondaryLittle seriousSeriousScore0∼22∼44∼66∼88∼10Table 2
Fuzzy membership degree of quantitative indicators.
NoQuantitative indexFuzzy membership degree1Mining depthWhen the mining depth is greater than 800, the membership degree is 0, less than 200 is 1, and the middle is 0.5.2Coal seamWhen 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.3Overburden structureWhen 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.4Subsidence before fillingWhen 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.5Nonbackfilling accountWhen there is no subsidence, the membership degree is taken as 1; when it is greater than 60 mm, it is taken as 06Floating coalWhen there is no compression, the membership degree is 1; when the compression is greater than 50 mm, the membership degree is 0.7Rock compressionWhen there is no compression, the membership degree is 1; when the compression is greater than 50 mm, the membership degree is 0.8Volume shrinkageWhen there is no shrinkage, the membership degree is 1; when it is greater than 30 mm, it is 0.9Volume compressionWhen there is no compression, the membership degree is 1; when the compression is greater than 50 mm, the membership degree is 0.10Strength of filling bodyWhen 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.11Slicing heightWhen 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 scaleMeaning1The weights ofui and uj are the same3The weight ofui is greater than uj5The weight ofui is larger than that of uj7The weight ofui is stronger than that of uj9The weight ofui is significantly higher than that of uj2, 4, 6, 8The median of 1, 3, 5, 7, and 9ReciprocalIf the ratio ofui and uj is aij, the ratio of uj and ui is 1/aijTable 4
RI value table.
nRI102030.5240.8951.1261.2671.3681.4191.46101.49111.52121.54131.56141.58151.59Table 5
Evaluation guidelines.
Stability gradesMembership degreeMeaning of each gradeI0.8∼1.0Good stabilityII0.6∼0.8Fairly good stabilityIII0.4∼0.6General stabilityIV0.2∼0.4Relatively poor stabilityV0∼0.2Poor stability


### 2.2.3. Association Analysis

(1)
Determine the matter element to be evaluated(a)
Matter element to be evaluated(4)R=pc1v1c2v2⋯⋯cnvn,wherevi (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 obtainedOrder(5)R0j=N0j,C,V0j=N0jc1V0j1c2V0j2⋯⋯cnV0jn=N0jc1a0j1,b0j1c2a0j2,b0j2⋯⋯cna0jn,b0jn,whereR0j is the matter element and N0j is the j-th evaluation category; ci (i = 1, 2, …, n) is the i-th evaluation index; V0ji (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 V0ji=a0ji,b0ji.Order(6)RP=P,C,VP=Pc1VP1c2VP1⋯⋯cnVPn=Pc1aP1,bP1c2aP2,bP2⋯⋯cnaPn,bPn,whereP is the whole of the evaluation categories; VPi (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(7)Kjvi=−ρvi,v0ji,V0jib0ji−a0ji,vi∈V0ji,ρvi,v0ji,V0jiρvi,VPi−ρvi,V0ji,vi∉V0ji,whereρvi,v0ji,V0ji is the lateral distance between point vi and interval v0ji with V0ji 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; ρvi,V0ji and ρvi,VPi are the lateral distances between point vi and interval V0ji and VPi, which are calculated according to equations (10) and (11) respectively.(8)ρvi,v0ji,V0ji=a0ji−vi,vi≤a0ji,vi−b0ji,vi≥v0ji,b0ji−v0ji⋅vi−a0jia0ji−v0ji,vi∈a0ji,v0ji,(9)ρvi,v0ji,V0ji=a0ji−vi,vi≤v0ji,vi−b0ji,vi≥b0ji,a0ji−v0ji⋅b0ji−vib0ji−v0ji,vi∈b0ji,v0ji,(10)ρvi,V0ji=vi−a0ji+b0ji2−b0ji−a0ji2,(11)ρvi,VPi=vi−aPi+bPi2−bPi−aPi2.(3)
Determination of evaluation index weight for key factors of surface subsidence in paste backfilling miningThe analytic hierarchy process is selected here to determine the weight.Let the elementC of the upper layer be the criterion, let the elements of the next layer be u1, u2, …, un, and let the relative weights of each element to criterion C be ω1, ω2, …, ω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 matrixAccording to the scale of 1–9, the importance degree of elementsui 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 matrixA = (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 methodSolve the eigenvalue of judgment matrixA; that is, the following formula is satisfied:(12)AW=λmaxW,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)(13)CI=λmax−nn−1,wheren is the order of matrix A, that is, the number of evaluation indexes.(ii)
Find the corresponding average random consistency index RITable4 gives the average random consistency index of order 1–15.(iii)
Calculate the consistency ratio CR(14)CR=CIRI.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 degreeRelevance degreeK 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.(15)Kjp=∑i=1nωiKjvi.(5)
Determine the evaluation levelAccording 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 Table5.Therefore, if(16)Kj0p=maxj0∈1,2,…,mKjp,then the matter elementp belongs to grade j0.


## 2.2.1. Fuzzy Matter Element

The 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–22]. If N is characterized by c1, c2, …, cn and the corresponding fuzzy values v1, v2, …, vn, then R is called n-dimensional fuzzy matter element, R=N,c,v. 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 ofRmn is rewritten as the value of fuzzy matter element, it is called the n-dimensional compound fuzzy matter element of m objects, recorded as(1)Rmn=M1M2⋅⋅⋅Mmc1v11v21⋅⋅⋅vm1c2v12v22⋅⋅⋅vm2⋅⋅⋅⋅⋅⋅⋅⋅⋅⋅⋅⋅⋅⋅⋅cnv1nv2n⋅⋅⋅vmn,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 Figure1 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 IndexThe 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 Table1.As shown in the following formula, the fuzzy membership degree is calculated according to the characteristics of each index.(2)μij=μijmaxμij,the bigger, the better,μij=minμijμij,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.The comprehensive membership degree of the first level qualitative index is calculated as follows:(3)μ=∑i=1kωij⋅μij,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 IndexThe determination method is shown in Table2.Table 1
Point system.
Influence degreeSlightLighterSecondaryLittle seriousSeriousScore0∼22∼44∼66∼88∼10Table 2
Fuzzy membership degree of quantitative indicators.
NoQuantitative indexFuzzy membership degree1Mining depthWhen the mining depth is greater than 800, the membership degree is 0, less than 200 is 1, and the middle is 0.5.2Coal seamWhen 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.3Overburden structureWhen 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.4Subsidence before fillingWhen 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.5Nonbackfilling accountWhen there is no subsidence, the membership degree is taken as 1; when it is greater than 60 mm, it is taken as 06Floating coalWhen there is no compression, the membership degree is 1; when the compression is greater than 50 mm, the membership degree is 0.7Rock compressionWhen there is no compression, the membership degree is 1; when the compression is greater than 50 mm, the membership degree is 0.8Volume shrinkageWhen there is no shrinkage, the membership degree is 1; when it is greater than 30 mm, it is 0.9Volume compressionWhen there is no compression, the membership degree is 1; when the compression is greater than 50 mm, the membership degree is 0.10Strength of filling bodyWhen 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.11Slicing heightWhen 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 scaleMeaning1The weights ofui and uj are the same3The weight ofui is greater than uj5The weight ofui is larger than that of uj7The weight ofui is stronger than that of uj9The weight ofui is significantly higher than that of uj2, 4, 6, 8The median of 1, 3, 5, 7, and 9ReciprocalIf the ratio ofui and uj is aij, the ratio of uj and ui is 1/aijTable 4
RI value table.
nRI102030.5240.8951.1261.2671.3681.4191.46101.49111.52121.54131.56141.58151.59Table 5
Evaluation guidelines.
Stability gradesMembership degreeMeaning of each gradeI0.8∼1.0Good stabilityII0.6∼0.8Fairly good stabilityIII0.4∼0.6General stabilityIV0.2∼0.4Relatively poor stabilityV0∼0.2Poor stability


## 2.2.3. Association Analysis

(1)
Determine the matter element to be evaluated(a)
Matter element to be evaluated(4)R=pc1v1c2v2⋯⋯cnvn,wherevi (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 obtainedOrder(5)R0j=N0j,C,V0j=N0jc1V0j1c2V0j2⋯⋯cnV0jn=N0jc1a0j1,b0j1c2a0j2,b0j2⋯⋯cna0jn,b0jn,whereR0j is the matter element and N0j is the j-th evaluation category; ci (i = 1, 2, …, n) is the i-th evaluation index; V0ji (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 V0ji=a0ji,b0ji.Order(6)RP=P,C,VP=Pc1VP1c2VP1⋯⋯cnVPn=Pc1aP1,bP1c2aP2,bP2⋯⋯cnaPn,bPn,whereP is the whole of the evaluation categories; VPi (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(7)Kjvi=−ρvi,v0ji,V0jib0ji−a0ji,vi∈V0ji,ρvi,v0ji,V0jiρvi,VPi−ρvi,V0ji,vi∉V0ji,whereρvi,v0ji,V0ji is the lateral distance between point vi and interval v0ji with V0ji 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; ρvi,V0ji and ρvi,VPi are the lateral distances between point vi and interval V0ji and VPi, which are calculated according to equations (10) and (11) respectively.(8)ρvi,v0ji,V0ji=a0ji−vi,vi≤a0ji,vi−b0ji,vi≥v0ji,b0ji−v0ji⋅vi−a0jia0ji−v0ji,vi∈a0ji,v0ji,(9)ρvi,v0ji,V0ji=a0ji−vi,vi≤v0ji,vi−b0ji,vi≥b0ji,a0ji−v0ji⋅b0ji−vib0ji−v0ji,vi∈b0ji,v0ji,(10)ρvi,V0ji=vi−a0ji+b0ji2−b0ji−a0ji2,(11)ρvi,VPi=vi−aPi+bPi2−bPi−aPi2.(3)
Determination of evaluation index weight for key factors of surface subsidence in paste backfilling miningThe analytic hierarchy process is selected here to determine the weight.Let the elementC of the upper layer be the criterion, let the elements of the next layer be u1, u2, …, un, and let the relative weights of each element to criterion C be ω1, ω2, …, ω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 matrixAccording to the scale of 1–9, the importance degree of elementsui 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 matrixA = (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 methodSolve the eigenvalue of judgment matrixA; that is, the following formula is satisfied:(12)AW=λmaxW,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)(13)CI=λmax−nn−1,wheren is the order of matrix A, that is, the number of evaluation indexes.(ii)
Find the corresponding average random consistency index RITable4 gives the average random consistency index of order 1–15.(iii)
Calculate the consistency ratio CR(14)CR=CIRI.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 degreeRelevance degreeK 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.(15)Kjp=∑i=1nωiKjvi.(5)
Determine the evaluation levelAccording 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 Table5.Therefore, if(16)Kj0p=maxj0∈1,2,…,mKjp,then the matter elementp belongs to grade j0.


## 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 Table6.Table 6
Weight of evaluation indicators.
Criterion layerWeight coefficientQuantitative/primary qualitative indicatorsWeight coefficientSecondary qualitative indicatorsWeight coefficientMining geological conditions0.15Mining depth0.20Original rock stress0.54Tectonic stress0.46Coal seam0.21Coal seam dip angle0.48Coal seam thickness0.52Surface water drainage0.23Overburden structure0.36Strength of roadway0.31Key layer location0.34Topsoil thickness0.35Approach of roof and floor0.28Subsidence before filling0.36Roof subsidence0.66Floor heave0.27Stress expansion0.07Nonbackfilling account0.47Floating coal0.12Rock compression0.05Filling body0.32Volume shrinkage0.17Bleeding rate0.45Water content0.55Volume compression0.35Bulk concentration0.53Volumetric strain0.47Strength of filling body0.48Poisson’s ratio0.49Elastic modulus0.51Mining and filling technology0.17Mining technology0.33Mining speed0.36Support mode0.38Coal pillar strength0.26Filling technology0.47Filling step0.40Unsupported roof span0.42Filling sequence0.18Slicing height0.20Overburden retained voids0.08Mining method0.38Strip mining0.67Longwall mining0.33Stratigraphic structure0.34Surface structure0.28Through 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.


## 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 Table6.Table 6
Weight of evaluation indicators.
Criterion layerWeight coefficientQuantitative/primary qualitative indicatorsWeight coefficientSecondary qualitative indicatorsWeight coefficientMining geological conditions0.15Mining depth0.20Original rock stress0.54Tectonic stress0.46Coal seam0.21Coal seam dip angle0.48Coal seam thickness0.52Surface water drainage0.23Overburden structure0.36Strength of roadway0.31Key layer location0.34Topsoil thickness0.35Approach of roof and floor0.28Subsidence before filling0.36Roof subsidence0.66Floor heave0.27Stress expansion0.07Nonbackfilling account0.47Floating coal0.12Rock compression0.05Filling body0.32Volume shrinkage0.17Bleeding rate0.45Water content0.55Volume compression0.35Bulk concentration0.53Volumetric strain0.47Strength of filling body0.48Poisson’s ratio0.49Elastic modulus0.51Mining and filling technology0.17Mining technology0.33Mining speed0.36Support mode0.38Coal pillar strength0.26Filling technology0.47Filling step0.40Unsupported roof span0.42Filling sequence0.18Slicing height0.20Overburden retained voids0.08Mining method0.38Strip mining0.67Longwall mining0.33Stratigraphic structure0.34Surface structure0.28Through 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.


## 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.



---

*Source: 6634331-2021-06-07.xml*