Overlying Strata Movement Laws Induced by Longwall Mining of Deep Buried Coal Seam with Superhigh-Water Material Backfilling Technology

School of Mines, State Key Laboratory of Coal Resources and Safe Mining, Key Laboratory of Deep Coal Resource Mining, Ministry of Education of China, China University of Mining and Technology, Xuzhou 221116, China Key Laboratory of Mine Geological Hazards Mechanism and Control, Xi’an 710054, China Shandong Yineng Coal Mine Co., Ltd., Jining 272511, China School of Environment Science and Spatial Informatics, China University of Mining and Technology, Xuzhou 221116, China


Introduction
e " ree Under Mining" (mining under the buildings, the railways, and the waters) coal reserves have been reached 13.79 billion tons in China, and the coal resources under the buildings are over 8.76 billion tons.
e coal mines located in the central and eastern China mining areas are mostly accompanied by the congregate village buildings.As a result, the mining of the " ree Under Mining" coal resources is restricted by the protection of surface ecological environment [1,2].Additionally, with the tens or almost one hundred years of exploitation, many central and eastern coal mines have entered into deep buried coal seams mining (elevation over 800 m), which leads to a series issues, such as high geo-stress, high water pressure, high ground temperature, and high gas pressure [3][4][5][6].Many in situ experiments show that longwall mining with backfilling technology is an efficient way to prevent the underground disasters and protect the surface ecological environment [7].However, the field application effects and reliabilities are affected by many factors, including backfilling strength and backfilling rates.
Domestic and overseas scholars have carried out a series of studies on the backfilling technology and overlying strata movement laws.Mkadmi et al. [8] used the SIGMA/W finite element code to simulate the characteristics of the rock-backfill interface influence on the magnitude of the stresses in backfilled stopes.Mohamed et al. [9] simulated the influence of filling gables on the long-term stability of gypsum cylinders in order to prevent the collapse of surface buildings and infrastructure.Kostecki and Spearing [10] used the FLAC 3D software to analyze the plastic flow characteristics of coal pillars and compared the shear strength, tensile strength, and the influence of mine space constraints on the strength of coal pillars.In China, hydraulic backfilling, solid waste backfilling, paste-like backfilling, and superhigh-water material (SHWM) backfilling have been employed in many underground coal mines.e field observations indicate that (1) the hydraulic backfilling technology has a good surface subsidence control effect, but it has no general applications due to complicated system and high cost [11].(2) e backfilling rate of the solid waste backfilling technology has increased to 85%, which benefits long-term stability for surface subsidence.However, it cannot meet the requirement of large area longwall mining due to high investment, lack of solid waste material [12].(3) Paste-like backfilling technology has high security, high mining rate, and good surface subsidence control effect, and is environment friendly, which also has many disadvantages, such as high initial investment and high production cost [13].(4) Superhigh-water material backfilling technology has many prominent advantages besides low initial investment, low production cost, simple backfilling process, easy operation, and good backfilling material mechanics, which are of great help to high-efficient mining and effective surface subsidence control [14,15].For the research of overlying movement laws, Li et al. [16] applied the elastic foundation thin plate theory to presenting that the elastic foundation coefficient is the key factor affecting the movement and deformation of overburden during backfilling mining.Jia et al. [17] carried out a physical simulation experiment of superhigh-water material (SHWM) backfilling technology, which illustrated the laws of overburden migration and surface subsidence.Yang et al. [18] analyzed the influence of different backfilling strength and backfilling rates on the overburden movement and deformation.
e research studies mentioned above can shed light on revealing the overlying strata movement laws of longwall mining in deep buried coal seams.
Based on the typical geological and production conditions of AB Coal Mine, the theoretical analysis, numerical simulation, and in situ measurement were comprehensively applied to uncover the overlying strata movement laws of longwall mining in deep buried coal seams with SHWM backfilling technology.e reasonable parameters, including backfilling rates and water volume fraction, are determined to guarantee the safe and efficient mining of deep buried coal resources and protection of the surface ecological environment.

Engineering Background
AB Coal Mine is located in Shandong Province, an eastern China mining area.e No. 3 coal seam, mined in the CG1302 working face, has an average thickness of 3.0 m. e geological structure is simple in the longwall panel, the advancing length of the CG1302 working face is 1030 m, and the inclination length is 110 m. e ground elevation is about +47 m; the underground elevation is about −776.7 m, and thus, the buried depth is over 820 m, which means high geo-stress, high water pressure, high ground temperature, and high gas pressure increasingly emerge during longwall mining.e coal seam comprehensive histogram and UDEC numerical simulation model are shown in Figure 1. e mechanical parameters of the rock layer are illustrated in Table 1.
More than 1790 village buildings are concentrated on the surface of AB Coal Mine, which is a typical " ree Under Mining."To achieve the goal of the safe and efficient mining of deep coal resources and the protection of the surface ecological environment, the longwall mining with SHWM backfilling technology was applied in the CG1302 working face.When the working face advances, the backfilling package will be hung up and injected with mixed SHWM in the goaf immediately behind the hydraulic support.After the mixed SHWM solidified, the mined out roof will be supported by the SHWM backfilling and the left coal pillars.e backfilling technology procedure and field application are shown in Figure 2.

Theoretical Analysis
e roof, hanging in the goaf after the working face advancing, has a load composed by both self-weight and the overlying strata load.e interaction between the overlying strata forms a composite beam structure, according to the combined beam theory [19], and the overburden load is where (q n ) 1 is the load that acts on the first layer with the consideration of the n th layer, MPa; E n is the elastic modulus of the n th layer, MPa; H n is the thickness of the n th layer, m; and c n is the bulk density of the n th layer, MN/m 3 .According to the combined beam theory, the load that acts on the immediate roof is calculated to be 0.141 MPa, and the load that acts on the main roof is 0.289 MPa.Before the first fracturing of the roof, the coal wall behind the open-off cut and in front of the working face forms a fixed beam structure, as shown in Figure 3.
According to the elastic mechanics theory, the maximum shear stress emerges in the center of the end-fixed beam.According to the maximum tensile stress criterion of the rock, the length before the fracture is [20] When the roof is fractured, the safety factor is n, and the length before the fracture is 2

Advances in Civil Engineering
where L s is the roof stability length, m; [σ] is the roof tensile strength, MPa; n is the safety factor; and q is the load acts on the roof, MPa.
When the roof beam fracturing, the maximum bending deflection equation for the center section is where ω max is the maximum bending deflection; I is the moment of inertia, m 4 ; I � bh 3 /12, b is the beam width, m and h is the beam height, m; and l is the length of the hanging roof, m.
After the first fracture of the roof strata, the beam structure is in a state of one end fixed and another end hanged, forming a cantilever beam structure.e safe length of the cantilever beam before fracture is [21] where l is the period weighting length of roof strata, m.According to the mechanical parameters of coal and rock strata of the CG1302 working face, most of the rock strata have large thickness and damage due to brittle failure.erefore, the safety factor n can be determined to be 1.2.As a result, the first roof weighting lengths of the immediate roof and the main roof are 29.7 m and 63.0 m, and the period roof weighting lengths are 9.5 m and 20.3 m, respectively.
In addition, with the consideration of the maximum vertical bending values of the immediate roof and the main roof are 0.26 m and 0.48 m, and the residual expansion coefficient of the rock is 1.15 and the maximum bending value of the main roof can reach 0.48 m and become fracture while the backfilling height reaches 1.53 m. e critical backfilling rate to prevent main roof fracturing is 51.0%, and the immediate roof will fracture if the backfilling rate is less than 91.3%.

Establishment of the Numerical Simulation
Model. e software of UDEC (Universal Distinct Element Code) was applied to analyzing the overlying strata movement laws.As shown in Figure 1, the numerical model has a length and a height of 240 m and 120 m, respectively.According to the geological conditions, the upper boundary stress is set to 17.86 MPa to simulate the weight of the overlying strata.e    To predict the surface subsidence under di erent lling rates, a survey line, used to monitor the vertical displacement, was set on the top of the UDEC numerical simulation model.As shown in Figure 5, in the condition of nonback lling, the maximum vertical displacement was above 0.95 m.With the increase of lling rate from 40% to 90%, the vertical displacement decreased from 0.76 m to 0.22 m.It can be predicted that the surface subsidence will be e ectively  2) When zones rst emerge in the overlying strata around the working face, the immediate roof, the main roof, and the overlying strata start to bend and fracture in proper order.Additionally, the plastic failure zone formed a shape similar to "reverse funnel" above the goaf.(3) In the condition of back lling, the plastic zones have no signi cant change when the working face advances from 35 m to 70 m; the plastic zones only occurred within the immediate roof, and the plastic damage degree was not large enough to break the main roof.

Vertical Stress Evolutions Mechanism.
As shown in Figure 7, (1) a "butter y wings" shape vertical stress accumulation area was formed both in front of the working face and behind the open-o cut, and the stress-relaxed area was  MPa when the lling rates were 90%, 70%, 40%, and 0, respectively.As the lling rate increases, the peak value of the vertical stress in front of the coal wall decreases gradually.It can reveal that the better control of both coal wall splitting and roadway surrounding rock deformation would be achieved through increasing the lling rate.As shown in Figure 8(c), the peak values of the vertical stress were 39.7 MPa, 43.6 MPa, 45.3 MPa, and 49.7 MPa when the water volume fractions were 91%, 93%, 95%, and 97%, respectively.Simultaneously, there were only low peak stress changes while the water volume fraction varies from 93% to 95%.Considering the superhigh-water material costs and the pipeline delivery e ciency, the proper water volume fraction was determined to be 95%.

Ground Pressure Measurements.
e lling rate of 90% and the water volume fraction of 95% were used in the exploitation of the CG1302 working face.During the longwall mining with SHWM back lling technology, a real-time digital ground pressure monitor was set up in the front column of the hydraulic support, which can capture the changes of the load acted on the hydraulic support.For example, the support resistance observations of the 53 # hydraulic support are shown in Figure 9.According to the eld measurement data, the average and maximum working resistances of the hydraulic support were 19.7 MPa and 36.0MPa, and the roof weighting criterion was determined to 28.1 MPa. Figure 9 also indicates that the average period roof weighting length of the immediate roof is 9.6 m, which is verifying the reasonable of the theoretical analysis (9.5 m).With the application of the back lling technology with a lling rate of 90%, only the immediate roof has been fractured and caved into the goaf, and the dynamic load coe cient of the working face was low enough for the safe mining of the deep buried coal seam.

Surface Subsidence Measurements.
To evaluate the e ect of SHWM back lling technology for controlling the surface subsidence, a total of 8 observation lines were arranged in the corresponding ground positions of CG1301 and CG1302.ere are 17 measuring points in the survey lines, from C1 to C17, respectively.

Advances in Civil Engineering
As shown in Figure 10, the maximum cumulative dynamic subsidence of the C line was less than 200 mm, which is consistent with the numerical simulation results.e surface subsidence measurements results indicate that the village building can be kept in a stable state even a ected by the underground longwall mining-induced overlying movements.It can be concluded that the goal of the safe and e cient mining of deep buried coal resources and protection of surface ecological environment would be realized for the " ree Under Mining" coal resources with the application of SHWM back lling technology and reasonable lling parameters design.Advances in Civil Engineering

Conclusions
(1) High geo-stress, high water pressure, high ground temperature, and high gas pressure increasingly emerge during longwall mining in deep buried coal seams.Longwall mining with superhigh-water material back lling technology is one of the high-e cient ways to prevent the underground disasters and protect the surface ecological environment.(2) e software of UDEC was applied to analyze the overlying strata movement laws.In the nonback lling conditions, the large concentrated vertical stress means a huge elastic deformation energy was involved around the working face, which was likely to cause rock burst accidents and threaten the miners' life.In contrast, the back lling played a signi cant e ect to decrease the stress concentrate rate around the goaf and prevent the rock burst disasters.(3) With the comparisons of the vertical displacement movements and the vertical stress distributions, the reasonable parameters were determined, such as the lling rate of 90% and the water volume fraction of 95%, which were carried out in the eld experiments.
(4) According to the eld measurement data, the underground dynamic load was low enough to enable the safe mining, and the village building can be kept in a stable state with the application of the back lling technology.e aim of the safe and e cient mining of deep buried coal resources and protection of surface ecological environment would be achieved for exploiting the " ree Under Mining" coal resources.
Data Availability e research data used to support the ndings of this study are available from the corresponding author upon request.
simulated longwall mining face advancing from 50 to 190 m with a total mining length of 140 m.e ZC7000/19/40 type hydraulic support is used in the CG1302 working face, the supporting strength is 0.67 MPa, and the support length is 4 m.e simulation contents include the vertical displacement movement laws, the plastic zone development characteristics, and the vertical stress evolutions mechanism, which vary with di erent roof manage modes, including roof caving and back lling.

4. 2 .
Numerical Simulation Results 4.2.1.Vertical Displacement Movement Laws.As shown in Figure 4, (1) the immediate roof is bended when the working face advances 30 m, and then fractured and caved into goaf when the working face advances 35 m. (2) e main roof is fractured when the working face advances 60-70 m. (3) After the back lling in the goaf, the immediate roof is not fractured when the working face advances 35 m and 70 m, and the vertical displacement is less than 0.3 m.

Figure 8 :
Figure 8: Concentrated vertical stress curves vary with excavation length and lling rate.(a) Di erent excavation length.(b) Di erent lling rate.(c) Di erent volume fraction of water.

Figure 9 :Figure 10 :
Figure 9: Support resistance change curves during the advancing of the back lling face.

Table 1 :
Mechanical parameters of coal seam and overlying strata.