Abstract
A novel and environmental-friendly backfill mining method known as upward slicing longwall-roadway cemented backfill (USLCB) technology has recently been proposed and successfully applied in mines extracting extra-thick coal seams located under sensitive areas. This paper studies the effects USLCB had on roof movement and failure behavior using the mechanical analysis approach. The application of USLCB in the Gonggeyingzi Mine is taken as a case study with roof movement behavior being monitored over a single mining cycle, as well as over multiple mining cycles of different coal slices. In addition, backfill performance requirements to prevent roof failures where USLCB is implemented are studied. The results show that the deflection curves of the roof at the end of each mining cycle during mining the first and the six slices are symmetrical, but they change from asymmetrical to symmetrical during the mining progresses of the second slice to the fifth slice. The final state of roof movement after the first slice, and through to the fifth slice, displays an obvious “flat bottom” pattern in the middle of the deflection curve. The roof movement during the removal of the top slice is noticeably different from other slices. The results also show that the requirements of the elastic modulus, as well as the strength of the backfill, increase as the number of mined slices increases from 1 to 5, but the requirements drop sharply for mining the top slice.
1. Introduction
The process of mining can be defined as economically recovering valuable minerals or other materials from the earth’s crust either at the surface or underground [1, 2]. Thick coal seams whose thickness is over 3.5 m comprise approximately half of the world’s coal reserves and 70–80% of thick coal seams are extracted using underground mining methods [3, 4]. In underground mining, the underground voids created as a result of extracting the valuable materials can present serious environmental challenges and are expected to be managed and remedied by filling them through a process known as backfilling technology [5–7]. Backfilling technology can control ground movement and prevent strata fracture and surface subsidence, as well as serve as an opportunity to dispose waste materials created during the extraction process [8, 9]. An important backfill technology is cemented paste backfill [10–14], which is widely used across the industry for filling excavated stopes with a mixture of mine waste materials such as tailings, slags, and water [15–18]. Although cemented backfill has been widely used in industry because of its advantages, such as high backfilling ratio and flexible distribution systems [19, 20], the technology has been rarely applied in operations that are extracting extra-thick coal seams.
There are currently two main backfill mining methods used for extracting the majority of thick coal seams, including single-pass backfill mining and multislice backfill mining [3, 4, 21, 22]. The major constraint with applying the single-pass backfill mining method in extra-thick coal seams is the limiting factor of the maximum mining height [23, 24]. The multislice backfill mining method can provide a viable alternative to the previously stated methods for mining extra-thick coal seams as the backfill can support the overlying strata and prevent fractures and failures from developing [25, 26].
A novel and environmental-friendly mining and backfilling method known as upward slicing longwall-roadway cemented backfill (USLCB) technology to be proposed by Deng et al. and was initially published in 2017 [14]. The technical principles of the USLCB are shown in Figure 1 and described as follows [14]. The cemented backfill is prepared in a backfill plant at the surface and is transported underground to stopes through a network of backfilling pipelines [7, 27–29]. The extra-thick coal seam is divided into several slices with the deepest slice being mined first and the resulting void being filled with backfill material [30, 31]. Once the backfill material hardens [32, 33], it is used as a floor to facilitate the mining of the subsequent overlying slices as the entire coal seam is mined out progressively from bottom to top. The general longwall mining roadway system, including the preparation of roadways and gateways, is constructed prior to mining each slice. The roadheader is used to excavate crosscuts from the headgate to the tailgate, and after one crosscut is driven, it is filled with cemented backfill while the next crosscut is driven concurrently. A mining cycle contains 7 mining steps, and the mining sequence has been shown in Figure 2. After a mining cycle is finished, the coal pillar between first and second backfilling sections is totally mined out. The driving of crosscuts and the backfilling processes are repeated until the entire slice is mined out and backfilled.


The USLCB technology was applied to extract an extra-thick coal seam with a thickness of 21 meters located under aquifers to protect the water, and it achieved positive application results [14]. However, although USLCB has indicated it can be a successful technology, the roof movement and failure behavior of USLCB need to be studied further. In this paper, a mechanical analysis approach is used to study roof movement behavior in mining extra-thick coal seams using USLCB technology. The application of this technology in the Gonggeyingzi Coal Mine is taken as a case study, with the roof movement behavior monitored and analyzed over a single mining cycle, over multiple mining cycles, and within each pass to excavate the different slices. In addition, the performance requirements of backfill material to prevent roof failures during the mining process of USLCB are calculated.
2. The Roof Movement of Each Slice in the Mining Process of USLCB
A mechanical analysis approach [34–36] is used to study the roof movement characteristics during the mining of each slice when extracting an extra-thick coal seam using USLCB technology. Nonhomogeneous elastic foundation beam models are developed in which the roof during the USLCB mining process is simplified as a beam. The mining process of USLCB is also considered in the mechanical models. It should be noted that the term roof is used to refer to different things. It refers to the stratum above the coal seam during the mining of the top slice. However, during the mining of the bottom slices, the roof refers to the remaining coal located above the extracted slice.
The coal body, backfilling body, and the combination of coal and backfilling are sectionally simplified as nonhomogeneous elastic foundations in the mechanical models [37–39]. According to the specific conditions of mining the extra-thick coal seam using USLCB technology, two mechanical models are established: a mechanical model for mining the first slice, located at the bottom, and another mechanical model for mining the remaining overlying slices.
2.1. The Roof Movement Behavior of the First Slice
2.1.1. The Mechanical Model of Mining the First Slice
The process of USLCB technology is that mined-out space is immediately filled with backfill material so that overlying strata are supported by the backfilling body and the coal pillars. The foundation structure of mining the first slice is shown in Figure 2. Therein, Ls represents the total length of the mining area; L is the length of each mining cycle; n is the number of completed mining cycles; nL is the length of mined range; ns represents the individual steps within the mining cycle; a represents the width of each excavation step; l1 and l3 denote the width of first and second backfilling sections in a mining cycle, respectively; l2 is the width of the coal pillar between the first and second backfilling sections in a mining cycle; and l4 is the width of the coal pillar after the second backfilling section in a mining cycle.
During the mining cycling, the values of l1, l2, l3, and l4 are constantly changing as the mining and backfilling operations undertaken; however, the sum of these lengths is fixed, as shown in the following equation:
There are 7 excavation steps in each mining cycle, so the value of ns is from 1 to 7. The sizes of the crosscuts, pillars, and backfilling areas represented by the values l1, l2, l3, and l4 change as each step of the mining cycle is completed, and how these values change with each step of the mining cycle is summarized in Table 1.
The mechanical model of mining the first slice is shown in Figure 3, in which the elastic foundation beam is divided into 5 sections. Therein, L0 and L4 represent the range of the stress-concentrated area in the boundary pillars; L1 is the range from the left mining limit boundary to the first coal pillar which means the width of past mining cycles plus the width of the first backfilling section in current mining cycle; L2 is the width of the coal pillar between the first and second backfilling sections in present mining cycle; and L3 represents the width of the second backfilling section in present mining cycle which begins after the first coal pillar and ends at the second coal pillar.

The parametric correspondence between the mechanical model of USLCB and the mining process of USLCB is shown in the following equation:
2.1.2. The Mechanical Analysis of Roof Deflection during the Mining of the First Slice
Taking the point O of the beam as the origin, a coordinate system is established with the X-axis as the direction of mining and the Y-axis as the vertically downward direction. These axes are shown in Figure 3. Therein, and are the elastic foundation coefficients of coal and the backfilling body, respectively. The deflection calculation process of the beam is shown as follows:(1)Section L0: When , the differential equation for deflection of the roof is shown in the following equation: The load on the coal pillar can be expressed as follows: where is the dimensionless stress concentration factor. The deflection of the roof section L0 is shown in the following equation: where is the characteristic coefficient, in which is the elastic foundation coefficient of coal body, N/m3; E is the elastic modulus of the roof, Pa; and I is the moment of inertia of the roof section, m4. As “x” tends towards negative infinity, the deflection of the roof is constant. The roof can be considered as a semi-infinite body, and this allows equation (5) to be simplified as follows:(2)Section L1: When , the differential equation for deflection of the roof is shown in the following equation: The deflection of the roof in section L1 can be obtained and shown in the following equation: where is the characteristic coefficient, in which is the elastic foundation coefficient of backfilling body, N/m3.(3)Section L2: When , the differential equation for deflection of the roof is shown in the following equation: where is the dimensionless stress concentration factor. The deflection of the roof section L2 can be obtained and shown in the following equation: where is the elastic foundation coefficient of coal body, N/m3.(4)Section L3: When , the differential equation for deflection of the roof is shown in the following equation: The deflection of the roof in section L3 can be expressed as follows: where is the elastic foundation coefficient of backfilling body, N/m3.(5)Section L4: When , the differential equation for deflection of the roof is shown in the following equation: The load on the coal pillar can be expressed as follows: where is the dimensionless stress concentration factor. The deflection of the roof section L4 is shown in the following equation: where is the elastic foundation coefficient of coal body, N/m3. As “x” tends towards positive infinity, the deflection of the roof is constant. The roof can be considered as a semi-infinite body, and equation (15) can be simplified as follows: The rotation , bending moment , and shearing force of an arbitrary cross section of the beam can be calculated by using the following equation: The boundary conditions and the continuity conditions of the mechanical model are shown in equation (18). All the parameters of the model can be solved according to the following conditions:
2.2. Roof Movement Behavior of the Remaining Upper Slices
2.2.1. The Mechanical Model of Mining the Remaining Upper Slices
When mining the first slice using the USLCB method, the floor is composed of the undisturbed original ground. However, after the first slice is mined, and subsequently backfilled, the remaining upper slices will have the backfilling body serving as their floor for mining activities. Therefore, the foundation structure for mining these remaining slices will be different from the first slice, as shown in Figure 4.

The mechanical model of mining the remaining upper slices is shown in Figure 5, in which the elastic foundation beam is divided into 6 sections. Therein, L5 and L10 represent the range of the stress-concentrated area in the boundary pillars; L6 is the range that is a combination of past mining cycles and the first backfilling section in current mining cycle; L7 is the range of the coal pillar between the first and second backfilling sections in current mining cycle; L8 represents the range of the second backfilling section in current mining cycle; and L9 is the range of the coal pillar between the second backfilling sections in current mining cycle and the mining boundary of the previous slice.

The parametric correspondence between the mechanical model and applying USLCB to mining the remaining slices is shown in the following equation:
2.2.2. Mechanical Model for Calculating Deflection When Mining the Remaining Upper Slices
The mechanical model to determine deflection when mining the remaining upper slices that follow the mining of the first slice is similar to how deflection was calculated when mining the first slice, and the deflection calculation process of the beam is shown as follows:(1)Section L5: When , the differential equation for deflection of the roof is shown in the following equation: The load on the coal pillar can be expressed as follows: where is the dimensionless stress concentration factor. The deflection of the roof in section L5 is shown in the following equation: where is the characteristic coefficient, in which (i = 2, …, 6) is the elastic foundation coefficient of coal body, N/m3; E is the elastic modulus of the roof, Pa; and I is the moment of inertia of the roof section, m4. As “x” tends towards negative infinity, the deflection of the roof is constant. The roof can be considered as a semi-infinite body, so and , and equation (22) can be simplified as follows:(2)Section L6: When , the differential equation for deflection of the roof is shown in the following equation: The deflection of the roof in section L6 can be expressed as follows: where is the characteristic coefficient, in which (i = 2, …, 6) is the elastic foundation coefficient of the backfilling body, N/m3.(3)Section L7: When , the differential equation for deflection of the roof is shown in the following equation: The deflection of the roof in section L7 can be obtained and expressed as follows: where is the characteristic coefficient, in which (i = 2, …, 6) is the elastic foundation coefficient for the combination of in situ coal and the backfilling body, N/m3.(4)Section L8: When , the differential equation for deflection of the roof is shown in the following equation: The deflection of the roof in section L8 can be expressed as follows: where is the characteristic coefficient, in which (i = 2, …, 6) is the elastic foundation coefficient of the backfilling body, N/m3.(5)Section L9: When , the differential equation for deflection of the roof is shown in the following equation: The deflection of the roof section L9 can be obtained and shown in the following equation: where is the characteristic coefficient, in which (i = 2, …, 6) is the elastic foundation coefficient of the combination of backfilling and coal body, N/m3.(6)Section L10: When , the differential equation for deflection of the roof is shown in the following equation: The load on the coal pillar can be expressed as follows: where is the dimensionless stress concentration factor. The deflection of the roof section L10 is shown in the following equation: where is the characteristic coefficient, in which (i = 2, …, 6) is the elastic foundation coefficient of the coal body, N/m3. As “x” tends towards positive infinity, the deflection of the roof is constant. The roof can be considered as a semi-infinite body, so and , and equation (34) can be simplified as follows: The boundary conditions and the continuity conditions of the mechanical model are shown in equation (36). All the parameters of the model can be solved according to the following conditions and equation (17):
2.3. Roof Failure Criteria When Mining with USLCB Technology
The first strength theory, also known as the maximum tensile stress theory, is chosen as the failure criteria of the roof stratum. From the elastic foundation theory, the maximum tensile stress of the roof can be obtained from the maximum bending moment Mmax, which can be determined from the bending moment function M(x), as shown in the following equation:where h is the thickness of the roof, m.
Therefore, the failure criteria of the roof stratum can be expressed as follows:where is the maximum allowable tensile stress of the roof, MPa.
3. Results and Discussion: A Case Study
The USLCB mining technology was applied at the Gonggeyingzi Coal Mine to exploit a 21 m thick coal seam that is located 52.3 m under multiple aquifers. Specific details about the location and geological conditions of this mine, as well as the mining technology implemented there, have been published in a previous paper entitled “Upward slicing longwall-roadway cemented backfilling technology for mining an extra-thick coal seam located under aquifers: a case study” [14]. The entire 21 m coal seam in this mine is divided equally into 6 slices, where the thickness of each slice is 3.5 m. The width of each excavation step is 5 m, the vertical in situ stress q0 is 2.5 MPa, and the elastic modulus of the backfilling body is 10 MPa. To study the roof movement of mining different slices of the extra-thick coal seam with USLCB technology, the specific parameters of each unique mining condition in the mine are substituted into the deflection equations seen throughout Section 2, and the results are shown as follows.
3.1. Roof Movement Behavior inside a Single Mining Cycle
The process of extracting the first slice during the first mining cycle is used to explain the changes and development of roof movement and deformation behavior inside a single mining cycle. The deflection curves of the roof after each mining step during this process are shown in Figure 6.

(a)

(b)
Figure 6 shows that, over a single mining cycle, the deflection curves of the roof can be divided into three types, i.e., the small range unimodal deflection curve of step 1-1-1, the bimodal deflection curves of steps 1-1-2 through to 1-1-6, and the large range unimodal deflection curve of step 1-1-7.(1)The small range unimodal deflection curve occurs during the first step (step 1-1-1) of the mining cycle and signifies a mining range that is a single backfilling crosscut. The maximum subsidence of the roof is 11.02 mm, and it appears in the middle of the deflection curve.(2)Bimodal deflection curves occur through the second to the sixth mining steps (step 1-1-2 through to step 1-1-6) when there is a coal pillar in the middle of the mining range, which separates the active mining areas. When the mining range is symmetrical, the deflection curve is also symmetrically distributed. The two peaks appear above the two active parts of the mining range, respectively. The distance between two peaks varies between 25 and 30 m, and the maximum subsidence of the roof increases from 11.26 mm to 13.75 mm gradually.(3)The large range unimodal deflection curve occurs during the last mining step (step 1-1-7) when the mining range is seven backfilling crosscuts and coal pillar in the middle has been removed. Although the shape of the deflection curve is similar to the curve during the first mining step, the subsidence range of the roof is much larger compared to the first mining step. The maximum subsidence of the roof appears in the middle of the deflection curve and it significantly increases to 25.22 mm.
3.2. Roof Movement Behavior during Multiple Mining Cycles
The first four mining cycles of each slice are used as examples to display different roof movement and deformation behavior during multiple mining cycles. The roof deflection curves of each slice at the end of the first four mining cycles are calculated and shown in Figure 7.

(a)

(b)

(c)

(d)

(e)

(f)
Figure 7(a) shows that all the roof deflection curves at the end of each mining cycle for the first slice all have a symmetrical distribution. As the number of mining cycles increases, both the range of roof subsidence and the maximum subsidence increase gradually. The maximum subsidence appears in the middle of the deflection curves. From the first to the fourth mining cycle of the first slice, the maximum subsidence of the roof is 25.22 mm, 45.48 mm, 52.86 mm, and 53.00 mm, respectively. After the third mining cycle of each slice, the maximum subsidence of the roof is almost achieved. Further subsidence after the third slice, resulting from an increase in the mining range, is relatively small, and only the range experiencing the near-maximum subsidence increases. After the fourth mining cycle, an obvious “flat bottom” pattern appeared in the middle of the deflection curve.
Figures 7(b)–7(e) show that the subsidence range of the roof does not significantly change during the mining of the second slice to the fifth slice, but the shape of the deflection curve changes significantly relative to the first slice. During the mining of the second slice to the fifth slice, the shape of the roof deflection curve changes from asymmetry (the end of the first three cycles) to symmetry (the end of the fourth cycle), and the maximum subsidence of the roof increases with its position moving gradually from the left side to the middle. Similar to the first slice, after the third mining cycle, there is a clear “flat bottom” pattern in the middle of the deflection curve of the roof, and the maximum subsidence is nearly achieved as there are no longer significant increases in subsidence with increases in the mining range.
Figure 7(f) shows that the deflection curve of the roof during the mining process of the sixth slice is obviously different from other slices. The main reason is that the roof for this final slice is the sandstone stratum overlying the entire extra-thick coal seam, whereas during the mining of previous slices the roof was overlying slices of coal that would eventually be extracted. After the fourth mining cycle is completed, the “flat bottom” pattern in the middle of the roof deflection curve is not particularly obvious. From the first to the fourth mining cycle, the maximum subsidence of the roof occurring at the end of each cycle is 433.51 mm, 480.56 mm, 510.58 mm, and 514.57 mm, respectively.
3.3. Comparing Roof Movement Behavior during the Mining of Different Slices
In order to compare the roof movement behaviors seen when mining the different slices, the roof deflection curve during the sixth step in the first mining cycle is taken as a baseline example. It is compared with the roof deflection curves seen in the sixth step of subsequent slices. The results are plotted in Figure 8.

Figure 8 shows that, during the sixth step of the first mining cycle, the subsidence range of the first slice is much smaller than that of subsequent slices. An explanation is that mining the first slice is the initial disturbance to the stope, and mining subsequent slices are simply modifying an area that has already been disturbed with the roof being subjected to secondary and even numerous mining operations. The deflection curves of the roof in this sixth step of the first to fifth slices are bimodal curves, while it is a unimodal curve in the sixth slice. The reason for the occurrence of the bimodal curve is that there is a coal pillar in the middle of the mining range which can support the roof, but this supporting effect is relatively weaker during the mining of upper slices.
To compare the roof movement behavior of different slices during the mining process, the deflection curves of the roof after each mining cycle are calculated and plotted as shown in Figure 9.

(a)

(b)

(c)

(d)
Figure 9 shows that the deflection curve of the roof fluctuates less and is more symmetrical at the end of each mining cycle of the first and sixth slices, and the shape of the curve is noticeably different from other slices. The shape of the roof deflection curves is similar at the ends of each mining cycle of the second, third, fourth, and fifth slices, and as more slices are mined, increased roof subsidence is seen.
3.4. The Requirement of Backfill Performance to Prevent Roof Failure
The failure criteria of the roof are presented in Section 2.3. The allowable tensile stress and bending moment of the roof are different for each slice because of the different roof conditions in the mining process, as shown in Table 2.
The elastic moduli for different types of backfilling materials are brought into the mechanical model to calculate the maximum bending moment of the roof during the mining process of each slice. Nonlinear regression is used to obtain a quantitative relationship between the elastic modulus of the backfilling material and the maximum bending moment of the roof for mining different slices. Combined with the allowable bending moment of the roof for each slice, the requirement for the elastic modulus of the backfilling material to ensure the roof does not break can be obtained. Then, the requirements for the strength of the backfilling material can be determined by measuring roof deflection and seeing when it reaches the critical failure state. The results of determining the performance requirement of cemented backfill are shown in Figure 10.

(a)

(b)
Figure 10(a) shows that the quantitative relationship between the maximum bending moment of the roof and the elastic modulus of the backfilling material can be described effectively by the power function. As the elastic modulus of the backfilling material increases, the maximum bending moment of the roof shows an exponential decay. Due to the different conditions of the roof, the maximum bending moment of the roof during the sixth slice is obviously larger than that of the other slices with the same backfilling material.
Figure 10(b) shows that the elastic modulus requirement for the backfill material increases from 2.84 MPa to 29.44 MPa as the number of slices mined increases from 1 to 5, but the requirement drops to 1.86 MPa after the sixth slice is mined. The compressive strength requirement of the backfill varies slightly within the range of 0.148 to 0.153 MPa when the first three slices are mined. However, the strength requirement rapidly increases to 0.219 MPa when the fifth slice is mined and decreases to 0.202 MPa when the sixth slice is mined. The significant change in backfill requirement for the sixth slice shown in the calculation results is due to the different roof conditions compared to other slices. In this case, of all the slices, it is the fifth slice that requires the strongest backfill material for mining the entire extra-thick coal seam with USLCB technology. The performance requirement of cemented backfill calculated by this model is the final condition, which means that the final strength and elastic modulus of the cemented backfill should meet this requirement theoretically. Taking into account the safety margins, the performance of the cemented backfill with a curing time of 7 days can be selected to meet this requirement for the industry application.
4. Conclusions
In this paper, a mechanical analysis approach is used to study roof movement and roof control in mining extra-thick coal seam using USLCB technology. The main findings and conclusions are as follows:(1)In a single mining cycle of the first slice, the deflection curves of the roof can be divided into three types, i.e., the small range unimodal deflection curve in the initial step, the bimodal deflection curves of the middle steps with two active areas, and the large range unimodal deflection curve in the final step when the support pillar is removed. The shape of the curve is significantly affected by the mined area.(2)All of the roof deflection curves that occur at the end of each mining cycle during the extraction of the first slice are symmetrical in their distribution. Also, there is an obvious “flat bottom” pattern occurring in the middle of the deflection curve after the fourth mining cycle is completed.(3)During the mining of the second slice to the end of the fifth slice, the shape of the roof deflection curve changes from asymmetrical to symmetrical. After the third mining cycle is completed during the mining of these slices, there is a clear “flat bottom” pattern in the middle of the roof deflection curve. The roof movement during the mining of the top slice is noticeably different from roof movement during the mining of other slices. This difference is due to that the roof lithology of the top slice is different from previous slices whose immediate roof is actually the remaining coal seam.(4)The relationship between the maximum bending moment of the roof and the elastic modulus of the backfilling material can be well described by the power function. The requirements of the backfill material’s elastic modulus and strength increase as the number of mined slices increases from 1 to 5, but the requirements drop sharply when the sixth slice is mined.
Data Availability
The data used to support the findings of this study are included within the article.
Conflicts of Interest
The authors declare that they have no conflicts of interest.
Acknowledgments
This work was supported by the National Key R&D Program of China (Grant no. 2018YFC0604701), National Natural Science Foundation of China (Youth Program) (Grant no. 51804308), the Research Fund of the State Key Laboratory of Coal Resources and Safe Mining, CUMT (Grant no. SKLCRSM18KF022), and the Natural Science Foundation of Heilongjiang Province of China (Grant no. LH2019E118). The authors would like to thank those researchers who have made great contributions to the project. Special thanks go to the Gonggeyingzi Coal Mine for providing everything for the field application studied in this research.