Research on the Roof Advanced Breaking Position and Inﬂuences of Large Mining Height Working Face in Shallow Coal Seam

: The large mining height (LMH) in shallow coal seam has been widely applied in the Shenfu coalﬁeld, China. The dynamic load is obvious, and the rib spalling is serious when the LMH working face concerns roof weighting. The advanced breaking position of the roof a ﬀ ects the strength of the ground pressure when the roof is broken. Firstly, based on a large number of actual measurements and physical simulation experiments, the rock formation in the fall zone, where the articulated structure cannot be articulated between the coal seam and the main roof, is called the equivalent immediate roof (EIR). When the mining height increases, the thickness of the EIR increases non-linearly. Next, based on the theory of “elastic foundation beam”, a mechanical model for the advanced breaking of the roof is established in shallow coal seam, and the calculation equation for the advanced breaking position of the roof is given; then, designed and carry out boreholes of the no. 22201 working face in the Zhangjiamao Coal Mine. The theoretical calculation of key strata results (5.6–6.9 m) are in the range of ﬁeld measurement results (5–8 m). According to the ﬁeld measurement results, the roof movement of the LMH working face is ahead of the roof weighting. Finally, we deﬁne the thickness of EIR and the mining height ratio as the immediate mining ratio k i , which a ﬀ ects the degree of ﬁlling of the goaf and determines the structural form of the main roof. When the k i is small, the goaf is fully ﬁlled; when the k i is large, the goaf is fully ﬁlled. Under the same conditions, di ﬀ erent ﬁlling rate conditions will form di ﬀ erent roof structures. Results of this research can be helpful to control roof weighting and provide early warning of possible safety problems related to the LMH working face in shallow coal seam.


Introduction
Shenfu-Dongsheng coalfield in China is one of the eight largest coalfields in the world. The main shallow coal seam in this area is Jurassic in age. The thickness of coal seam mining is 3.5-8 m. Most of the working faces adopt the mining method of LMH and full mining height [1,2]. The roof structure will show the different strength of ground pressures before and after the fracture [3][4][5][6][7]. The roof fracture intervals depend on different boundary conditions. The fracture position and rotation angle of the main roof will affect the stability of the roof [8,9]. The structural model of roof breaking was established in shallow coal seam [10,11]. Meanwhile, as the mining height increases, mining directly affects the ground surface [12]. During the period of incoming weighting, there are obvious steps sinking, dynamic load, and serious phenomenon of rib spalling [13][14][15][16][17][18]. Based on a large number of field measurement analyses, the evolution characteristics of the roof structure at 4-7 m LMH working faces are analyzed. This, combined with the physical simulation and numerical calculation, the mechanism of the upward movement of the hinge point of the roof structure (caused by the increase of the roof thickness) and the definition of the equivalent immediate roof (EIR) of the large mining height working face is proposed [19][20][21]. The immediate roof of the high mining height is divided into three types, and the roof control mechanical model of the stope is given [22,23]. Based on the characteristics of the underground pressure at no. 22303 working face in the Bulianta coal mine, a structure model of the main roof cantilever beam with a 7.0 m LMH was proposed [24][25][26]. A short cantilever beam-hinged rock beam structure of a fully mechanized mining roof with a large mining height is proposed to determine the working resistance of the support [27]. Other scholars have also studied the roof structure of the LMH stope [28][29][30][31][32]. However, there are few studies on the roof advanced breaking position of LMH working face in shallow coal seam.
In this paper, the research background concerns the 22201 working face in the Zhangjiamao Coal Mine in China, adopting a combination of theoretical analyses and field measurement methods to research on the advanced breaking position of the shallow-buried and high-mining working face. The present research is aimed at providing the basis for stable roof structure and optimize the choice of stope support and roadway advanced support methods.

Proposition of EIR
Based on a large number of field measurements and physical simulation experiments of shallow buried and LMH working faces, authors [19][20][21] found that, as the mining height increases, some of the overlying strata that originally belonged to the geological main roof also caved. This phenomenon presents the function of immediate roof. Therefore, we propose that all rock layers that do not form a hinged structure between the main roof hinge layer and the coal seam are called "EIR". During normal mining, the geological immediate roof caves followed by excavating, and the lower part of the geological main roof is cantilever state. According to the statistics of EIR with 4-7 m different high-mining working faces in five typical coal mines in the Yushenfu mining area, the results are shown in Table 1. As the mining height increases, the thickness of EIR and k i (k i is the thickness of EIR and the mining height ratio as the immediate mining ratio) increase non-linearly, as shown in Figure 1. Using the single variable method, we change the mining height on the same 5 m physical model experimental platform [16], in order to obtain the thickness of EIR at different mining heights. The results are shown in Figure 2. According to the Figures 1 and 2, we found that the thickness of EIR showed an increasing trend with the increase of the mining height, which caused the horizon of the main roof hinge structure rise; the caved roof filled the goaf area differently, which caused the different hinged structures in the overlying strata [19,20]. Many authors have proposed different roof structure models based on different mining geological conditions [29,30]. S h=11m Sh=13.5m Sh=17.5m

Mechanical Model of Key Stratum Advanced Breaking Position
The research found that when mining in shallow-buried and LMH working faces, the first and periodic breaks are all one-time full-thick breaks, the ground pressure appeared severely, the surface has obvious step sinking [1]. The structural effect of the roof was revealed, and the law of load transfer in thick loose layer of overburden is proposed [33][34][35]. The roof periodic breaking of the LMH working face in shallow coal seam can be divided into three stages: ① advanced cracks appear ahead of the working face; ② advanced cracks continue to expand; ③ the working face advances to the position of the advanced crack nearly, the roof sinks and forms a hinge structure.
The advanced breaking mechanics model of the key strata is established, as shown in Figure 3. The broken interval of the roof strata is mainly composed of the rock beam overhang length L and the advanced breaking distance x. The bending moment (M0) and shear stress (Q0) of key strata above the coal wall in working face. The mining height of coal seam (M) and the height of EIR (Σh) Using the single variable method, we change the mining height on the same 5 m physical model experimental platform [16], in order to obtain the thickness of EIR at different mining heights. The results are shown in Figure 2. According to the Figures 1 and 2, we found that the thickness of EIR showed an increasing trend with the increase of the mining height, which caused the horizon of the main roof hinge structure rise; the caved roof filled the goaf area differently, which caused the different hinged structures in the overlying strata [19,20]. Many authors have proposed different roof structure models based on different mining geological conditions [29,30].  Using the single variable method, we change the mining height on the same 5 m physical model experimental platform [16], in order to obtain the thickness of EIR at different mining heights. The results are shown in Figure 2. According to the Figures 1 and 2, we found that the thickness of EIR showed an increasing trend with the increase of the mining height, which caused the horizon of the main roof hinge structure rise; the caved roof filled the goaf area differently, which caused the different hinged structures in the overlying strata [19,20]. Many authors have proposed different roof structure models based on different mining geological conditions [29,30]. S h=11m Sh=13.5m Sh=17.5m

Mechanical Model of Key Stratum Advanced Breaking Position
The research found that when mining in shallow-buried and LMH working faces, the first and periodic breaks are all one-time full-thick breaks, the ground pressure appeared severely, the surface has obvious step sinking [1]. The structural effect of the roof was revealed, and the law of load transfer in thick loose layer of overburden is proposed [33][34][35]. The roof periodic breaking of the LMH working face in shallow coal seam can be divided into three stages: ① advanced cracks appear ahead of the working face; ② advanced cracks continue to expand; ③ the working face advances to the position of the advanced crack nearly, the roof sinks and forms a hinge structure.
The advanced breaking mechanics model of the key strata is established, as shown in Figure 3. The broken interval of the roof strata is mainly composed of the rock beam overhang length L and the advanced breaking distance x. The bending moment (M0) and shear stress (Q0) of key strata above the coal wall in working face. The mining height of coal seam (M) and the height of EIR (Σh)

Mechanical Model of Key Stratum Advanced Breaking Position
The research found that when mining in shallow-buried and LMH working faces, the first and periodic breaks are all one-time full-thick breaks, the ground pressure appeared severely, the surface has obvious step sinking [1]. The structural effect of the roof was revealed, and the law of load transfer in thick loose layer of overburden is proposed [33][34][35]. The roof periodic breaking of the LMH working face in shallow coal seam can be divided into three stages: 1 advanced cracks appear ahead of the working face; 2 advanced cracks continue to expand; 3 the working face advances to the position of the advanced crack nearly, the roof sinks and forms a hinge structure.
The advanced breaking mechanics model of the key strata is established, as shown in Figure 3. The broken interval of the roof strata is mainly composed of the rock beam overhang length L and the advanced breaking distance x. The bending moment (M 0 ) and shear stress (Q 0 ) of key strata above the Energies 2020, 13, 1685 4 of 15 coal wall in working face. The mining height of coal seam (M) and the height of EIR (Σh) as elastic foundations, the roof strata as a beam located on the elastic foundation. The support effect is offset by balancing the weight of the EIR. The broken interval is l, the overhanging section of the roof strata on the elastic foundation is reduced to R = 0.5ql by the force of the key strata, the overlying rock force is regarded as the uniformly distributed load q.

Energies
4 of 15 as elastic foundations, the roof strata as a beam located on the elastic foundation. The support effect is offset by balancing the weight of the EIR. The broken interval is l, the overhanging section of the roof strata on the elastic foundation is reduced to R=0.5ql by the force of the key strata, the overlying rock force is regarded as the uniformly distributed load q.  For the key stratum beam of unit width, the constitutive equation is: where E is the elastic modulus of the rock stratum of the main roof, MPa; I is the moment of inertia of the key stratum of roof, m 4 ; y is the vertical displacement of the beam, m; q(x) is the load of overlying strata on the key stratum beam of roof, MPa; p(x) is the force of the lower stratum on the roof key stratum, MN. According to the elastic foundation beam theory, there are: where k is the EIR and coal seam elastic foundation coefficient, MN/m 3 . According to Figure 3, the elastic foundation coefficient of the key stratum of the roof in the LMH working face is: where M is the mining height, m; Em is the elastic modulus of the coal seam, MPa; Σh is the thickness of EIR, m; E1 is the elastic modulus of EIR, MPa.
Then, the constitutive equation of the rock beam in the key stratum of the roof can be transformed into: In Equation (4), the undetermined coefficients C1, C2, C3, and C4 can be determined according to the boundary conditions. When the x→∞, e -ax →0, e ax →∞, so the C3 = C4 = 0, substituting the coefficient into the above equation can be obtained: Calculating Equation (5), the results are: For the key stratum beam of unit width, the constitutive equation is: where E is the elastic modulus of the rock stratum of the main roof, MPa; I is the moment of inertia of the key stratum of roof, m 4 ; y is the vertical displacement of the beam, m; q(x) is the load of overlying strata on the key stratum beam of roof, MPa; p(x) is the force of the lower stratum on the roof key stratum, MN. According to the elastic foundation beam theory, there are: where k is the EIR and coal seam elastic foundation coefficient, MN/m 3 . According to Figure 3, the elastic foundation coefficient of the key stratum of the roof in the LMH working face is: where M is the mining height, m; E m is the elastic modulus of the coal seam, MPa; Σh is the thickness of EIR, m; E 1 is the elastic modulus of EIR, MPa. Then, the constitutive equation of the rock beam in the key stratum of the roof can be transformed into: When the α = 4 k 4EI , the general solution of Equation (1) is: In Equation (4), the undetermined coefficients C 1 , C 2 , C 3 , and C 4 can be determined according to the boundary conditions. When the x→∞, e −ax →0, e ax →∞, so the C 3 = C 4 = 0, substituting the coefficient into the above equation can be obtained: Energies 2020, 13, 1685

of 15
Calculating Equation (5), the results are: According to the theory of elasticity, the bending moment M 0 and shear stress Q 0 of the key stratum beam at the rib in Figure 3 are: Substituting M 0 and Q 0 into Equations (7) and (8), the results are: Taking Equation (10) into Equation (6), the deflection deformation equation of the key stratum of the roof is: The maximum bending moment M max of the key beam in the main roof and the distance x ahead of the rib are: where x is the maximum bending moment position of the rock beam in the key stratum of the roof (the distance of rib advanced breaking), m. According to the maximum tensile stress criterion, the judgment conditions for the breaking of the rock stratum in the key stratum of the main roof are: where L is the overhang distance of the key strata, m; x is the leading breaking distance of the key stratum, m; q is the load of the overlying strata on the key strata, MPa; h is the thickness of the key strata, m; R t is key strata tensile strength, MPa.

Mechanical Model of EIR Advanced Breaking Position
During the advancement of the working face, it is necessary to provide a certain initial force through the support to control the movement of the immediate roof. The coal seam is used as the elastic foundation; the EIR is regarded as the elastic foundation beam located on the coal seam. The mechanical model of the EIR elastic foundation beam is established as shown in Figure 4. The EIR deflection and coal seam load under the supporting force of the support can be simplified by using the semi-infinite load beam theory.  In Figure 4, the EIR suspension elongation is l, and the length of the support roof beam is lk . The bending moment (M0) and shear stress (Q0) of EIR above coal wall in working face. The mining height of coal seam (M) and the height of EIR (Σh) as elastic foundations, the load of EIR is q1, the support strength of the support is q2.
According to the actual measurement [19][20][21], it is generally 1.5-2 times lk (lk is the length of the roof beam of the support). To ensure safety, take the l = 2lk , it can be obtained the M0 and Q0 at the coal wall (x = 0) are: where q1 is the load of EIR, MPa; q2 is the support strength of the support, MPa; lk is the length of the support roof beam, m.
Substituting Equations (16) and (17) and related parameters of EIR into Equations (7) and (8), we get: where E1 is the EIR elastic modulus, MPa; I1= h 3 /12 is the EIR inertia moment, m 4 ; km is the elastic foundation coefficient of the coal seam. Substituting Equation (18) into Equations (6) and (7), the EIR deflection y', Mx, and Qx equations can be obtained: The deformation of q1/km has already occurred when the working face excavated, so it can be ignored. Both sides of the Equation (19) multiplying km, according to the elastic foundation beam theoretical, the EIR load on the coal seam is: Equation (20) shows the maximum bending moment in the EIR occurs at Qx = 0, the EIR advances the coal wall breaking distance (x') is: In Figure 4, the EIR suspension elongation is l, and the length of the support roof beam is l k . The bending moment (M 0 ) and shear stress (Q 0 ) of EIR above coal wall in working face. The mining height of coal seam (M) and the height of EIR (Σh) as elastic foundations, the load of EIR is q 1 , the support strength of the support is q 2 .
According to the actual measurement [19][20][21], it is generally 1.5-2 times l k (l k is the length of the roof beam of the support). To ensure safety, take the l = 2l k , it can be obtained the M 0 and Q 0 at the coal wall (x = 0) are: where q 1 is the load of EIR, MPa; q 2 is the support strength of the support, MPa; l k is the length of the support roof beam, m. Substituting Equations (16) and (17) and related parameters of EIR into Equations (7) and (8), we get: where E 1 is the EIR elastic modulus, MPa; I 1 = h 3 /12 is the EIR inertia moment, m 4 ; k m is the elastic foundation coefficient of the coal seam. Substituting Equation (18) into Equations (6) and (7), the EIR deflection y , M x , and Q x equations can be obtained: The deformation of q 1 /k m has already occurred when the working face excavated, so it can be ignored. Both sides of the Equation (19) multiplying k m , according to the elastic foundation beam theoretical, the EIR load on the coal seam is: Equation (20) shows the maximum bending moment in the EIR occurs at Q x = 0, the EIR advances the coal wall breaking distance (x ) is: Maximum bending stress of equivalent direct overhanging beam: In the equation, Z = b Σh 2 /6 is the section coefficient. The width b = 1 m in this model, then Substituting Equation (26) into Equation (25), the equivalent maximum direct tensile stress is: The maximum tensile stress at the upper end of the equivalent direct overhanging beam exceeds the ultimate tensile stress.
where σ t is the uniaxial tensile strength of the EIR rock formation, MPa. Substituting M max into Equation (29), we can obtain the EIR rock strata fracture discriminant equation: The Equation (30) shows that the stability of the EIR cantilever beam structure of the large mining height working face is not only related to its size, weight, and properties, but also affected by its upper load and support force of the bracket.

Field Measurement Method
The 22201 working face is the first ever performed in the 2-2 coal seam in the Zhangjiamao Coal Mine. The average buried depth of the 2-2 coal seam is about 95 m, the mining height is 6 m, and the hydraulic support model is ZY12000/28/63D. In the middle of the 22201 working face, the cores are drilled to the roof coal and overlying strata mainly measured the compressive strength, tensile strength, elastic modulus, and Poisson's ratio of coal rock mass, experimental results, as shown in Table 2 [32]. In order to observe the whole process of roof advanced breaking position of the 22201 working face in mining. Three observing stations were fixed along the auxiliary headgate gateway (AHG) and headgate gateway (HG), three stations are located in 37 m, 67 m, and 97 m in front of the setup Energies 2020, 13, 1685 8 of 15 entry. Every observed station is including three drilling holes, which are drilled into the roof as shown in Figure 5a. By the drilling peep and the mechanical multi-point displacement meter, the breaking position of the roof is observed [32].
According to the drilling angle of every drill as shown in Figure 5b, the broken position of any point in space, the roof with horizon of broken location and the suspension distance with AHG can be calculated. In order to observe the whole process of roof advanced breaking position of the 22201 working face in mining. Three observing stations were fixed along the auxiliary headgate gateway (AHG) and headgate gateway (HG), three stations are located in 37 m, 67 m, and 97 m in front of the setup entry. Every observed station is including three drilling holes, which are drilled into the roof as shown in Figure 5a. By the drilling peep and the mechanical multi-point displacement meter, the breaking position of the roof is observed [32].
According to the drilling angle of every drill as shown in Figure 5b, the broken position of any point in space, the roof with horizon of broken location and the suspension distance with AHG can be calculated.    roof distance l k = 5.6 m, support strength q 2 = 0.9 MPa, initial support force 7916 kN, rated working resistance is 12,000 kN.

Calculation Results of the Theoretical Model
According to Equation (3), the elastic foundation coefficient of the key stratum of the roof is: Substituting I = h 3 /12 and related parameters into Equation (3), α can be calculated: Based on the field measurements [28], the average overhang distance of the key strata L = 6-8 m, the range of x can be obtained from Equation (13): x = 5.6-6.9 m Substituting parameters x = 5.6 m into Equation (15), we get 6.44 MPa > 1.14 MPa; when x = 6.9, we get 3.06 MPa > 1.14 MPa. Results indicate that the leading layer of the main roof has broken in advance during the working face, the distance in advance is 5.6-6.9 m.
(2) Calculation Result of EIR Fracture The elastic foundation coefficient of the coal seam is: Substituting I 1 = h 1 3 /12 and related parameters into Equation (18) can be obtained: The overburden strata acts on the EIR load q 1 = 0.85 MPa, the support loads the roof load of 0.9 MPa, so the x can be obtained from Equation (23).
The calculations show that the main roof and the EIR of the high-mining working face have advanced fracture phenomena in shallow coal seam. The theoretical calculation of key strata results (5.6-6.9 m) are in the range of field measurement results (5-8 m).

Field Measurement Results
(1) First Caving of the Main Roof When the working face advanced to 32 m, the main roof occurred first weighting, the first station located 5 m in front of the rib is not obviously broken. However, obvious movement (about 0.8 m displacement) was measured on the roof of the horizon of 0-10 m, which indicates the movement of roof strata was ahead of the mining working face weighting, as shown in Figure 6. Based on the actual measurement, the position of roof broken is 5 m ahead of the rib for the first weighting. The height of the caving zone is 10 m, and the hinged roof is the original bedrock on the horizon of 10-20 m. According to [1], combined with the above-mentioned measured parameters, the formation of the roof "asymmetric three-hinged arch" structure during the first caving of the main roof is shown in Figure 7.  Based on the actual measurement, the position of roof broken is 5 m ahead of the rib for the first weighting. The height of the caving zone is 10 m, and the hinged roof is the original bedrock on the horizon of 10-20 m. According to [1], combined with the above-mentioned measured parameters, the formation of the roof "asymmetric three-hinged arch" structure during the first caving of the main roof is shown in Figure 7. Based on the actual measurement, the position of roof broken is 5 m ahead of the rib for the first weighting. The height of the caving zone is 10 m, and the hinged roof is the original bedrock on the horizon of 10-20 m. According to [1], combined with the above-mentioned measured parameters, the formation of the roof "asymmetric three-hinged arch" structure during the first caving of the main roof is shown in Figure 7.  In the same way, when the main roof is first period caving, the roof is broken at 5-8 m ahead of the rib, and the roof with horizon of 16-20 m, as shown in Figure 9a. When the main roof is second period caving, a large crack occurred at the roof with horizon of 25 m and by 5-8 m ahead of the rib, as shown in Figure 9b. It shows that from the first weighting to the periodic weighting, the hinged level of the roof plate gradually moved from 10-20 m to 16-30 m. During the period weighting, the height of the caving zone is about 16 m, the roof hinge layer is at the level of 16-30 m, the periodical weighting interval is 13 m. The main roof is broken to form a periodic step rock beam structure, and the periodic step is 1.7 m.

The "Space-Time" Relationship between the Roof Structure Movement and Ground Pressure
The results of roof advanced breaking position obtained from the displacement and the drilling peep of the roof in the AHG, which can explain the "space-time" relationship between the roof structure formed by the periodic breaking and the ground pressure.
(1) Before weighting: as the working face advances, the main roof bends at the 5-8 m ahead of coal wall, and the main roof reaches the limit state, also obvious cracks are observed in the drilling results. Meanwhile, the EIR has not broken; the EIR can support the main roof. The working face is not weighting, and the main roof remains stable, as shown in Figure 10a.
(2) Starting weighting: when the working face advances to 2-4 m from the advanced breaking line, the main roof structure M block starts to rotate, and the advanced breaking line extends downward to the EIR, so the stability of the immediate roof reduced. The same time, the support capacity of the structure and the stability of the roof structure are reduced, and the pressure of the support on the measured in the 22201 working face rises and starts to come [32], as shown in Figure  10b.

The "Space-Time" Relationship between the Roof Structure Movement and Ground Pressure
The results of roof advanced breaking position obtained from the displacement and the drilling peep of the roof in the AHG, which can explain the "space-time" relationship between the roof structure formed by the periodic breaking and the ground pressure.
(1) Before weighting: as the working face advances, the main roof bends at the 5-8 m ahead of coal wall, and the main roof reaches the limit state, also obvious cracks are observed in the drilling results. Meanwhile, the EIR has not broken; the EIR can support the main roof. The working face is not weighting, and the main roof remains stable, as shown in Figure 10a.
(2) Starting weighting: when the working face advances to 2-4 m from the advanced breaking line, the main roof structure M block starts to rotate, and the advanced breaking line extends downward to the EIR, so the stability of the immediate roof reduced. The same time, the support capacity of the structure and the stability of the roof structure are reduced, and the pressure of the support on the measured in the 22201 working face rises and starts to come [32], as shown in Figure  10b.
(3) Coming time: when the working face is advanced near the advanced breaking line of the main roof, the main roof fracture intensifies and a large turning subsidence occurs (the angle of rotation is about 8°), which is EIR internal crack expansion. Therefore, the stability of the main roof structure reduces, the dynamic load increases, and the main roof structure sinks into steps, which cause the working surface to be pressed, as shown in Figure 10c.
(4) After weighting: when the working face is advanced about 4-6 m across the break line position, the supporting face of the working face enters into the new overhang EIR, the ground pressure decreases, and the weighting of working face comes to an end. At this moment, the post-frame fractured main roof rock block slipped onto the falling goaf. Due to the mining height of

The "Space-Time" Relationship between the Roof Structure Movement and Ground Pressure
The results of roof advanced breaking position obtained from the displacement and the drilling peep of the roof in the AHG, which can explain the "space-time" relationship between the roof structure formed by the periodic breaking and the ground pressure.
(1) Before weighting: as the working face advances, the main roof bends at the 5-8 m ahead of coal wall, and the main roof reaches the limit state, also obvious cracks are observed in the drilling results. Meanwhile, the EIR has not broken; the EIR can support the main roof. The working face is not weighting, and the main roof remains stable, as shown in Figure 10a.
(2) Starting weighting: when the working face advances to 2-4 m from the advanced breaking line, the main roof structure M block starts to rotate, and the advanced breaking line extends downward to the EIR, so the stability of the immediate roof reduced. The same time, the support capacity of the structure and the stability of the roof structure are reduced, and the pressure of the support on the measured in the 22201 working face rises and starts to come [32], as shown in Figure 10b.
(3) Coming time: when the working face is advanced near the advanced breaking line of the main roof, the main roof fracture intensifies and a large turning subsidence occurs (the angle of rotation is about 8 • ), which is EIR internal crack expansion. Therefore, the stability of the main roof structure reduces, the dynamic load increases, and the main roof structure sinks into steps, which cause the working surface to be pressed, as shown in Figure 10c.

The Relationship between Filling Rate of EIR and Roof Structure
Due to the different mining heights, different thicknesses of EIR will be formed, and the extent of filling of the goaf will also be different. Different advance breaking positions will cause different loads above the support; it will lead to different roof structure and ground pressure.
According to Figure 11, the height of the interspace between the EIR and the main roof is: where hc is the filling height of EIR, m; Σhi is the thickness of EIR, m; M is the mining height, m.
Σhi M hc  Figure 11. Sketch-profile evidencing variations in the relationship between kp, ki and filling rate.
The fill rate of the EIR is: (32) where δ is the EIR filling rate, %.
The filling height of the EIR collapse is:

The Relationship between Filling Rate of EIR and Roof Structure
Due to the different mining heights, different thicknesses of EIR will be formed, and the extent of filling of the goaf will also be different. Different advance breaking positions will cause different loads above the support; it will lead to different roof structure and ground pressure.
According to Figure 11, the height of the interspace between the EIR and the main roof is: where h c is the filling height of EIR, m; Σh i is the thickness of EIR, m; M is the mining height, m.

The Relationship between Filling Rate of EIR and Roof Structure
Due to the different mining heights, different thicknesses of EIR will be formed, and the extent of filling of the goaf will also be different. Different advance breaking positions will cause different loads above the support; it will lead to different roof structure and ground pressure.
According to Figure 11, the height of the interspace between the EIR and the main roof is: where hc is the filling height of EIR, m; Σhi is the thickness of EIR, m; M is the mining height, m. The fill rate of the EIR is: where δ is the EIR filling rate, %.
The filling height of the EIR collapse is: Figure 11. Sketch-profile evidencing variations in the relationship between k p , k i and filling rate.
The fill rate of the EIR is: where δ is the EIR filling rate, %. The filling height of the EIR collapse is: Then, the thickness of EIR and mining height ratio is Substituting Equations (33) and (34) into Equation (32), we can get: According to Equation (35), if the k p is constant, when the k i increases, the fill rate also increases. The changing relationship between k p , k i and filling rate shown in Figure 12.
Substituting Equations (33) and (34) into Equation (32), we can get: According to Equation (35), if the kp is constant, when the ki increases, the fill rate also increases. The changing relationship between kp, ki and filling rate shown in Figure 12.

Conclusions
(1) Based on the theory of the "elastic foundation beam", mechanical models for the key strata and EIR advanced breaking position of large mining height working face in shallow coal seam are established. According to the field measured physical and mechanical parameters of the 22201 LMH working face, the theoretical calculation of key strata results (5.6-6.9 m) are in the range of field measurement results (5-8 m).
(2) In mining practice, the roof with horizon of 0-10 m has detected advanced movement displacement in the 22201 LMH working face. It is shown that the advanced breaking of the LMH working face in shallow coal seam can be used as an indication of roof weighting, providing early warning for the site of large-mining work face, and ensuring the safety of the mining.
(3) Under the condition of different filling rates, different EIR have different filling degrees of goaf. The degree of filling in the goaf determines the shape and stability of the main roof structure of the stope in shallow coal seam. Different roof structures will lead to different degrees of weighting.