Key Stratum Structure and Support Working Resistance of Longwall Face with Large Mining Height in the Shallow Coal Seams, China

The fully mechanized mining with large mining height is the main method for high yield and eﬃcient coal mining in China. The key stratum structure (KSS) is the basis of revealing the mechanism of roof weighting and determination of support working resistance of the longwall face with large mining height (LFLMH) in the shallow coal seam. The height of the caving zone at LFLMH is large, the thick immediate roof forms the “short cantilever beam” structure commonly, and the hinge layer of the overlying key stratum will move upward to the higher position. The “high position oblique step voussoir beam” structure of single-key stratum (SKS) and “oblique step voussoir beam and voussoir beam” structure of double-key stratum (DKS) in the shallow coal seam were proposed with physical simulation and Universal Distinct Element Code (UDEC). The analysis of the KSS and numerical simulation reveals the mechanism of strong roof weighting at the SKS longwall face and large-small alternate periodic weighting at the DKS longwall. It is concluded that the large static load caused by the “equivalent immediate roof (EIR)” is the basic load, and the instability load of the KSS is the additional dynamic load of support. Besides, the calculation methods of the reasonable support working resistance at LFLMH were obtained and veriﬁed with engineering applications.


Introduction
e reserves of thick coal seams contribute about 44% of the country's total coal reserves, and the coal thereof produced accounts for approximately 45% of the total raw coal output in China [1,2], indicating that the safe and efficient mining of thick coal seams is critical to Chinese energy supply. Shenfu-Dongsheng Coalfield, as depicted in Figure 1, is one of the largest coal production coalfields worldwide, and there are a large number of coal seams characterized by the shallow cover and extreme thickness in the northwest of China [3]. Fully mechanized mining with large mining height technology with the advantages of high coal recovery rate and low roadway driving rate has been the common method for high-efficiency coal mining with the thickness of coal seam ranging from 3.5 to 8 m in China [1,[4][5][6]. China is also the country with the most coal mine accidents in the world [7,8], including the coal wall falls, roof falls, coal burst, and support failure accidents, affecting coal production and its efficiency seriously [9][10][11][12]. In 2017, 375 people died of coal mine disasters, where 32% are roof accidents, indicating that it is quite essential to ensure the roof stability for the safety of mining workers and equipment [8]. It is an essential issue to find out the mechanism of roof instability and determine the effective support working resistance for the longwall face in the shallow thick seams [13].
With the development of the Shendong Coalfield, Chinese scholars have conducted systematic studies on the KSS and the theory of ground control in the shallow coal seam. e ''key stratum theory'' [14] provides a new approach to study the KSS and mine pressure behaviors, where the stratum that controls the movement of part or all of the overlying strata was defined as the key stratum (KS). Huang has proposed the "asymmetric three-hinged arches" structure of the roof in the first weighting [15], the "short voussoir beam" and "step voussoir beam" structures of the roof in the periodic weighting [16], and the load transfer theory of unconsolidated layers [17]. All these provide the theoretical basis for the roof control in shallow coal seam mining.
In recent years, with the rapid development of large mining height mining equipment and technology, the maximum mining height has reached 8.8 m in Shendong Coalfield, China, the mine pressure behaviors of the longwall face are strong, and the roof control is complicated. e relationship between support and surrounding rock at LFLMH has become a study hotspot. Ju and Xu have analyzed the structural characteristics of key stratum and strata behavior of the longwall face with super large mining height [18] and proposed the "cantilever beam and voussoir beam" structure, which leads to the alternating roof pressure [19]. Yan et al. [20] have proposed the "short cantilever beam and articulated voussoir beam" structure at LFLMH and found that the fracture interval of the roof has a significant influence on the support working resistance. Based on the statistical analysis for field measurements of support load at the longwall faces with 3.5 to 5.5 m mining height, Gong [21] has revealed that the influence of the main roof structure on the support load decreases, and the influence of the immediate roof on the support load increases. Huang et al. have proposed the concept of EIR [22] and the "oblique step voussoir beam" structure at LFLMH [23] and studied the crack development characteristics under different mining conditions in the shallow coal seam [24]. e abovementioned studies have promoted the development of roof control at LFLMH. However, the effect of filling status of the goaf on SKS and DKS structures in the shallow coal seam still needs to be studied because the immediate roof has a significant impact on the stability of the KSS and roof control at LFLMH. A clear scientific understanding of the KSS and roof weighting mechanism of LFLMH in the shallow coal seam is still lacking.
is study constructed the mechanical models of SKS and DKS roof structures at LFLMH and revealed the mechanism of roof weighting of the longwall face and proposed the calculation methods of the reasonable support working resistance. We wish that the study's results are used as a reference for roof control at LFLMH in the shallow coal seam.

Caving Characteristics of the Immediate
Roof at LFLMH e height of the caving zone at LFLMH increases linearly with the increase of mining height [13,[25][26][27]; the hinge layer of overlying strata moves upward in the higher position. In this part, the caving characteristics of the immediate roof were studied using the physical simulation experiment.

Geological Conditions.
e physical simulation experiment was based on the mining conditions of No. 5-2 seam in Zhangjiamao coal mine, located in Yulin, Shaanxi, China.
e No. 5-2 seam belongs to a typical shallow coal seam with average depth and thickness of 120 m and 6.1 m. e roof of No. 5-2 seam is mainly made up of the siltstone, the fine grained sandstone, and the overlying strata with thin bedrock and thick unconsolidated layers. e composite histogram of No. 5-2 seam is shown in Figure 2.

Physical Simulation Analysis.
Physical simulation is one of the most effective techniques in research on rock strata behavior [28][29][30], which has been applied to simulate longwall mining-induced subsidence and other related problems [31][32][33]. Due to the difference in scale between coal fields and laboratory experiments, proper model materials must be selected to maintain meaningful physical proportions to the field conditions. e properties of model material need to be tested to satisfy the similarity theory principles [32,34]. According to similarity theory and field conditions, the physical simulation model was built with the dimensions of 5 m × 0.2 m × 1.2 m, the geometric ratio of the simulation model was determined as 1 : 50, the bulk density ratio is 1 : 1.5, and the strength ratio is 1 : 75; the overlying strata that the model has not simulated were replaced by the additional loads on the top of the model. e caving characteristics of the immediate roof were simulated with the mining height of 4 m, 5 m, 6 m, and 7 m on the same model, as shown in Figure 3. e material ratio of the physical simulation model is listed in Table 1 based on the rock mechanics property tests.
Notes: A, B, C, and D are the longwall faces simulation areas with the mining height of 4 m, 5 m, 6 m, and 7 m, respectively. e physical simulation results show that the thickness of the immediate roof at LFLMH is large. e immediate roof will not cave immediately after coal mining and form the "short cantilever beam" structure ( Figure 4(a)). With the continued advance of the longwall face, the immediate roof will cave (Figure 4(b)). Some lower main roofs cannot form the hinged structure, but they can cave and fill the goaf to act as the "immediate roof." e strata below the KS that cannot form the hinged structure is defined as the EIR [22].
According to the physical simulation results, when the mining heights are 4 m, 5 m, 6 m, and 7 m, the thicknesses of EIR is 10 m, 15 m, 21 m, and 26 m, respectively, as shown in Figure 5. e thickness of the EIR linearly correlates with mining height.

Classification of the EIR.
When the thick EIR fractures, a large static load will be generated on the support, and the load is needed to be controlled by support setting load. e different filling status of the goaf will cause different stability of KSS and the strength of the roof weighting at the longwall face.
Based on the physical simulation results, the field observations of step subsidence [3], and the KS theory [14], the EIR can be divided into "fully filling type" and "general filling type." e "step voussoir beam" structure can be formed in general filling type, and it is the universal type of overlying strata during the shallow coal seam mining.
(1) Fully Filling Type. e thickness of the EIR is over 3 to 3.3 times the mining height. Under these conditions, the rotation space of the main roof structure is small; thus, the main roof can form the "voussoir beam" structure. (2) General Filling Type. When the EIR cannot fully fill the goaf, the main roof structure will present step subsidence and form the "step voussoir beam" structure. Such thickness of EIR is generally less than 3 times the mining height, which is a universal condition at LFLMH. e "step voussoir beam" structure is prone to sliding instability [16]; the roof pressure of the "step voussoir beam" structure is larger than the roof pressure of the "voussoir beam" structure. erefore, the determination of support working resistance should be based on this KSS.

Numerical Model Establishment.
e No. 15201 LFLMH (LW15201) at Zhangjiamao coal mine and the No. 22303 LFLMH (LW22303) at Bulianta coal mine are taken as the background to study the overlying strata structure and mechanism of roof weighting at LFLMH in the shallow coal seam. e UDEC numerical models were established with the Mohr-Coulomb failure criterion to study the overlying strata structure during the advance of the longwall face. e lateral boundaries are fixed in displacement at the horizontal direction, the bottom boundary is fixed in displacement at the horizontal and vertical directions, and the top boundary is free. e dimensions of the two models were 500 m in length × 120 m in height and 500 m in length × 172 m height, respectively. e 100 m boundary coal pillars were set at the left and right ends of the models. e length of the simulation mining region is 300 m in the middle of the model. e reasonable parameter setting and meshing are crucial to ensure the accuracy of the simulation results. Mohammad et al. [35] have suggested that the stiffness of the numerical model should be equal to 0.47 times the average stiffness value of the laboratory tests. Cai et al. [36] have suggested that the elastic modulus of coal and rocks should be 0.1 to 0.25 times that of the laboratory tests, while Poisson's ratio of coal and rocks is assumed to be 1.2 to 1.4 times that of the laboratory tests. erefore, the elastic modulus and tensile strength of coal and rocks are equal to    Tables 2 and 3. e grids in this numerical model were divided reasonably, referring to the UDEC numerical simulation of overlying strata movement at longwall faces [37,38].

Simulation Results Analysis.
e simulation results of LW15201 at Zhangjiamao coal mine show that the overlying strata of the longwall face can form the SKS roof structure with the average breaking angle of the KS being 70°. e SKS roof structure is prone to sliding, making the overlying load fully transfer and the roof pressure of the longwall face increase rapidly. Figure 6 shows the "high position oblique step voussoir beam" SKS roof structure.   e simulation results of LW22303 at Bulianta coal mine show that the overlying strata of the longwall face can form the DKS roof structure, the lower KS forms the "lower oblique step voussoir beam" structure, and the upper KS forms the "upper voussoir beam" structure.
ere are two main types in the DKS roof structure: (1) when the DKS fractures synchronously (Figure 7(a)), the support load of the longwall face caused by the overlying load fully transfer is large, forming large periodic weighting of the longwall face; (2) when the upper KS fractures later than the lower KS (Figure 7(b)), the strata below the upper KS can form bedding separation, the support load of the longwall face is mainly caused by the weight of the lower KS, forming small periodic weighting of the longwall face. During the coal mining process, the large and small roof weighting of the longwall face alternate.

"High Position Oblique
Step Voussoir Beam" Structure. For the typical shallow coal seam [39], the thickness of the bedrock is small, and there is only SKS structure in the main roof. According to the field observations of step subsidence [3], the KS theory [14], the physical simulation experiment ( Figure 8) [22,23], and the numerical simulation results, the filling status of the goaf at LFLMH is general filling and the thickness of EIR is large. us, the hinge point position of the main roof will move higher, and the overlying strata will form the "high position oblique step voussoir beam" structure, as shown in Figures 8 and 9.
In Figure 9, h 1 is the thickness of the KS, m; h i is the thickness of the EIR, m; m is the mining height, m; h 2 is the thickness of the "short cantilever beam," m; h 3 is the The ratio of the thickness of the EIR to mining height y = -11.7 + 5.4x   structure; R 1 is the force on the EIR from the block M, kN; R 2 is the weight of the "short cantilever beam," kN; R 3 is the weight of easily caving EIR, kN; P m is the support load, kN; A, B, and C are the hinge point of the key blocks. P m can be calculated using the following formula: where R 1 , R 2 , and R 3 are calculated with support width b.

Stability Analysis of "High Position Oblique
Step Voussoir Beam" Structure. Based on the "step voussoir beam" structure model [16], the mechanical model of "oblique step" key blocks can be described as in Figure 10. e key block N falls entirely on the caving rocks, and the key block M is supported by the key block N at point C. At this point, the key block N is compacted; it is desirable that R N ≈ P 2 , Q B ≈ 0. Taking the torque for the key block M at point C, ΣM C � 0: where T is the horizontal squeezing force, kN; β is the breaking angle of the KS,°; θ is the rotation angle of the key block M,°; W is the subsidence height of the key block N, m; a is the height of contact surface, m; P 1 is the sum of the weight of the key block M and the overlying load, kN; L is the length of the key blocks M and N, m; d is the step height of the key blocks M and N, m; Q A is the shear force at point A, kN.
Taking the entire vertical force ΣF y � 0, the equilibrium equation is as follows: where Q B is the shear force at point B, kN; R N is the support force for the key block N, kN; and P 2 is the sum of the weight of the key block N and the overlying load, kN. Substituting R N ≈ P 2 and Q B ≈ 0 into equations (2) and (3), T can be derived, as follows: According to the "S-R" stability analysis [40], the "high position oblique step voussoir beam" structure is prone to forming sliding instability, and the condition for preventing the KSS sliding instability is where tan φ is the friction coefficient between the blocks and R 1 is the support force for maintaining the stability of the block M, kN. rough equations (4) to (6), the supporting force for maintaining the stability of "high position oblique step voussoir beam" structure is where P 1 consists of the weight of the key block M (R G ) and the weight of the load layer (R Z ) [41]: where b is the support width, m; c is the bulk density of the bedrock, kN/m 3 ; K G (≤1) is the load transfer coefficient [41]; h Z is the thickness of the load layer, m; c Z is the bulk density of the load layer, kN/m 3 . rough equations (7) to (10), R 1 can be calculated using the following formula:

Mining direction
High position oblique step voussoir beam

Determination of Support Working Resistance at SKS Longwall Face.
e composition of the support load of the longwall face is as shown in Figure 11. e weight of the thick EIR (R 2 + R 3 ) is the basic static load, and the instability load of SKS structure (R 1 ) is an additional dynamic load of the support.
According to the geometric conditions in Figure 11, R 2 and R 3 can be calculated using the following formula, respectively: where l is the length of "short cantilever beam," m; l k is the roof control distance of support, m; and α is the breaking angle of the EIR,°.
Based on the mine pressure theory [40], the subsidence height of the key block N can be calculated using the following formula: where K p is the bulking coefficient of the rock mass in the goaf, set to 1.3. According to the experimental results [42], take the value of tanφ as 0.5; because the values of θ and a are small, the values of θ and 0.5a can be set to 0. Substituting the above parameters into equations (1) and (11) to (14), meanwhile considering the actual efficiency of support, the reasonable support working resistance at SKS longwall face can be calculated using the following formula: where μ is the real efficiency of support, which is set to 0.9.

Mining direction
High position oblique step voussoir beam Advances in Civil Engineering 7

DKS Structure.
When the thickness of the bedrock is large, the multi-key strata may exist in the main roof. According to the KS theory, the KS is close to the longwall face, which has a significant impact on the support load. If the spacing of the DKS is small, two KSs interact and generate the DKS structure effect. e spacing that can cause the DKS structure effect is related to the breaking angle of the KS and the periodic weighting interval, as shown in Figure 12.
When the front hinge joint of the upper KS acts on the lower-key block M 1 , the load of the upper KS can be transferred to the longwall face, according to the geometric conditions in Figure 12; h j·max can be calculated using the following formula: where h j·max is the maximum spacing of the DKS, m, and L 1 is the length of the lower-key block (approximated to be periodic weighting interval), m. erefore, the strength of the periodic weighting of the longwall face is affected if the spacing h j between lower KS and upper KS satisfies the following formula: Since the periodic weighting interval is generally 10 to 20 m, and the breaking angle of the lower KS is about 70°at LFLMH of the shallow coal seam in Shendong mining area [43], according to equation (17), the overlying strata can form the DKS structure effect when h j ≤ 27 m; the upper key block M 2 can influence the stability of the lower KS. Figure 13 presents the physical simulation experiment results of the longwall mining with large mining height. When the caving EIR does not fully fill the goaf, the key blocks M 1 and N 1 will present "step subsidence," and the lower KS will form the "oblique step voussoir beam"  Advances in Civil Engineering structure. At this point, the upper KS is supported by the caving EIR, the lower KS, and the weak strata. When the caving EIR reaches fully filling status, the rotation angle of the key block M 2 is small, and the upper KS will form the "voussoir beam" structure. erefore, the overlying strata of the longwall face can form the DKS structure of the "oblique step voussoir beam and voussoir beam."

Large Periodic Weighting.
When the DKS fractures synchronously, the overlying strata of the longwall face will move fully. e instability load of the upper KS has a dynamic impact on the lower-key block M 1 , and the support load of the longwall face is large, forming a large period weighting. According to the field measurements and physical simulation [43], the large period weighting interval is two times the small period weighting interval, and the roof control of the longwall face is the most challenging work during the large period weighting.

Small Periodic Weighting.
When the lower KS is unstable, and the upper KS is stable, the strata below the upper KS will form bedding separations. e upper KS will obstruct the load transfer of overlying strata, and the support load is mainly made up of the weight of the lower-key block M 1 and the weak strata. At this point, the support load of the longwall face is small and forms a small periodic weighting of the longwall face.

Analysis of the KSS during the Large Periodic Weighting.
When the DKS fractures synchronously, the load of overlying strata will fully act on the longwall face, the support load reaches the largest, and the longwall face can form the large periodic weighting; determining the reasonable support working resistance during the coal mining process should be based on this KSS; the overlying strata structure is shown in Figure 14.
In Figure 14, h 1 is the thickness of the lower KS, m; h 4 is the thickness of the weak strata between the upper KS and lower KS and h 5 is the thickness of the upper KS, m; R 1 is the force on the EIR from the block M 1 , kN; R 4 is the sum of the weight of the block M 1 and weak strata, kN; R 5 is the force on the weak strata from the block M 2 , kN; P 4 is the sum of the weight of the key block M 2 and the overlying load, kN; L 1 and L 2 are the lengths of the key blocks M 1 and M 2 , respectively, m; W 1 and W 2 are the subsidence heights of the key blocks N 1 and N 2 , respectively, m; and θ is the rotation angle of the key block M 1 ,°.
As shown in Figure 14, the support load is mainly made up of the basic load provided by the static load of the EIR and the additional load offered by the instability dynamic load of the "lower oblique step voussoir beam" structure. e effect of upper KS is reflected by its instability load acting on lower KS and leads to an increase in the support load. e composition of the support load also by R 1 , R 2 , and R 3 can be calculated using equation (1). e load acting on the EIR from the key block M 1 can be calculated using the following formula: where P 3 is the sum of the weight of the key block M 1 and the overlying load, kN.
Since the DKS fractures synchronously, the instability load of the upper KS and the weight of weak strata will simultaneously act on the lower KS. erefore, P 3 includes the sum of the weight of the key block M 1 and the weak strata (R 4 ), as well as the load acting on the weak strata from the key block M 2 (R 5 ). R 4 can be calculated using the following formula:

Advances in Civil Engineering
According to Huang [43], R 5 can be calculated using the following formula: where φ is the friction between the blocks (°), and P 4 can be determined by the calculation method for P 1 . P 3 can be calculated using equations (19) and (20): rough equations (1), (12), (13), (18), and (21), the support working resistance during the large periodic weighting can be calculated using the following formula:

Analysis of the KSS during the Small Periodic Weighting.
When the lower KS is unstable, and the upper KS is stable, upper and lower KS form delamination. Figure 15 shows the overlying strata structure of the longwall face. e load transfer of the overlying strata is obstructed, and the longwall face can form the small periodic weighting. At this point, the support load is mainly made up of the instability load of the lower KS and the weight of the EIR.
During the small periodic weighting of the longwall face, the load of the lower-key block M 1 is the only weight of the weak strata, and P 3 can be calculated using the following formula: Combining equations (1), (12), (13), (18), and (23), the support working resistance during the small periodic weighting can be calculated using the following formula:  overlying strata load during the large periodic weighting, and the values of W 1 and W 2 can be calculated by equation (14). Considering the real efficiency of support, the reasonable support working resistance at DKS longwall face can be calculated using the following formula:

Calculation of Support Working Resistance at SKS
Longwall Face. According to the borehole column histogram and KS identification of LW22406 at Halagou coal mine, an SKS structure is found in overlying strata [44]. e correctness of the calculation formula of support working resistance at the SKS longwall face can be verified with the parameters of LW22406. Based on the geological data, mining conditions, and the simulation experiment results, the following parameters are determined: the mining height m is 5.2 m, the thickness h 1 of the KS is 10 m, the bulk density c of the bedrock is 25 kN/m 3 , the bulk density c Z of the load layer is 25 kN/m 3 , the length L of the key block (approximated to be periodic weighting interval) is 12.1 m, the thickness h z of the load layer is 10 m, the thickness h i of the immediate roof is 15 m, the length l of the "short cantilever beam" is 10.5 m, and the values of h 2 and h 3 are 13 m and 2 m, respectively. e LW22406 was equipped with ZY12000/25.5/55D-type supports, and the support width b is 1.75 m, the roof control distance l k of support is 5.0 m, the breaking angle β of the KS is 65°on average, the breaking angle α of the lower EIR is 60°on average, according to engineering experience, the value of K G is 0.38.
Substituting the above parameters into equation (15), the reasonable support working resistance of the longwall face can be calculated as 10620 kN. According to the monitoring data of support load at LW22406 (Figure 16), the average support working resistance during the periodic weighting is 84% to 92% of rated support working resistance (12000 kN), which is very close to the support working resistance in the field practice, indicating that the calculation formula can provide some guidance for support selection.

Calculation of Support Working Resistance at DKS Longwall Face.
rough the KS identification with three borehole column histograms of LW22303 at Bulianta coal mine in Shendong mining area, the DKS structure is found in the overlying strata [45]. e correctness of the calculation formula of support working resistance at the DKS longwall face can be verified with the parameters of LW22303. Based on the geological data, mining conditions, and the simulation experiment results, the following parameters are determined: the mining height m is 6.8 m, the bulk density c of the bedrock is 25 kN/m 3 , the length L 1 of the lower-key block is 13.2 m, the length L 2 of the upper key block is 25 m, the thickness h i of the immediate roof is 17.3 m, the length l of the "short cantilever beam" is 11.8 m, the values of h 1 to h 5 are 10.2 m, 14 m, 3.3 m, 4.6 m, and 8.7 m, respectively, and the value of P 4 is 1725 kN. e LW22303 was equipped with ZY16800/32/70D-type supports, and the support width b is 2.05 m, the roof control distance l k of support is 6.62 m, the breaking angle β of the KS is 65°on average, the friction angle φ between the blocks is about 27°, and the breaking angle α of the lower EIR is 60°on average.

Advances in Civil Engineering
Substituting the above parameters into equation (25), the reasonable support working resistance of the longwall face can be calculated as 15860 kN. According to the field measurements, the average support working resistance during the coal mining process is shown in Figure 17. e support working resistance shows the large-small alternate phenomenon, which verifies the correctness of the KSS analysis, and the theoretical calculation result is very close to the support working resistance during the large periodic weighting. e rated support working resistance of supports in the field is 16800 kN, which meets the requirements for roof control.

Conclusions
Based on the physical simulation experiment and the numerical simulation, the caving characteristics and the cantilever fracture pattern of the EIR at LFLMH were analyzed, and the SKS and DKS roof structures models in the shallow coal seam were constructed.
rough the mechanical analysis, the mechanism of roof weighting at LFLMH was revealed, and the calculation formulas of the reasonable support working resistance of the longwall face were proposed. e conclusions are as follows: (1) e roof caving space increases at LFLMH in the shallow coal seam. Some lower main roof is challenging to form the hinged structure, and it can fill the goaf. e thickness of the immediate roof increases largely and presents a "short cantilever beam" structure commonly. e main roof can form a high position hinged structure. e studies show the ratios of the immediate roof static load to the increase in roof pressure. e value of the dynamic load is small, and the overall roof pressure of the longwall face is large.
(2) According to the filling status of the EIR to goaf, the EIR can be divided into "fully filling type" and "general filling type." e fully filling of goaf can form the stable "voussoir beam" structure with relatively small roof weighting of the longwall face. e general filling of goaf can form the "high position oblique step voussoir beam" structure with strong roof weighting of the longwall face, and this KSS is the universal form at LFLMH in the shallow coal seam. (3) When the overlying strata only have SKS of LFLMH in the shallow coal seam, the main roof can form the "high position oblique step voussoir beam" structure, and the sliding of the KSS is the root cause to roof weighting of the longwall face that increases rapidly. When there is DKS in overlying strata, the "oblique step voussoir beam, and voussoir beam" structure can be formed. e instability of the DKS synchronously causes large periodic weighting, and the instability of the lower KS causes small period weighting. During the coal mining process, the large and small roof weighting of the longwall face alternate. Besides, the calculation formulas of the reasonable support working resistance were proposed, which can provide a reference for the support selection and roof control at LFLMH in the shallow coal seam.

Data Availability
e data used to support the findings of this study are included within the paper. Other data used to support the findings of this study are available from the corresponding author upon request.

Conflicts of Interest
e authors declare that there are no conflicts of interest regarding the publication of this paper.   Figure 17: e average support working resistance during the coal mining process at LW22303.