Analysis on Filling Ratio and Shield Supporting Pressure for Overburden Movement Control in Coal Mining with Compacted Backfilling

Since the weight of overburden is sustained by both the backfill body and the unmined solid coal in coal mining with compacted backfilling (CMCB) panels, the stress and deformation characteristics of the surrounding rocks in coal mining are radically changed. The overburden movement control mechanism by coordinating with backfill body and shield in CMCB was studied systematically in this paper. Based on the analysis of deformational and structural characteristics of surrounding rock in CMCB panels, the methods of theoretical analysis, numerical simulation and engineering test are employed. The results show that the fracture of the main roof is mainly controlled by the filling ratio φ and is non-correlated to the shield supporting pressure p. However, p has a significant control effect on the deflection of roof within the shield canopy length, and adversely affects the filling ratio. With the increase of the filling ratio of the gob, the maximum sagging of the immediate and the main roofs, the peak front and the influence range of the abutment pressures are gradually reduced. Correspondingly, the stable period of internal pressure of backfill body in the gob is shortened. Engineering practice shows that the sagging of the gob roof, the distribution of the abutment pressure, the distribution of the internal pressure in the backfill body, and the ground surface sagging results obtained by the in-situ measurement are approximately corresponding to the theoretical analysis and numerical simulation results.


Introduction
With the development of mining technology, several methods of refilling or backfilling the developed voids, such as the hydraulic backfill, cemented paste backfill, rock backfill and silica alumina-based backill methods have been proposed [1].Hydraulic backfill uses water as the transportation medium to convey the hydraulic backfill materials, such as waste tailings, water hydrophilic slag, mountain sand, river sand, and crushing sand to fill underground mined voids like stopes [2].Hydraulic backfill can be placed cemented or non-cemented, of which the latter one is one of the cheapest methods in mines when small size waste particles are available [3].Cemented paste backfill (CPB) is a non-homogenous material made by mixing waste tailings, water, and cement.CPB is one of new and fast growing waste management and backfill technologies used in the mining industry, which creates an economical way of providing support and safe working place in underground mining development beside its environmental benefits [4,5].Rock backfill can be described as a technology for transportation of backfill forming components such as stone, gravel, soil, industrial solid waste Energies 2017, 10, 31 2 of 23 using manpower, gravity or machinery equipment to fill underground mined voids and production of compressed backfill body.This method decreases waste material on the surface and expands usable land on the ground, decreases environmental pollution by transferring waste rocks to deeper levels out of rain contact, increases stability of mined area and decreases land subsidence and rock bursts due to stress pattern changes [6].
Coal mining with compacted backfilling (CMCB) technology is a "green mining" technology, which is proposed for the problems of coal seams under structures, under underground aquifers and under infrastructures (3U coal seam), and solvents of the gangue in coal mining [7,8].Solid wastes, such as gangue, aeolian sand and open pit mine slags are backfilled into the gob.Simultaneously, the purpose of replacing the 3U coal seam, controlling the overburden movement, and protection of surface buildings are achieved [9,10].In CMCB panels, the backfill body serves as a supporting structure and sustains most of the load from the overburden above the gob.It changes the characteristics of stress redistribution in surrounding rock, and restricts the crack widely development and large area caving of roof.Different from the traditional coal caving mining method, the weight of overburden is carried by both the backfill body and the unmined solid coal in CMCB, and the stress distribution and deformation characteristics of surrounding rocks in coal mining are changed radically [11][12][13].
To date, most theoretical studies have focused on the characters and features of overburden movement in CMCB.The equivalent mining height theory of backfill mining [7] was proposed for the key stratum foundation beam theory [14] and the elastic foundation beam of the roof theory [15], etc.In CMCB, the backfilling degree of the gob can be measured by the filling ratio, which is the value of the final height of the compacted backfill body to the mining height [16].These studies provided useful guidance in understanding the deformational characteristics of overburden at the middle of the panel, but still have their limits.Therefore, an elastic-thin-plate mechanical model is not only established, but also solved in this paper.Meanwhile, overburden movement control mechanism by coordination with backfill body and shield in CMCB were also studied systematically by using numerical simulations.Thereafter, the results of the elastic-thin-plate mechanical model derivation are verified by an engineering practice.The conclusions could serve as reference and guidance for further study on the theory of the overburden movement in CMCB.

Overburden Deformational Characteristics in Traditional Caving Method
In the traditional caving method, periodic fracturing and caving of the immediate and the main roofs were observed with the advancement of the working face.When the gob area was large enough, the resulted disturbance due to the underground excavation may reach the ground surface, creating a surface sagging basin whose area was greater than the gob area.The overburden can be divided into "three horizontal zones" and "three vertical zones" based on their movement and failure characteristics: the coal wall abutment pressure affecting zone, the separation zone and the recompaction zone in horizontal direction; and the caving zone, the fractured zone and the continuous deformation zone in the vertical direction.Figure 1 shows a diagram of the spatial relationship of the "three horizontal zones" and the "three vertical zones" [17].

Overburden Deformation Characteristics in Coal Mining with Compacted Backfilling (CMCB) Method
In CMCB panels, backfill body such as gangue was timely filled into the gob.After compacting by the compactor of the shield [9], the backfill body would have full contact with the roof, thus, it can be seen as a supporting structure.Once in place, it will sustain the load of overburden above the gob, and further, restrict the sagging of the immediate and the main roofs, and prevent them from fracturing and caving.The overburden movement is featured by bending, deflection and local fractures, without caving zone.Consequently, the deformation behaviors of surrounding rocks were totally changed [18], as seen in Figure 2. Features of overburden strata movement in coal mining with compacted backfilling (CMCB) method [11].

Establishment of Elastic-Thin-Plate Mechanical Model of Overburden Strata Movement in CMCB
According to the Winkler hypothesis [19], if a plate were placed on a continuous elastic foundation and bent under the load [q(x,y)] perpendicular to the plate, the contact pressure, R, is proportional to the deflection, w, given the condition that the deflection is far smaller than the thickness of the plate:

R kw =
(1) where k refers to the elastic foundation coefficient.The load on a point on the plate, Q, can be calculated as follows: ( , ) Q q x y kw = − (2) The deflection of the elastic foundation would be consistent with the following differential equation:

Overburden Deformation Characteristics in Coal Mining with Compacted Backfilling (CMCB) Method
In CMCB panels, backfill body such as gangue was timely filled into the gob.After compacting by the compactor of the shield [9], the backfill body would have full contact with the roof, thus, it can be seen as a supporting structure.Once in place, it will sustain the load of overburden above the gob, and further, restrict the sagging of the immediate and the main roofs, and prevent them from fracturing and caving.The overburden movement is featured by bending, deflection and local fractures, without caving zone.Consequently, the deformation behaviors of surrounding rocks were totally changed [18], as seen in Figure 2.

Overburden Deformation Characteristics in Coal Mining with Compacted Backfilling (CMCB) Method
In CMCB panels, backfill body such as gangue was timely filled into the gob.After compacting by the compactor of the shield [9], the backfill body would have full contact with the roof, thus, it can be seen as a supporting structure.Once in place, it will sustain the load of overburden above the gob, and further, restrict the sagging of the immediate and the main roofs, and prevent them from fracturing and caving.The overburden movement is featured by bending, deflection and local fractures, without caving zone.Consequently, the deformation behaviors of surrounding rocks were totally changed [18], as seen in Figure 2.

Establishment of Elastic-Thin-Plate Mechanical Model of Overburden Strata Movement in CMCB
According to the Winkler hypothesis [19], if a plate were placed on a continuous elastic foundation and bent under the load [q(x,y)] perpendicular to the plate, the contact pressure, R, is proportional to the deflection, w, given the condition that the deflection is far smaller than the thickness of the plate:

R kw =
(1) where k refers to the elastic foundation coefficient.The load on a point on the plate, Q, can be calculated as follows: ( , ) The deflection of the elastic foundation would be consistent with the following differential equation:

Establishment of Elastic-Thin-Plate Mechanical Model of Overburden Strata Movement in CMCB
According to the Winkler hypothesis [19], if a plate were placed on a continuous elastic foundation and bent under the load [q(x,y)] perpendicular to the plate, the contact pressure, R, is proportional to the deflection, w, given the condition that the deflection is far smaller than the thickness of the plate: where k refers to the elastic foundation coefficient.The load on a point on the plate, Q, can be calculated as follows: The deflection of the elastic foundation would be consistent with the following differential equation: where D refers to the flexural rigidity of the plate, D = Eh 3 12(1−v 2 ) ; E refers to the elastic modulus of the plate; υ refers to Poisson ratio; h refers to thickness of the plate.
In CMCB, the typical length of panel, l, was around 100-150 m, the advancement length, L, was around 1000-2000 m, and the thickness, h, of the main roof was about 3-20 m.Thus the value of h/l can be estimated to be between 1/50 and 1/5.According to the theory of elastic-thin-plate [20,21], the main roof can be categorized as an elastic-thin-plate.
The main roof thus was regarded as an elastic-thin-plate with fixed boundary conditions at the four edges of the plate.The loads on the elastic-thin-plate, q, which is a load from the overburden, is assumed to be uniformly distributed.The working pressure of the shield on the main roof could be simplified as a uniform pressure, p.The immediate roof and backfill body can be seen as an elastic foundation.The pressure delivered to the main roof from this elastic foundation is related to the deflection of the main roof and can be calculated using Equation (1).The final simplified model for the main roof in CMCB is illustrated in Figure 3. ( 2 ) ( , ) w w w D q x y k w x x y y where D refers to the flexural rigidity of the plate, ; E refers to the elastic modulus of the plate; υ refers to Poisson ratio; h refers to thickness of the plate.
In CMCB, the typical length of panel, l, was around 100-150 m, the advancement length, L, was around 1000-2000 m, and the thickness, h, of the main roof was about 3-20 m.Thus the value of h/l can be estimated to be between 1/50 and 1/5.According to the theory of elastic-thin-plate [20,21], the main roof can be categorized as an elastic-thin-plate.
The main roof thus was regarded as an elastic-thin-plate with fixed boundary conditions at the four edges of the plate.The loads on the elastic-thin-plate, q, which is a load from the overburden, is assumed to be uniformly distributed.The working pressure of the shield on the main roof could be simplified as a uniform pressure, p.The immediate roof and backfill body can be seen as an elastic foundation.The pressure delivered to the main roof from this elastic foundation is related to the deflection of the main roof and can be calculated using Equation 1.The final simplified model for the main roof in CMCB is illustrated in Figure 3.The establishment of the coordinate system is as shown in Figure 4. Here, l refers to the length of panel, L refers to the advancement length of the panel, L' refers to the shield length, q refers to the uniform load of overburden on the main roof, p refers to the supporting pressure of the shield, k refers to the elastic foundation coefficient and ω refers to the deflection of main roof.The establishment of the coordinate system is as shown in Figure 4. Here, l refers to the length of panel, L refers to the advancement length of the panel, L refers to the shield length, q refers to the uniform load of overburden on the main roof, p refers to the supporting pressure of the shield, k refers to the elastic foundation coefficient and ω refers to the deflection of main roof.Based on the principle of equivalent substitution [22], the fixed boundary conditions at the four edges of the proposed model can be replaced by the hinged boundary with an unknown distribution of bending moment These conditions will stay true in the equivalent model.The simplified mechanical model after the equivalent substitution is also shown in Figure 4.

Solving of the Elastic-Thin-Plate Model
According to the simplified elastic-thin-plate mechanical model, the double sine series selected by Navier [23,24] was used to express the deflection equation of the plate: Here, the coefficient amn could be calculated by the principle of virtual displacement [25].The first derivative of deflection along y = 0, y = l and x = 0, x = L was calculated as follows: The main roof, shield, and elastic foundation were considered as an integrated system.Using the principle of virtual displacement, the sum of the work done by the uniform load of the overburden, the moment acting on the hinged boundaries on the virtual displacement, and the increment of deformation energy of the entire system equals to zero.
Here, the bending deformation energy dV1 at a point in the plate was calculated as follows: These conditions will stay true in the equivalent model.The simplified mechanical model after the equivalent substitution is also shown in Figure 4.

Solving of the Elastic-Thin-Plate Model
According to the simplified elastic-thin-plate mechanical model, the double sine series selected by Navier [23,24] was used to express the deflection equation of the plate: Here, the coefficient a mn could be calculated by the principle of virtual displacement [25].The first derivative of deflection along y = 0, y = l and x = 0, x = L was calculated as follows: Energies 2017, 10, 31 6 of 23 The main roof, shield, and elastic foundation were considered as an integrated system.Using the principle of virtual displacement, the sum of the work done by the uniform load of the overburden, the moment acting on the hinged boundaries on the virtual displacement, and the increment of deformation energy of the entire system equals to zero.
Here, the bending deformation energy dV 1 at a point in the plate was calculated as follows: ]dxdy (7) where M x , M y , M xy refer to bending moments, which can be expressed as follows.
The deformation energy, V 2 , of the elastic foundation can be calculated by Equation (10).
The deformation energy, V 3 , of the shield can be calculated by Equation (11).
The total deformation energy of the system can be calculated by when a mn V was increased by δa mn , the increment of the deformation energy of the system was calculated as follows: The virtual work, δT, done by the uniform load q of overburden, and the bending moments M(x)| y = 0 y = l and M(y)| x = 0 x = L on the hinged boundary can be expressed as follows.
In Equation ( 13), E m and F n refer to the coefficients in terms of that M(x)| y = 0 were expanded into a Fourier series [26]: According to the virtual displacement principle and by reorganizing equation, the expression for a mn could be obtained: After substituting Equation (15) into Equation ( 5), the expression for the deflection becomes: Partially differentiate Equation ( 16) by y and by x, the partial derivative of deflection on boundary x = 0 and y = 0 can be expressed by Equation ( 17): Energies 2017, 10, 31 8 of 23 Since the conditions shown in Equation ( 4) would be maintained in the simplified model: Equating the right hand side of Equation ( 17) with zero and getting: In Equation ( 19), the Fourier coefficients E 1 , E 3 , E 5 , E 7 , . . .and F 1 , F 3 , F 5 , F 7 , . . .could be solved through the successive approximation method.
The relationship between the mining height, the filling ratio and the final sagging of the roof can be described by Equation ( 20): where h refers to the mining height, ϕ refers to the filling ratio, and ∆ refers to the final sagging of roof.
The relationship between the elastic foundation coefficient and the final sagging of the roof is shown as follows: where σ 0 refers to in-situ stress, and k refers to the elastic foundation coefficient.The relationship between the elastic foundation coefficient and filling ratio, ϕ, was then deduced from Equations ( 20) and ( 21): In a certain geological condition, σ 0 and h are constant, thus, the elastic foundation coefficient is mainly related to the filling ratio of the gob in CMCB.
Substitute k in Equation ( 16) with the expression in Equation ( 22), the deflection of the plate can now be expressed as follows: Energies 2017, 10, 31 From Equation ( 23), the expression of the bending moments on the plate, M x and M y , can be derived:

Critical Conditions of Main Roof Failure
The maximum stress component within the plate: Through analysis of Equations ( 24) and ( 25), it was found that the maximum bending moment was on the fixed side.The bending moment can be expressed in Fourier series and the expression is shown in Equation ( 14) in previous sections.
M x,max and M y,max were replaced in Equation ( 26) with the expression, as shown in Equation ( 14).The expression for maximum tensile stress within the plate can thus be derived.The main roof would stay stable if the maximum tensile stress in the plate is less than the tensile strength of the main roof:

Application of Analysis Result on an Actual Case
From the deflection and maximum tensile stress equations of the main roof in CMCB, the main factors influencing the deflection and maximum tensile stress of main roof include: filling ratio, ϕ, the length of the panel, l, the advancement length, L, the elastic modulus of the main roof, E, load of overburden above the main roof, q, and the working pressure of the shield, p.Under a certain mining geological condition l, L, h, E, q are approximately constant, so the filling ratio ϕ and the working pressure of the shield p are the most important factors influencing the movement of the roof.
In this section, we will analyze the cooperative roof control effect of the filling ratio and the working pressure of the shield in CMCB, using the field data from Panel 13080 of Pingdingshan No. 12 Coal Mine.In this field case, the length of the panel, l, was 100 m; the advancement length, L, was 350 m; the thickness of the coal seam was 3.3 m; the dip angle of the coal seam was found to be 0 • -4 • .The main roof was a sandstone layer with a thickness of 31.5 m.The elastic modulus Energies 2017, 10, 31 10 of 23 E of the sandstone was 28 GPa, and the Poisson's ratio ν was 0.25, while the uniform load q from the overburden was 4.0 MPa (the depth of cover between main roof and main key stratum, H, was 156.9 m, and the average unit weight γ of the overburden strata was 25 kN/m 3 ) [27,28], The shield canopy length L was 8 m.
(1) Effect of filling ratio and working pressure of shield on main roof stability The first two terms of the Fourier expansion were used to express the bending moments on the fixed boundaries, M(x)| y = 0 . The solving of E 1 , E 3 , F 1 and F 3 was done by the mathematical tool software, Maple.Replace the E 1 , E 3 , F 1 and F 3 in Equation ( 26) and get the expression for maximum tensile stress in the main roof, σ t,max , as a function of shield supporting pressure, p, and the filling ratio, ϕ: The relationship between the maximum tensile stress in the main roof and the shield supporting pressure and filling ratio drawn by the Maple software was shown in Figure 5.
As shown in Figure 5, when the filling ratio of the gob remains constant, increasing or decreasing the shield supporting pressure, p, will have little influence on the maximum tensile stress of the main roof.Therefore, the influence of the shield supporting pressure, p, on the maximum tensile stress of the main roof could be neglected under a determined filling ratio.The maximum tensile stress of the main roof in Panel 13080 was determined to be 6 MPa.From Figure 5, it could be concluded that the filling ratio of the gob should be greater than 82%, i.e., when ϕ ≥ 82%, the failure of the main roof is unlikely to occur.The first two terms of the Fourier expansion were used to express the bending moments on the fixed boundaries,  26) and get the expression for maximum tensile stress in the main roof, σt,max, as a function of shield supporting pressure, p, and the filling ratio, φ: The relationship between the maximum tensile stress in the main roof and the shield supporting pressure and filling ratio drawn by the Maple software was shown in Figure 5.
As shown in Figure 5, when the filling ratio of the gob remains constant, increasing or decreasing the shield supporting pressure, p, will have little influence on the maximum tensile stress of the main roof.Therefore, the influence of the shield supporting pressure, p, on the maximum tensile stress of the main roof could be neglected under a determined filling ratio.The maximum tensile stress of the main roof in Panel 13080 was determined to be 6 MPa.From Figure 5, it could be concluded that the filling ratio of the gob should be greater than 82%, i.e., when φ ≥ 82%, the failure of the main roof is unlikely to occur.(2) Deflection of roof within the shield canopy length before backfilling Since the deformation of the roof is not reversible, the deflection of the roof within the shield canopy length before backfilling is critical for the filling ratio.Based on the theory of the plate, the maximum deflection of the roof within the shield canopy could be observed at x = 8 m, y = 50 m, i.e., the middle of the panel and the end of the shield canopy.
To solve the maximum deflection, we substituted φ = 82%, E1, E3, F1, F3, x = 8 m and y = 50 m into Equation (23).The relationship between the maximum deflection wmax and the shield supporting pressure p was as follows: +5.87 10 0.119 0.82 10 Based on Equation ( 29), the curve of the maximum deflection of main roof versus shield supporting pressure was drawn using the Maple software, as shown in Figure 6.(2) Deflection of roof within the shield canopy length before backfilling Since the deformation of the roof is not reversible, the deflection of the roof within the shield canopy length before backfilling is critical for the filling ratio.Based on the theory of the plate, the maximum deflection of the roof within the shield canopy could be observed at x = 8 m, y = 50 m, i.e., the middle of the panel and the end of the shield canopy.
To solve the maximum deflection, we substituted ϕ = 82%, E 1 , E 3 , F 1 , F 3 , x = 8 m and y = 50 m into Equation (23).The relationship between the maximum deflection w max and the shield supporting pressure p was as follows: Based on Equation ( 29), the curve of the maximum deflection of main roof versus shield supporting pressure was drawn using the Maple software, as shown in Figure 6.The following conclusions could be drawn based on Figure 6: (1) When the filling ratio of the gob was constant, the relationship between the maximum deflection of the main roof within shield canopy length before backfilling, and the shield supporting pressure can be represented by a hyperbolic curve, which can be expressed by function y = 54.10781+ 2809.71225× e −x/0.16579 ; (2) The deflection of the main roof within the shield canopy length before backfilling decreases apparently with the increase of the shield supporting pressure.When the shield supporting pressure increases from 0.4 MPa to 0.9 MPa, the deflection of the main roof decreases by 77%, dropping from 350 mm to 80 mm, which would highly increase the backfilling height, and promote the filling ratio to a large extent.While the shield supporting pressure increases from 0.9 MPa to 1.8 MPa (implying high investment of shield), the decrease in roof deflection is relatively small, dropping from 80 mm to 54 mm, indicating a low increase of backfilling height and filling ratio.Based on this observation, it can be concluded that shield supporting pressure has a significant influence on the roof sagging within the shield canopy length, and adversely affect the filling ratio.

The Numerical Simulation of the Overburden Movement in CMCB
The previous results showed that the filling ratio φ and the shield supporting pressure p are the key factors influencing the movement of the roof in CMCB.And the filling ratio is the critical factor determining the fracture of the main roof.Therefore, this chapter demonstrated the movement characteristics of overburden with the filling ratios of 50%, 65%, 70%, 75%, 80%, 85%, 90%, employing numerical simulations.

Establishment of Numerical Model
The FLAC 3D numerical model is established based on the engineering geological conditions of Panel 13080 employed CMCB of Pingdingshan No. 12 Coal Mine, with a size of 550 m in length, 200 m in width, and 107 m in height, as shown in Figure 7.The following conclusions could be drawn based on Figure 6: (1) When the filling ratio of the gob was constant, the relationship between the maximum deflection of the main roof within shield canopy length before backfilling, and the shield supporting pressure can be represented by a hyperbolic curve, which can be expressed by function y = 54.10781+ 2809.71225× e −x/0.16579 ; (2) The deflection of the main roof within the shield canopy length before backfilling decreases apparently with the increase of the shield supporting pressure.When the shield supporting pressure increases from 0.4 MPa to 0.9 MPa, the deflection of the main roof decreases by 77%, dropping from 350 mm to 80 mm, which would highly increase the backfilling height, and promote the filling ratio to a large extent.While the shield supporting pressure increases from 0.9 MPa to 1.8 MPa (implying high investment of shield), the decrease in roof deflection is relatively small, dropping from 80 mm to 54 mm, indicating a low increase of backfilling height and filling ratio.Based on this observation, it can be concluded that shield supporting pressure has a significant influence on the roof sagging within the shield canopy length, and adversely affect the filling ratio.

The Numerical Simulation of the Overburden Movement in CMCB
The previous results showed that the filling ratio ϕ and the shield supporting pressure p are the key factors influencing the movement of the roof in CMCB.And the filling ratio is the critical factor determining the fracture of the main roof.Therefore, this chapter demonstrated the movement characteristics of overburden with the filling ratios of 50%, 65%, 70%, 75%, 80%, 85%, 90%, employing numerical simulations.

Establishment of Numerical Model
The FLAC 3D numerical model is established based on the engineering geological conditions of Panel 13080 employed CMCB of Pingdingshan No. 12 Coal Mine, with a size of 550 m in length, 200 m in width, and 107 m in height, as shown in Figure 7.The final monitoring data sets of the displacement were obtained after the stabile movement of the overburden.The final sagging curves of the immediate and the main roofs with different filling ratios versus the distance to the working face are shown in Figures 9 and 10.The maximum sagging of immediate roof and main roof versus different filling ratios are shown in Figure 11.The final monitoring data sets of the displacement were obtained after the stabile movement of the overburden.The final sagging curves of the immediate and the main roofs with different filling ratios versus the distance to the working face are shown in Figures 9 and 10.The maximum sagging of immediate roof and main roof versus different filling ratios are shown in Figure 11.The final monitoring data sets of the displacement were obtained after the stabile movement of the overburden.The final sagging curves of the immediate and the main roofs with different filling ratios versus the distance to the working face are shown in Figures 9 and 10.The maximum sagging of immediate roof and main roof versus different filling ratios are shown in Figure 11.(1) The final sagging curves of the immediate and main roofs are a bowl shaped along the advancing direction of panel.The maximum sagging of the immediate and main roofs both appeared in the middle of the gob, where the platform area appeared.The final sagging of the immediate and the main roofs gradually decreases from the middle of the gob to the ends and is close to zero at the position of the open-off cut and the stopping line.(2) The maximum sagging of immediate roof decreases with the increase of the filling ratio of gob.
When the filling ratio is 50%, 65%, 75%, 85%, 90%, the maximum sagging of immediate roof is 1625 mm, 1198 mm, 876 mm, 501 mm, 315 mm, respectively.(3) The maximum sagging of main roof decreases with the increase of the filling ratio of the gob.
When the filling ratio is 50%, 65%, 75%, 85%, 90%, the maximum sagging of main roof is 1338 mm, 1023 mm, 609 mm, 263 mm, 163 mm, respectively.When the filling ratio increased from 75% to 85%, the maximum sagging is severely decreased from 609 mm to 263 mm.It concludes that the main roof could remain continuously deformed without breaking when the filling ratio is above 85%, which is consistent with the elastic-thin-plate model.

Distributions of Abutment Pressure and Backfill Body Internal Pressure under Different Filling Ratios
(1) The distribution characteristics of abutment pressure With the working face advancing a certain distance, the monitoring data of abutment pressure could be illustrated into a curve versus distance to the working face.The abutment pressure distributions under different filling ratios are shown in Figure 12.  (1) The final sagging curves of the immediate and main roofs are a bowl shaped along the advancing direction of panel.The maximum sagging of the immediate and main roofs both appeared in the middle of the gob, where the platform area appeared.The final sagging of the immediate and the main roofs gradually decreases from the middle of the gob to the ends and is close to zero at the position of the open-off cut and the stopping line.(2) The maximum sagging of immediate roof decreases with the increase of the filling ratio of gob.
When the filling ratio is 50%, 65%, 75%, 85%, 90%, the maximum sagging of immediate roof is 1625 mm, 1198 mm, 876 mm, 501 mm, 315 mm, respectively.(3) The maximum sagging of main roof decreases with the increase of the filling ratio of the gob.
When the filling ratio is 50%, 65%, 75%, 85%, 90%, the maximum sagging of main roof is 1338 mm, 1023 mm, 609 mm, 263 mm, 163 mm, respectively.When the filling ratio increased from 75% to 85%, the maximum sagging is severely decreased from 609 mm to 263 mm.It concludes that the main roof could remain continuously deformed without breaking when the filling ratio is above 85%, which is consistent with the elastic-thin-plate model.

Distributions of Abutment Pressure and Backfill Body Internal Pressure under Different Filling Ratios
(1) The distribution characteristics of abutment pressure With the working face advancing a certain distance, the monitoring data of abutment pressure could be illustrated into a curve versus distance to the working face.The abutment pressure distributions under different filling ratios are shown in Figure 12. (1) The final sagging curves of the immediate and main roofs are a bowl shaped along the advancing direction of panel.The maximum sagging of the immediate and main roofs both appeared in the middle of the gob, where the platform area appeared.The final sagging of the immediate and the main roofs gradually decreases from the middle of the gob to the ends and is close to zero at the position of the open-off cut and the stopping line.
(2) The maximum sagging of immediate roof decreases with the increase of the filling ratio of gob.
When the filling ratio is 50%, 65%, 75%, 85%, 90%, the maximum sagging of immediate roof is 1625 mm, 1198 mm, 876 mm, 501 mm, 315 mm, respectively.(3) The maximum sagging of main roof decreases with the increase of the filling ratio of the gob.
When the filling ratio is 50%, 65%, 75%, 85%, 90%, the maximum sagging of main roof is 1338 mm, 1023 mm, 609 mm, 263 mm, 163 mm, respectively.When the filling ratio increased from 75% to 85%, the maximum sagging is severely decreased from 609 mm to 263 mm.It concludes that the main roof could remain continuously deformed without breaking when the filling ratio is above 85%, which is consistent with the elastic-thin-plate model.

Distributions of Abutment Pressure and Backfill Body Internal Pressure under Different Filling Ratios
(1) The distribution characteristics of abutment pressure With the working face advancing a certain distance, the monitoring data of abutment pressure could be illustrated into a curve versus distance to the working face.The abutment pressure distributions under different filling ratios are shown in Figure 12.The following conclusions could be indicated based on Figure 12: under different filling ratios, the peak front of abutment pressure appears in the range of 10-20 m from the working face.When the filling ratio is 50%, the peak front value of the abutment pressure is 16 MPa with an influence range 25 m, while the filling ratio 85% with the peak front 11.5 MPa and an influence range 16 m.Thus, the filling ratio of gob has a significant influence on the value of abutment pressure and influence range, i.e., a higher the ratio results in a lower peak front and a smaller influence range.Additionally, when the filling ratio increases from 80% to 85%, the peak front and influence range decrease at a large extent, indicating that the main roof remains continuously deformed when the filling ratio is above 85%, without broken.
(2) The distribution characteristics of backfill body internal pressure According to data from the stress measuring line set at the middle of the backfill body along the advance direction of the panel, the curves of the backfill body internal pressure under different ratios versus the working face advancing distance are illustrated in Figure 13.The stable periods of backfill body internal pressure under different ratios are shown in Figure 14.
Energies 2017, 10, 31 15 of 23 The following conclusions could be indicated based on Figure 12: under different filling ratios, the peak front of abutment pressure appears in the range of 10-20 m from the working face.When the filling ratio is 50%, the peak front value of the abutment pressure is 16 MPa with an influence range 25 m, while the filling ratio 85% with the peak front 11.5 MPa and an influence range 16 m.Thus, the filling ratio of gob has a significant influence on the value of abutment pressure and influence range, i.e., a higher the ratio results in a lower peak front and a smaller influence range.Additionally, when the filling ratio increases from 80% to 85%, the peak front and influence range decrease at a large extent, indicating that the main roof remains continuously deformed when the filling ratio is above 85%, without broken.
(2) The distribution characteristics of backfill body internal pressure According to data from the stress measuring line set at the middle of the backfill body along the advance direction of the panel, the curves of the backfill body internal pressure under different ratios versus the working face advancing distance are illustrated in Figure 13.The stable periods of backfill body internal pressure under different ratios are shown in Figure 14.The following conclusions could be implied based on Figures 13 and 14: (a) The backfill body internal pressure gradually increased with the increase of the advancing distance of the working face.The increase pressure curves could be divided into three phases, i.e., initial pressure phase, accelerated increasing phase, and the stable pressure phase.In the initial pressure phase, the backfill body internal pressure slowly increases when the advancing The following conclusions could be indicated based on Figure 12: under different filling ratios, the peak front of abutment pressure appears in the range of 10-20 m from the working face.When the filling ratio is 50%, the peak front value of the abutment pressure is 16 MPa with an influence range 25 m, while the filling ratio 85% with the peak front 11.5 MPa and an influence range 16 m.Thus, the filling ratio of gob has a significant influence on the value of abutment pressure and influence range, i.e., a higher the ratio results in a lower peak front and a smaller influence range.Additionally, when the filling ratio increases from 80% to 85%, the peak front and influence range decrease at a large extent, indicating that the main roof remains continuously deformed when the filling ratio is above 85%, without broken.
(2) The distribution characteristics of backfill body internal pressure According to data from the stress measuring line set at the middle of the backfill body along the advance direction of the panel, the curves of the backfill body internal pressure under different ratios versus the working face advancing distance are illustrated in Figure 13.The stable periods of backfill body internal pressure under different ratios are shown in Figure 14.The following conclusions could be implied based on Figures 13 and 14: (a) The backfill body internal pressure gradually increased with the increase of the advancing distance of the working face.The increase pressure curves could be divided into three phases, i.e., initial pressure phase, accelerated increasing phase, and the stable pressure phase.In the initial pressure phase, the backfill body internal pressure slowly increases when the advancing The following conclusions could be implied based on Figures 13 and 14: (a) The backfill body internal pressure gradually increased with the increase of the advancing distance of the working face.The increase pressure curves could be divided into three phases, i.e., initial pressure phase, accelerated increasing phase, and the stable pressure phase.In the initial pressure phase, the backfill body internal pressure slowly increases when the advancing distance of the working face increases, which indicates that the overburden subjected to the supports of front solid coal and rear backfill body bends slowly to the gob.Within the accelerated increasing phase, while the advancing distance of the face and span of the gob increase, the overburden accelerately sags to the gob, leading to the rapidly increase of the backfill body internal pressure.After the accelerated increasing phase, the internal pressure goes into a stable pressure phase, in which the internal pressure gradually approaches to the in-situ stress and the movement of the overburden reaches a new equilibrium.(b) When the filling ratio is 50%, 65%, 70%, 75%, 80%, 85%, 90%, the advancing distance of the initial pressure phase is 68 m, 59 m, 52 m, 46 m, 43 m, 30 m, 20 m, respectively, and the advancing distance at which the internal pressure reaches stable is 105 m, 82 m, 78 m, 74 m, 72 m, 70 m, 69 m.Thus, the bigger the filling ratio is, the faster the backfill body internal pressure reaches stable period.

Background
The underground layout of Panel 13080 employed CMCB of Pingdingshan No. 12 Coal Mine, and buildings location on or nearby the ground of the panel are shown in Figure 15.There are two villages and two mine surface plants located on or nearby the surface of Panel 13080.Donggaohuang Village lies on the west surface of Panel 13080, while Mifengwang Village lies on the east.Mine surface plant of No. 12 lies on the north surface of Panel 13080, while mine surface plant of Hezhuang Coal Mine lies on the south.The backfill material is from gangue dump of No. 12 Coal Mine.According to the former analysis, the designed filling ratio is 85%.
Energies 2017, 10, 31 16 of 23 distance of the working face increases, which indicates that the overburden subjected to the supports of front solid coal and rear backfill body bends slowly to the gob.Within the accelerated increasing phase, while the advancing distance of the face and span of the gob increase, the overburden accelerately sags to the gob, leading to the rapidly increase of the backfill body internal pressure.After the accelerated increasing phase, the internal pressure goes into a stable pressure phase, in which the internal pressure gradually approaches to the in-situ stress and the movement of the overburden reaches a new equilibrium.(b) When the filling ratio is 50%, 65%, 70%, 75%, 80%, 85%, 90%, the advancing distance of the initial pressure phase is 68 m, 59 m, 52 m, 46 m, 43 m, 30 m, 20 m, respectively, and the advancing distance at which the internal pressure reaches stable is 105 m, 82 m, 78 m, 74 m, 72 m, 70 m, 69 m.Thus, the bigger the filling ratio is, the faster the backfill body internal pressure reaches stable period.

Background
The underground layout of Panel 13080 employed CMCB of Pingdingshan No. 12 Coal Mine, and buildings location on or nearby the ground of the panel are shown in Figure 15.There are two villages and two mine surface plants located on or nearby the surface of Panel 13080.Donggaohuang Village lies on the west surface of Panel 13080, while Mifengwang Village lies on the east.Mine surface plant of No. 12 lies on the north surface of Panel 13080, while mine surface plant of Hezhuang Coal Mine lies on the south.The backfill material is from gangue dump of No. 12 Coal Mine.According to the former analysis, the designed filling ratio is 85%.

Monitoring Layouts and Instruments
(1) Layout of roof sagging, pressure monitoring and instruments With the advancement and backfilling of Panel 13080, the gob roof sagging, abutment pressure and backfill body internal pressure were monitored.The monitoring layouts are shown in Figure 16, and the monitoring instruments are shown in Figure 17, with parameters shown in Table 2.As shown in Figure 16, #1 to #5 monitors for sagging of gob roof and #1 to #5 monitors for the backfill body internal pressure were installed when the working face advanced 20 m; #6 to #10 monitors for sagging of gob roof and #6 to #10 monitors for the backfill body internal pressure were installed when the working face advanced 40 m.

Monitoring Layouts and Instruments
(1) Layout of roof sagging, pressure monitoring and instruments With the advancement and backfilling of Panel 13080, the gob roof sagging, abutment pressure and backfill body internal pressure were monitored.The monitoring layouts are shown in Figure 16, and the monitoring instruments are shown in Figure 17, with parameters shown in Table 2.As shown in Figure 16, #1 to #5 monitors for sagging of gob roof and #1 to #5 monitors for the backfill body internal pressure were installed when the working face advanced 20 m; #6 to #10 monitors for sagging of gob roof and #6 to #10 monitors for the backfill body internal pressure were installed when the working face advanced 40 m.(2) Surface Sagging Monitoring Layout and instrument According to the design principles of surface sagging observation [22], a surface sagging monitoring line was set along the dip direction of the panel, considering the distribution of the buildings and roads.As shown in Figure 18, the surface sagging monitoring line was 1025 m, with 3 control observation (R1, R2, R3, 50 m of interval) and 36 dynamic observation (25 m of interval).The robotic total station was employed in the surface sagging monitoring.(2) Surface Sagging Monitoring Layout and instrument According to the design principles of surface sagging observation [22], a surface sagging monitoring line was set along the dip direction of the panel, considering the distribution of the buildings and roads.As shown in Figure 18, the surface sagging monitoring line was 1025 m, with 3 control observation (R1, R2, R3, 50 m of interval) and 36 dynamic observation (25 m of interval).The robotic total station was employed in the surface sagging monitoring.(2) Surface Sagging Monitoring Layout and instrument According to the design principles of surface sagging observation [22], a surface sagging monitoring line was set along the dip direction of the panel, considering the distribution of the buildings and roads.As shown in Figure 18, the surface sagging monitoring line was 1025 m, with 3 control observation (R 1 , R 2 , R 3 , 50 m of interval) and 36 dynamic observation (25 m of interval).The robotic total station was employed in the surface sagging monitoring.The backfill body mainly subjected to the overburden load.The maximum sagging of #3 and #8 were 340 mm, with the filling ratio of 89.5%, which was consistent with the theoretical results.

Monitoring Results and Analysis
(2) Abutment pressure The abutment pressure in different depth of solid coal was monitored by setting the borehole stress meters in front of working face, shown in Figure 20.(2) Abutment pressure The abutment pressure in different depth of solid coal was monitored by setting the borehole stress meters in front of working face, shown in Figure 20.The following conclusions could be drawn based on Figure 20: (a) The abutment pressure distribution of CMCB panels is similar to the caving mining panels, with 4 zones along the advancement of the face, i.e., de-stress zone (0-5 m), stress rapid increase zone (5-15 m), stress decrease zone (15-25 m) and in-situ zone (after 25 m).(b) However, the peak front and the influence range of the abutment pressures are gradually reduced compared to caving mining method.Furthermore, the maximum stress concentration factor is only 1.3, which is lower than the concentration factor in caving mining panels.(c) When the drilling depth is 15 m, the monitoring results of abutment pressures are most close to the numerical simulation results.But the peak value of abutment pressures are slightly lower in numerical simulation, and the distance of peak position to the working face is large.It is considered that the numerical simulation cannot completely reproduce all the complicated mining geological conditions in engineering practice, but the error is still within the acceptable range.
(  The curves of maximum surface sagging from the data through field measurement and prediction software MSPS [29] with a filling ratio of 0%, 85% were illustrated, as shown in Figure 22 and Table 3.The curves of maximum surface sagging from the data through field measurement and prediction software MSPS [29] with a filling ratio of 0%, 85% were illustrated, as shown in Figure 22 and Table 3.As show in Figure 22 and Table 3, after the mined-out of Panel 13080, the maximum value of surface sagging was 173 mm, and the maximum horizontal deformation is 0.2 mm/m.According to the Building, Railway and Water bodies Pillars Remain and Mining Act of 2000, damage level of the buildings on the ground is I level (very slightly) without the need for reinforcement or repair, and it's in consistent with the theoretical results.

Conclusions
(1) The elastic-thin-plate mechanical model for overburden in CMCB was established according to the deformation characteristics of the surrounding rocks.The deflection equation, bending moment equation, and maximum tensile stress equations associated with the main roof in CMCB were derived from virtual displacement principle.The critical conditions of main roof failure were also obtained.(2) By theoretical analysis, the main roof stability is primarily controlled by the filling ratio φ and is non-correlated to the shield supporting pressure p.The relationship between the maximum sagging of main roof within the shield canopy and shield supporting pressure was given, which shows that the support pressure has a significant control effect on the sagging of the roof within the shield canopy length, and adversely affects the filling ratio.As show in Figure 22 and Table 3, after the mined-out of Panel 13080, the maximum value of surface sagging was 173 mm, and the maximum horizontal deformation is 0.2 mm/m.According to the Building, Railway and Water bodies Pillars Remain and Mining Act of 2000, damage level of the buildings on the ground is I level (very slightly) without the need for reinforcement or repair, and it's in consistent with the theoretical results.

Conclusions
(1) The elastic-thin-plate mechanical model for overburden in coal mining with compacted backfilling (CMCB) was established according to the deformation characteristics of the surrounding rocks.The deflection equation, bending moment equation, and maximum tensile stress equations associated with the main roof in CMCB were derived from virtual displacement principle.The critical conditions of main roof failure were also obtained.(2) By theoretical analysis, the main roof stability is primarily controlled by the filling ratio ϕ and is non-correlated to the shield supporting pressure p.The relationship between the maximum sagging of main roof within the shield canopy and shield supporting pressure was given, which shows that the support pressure has a significant control effect on the sagging of the roof within the shield canopy length, and adversely affects the filling ratio.(3) The numerical simulation results show that the maximum sagging of immediate roof and main roof decreases with the increase of the filling ratios of gob.When the filling ratio increased from 75% to 85%, the maximum sagging is decreased from 609 mm to 263 mm, respectively.With the increase of the filling ratio of gob, the peak value and the influence range of the abutment pressure are gradually reduced, and the stability period of backfill body internal pressure is shortened gradually.(4) Filling ratio of working face was determined to be 85% according to the theoretical analysis and the numerical calculation results.The abutment pressure, the sagging of roof, the backfill body internal pressure, surface sagging were monitored and analysed.Monitoring results show that the measured data are consistent with the theoretical analysis and numerical simulation results.

Figure 3 .
Figure 3. Elastic-thin-plate mechanical model of main roof in CMCB.(a) 3-D view of model; (b) Section view of model.

Figure 3 .
Figure 3. Elastic-thin-plate mechanical model of main roof in CMCB.(a) 3-D view of model; (b) Section view of model.

Figure 4 .
Figure 4. Simplified elastic-thin-plate mechanical model of main roof.
on those boundaries if the deformation character at those boundaries would maintain.Bending moment 0 ( ) y y l M x = = distribution in both y = 0, y = l sides changes with the x location of each point.Bending moment 0 x = 0, x = l sides changes with the y location of each point.For fixed boundary condition, the deflection and the first derivative of the deflection at the direction normal to the boundary would be zero, as shown in Equation (4):

Figure 4 .
Figure 4. Simplified elastic-thin-plate mechanical model of main roof.

.
The solving of E1, E3, F1 and F3 was done by the mathematical tool software, Maple.Replace the E1, E3, F1 and F3 in Equation (

Figure 5 .
Figure 5. Relationship between the maximum tensile stress and the shield supporting pressure and filling ratio.

Figure 5 .
Figure 5. Relationship between the maximum tensile stress and the shield supporting pressure and filling ratio.

23 Figure 6 .
Figure 6.Relationship between the maximum deflection and the shield supporting pressure.

Figure 6 .
Figure 6.Relationship between the maximum deflection and the shield supporting pressure.

Figure 8 .Figure 8 .
Figure 8. Monitoring points in the profile numerical model.

Figure 9 .
Figure 9.The final sagging curves of the immediate roof with different filling ratios.

Figure 10 .
Figure 10.The final sagging curves of the main roof with different filling ratios.

Figure 9 .
Figure 9.The final sagging curves of the immediate roof with different filling ratios.

Figure 9 .
Figure 9.The final sagging curves of the immediate roof with different filling ratios.

Figure 10 .
Figure 10.The final sagging curves of the main roof with different filling ratios.

Figure 10 .
Figure 10.The final sagging curves of the main roof with different filling ratios.

Figure 11 .
Figure 11.The maximum sagging of the immediate and the main roofs with the filling ratios.

Figure 12 .
Figure 12.The abutment pressure distributions under different filling ratios with the advancement distance of the working face.

18 Figure 11 .
Figure 11.The maximum sagging of the immediate and the main roofs with the filling ratios.

Energies 2017, 10 , 31 14 of 23 Figure 11 .
Figure 11.The maximum sagging of the immediate and the main roofs with the filling ratios.

Figure 12 .
Figure 12.The abutment pressure distributions under different filling ratios with the advancement distance of the working face.

18 Figure 12 .
Figure 12.The abutment pressure distributions under different filling ratios with the advancement distance of the working face.

Figure 13 .
Figure 13.The distribution of backfill body internal pressure under different filling ratios.

Figure 14 .
Figure 14.The distance from working face when backfill body internal pressure is stable under the condition of different filling ratios

Figure 13 .
Figure 13.The distribution of backfill body internal pressure under different filling ratios.

Figure 13 .
Figure 13.The distribution of backfill body internal pressure under different filling ratios.

Figure 14 .
Figure 14.The distance from working face when backfill body internal pressure is stable under the condition of different filling ratios

Figure 14 .
Figure 14.The distance from working face when backfill body internal pressure is stable under the condition of different filling ratios.

Figure 15 .
Figure 15.The surface-underground contrast layout of panel 13080.

Figure 15 .
Figure 15.The surface-underground contrast layout of panel 13080.

Figure 17 .
Figure 17.Mine pressure monitoring instruments.(a) Monitor for sagging of gob roof; (b) Monitor for backfill body internal pressure; (c) Borehole stress meter.

( 1 )
The sagging of the gob roof Taken data collected from #3 roof sagging monitor (in the middle of the face line, 20 m away from the open-cut) and #8 roof sagging monitor (in the middle of the face line, 40 m away from the open-cut) for instance, the curves of the sagging of the gob roof versus advancement distance of the working face are as shown in Figure 19.

Figure 18 .
Figure 18.The sketch of the surface sagging observation line layout.

( 1 )
The sagging of the gob roof Taken data collected from #3 roof sagging monitor (in the middle of the face line, 20 m away from the open-cut) and #8 roof sagging monitor (in the middle of the face line, 40 m away from the open-cut) for instance, the curves of the sagging of the gob roof versus advancement distance of the working face are as shown in Figure 19.Energies 2017, 10, 31 19 of 23 5.3.Monitoring Results and Analysis (1) The sagging of the gob roof Taken data collected from #3 roof sagging monitor (in the middle of the face line, 20 m away from the open-cut) and #8 roof sagging monitor (in the middle of the face line, 40 m away from the open-cut) for instance, the curves of the sagging of the gob roof versus advancement distance of the working face are as shown in Figure 19.

Figure 19 .
Figure 19.The sagging curves of the gob roof versus advancement distance of the working face.

Energies 2017, 10 , 31 19 of 23 Figure 19 .
Figure 19.The sagging curves of the gob roof versus advancement distance of the working face.The following conclusions could be drawn based on Figure 19: (a) The curves of the gob roof sagging versus the working face advancement could be divided into three stages, i.e., the slowly sagging stage (0-20 m from the working face), the accelerated sagging stage (20-65 m from the working face) and the stable sagging stage (above 65 m from the working face), representing the three stages of the gob roof movement.(b) In the slowly sagging stage (0-25 m from the working face), since the distance advancement of the working face was short, the pillar at the open-cut, solid coal subjected to the main load of the overburden, thus, the roof sagged slowly.(c) With the distance advancement of the working face continue increasing, the gob roof acceleratedly sagged mainly bending as an intact stratum, corresponding to the accelerated sagging stage (20-65 m from the working face).(d) After the face advancing 65 m, the sagging of the gob roof approximately verged to a stable stage.The backfill body mainly subjected to the overburden load.The maximum sagging of #3 and #8 were 340 mm, with the filling ratio of 89.5%, which was consistent with the theoretical results.

Figure 20 .
Figure 20.The abutment pressure in different depth of solid coal with the distance to working face.The following conclusions could be drawn based on Figure 20: (a) The abutment pressure distribution of CMCB panels is similar to the caving mining panels, with 4 zones along the advancement of the face, i.e., de-stress zone (0-5 m), stress rapid increase zone (5-15 m), stress decrease zone (15-25 m) and in-situ zone (after 25 m).(b) However, the peak front and the influence range of the abutment pressures are gradually

Figure 20 .
Figure 20.The abutment pressure in different depth of solid coal with the distance to working face.

23 ( 3 )
) Backfill body internal pressure Taken data collected from #3 pressure meter (in the middle of the face line, 20 m away from the open-cut) and #8 pressure meter (in the middle of the face line, 40 m away from the open-cut) for instance, the curves of the backfill body internal pressure versus advancement distance of the working face are shown in Figure 21.Energies 2017, 10, 31 20 of Backfill body internal pressure Taken data collected from #3 pressure meter (in the middle of the face line, 20 m away from the open-cut) and #8 pressure meter (in the middle of the face line, 40 m away from the open-cut) for instance, the curves of the backfill body internal pressure versus advancement distance of the working face are shown in Figure 21.

Figure 21 .
Figure 21.The curves of backfill body internal pressure of monitoring points #3 and #8.The following conclusions could be drawn based on Figure 21: (a) The backfill body internal pressure curves versus the advancement distance of the working face could be divided into three stages: the initial stress stage (0-20 m from the working face), the stress rapidly increasing stage (20-65 m from the working face) and the stress stability stage (after 65 m from the working face), in correlation with the three stages of gob roof sagging, i.e., the slowly sagging stage, the accelerated sagging stage and the stable sagging stage.(b) In initial stress stage, the backfill body internal pressure varied from 1.8 MPa to 2.5 MPa, while in the stress rapidly increasing stage, it increase rapidly from 2.5 MPa to approximate 8.5 MPa.It indicated that the main carriers of the overburden load changed from the solid coal of the face and the open-cut pillars in the initial stage to backfill body in the stress rapidly increasing stage.(c) After 65 m, the backfill body internal pressure was gradually stabilized with the distance advancement of the working face, with a maximum value of approximate 8.6 MPa which is close to the in-situ stress.It shows that the movement of the overburden at this stage was close to stable under the support of the backfill body.(4) Measurement of surface sagging

Figure 21 .
Figure 21.The curves of backfill body internal pressure of monitoring points #3 and #8.

Table 1 .
Mechanical parameters of coal and rock strata.

Table 1 .
Mechanical parameters of coal and rock strata.

Table 1 .
Mechanical parameters of coal and rock strata.

Table 2 .
Selection of mine pressure monitoring instruments.

Table 2 .
Selection of mine pressure monitoring instruments.

Table 2 .
Selection of mine pressure monitoring instruments.

Table 3 .
Data of surface movement prediction and field measurement of Panel 13080.

Table 3 .
Data of surface movement prediction and field measurement of Panel 13080.