Abstract

Extremely close coal seam groups are widely distributed in China, and the main mining method is downward mining. In the downward mining process of extremely close coal seam groups, the violent movement of overlying strata will cause the redistribution of surrounding rock stress. It not only produces stress concentration on the pillar but also causes the roof of the lower coal seam to be broken and difficulty in maintaining the mining roadway. In this study, the physical similitude modeling method and field observations were used to study the breakage and migration law of overlying strata in the downward mining of extremely close coal seams. Results show that in the process of mining upper coal seam, the first weighting step of the main roof is 37.5 m and the periodic weighting step is 12.5 m. The occurrence of strata separation is beneficial to the prediction of roof weighting. When the working face advances to 25 m, the rock stratum approximating a parallelogram of height 5 m does not collapse, and the working face is relatively dangerous. When mining the lower coal seam, the overall pressure of the working face is large, but the periodic weighting of the working face is not obvious. The first collapse step of the immediate roof is 15 m. When mining the upper and lower coal seams, the subsidence of the monitoring point increases significantly at 17.5 and 15 m, respectively. The roof collapse of the lower coal seam occurs 2.5 m ahead of that of the upper coal seam. The hydraulic value of the support, roof fall height, and sloughing depth in the entire working face reach the maximum at the coal pillar, and the extreme points at the coal pillar are relatively concentrated. This research provides some guidance for the safe and efficient mining of extremely close coal seams in the future.

1. Introduction

China produces and consumes the most amount of coal in the world [1]. The proportion of coal in primary energy production and consumption is approximately 70%. According to the British Petroleum Statistical Review of World Energy 2018 program, China’s coal output reached 3.523 billion tons in 2017, accounting for 45.6% of the world’s total output [2]. Therefore, coal is important for China’s economic development [3]. Recently, with large-scale mining, coal seams with good geological conditions have been gradually exhausted [4]. Owing to the large proportion of extremely close coal seams in China, the mining of extremely close coal seams is becoming more common for improving the utilization of coal resources [5]. Extremely close coal seams are closely spaced and interact with each other during mining. With the decrease in spacing, the interaction between coal seams will gradually increase [6]. When the distance of coal seams is extremely close, the roof integrity of the lower coal seam will be damaged by the mining of the upper coal seam. The area above the roof is the caving zone formed by the collapsed immediate roof [7]. Moreover, the remaining pillars in the upper coal seam can easily cause stress concentration [8, 9]. Consequently, the roof structure and stress environment in the mining area of the lower coal seam changed, and many new mine pressure phenomena occurred during the mining of extremely close coal seams [10].

Many engineering examples show that the violent movement of overlying strata will induce serious air leakages and water and gas accidents [11, 12]. Factors such as mining thickness, burial depth, and dip angle of coal seam are closely related to the law of overlying strata movement [13]. Because the process of strata movement is complex and no theoretical analysis method exists to satisfy engineering practices [14], physical similarity simulation is still the main method to study strata movement. This method overcomes the invisibility of mine pressure and overlying strata movement in field production, reflects the mechanical phenomenon visually, and can simulate the entire process in a short time. Based on physical similarity simulations, Huang et al. [15] studied the characteristics of overlying strata movement and strata behavior law in fully mechanized coal mining and backfilling longwall faces. Yuan et al. [16] proposed a new method for a similar material simulation experiment of steeply inclined upper protective layer mining and successfully applied it to the Nantong mining district. Niu et al. [17] constructed a similar physical model of coal rock to verify that a new method could be applied for the monitoring and early warning of coal and rock dynamic disasters. Zhang et al. [18] discussed the roof movement law of a fully mechanized mining face under a large dip angle through physical similarity simulations. Based on the engineering background of the Wuhushan coal mine, the law of overlying strata breakage and migration in the downward mining of extremely close coal seams was studied using the physical similitude modeling method. The current study can provide important guidance for the safe and efficient mining of extremely close coal seams in the future.

Recently, experts and scholars have performed relevant research and exploration on the mining system and safety technology of extremely close coal seams. Based on a mechanical model and FLAC 3D numerical simulations, Wu et al. [19] studied the stress distribution under a coal pillar and optimized the roadway layout. Based on the in situ monitoring of overburden failure, Ning et al. [20] proposed a statistical formula for predicting the maximum height of overburden failure induced by extremely close coal seam mining. Zhang et al. [21] discussed the relationship between pillar size and roadway stability, incorporating a strain softening model for pillars and a double yield model for goaf material. By considering horizontal, vertical, and tangential stresses, Yan et al. [22] and Yuan et al. [23] proposed a new method for calculating the stress distribution of coal pillars. Based on experimental research and the UDEC software, Zhang et al. [24] analyzed the relationship between mining sequence under water body and overburden failure degree. Zhang et al. [25] incorporated geotechnical considerations for concurrent pillar recovery in extremely close coal seams, where mining sequence, panel layout, and pillar size were considered. Based on the observations of surface subsidence and three-dimensional simulations, Yu et al. [26] and Zhu et al. [27] studied the relationship between upper coal pillar and lower working face. Based on the floor failure mechanics model, Zhang et al. [28] proposed a new method to monitor floor failure depth and successfully applied it to the Caocun coal mine in China. Liu et al. [29] deduced a formula for the analysis of floor stress distribution and roadway position in extremely close coal seams. Aiming at the large deformation and destruction of roadway in extremely close coal seams, Li et al. [30] proposed an asymmetric support scheme, which has been successfully applied in other mines. Based on the finite element method, Ghabraie et al. [31] and Khanal et al. [32] developed a new method that can accurately simulate the collapse of overlying strata and surface subsidence during multiseam mining. Based on the law of gas occurrence and outburst characteristics, Wang et al. [33] and Konicek and Schreiber [34] studied the sequence of coal seam mining, key protective seam mining technology, and gas control measures.

As described above, scholars primarily focused on the layout of mining roadways in the lower coal seam, gas control, mining sequence, and stress distribution of pillars and floor in the upper coal seam. However, studies regarding the breakage and migration law of overlying strata by physical similarity simulations are rare. Therefore, the law of strata breakage and migration must be studied to realize the safe and efficient mining of extremely close coal seams.

2. Engineering Background

2.1. Mining and Geological Condition

The Wuhushan coal mine, located in Wuhai city, Inner Mongolia autonomous region, China (Figure 1), covers a mining area of 12.6 km2. Coals 9 and 10 are extremely close coal seams, with a 0.45–5.02 m layer of sandy mudstone in the middle. The inclination and strike length of the working face are 130 and 400 m, respectively. A fully mechanized mining method was adopted. The average thickness of coal 9 was 3.2 m. The rock strata above coal 9 were mudstone of average thickness 9.4 m and medium sandstone of average thickness 6.0 m, while the rock stratum below coal 9 was sandy mudstone of average thickness 2.0 m. The average thickness of coal 10 was 2.2 m. Its roof was also the floor of coal 9, and the rock stratum below coal 10 was siltstone of average thickness 5.4 m. Figure 2 shows the generalized stratigraphy column.

2.2. Experiments for Determining Rock Mechanical Properties

To understand the rock mechanical properties better, the coal and rock mass in the field were processed into a certain shape using the ZS-100 fully automatic drilling machine, SCM200 double-end grinder, HJD-150A concrete sawing machine, and SC200 automatic core-taking machine. Figure 3 shows the coal and rock samples for use in experiments. Figure 4 shows the processing equipment. Uniaxial compression, splitting, and shear strength tests were performed on the samples to determine the mechanical parameters of coal and rock mass [35], as shown in Table 1.

3. Similar Material Simulation

3.1. Similarity Theory

A similar material simulation was performed based on the similarity theory. Geometric, time, and dynamic similarities must be considered between the model and prototype. Based on [18] and “dimensional analysis,” the dynamic similarity rate is presented as shown in equation (4). Meanwhile, Ren et al. [36] indicated that the dynamic similarity requires the force of the model and prototype at the corresponding point and time to be at a certain proportion to each other, and the main characteristics of force are reflected by compressive strength and bulk density in the experiment. Therefore, the compressive strength can be described as the dynamic similarity rate.

The geometric similarity rate of the model iswhere CL refers to the length similarity constant and Lm and Lp are the lengths of the similar material simulation model and prototype, respectively.

The time similarity rate of the model iswhere CT is the time similarity constant and Tm and Tp are the time of the similar material simulation model and prototype, respectively.

The density similarity rate of the model iswhere is the density similarity constant and ρm and ρp are the densities of the similar material simulation model and prototype, respectively.

The dynamic similarity rate of the model iswhere Cσ is the strength similarity constant and σm, σp, γm, and γp are the compressive strengths and bulk densities of the similar material simulation model and prototype, respectively.

According to the dynamic similarity rate formula, the compressive strength and bulk density of the strata in the model and prototype can be obtained (Table 2).

3.2. Overall Design of Physical Similarity Model

Based on the actual geological data of the fully mechanized mining face of the Wuhushan coal mine, fine sand, lime, and gypsum were selected as similar materials. The size of the test bench was 1800 mm (length) × 160 mm (width) × 1300 mm (height), and the plane stress model was adopted. The model building process is presented as follows [37]: (1) based on the calculation in Table 3, sand, lime, and gypsum were weighed and combined in a mixer. (2) The mixed material was paved evenly and compacted to maintain the required bulk density. Subsequently, mica powder was sprinkled on the strata to clarify the model bedding. (3) The other strata of the model followed the preceding steps until the required height was reached. (4) The weight of the overlying strata above the model was determined by adding the counterweight. (5) The model was dried naturally for five days.

Because the thickness and strength of the floor of coal 10 will not significantly affect the test, they can be simplified during building. The average height was 200 m from the actual working face to the surface. The thickness of the simulated overlying strata was 45.2 m, and the remaining height of 154.8 m was generated by the simulated pressure. The total height of this test was 1148 mm. The total excavation length was 1000 mm, and the length of each excavation was 50 mm. During building, the actual size of the strata should be adhered strictly. The specific amount of similar material is shown in Table 3.

3.3. Layout of Monitoring Points

To acquire the displacement variation of the overlying strata, displacement monitoring points were evenly arranged. An electronic theodolite with high precision was used to measure strata movement during mining. As shown in Figure 5, the monitoring points were evenly arranged above the roof of coal 10. Six detection lines were arranged in the model, namely, 2, 12, 22, 32, 42, and 52 cm from Coal 9. A total of 11 monitoring points were set on each line. A 15 cm × 10 cm grid design was adopted. To acquire data accurately during mining, a data collector was used to record the pressure data automatically; subsequently, the data are transmitted to a computer for analysis, as shown in Figure 6.

4. Mining Result Analysis of Upper Coal Seam

4.1. Mine Pressure Appearance Law

The open-off cut of the working face is 7.5 m. When the working face advances to 17.5 m, the first collapse of the immediate roof occurs. The collapse height is 2.5 m. As shown in Figure 7, when the working face advances to 25 m, the mining-induced fracture will not extend to the main roof. The collapsed strata form a two-part masonry beam articulated structure. The collapse height is 5 m, which is approximately 1.5 times the mining height. The upper minimum collapse range is 10 m. Because of its self-stabilization ability, the strata approximating a parallelogram did not collapse. However, when the working face advances to 27.5 m, the first collapse of the main roof occurs. The working face encounters the first weighting of the main roof. As shown in Figure 8, the overlying strata begin to separate when the working face advances to 32.5 m. The first periodic weighting of the main roof occurs when the working face advances to 37.5 m (Figure 9). The second periodic weighting of the main roof occurs when the working face advances to 50 m (Figure 10). The average periodic weighting step is 12.5 m.

From the discussion above, it is clear that with the advance of the working face, the first collapse of the immediate roof and the first and periodic weighting of the main roof will occur. Finally, the overlying strata will collapse in a large area. When the periodic weighting of the main roof occurs, the collapsed roof will exhibit a specific regularity. The collapse length of the overlying strata is the same, which is approximately equal to the periodic weighting step. From the working face up, the collapsed strata become more orderly and a stable articulated structure can be formed easily. This is because the strength of the first collapsed strata is low, and with the advance of the working face, the collapsed strata are gradually crushed by the strata above.

4.2. Analysis of Roof Subsidence

As shown in Figure 11, the subsidence of the overlying strata shows a certain regularity. The maximum and uniform subsidence is line 1, which is 2 cm from coal 9. The variation range is between 55 and 64 mm, which is close to coal seam thickness. The maximum subsidence of lines 2, 3, 4, 5, and 6 is 48, 45, 42, 36, and 26 mm, respectively. This shows that the subsidence of the overlying strata decreases with the increase in distance from the coal seam. This is because with the increase in distance from the coal seam, the probability of interaction increases between the collapsed overlying strata. In addition, some stable structures may be formed between the large strata. Consequently, the space between the strata and the dilatancy coefficient increase. Ultimately, the subsidence of the overlying strata is reduced. When the advancement distance of the working face is 17.5 m, the subsidence of the monitoring point increases significantly. This phenomenon is caused by a sudden roof caving in the gob, which is consistent with the physical similarity simulation results.

5. Mining Result Analysis of Lower Coal Seam

5.1. Mine Pressure Appearance Law

When the lower coal seam was mined, the overlying strata and the roof of the upper coal seam collapsed and recemented. Owing to the mining activity of the upper coal seam, the roof of coal 10 was damaged and generated some microfractures; additionally, the roof strength was low. Therefore, the first collapse of the immediate roof occurred at 15 m during mining. The roof collapse of the lower coal seam was 2.5 m ahead of that of the upper coal seam. In the subsequent mining process, no obvious periodic weighting of the main roof occurred, the roof falls with mining, and no obvious structure was formed, as shown in Figure 12.

5.2. Analysis of Roof Subsidence

As shown in Figures 11 and 13, the subsidence of the overlying strata has little effect on the outside of the mining area. However, the subsidence of the overlying strata above the mining area changed significantly. Most of the subsidence was concentrated between 80 and 104 mm. The subsidence of lines 1, 2, and 3 was approximately equal to the thickness of Coal 10. Owing to effect of mining, the original structure of the overlying strata was destroyed and the subsidence of the upward detection lines increased. Ultimately, the collapsed strata’s dilatancy coefficient would be reduced, and the rock mass further compacted. The increase in the sinking point of lines 4, 5, and 6 exceeded the coal seam thickness after mining Coal 10, and the maximum increase could reach to 227.9%. This was caused by the decrease in the coefficient of fragmentation and expansion. Meanwhile, it could be attributed to Coal 10 being directly excavated before the strata movement had stopped completely. As shown in Figure 13, when the advancement distance of the working face is 15 m, the subsidence of the monitoring point increases sharply and the roof collapses in a large area. The subsidence curve is approximately symmetric, with an unstable area on both sides and no obvious periodic weighting area in the middle.

6. Field Observation

6.1. Layout of Stations

To understand the law of mine pressure in extremely close coal seams, the mine pressure in 1001 working face of the Wuhushan coal mine was observed. Eight stations were arranged in the working face. The stations were densely distributed under the coal pillar and evenly distributed in other places, as shown in Figure 14. The stations were located at hydraulic support nos. 6, 20, 24, 28, 31, 56, 81, and 106.

6.2. Analysis of Observation Results

The collected hydraulic information of all supports was divided into four regions: the upper, middle, coal pillar, and lower regions. The pressure values of the hydraulic support and the change characteristics of the roof fall and sloughing in each region were considered, as shown in Figure 15.

From the data, it can be concluded that the first collapse steps of the immediate roof in the upper, middle, coal pillar, and lower regions of the working face are 16, 16, 14.5, and 14.3 m, respectively. Comprehensive analysis shows that the average first collapse step of the immediate roof in 1001 working face is 15 m and no obvious periodic weighting is shown, which is consistent with the physical similarity simulation results. The hydraulic value of the support, roof fall height, and sloughing depth in the entire working face reached the maximum at the coal pillar, and the extreme points at the coal pillar were relatively concentrated. Furthermore, maximum points appeared at the upper and lower regions, but the entire working face was not as large as the coal pillar. The hydraulic value of the working face was generally large, roof fall and sloughing occurred occasionally, and preventive measures must be improved.

7. Conclusions

In this study, the physical similitude modeling method was used to study the breakage and migration law of overlying strata in the downward mining of extremely close coal seams, which was verified by field observations in the working face. The conclusions are as follows:(1)In the process of mining upper coal seam, the first weighting step of the main roof was 37.5 m, and the periodic weighting step was 12.5 m. The occurrence of strata separation was beneficial to the prediction of roof weighting.(2)When the working face advanced to 25 m, the rock stratum approximating a parallelogram of height 5 m did not collapse, and the working face was relatively dangerous.(3)When mining the lower coal seam, the overall pressure of the working face was large, but the periodic weighting of the working face was not obvious. The first collapse step of the immediate roof was 15 m.(4)When mining the upper and lower coal seams, the subsidence of the monitoring point increased significantly at 17.5 and 15 m, respectively. The roof collapse of the lower coal seam was 2.5 m ahead of that of the upper coal seam.(5)The hydraulic value of the support, roof fall height, and sloughing depth in the entire working face reached the maximum at the coal pillar, and the extreme points at the coal pillar were relatively concentrated. The hydraulic value of the working face was generally large, roof fall and sloughing occurred occasionally, and preventive measures must be improved.

Data Availability

The data used to support the findings of this research are included within the paper.

Conflicts of Interest

The authors declare that they have no conflicts of interest.

Acknowledgments

This work was supported by the National Natural Science Foundation of China (No. 51974317), the Yue Qi Distinguished Scholar Project (800015Z1138), China University of Mining and Technology, Beijing, and the Fundamental Research Funds for the Central Universities (800015J6).