4.1 Model construction and simulation
We performed a numerical simulation using the FLAC3D software to study the surrounding rock deformation and migration laws in the working face of the longwall-roadway CPB mining in the Nantun Coal Mine. The Mohr-Coulomb model was used for characterizing the failure of the coal-rock masses. According to the actual geological conditions of Nantun Coal Mine, the numerical model illustrated in Fig. 10 had a height of 273. 2 m, a strike length of 300 m, and a dip length of 172 m. First, the numerical model was balanced, with the horizontal displacements of the four vertical planes restrained in the normal direction. The vertical displacement of the model's bottom was set to 0. The ela
Table 1 Physical and mechanical characteristic parameters of the roof, floor, and coal seam
Mechanical
properties
Lithology
|
Thickness
/m
|
Bulk modulus
/GPa
|
Shear modulus
/GPa
|
Density
/t·m-3
|
Tensile strength
/MPa
|
Cohesive force
/MPa
|
Internal friction angle
/(°)
|
Mudstone
|
16
|
7.8
|
2.6
|
2.2
|
3
|
4
|
26
|
Coarse sandstone
|
21
|
17.5
|
8.6
|
2.6
|
3.2
|
3.0
|
26
|
Sandy mudstone
|
30
|
12.8
|
4.6
|
2.5
|
3.2
|
2.4
|
27
|
Coarse sandstone
|
28
|
17.5
|
8.6
|
2.6
|
3.2
|
3.0
|
26
|
Medium sandstone
|
25
|
23.3
|
13.2
|
2.6
|
3.1
|
8.3
|
26
|
Fine conglomerate
|
25
|
18.5
|
12.6
|
2.6
|
2.5
|
2.26
|
26
|
Medium sandstone
|
28.8
|
23.3
|
13.2
|
2.6
|
3.1
|
8.3
|
26
|
Fine conglomerate
|
20
|
20.5
|
12.6
|
2.6
|
2.5
|
2.26
|
35
|
Siltstone
|
7.2
|
20.5
|
6.2
|
2.6
|
1.1
|
1.26
|
30
|
Fine sandstone
|
7.3
|
18.5
|
12.6
|
2.6
|
2.5
|
2.26
|
26
|
Mudstone
|
9.8
|
17.8
|
4.6
|
2.7
|
3
|
4
|
33
|
Medium sandstone
|
10.8
|
23.3
|
13.2
|
2.6
|
3.1
|
8.3
|
26
|
Fine conglomerate
|
3.05
|
18.5
|
12.6
|
2.6
|
2.5
|
2.26
|
26
|
Medium sandstone
|
4.66
|
23.3
|
13.2
|
2.6
|
3.1
|
8.3
|
26
|
Siltstone
|
4.05
|
20.5
|
1.2
|
2.6
|
1.1
|
1.26
|
26
|
3upper coal seam
|
4.73
|
8.6
|
0.9
|
1.38
|
1.8
|
1.96
|
26
|
Interbedding of siltstone and fine sandstone
|
2.2
|
18.5
|
12.6
|
2.6
|
2.5
|
2.26
|
26
|
3lower coal seam
|
3.0
|
8.6
|
0.9
|
1.38
|
1.8
|
1.96
|
26
|
Fine sandstone
|
5.41
|
18.5
|
12.6
|
2.6
|
2.5
|
2.26
|
26
|
Siltstone
|
6.15
|
20.5
|
1.2
|
2.6
|
1.1
|
1.26
|
26
|
Mudstone
|
6.05
|
17.8
|
4.6
|
2.7
|
3
|
4
|
33
|
The numerical simulation of the Nantun Coal Mine was divided into four stages: (Ⅰ) calculate the initial stress state caused by the gravity of the overlying strata, (Ⅱ) excavate the 3upper coal seam, (Ⅲ) fill the 3upper coal seam, (Ⅳ) excavate the 3lower coal seam, (Ⅴ) fill the 3lower coal seam. We were particularly concerned with the changes in the stress state of the surrounding rocks and the changes in the displacement state of the overlying strata in the working face before and after the filling of the 3upper and 3lower coal seams.
4.2 Numerical simulation results
The 3upper and 3lower coal seams were excavated by longwall-roadway CPB mining. A numerical simulation was performed to analyze the variation of overlying strata displacement and surrounding rock stress distribution before and after filling the 3lower and 3upper coal seams with CPB.
4.2.1 Comparative analysis of surrounding rock displacement
The variations in overlying strata displacement and ground surface displacement were obtained by numerical simulation under the backfilling mining of the 3upper and 3lower coal seams of the Nantun Coal Mine. Figure 11 shows the nephograms of overlying strata displacement at the end of the excavation of the 3upper coal seam, filling of the 3upper coal seam, excavation of the 3lower coal seam, and filling of the 3lower coal seam. The displacement changes in the nephograms were analyzed for the 3upper and 3lower coal seams of Nantun Coal Mine along the profile at different stages of mining. As shown in Fig. 11, the monitoring lines were arranged in the immediate roof and ground surface after excavating and filling the 3upper and 3lower coal seams. The displacement was monitored throughout the entire excavation and filling process. The control effect of the CPB on overlying strata and ground surface movements were analyzed.
As shown in Fig. 12 (a), when the excavation of the 3upper and 3lower coal seams was completed, the maximum displacement of the immediate roof was 3.02 m and 3.49 m, respectively. When the filling of the goaf in the 3upper and 3lower coal seams was completed, the maximum subsidence of the immediate roof was 0.14 m and 0.24 m, respectively. They were reduced by 95.3% and 93.12% compared with those under the caving mining, respectively. The variation of the vertical displacement of the immediate roof under caving mining could be described using a U-shaped curve, with smaller values at the two ends and larger values in the middle. The maximum displacement occurred in the center of the goaf. As shown in Fig. 12 (b), the ground surface subsidence was generally small. The maximum ground surface subsidence was 1.62 m and 0.08 m after the excavation and filling of the 3upper coal seam, respectively. After the backfilling mining, the ground surface subsidence decreased by 95.1%. The maximum ground surface subsidence was 240 mm and 120 mm after the excavation and filling of the 3lower coal seam, respectively, both being lower than 300 mm. After the backfilling mining, the ground surface subsidence decreased by 95%. The nephograms of vertical displacement of the surrounding rocks are shown in Fig. 11. The variation of displacement of the overlying strata in the working face was approximately symmetrical about the midline of the advance distance. The vertical displacement of the overlying strata was the largest when approaching the midpoint of the advance distance of the working face. The vertical displacement of the overlying strata decreased significantly after goaf filling.
CPB was used to fill the space taken up by the excavated coal. Besides, CPB itself has a specific strength and can support the roof load upon roof weighting induced by backfilling mining. Hence, the roof integrity is preserved, and the overlying strata damage is put under proper control. As a result, roof caving was lessened, which prevented the upward propagation of failure in the rock strata, thereby reducing ground surface deformation. Given the above, it is feasible to control ground surface deformation by backfilling mining.
4.2.2 Comparative analysis of surrounding rock stress
We compared the nephograms of overlying strata stress distribution under the excavation and filling of the 3upper and 3lower coal seams to investigate the weighting laws in the backfilling mining stope, as shown in Fig. 13. Measuring lines were arranged in the immediate roof and the main roof. The vertical stress curves measured along the horizontal distance are shown in Fig. 14.
The surrounding rock stress was redistributed in the stope, as shown in Fig. 13. There were stress concentration regions in the coal pillars in the cut hole and on the two sides of the stop line. Due to the mining-induced influence, the roof above the working face was considerably unloaded. This area remained within the stress relief zone. Tensile stress occurred in the middle of the goaf. The degree of stress concentration in the coal pillars decreased after the backfill mining of the 3upper coal seam was completed. The peak stress in the immediate roof during caving mining was 15.8 MPa, and it decreased to 6.72 MPa after the backfilling. The point with the peak stress made a rightward shift towards the upper margin of the coal pillar. As shown in Fig. 14, the peak stress in the immediate roof was 17.3 MPa, and 7.34 MPa after the caving mining and the backfilling mining in the 3lower coal seam, respectively. The peak stress in the immediate roof was 15.0 MPa and 5.72 MPa, respectively, indicating that the stress in the immediate roof decreased significantly under backfilling mining than under caving mining. The stress curve of the immediate roof was symmetrical about the midline of the advanced distance. Under caving mining, the stress above the coal pillars on the two sides of the goaf first increased, then decreased, and then increased again, showing minor fluctuations. The stress in the roof above the goaf remained approximately unchanged, consistently being tensile.
Finally, a stress shell was formed in the coal pillars outside the cut hole and the stop line. The stress in the goaf roof finally evolved into tensile stress. The roof was not broken or caved in due to the support of CPB under backfilling mining. The overlying strata tended to move altogether, with a small overall displacement. The rock bursts were inconspicuous in the stope. Therefore, CPB mining could effectively control the stress in overlying strata, reducing the disturbance caused by slice mining of the thick coal seam.
4.3 Prediction of overlying strata deformation
4.3.1 Model for predicting overlying strata deformation
(1) Excavation plan for the coal seam
The ground surface subsidence prediction for the 3upper and 3lower coal seams consisted of the prediction of four cumulative ground surface subsidences: ground surface subsidence upon excavation of the 3upper coal seam, cumulative ground surface subsidence upon backfilling of the 3upper coal seam, cumulative ground surface subsidence upon excavation of the 3lower coal seam, and cumulative ground surface subsidence upon backfilling of the 3lower coal seam. The basic parameters used for predicting ground surface subsidence induced by the excavation of the 3lower and the 3upper coal seams are tabulated in Table 2. Table 3.
Table 2 Basic parameters for ground surface subsidence prediction
Name of the working face
|
Angle of strike
|
Mining thickness
(mm)
|
Coal seam dip angle (°)
|
Subsidence coefficient
|
Tangent of the main influence angle
|
Horizontal movement coefficient
|
Propagation angle (°)
|
Ground surface subsidence induced by caving mining of the 3lower coal seam
|
C13 upper 01/02
|
10
|
4730
|
5
|
0.5
|
1.4
|
0.27
|
88
|
C13 lower 01/02
|
10
|
3000
|
5
|
0.5
|
1.4
|
0.27
|
88
|
Ground surface subsidence induced by the excavation of the 3lower coal seam
|
C13 upper 01/02
|
10
|
236.5
|
5
|
0.5
|
1.4
|
0.27
|
88
|
C13 lower 01/02
|
10
|
150
|
5
|
0.5
|
1.4
|
0.27
|
88
|
Table 3 Coordinate parameters of working faces
1X
|
1Y
|
2X
|
2Y
|
3X
|
3Y
|
4X
|
4Y
|
5X
|
5Y
|
6X
|
6Y
|
H
|
C13upper (lower) 01 working face
|
3916387.6543
|
39532404.9715
|
3916387.3679
|
39532524.9715
|
3917012.2816
|
39532524.9715
|
3917012.2816
|
39532805.5888
|
3917120.6483
|
39532798.8176
|
3917122.8205
|
39532404.9715
|
350
|
C13upper (lower) 02 working face
|
3916387.8696
|
39532284.9715
|
3916387.6292
|
39532400.1427
|
3917122.8206
|
39532400.1427
|
3917122.8206
|
39532284.9715
|
-
|
-
|
-
|
-
|
350
|
4.3.2 Comparison of predicted ground surface deformations
(1) Excavation of 3upper coal seam
According to the prediction methodology, the subsidence was first predicted for the excavation and backfilling of the 3upper coal seam. The basic parameters in the different backfilling working faces shown in Tables 2 and 3 provided important data for predicting the ground surface subsidence. Fig. 15 and 16 show the ground surface subsidence curve vs. position of the working face and the 3D diagrams of mining-induced ground surface subsidence under the excavation and backfilling of the 3upper coal seam using the two mining methods, respectively.
As shown in Fig. 15 and 16, the red box indicates the scope of the backfilling working face. The prediction of ground surface subsidence directly above the backfilling working face was the priority. The ground surface subsidence curves in Fig. 15 show that the maximum ground surface subsidence was 1620 m when the caving mining of the 3upper coal seam was completed. The maximum ground surface subsidence was only 80 mm when backfilling mining was applied to the excavation of the 3upper coal seam. CPB mining had a good control effect for ground surface deformation under the backfilling mining of the working face. As shown by the 3D diagrams of ground surface subsidence under the excavation of the 3upper coal seam in Fig. 16, the ground surface subsidence was the largest around the corner of the CT3 upper (lower) 01 backfilling working face. This was primarily because the mining area was larger around the corners of the CT3 upper (lower) 01 backfilling working face and the CT3 upper (lower) 02 backfilling working face. As a result, the ground surface subsidence was substantial.
(2) Excavation of the 3upper coal seam and 3lower coal seam
We then studied the influence of CPB mining on overlying strata movement under multi-slice mining of the thick coal seam. Following the mining of the 3upper coal seam, we investigated the influence of subsequent CPB mining of the 3lower coal seam on the cumulative ground surface subsidence. Fig. 17 and 18 show the ground surface subsidence curves and the 3D diagrams of ground surface subsidence under the caving mining and CPB mining of the 3upper and 3lower coal seams.
As shown in Fig. 17, the cumulative ground surface subsidence under the caving mining of the 3upper and 3lower coal seams reached 2400 mm. By contrast, the cumulative ground surface subsidence reached 120 mm when the 3upper coal seam was first excavated by CPB mining, followed by CPB mining of the 3lower coal seam. The surface subsidence induced by CPB mining was 1/20 that of the cumulative ground surface subsidence caused by caving mining, showing a 20-fold reduction. CPB mining could effectively control surface subsidence caused by multi-slice mining of the thick coal seam, offering protection for the buildings above the ground.
As shown in Fig. 18, caving mining of two-slice coal and CPB mining had a similar influence on ground surface subsidence. The ground surface subsidence was the largest around the corner of the backfilling working face using either method. According to the 3D diagram of ground surface subsidence induced by backfilling mining of the coal mine below buildings, railways, and water bodies, the variation in ground surface subsidence could be described by an “L-shaped” curve. The maximum ground surface subsidence was found at the corners. This was because the mining area was larger around the corner of the working face, resulting in a larger influence scope and more severe ground surface subsidence. Besides, multi-slice mining resulted in greater overlying strata movement than single-slice mining.