Deformation Behavior of Hard Roofs in Solid Backfill Coal Mining Using Physical Models

Solid backfill coal mining technology has been widely applied in coal seams that are at risk of hard roof. Using actual measured strain–stress curves of the backfill body and the similarity theory, this study designed and employed four experimental models for physical simulation, corresponding to roof-controlled backfilling ratios of 0%, 40%, 82.5% and 97% using the geological conditions of Face No. 6304 in the Jining No. 3 coal mine—a solid backfill coal mining face under a hard roof. A non-contact strain measurement system and pressure sensors were used to monitor the deformation of the overlying strata and changes in abutment stress ahead of the face during mining of the models for varying roof-controlled backfilling ratios. The results indicated that the solid backfill body was able to support the roof. As the roof-controlled backfilling ratio was increased, the maximum subsidence of the roof and the maximum height of the cracks decreased. When the roof-controlled backfilling ratio was 82.5% or higher, the working face did not display any obvious initial fractures or periodic fractures, and both the value and the impact range of the abutment stress ahead of the face decreased.


Introduction
A hard roof [1][2][3] refers to the thick strata above a coal seam or above a thin immediate roof.Such a roof generally has high strength and low caving immediately after each mining operation.With the continuing advance of the longwall face and the increase of the overhang area, the stress in the hard roof eventually exceeds its ultimate strength and the roof caves suddenly, which causes a rapid release of the energy stored in the roof and the coal seam.This often results in coal mine catastrophic dynamic events, such as wind blasting, causing serious damage to mining equipment, significant mining delays, and sometimes casualties [4][5][6][7].The hard roof is one of the main causes of a coal mine burst at the longwall face [8,9].It is of great importance to note that in China, coal reserves under hard roofs account for about one-third of total reserves; moreover, nearly 40% of fully mechanized coal mining panels have hard roofs, and more than 50% of mining areas suffer from problems associated with hard roofs [10,11].
In recent years, backfill coal mining [12,13] has become a popular and widely applied mining method for the safe extraction of coal resources trapped under buildings, railways, and water bodies.This technique fills the goaf area with solid waste, such as gangue, to effectively control the movement of both the roof and the overlying strata.Numerous case studies [14][15][16] have demonstrated that backfill coal mining can effectively control the movement of both the roof and the overlying strata.In comparison with the traditional longwall caving method, backfill coal mining can greatly reduce the abutment stress around the excavation surface.For these reasons, Zhou and Zhang [17,18] proposed a solid backfill method for controlling hard-roof-induced dynamic disasters; their method has been verified both in theory and in practice.However, different backfill systems have different abilities to control the deformation of the hard roof and different pressure profiles in the working face.Based on the geological conditions of the 6304-1 working face in the Jisan coal mine, this paper compares the deformation of the hard roof in the working face for four backfill scenarios using physical simulation models (Compared with numerical modelling, physical modelling results are more reliable because numerical modelling has difficulties in parameter selection).The results can provide basic data to inform the design of backfill parameters, the selection of appropriate backfill support, and roadway support.

Principle of Solid Backfilling
The backfilling system was developed to handle solid backfill materials as an integral part of the original fully mechanized coal mining system [19,20].Backfill coal mining technology consists of five main production systems.Its coal transport system, auxiliary transport system, and ventilation system are the same as those used in traditional fully mechanized coal mining.The two additional systems that backfill coal mining employs are a solid waste transport system and a backfill system.The solid waste transport system transports the solid waste from the surface to the underground working face.The backfill system places the dense backfill material in the gob area.
The backfill hydraulic supports [21,22] are key to the success of solid backfilling.As shown in Figure 1, its front canopy supports the roof, providing a safe space to operate mining machines, and its back canopy offers the space required to transport the backfill material, and to dump and compact it in the gob.The backfilling scraper conveyor, used for transporting the solid backfill materials, is hung below the rear roof canopy.Another important piece of equipment in the rear part of the backfill hydraulic support is the tamping arm, which can provide a pressure of 2 MPa to push the backfill materials into the gob and compact them to a sufficient density to effectively support the roof [11].
Energies 2017, 10, 557 2 of 20 overlying strata.In comparison with the traditional longwall caving method, backfill coal mining can greatly reduce the abutment stress around the excavation surface.For these reasons, Zhou and Zhang [17,18] proposed a solid backfill method for controlling hard-roof-induced dynamic disasters; their method has been verified both in theory and in practice.However, different backfill systems have different abilities to control the deformation of the hard roof and different pressure profiles in the working face.Based on the geological conditions of the 6304-1 working face in the Jisan coal mine, this paper compares the deformation of the hard roof in the working face for four backfill scenarios using physical simulation models (Compared with numerical modelling, physical modelling results are more reliable because numerical modelling has difficulties in parameter selection).The results can provide basic data to inform the design of backfill parameters, the selection of appropriate backfill support, and roadway support.

Principle of Solid Backfilling
The backfilling system was developed to handle solid backfill materials as an integral part of the original fully mechanized coal mining system [19,20].Backfill coal mining technology consists of five main production systems.Its coal transport system, auxiliary transport system, and ventilation system are the same as those used in traditional fully mechanized coal mining.The two additional systems that backfill coal mining employs are a solid waste transport system and a backfill system.The solid waste transport system transports the solid waste from the surface to the underground working face.The backfill system places the dense backfill material in the gob area.
The backfill hydraulic supports [21,22] are key to the success of solid backfilling.As shown in Figure 1, its front canopy supports the roof, providing a safe space to operate mining machines, and its back canopy offers the space required to transport the backfill material, and to dump and compact it in the gob.The backfilling scraper conveyor, used for transporting the solid backfill materials, is hung below the rear roof canopy.Another important piece of equipment in the rear part of the backfill hydraulic support is the tamping arm, which can provide a pressure of 2 MPa to push the backfill materials into the gob and compact them to a sufficient density to effectively support the roof [11].

Basic Parameters of the Working Face
Backfill coal mining technology has been applied in Panel No. 6304-1 of the Jisan coal mine, where the No. 3 coal seam of the Shanxi formation is the main mineable coal seam.The average thickness of the No. 3 coal seam is 3.5 m, its average dip is 2.5°, its density is 1.36 t/m 3 , and its average depth is 650 m.It has no immediate roof.The main roof is a fine sandstone of 32.5-49.75m thick (with an average thickness of 41.6 m), the uniaxial compressive strength is 84.4 MPa, the tensile strength is 13.5 MPa, and the elastic modulus is 17.0 GPa.It can be regarded as a hard roof.The panel, No. 6304,

Experimental Installation
Model supports measuring 2.30 m × 0.30 m × 1.80 m were used to perform the experiments, with pressure water bags used as the compensatory loads.Changes in the abutment stress ahead of the working face were monitored using pressure sensors arranged inside the coal seam and a TS3890 static strain indicator (Yangzhou Thai Electronics Co., Ltd., Yangzhou, China).A non-contact strain measurement system was used to monitor displacement variations in the overlying strata.The system used an advanced 2D digital image correlation algorithm, which can measure arbitrary displacement and deformation from 500 micro strain to more than five times strain.The size of the observed model is from 1 mm to 10 m.This monitoring system uses an industrial digital camera (29 million pixels), optical lens, lighting, matching notebook, Flycapture software, and Vic-2D (Beijing Deanwell Technology Co., Ltd., Beijing, China).The Flycapture software relies on the Canadian PointGrey industrial digital camera.Vic-2D can realize the data-output constantly and deal with the data.The arrangement of the measurement and data processing systems is shown in Figure 3.

Experimental Installation
Model supports measuring 2.30 m × 0.30 m × 1.80 m were used to perform the experiments, with pressure water bags used as the compensatory loads.Changes in the abutment stress ahead of the working face were monitored using pressure sensors arranged inside the coal seam and a TS3890 static strain indicator (Yangzhou Thai Electronics Co., Ltd., Yangzhou, China).A non-contact strain measurement system was used to monitor displacement variations in the overlying strata.The system used an advanced 2D digital image correlation algorithm, which can measure arbitrary displacement and deformation from 500 micro strain to more than five times strain.The size of the observed model is from 1 mm to 10 m.This monitoring system uses an industrial digital camera (29 million pixels), optical lens, lighting, matching notebook, Flycapture software, and Vic-2D (Beijing Deanwell Technology Co., Ltd., Beijing, China).The Flycapture software relies on the Canadian PointGrey industrial digital camera.Vic-2D can realize the data-output constantly and deal with the data.The arrangement of the measurement and data processing systems is shown in Figure 3.

Design of Model Parameters
Each model was designed to consist of nine strata including the coal seam, hard roof, immediate roof, and partially overlying strata.The model size was 2.30 m × 0.30 m × 1.23 m (L × W × H).The real height on site was 184.80 m.The rock strata that were not laid were approximately 475 m thick.According to the stress similarity ratio, a compensatory load of 0.0475 MPa was required, which was applied using pressure water bags.
The height, strength, and material mixing proportion for each rock stratum in the models were calculated based on the site test results of the actual lithology and the principle of similarity.They are shown in Table 1.
The rock masses in the models were laid from bottom to top.Before the placement of each stratum, sand, calcium carbonate, plaster, and water were mixed in a blender with the mass of each component as shown in Table 2.After mixing, the materials were spread evenly and smoothly on the model support and tamped by a special hammer to create the appropriate stratum height.Next, a layer of mica plate with a thickness not more than 1 mm was laid evenly over the stratum before the placement of the next stratum.Once completed, the surface of the model was monitored by the speckle, as shown in Figure 4.The model is laid on the model support which can provide boundary

Determination of the Model Simulation Scale
Based on similarity criteria including geometrical conditions, ambient conditions, boundary conditions, and initial conditions, and considering the geological conditions of the actual working face and the size of the model supports, the geometric similarity ratio C l of the model was determined to be 150:1.As the volume-weight of actual rocks is about 25 kN/m 3 , and the materials used to establish the model have an average volume-weight of 15 kN/m 3 , the volume-weight similarity ratio of the model was determined to be 1.667:1.Calculations indicated that C t = √ c l = 12.25:1, and that the stress similarity ratio C p = C r C l = 250:1.

Design of Model Parameters
Each model was designed to consist of nine strata including the coal seam, hard roof, immediate roof, and partially overlying strata.The model size was 2.30 m × 0.30 m × 1.23 m (L × W × H).The real height on site was 184.80 m.The rock strata that were not laid were approximately 475 m thick.According to the stress similarity ratio, a compensatory load of 0.0475 MPa was required, which was applied using pressure water bags.
The height, strength, and material mixing proportion for each rock stratum in the models were calculated based on the site test results of the actual lithology and the principle of similarity.They are shown in Table 1.
The rock masses in the models were laid from bottom to top.Before the placement of each stratum, sand, calcium carbonate, plaster, and water were mixed in a blender with the mass of each component as shown in Table 2.After mixing, the materials were spread evenly and smoothly on the model support and tamped by a special hammer to create the appropriate stratum height.Next, a layer of mica plate with a thickness not more than 1 mm was laid evenly over the stratum before the placement of the next stratum.Once completed, the surface of the model was monitored by the speckle, as shown in Figure 4.The model is laid on the model support which can provide boundary constraints, and the model can simulate the movement and deformation of the coal seam, hard roof, and overlying strata.

40%
The height of mining is 3500 mm, the height of backfill is 1800 mm, no compaction.
The height of mining is 24 mm, the height of backfill is 12 mm, initial loading stress of backfill body is 0 MPa.
The height of mining is 3500 mm, the height of backfill is 3600 mm, no compaction.
The height of mining is 24 mm, the height of backfill is 24 mm, initial loading stress of backfill body is 0 MPa.

97%
The height of mining is 3500 mm, the height of backfill is 3600 mm, compaction strength is 2 MPa.
The height of mining is 24 mm, the height of backfill is 24 mm, initial loading stress of backfill body is 8 kPa.
constraints, and the model can simulate the movement and deformation of the coal seam, hard roof, and overlying strata.

40%
The height of mining is 3500 mm, the height of backfill is 1800 mm, no compaction.
The height of mining is 24 mm, the height of backfill is 12 mm, initial loading stress of backfill body is 0 MPa.
The height of mining is 3500 mm, the height of backfill is 3600 mm, no compaction.
The height of mining is 24 mm, the height of backfill is 24 mm, initial loading stress of backfill body is 0 MPa.

97%
The height of mining is 3500 mm, the height of backfill is 3600 mm, compaction strength is 2 MPa.
The height of mining is 24 mm, the height of backfill is 24 mm, initial loading stress of backfill body is 8 kPa.

Design of Scenarios for Comparison
Four groups of physical simulation experiments with different backfilling states were performed.The mining and backfilling states of each group are shown in Table 2.   2.
Energies 2017, 10, 557 6 of 20 For convenience of description, the term "roof-controlled backfilling ratio" was devised to describe the backfilling effects [23,24].It refers to the ratio of the backfill body height to the mining height under the condition that the stress to which the backfill body is subjected is due to the weight of the main roof and the overlying soft strata.This index directly reflects the control effects of a backfill body on the roof.Specific to Face No. 6304-1, the gravity sum of the hard roof and the weak strata that deform with the hard roof was 2.27 MPa and the roof-controlled backfilling ratio in this case was referred to as the ratio of the backfill body height and the mining height when the stress endured by the backfill body was 2.27 MPa.

Backfill Mining and Monitoring of the Models
Mining commenced once the models had been laid and dried for a period of 10 days.On the actual site, each shift lasted 8 h.In each 24-hour period there were two shifts for mining and one shift for maintenance, and a footage of about 3600 mm was achieved.Based on the advance rate of two days, a cycle of mining (50 mm) and backfilling (50 mm) of the model every 200 min was established, until the stop line of the face was reached.To eliminate the boundary effect, both the start line and the stop line were 250 mm from the boundary.The mining distance was 1800 mm, for which 36 cycles of mining were required.
The non-contact strain measurement system collected the displacement data of the model every 2 s-to monitor the displacement of each rock stratum, before and after each cycle of mining and backfilling.Six pressure sensors were arranged in the coal seam to assist the strain measurement system and to monitor the changes of abutment stress ahead of the face during the mining process.The model mining and backfilling sequence and the distribution of pressure sensors are shown in Figure 5.
Energies 2017, 10, 557 6 of 20 For convenience of description, the term "roof-controlled backfilling ratio" was devised to describe the backfilling effects [23,24].It refers to the ratio of the backfill body height to the mining height under the condition that the stress to which the backfill body is subjected is due to the weight of the main roof and the overlying soft strata.This index directly reflects the control effects of a backfill body on the roof.Specific to Face No. 6304-1, the gravity sum of the hard roof and the weak strata that deform with the hard roof was 2.27 MPa and the roof-controlled backfilling ratio in this case was referred to as the ratio of the backfill body height and the mining height when the stress endured by the backfill body was 2.27 MPa.

Backfill Mining and Monitoring of the Models
Mining commenced once the models had been laid and dried for a period of 10 days.On the actual site, each shift lasted 8 h.In each 24-hour period there were two shifts for mining and one shift for maintenance, and a footage of about 3600 mm was achieved.Based on the advance rate of two days, a cycle of mining (50 mm) and backfilling (50 mm) of the model every 200 min was established, until the stop line of the face was reached.To eliminate the boundary effect, both the start line and the stop line were 250 mm from the boundary.The mining distance was 1800 mm, for which 36 cycles of mining were required.
The non-contact strain measurement system collected the displacement data of the model every 2 s-to monitor the displacement of each rock stratum, before and after each cycle of mining and backfilling.Six pressure sensors were arranged in the coal seam to assist the strain measurement system and to monitor the changes of abutment stress ahead of the face during the mining process.The model mining and backfilling sequence and the distribution of pressure sensors are shown in Figure 5.

Tests of the Deformation Behavior of the Solid Backfill Body
A solid backfill body is subjected to confined compression after the roof in the mined-out area becomes deformed.The degree of deformation of the solid backfill body varies with the stress imposed on it by the roof.The coal gangue used by the solid backfill coal mining face in the Jining No. 3 coal mine was adopted for testing purposes.The universal material testing system and specially-made compact steel cylinder shown in Figure 6 were used to test the coal gangue's deformation behavior under confined compression [11].

Tests of the Deformation Behavior of the Solid Backfill Body
A solid backfill body is subjected to confined compression after the roof in the mined-out area becomes deformed.The degree of deformation of the solid backfill body varies with the stress imposed on it by the roof.The coal gangue used by the solid backfill coal mining face in the Jining No. 3 coal mine was adopted for testing purposes.The universal material testing system and specially-made compact steel cylinder shown in Figure 6 were used to test the coal gangue's deformation behavior under confined compression [11].Two sets of tests were performed to simulate two backfilling states of the coal gangue: no tamping, and tamping by a force of 2 MPa.During the test involving no tamping force, stress was directly applied to the coal gangue inside the steel cylinder until it reached 16.5 MPa, which was the stress to which the primary rock in Face No. 6304 was subjected.The strain-stress curve of the coal gangue during the compaction process was recorded.During the test involving tamping force, a stress of 2 MPa was initially applied to the coal gangue inside the steel cylinder, before being unloaded to 0 MPa.Stress was then reapplied until it reached 16.5 MPa.The strain-stress curve of the coal gangue during this process was also recorded.

Design of Solid Backfill Analogs for Different Backfilling Scenarios
To precisely simulate the roof deformation and stress changes of the face during solid backfill coal mining, it was necessary to design a suitable solid backfill analog that had a similar strain-stress relationship to that of the actual solid backfill body.Based on past test experiences, sponges and plain paper were used as the solid backfill analog.Tests have revealed that after the sponge is precompressed, its strain-stress curve is closest to that of an actual solid backfill body.The sponge was therefore pre-compressed by veneers and Scotch tape and fixed at an appropriate height, as shown in Figure 7.In view of the backfilling heights of 1.2 cm and 2.4 cm required by the tests, different combinations of sponges and paper were designed, which resulted in 13 different variations, as shown in Table 3.
An MTS electro-hydraulic servo tester (MTS Systems Co., Eden Prairie, MN, USA) was used to perform uniaxial compression tests on each group.The strain-stress curve of each group was obtained by providing stress loading ranging from 0 kPa to 66 kPa (similar pressure stress to 16.5 MPa).Comparisons to the strain-stress curves of the coal gangue mentioned in Section 3.1 indicate that the strain-stress curves of groups 2, 7 and 10 were the closest.Therefore, those groups were used Two sets of tests were performed to simulate two backfilling states of the coal gangue: no tamping, and tamping by a force of 2 MPa.During the test involving no tamping force, stress was directly applied to the coal gangue inside the steel cylinder until it reached 16.5 MPa, which was the stress to which the primary rock in Face No. 6304 was subjected.The strain-stress curve of the coal gangue during the compaction process was recorded.During the test involving tamping force, a stress of 2 MPa was initially applied to the coal gangue inside the steel cylinder, before being unloaded to 0 MPa.Stress was then reapplied until it reached 16.5 MPa.The strain-stress curve of the coal gangue during this process was also recorded.

Design of Solid Backfill Analogs for Different Backfilling Scenarios
To precisely simulate the roof deformation and stress changes of the face during solid backfill coal mining, it was necessary to design a suitable solid backfill analog that had a similar strain-stress relationship to that of the actual solid backfill body.Based on past test experiences, sponges and plain paper were used as the solid backfill analog.Tests have revealed that after the sponge is pre-compressed, its strain-stress curve is closest to that of an actual solid backfill body.The sponge was therefore pre-compressed by veneers and Scotch tape and fixed at an appropriate height, as shown in Figure 7. Two sets of tests were performed to simulate two backfilling states of the coal gangue: no tamping, and tamping by a force of 2 MPa.During the test involving no tamping force, stress was directly applied to the coal gangue inside the steel cylinder until it reached 16.5 MPa, which was the stress to which the primary rock in Face No. 6304 was subjected.The strain-stress curve of the coal gangue during the compaction process was recorded.During the test involving tamping force, a stress of 2 MPa was initially applied to the coal gangue inside the steel cylinder, before being unloaded to 0 MPa.Stress was then reapplied until it reached 16.5 MPa.The strain-stress curve of the coal gangue during this process was also recorded.

Design of Solid Backfill Analogs for Different Backfilling Scenarios
To precisely simulate the roof deformation and stress changes of the face during solid backfill coal mining, it was necessary to design a suitable solid backfill analog that had a similar strain-stress relationship to that of the actual solid backfill body.Based on past test experiences, sponges and plain paper were used as the solid backfill analog.Tests have revealed that after the sponge is precompressed, its strain-stress curve is closest to that of an actual solid backfill body.The sponge was therefore pre-compressed by veneers and Scotch tape and fixed at an appropriate height, as shown in Figure 7.In view of the backfilling heights of 1.2 cm and 2.4 cm required by the tests, different combinations of sponges and paper were designed, which resulted in 13 different variations, as shown in Table 3.
An MTS electro-hydraulic servo tester (MTS Systems Co., Eden Prairie, MN, USA) was used to perform uniaxial compression tests on each group.The strain-stress curve of each group was obtained by providing stress loading ranging from 0 kPa to 66 kPa (similar pressure stress to 16.5 MPa).Comparisons to the strain-stress curves of the coal gangue mentioned in Section 3.1 indicate that the strain-stress curves of groups 2, 7 and 10 were the closest.Therefore, those groups were used In view of the backfilling heights of 1.2 cm and 2.4 cm required by the tests, different combinations of sponges and paper were designed, which resulted in 13 different variations, as shown in Table 3.
An MTS electro-hydraulic servo tester (MTS Systems Co., Eden Prairie, MN, USA) was used to perform uniaxial compression tests on each group.The strain-stress curve of each group was obtained by providing stress loading ranging from 0 kPa to 66 kPa (similar pressure stress to 16.5 MPa).Comparisons to the strain-stress curves of the coal gangue mentioned in Section 3.1 indicate that the strain-stress curves of groups 2, 7 and 10 were the closest.Therefore, those groups were used as solid backfill analogs for the scenarios with roof-controlled backfilling ratios of 40%, 82.5%, and 97%. Figure 8 displays the results of the strain-stress curves for comparison.as solid backfill analogs for the scenarios with roof-controlled backfilling ratios of 40%, 82.5%, and 97%. Figure 8 displays the results of the strain-stress curves for comparison.When performing the physical simulation experiments, the solid backfill analogs were made into cuboids measuring 300 mm × 50 mm (L × W).For every 50 mm of mining, a corresponding solid backfill analog was immediately placed inside the mined-out space, and the mining and backfilling continued like this in sequence.

Scenario 1: Roof-Controlled Backfilling Ratio of 0
When the roof-controlled backfilling ratio was 0, i.e., the caving method was used, the hard roof displayed obvious characteristics of initial fractures and periodic fractures.In order to directly correspond to the actual face, the data in the similarity model had to be converted from the data from the actual working face using a similarity ratio.The following is a staged description of the deformation behavior for Scenario 1: (1) In the initial phase of mining, the roof displayed little deformation.When the mining reached a distance of 60 m from the start line, the maximum displacement of the roof was only 0.11 m, and no obvious cracks were observed.(2) When the mining reached a distance of 127.5 m, obvious separation was observed between the hard roof and the No. 6 stratum-mudstone and sandstone interbedding; the maximum When performing the physical simulation experiments, the solid backfill analogs were made into cuboids measuring 300 mm × 50 mm (L × W).For every 50 mm of mining, a corresponding solid backfill analog was immediately placed inside the mined-out space, and the mining and backfilling continued like this in sequence.

Scenario 1: Roof-Controlled Backfilling Ratio of 0
When the roof-controlled backfilling ratio was 0, i.e., the caving method was used, the hard roof displayed obvious characteristics of initial fractures and periodic fractures.In order to directly correspond to the actual face, the data in the similarity model had to be converted from the data from the actual working face using a similarity ratio.The following is a staged description of the deformation behavior for Scenario 1:  (4) When the mining reached a distance of 217.5 m, the hard roof experienced a periodic fracture, with a fracture interval of 60 m.The overlying strata were affected by the mining activities; the No. 5 stratum experienced a maximum subsidence of 2.2 m, as shown in Figure 10.(5) When the mining reached a distance of 270 m (the stop line), a second periodic fracture occurred in the hard roof, with a fracture interval of 52 m.The hard roof subsided by a maximum amount of 3.42 m, and cracks were seen developing across the whole model height, as shown in Figure 11.(4) When the mining reached a distance of 217.5 m, the hard roof experienced a periodic fracture, with a fracture interval of 60 m.The overlying strata were affected by the mining activities; the No. 5 stratum experienced a maximum subsidence of 2.2 m, as shown in Figure 10.(4) When the mining reached a distance of 217.5 m, the hard roof experienced a periodic fracture, with a fracture interval of 60 m.The overlying strata were affected by the mining activities; the No. 5 stratum experienced a maximum subsidence of 2.2 m, as shown in Figure 10.(5) When the mining reached a distance of 270 m (the stop line), a second periodic fracture occurred in the hard roof, with a fracture interval of 52 m.The hard roof subsided by a maximum amount of 3.42 m, and cracks were seen developing across the whole model height, as shown in Figure 11.(5) When the mining reached a distance of 270 m (the stop line), a second periodic fracture occurred in the hard roof, with a fracture interval of 52 m.The hard roof subsided by a maximum amount of 3.42 m, and cracks were seen developing across the whole model height, as shown in Figure 11.

Scenario 2: Roof-Controlled Backfilling Ratio of 40%
When the roof-controlled backfilling ratio was 40%, i.e., the mining height was 3.5 m and the backfilling height was 1.75 m, the hard roof still displayed obvious characteristics of initial fractures and periodic fractures.The deformation behavior can be summarized by the following five stages: (1) Before the mining reached a distance of 135 m, the roof displayed some degree of deformation but did not show obvious cracks or separations.(2) When the mining reached a distance of 135 m, at a distance of 28 m from the start line, a crack about 11 m high occurred at the lower part of the roof.The crack did not penetrate through the roof.At this time, the roof subsided by a maximum amount of 0.51 m and was separated from the No. 6 stratum.(3) When the mining reached a distance of 165 m, the hard roof experienced its initial fracture.
Obvious cracks were observed above the start line, above the middle of the mined-out area, and above the working face.Due to uneven laying of the model, the fracture of the roof was somewhat stratum-wise.Above the start line and middle of the mined-out area, cracks did not penetrate through the rock mass but developed from the vertical direction toward the horizontal direction.At this time the maximum subsidence of the roof was 2.41 m, as shown in Figure 12.

Scenario 2: Roof-Controlled Backfilling Ratio of 40%
When the roof-controlled backfilling ratio was 40%, i.e., the mining height was 3.5 m and the backfilling height was 1.75 m, the hard roof still displayed obvious characteristics of initial fractures and periodic fractures.The deformation behavior can be summarized by the following five stages: (1) Before the mining reached a distance of 135 m, the roof displayed some degree of deformation but did not show obvious cracks or separations.(2) When the mining reached a distance of 135 m, at a distance of 28 m from the start line, a crack about 11 m high occurred at the lower part of the roof.The crack did not penetrate through the roof.At this time, the roof subsided by a maximum amount of 0.51 m and was separated from the No. 6 stratum.(3) When the mining reached a distance of 165 m, the hard roof experienced its initial fracture.
Obvious cracks were observed above the start line, above the middle of the mined-out area, and above the working face.Due to uneven laying of the model, the fracture of the roof was somewhat stratum-wise.Above the start line and middle of the mined-out area, cracks did not penetrate through the rock mass but developed from the vertical direction toward the horizontal direction.At this time the maximum subsidence of the roof was 2.41 m, as shown in Figure 12.  (5) When the mining reached a distance of 270 m (the stop line), a second periodic fracture occurred in the hard roof, with a fracture interval of 45 m.By this time cracks were developing across the entire model height, as shown in Figure 14.

Scenario 3: Roof-Controlled Backfilling Ratio of 82.5%
When the roof-controlled backfilling ratio was 82.5%, i.e., both the backfilling height and the mining height were 3.5 m, the hard roof fractured but did not display obvious weighting phenomena.During the entire backfill mining process, the rock mass did not show any caving zones; only cracking zones and curve subsidence zones.The deformation behavior can be summarized by the following three stages:  (3) No obvious new cracks or separations occurred until the completion of mining in the working face.Once the deformation of the hard roof and the overlying strata had become stable, the existing crack developed to a height of 63 m, with the roof subsiding by a maximum amount of 1.13 m, as shown in Figure 16.When the roof-controlled backfilling ratio was 82.5%, i.e., both the backfilling height and the mining height were 3.5 m, the hard roof fractured but did not display obvious weighting phenomena.During the entire backfill mining process, the rock mass did not show any caving zones; only cracking zones and curve subsidence zones.The deformation behavior can be summarized by the following three stages:

Scenario 3: Roof-Controlled Backfilling Ratio of 82.5%
When the roof-controlled backfilling ratio was 82.5%, i.e., both the backfilling height and the mining height were 3.5 m, the hard roof fractured but did not display obvious weighting phenomena.During the entire backfill mining process, the rock mass did not show any caving zones; only cracking zones and curve subsidence zones.The deformation behavior can be summarized by the following three stages:  (3) No obvious new cracks or separations occurred until the completion of mining in the working face.Once the deformation of the hard roof and the overlying strata had become stable, the existing crack developed to a height of 63 m, with the roof subsiding by a maximum amount of 1.13 m, as shown in Figure 16.(3) No obvious new cracks or separations occurred until the completion of mining in the working face.
Once the deformation of the hard roof and the overlying strata had become stable, the existing crack developed to a height of 63 m, with the roof subsiding by a maximum amount of 1.13 m, as shown in Figure 16.

Scenario 4: Roof-Controlled Backfilling Ratio of 97%
When the roof-controlled backfilling ratio was 97%, i.e., both the backfilling height and the mining height were 3.5 m, the backfill body was compressed by an initial stress of 2 MPa.During the modeled mining process, the hard roof did not fracture and only one vertical crack developed.The overlying strata did not show any caving zones; only cracking zones and curve subsidence zones.The deformation behavior can be summarized by the following three stages: (1) Before a distance of 232.5 m, the roof only showed curve subsidence, with no obvious separations or cracks.The maximum subsidence of the roof was 0.43 m at a distance of 150 m. (2) When the mining reached a distance of 232.5 m, at a distance of 125.5 m from the start line, a crack about 10 m high occurred at the lower part of the roof.The crack did not penetrate through the roof.At this time the roof subsided by a maximum amount of 0.77 m. (3) As the mining advanced, the crack developed higher, but no new cracks or separations developed.After the completion of mining in the working face, the existing crack grew to a maximum of 32 m but still did not penetrate through the roof.The maximum subsidence of the roof was 0.81 m, as shown in Figure 17.

Scenario 4: Roof-Controlled Backfilling Ratio of 97%
When the roof-controlled backfilling ratio was 97%, i.e., both the backfilling height and the mining height were 3.5 m, the backfill body was compressed by an initial stress of 2 MPa.During the modeled mining process, the hard roof did not fracture and only one vertical crack developed.The overlying strata did not show any caving zones; only cracking zones and curve subsidence zones.The deformation behavior can be summarized by the following three stages: (1) Before a distance of 232.5 m, the roof only showed curve subsidence, with no obvious separations or cracks.The maximum subsidence of the roof was 0.43 m at a distance of 150 m. (2) When the mining reached a distance of 232.5 m, at a distance of 125.5 m from the start line, a crack about 10 m high occurred at the lower part of the roof.The crack did not penetrate through the roof.At this time the roof subsided by a maximum amount of 0.77 m. (3) As the mining advanced, the crack developed higher, but no new cracks or separations developed.
After the completion of mining in the working face, the existing crack grew to a maximum of 32 m but still did not penetrate through the roof.The maximum subsidence of the roof was 0.81 m, as shown in Figure 17.

Effect of the Solid Backfill Body on Roof Deformation
Deformation characteristics of the hard roof were summarized, compared, and analyzed based on the four groups of experiments with different roof-controlled backfilling ratios.The results are shown in Table 4.At the bottom of the roof, 37 evenly spaced points were selected to obtain the final roof subsidence curves after the completion of mining in each scenario.These curves are displayed in Figure 18.   4. At the bottom of the roof, 37 evenly spaced points were selected to obtain the final roof subsidence curves after the completion of mining in each scenario.These curves are displayed in Figure 18.From Table 4 and Figure 18, the following becomes apparent: (1) The maximum subsidence of the hard roof decreased with increasing roof-controlled backfilling ratio.When the roof-controlled backfilling ratio was 97% in comparison with the caving method (0%), the maximum subsidence of the hard roof dropped from 3.42 m to 0.81 m; a decrease of 76.3%.As the roof-controlled backfilling ratio increased, the maximum height of the cracks that developed decreased; the maximum was only 32 m for the scenario with a backfilling ratio of 97%.
(2) When the roof-controlled backfilling ratios were 0% and 40%, the hard roof did not make contact with the backfill body prior to the occurrence of fractures, and displayed obvious initial and periodic fractures.Meanwhile, the mining distance between two fractures was only the distance of one cycle of mining.Due to the uneven laying of the rock strata and the error in the mining time, it was considered that the fracture intervals under these two scenarios were the same and that the mining distance at which the first roof crack appeared under both scenarios was also the same.(3) When the roof-controlled backfilling ratios were 82.5% and 97%, the hard roof maintained contact with the backfill body the entire time.As mining progressed along the working face, the roof displayed successional curve subsidence, but with no sudden fractures.Only one crack developed and the mining distance at which the crack appeared was 150 m and 232.5 m for ratios of 82.5% and 97%, respectively, compared to a distance of 127.5 m with the caving method.
The above results indicate that, with an increase of the backfilling ratio, the solid backfill body effectively controlled the roof's deformation, and the final subsidence of the roof decreased with increasing backfilling ratio.However, the backfilling ratio needed to reach a certain level before the roof was totally impervious to fractures.The critical factor was whether the roof contacted the backfill body prior to the occurrence of the fracture.When the backfilling ratio reached 82.5% and higher, the roof no longer displayed any initial fractures or periodic fractures.

Effect of the Solid Backfill Body on Distribution of Abutment Stress Ahead of the Face
During mining in the four scenarios, the data from the pressure sensors in the coal seam were recorded.The stable values from each mining and backfilling cycle were used to derive distribution curves of the abutment stress ahead of the face for different backfilling ratios.The results are shown in Figure 19 and in Table 5. Roof-controlled backfilling ratio 0% 40% 82.5% 97% From Table 4 and Figure 18, the following becomes apparent: (1) The maximum subsidence of the hard roof decreased with increasing roof-controlled backfilling ratio.When the roof-controlled backfilling ratio was 97% in comparison with the caving method (0%), the maximum subsidence of the hard roof dropped from 3.42 m to 0.81 m; a decrease of 76.3%.As the roof-controlled backfilling ratio increased, the maximum height of the cracks that developed decreased; the maximum was only 32 m for the scenario with a backfilling ratio of 97%.(2) When the roof-controlled backfilling ratios were 0% and 40%, the hard roof did not make contact with the backfill body prior to the occurrence of fractures, and displayed obvious initial and periodic fractures.Meanwhile, the mining distance between two fractures was only the distance of one cycle of mining.Due to the uneven laying of the rock strata and the error in the mining time, it was considered that the fracture intervals under these two scenarios were the same and that the mining distance at which the first roof crack appeared under both scenarios was also the same.(3) When the roof-controlled backfilling ratios were 82.5% and 97%, the hard roof maintained contact with the backfill body the entire time.As mining progressed along the working face, the roof displayed successional curve subsidence, but with no sudden fractures.Only one crack developed and the mining distance at which the crack appeared was 150 m and 232.5 m for ratios of 82.5% and 97%, respectively, compared to a distance of 127.5 m with the caving method.
The above results indicate that, with an increase of the backfilling ratio, the solid backfill body effectively controlled the roof's deformation, and the final subsidence of the roof decreased with increasing backfilling ratio.However, the backfilling ratio needed to reach a certain level before the roof was totally impervious to fractures.The critical factor was whether the roof contacted the backfill body prior to the occurrence of the fracture.When the backfilling ratio reached 82.5% and higher, the roof no longer displayed any initial fractures or periodic fractures.

Effect of the Solid Backfill Body on Distribution of Abutment Stress Ahead of the Face
During mining in the four scenarios, the data from the pressure sensors in the coal seam were recorded.The stable values from each mining and backfilling cycle were used to derive distribution curves of the abutment stress ahead of the face for different backfilling ratios.The results are shown in Figure 19 and in Table 5.The following observations can be derived from Figure 19 and Table 5: (1) When the backfilling ratio was 40%, the roof did not contact the backfill body prior to the occurrence of the fracture.The distribution of abutment stress ahead of the face under this scenario was basically the same as that under the caving method (backfilling ratio of 0%).The peak abutment stress ahead of the face reached approximately 43 MPa, with an impact range of about 45 m. (2) When the backfilling ratio reached 82.5% and 97%, the abutment stress ahead of the face decreased to 29 MPa and 23.5 MPa, respectively; a decrease of 32.7% and 45.5%, respectively, when compared to the caving method.The impact range also dropped significantly, reducing from 45 m to 30 m and 22.5 m with the backfilling ratios of 82.5% and 97%, respectively.

Mining and Geological Conditions of the Test Site
The test site is located to the north of the main and auxiliary shafts.It was 250 m wide by 2200 m long and three panels were planned in the pilot area.The first panel, No. 6304-1, was designed to be 80 m wide by 548 m long, with a recoverable reserve of 230,000 tons; the average panel depth was about 660 m.The uniaxial compressive strength of No. 3 coal is 11.2 MPa, and its density is 1.36 t/m 3 .The height of the coal seam is about 3.5 m [11].The geological features of the roof and floor are shown in Table 6.The following observations can be derived from Figure 19 and Table 5: (1) When the backfilling ratio was 40%, the roof did not contact the backfill body prior to the occurrence of the fracture.The distribution of abutment stress ahead of the face under this scenario was basically the same as that under the caving method (backfilling ratio of 0%).The peak abutment stress ahead of the face reached approximately 43 MPa, with an impact range of about 45 m. (2) When the backfilling ratio reached 82.5% and 97%, the abutment stress ahead of the face decreased to 29 MPa and 23.5 MPa, respectively; a decrease of 32.7% and 45.5%, respectively, when compared to the caving method.The impact range also dropped significantly, reducing from 45 m to 30 m and 22.5 m with the backfilling ratios of 82.5% and 97%, respectively.

Mining and Geological Conditions of the Test Site
The test site is located to the north of the main and auxiliary shafts.It was 250 m wide by 2200 m long and three panels were planned in the pilot area.The first panel, No. 6304-1, was designed to be 80 m wide by 548 m long, with a recoverable reserve of 230,000 tons; the average panel depth was about 660 m.The uniaxial compressive strength of No. 3 coal is 11.2 MPa, and its density is 1.36 t/m 3 .The height of the coal seam is about 3.5 m [11].The geological features of the roof and floor are shown in Table 6.

Measurement of the Ratio of Backfilled Mass to Mined Mass
In order to monitor the backfilling quality of the panel, the mass of backfilled waste rock and mined coal is monitored every day.The backfilled-to-mined coal mass ratio is used to evaluate the backfilling effect.The mass ratio of backfilled to mined coal τ is related to the roof-controlled backfilling ratio ϕ k by: (1) where g C = 16.46 kN/m 3 is the bulk density of the backfilled body obtained under a pressure of 2.27 MPa and c C = 13.33 kN/m 3 is the bulk density of the coal.
Based on monitoring data of the amount of backfilled waste rock and coal mined during the period from 15 March 2011 to 15 September 2011, the mass ratio of backfilled to mined coal of Panel No. 6304-1 is shown in Figure 19.
Figure 20 shows that the mass ratio was mostly greater than the design value of 0.93, with an average of 1.19.Through back-calculation, the average roof-controlled backfilling ratio φk was actually 96.4%.This result is close to the fourth scheme (97%) in the physical similarity simulation [11].

Measurement of Roof Subsidence
During mining, a roof deformation monitoring device was installed at the rear end of the shield support.It was connected via cables to a data acquisition instrument near the take-off location of the panel, so that changes in roof subsidence were monitored in real time.Three monitoring stations were deployed at the locations of 10, 30 and 60 m from the cutting hole.Each station had five monitoring devices uniformly distributed in the backfilling area.Taking the measurements at the 60 m location for instance, the largest roof subsidence was obtained, as shown in Figure 21 [11].

Measurement of Roof Subsidence
During mining, a roof deformation monitoring device was installed at the rear end of the shield support.It was connected via cables to a data acquisition instrument near the take-off location of the panel, so that changes in roof subsidence were monitored in real time.Three monitoring stations were deployed at the locations of 10, 30 and 60 m from the cutting hole.Each station had five monitoring devices uniformly distributed in the backfilling area.Taking the measurements at the 60 m location for instance, the largest roof subsidence was obtained, as shown in Figure 21 [11].It can be known that the maximum subsidence of hard roof is 197.8 mm, and the roof subsidence is only 5.7% of the mining height (3.5 m).If the hard roof is broken, the subsidence value should be much more than this measured value.It can be concluded that the roof has not been broken at the position of the 60 m distance from the open-off cut because of the supporting function of the backfill body.As a result, the physical similarity simulations are verified.

Measurement of Advanced Abutment Pressure
Three monitoring stations were deployed in both the maingate and tailgate of Panel No. 6304-1.At each monitoring point, four different drill hole depths (i.e., 3 m, 5 m, 10 m, and 15 m) were designed with a total of 24 borehole stress meters in total.Taking the measurements of the station 120 m from the setup room as an example, the front abutment pressure in the coal block at different depths is shown in Figure 22 [11].Figure 22 shows that the backfilled body restricted the deformation of the hard roof, mitigating the effects of mining.The stress concentration factor of the front abutment pressure was 1.44 and the influence distance was about 28 m.This result is close to the fourth scheme (97%) in the physical similarity simulation.

Conclusions
The following conclusions can be drawn from this study: (1) By using similarity theory and test results of the mechanical properties of backfill materials, physical simulation experiments were performed for four scenarios of varying roof-controlled backfilling ratios (0%, 40%, 82.5% and 97%) to study the deformation behavior of the hard roof.A combination of pre-compressed sponges and paper was used as a backfill body analogue, It can be known that the maximum subsidence of hard roof is 197.8 mm, and the roof subsidence is only 5.7% of the mining height (3.5 m).If the hard roof is broken, the subsidence value should be much more than this measured value.It can be concluded that the roof has not been broken at the position of the 60 m distance from the open-off cut because of the supporting function of the backfill body.As a result, the physical similarity simulations are verified.

Measurement of Advanced Abutment Pressure
Three monitoring stations were deployed in both the maingate and tailgate of Panel No. 6304-1.At each monitoring point, four different drill hole depths (i.e., 3 m, 5 m, 10 m, and 15 m) were designed with a total of 24 borehole stress meters in total.Taking the measurements of the station 120 m from the setup room as an example, the front abutment pressure in the coal block at different depths is shown in Figure 22 [11].It can be known that the maximum subsidence of hard roof is 197.8 mm, and the roof subsidence is only 5.7% of the mining height (3.5 m).If the hard roof is broken, the subsidence value should be much more than this measured value.It can be concluded that the roof has not been broken at the position of the 60 m distance from the open-off cut because of the supporting function of the backfill body.As a result, the physical similarity simulations are verified.

Measurement of Advanced Abutment Pressure
Three monitoring stations were deployed in both the maingate and tailgate of Panel No. 6304-1.At each monitoring point, four different drill hole depths (i.e., 3 m, 5 m, 10 m, and 15 m) were designed with a total of 24 borehole stress meters in total.Taking the measurements of the station 120 m from the setup room as an example, the front abutment pressure in the coal block at different depths is shown in Figure 22 [11].Figure 22 shows that the backfilled body restricted the deformation of the hard roof, mitigating the effects of mining.The stress concentration factor of the front abutment pressure was 1.44 and the influence distance was about 28 m.This result is close to the fourth scheme (97%) in the physical similarity simulation.

Conclusions
The following conclusions can be drawn from this study: (1) By using similarity theory and test results of the mechanical properties of backfill materials, physical simulation experiments were performed for four scenarios of varying roof-controlled backfilling ratios (0%, 40%, 82.5% and 97%) to study the deformation behavior of the hard roof.A combination of pre-compressed sponges and paper was used as a backfill body analogue, Figure 22 shows that the backfilled body restricted the deformation of the hard roof, mitigating the effects of mining.The stress concentration factor of the front abutment pressure was 1.44 and the influence distance was about 28 m.This result is close to the fourth scheme (97%) in the physical similarity simulation.

Conclusions
The following conclusions can be drawn from this study: (1) By using similarity theory and test results of the mechanical properties of backfill materials, physical simulation experiments were performed for four scenarios of varying roof-controlled backfilling ratios (0%, 40%, 82.5% and 97%) to study the deformation behavior of the hard roof.A combination of pre-compressed sponges and paper was used as a backfill body analogue, which precisely simulated the deformation process that the backfill body experienced under the stress of the roof for each backfilling scenario.

Figure 1 .
Figure 1.Backfilling hydraulic support at the working face.

Figure 1 .
Figure 1.Backfilling hydraulic support at the working face.

2. 2 .
Basic Parameters of the Working Face Backfill coal mining technology has been applied in Panel No. 6304-1 of the Jisan coal mine, where the No. 3 coal seam of the Shanxi formation is the main mineable coal seam.The average thickness of the No. 3 coal seam is 3.5 m, its average dip is 2.5 • , its density is 1.36 t/m 3 , and its average depth is 650 m.It has no immediate roof.The main roof is a fine sandstone of 32.5-49.75m Energies 2017, 10, 557 3 of 20 thick (with an average thickness of 41.6 m), the uniaxial compressive strength is 84.4 MPa, the tensile strength is 13.5 MPa, and the elastic modulus is 17.0 GPa.It can be regarded as a hard roof.The panel, No. 6304, was divided into three working faces, with a face length of 80 m and advancing distance of 518 m.The physical simulation experiment was mainly based on the specific conditions of Panel No. 6304-1, whose working face layout is shown in Figure 2. Energies 2017, 10, 557 3 of 20 was divided into three working faces, with a face length of 80 m and advancing distance of 518 m.The physical simulation experiment was mainly based on the specific conditions of Panel No. 6304-1, whose working face layout is shown in Figure 2.

Figure 3 .
Figure 3. Non-contact strain measurement system.(a) Monitoring equipment and its placement; (b) Data processing system.2.3.2.Determination of the Model Simulation ScaleBased on similarity criteria including geometrical conditions, ambient conditions, boundary conditions, and initial conditions, and considering the geological conditions of the actual working face and the size of the model supports, the geometric similarity ratio Cl of the model was determined to be 150:1.As the volume-weight of actual rocks is about 25 kN/m 3 , and the materials used to establish the model have an average volume-weight of 15 kN/m 3 , the volume-weight similarity ratio of the model was determined to be 1.667:1.Calculations indicated that Ct =

Figure 3 .
Figure 3. Non-contact strain measurement system.(a) Monitoring equipment and its placement; (b) Data processing system.

2. 5 .
Backfill Mining of the Models and Design of Monitoring Scheme 2.5.1.Design of Scenarios for Comparison Four groups of physical simulation experiments with different backfilling states were performed.The mining and backfilling states of each group are shown in Table

Figure 5 .
Figure 5.The model mining and backfilling sequence and the distribution of pressure sensors.

Figure 5 .
Figure 5.The model mining and backfilling sequence and the distribution of pressure sensors.

Figure 6 .
Figure 6.The universal material testing system and specially-made compact steel cylinder.

Figure 7 .
Figure 7. Diagram and field photo of the material simulating the backfill body.

Figure 6 .
Figure 6.The universal material testing system and specially-made compact steel cylinder.

Figure 6 .
Figure 6.The universal material testing system and specially-made compact steel cylinder.

Figure 7 .
Figure 7. Diagram and field photo of the material simulating the backfill body.

Figure 7 .
Figure 7. Diagram and field photo of the material simulating the backfill body.

Figure 8 .
Figure 8. Comparisons of strain-stress curves between similar materials and backfilling bodies.

Figure 8 .
Figure 8. Comparisons of strain-stress curves between similar materials and backfilling bodies.

( 1 ) 9 .
In the initial phase of mining, the roof displayed little deformation.When the mining reached a distance of 60 m from the start line, the maximum displacement of the roof was only 0.11 m, and no obvious cracks were observed.(2)When the mining reached a distance of 127.5 m, obvious separation was observed between the hard roof and the No. 6 stratum-mudstone and sandstone interbedding; the maximum displacement in the vertical direction was 0.29 m.At a distance of 67 m from the start line, a crack about 14 m high occurred at the lower part of the hard roof.The crack did not penetrate through the roof.(3) When the mining reached a distance of 157.5 m, the hard roof experienced its initial fracture.Penetrating cracks occurred above the start line, above the middle of the mined-out area, and above the working face.The roof subsided by a maximum amount of 2.91 m.At this time, the No. 6 stratum subsided and was separated from the No. 5 stratum above, as shown in Figure Energies 2017, 10, 557 9 of 20 displacement in the vertical direction was 0.29 m.At a distance of 67 m from the start line, a crack about 14 m high occurred at the lower part of the hard roof.The crack did not penetrate through the roof.(3) When the mining reached a distance of 157.5 m, the hard roof experienced its initial fracture.Penetrating cracks occurred above the start line, above the middle of the mined-out area, and above the working face.The roof subsided by a maximum amount of 2.91 m.At this time, the No. 6 stratum subsided and was separated from the No. 5 stratum above, as shown in Figure 9.

Figure 9 .
Figure 9. Deformation and displacement distribution cloud graph of overlying strata when the mining had advanced to 157.5 m with a roof-controlled backfilling ratio of 0 (caving method).

Figure 10 .
Figure 10.Deformation and displacement distribution cloud graph of overlying strata when the mining had advanced to 217.5 m with a roof-controlled backfilling ratio of 0 (caving method).

Figure 9 .
Figure 9. Deformation and displacement distribution cloud graph of overlying strata when the mining had advanced to 157.5 m with a roof-controlled backfilling ratio of 0 (caving method).
vertical direction was 0.29 m.At a distance of 67 m from the start line, a crack about 14 m high occurred at the lower part of the hard roof.The crack did not penetrate through the roof.(3)When the mining reached a distance of 157.5 m, the hard roof experienced its initial fracture.Penetrating cracks occurred above the start line, above the middle of the mined-out area, and above the working face.The roof subsided by a maximum amount of 2.91 m.At this time, the No. 6 stratum subsided and was separated from the No. 5 stratum above, as shown in Figure9.

Figure 9 .
Figure 9. Deformation and displacement distribution cloud graph of overlying strata when the mining had advanced to 157.5 m with a roof-controlled backfilling ratio of 0 (caving method).

Figure 10 .
Figure 10.Deformation and displacement distribution cloud graph of overlying strata when the mining had advanced to 217.5 m with a roof-controlled backfilling ratio of 0 (caving method).

Figure 10 .
Figure 10.Deformation and displacement distribution cloud graph of overlying strata when the mining had advanced to 217.5 m with a roof-controlled backfilling ratio of 0 (caving method).

Figure 11 .
Figure 11.Deformation and displacement distribution cloud graph of overlying strata when the mining had advanced to 270 m with a roof-controlled backfilling ratio of 0 (caving method).(a) Deformation of overlying strata; (b) Displacement distribution cloud graph of overlying strata.

Figure 11 .
Figure 11.Deformation and displacement distribution cloud graph of overlying strata when the mining had advanced to 270 m with a roof-controlled backfilling ratio of 0 (caving method).(a) Deformation of overlying strata; (b) Displacement distribution cloud graph of overlying strata.

Figure 12 .( 4 )
Figure 12.Deformation and displacement distribution cloud graph of overlying strata when the mining had advanced to 165 m with a roof-controlled backfilling ratio of 40%.(4) When the mining reached a distance of 225 m, the hard roof experienced a periodic fracture, with a fracture interval of 60 m.The fractured roof weighed completely on the backfill body, and penetrating cracks occurred above the working face.At this time, the separations and cracks above the middle of the mined-out area gradually closed.The No. 6 stratum displayed vertical penetrating cracks and was separated from the No. 5 stratum, which also subsided by a maximum amount of 1.5 m, as shown in Figure 13.

Figure 13 .
Figure 13.Deformation and displacement distribution cloud graph of overlying strata when the mining had advanced to 225 m with a roof-controlled backfilling ratio of 40%.

Figure 12 .( 4 ) 20 Figure 12 .
Figure 12.Deformation and displacement distribution cloud graph of overlying strata when the mining had advanced to 165 m with a roof-controlled backfilling ratio of 40%.

Figure 13 .Figure 13 . 20 Figure 12 .( 4 )
Figure 13.Deformation and displacement distribution cloud graph of overlying strata when the mining had advanced to 225 m with a roof-controlled backfilling ratio of 40%.(5)When the mining reached a distance of 270 m (the stop line), a second periodic fracture occurred in the hard roof, with a fracture interval of 45 m.By this time cracks were developing across the entire model height, as shown in Figure14.

Figure 13 .( 5 )Figure 14 .Figure 14 .
Figure 13.Deformation and displacement distribution cloud graph of overlying strata when the mining had advanced to 225 m with a roof-controlled backfilling ratio of 40%.(5) When the mining reached a distance of 270 m (the stop line), a second periodic fracture occurred in the hard roof, with a fracture interval of 45 m.By this time cracks were developing across the entire model height, as shown in Figure 14.

( 1 )
Before a distance of 150 m, the roof displayed successional curve subsidence with mining in the face.The maximum subsidence of the roof was only 0.13 m at a distance of 75 m.(2) When the mining reached a distance of 150 m, at a distance of 67 m from the start line, a crack about 13 m high occurred at the lower part of the roof.At this time the roof subsided by a maximum amount of 0.53 m, as shown in Figure 15.

Figure 15 .
Figure 15.Deformation and displacement distribution cloud graph of overlying strata when the mining had advanced to 150 m with a roof-controlled backfilling ratio of 82.5%.

Figure 14 .
Figure 14.Deformation and displacement distribution cloud graph of overlying strata when the mining had advanced to 270 m with a roof-controlled backfilling ratio of 40%.(a) Deformation of overlying strata; (b) Displacement distribution cloud graph of overlying strata.

( 1 )Figure 14 .
Figure 14.Deformation and displacement distribution cloud graph of overlying strata when the mining had advanced to 270 m with a roof-controlled backfilling ratio of 40%.(a) Deformation of overlying strata; (b) Displacement distribution cloud graph of overlying strata.

( 1 )
Before a distance of 150 m, the roof displayed successional curve subsidence with mining in the face.The maximum subsidence of the roof was only 0.13 m at a distance of 75 m.(2) When the mining reached a distance of 150 m, at a distance of 67 m from the start line, a crack about 13 m high occurred at the lower part of the roof.At this time the roof subsided by a maximum amount of 0.53 m, as shown in Figure 15.

Figure 15 .
Figure 15.Deformation and displacement distribution cloud graph of overlying strata when the mining had advanced to 150 m with a roof-controlled backfilling ratio of 82.5%.

Figure 15 .
Figure 15.Deformation and displacement distribution cloud graph of overlying strata when the mining had advanced to 150 m with a roof-controlled backfilling ratio of 82.5%.

Figure 16 .
Figure 16.Deformation and displacement distribution cloud graph of overlying strata when the mining had advanced to 270 m with a roof-controlled backfilling ratio of 82.5%.(a) Deformation of overlying strata; (b) Displacement distribution cloud graph of overlying strata.

Figure 16 .
Figure 16.Deformation and displacement distribution cloud graph of overlying strata when the mining had advanced to 270 m with a roof-controlled backfilling ratio of 82.5%.(a) Deformation of overlying strata; (b) Displacement distribution cloud graph of overlying strata.

Figure 17 .
Figure 17.Deformation and displacement distribution cloud graph of overlying strata when the mining had advanced to 270 m with a roof-controlled backfilling ratio of 97%.(a) Deformation of overlying strata; (b) Displacement distribution cloud graph of overlying strata.

Figure 17 . 5 . 5 . 1 .
Figure 17.Deformation and displacement distribution cloud graph of overlying strata when the mining had advanced to 270 m with a roof-controlled backfilling ratio of 97%.(a) Deformation of overlying strata; (b) Displacement distribution cloud graph of overlying strata.

Figure 20 .
Figure 20.Variation of the mass ratio of backfilled waste to mined coal.

Figure 20 .
Figure 20.Variation of the mass ratio of backfilled waste to mined coal.

Figure 21 .
Figure 21.Measured roof subsidence in the backfilling area.

Figure 22 .
Figure 22.Measured abutment pressures at the 120 m station.

Figure 21 .
Figure 21.Measured roof subsidence in the backfilling area.

Figure 22 .
Figure 22.Measured abutment pressures at the 120 m station.

Figure 22 .
Figure 22.Measured abutment pressures at the 120 m station.

Table 1 .
Thickness and mechanical parameters of coal and rock mass.

Table 2 .
Backfill mining status of each simulation program.

Table 1 .
Thickness and mechanical parameters of coal and rock mass.

Table 2 .
Backfill mining status of each simulation program.

Table 3 .
Various combinations of sponge and paper.

Table 3 .
Various combinations of sponge and paper.

Table 4 .
Roof deformation comparison for different roof-controlled backfilling rates.

Table 4 .
Roof deformation comparison for different roof-controlled backfilling rates.

Table 5 .
Comparison of abutment stress distribution ahead of the face for different roof-controlled backfilling ratios.

Table 5 .
Comparison of abutment stress distribution ahead of the face for different roof-controlled backfilling ratios.

Table 6 .
No. 3 coal seam roof and floor.