Abstract

Roof-breaking characteristics and ground pressure behavior of the coalface are instrumental in guiding deep Jurassic coal seam mining, in particular in the Shaanxi and Inner Mongolia regions of China. A thick-plate mechanical model (TPMM) of the main roof was developed and applied to the case study of 21102 first-mined coalface (FMC) of the Hulusu Coal Mine (HCM) in the Hujirt Mining Area (HMA), China. A theoretical analysis performed via the developed model revealed that the first and periodic breaking intervals of the main roof were 40.6 and 25.0 m, respectively. The roof failure occurred in the tensile mode, was controlled by the internal stress in the rock strata, and started from the center of the long side with the fixed support in the goaf. The field measurement of roof weighting was also performed for the coalface advance from zero to 400 m. The measurement results showed that the first weighting average interval was 41.4 m, and the average interval of periodic weighting was 22.0 m, which agreed with the theoretical calculation and proved the proposed model’s feasibility. Finally, the frequency distribution features of the hydraulic support working resistance in the FMC were analyzed statistically. The results showed that the ZY10000-16/32D supports could adapt to the mining geological conditions of the FMC. However, the margin of the rated working resistance of supports was still small. Thus, roof management enhancement during the mining process was strongly recommended. These research findings could offer theoretical guidance for safe and high-efficiency production in the coal mines under similar geological conditions.

1. Introduction

In recent years, the development of coal resources in China has rapidly shifted to Northwestern China, especially the Inner Mongolia Autonomous Region [1]. From January to December 2020, China’s northwestern coal mining areas accounted for about 70% of the total national coal output and became a forepost of China’s coal industry upgrading [25]. A good example is the Ordos City in the Inner Mongolia Autonomous Region, in which coal production in 2018, 2019, and 2020 amounted to 616, 679, and 640 million tons, respectively.

The Dongsheng Coalfield in the Ordos City extends for about 100 km in the north-south direction and has the largest width in the east-west direction of 100 km, covering an area of about 8790 km2 [6, 7]. The main coal-bearing seams in the Dongsheng Coalfield belong to the Middle Jurassic Yanan Formation, which lithology consists primarily of medium sandstone, siltstone, mudstone, and coal seam [8]. The coal seam’s burial depth generally drops from north to south but increases from east to west. The coal seams mined at an early stage are generally typical shallow-buried (below 300 m) with simple geological conditions, which have been investigated theoretically and subjected to field studies. Abundant reports have been published concerning roof-breaking characteristics, mining-induced fissure development, and ground pressure behavior [918]. These findings lay a theoretical foundation and offer ready technologies for the safe and high-efficiency mining of shallow-buried coal seams at an early stage. However, due to the variabilities of the burial depth of the coal seam and strata structure, the first-mined coal seams in some newly developed large-scale coal bases in the Dongsheng Coalfield have already entered the stage of deep mining [19]. One example of this is the HMA, a state planning mining area located in the midwestern part of the Dongsheng Coalfield. The main coal seams in the HMA are #2 and #3 coal seams in the upper Jurassic Yanan Formation, with a depth above 600 m. The available roof-breaking characteristics and ground pressure behavior in mining shallow-buried coal seams can no longer offer guidance for deep coal resource development in this region. It is necessary to carry out in-depth investigations into the occurrence conditions of the Jurassic deep-buried coal seam in this mining area.

The HMA is located on the border of Inner Mongolia and Shaanxi. The total planning area is about 2161 km2. According to preliminary planning, this mining area is subdivided into seven mine fields, two exploration areas, and one prospective area. The total planning area is 63 Mt/a. HCM, located in the mid-southern part of the mining area, is one of the first developed in this region and has a design production capacity of 13.0 Mt/a. The burial depth of the mined coal seam is 650-900 m. The HCM is a typical newly built deep coal mine in Inner Mongolia and Shaanxi. The increased mining scope was accompanied by a series of problems: large intensity of overburdened rock activities in the stope, intense ground pressure at the coalface, and large deformation of the roadway-surrounding rocks [20]. However, there are only a few studies on roof breaking and ground pressure in this mining area [2125].

Given this, we elaborated a TPMM for the main roof breakage in the 21102 FMC of the HCM in the HMA under the actual geological conditions. A theoretical analysis of the main roof’s first and periodic breaking was performed using the established model. On this basis, we conducted a field measurement to study ground pressure behavior at the FMC. The applicability of the hydraulic supports used for this FMC was evaluated. Our research findings offer theoretical guidance for the normal mining of other coalfaces in the HCM. Moreover, these results can be utilized to guide the safe and high-efficiency coal production under similar geological conditions in the HMA.

2. Mining Technical Conditions for the FMC

The 21102 coalface of the HCM, the FMC in the entire mine field, is located in the first panel of the #2-1 coal seam. There are three roadways on the east wing of the mine north of the coalface. In the south, there is the fifth panel of the coal seam. The length of the coalface along the dip direction is 320 m, and that along the strike direction is 4150 m. The stoping area is 1,327,900 m2. The ground surface above the coalface is hilly sand land, with an elevation of +1304 to +1328 m. The elevation of the underground coal seam roof is +673 to +694 m. The burial depth of the coalface is 626.17 to 647.91 m. The coalface under study is a typical deep-lying coalface in this region. The coal seam mined in the coalface is the #2-1 coal seam of the Jurassic Yanan Formation, with a density of about 1.31 t/m3. The coal seam thickness is 1.7-3.0 m, with an average of 2.55 m, indicating a medium thickness. The coal seam is of a simple structure with slight undulation and a dip angle of 1°-3°. This coal seam is near-horizontal.

The geological structure of this coalface is also simple. The false roof of the coal seam is mainly composed of sandy mudstone, with a thickness of 0.2-0.5 m. The immediate roof comprises sandy mudstone and silty sandstone, with a thickness of 2.43-7.8 m. The main roof is composed of fine sandstone and medium sandstone, with a thickness of 13.30-23.3 m. The immediate floor comprises sandy mudstone and silty sandstone, with a thickness of 3.4-8.9 m. The stratigraphic synthesis column of the study area is shown in Figure 1. In this coalface, the full-thickness longwall mining method along the strike (with an average mining height of 2.55 m) is performed. The goaf roof is managed by using the full caving method. The hydraulic supports used in the coalface are ZY10000-16/32D two-column shield-type hydraulic supports manufactured by China Coal Beijing Coal Mining Machinery Co., Ltd. The rated working resistance of the supports is 10,000 kN, and the distance between the centers of the two supports is 1.75 m.

3. Mechanical Analysis of Breaking Characteristics of the Thick Sandstone Main Roof in the FMC

A goaf is formed along with the constant advance of the coalface from the open-off cut. The first breaking and weighting of the main roof occur when the roof overhang is too big. As the coalface continues to advance after the first breaking, the periodic weighting of the main roof occurs. Since the main roof of the FMC is composed of hard thick sandstone, the conventional beam model may not apply to the theoretical analysis of the roof-breaking characteristics [26]. According to the field measurement of ground pressure in the deep-buried mine near the HCM, we considered the influence of internal transverse shear stress of the plate based on the judgment criteria of the plate model. We developed the Reissner thick-plate model for the breaking of the thick sandstone main roof and performed a theoretical analysis of stress distribution in the main roof at different mining stages in the coalface. The inducing factors of the first and periodic breaking of the main roof were determined. The first and periodic weighting intervals were calculated for this coalface.

3.1. Establishment of the TPMM for Main Roof Breaking

Since the 21102 coalface is the FMC, the main roof can be simplified into a thick rectangular plate that has fixed support on the four sides upon the moment of the first breaking. Under this assumption, we built the TPMM for the first breaking of the main roof, as shown in Figure 2(a). After the first breaking of the main roof, broken rock masses in the goaf form a three-hinged arch. At this moment, the main roof can be reduced to a thick rectangular plate with simple support on one side and fixed support on the other three sides. Thus, the TPMM for periodic breaking of the main roof was built, as shown in Figure 2(b), where and are the short and long sides of the thick rectangular plate, respectively; is its thickness; and is the load acting on the upper part of the thick rectangular plate.

According to the Reissner plate theory, the basic equations can be written as follows [2729]: where is the bending stiffness of the main roof, (N·m), is the elastic modulus of the main roof (GPa), is the deflection of the main roof (m), is the Laplace operator, is the loading acting on the upper part of the main roof (MPa), is Poisson’s ratio of the main roof, is the main roof thickness (m), and is the stress function within the main roof.

The expressions of the shear stress and the bending moment for the main roof can be derived as follows:

3.2. Mechanical Analysis of the First Breaking Characteristics of the Main Roof
3.2.1. A Theoretical Analysis of the Mechanical Model Describing the First Breaking of the Main Roof

The governing equations for the fixed support of the four sides before the first breaking of the main roof in the FMC have the following forms:

Using the fixed support boundary conditions for the four sides of the main roof, the following deflection equation that satisfied the boundary conditions was derived:

Substituting Equation (8) and (the uniformly distributed load applied to the upper part of the main roof) into Equation (1), we can calculate the coefficient by the Matlab software program:

Substituting Equation (9) into Equation (8), we obtain the deflection function under the fixed support boundary conditions for the four sides of the main roof:

Substituting Equation (10) into Equations (3)–(6), the shear stress function and the bending moment function of the main roof are derived as follows: where .

According to the theory of elasticity and thick-plate theory, the stress in the main roof is related to the bending moment as . From this, we can derive the stress function of the main roof:

3.2.2. First Breaking Interval of the Main Roof and Its Mechanical Characteristics

According to the stress function calculated in Section 3.2.1, at and , both reach their peak negative values (that is, are tensile stresses, ). Both are located at the centers of the long sides in the TPMM. are given, respectively, as follows:

If in the main roof reaches the tensile strength , the long side with fixed support in the TPMM undergoes the first tensile failure. According to the geological conditions for the FMC, the basic parameters are substituted into the above formulas to estimate the first breaking interval of the main roof as  m.

Below, we analyze the mechanical characteristics of the main roof upon first weighting in the FMC by plotting the deflection and stress distribution of the main roof upon the first breaking using the Matlab software program based on the parameters calculated above. The results are shown in Figures 35.

As shown in Figure 3, the deflection of the main roof under fixed support on the four sides had a symmetric distribution about the goaf center. The largest deflection occurred in the center, and the deflection decreased gradually towards the periphery. The largest defection at the center of the goaf () reached 1.25 mm.

As shown in Figure 4, the stress in the main roof under the condition of fixed support on the four sides had an axisymmetric distribution about the goaf center with coordinates (), where the stress in the main roof reached its maximum ( MPa), indicating the largest compressive stress at this position. At the centers of the two long sides in the goaf, that is, at and , the stress in the main roof reached the peak negative value ( MPa), indicating the largest tensile stress in this position. Besides, the tensile stress already reached the ultimate tensile strength of the main roof. Therefore, tensile failure of the main roof occurred at the centers of the goaf’s long sides.

As shown in Figure 5, stresses in the main roof under the boundary condition of fixed support on four sides also had an axisymmetric distribution about the goaf center. At the center of the goaf , the stress in the main roof reached the largest positive value ( MPa), indicating the largest compressive stress at this position. At the centers of the two long sides in the goaf, i.e., and , the stress in the main roof reached the largest negative value ( MPa), indicating the largest tensile stress in this position. However, the tensile stress did not reach the ultimate tensile strength of the main roof. Given the above, the failure of the main roof was not caused by the stress .

3.3. Mechanical Analysis of the Periodic Breaking Characteristics of the Main Roof
3.3.1. Theoretical Analysis of the Mechanical Model for the Periodic Breaking of the Main Roof

The governing equation for the boundary condition of fixed support on three sides and simple support on one side before the periodic breaking of the main roof in the FMC can be written as follows:

Under the boundary condition of fixed support on three sides and simple support on one side, we derive the following deflection equation that satisfies the boundary conditions:

Substituting Equation (15) and into Equation (1), we can calculate the coefficient by the Matlab software program:

Substituting Equation (16) into Equation (15), we can derive the deflection function under the condition of fixed support on the three sides and simple support on one side for the main roof:

Substituting Equation (17) into Equations (3)–(6), we can derive the shear stress function and the bending moment function of the main roof as follows: where .

Substituting the abovementioned bending moment function into , we can obtain the corresponding stress function:

3.3.2. Periodic Breaking Interval of the Main Roof and Its Mechanical Characteristics

According to the stress function calculated in Section 3.3.1, at , both both reach the largest negative value (that is, are tensile stresses, ), indicating that both are located at the centers of the long sides with fixed support in the TPMM. Values of are derived as follows:

Thus, when in the main roof reached the tensile strength , tensile failure occurred periodically at the long sides with fixed support in the TPMM. By substituting the basic parameters of the FMC into the above equations, we estimated the periodic breaking interval of the main roof as  m.

Similarly, we plotted the deflection and stress distributions of the main roof upon periodic breaking using the Matlab software program (Figures 68). Then, the mechanical characteristics of the main roof upon periodic weighting in the FMC were analyzed.

As shown in Figure 6, the deflection of the main roof had an axisymmetric distribution about the goaf center () under the boundary condition of fixed support on the three sides and simple support on one side. We observed the largest deflection at the center, with a gradual decrease towards the periphery. In the position (, ) in the goaf, the main roof deflection reached a maximum of 0.36 mm.

As shown in Figure 7, the stress in the main roof had an axisymmetric distribution about the goaf center () under the boundary condition of fixed support on the three sides and simple support on one side. At the center of the long side with fixed support in the goaf , the stress in the basic stress reached the peak negative value ( MPa), indicating the largest tensile stress in this position. At this moment, the tensile stress already reached the ultimate tensile strength of the main roof. Therefore, tensile failure of the main roof occurred at the centers of the long sides with fixed support in the goaf.

As shown in Figure 8, both stresses and in the main roof had an axisymmetric distribution about the goaf center () under conditions of fixed support on the three sides and simple support on one side. At (, ) in the goaf, stress reached the largest positive value ( MPa), indicating the largest compressive stress in this position. At the center of the goaf long side , the stress in the main roof reached the largest negative value ( MPa), indicating the largest tensile stress in this position. However, the tensile stress did not reach the ultimate tensile strength in the main roof at this moment. It could be inferred that the stress caused no failure in the main roof.

4. Field Measurement of Ground Pressure Behavior in the FMC

In this section, the working resistance of hydraulic supports in the FMC was studied during the stoping period. The working resistance variations of the hydraulic supports were characterized. We analyzed the first and periodic weighting patterns in the coalface and verified the feasibility of the proposed TPMM for the main roof breaking. On this basis, we evaluated the applicability of the hydraulic supports used in the coalface, offering valuable clues for roof management in the coalface.

4.1. Field Measurement Plan for Ground Pressure

A total of 25 measurement stations were installed along the dip direction of the coalface to analyze ground pressure behavior at different positions in the roof (Figure 9). The measurement stations were located at the following hydraulic supports: #4, #8, #15, #23, #30, #38, #46, #53, #60, #68, #76, #83, #92, #98, #106, #113, #120, #128, #136, #143, #150, #158, #166, #174, and #183. The YHY60W(A) support pressure gauges for mines developed by the Xuzhou Shengneng Technology Co., Ltd. were installed in the corresponding measurement stations for continuous monitoring of the working resistance of hydraulic supports.

4.2. Analysis of the First Weighting Characteristics in the Coalface

Generally speaking, weighting in the coalface can be more accurately performed using the circulation terminal resistance (CTR) of hydraulic supports. However, after installing and debugging the special ground pressure-monitoring system, the FMC had already advanced by 140 m. Therefore, the CTR could not be timely and accurately captured by the electrohydraulic control system embedded in the hydraulic supports. For this reason, we analyzed the first weighting of the main roof in the coalface according to the actual roof caving, time-weighted resistance (TWR) of the hydraulic supports, and daily maximum resistance (DMR) of the hydraulic supports.

During the primary mining of the FMC, when the coalface advanced from the open-off cut for about 12 m, the immediate roof in the middle of the coalface collapsed. As the coalface further advanced to 20 m, the entire immediate roof had collapsed, filling the goaf. When the coalface advanced to about 40 m, clear breaking sounds could be heard from the main roof. Thus, the occurrence of the first weighting of the coalface was confirmed. During the primary mining of the coalface, the overall distribution of the daily average working resistance of the hydraulic supports is monitored and plotted in Figure 10. From the coalface’s head to its tail, the low-stress region, high-stress region, low-stress region, high-stress region, and low-stress region appeared successively. This was a saddle-shaped distribution, which was typical of the ultralong coalface.

The judgment criteria of time-weighted resistance (JCTWR) and the judgment criteria of daily maximum resistance (JCDMR) of some hydraulic supports in the coalface were calculated for the first weighting, as shown in Table 1. The variation curves of the working resistance of hydraulic supports are presented in Figure 11 (the red vertical line indicates the position of the first weighting in the coalface).

According to Table 1 and Figure 11, the first weighting interval of the main roof in the FMC was 38.5-45.9 m, with an average of 41.4 m, which agreed with the field observation. The main roof of the FMC was composed of hard and thick sandstone. Deep-hole presplit blasting has been performed for the open-off cut roof to prevent roof disaster caused by the extensive roof overhang. Therefore, the first weighting of the coalface had low intensity, and the safety valves of the hydraulic supports did not open extensively. The dynamic load factor of the hydraulic supports during the first weighting of the coalface was 1.04-1.18, with an average of 1.08, which values were relatively small.

4.3. Analysis of the Periodic Weighting Characteristics in the Coalface

Field measurement data were collected from the special ground pressure monitoring system installed in the hydraulic supports. Periodic weighting was analyzed as the coalface advanced from 140 to 400 m. The judgment criteria of time-weighted resistance (JCTR) and the judgment criteria of circulation terminal resistance (JCCTR) of some hydraulic supports were calculated; the results are shown in Table 2. The variation curves of the working resistance of some hydraulic supports are shown in Figure 12 (the red vertical line indicates the position of the periodic weighting in the coalface).

As shown in Table 2 and Figure 12, as the FMC advanced from 140 to 400 m (i.e., by a distance of 260 m), 10-13 events of periodic weighting occurred in total. The minimum, average, and maximum weighting intervals were 11.6, 22.0, and 38.7 m. The dynamic load factor of the hydraulic supports during the periodic weighting of the coalface was 1.03-1.09, with an average of 1.06. This dynamic load factor was relatively small.

4.4. Applicability Evaluation of the Hydraulic Supports Used

The operational resistance data of the hydraulic supports were analyzed as the FMC advanced from 0 to 400 m. The distribution of the daily average working resistance of the hydraulic supports is shown in Figure 13. From the head to the tail of the coalface, hydraulic supports’ daily average working resistance presented a “saddle-shaped” distribution similar to the primary mining of the coalface (low-stress region, high-stress region, low-stress region, high-stress region, and low-stress region). Along the advanced direction of the coalface, there were three types of regions, classified by their daily average working resistance, namely, low-stress region, transition region, and high-stress region. (1)Low-Stress Region (0-100 m). The overburden failure in the coalface was extremely insufficient within this region. The daily average working resistance of hydraulic supports was generally below 7000 kN and changed little. No significant coal wall spalling was observed during the coalface advance.(2)Transition Region (100-300 m). Within this region, the degree of overburden failure in the coalface grew with the coalface advance. The daily average working resistance of hydraulic supports gradually increased.(3)High-Stress Region (Beyond 300 m). The overburden failure in the coalface already became sufficient within this region. The daily average working resistance of hydraulic supports was generally above 8000 kN. Significant coal wall spalling was observed during the coalface advance.

The frequency distribution of working resistance usually characterizes the operational performance of hydraulic supports. A reasonable working resistance distribution of hydraulic supports is a near-normal distribution [30]. A higher share of low working resistance usually indicates a lower efficiency of the hydraulic supports. Besides, the rated working resistance margin is large under the above situation. By contrast, a higher share of the high working resistance indicates that the hydraulic supports mostly work at a high load. The rated working resistance is lower and cannot meet the actual support demand on-site. The frequency distribution of the working resistance of hydraulic supports was analyzed statistically as the FMC advanced from 0 to 400 m. It yielded 0.4% for the 0-2000 kN interval, 0.2% for the 2000-4000 kN interval, 1.1% for the 4000-6000 kN interval, 28.2% for the 6000-8000 kN interval, 65.3% for the 8000-10000 kN interval, and 4.8% for the interval above 10000 kN. According to these results, high working resistance (above 8000 kN) had the dominating share of 70.1%. This indicated that the hydraulic supports generally worked at high loads with a small margin of the rated working resistance.

5. Conclusions

(1)The proposed TPMM of the main roof breaking predicted that the first and periodic weighting intervals of the main roof caused by the internal stress in the rock strata were 40.6 and 25.0 m, respectively. The tensile failure first started in the goaf from the center of the long side with fixed support(2)Field measurements of ground pressure showed that the first weighting interval of the main roof in the FMC was 38.5-45.9 m, with an average of 41.4 m. The periodic weighting interval was 11.6-38.7 m, with an average of 22.0 m. The discrepancy between the model predictions and field data was quite small. The above results indicated that the TPMM was applicable to the particular geological conditions of the FMC(3)The estimation of the daily average working resistance of hydraulic supports exhibited a saddle-shaped pattern along the coalface dip direction, with the following sequence order: low-stress region, high-stress region, low-stress region, high-stress region, and low-stress region, judging by the daily average working resistance. This saddle-shaped distribution was typical of the ultralong coalface. Along the advanced direction of the coalface, the daily average working resistance could be classified into three types of regions by the low-stress region, transition region, and high-stress region(4)Analysis of the working resistance frequency distribution revealed that the ZY10000-16/32D two-column shield-type hydraulic supports could adapt to the geological conditions under which the FMC was mined. However, the hydraulic supports generally experienced high loads, with an insufficient margin of the rated working resistance (i.e., 10000 kN). Therefore, enhanced roof management in the presence of special geological structures (e.g., faults and goaf passages) during the coalface mining is strongly recommended

Data Availability

The data used to support the findings of this study are included within the article.

Conflicts of Interest

The authors declare that there are no conflicts of interest regarding the publication of this paper.

Acknowledgments

This work was financially supported by the National Natural Science Foundation of China (No. 51874278), the “Six Talent Peaks” Project of Jiangsu Province (No. JNHB-087), the Postgraduate Research and Practice Innovation Program of Jiangsu Province (Nos. KYCX21_2339 and KYCX21_2247), and the Independent Research Project of State Key Laboratory of Coal Resources and Safe Mining (No. SKLCRSM2020X04). We wish to thank the HCM in HMA, China, for supporting field help to conduct this important study.