Abstract
The fracturing of hard roofs in different layers would result in complex ground pressure on the working face, such as supports collapsed and severe roadway deformation. However, the mechanism of the ground pressure induced by hard roof fracturing in different layers is still unclear. In the paper, a physical model of a 20 m extrathick coal seam mined with hard roofs existing was established based on the physical simulation similarity criterion. The overburden fracturing structure, abutment stress distribution, and failure characteristics of the coal body were monitored by a noncontact strain measurement system and resistance strain gauges, to reveal the mechanism of ground pressure induced by hard roof fracturing. Furthermore, on-site measurement was used to monitor and analyze the ground pressure affected by hard roofs in different levels. The results provide a theoretical basis for the control of ground pressure in extrathick coal seam mining with hard roofs.
1. Introduction
In coal seam mining, when hard roofs exist in the overburden, the breaking span is usually huge due to their great strength, thus would result in a strong ground pressure in the mining panel such as support failure and roadway deformation. Through numerical simulation and field measurement, He et al. and Wang et al. found that the large suspended area of hard roofs would easily cause a stress concentration and failure in advanced coal rock mass [1, 2]. When hard roof fracturing, the intensity of the energy release is higher due to the large breaking span. Ning et al. [3] studied the fracture energy of thick hard roofs by means of microseismic monitoring. The results showed that the high energy released by the breaking of thick and hard roofs is the main reason causing strong ground pressure. Bednarek and Majcherczyk [4] discussed the rock mass characteristics which influence the choice of support. Zhao et al. [5] studied the fracture characteristic of an extrathick and hard roof based on the long beam theory. Li et al. [6] studied the periodic breaking span of a thick and hard roof based on the Vlasov plate theory and the strain energy distribution characteristics of coal seam under a thick hard roof. Shen et al. [7] revealed the strong response for an entry influenced by overlying hard roof and the influence of the hard roof fracture position on the entry. Ju and Xu [8] analyzed the structural characteristics of overlying hard strata and ground pressure in the panel after a 7 m-thick coal seam mining. Xia et al. [9] studied the characteristics and mechanism of ground pressure in the mining panel under the combined action of hard roofs and a coal pillar. It was found in the above research that, because of the high strength and large overhang of hard roofs, the stress concentration is obviously on the hard roofs, which has a significant impact on the working face and roadway in the mining panel.
During the mining of an extrathick coal seam under hard roofs, due to the large mining thickness, the migration range of overburden is wide. The results showed that the failure height could reach 200 to 350 m during a 14 to 20 m-thick coal seam mining [10–12]. Field monitoring showed that the fracturing of hard roofs in a large space frequently resulted in the occurrence of ground pressure on the working face with a different strength and manifestation [13, 14]. The ground pressure on the working face shows the characteristics of “long and short duration”, of which the “long duration” interval was 30 to 60 m and the “short duration” interval was 12 to 20 m, the pressure intensity increased when the “long duration” occurred. However, there was no obvious regularity in the occurrence of strong ground pressure such as support failure or severe roadway deformation [15, 16]. Li et al. [17] found that the huge mining thickness of an extrathick coal seam results in a larger mobile space of high-level hard roofs, and the sliding instability of the high-level hard roofs would cause strong ground pressure. Singh et al. [18] studied the strata movement during underground mining of a thick coal seam. According to Li et al.’s research [19], the rotational movement of key strata in the overburden had a direct impact on the supports in the working face, and the dynamic loading induced by the instability of the Voussoir beam structure was analyzed. Xie and Xu [20] analyzed the influence law of different thickness and level of hard roofs on the peak value and influencing range of the abutment stress. Yu et al. and Chen et al. conducted field measurement on the deformation characteristics and stress distribution law of a roadway in the process of mining an extrathick coal seam under hard roofs [21, 22]. Guo et al. [23] studied the relationship between the support resistance and the overhanging distance of hard roofs, and criteria for reasonable selection of support resistance were given. Mondal et al. [24] monitored the strata behavior in the destressed zone of a shallow Indian longwall panel with hard sandstone cover using mine-microseismicity and borehole televiewer data. However, other scholars also invest in the hard roof research [25, 26]; the results show that the occurrence of the hard roof has a direct impact on the ground pressure in the working face and is complex.
In terms of the mining conditions of extrathick coal seams with hard roofs, failure and instability of hard roofs will occur gradually within a large space, which may further result in frequent occurrence of ground pressure of different strengths on the working face. Most of the aforementioned scholars focused on the structural characteristics of the overburden or worked from the perspective of mining pressure development but did not do in-depth research on the mechanism of ground pressure induced by the fracturing of hard roofs at different levels. According to the key stratum theory [27, 28], the lower level key stratum is close to the coal seam, and the breaking span is small. With increase in the occurrence level, the breaking span and the strength of the high-level key stratum increase accordingly. Therefore, due to the influence of the distance from the coal seam and the breaking span, the failure and instability of key strata at the low and high levels may have different effects on the ground pressure. In addition, the large breaking span and instability of high-level hard roofs may cause the occurrence of strong ground pressure on the working face and could have interaction with the lower key layers. A physical simulation experiment could reflect the structural characteristics of the overburden and its influence on the abutment stress of a mining panel directly and was considered an effective means to study the overburden structure and ground pressure in coal seam mining [29, 30]. Based on the exploitation of extrathick coal seams with hard roofs, for this paper, the research method of physical simulation was adopted to study the fracture law of key strata and the strength of ground pressure induced by the breakage and instability of strata in different levels.
2. Experimental Model
2.1. Geological Background
The carboniferous #3–5 extrathick coal seams were mainly mined in Tashan coal mine of the Datong mining area, with a thickness of 14 to 20 m. The method of top coal caving mining was adopted. The buried depth of the coal seam is 400 to 800 m, and the overburden is covered with multilayer hard roofs with the compressive strength of 60 to 120 MPa. Due to the large mining thickness of the coal seam, the fracturing of multiple layers in a large space resulted in a frequently strong ground pressure, accompanied by the supports crashed. The floor heave was severe within 10 to 40 m in advanced roadway, and the maximum floor heave was up to 0.8 m. The roof subsidence was up to 0.6 m, the shotcrete of the two sides of the roadway was seriously cracked, and the advanced individual props were seriously damaged, as shown in Figure 1.
Taking the #8216 working face in Tashan coal mine as an example, the Carboniferous #3–5 coal seam was mined, the average thickness of the coal seam was 16 m, and the buried depth and inclination of the coal seam was 418 to 522 m and 1 to 3°, respectively. The length and mining distance of the working face was 230 and 1500 m, respectively, and the coal seam was covered with multiple hard roofs.
Based on the statistics of the strong ground pressure during the mining process, as shown in Table 1, it can be found that the pressure on the working face was relatively strong within the initial mining range of 214 m. The roadway deformation and hydraulic props collapsed were serious in advanced 35 m, the support resistance increased obviously and was even accompanied by safety valve opening. In the later mining process, the strength of the ground pressure was relatively weakened, the influence range was reduced to 10 m in an advanced roadway, and the deformation was also reduced.
2.2. Basic Parameters of the Model
In order to study the effect of key stratum fracturing on the ground pressure in the working face, the method of physical simulation was used in the laboratory. Taking the #8101 working face as the background, the Carboniferous #3–5 coal seam was mined. The thickness and burial depth of the coal seam was 20 and 470 m, respectively. The coal seam dip was 1 to 3°, and the length and mining distance of the working face was 230 and 1500 m, respectively. The frame size of the physical model in the laboratory was m (). The geometric similitude ratio of the designed model was 150 : 1; the actual height of the model was 1.47 m, which simulated a height of 220 m. Materials including sand, calcium carbonate, and gypsum were used to simulate the overlying rock. The weight ratio, time ratio, and stress ratio were 1.667 : 1, 12.25 : 1, and 250 : 1, respectively. The actual working time of the working face is 16 hours per day, and the maintenance work time is 8 hours. The advanced distance is about 4 m per day. Therefore, the actual working time and excavation distance on the model can be calculated to be every 78 min and 2.6 cm each time. The thickness of the overlying unlaid strata was 272.65 m; then, the weight of the overburden rock is 6.816 MPa, and the magnitude of the compensation stress added to the upper part of the similar model is 0.027265 MPa. The basic parameters of the model are shown in Table 2.
Table 3 shows the physical and mechanical parameters of the coal and rock mass [31]. According to the similarity ratio, the matching parameters of each rock formation in the model are shown in Table 3. The rock layers No.32, No.27, No.22, No.16, and No.9 above the coal seam are calculated to be the key strata based on the key stratum theory [27, 28], which are the emphasis to be studied in this test.
2.3. Monitoring System
The monitoring system schematic is shown in Figure 2. Taking five key strata as the research objects (17, 45, 75, 107, and 146 m vertically above the coal seam), the thickness of each key layer was 9.44, 9.1, 10.12, 12.2, and 12.9 m, respectively. In order to monitor the influence of each key layer breaking on the distribution characteristics of the abutment stress in the coal body, a total of seven strain monitors were arranged at intervals of 30 cm in the coal seam. The first measurement point was 40 cm from the open-off cut of the working face. In order to truly reflect the structural characteristics of the overburden, a small simulated hydraulic support was placed on the working face.
(a) Physical model
(b) Strain monitor by noncontact Vic-2D
(c) Monitoring scheme of abutment stress
A noncontact strain monitoring system (Vic-2D) was used to monitor the overburden displacement timely, as seen in Figure 2(b). Black flecks were randomly sprayed onto the model surface, and a camera was used to capture movement of the flecks in real time as the overburden moved. The displacement of the flecks was obtained by postprocessing software and was inverted to the displacement-changes law of the overlying strata.
3. Analysis of Experimental Results
3.1. Structural Characteristics of the Overlying Strata
The mined interval of the simulated working face was 5 cm every 30 minutes. The structural characteristics and overburden displacement variation during hard roof fracturing at different levels were obtained, as shown in Figure 3.
(a) Fracturing of immediate roof
(b) Fracturing of KS1
(c) Fracturing of KS2
(d) Fracturing of KS3
(e) Fracturing of KS4
(f) Fracturing of KS5
When the working face was mined 60 m, the roof broke for the first time, with a breaking span of 42 m. The length of the broken block above the support was 7 m, and its weight was borne by the support. After the roof broke for the first time, the overlying unbroken strata flexed and sank by 20 cm, as shown in Figure 3(a).
When the working face was mined for 105 m, KS1 broke for the first time, with a breaking span of 82.5 m. The fracturing of KS1 caused synchronous rotation and collapse of lower strata and resulted in vertical variation of 0.42 m in the support. This increased the working resistance of the support, as shown in Figure 3(a).
When the working face was mined for 140 m, KS2 broke for the first time, and the initial fracturing step of KS2 was up to 105 m. The fracturing of KS2 caused synchronous failure and instability of the underlying strata, resulting in rotary movement, as shown in Figure 3(c). The intension was low during the KS2 breakage and the KS1 structure remained stable, which reduced the strength of the KS2 breakage. According to Figure 3(c), we could see that the vertical displacement of the support was only 0.22 m and the rock layer above the support showed the structural characteristics of a composite cantilever beam.
When the working face was mined 180 m, KS3 broke for the first time, with a breaking span of 140 m and a rotary subsidence of 7.5 m. The thickness of KS3 was 10.12 m. Due to its large thickness, breaking span, and rotary subsidence, the energy release intensity was relatively high during the breakage of KS3, thus resulting in synchronous rotary movement of the underlying strata. As shown in Figure 3(d), the drop-off in the support reached 0.9 m.
When the working face was mined to 210 m, KS4 in the high and far field broke for the first time. The breaking step distance was up to 170 m and the thickness of KS4 was 12.2 m. The high intensity of the KS4 fracturing resulted in the synchronous rotation of KS1, KS2, and KS3. The rotary movement of this large rock structure acted on the cantilever beam structure above the support, resulting in a 2.5 m vertical displacement of the cantilever beam structure and intense pressure on the support, as shown in Figure 3(e).
When the working face was mined for 225 m, KS5 broke for the first time. The breaking step distance was 175 m, and the thickness of KS5 was 12.9 m. In the case of the KS5 fracturing, although the intensity was relatively high, due to its long distance from the coal seam, the fracturing of KS5 did not cause instability rotation in all the underlying strata. As shown in Figure 3(f), the cantilever structure of KS1 maintained its own stability, thus having a certain protective effect on the working face. However, due to the high breaking intensity of KS5, the vertical displacement at the immediate roof was still up to 0.48 m, and the breaking impact on the support was slightly lower than that of KS4.
Based on the comprehensive analysis of the overlying strata fracturing, it was found that due to the large mining thickness of the coal seam and the wide range of overlying rock migration, the key strata of KS1 to KS5 all broke, and the key strata at different levels had different influences on the working face. The fracturing of KS1 and KS2 in the near field formed a cantilever beam and masonry beam structure, which resulted in a slight ground pressure on the working face. When the far field KS4 broke, due to its great thickness and wide breaking span of 170 m, the ground pressure on the working face was most intense. The intension of the KS5 fracturing was also high with a breaking span of 175 m, but due to the long distance from the coal seam and the stable structure in the lower strata, the effects of KS5 fracturing on the working face was reduced compared to KS4, and the working face showed slightly strong ground pressure.
3.2. Variation of Abutment Stress in the Coal Body
The above analysis was made based on the overburden structural characteristics and displacements when the key strata broke at different layers. The influence of the occurrence and fracturing of key strata on the stress distribution in the coal body could be gained based on the monitored strain value; the calculation process is shown in the following. The strain gauges in the coal seam were numbered 1-7 along the mining direction of the working face, as shown in Figure 2(c). Four strain gauges were selected to decide the elasticity modulus, as shown in Table 4. Thus, based on the elasticity modulus and monitored strain, the stress variation at each monitor point were gained, as shown in Figure 4.
(a) Main roof fracturing for the first time
(b) KS1 fracturing for the first time
(c) KS2 fracturing for the first time
(d) KS3 fracturing for the first time
(e) KS4 fracturing for the first time
(f) KS5 fracturing for the first time
The abutment stress distribution during the first fracturing of the main roof is shown in Figure 4(a), where the No.1 monitoring point is next to the working face. After the main roof breaks, the stress at the No.1 monitoring point decreases from 24 to 14 MPa, indicating that the stress concentration on the coal wall is reduced to a certain extent. When KS1 broke for the first time, the No.2 monitoring point is close to the working face. Due to the rotation of KS1, the stress at the coal wall increases from 24 to 28 MPa, as shown in Figure 4(b). During the first breaking of KS2, the No.3 monitoring point is 15 m from the working face, and the strata fracturing has little impact on the abutment stress, as shown in Figure 4(c). When KS3 first broke, the No.4 monitoring point is 15 m from the working face. As a result of KS3 breaking, the stress value at No.4 monitoring point increases from 36 to 43 MPa. Meanwhile, at monitoring points No.5 and No.6, the stress decreases slightly, as shown in Figure 4(d), indicating that the stress concentration in the advanced 60 m is reduced. When KS4 breaks, monitoring points No.5 and No.6 are 30 m and 75 m, respectively, away from the working face. The stress values at monitoring points #5 and #6 are reduced from 31 and 34 MPa to 29 and 32 MPa, respectively, as shown in Figure 4(e). When KS5 broke for the first time, monitoring points No.5 and No.6 is 15 m and 60 m, respectively, away from the working face. The stress change law is similar to that for KS4, as shown in Figure 4(f).
The influence of hard roof fracturing at different layers on the peak stress and the range of the abutment stress in the advanced coal body were obtained statistically (see Table 5).
As shown in Table 5, with the development of the overburden caving height, the abutment stress and the influence range in advanced coal seam increase accordingly. The peak stress in the advanced coal body increases slightly (30 to 34 MPa), and the advanced range of influence increases greatly. Before the critical fracturing of KS4, the range of influence of the advanced abutment stress reaches 120 m.
When the key strata broke, the strata rotations squeezed the coal body near the working face, causing the stress on the coal body to rise. Affected by the strain monitoring layout interval, it was impossible to accurately calculate the range of influence in the advanced coal body affected by the fracturing of each key layer, but we could still find from Table 5 that the fracturing of KS4 resulted in a strong compression zone of 15–30 m in advance of the working face. With the increase of the distance between the key stratum and coal seam, the buffer action of the underlying strata was strengthened, and the strong compression area in the advanced coal body was weakened accordingly (e.g., KS5). As the distance between the key stratum and the coal seam decreased, the breaking span and the range of the advanced compression zone decreased accordingly (e.g., KS3).
3.3. The Advanced Failure Zone
Based on the noncontact strain monitoring system, the strain variation characteristics of the advanced coal body were gained through the strain measurement line to further study the range and degree of influence on strong compression areas caused by the fracturing of each key layer.
Strain measurement lines were arranged at the interface of the coal seam and immediate roof to obtain strain variation of the coal body after the fracturing of KS1 to KS5, as shown in Figure 2; the results are shown in Figure 5. The abscissa represents the location of monitoring points along the mining direction, the leftmost side of the model is the coordinate zero point, and the coordinate axis increases along the mining direction of the working face accordingly. The length of the model is 375 m.
(a) Strain of coal body during KS1 fracturing
(b) Strain of coal body during KS2 fracturing
(c) Strain of coal body during KS3 fracturing
(d) Strain of coal body during KS4 fracturing
(e) Strain of coal body during KS5 fracturing
The boundary pillar of 15 m was retained in the model. When the working face was mined to 105 (corresponds to abscissa m), KS1 broke, and a statistical strain curve of the coal body before and after KS1 was obtained, as shown in Figure 5(a). The simulated support is 15 m long. From Figure 5(a), we can find that the coal body is compressed in the advanced 5 m, that is, the impact range of the KS1 fracturing on the advanced coal body is about 5 m. Similarly, the strain curve during KS2 fracturing is obtained as shown in Figure 5(b). Compared with KS1, the compressed zone in the advanced coal body caused by KS2 fracturing slightly increases, reaching 6 m. As shown in Figure 5(c), the compressed range in the advanced coal body during the fracturing of KS3 reaches 10 m. According to the coal body failure criterion obtained from the test, the fracturing of KS3 resulted in 3 m of advanced coal body failure.
When KS4 broke, due to its high breaking strength and strong impact, the strain on the coal body was greatly increased, as shown in Figure 5(d). The maximum strain variable in the advanced coal body is up to 0.07, which could easily cause the crash of supports and the advanced influence range of 20 m. The coal body is severely compressed in the advanced 20 m. Similar to KS4, the action of KS5 fracturing on the advanced coal body is also relatively serious, but the range of influence and action intensity are lower than that of KS4, as shown in Figure 5(e). The range of influence from the KS5 fracturing on the advanced coal body is 16 m.
4. Discussion
According to Figure 5, the range of influence and degree of overburden-key-strata fracturing in the advanced coal body are statistically obtained as shown in Figure 6. It can be seen from Figure 6 that, with an increase of the occurrence level of key strata, the impact strength of each key layer fracturing on the coal body increases accordingly. The fracturing of KS4 has the strongest effect on the advanced coal body, with the maximum strain variation of the coal body reaching 0.036 and the influence range reaching 20 m. With further increase of the distance between the key layer and the coal seam, taking KS5 as an example, the action strength of KS5 on the advanced coal body decreases, the change of strain on the coal body is 0.013, and the range of influence decreases to 16 m. As the distance between the key strata and the coal seam gets closer, the intensity of the key strata fracturing will also decrease due to decrease of the breaking span of the key strata. Taking KS3 as an example, the maximum change in the strain on the advanced coal body caused by the KS3 fracturing is 0.005, and the range of influence is 10 m.
Combined with the above studies, we can find that the impact strength of the fracturing of key strata on the mining panel pressure is directly related to the breaking span, occurrence location, and distribution of the underlying strata. The physical simulation studies show that KS4 fracturing had the strongest effect on the mining panel pressure, followed by KS5. The ratio of the distance between KS4 and KS5 from the coal seam and the coal seam thickness are 5.35 : 1 and 7.3 : 1, respectively.
From Figure 6, we can see that the compression induced in the coal body by the fracturing of KS4 and KS5 are greater than that of KS1–KS3. This is mainly because during the fracturing of KS1 to KS3, the lower bed-separation space and key layer rotation angle are big. Although the impact velocity of the key strata is larger, the breaking step of key strata is relatively small. In addition, the key strata are more likely to act on the goaf caved rock mass, as shown in Figure 7(a), which has a weak impact on the mining panel pressure.
(a) Ground pressure during near field strata breaking
(b) Ground pressure during far field strata breaking
As the overburden fracture height develops, the bed-separation space under the far field key strata decreases, and the impact velocity after the breaking of key strata decreases accordingly. However, the kinetic energy after a far field key strata break is huge due to its great thickness, high strength, and large breaking interval, coupled with small rotary space. The kinetic energy of strata fracturing mostly acts on the lower overburden, resulting in a synchronous rotary movement of the lower overburden rock strata. This causes strong ground pressure in the advanced coal body and working face, as shown in Figure 7(b).
With further increase of the vertical distance between the key strata and the coal seam, the bed-separation space is further reduced, and the energy transferred from the key strata fracturing to the coal seam is weakened accordingly. Strata fracturing does not necessarily cause instability and synchronous movement of the underlying strata. According to the above physical simulation analysis, in this geological condition, the KS4 fracturing (in the far field) has the most serious impact on the working face, followed by KS5. The vertical distance from the coal seam of KS4 and KS5 are 107 and 146 m, respectively, and the ratio of the vertical distance to the coal-seam thickness is 5.35 and 7.3, respectively.
At the same time, we also used the method of field measurement to study the ground pressure in the mining of extrathick coal seams with hard roofs [13, 15]. The thickness of the coal seam was 19 m, and the method of top coal caving mining was adopted. The strata movement measurement points were arranged in key strata at different levels (22, 51, and 104 m away from the coal seam). The thickness of the three key strata was 12, 9.8, and 23 m from the bottom up. Meanwhile, the resistance characteristics of the working face support were recorded in real time, as shown in Figure 8. The monitoring results showed that the support resistance in the working face increased with fracturing of the two key strata, which were 22 and 51 m away from the coal seam at the low level, and the dynamic load coefficient of the support was 1.15 and 1.34, respectively. The pressure duration was 7 h and 16 h, respectively, and the working face had no obvious indications of strong ground pressure. When the 23 m-thick key layer (which was 104 m from the coal seam) broke, the #35–95 supports in the working face were crushed, and the dynamic load coefficient of the supports reached 1.54. The pressure duration reached 43 h, and the ratio of the distance between the highest key layer and the coal seam, to the coal seam thickness, was 5.47.
Thus, due to the great thickness of the extrathick coal seam, the range of the overburden movement was huge and the key strata at all different levels broke. Among them, the hard roof in the far field had the largest breaking span and highest energy intensity. This resulted in the synchronous rotation of the underlying strata, which is the main reason of the occurrence of strong ground pressure on the working face. Based on the above studies, it was found that thick, hard strata, with the ratio of the distance from the coal seam to the coal seam thickness between 5.3 and 7.3, have the most serious impact on the working face. That is, when a 20 m-thick coal seam is mined, the thick and hard roof, which is within the range of 106 m to 146 m away from the coal seam, has the greatest impact on the mining pressure of the working face. The results provide a theoretical basis for guiding the selection of the hard roof control layer in a large space and improve the reliability of the hard roof control. The results show that the high-level hard strata are the main factor causing strong ground pressure in the stope, which breaks through our idea of control ground pressure by low-level rock weakening and improves the reliability and accuracy of ground pressure control. For the control of high-level hard strata, traditional technical means cannot reach that height. Thus, our team developed a technical method of ground fracturing to control high-level hard strata and achieved success. This will be introduced in the follow-up researches.
5. Conclusions
(1)Based on physical simulation criterion, a physical model for the mining of an extrathick coal seam with hard roofs was established and the methods of noncontact strain monitor and abutment stress monitoring were used to study the structure characteristics and the corresponding ground pressure of hard roofs at different levels(2)The research showed that the lower key strata mostly collapse into goaf, thus having less impact on the working face and advanced coal body. During the thick, hard key stratum in the far field breaking, the large fracturing span and high strength result in instability and synchronous rotation of the underlying strata, thus causing the supports crashed and advanced coal body damaged, which is the main reason for the occurrence of strong ground pressure(3)Combined with the physical simulation and field measurement studies, it is found that the characteristics of ground pressure vary with the levels at which the key strata occur. During the mining of an extrathick coal seam with hard roofs, thick and hard strata with a ratio (of the distance from the coal seam to seam thickness) of 5.3–7.3 have the greatest influence on the working face
Data Availability
The data and material are transparent and included in the paper.
Conflicts of Interest
The authors declare that there is no conflict of interest.
Acknowledgments
We are grateful Yun Zhang, Lixin Lan, Baiyi Li, and Hao Yan for their assistance in physical simulation experiment. This work was supported by the State Key Research Development Program of China (grant numbers 2018YFC0604500, 2018YFC0604506), China Postdoctoral Science Foundation (grant number 2019M651080), applied basic research project of Shanxi Province (grant number 201901D211030), Scientific and Technological Innovation Programs of Higher Education Institutions in Shanxi (STIP) (grant number 2019L0208), Major Program in Shanxi Province (grant number 20191101015), and Independent Research Project of State Key Laboratory of Coal Resources and Safe Mining, CUMT (SKLCRSM18X006).