Abstract

Fully mechanized caving mining in extrathick coal seams is currently the most efficient mining method. To figure out the mining-induced stress in extrathick coal seams which is significantly different from the normal ones is essential for mining safety. This research is based on the No. 3106 working face of Lu-Zi-Gou coal mine in China, which extracts extrathick coal seams by using the comprehensive mechanized top-coal caving method. According to the elastoplastic theory, a multiforce superposition mining-induced stress model is constructed in this research. In this model, the mining-induced stress is generated by the side coal seam of the stope, including the load formed by the broken roof blocks in the caving zone, the unbroken rock layers in the side of the goaf, and the rock layers in the bending zone. To verify the model, electromagnetic wave CT detection and borehole endoscopy are used to test the coal’s electromagnetic wave attenuation and internal joints. The results indicate that the mining of extrathick coal seams forms internal and external stress fields in the side coal seams of the stope. As the depth increases, the mining-induced stress first increases and then decreases within the two areas. According to the above conclusions, the width of coal pillars and the location of roadways between adjacent working faces can be reasonably determined.

1. Introduction

In Northwest China, there are vast reserves of coal resources. Some coal seams’ thickness exceeds 20 m, which is an extrathick coal seam [1]. The extrathick coal seam provides geological conditions for maximizing coal mining efficiency. The mining of extrathick coal seams is an efficient way of extracting coal resources. In the past ten years, mines excavating extrathick coal seams have occupied an increasingly important position in China’s energy system, such as Tashan and Tong-Qi coal mines in Shanxi, Huating coal mine in Gansu, and Bu-Lian-Gou coal mine [2, 3]. A working face of the Bu-Lian-Ta coal mine in Inner Mongolia for mining extrathick coal seams has set a world record of 14 million tons of raw coal output per year [4]. The length of the working face in the extrathick coal seam is between 150 m and 300 m, and the fully mechanized top-coal caving mining method is mostly utilized [5]. However, the mining of extrathick coal seams will produce oversized goaf. In addition, the range of the influence of mining-induced stress, the maximum pressure, and the damage depth of the surrounding rock is much more extensive than that in regular coal seam mining [6]. Under the influence of mining-induced stress, the surrounding rock in the roadway, working face, chamber, and other areas deform massively which is hard to control. The above hazards seriously restrict the safe and efficient production of extrathick coal seams.

Due to the lack of long-term, stable, and effective technology and equipment for monitoring mining-induced stress and surrounding rock damage [7], physical simulation, numerical simulation, and theoretical modeling are often used to study the distribution and evolution of mining-induced stress under the condition of comprehensive top-coal caving mining in the extrathick coal seam [8]. The research results show that the mining-induced stress forms internal and external stress fields in the “caved zone” and “bending subsidence zone,” which bear the load generated by the overlying strata and the broken roof block in the goaf, respectively. The “internal” stress field is the primary source of mining-induced stress [9, 10]. It may also increase the ineffectiveness of supports and unpredictability at the goaf edges which may result in risky caving in-bye in the top-caving mining [11]. Kumar et al. and Mandal et al. conducted field research on Zero coal seams in Kajora area of which the thickness was more than 7.0 m. It showed that the mining of extrathick coal seams caused different surface subsidence and crack distribution [12, 13]. He et al. [14], Zhao et al. [15], and Zhu et al. [16] studied the internal and external stress fields under the mining of extrathick coal seams. They obtained the distribution and evolution of high-stress concentration in the internal stress field in the mining plane.

The cross-hole electromagnetic wave CT detection is one of the on-site geophysical testing methods which can be used to monitor mining-induced stress. The electromagnetic waves are generated and received by the devices in the two boreholes, and the properties of the rock strata are determined based on the attenuation of the electromagnetic wave [17, 18]. The electromagnetic wave CT technology has been widely used because it has advantages of high detection accuracy, fast speed, and convenient equipment operation, including exploration and examination of the coal mine geological structure, distribution of groundwater, damage area of the surrounding rock, stability of the intersection coal pillar, and other geological and disaster information [19, 20].

At present, the research of the distribution and evolution of mining-induced stress needs to be further studied under the condition of fully mechanized top-coal caving in the extrathick coal seam. Therefore, this research explores the use of cross-hole electromagnetic wave CT detection to examine the mining-induced stress of extrathick coal seams. In addition, in the process of fully mechanized caving mining in extrathick coal seams, the evolution process of mining-induced stress is analyzed.

2. Distribution of Mining-Induced Stress in the Extrathick Coal Seam

2.1. Geology and Mining Conditions

The Lu-Zi-Gou coal mine is located in Shanxi Province, China, with a coal output of 3 Mt/a. No. 3106 working face is one of the prominent working faces of this mine. It excavates a composite coal seam which included three coal seams, Nos. 2#, 3#, and 5#. The average thickness of the coal seam in this area is 25 m, the maximum value is 29 m, the inclination angle is 3°, and the buried depth is 350 m. The mining height is 25 m. No. 3106 working face adopts the comprehensive top-coal caving method, with a dip length of 150 m and the strike length of 1200 m. Hydraulic fracturing and blasting are used to weaken the top coal above the hydraulic support to improve the recovery rate of coal resources. The width of the mining roadway is 4 m. The adjacent working face of No. 3106 is No. 3107. There is a 30 m-wide coal pillar between these two working faces (Figure 1).

2.2. Mechanical Modeling

After the extrathick coal seam is mined, the breaking and moving boundary of the roof can be approximated as a straight line obliquely crossing the coal seam with angles and (Figure 2). In the vertical direction, according to rock strata’s movement, deformation, and fracture characteristics above the goaf, it can be divided into the caved, fracture, and bending zones [21]. The mining-induced stress formed at point M in the side coal seam of the stope is caused by the superposition of multiple forces. In the caved zone, the broken roof in the goaf attaches and applies a load on the unbroken rock strata, which are at the side of the stope. The lateral intact rock strata of the stope deflect under the action of its weight. The two ends of the curved zone are, respectively, located above the lateral coal body in the goaf and are firmly supported. The broken roof in the goaf, the unbroken roof on the side of the goaf, and the rock strata in the bending zone together form the mining-induced stress. The stress at point M can be analyzed based on stress superposition. According to the above analysis, a multiforce superposition mining-induced stress model (MSMS model) can be established. The MSMS model is based on elastic beams, and the force analysis is based on theoretical mechanics.

For point M in the side coal seam of the stope, the mining-induced stress is composed of the weight stress of the overlying strata, the additional stress which is generated by the broken roof in the goaf, and the stress which is produced by the bending zone:

According to the horizontal location of point M in the coal seam, can be calculated by the following equation:

The load generated by the broken roof in the side of the goaf can be calculated through the following formula [10, 11]:

The stress generated by Δqi at point M is , and the mechanical model to calculate is shown in Figure 3.

In Figure 3, , , and can be calculated through the following formula:

The vertical stress component generated by at point M is

After calculus and superposition calculation, the following formula can be obtained:

By superimposing the stress which is generated by all lateral unbroken rock layers, the calculation equation of is

The ends of the rock stratum in the bending zone are fixed. Near the fulcrum, the loads generated by deflection are as shown in the following formula:

The stress generated by at point M can be calculated as follows:

The calculation equation of can be obtained by accumulating the stress generated by the rock strata in the bending zone:

Bringing equations (5) and (8) into equation (1), the mining-induced stress calculation formula at point M can be obtained:

2.3. Mining-Induced Stress Distribution

Based on the geology and mining situations of the Lu-Zi-Gou coal mine, the values of the parameters in the MSMS model can be obtained (Table 1).

By bringing the values in Table 1 into the MSMS model, the distribution curve of mining-induced stress in the coal seam can be calculated (Figure 4). It shows that the range of mining impact area reaches 60 m. Moreover, there are two stress fields within this area: internal stress field (ISF) and external stress field (ESF). In the two fields, as the buried depth increases, the stress first rises and then decreases. The width of the internal stress field is 0–32 m, and the maximum value is 23 MPa which is 16 m away from the goaf. The range of the external stress field is 32–60 m with the peak value of 34 MPa, which is 44 m to the goaf.

According to the distribution of mining-induced stress, the internal and external stress fields of fully mechanized top-coal caving mining in extrathick coal seams can be divided into five areas. From the edge of the goaf to the inside of the coal seam, they are the low-stress area in ISF (I), the high-stress area in ISF (II), low-stress area at the junction (III), the high-stress area in ESF (IV), and the low-stress area in ESF (V). The distribution of mining-induced stress and regional division provide a basis for reasonably determining the location of adjacent roadways and the width of coal pillars. Area III is the best scope for the adjacent working faces’ mining roadways to minimize coal loss. Therefore, the width of coal pillars between No. 3106 and No. 3107 working faces is 30 m.

3. Field Test

3.1. Principle of Electromagnetic Wave CT Detection

The electromagnetic wave CT detection technology is based on the geometric ray theory. The electrical difference between the detection target and the surrounding medium is used to calculate the propagation characteristics of electromagnetic waves between the two boreholes and the absorption by the medium. Furthermore, the position, shape, and size of the target body, the degree of damage of the rock mass, and the difference in stress appearance can be further determined. As the penetration depth increases, the electromagnetic wave continues to attenuate. The reasons that affect the attenuation include the diffusion of electromagnetic waves, the reflection and scattering on the media interface with different electrical properties, and the absorption of coal and rock media. The electromagnetic wave attenuation coefficient is used to describe this phenomenon quantitatively. According to Maxwell’s equation, ßs can be calculated [22]:

The following formula can calculate the electromagnetic field strength in the medium:in which E is the field strength of the receiving point, V/m; is the initial radiation constant, V/m; β is the absorption coefficient, which is the absorption value of electromagnetic wave per unit distance in the medium, dB/m; is the direction factor function of the electromagnetic wave generating and receiving probe; is the distance between the detection holes, m.

According to the change of βs in the test area, particular analysis software is used in the electromagnetic wave CT detection to invert the distribution of the electromagnetic wave on the fluoroscopic section. Due to differences in the absorption of the electromagnetic wave by different media, when the electromagnetic wave passes through cave and fracture zones, βs is much larger than that in the surrounding rock. Therefore, the field strength consumption in these areas is small, manifested as a negative anomaly. The distribution of βs in the electromagnetic wave image can help identify the cracks in the coal body, the stress of the coal body, and other information.

In addition, the resistivity of coal and rock mass is the dominant factor affecting electromagnetic wave attenuation. Moreover, the resistivity and βs are negatively correlated. That is, the higher the resistivity, the smaller βs. The electrical resistivity of the coal seam can be influenced by the joints, metamorphism, water and gas contents, and stress status. The factors can also indirectly change βs. Chen et al. [23] analyzed the change of coal resistivity during uniaxial compression. The results showed that the high-stress concentration closed the initial cracks and cavities inside the coal specimen and increased the density.

Nevertheless, the resistivity of the sample decreased with the increase of the load before the peak strength. When load approached and exceeded the compressive strength, the formation of new cracks led to an increase in resistivity. Therefore, when the stress concentration in the coal seam is higher, βs is larger. Instead, if the mining-induced stress in the coal seam is small or broken, βs decreases. Accordingly, the electromagnetic wave CT detection technology can be used to quantitatively analyze the coal’s stress concentration and degree of damage.

3.2. Preparation before the Test

The operation of electromagnetic wave CT detection requires a transmitting hole and a receiving hole. The electromagnetic wave is generated by generators installed in the transmitting hole. The probes collect the electromagnetic wave in the receiving borehole after passing through the coal and rock mass between the two boreholes (Figure 5).

JW-6 underground high-frequency electromagnetic wave CT system (Figure 6) is used in this study, composed of the receiver, transmitter, antenna, mounting rod, and image postprocessing software.

The site was in Roadway A, No. 3107 working face. In different areas, the electromagnetic wave distribution in the coal seam was tested. Three sets of parallel cross-holes which are 110 m in front of the working face (station I), 150 m (station II), and 400 m behind the working face (station III) (Figure 7), were constructed in Roadway A. The label “R” stands for receiving drilling, and “T” stands for transmitting drilling. The coordinate system can be established by taking the boundary of No. 3106 working face as the origin and the inclination of the coal seam as the positive direction. The widths of the coal pillars and Roadway A are 30 m and 4 m.

The field test results showed that the physical properties of the rock strata were affected by many factors, especially homogeneity, water content, and integrity. When the electromagnetic wave CT and borehole endoscopy are used to observe the surrounding rock undermining, it would be better to operate the test in a more homogeneous surrounding rock. Also, the inverse analysis of the observation results determines the key to the accuracy of the test. At the same time, repeated testing is also a high-precision method. There was no clear structure and water seepage in the coal seam where the stations were located. The length of the boreholes was 16 m, and the spacing of the holes was 9 m. The operating frequency of the electromagnetic wave CT detection system was set to 28 MHz. It adopted fixed-point receiving and completed the working measurement mode. The borehole endoscope was used to observe the development of cracks on the surface of the “1T” and “1R” boreholes in station I to obtain the initial state of the surrounding rock of Roadway A. Due to the impact of mining, the boreholes in stations II and III were at the risk of collapse. In this study, borehole endoscopes were not used to observe these two stations.

3.3. Results

Due to the increase in mining-induced stress, the joints in the coal seam close, the resistivity reduces, and the degree of electromagnetic wave attenuation enhances. The electromagnetic wave attenuation coefficient is positively related to the mining-induced stress. The higher the pressure in the coal seam, the greater the electromagnetic wave attenuation coefficient.

3.3.1. Station I

The surrounding rock where station I was located had not been affected by the mining of the working face, but the roadway excavation caused damage to the surrounding rock. It can be seen from the deterioration of the borehole surface (Figure 8) that, within the depth of  0-1 m, the wall of the borehole is rough, and there are many cracks. It means that the coal seam is seriously damaged. In the depth of 2–4 m, the wall of the borehole becomes much smoother. However, there are still some cracks developed in this area because of the influence of tunneling. When the depth exceeds 5 m, the borehole wall is complete, and there are few cracks.

The change of βs in the coal seam where station I is located is shown in Figure 9. Along the borehole, βs varies between 2.5 and 5 dB/m. In the range of 34–38 m, that is, within 4 m of the roadway, the electromagnetic wave attenuation of the surrounding rock is lower than that in other areas. βs in some areas close to the borehole is also relatively small, indicating that the coal seam in this area has a higher degree of fracture. The results are consistent with the conclusions obtained by drilling endoscopy. Overall, βs of the unmined coal seam is small, and the stress concentration is low.

3.3.2. Station II

For station II and station III, the cross-holes in the coal seam of No. 3107 working face were located in the external stress field, and the boreholes in the intersection coal pillar were within the internal stress field. βs in the area where the two stations are located is significantly increased. The electromagnetic wave image obtained through the No. 2 cross-hole (Figure 10(a)) shows that βs in this area is generally less than 5 dB/m in the range of 34–48 m. Among them, the low EWT zone in the coal seam in the shallow depth area near the roadway is broken due to the impact of excavation. In the range of 38–46 m, βs varies between 5 and 10 dB/m; within 46–50 m, the scope of βs is 7–22 dB/m, and βs reaches its maximum value at 48 m. In summary, in the external stress field, as the depth of the coal seam increases, βs in the coal seam increases rapidly, which means that the concentration degree of mining-induced stress rises. The electromagnetic wave CT image obtained by the No. 3 cross-hole in the internal stress field is shown in Figure 10(b). The range of 14–20 m is the high EWT zone, and the peak value is about 18 dB/m, which is lower than the result of the external stress field. In the range of 22–30 m, there is an apparent low EWT zone, and the variation range of βs is 5–12 dB/m.

3.3.3. Station III

Station III is farther from the work surface than station II. In this area, the movement and deformation of the overlying strata reached a stable state. The features of electromagnetic wave CT images obtained by No. 4 and No. 5 cross-holes in station III are similar to those of station II. There is a high-stress concentration area in the external stress field in the range of 46–50 m, and the maximum value of βs is 14 dB/m (Figure 10(c)). 14–18 m is a high-stress concentration area in the internal stress field, and the maximum value of βs is 18 dB/m (Figure 10(d)).

Combining the observation results of electromagnetic wave CT and the calculation results of the MSMS model, it can be seen that both methods have confirmed the existence of internal and external stress fields. For the Lu-Zi-Gou coal mine, the maximum mining-induced stresses of the internal and external stress fields are located at 14–16 m and 44–50 m in the lateral direction of the stope, respectively.

4. Conclusion

This research is aimed to obtain the distribution of the mining-induced stress in the extrathick coal seam. The elastoplastic analysis, borehole endoscopy, and electromagnetic wave CT detection were utilized. The main conclusions obtained in this study are as follows:(1)The mining-induced stress in the side coal seam of the stope is formed by multisources which are the broken roof block in the caving zone, the side unbroken rock strata in the goaf, and the rock stratum in the bending zone. Accordingly, the multiforce superposition mining-induced stress model was established, and the calculation formula for mining-induced stress in the coal seam was derived.(2)Based on the mining and geological situation of the Lu-Zi-Gou coal mine, the distribution of the mining-induced stress in the extrathick coal seam is calculated. The results show that the mining-induced stress has formed internal and external stress fields in the coal seam beside the goaf. As the horizontal burial depth increases in the two fields, the stress in both fields increases first and then decreases. The maximum value of mining-induced stress and its position in the internal and external stress fields are 23 MPa, 16 m, and 34 MPa, 44 m, respectively.(3)Electromagnetic wave CT detection technology and borehole endoscopy are adopted to conduct on-site tests. The value of βs of the unaffected coal seam is 2.5–5 dB/m. Due to the mining influenced, βs increases significantly. In the internal and external stress fields, the maximum value of βs is about 18–22 dB/m. The calculated results of the MSMS model are in good agreement with the field data received by electromagnetic wave CT.

Nomenclature

H:Burial depth of the coal seam
h:Mining height of the working face
hi:Thickness of each rock stratum
L:Length of the working face
α, β:Angle between the moving and breaking boundary of the roof and the horizontal direction
H1, H2, and H3:Height of the caved, fracture, and bending zones, respectively
Pi:Load applied on the unbroken rock strata by the broken roof in the goaf in the caved zone
∆qi:Load generated by the weight of unbroken rock strata
Fi:Load generated by the bending zone
σ:Mining-induced stress
:Weight stress of the overlying strata
σ1:Additional stress generated by the broken roof in the goaf
σ2:Stress produced by the bending zone
γ:Unit weight of rock strata
ßs:Electromagnetic wave attenuation coefficient
K:Resistivity, Ω/m
ω:Angular frequency, 1/s
μ:Magnetic permeability, H/m
ε:Dielectric constant of the medium, F/m.

Data Availability

The data are available on request due to privacy/ethical restrictions.

Conflicts of Interest

The authors declare that they have no conflicts of interest.

Acknowledgments

Support for this work was provided by the Natural Science Foundation of Hebei Province (E2020402042, E2020402041, and E2021402061), the Open Fund Project of State Key Laboratory (SKLCRKF20-04), and the National Natural Science Foundation of China (51804093).