Next Article in Journal
Prospective Smart Distribution Substation in Bangladesh: Modeling and Analysis
Next Article in Special Issue
Factors Influencing Pile Friction Bearing Capacity: Proposing a Novel Procedure Based on Gradient Boosted Tree Technique
Previous Article in Journal
Linking Sustainability, Embeddedness and Marketing Strategies: A Study on the Craft Beer Sector in Italy
Previous Article in Special Issue
The Effects of Rock Index Tests on Prediction of Tensile Strength of Granitic Samples: A Neuro-Fuzzy Intelligent System
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Stability Analysis of Paste Filling Roof by Cut and Fill Mining

1
Key Laboratory of Deep Coal Resource Mining, Ministry of Education of China, School of Mines, China University of Mining and Technology, Xuzhou 221116, China
2
People’s Government of Dongshan Town, Qilin District, Qujing City 655004, China
*
Authors to whom correspondence should be addressed.
Sustainability 2021, 13(19), 10899; https://doi.org/10.3390/su131910899
Submission received: 24 July 2021 / Revised: 25 September 2021 / Accepted: 26 September 2021 / Published: 30 September 2021
(This article belongs to the Special Issue Advances in Rock Mechanics and Geotechnical Engineering)

Abstract

:
Ensuring the stability of paste false rooves is an important issue in the study of the process of paste filling and slicing mining. Here, a mechanical model of a paste false roof is created to analyze its stability in the process of lower slicing mining in order to determine the minimum slicing thickness of the false roof. We use FLAC3D to simulate and analyze the influence of changes in paste false roof thickness on the stability of the roof. The quantitative functional relationship between the thickness and the subsidence of a false roof, and the optimal thickness of the artificial paste roof, is finally obtained by the development law of the plastic zone in the lower slicing face. The results show that when the thickness of the paste false roof is 3.2 m, the roof can maintain its self-stabilization state and ensure the normal mining of lower layers. Because the same thickness of the upper and lower layers is beneficial for mining replacement and equipment selection in different layered working faces, the optimal thickness of a paste false roof is determined to be 3.2 m.

1. Introduction

Paste filling mining can effectively control surface deformation and overlying rock movement and has become one of the most important technical means to recover coal resources under buildings (structures) [1]. Paste false roofs need to have sufficient strength and a stable bearing capacity. Mastering the instability law of paste filling and layered mining is a prerequisite for safe and efficient mining under pressure [2]. A series of studies have been carried out.
Some researchers established a mechanical model of the basic roof key blocks for paste filling mining and calculated the thickness of the cracked roof [3,4,5]. For example, the mechanical model of the roof rock beam for backfill mining and its differential equations were established. The thin plate theory was used to analyze the stress conditions of the false roof under the mine stope [6,7]. The “embedded beam” mechanical model for the overfilled roof of the roadway was established. The finite-length beam model was used to mechanically analyze the “filling body-direct roof” integral support structure [8,9,10].
At present, the existing research mainly analyzes the stability of the false roof as it relates to the backfill, regarding aspects such as the amount of underfilling, the compression rate of the filling, the filling rate, the mining method, and the overburden load [11,12]. There are few reports on the thickness, strength, and stability of the false roof layer during the layered filling and recovery operation under the paste filling and layered mining [13,14]. Because of the limited space, this study will focus on the influencing factor of the thickness of the paste top layer, and use the method of numerical simulation to study the effect of different thicknesses of the false top layer on its stability [15,16].
Under some mining geological conditions, to effectively control surface deformation, recover the “three-under” compressed coal resources, and digest gangue to protect the ecological environment, paste filling mining is the best solution [17]. At present, the paste filling method has been used to implement paste filling in multiple mines, such as the Taiping Coal Mine of Jining Mining Group, the Daizhuang Coal Mine of Shandong Energy Zi Mining Group, the Xima Coal Mine (Shenyang Province, China)., and the Xiaotun Coal Mine of Jizhong Energy Fengfeng Group. Paste filling mining the full height at one time means extracting the full height of the coal seam at one time and then implementing the filling [18]. The paste filling layer mining method means filling the upper layer first, and the upper layer filling body is used as the false roof for the lower layer. The practice of backfill mining shows that layered mining with paste filling has a better ground deformation control effect than one-time mining with paste filling, and the early strength requirement of the filling is lower than that of the layered mining, so the former has stronger industrial applicability than the latter.
This study takes the paste filling mining project at the E1302 face of the Gaohe Coal Mine of Shanxi Lu’an Group as the background. According to the mining geological conditions of the Gaohe Coal Mine, the paste filling layered coal mining method is selected to recover coal resources. Theoretical analysis combined with numerical simulation is applied to study the relationship between the thickness of the paste false roof and the stability of the false roof.

2. Project Overview

The E1302 paste filling face is located in the East first panel of the Gaohe Mine, Lu’an, with a length of 230 m and an advanced length of 430 m. The main mining seam is the No. 3 coal seam, and the average coal thickness is 6.4 m. The paste filling is used for two-layer down mining. The coal hardness coefficient is 0.7, the average inclination angle of the coal seam is 8°, and the roof of the No. 3 coal seam is mudstone, sandy mudstone, siltstone, and partly sandstone. The floor is black mudstone, sandy mudstone, and dark-gray siltstone. The lithological characteristics of the roof and floor rocks are shown in Table 1.

3. Analysis of the Mechanism of the Instability Thickness of the Paste False Roof

3.1. Calculation of False Roof Load

(1) Direct top load Q1
Q1 = h1L1γ1
Formula:
h1—the thickness of the direct roof, m.
L1—suspended ceiling distance, m.
γ1—Volume force, MN/m3.
If the suspended ceiling distance is equal to the controlled ceiling distance, then:
Q1 = h1γ1 = 0.216 MPa
(2) Basic top load Q2
According to the composite beam theory, the high-strength rock layer and the upper, softer rock layer form a composite beam structure, and the load formed by the upper, softer rock layer and the lower key layer can be solved using composite beam theory.
Considering the influence of n layers of overlying rock on the bottom layer, replace (q1)x with (qn)1, and then:
q n 1 = E 1 h 1 3 γ 1 h 1 + γ 2 h 2 + + γ n h n E 1 h 1 3 + E 2 h 2 3 + + E n h n 3
Using this formula to calculate the load formed on the bottom rock beam, it must be calculated gradually from bottom to top, that is, the sequence is (q1)1, (q2)1, ..., (qn)1; when the calculation is (qi + 1)1< (qi), then the calculation stops at one o’clock. It shows that the i + 1 layer does not generate additional load on the bottom layer. At the time (qn)1 = (qi)1, it can be judged that the n + 1th layer does not affect the load of the first layer, that is, the strength of the nth layer on top of the rock. It is higher and can support the overburdened load without load on the basic roof. After calculating (q6)1 < (q5)1, the sixth layer of siltstone above the basic roof has the characteristics of large thickness and high strength, and the overlying load will not be transmitted downwards. Therefore, the load on the basic roof is Q2= (q5)1 = 0.443 MPa.
(3) Calculation of self-load of paste false top
The load of the paste false top is considered to be safe, and it is calculated according to the larger thickness of 5.2 m:
Q3 = hmaxγ3 = 0.095 MPa
In summary, the false roof bearing load q is calculated as:
q = Q1 + Q2 + Q3 = 0.754 MPa

3.2. Theoretical Analysis of False Roof Instability Thickness

There will be no collapse zone on the roof above the mining face with the paste filling, but only a fracture zone and curved subsidence zone [19]. The practice has shown that the roof control effect of paste filling mining is generally affected by factors such as the amount of roof and floor moving before filling, the amount of underfilling, the compression rate, and the filling rate of the filling body. With the continuous improvement of the filling process and the development of key filling equipment, paste filling technology has become more mature. The amount of top and bottom plate moving before filling, the amount of underfilling, and the compression rate of the filling body can be controlled within a small range, while still meeting the requirements [20,21]. The filling body fills the mined-out area and controls the roof in time, and therefore effectively controls the roof subsidence and the filling requirements of the surface collapse. At present, the commonly used mechanical models of the instability mechanism of the paste false top are the thin “plate” theoretical model and the simply supported “beam” theoretical model [22]. Both theoretical models assume that the paste false top is continuous and homogeneous, and that it was a linear elastic body before yield failure. Since the length of the filling working face in the inclined direction is much longer than the span in the strike direction, the roof in the middle of the working face can be regarded as a beam of countless unit widths inserted into the front and rear rock masses, so it will be filled. The artificial top is considered to be a simply supported beam for research purposes [22]. The load of the simply supported beam and above is transmitted from the beam to the front and rear of the stope before it breaks. At this time, the overburden load is regarded as the uniform load q, the false roof thickness h, and the rock beam length L, as shown in Figure 1. The mechanical model analyzes its force [23].
From elastic mechanics, we know that the bending stress in the beam σx is caused by the bending moment, the extrusion stress σy is caused by the overlying load q, and τxy is caused by the shear stress. Assuming that σy is only a function of y, we have:
σ y = f y
The Airy stress function is:
σ x = 2 Φ y 2 f x x ,     σ y = 2 Φ x 2 f y y ,     τ xy = 2 Φ x y
Bringing the formula into the Airy stress function Formula (7), we obtain:
Φ = x 2 2 f y + xf 1 y + f 2 y
Solve the stress function from the compatibility equation, put the Equation (8) into the compatibility equation, and simplify it to obtain:
4 Φ x 4 + 2 4 Φ x 2 y 2 + 4 Φ x 4 = 0
Omit the primary term and constant term, and further simplify to obtain the stress function:
Φ = x 2 2 A y 3 + By 2 + C y + D + x E y 3 + F y 2 + G y A 10 y 5 B 6 y 4 + H y 3 + K y 2
Solve the stress component from the stress function, and put the Formula (10) into the Airy stress function (7) to obtain the stress component:
σ x = x 2 2 6 A y + 2 B + x 6 E y + 2 F 2 A y 3 2 B y 2 + 6 H y + 2 K
σ y = A y 3 + B y 2 + C y + D
σ x y = x 3 A y 2 + 2 B y + C 3 E y 2 + 2 F y + G
Considering that the span of the beam is much larger than the thickness of the beam, the stress boundary conditions on the main boundary can be satisfied, and the integral stress boundary of the Saint-Venant principle is used to replace the minor part of the boundary so that the boundary conditions are approximately satisfied:
σ y y = h 2 = 0 ,     σ y y = h 2 = q ,     ( τ x y ) y = ± h 2 = 0
Therefore, first consider the main boundary conditions, combine the formula with the formula, and solve:
σ x = 6 q h 3 x 2 y + 4 q h 3 y 3 + 6 H y + 2 K
σ y = 2 q h 3 y 3 + 3 q 2 h y q 2
τ x y = 6 q h 3 x y 2 3 q 2 h x
Next, consider the left and right boundary conditions. Since the model is arranged symmetrically, one side can be considered. On the right side of the beam, there is no horizontal plane force, that is, σx = 0, so:
h 2 h 2 σ x x = L d y = 0
h 2 h 2 σ x x = L y d y = 0
Combining Equations (15)–(17) with conditional Equations (18) and (19), the stress component is solved as:
σ x = 6 q 17 h 3 ( L 2 x 2 ) y + q y h 4 y 2 h 2 3 85
σ y = q 2 L + y h × L 2 y h 2
τ xy = 6 q h 3 x   h 2 4 y 2
It can be seen from the stress expressions (20)–(22) that in a long beam whose length is far greater than its thickness, the bending stress σx is of the same order as q(L/h)2, which is the main stress; the shear stress τxy is the same as q(L/h)2. The same order of q(L/h) is the secondary stress; the extrusion stress σy is the same order of q, and this is secondary stress, which is not considered. Therefore, the failure of the false roof is mainly caused by the bending moment and shearing force.
Breaking of the false roof is caused by the bending moment:
The maximum bending stress on the beam reaches its maximum value at y = h/2, so:
σ max = 6 q L 2 x 2 17 h 2 + q 1 2 h 6 3 q 85 h 2
Omit the high-order infinitesimal terms, then:
σ max = 30 q L 2 x 2 3 q 85 h 2
If the tensile strength of the false top is FT, the critical conditions for the bending instability and breaking of the paste false top are considered as follows:
σ max σ = F T
Then there are:
30 q L 2 x 2 3 q 85 h 2 F T
Breaking of the false roof is caused by shearing force:
If the maximum bending stress on the beam reaches its maximum value on the neutral axis at y = 0, omitting the high-order infinitesimal terms, then:
( τ xy ) max = 3 q 2 h x
If the tensile strength of the false roof is FS, the critical condition of shear instability and fracture of the paste false roof is considered as follows:
( τ xy ) max τ xy = F s
Then there are:
3 q 2 h x F s
According to the formula, the minimum thickness of false roof hmin is obtained because it also satisfies:
h min 30 q L 2 x 2 3 q 85 F T
h min 3 qx 2 F s
This is the case considering that x = 0.5 L is the dangerous section of the false roof’s simply supported beam. According to the paste material mechanics experiment, the 28 d tensile strength of the paste false roof is 0.5 MPa, and the shear strength is 1.5 MPa. The false roof span L is 8 m, and the false roof bears the load of 0.754 MPa, and Formulas (30) and (31) are used to obtain:
h min 30 q L 2 x 2 3 q 85 F T = 447 q 85 F T = 2.81 m
Therefore, the comprehensive analysis shows that the thickness of the false roof layer should be greater than 3.01 m.

4. Numerical Simulation Analysis of the Instability Thickness of the Paste False Roof

4.1. Numerical Simulation Modeling of Layered Mining

The length of the E1302 working face is 230 m, the advance length is 430 m, the average thickness is 6.4 m, and the buried depth of the three coal seams is about 420 m. According to the actual filling and mining process of the working face, the geometric model of the stope structure of paste filling and layered mining is established as shown in Figure 2. By referring to the geometric model and combining the physical and mechanical parameters of the rock formation in Table 2, we can use FLAC3D [24,25] to establish a numerical simulation model as shown in Figure 3. Horizontal displacement constraints are imposed on the north and south boundaries of the model, and the bottom boundary of the model is fixed in the vertical direction. The top rock layer is buried at about 310 m. Therefore, a vertical downward uniform load of 7.75 MPa is applied to the top to simulate the overburden load. To eliminate the influence of the boundary of the model, a boundary coal pillar of 60 m is left on both sides of the working face along the advancing direction. The size of the model is 550 m × 300 m × 160 m in the order of length × width × height. In the process of analysis and calculation, the weight of each rock layer is considered in the vertical direction of the model, and the Mohr–Coulomb strength criterion is used to solve the constitutive relationship [26,27].

4.2. Numerical Simulation Schemes of Different False Roof Thicknesses

According to the actual mining and filling operation of the E1302 working face of Gaohe filling mining, the span of the roof control area of the working face is about 8 m. The upper stratified filling body is considered by the unidirectional force state, and the later strength requirement of the filling body is determined to be 5 MPa according to the Bieniawski paste late strength formula [28,29]. In the simulation, the controlled variable method is adopted to control the layer thickness as a single variable. According to the actual mining and filling operation, the compressive strength of the filling body false top is set to 5 MPa, and the collapse distance of the lower layer working face is 8 m. Other simulation variables are assumed to be consistent. The sum of the height of the false top layer and the lower layer working face of the model is the total thickness of the coal seam, which is 6.4 m. Based on the experience in the field, by changing the thickness of the false roof (2.2 m, 3.2 m, 4.2 m, 5.2 m), the distribution of the false roof plastic zone and the amount of roof subsidence have been analyzed. Previously, scholars have used the method of finite element numerical simulation to analyze the relationship between the displacement of the middle part of the false roof section and the instability and fracture of the false roof under the action of the weight and the upper rock formation according to the simulation results of the displacement field [10]. Some scholars also use FLAC to establish a numerical simulation model and analyze the critical conditions for roof failure and the instability of the working face under different filling steps by observing the plastic zone.
In this paper, the FLAC3D 5.0 numerical simulation method is used to study the roof of the working face through the development of the plastic zone, and the influence of the false roof thickness on the stability of the false roof is studied [30]. In the 3D geometric model, the blue survey lines and pink survey points are arranged at the bottom of the upper paste filling roof, and 19 survey points are set at equal intervals with an interval of 1 m. The survey line layout section is shown in Figure 4.

4.3. Simulation Results and Analysis of the Displacement of the Paste False Top

4.3.1. Correlation Analysis of False Roof Thickness and False Roof Subsidence

According to the different thickness conditions in the scheme, a simulation is carried out, and the displacement cloud diagram of the lower layer working face is shown in Figure 5.
It can be seen from Figure 5 that: when the thickness of the false roof is 2.2 m, the maximum sinking amount of the roof is 81 mm; when the thickness of the false roof is 3.2 m, the maximum sinking amount of the roof is 75 mm; and when the thickness of the false roof is 4.2 m, the sinking amount of the roof is at its maximum. When the false roof is 72 mm and the thickness of the false roof is 5.2 m, the maximum subsidence of the roof is 58 mm. It can be roughly judged that the maximum subsidence of the false roof is negatively correlated with the thickness of the false roof.
One can draw the false roof subsidence curve of each measurement point according to the cloud map displacement. As shown in Figure 6, the false roof subsidence area can be divided into three areas, A, B, and C, according to the distribution law of the subsidence curve; that is, the false roof subsidence boundary zone (zone A), the false roof subsidence transition zone (zone B), and the false roof subsidence central zone (zone C). Select different sections in each area of A, B, and C in Figure 6, and one can draw the false roof in sections. The relationship between the amount of subsidence and the thickness of the false roof is shown in Figure 7.
y = 2.83 x 3 27 x 2 + 87.21 x 170
To further quantitatively analyze the functional relationship between the false roof thickness and the false roof subsidence in the central area of false roof subsidence (zone C), the distance between the working face and the origin of the survey line is x, and when 52 m ≤ x ≤ 57 m, C is obtained the area function relationship is Formula (32), the degree of fit R2 = 0.989, and the fitted curve is shown in Figure 8.
From Figure 8 and Formula (32), it can be seen that in the central area of false roof subsidence (zone C), the thickness of the false roof is negatively correlated with the amount of false roof subsidence, and the amount of subsidence decreases with the increase of false roof thickness. The slope of the curve first gradually decreases and then increases, which means that when the false roof thickness is between 2.2 and 3.2 m, the subsidence speed of the false roof decreases with the increase in thickness. When the false roof thickness is between 3.2 and 4.2 m, the subsidence value of the false roof changes little. The sinking speed of the false roof increases with the increase of thickness after 4.2 m.
Similarly, quantitatively analyze areas B and C, and one can obtain the piecewise function of the correlation between the false roof thickness and the false roof subsidence:
y =   0.35 x 2 + 0.51 x 1.11                               R 2 = 0.999 ,       x 30 , 50 58 , 78 2.33 x 3 + 23 x 2 67.66 x + 20               R 2 = 0.999 ,       x 50 , 52 57 , 58 2.83 x 3 27 x 2 + 87.21 x 170         R 2 = 0.989 ,       x 52 , 57                                    

4.3.2. Analysis of the Development Law of Plastic Zone under Different False Top Thickness

The plastic zone distribution of the lower layer working face with different filling body false roof thickness is shown in Figure 9.
It can be seen from the figure that the failure rule of the plastic zone around the working face is the same. The false roof subsidence center area mainly has a tensile failure, and the subsidence boundary area has a shear failure. In the lower part of the subsidence transition area, the working face bottom plate is subjected to shear stress and tension at the same time. Plastic failure occurs under the influence of stress, and it is in a state of multi-directional stress failure.
When the thickness of the false roof is 2.2 m, tensile failure occurs mainly in the sinking center area of the false roof, while shear failure occurs in the sinking boundary area. When the false roof thickness is 3.2 m, there is a partial tensile failure within 1 m of the subsidence center area, and shear failure occurs in the subsidence boundary area. When the thickness of the false roof is 3.2~4.2 m, the deep part of the plastic zone does not increase significantly, but rather shear failure occurs in the sinking boundary area, and the plastic failure range does not increase significantly with the increase of the false roof thickness in the vertical and horizontal directions. When the thickness of the false roof is 4.2~5.2 m, the distribution area of the plastic zone around the working face is significantly reduced, and only less shear deformation occurs in the subsidence boundary area, and the false roof is in a stable state.
In summary, when the thickness of the false roof is less than 3.2 m, the extent of the plastic failure zone is significantly reduced as the thickness of the false roof increases; when the thickness of the false roof is 3.2–4.2 m, the thickness of the false roof increases, and in the plastic failure zone there is no significant expansion of the range; when the false roof thickness is greater than 4.2 m, the plastic failure zone will be further reduced with the increase of the false roof thickness. It should be pointed out that though the study is based on theoretical and numerical investigation, the results are relatively limited due to the selection of parameters. In the future, a field study will be conducted and then compared to the numerical study, which will be a guideline in this study.

5. Conclusions

(1) Based on the geological conditions of the Gaohe Mine, a mechanical model of the paste false top was established. By analyzing the two limit states of the paste false top in tensile failure and shear failure, it was concluded that the thickness of the upper layer should be at least 3.01 m for stability.
(2) The false roof subsidence area is divided into three areas: false roof subsidence boundary area, false roof subsidence transition area, and false roof subsidence central area, and the false roof thickness in the central area of false roof subsidence is obtained by the fitting curve. The cubic function relationship of the subsidence and the corresponding change trend of the false roof subsidence with the increase of thickness is analyzed through the slope of the curve.
(3) By analyzing the distribution of the plastic zone in the lower layered working face, it can be seen that when the false roof thickness is 2.2–3.2 m, the plastic failure range decreases with the increase of the false roof thickness; when the false roof thickness is 3.2–4.2 m, the plastic failure zone does not expand significantly with the increase of the false roof thickness; when the false roof thickness is 4.2–5.2 m, the plastic failure zone is further reduced with the increase of the false roof thickness. Since the total thickness of the top layer and the bottom layer is 6.4 m, the thickness of the top layer is too large or the thickness of the top layer is too small, which will increase the difficulty of mining, increase the requirements for support and equipment selection, and is not conducive to the work of upper and lower layers. Therefore, it is finally determined that the thickness of the false top layer of the paste filling is 3.2 m.

Author Contributions

Conceptualization by Q.C. and Y.S. (Yifeng Sun); methodology, Q.L.; investigation, Q.C.; writing—original draft preparation, Q.C. and Q.L.; writing—review and editing, Z.L., H.Z., and Y.S. (Yuantian Sun); supervision, Q.C.; funding acquisition, Q.C. and H.Z. All authors have read and agreed to the published version of the manuscript.

Funding

This research was supported by the projects of “the Fundamental Research Funds for the Central Universities (2020ZDPY0221, 2021QN1003)”, and the “National Natural Science Foundation of China (52174130, 52104106, 52174089)”.

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Data Availability Statement

The data presented in this study are available on request from the corresponding author. The data are not publicly available due to privacy.

Acknowledgments

The authors are grateful to Gaohe coal mine.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Zhou, H.; Hou, C.; Sun, X.; Qu, Q.; Chen, D. Solid waste paste filling for none-village-relocation coal mining. J. China Univ. Min. Technol. 2004, 33, 154–158, 177. [Google Scholar]
  2. Chang, Q.; Tang, W.; Xu, Y.; Zhou, H. Research on the width of filling body in gob-side entry retaining with high-water materials. Int. J. Min. Sci. Technol. 2018, 28, 519–524. [Google Scholar] [CrossRef]
  3. Du, X.J.; Feng, G.R.; Qi, T.Y. Roof stability analyses of water-preserved and water-stored coal mining with construction backfill mining. J. China Coal Soc. 2019, 44, 821–830. [Google Scholar]
  4. Xiuzhi, S.H.I.; Yonggang, G.O.U.; Xin, C.; Jia, L.U.O. Determination and application of the backfill thickness in the residual ore recovery of crown and sill pillar. J. Min. Saf. Eng. 2016, 33, 1080. [Google Scholar]
  5. Wu, A.; Zhang, A.; Wang, H.; Gongcheng, L.; Rong, K. The mechanical model for the paste false roof and finite element analysis. J. Min. Saf. Eng. 2017, 3, 98–101. [Google Scholar]
  6. Xikui, S.U.N.; Xiushan, L.I.; Xianyuan, S.H.I. Study on mine strata behavior with full backfill effect of paste backfill in coal mine. Coal Sci. Technol. 2017, 45, 48–53. [Google Scholar]
  7. Tang, W.; Sun, X.; Wang, H.; Shi, X.Y.; Li, X.Y.; Li, Y. Study on effect of overlying strata stability during strip coal pillar excavation with paste backfilling. Coal Sci. Technol. 2017, 45, 109–115. [Google Scholar]
  8. Sun, Y.; Li, G.; Zhang, J.; Xu, J. Failure Mechanisms of Rheological Coal Roadway. Sustainability 2020, 12, 2885. [Google Scholar] [CrossRef] [Green Version]
  9. Basarir, H.; Sun, Y.; Li, G. Gateway stability analysis by global-local modeling approach. Int. J. Rock Mech. Min. Sci. 2019, 113, 31–40. [Google Scholar] [CrossRef]
  10. Sun, Y.; Li, G.; Basarir, H.; Karrech, A.; Azadi, M.R. Laboratory evaluation of shear strength properties for cement-based grouted coal mass. Arab. J. Geosci. 2019, 12, 1–11. [Google Scholar] [CrossRef]
  11. Ma, D.; Kong, S.; Li, Z.; Zhang, Q.; Wang, Z.; Zhou, Z. Effect of wetting-drying cycle on hydraulic and mechanical properties of cemented paste backfill of the recycled solid wastes. Chemosphere 2021, 282, 131163. [Google Scholar] [CrossRef]
  12. Chen, J.; Liu, P.; Zhao, H.; Zhang, C.; Zhang, J. Analytical studying the axial performance of fully encapsulated rock bolts. Eng. Fail. Anal. 2021, 128, 105580. [Google Scholar] [CrossRef]
  13. Sun, Y.; Li, G.; Zhang, J. Investigation on jet grouting support strategy for controlling time-dependent deformation in the roadway. Energy Sci. Eng. 2020, 8, 1–8. [Google Scholar] [CrossRef] [Green Version]
  14. Sun, Y.; Li, G.; Zhang, J.; Qian, D. Experimental and numerical investigation on a novel support system for controlling roadway deformation in underground coal mines. Energy Sci. Eng. 2020, 8, 490–500. [Google Scholar] [CrossRef] [Green Version]
  15. Ma, D.; Duan, H.; Liu, J.; Li, X.; Zhou, Z. The role of gangue on the mitigation of mining-induced hazards and environmental pollution: An experimental investigation. Sci. Total Environ. 2019, 664, 436–448. [Google Scholar] [CrossRef] [PubMed]
  16. Chen, J.; Zhao, H.; He, F.; Zhang, J.; Tao, K. Studying the performance of fully encapsulated rock bolts with modified structural elements. Int. J. Coal Sci. Technol. 2021, 8, 64–76. [Google Scholar] [CrossRef]
  17. Sun, Y.; Li, G.; Zhang, J.; Qian, D. Stability Control for the Rheological Roadway by a Novel High-Efficiency Jet Grouting Technique in Deep Underground Coal Mines. Sustainability 2019, 11, 6494. [Google Scholar] [CrossRef] [Green Version]
  18. Ma, D.; Zhang, J.; Duan, H.; Huang, Y.; Li, M.; Sun, Q.; Zhou, N. Reutilization of gangue wastes in underground backfilling mining: Overburden aquifer protection. Chemosphere 2021, 264, 128400. [Google Scholar] [CrossRef]
  19. Chang, Q.; Sun, X.; Dong, X.; Shao, S. Stability analysis of cemented paste backfill false roof in highwall mining: A case study. Desalin. Water Treat. 2021, 219, 96–102. [Google Scholar] [CrossRef]
  20. Sun, Y.; Zhang, J.; Li, G.; Ma, G.; Huang, Y.; Sun, J.; Wang, Y.; Nener, B. Determination of Young’s modulus of jet grouted coalcretes using an intelligent model. Eng. Geol. 2019, 252, 43–53. [Google Scholar] [CrossRef]
  21. Sun, Y.; Li, G.; Zhang, N.; Chang, Q.; Xu, J.; Zhang, J. Development of ensemble learning models to evaluate the strength of coal-grout materials. Int. J. Min. Sci. Technol. 2021, 31, 153–162. [Google Scholar] [CrossRef]
  22. Lv, W.; Wu, Y.; Ming, L.; Yin, J. Migration law of the roof of a composited backfilling longwall face in a steeply dipping coal seam. Minerals 2019, 9, 188. [Google Scholar] [CrossRef] [Green Version]
  23. Fengshan, M.A.; Feng, L.I.U.; Jie, G.U.O.; Rong, L.U.; Huigao, G.U.O.; Yongyuan, K.O.U. Backfill stability analysis of steep dip mine using filling method. J. Eng. Geol. 2018, 26, 1351–1359. [Google Scholar]
  24. Abbasi, B.; Russell, D.; Taghavi, R. FLAC3D mesh and zone quality. In Proceedings of the FLAC/DEM Symposium, Hangzhou, China, 22–24 October 2013; Itasca Consulting Group: Minneapolis, MN, USA, 2013. [Google Scholar]
  25. Sears, M.M.; Slaker, B.; Rashed, G.; Winfield, J. Numerical model validation and analysis of a dipping limestone pillar using FLAC3D. In Proceedings of the 53rd US Rock Mechanics/Geomechanics Symposium, New York, NY, USA, 23 June 2019; OnePetro: New York, NY, USA, 2019. [Google Scholar]
  26. Sun, Y.; Li, G.; Zhang, J.; Sun, J.; Huang, J.; Taherdangkoo, R. New Insights of Grouting in Coal Mass: From Small-Scale Experiments to Microstructures. Sustainability 2021, 13, 9315. [Google Scholar] [CrossRef]
  27. Sun, Y.; Bi, R.; Chang, Q.; Taherdangkoo, R.; Zhang, J.; Sun, J.; Huang, J.; Li, G. Stability Analysis of Roadway Groups under Multi-Mining Disturbances. Appl. Sci. 2021, 11, 7953. [Google Scholar] [CrossRef]
  28. Wu, K.; Shao, Z.; Qin, S.; Wei, W.; Chu, Z. A critical review on the performance of yielding supports in squeezing tunnels. Tunn. Undergr. Sp. Technol. 2021, 115, 103815. [Google Scholar] [CrossRef]
  29. Wu, K.; Shao, Z.; Sharifzadeh, M.; Hong, S.; Qin, S. Analytical computation of support characteristic curve for circumferential yielding lining in tunnel design. J. Rock Mech. Geotech. Eng. 2021, 13, 1–13. [Google Scholar]
  30. Zhu, D.; Yan, E.; Hu, G.; Lin, Y. Revival deformation mechanism of Hefeng Landslide in the Three Gorges Reservoir based on FLAC3D software. Procedia Eng. 2011, 15, 2847–2851. [Google Scholar] [CrossRef] [Green Version]
Figure 1. The mechanical model of the paste false top.
Figure 1. The mechanical model of the paste false top.
Sustainability 13 10899 g001
Figure 2. Schematic diagram of the geometric model of the stope structure of the paste filling and layered mining.
Figure 2. Schematic diagram of the geometric model of the stope structure of the paste filling and layered mining.
Sustainability 13 10899 g002
Figure 3. Numerical simulation model of paste filling and layered mining.
Figure 3. Numerical simulation model of paste filling and layered mining.
Sustainability 13 10899 g003
Figure 4. The sectional view of the model measuring the point layout.
Figure 4. The sectional view of the model measuring the point layout.
Sustainability 13 10899 g004
Figure 5. Displacement cloud diagrams of different false roof thicknesses. (a) Displacement cloud image of false roof thickness of 2.2 m, (b) Displacement cloud map of false roof thickness of 3.2 m, (c) Displacement cloud image of false roof thickness 4.2 m, (d) Displacement cloud image of false roof thickness of 5.2 m.
Figure 5. Displacement cloud diagrams of different false roof thicknesses. (a) Displacement cloud image of false roof thickness of 2.2 m, (b) Displacement cloud map of false roof thickness of 3.2 m, (c) Displacement cloud image of false roof thickness 4.2 m, (d) Displacement cloud image of false roof thickness of 5.2 m.
Sustainability 13 10899 g005
Figure 6. The amount of false roof subsidence at different positions of the working face.
Figure 6. The amount of false roof subsidence at different positions of the working face.
Sustainability 13 10899 g006
Figure 7. The relationship between false roof thickness and subsidence in different regions positions of the working face.
Figure 7. The relationship between false roof thickness and subsidence in different regions positions of the working face.
Sustainability 13 10899 g007
Figure 8. The fitting curve of false roof subsidence–thickness in area C.
Figure 8. The fitting curve of false roof subsidence–thickness in area C.
Sustainability 13 10899 g008
Figure 9. Distribution of the plastic zone with different false roof thicknesses.
Figure 9. Distribution of the plastic zone with different false roof thicknesses.
Sustainability 13 10899 g009
Table 1. Lithological characteristics of the roof and floor of the No. 3 coal seam.
Table 1. Lithological characteristics of the roof and floor of the No. 3 coal seam.
Top and Bottom NameRock Stone
Name
Platts Hardness
Coefficient
Thickness/m
Old topFine-grained sandstone—Siltstone4.0412.42~ 13.65
Direct topSandy mudstone—Fine siltstone2.746.87~7.85
Mudstone1.221.82~2.13
Direct bottomMudstone—Siltstone2.712.96~3.05
Old bottomSandy mudstone—Fine-grained sandstone3.628.37~8.94
Table 2. Part of the physical and mechanical parameters of the model rock.
Table 2. Part of the physical and mechanical parameters of the model rock.
Rock Formation
Lithology
Basic Parameters
Layered Thick
Degree /m
Bulk Modulus the Amount/GPaShear Modulus/GPaElastic Modulus/GPaPoisson’s RatioCohesion/MPaInternal
Friction Angle/°
Density/kg·m−3Tensile Strength/MPa
Medium-grained sandstone22019.1113.1632.120.222.7225.1326502.7
Siltstone13.113.068.6021.160.232.0226.1026002
Sandy mudstone7.210.486.2915.730.251.1930.6725101.1
Mudstone1.97.655.5013.310.210.8931.0323200.8
3# coal seam6.41.370.671.730.290.5329.2414500.5
Mudstone1.17.655.5013.310.210.8931.0323200.8
Siltstone1.913.068.6021.160.232.0226.1026002
Fine-grained sandstone8.911.647.3218.160.241.3233.2124102.3
5 MPa filling body-0.080.130.20.270.3330.0019000.24
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Chang, Q.; Sun, Y.; Leng, Q.; Liu, Z.; Zhou, H.; Sun, Y. Stability Analysis of Paste Filling Roof by Cut and Fill Mining. Sustainability 2021, 13, 10899. https://doi.org/10.3390/su131910899

AMA Style

Chang Q, Sun Y, Leng Q, Liu Z, Zhou H, Sun Y. Stability Analysis of Paste Filling Roof by Cut and Fill Mining. Sustainability. 2021; 13(19):10899. https://doi.org/10.3390/su131910899

Chicago/Turabian Style

Chang, Qingliang, Yifeng Sun, Qiang Leng, Zexu Liu, Huaqiang Zhou, and Yuantian Sun. 2021. "Stability Analysis of Paste Filling Roof by Cut and Fill Mining" Sustainability 13, no. 19: 10899. https://doi.org/10.3390/su131910899

APA Style

Chang, Q., Sun, Y., Leng, Q., Liu, Z., Zhou, H., & Sun, Y. (2021). Stability Analysis of Paste Filling Roof by Cut and Fill Mining. Sustainability, 13(19), 10899. https://doi.org/10.3390/su131910899

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop