Abstract
In order to realize the safe and efficient mining of the short-distance isolated island working face under the deep goaf area, the 120502 isolated island working face of Liuzhuang Mine was taken as the engineering background. The method of combining numerical simulation and field measurement were used comprehensively to systematically simulate and study the spatial evolution of the stress field, plastic strain field, and fracture field of coal rock during the mining process. The leading support pressure and the vertical displacement of the roof in the overlapping section and noncoinciding section of the isolated working face and the goaf above were measured on site. The results are that the peak value of the advanced support pressure of the overlap section and the nonoverlapping section is 10 m before the coal wall of the working face; the advanced support pressure of the nonoverlapping section is 33.3 MPa, and the vertical displacement of the roof is 300 mm. The advanced support pressure and the vertical displacement of the roof in the noncoincidence section were significantly higher than those in the coincidence section of 18.2 MPa and 210 mm. The results are consistent with those predicted by numerical simulation. This provides theoretical support for the safe mining of the 120502 isolated island working face in Liuzhuang Mine and, at the same time, provides a reference for the study of similar working faces in other domestic mining areas.
1. Introduction
As the mining of coal resources in our country continues to move to the depths, due to the influence of geological conditions and the need for safe mining in the mines, it is inevitable that deep mine isolated islands are formed [1–4]. The coal rock masses in many mines have a certain tendency to impact, so the mining of isolated island working faces will encounter more complicated safety problems. Due to the open space on both sides of the isolated working face, the surrounding rock of the roadway is under a high-stress environment for a long time. As a result, the pressure of the mine becomes obvious, and the surrounding rock is deformed and damaged seriously [5–7]. It is difficult to support and maintain the two lanes and the open-off cut. The adjacent rock mass is in a rheological state after being destroyed and affected by dynamic pressure. Therefore, the behavior of mine pressure is complicated, which affects the safe advancement of the working face [8]. To this end, researchers have conducted a lot of research on the behavior of rock pressure on the isolated island working face and rock burst problems.
Liu et al. [9] and Lu and Shi [10] analyzed the support pressure distribution characteristics of the ultralong-isolated island working face and the rock pressure behavior when the island working face crosses the fault structure. Wang et al. [11] established a thin plate mechanics model and analyzed the distribution behavior of the roof bending moment of the isolated island working face. Wang et al. [12, 13] studied the behavior of rock pressure in the short-walled coal pillar working face of an isolated island. Zhao et al. [14] studied the behavior of rock pressure in a typical coal seam island face in western China. Yang et al. [15] comprehensively considered overlying rock movement, geological structure, mining and space, and other factors and classified isolated working faces based on scour prevention. The isolated island working face is divided into six types: full mining, insufficient mining, full-insufficient mining, three-dimensional, “recessive,” and compound type. In order to carry out the prevention and control of rock burst in different types of isolated island working faces, Dou [16, 17] proposed the detection and prediction and control technology of rock burst hazard in isolated island working face [18–20]. Wang et al. [21] analyzed the roof rupture height of the isolated island working face and the temporal and spatial distribution characteristics of microseismic events in the “square” area of the working face. Zhang et al. [22] used elastic plate theory to analyze the elastic energy value released during the fracture of the hard roof and used FLAC3D numerical simulation to divide the high-stress area during the mining face. At the same time, according to the fractal theory, the internal correlation between the time distribution and spatial distribution of microseismic events and rock bursts was analyzed.
The above research studies have analyzed the appearance behaviors of rock pressure in different types of isolated island working faces and the prevention and control of rock bursts and have achieved good results. However, there is little research on island working faces under some complex conditions, especially those under the goaf. Based on the engineering background of 120502 isolated island working face in Liu Zhuang Coal Mine, this paper analyzes the behavior of underground pressure on the isolated island working face under the goaf of the deep close coal seam group through numerical simulation and field measurement.
2. Project Overview
The research working face is 120502 isolated island working face of Liu Zhuang Coal Mine. 120502 isolated island working face is located in the first level and second mining area. The south of the working face is the goaf of 120503 working face, the north is near the goaf of 120501 working face, and the above is the 6-1 and 8 coal goaf. Also, the vertical distance between the overlying 6-1# coal seam goaf and the 5# coal seam is 18 m, which is an isolated island working face under the short-distance goaf. The specific spatial position relationship is shown in Figure 1.
The 120502 working face is a monoclinic structure with an average dip angle of 13.0°, the coal thickness of the working face is 0.98~5.58 m, and the average thickness is 4.05 m. 120502 working face 5 coal false roof is not developed, the direct roof is mudstone, the upper part contains a small amount of sandy content, and the average thickness is 8.73 m. The main roof is fine sandstone. The geological histogram of the coal roof and floor of the working face is shown in Figure 2. This working face is a low-gas working face, and the coal seam is easy to combust spontaneously.
3. Numerical Simulation
In order to simulate the stress and deformation characteristics of the overlying strata in different spatial layers during the mining process and to reproduce the spatial change process of the entire overlying strata movement and deformation, this paper uses FLAC3D software to carry out a numerical simulation.
3.1. Model Building
A numerical model is established on the fully mechanized mining face of the 120502 isolated island with main goafs on the left, right, and upper sides, The model is 550 m long, 360 m wide, and 200 m high. The Mohr–Coulomb constitutive model is adopted for coal and rock mass, and its physical and mechanical parameters are shown in Table 1.The characteristics of the top and bottom of the working face are shown in Figure 2. Displacement boundary conditions are applied around the model, and the upper part of the model is free boundary conditions. A vertical stress of 25.1 MPa is applied at the top of the model, a stress of 30.2 MPa is applied at the bottom, and the surrounding stress is 42.3 MPa According to the geological and ground stress test report of the working face, considering the elevation of the working face and the weight of the model, a vertical stress of 9 MPa was applied to the top of the model. The elastoplastic model is adopted, and the Mohr–Coulomb failure criterion is selected. The specific mechanical parameters of the main coal and rock layers are shown in Table 1.
3.2. Simulation Scheme
In the process of numerical simulation, in order to correctly simulate the stress, displacement, and fracture deformation characteristics of the surrounding rock of the stope caused by the stopping of the isolated working face, the model reaches the initial stress balance state under the given ground stress and boundary conditions. According to the sequence of the stope, the 120501, 120601, and 120503 working faces were excavated to construct a “three-sided” isolated island working face. The stopping of the 120502 mining working face is performed step by step. The stopping distances are, respectively, 80 m, 130 m, 180 m, and 240 m, and each excavation is designed to be 3000 steps. The distribution characteristics of stress, displacement, and plastic zone during the mining process are studied from the two directions of working face inclination and alignment direction. And, the distribution cloud map is quantified, so as to conduct in-depth research on the deformation characteristics of the overlying rock fractures during the stopping process of the isolated working face.
3.3. Simulation Result Analysis
3.3.1. Analysis of Evolution Characteristics of Mining Stress Field
After the model starts excavating, the stress distribution characteristics of the central coal rock mass in the inclined direction at the coal wall position at different advancing distances are shown in Figure 3. The pressure relief area in the inclined direction moves upward in an arch shape as the mining distance increases and gradually merges into one body. The coal pillar on the right side of the working face has a high degree of stress concentration, and the stress concentration factor has been increased from 2.54 to 4.41. Because the coal pillar on the left is in the pressure relief range of the upper goaf, and the stress concentration is relatively weak.
(a)
(b)
(c)
(d)
Figure 4 shows the vertical stress distribution cloud diagram in the middle of the goaf in the alignment direction when advancing at different distances. It can be seen from Figure 1 that the 120502 working face and the 120601 working face are staggered in the horizontal direction. At the square of the working face, affected by the mining of 5th coal, the vertical stress distribution range of the floor rock is enlarged. The concentrated stress inside the coal body at the open-off cut gradually shifted from the 6th coal to the 5th coal, showing that the concentrated stress at the end of the 6th coal weakened and the concentrated stress at the 5th coal end increased. The floor is affected by the upper coal mining, and the front of the working face is in the mutual influence area between the pressure relief zone and the concentrated stress zone, and the concentrated stress is not large. When the working face advances to 240 m, the pressure at the end of coal 6 is released. The stress is concentrated at the end of the lower layer coal 5, which increases the concentrated stress and increases the degree of damage.
(a)
(b)
(c)
(d)
The overburden strata have basically similar changes in the alignment and inclination directions. The stress concentration area and the stress reduction area are also basically the same, and they are continuously moved and expanded in an arch shape to form connections with other mining areas. The main goafs on the left, right, and above all have an impact on the deformation of the overlying rock in the goaf, but the upper main goaf has the greatest disturbance to it.
It can be seen from Figure 5 that, in the alignment direction, in the section overlapping with the upper goaf, the stress concentration areas are behind the open-off cut and in front of the working face. When the mining distance does not exceed 240 m, the distribution in front of the working face will first decrease and then increase. When the mining distance of the working face reaches 240 m, the distribution of “horizontal three areas” appears. The support pressure in front of the working face is increasing, and the distance between the peak support pressure and the coal wall is stable at about 10 m. In the noncoincident zone, it can be seen from Figure 6 that the stress concentration area remains unchanged, and the peak stress concentration is also equivalent under different advancing distances. Compared with the coincident section, the peak support pressure in the front of the working face under the same advance distance is larger. Comparing the two sections, it can be seen that the overburden stress distribution characteristics of the stope in the direction of the coincident and noncoincident sections are basically the same. The areas of stress concentration are both behind the open-off cut and in front of the working face, and the coal rock masses nearby are subject to greater stress, resulting in compression and destruction.
3.3.2. Analysis of Change Characteristics of Mining Displacement Field
Figure 7 reflects the vertical displacement change curve in the direction of the 120502 and 120601 working faces at different advancing distances in the overlapping sections. The direction of the vertical displacement in the section that coincides with the upper goaf in the alignment direction changes from downward to upward, and the vertical upward displacement peak remains stable. The change curve in the noncoincident zone is shown in Figure 8. The upper main golf area has a great influence on the vertical displacement. The peak value of the vertical displacement increases with the increase of the mining distance, and the vertical displacement of the overlying strata in the mined-out area is downward during the whole advancement process.
From the vertical displacement curve during the advancing process of the working face, it can be seen that the vertical displacement of the overburden under different advancing distances forms the largest sinking area in the middle of the mined-out area, resulting in stress concentration in the area behind the open-off cut and near the working front.
It can be seen from Figure 9 that, during the mining process of the working face, the horizontal displacement increases continuously with the advance of the working face in the alignment direction. The changing graphics show curved sinking and symmetrical, with the center of symmetry in the middle of the goaf. The horizontal displacement changes strongly behind the open-off cut and in the goaf area in front of the working face. The horizontal deformation of the overlying rock is relatively large, and the overlying rock is subjected to the combined action of horizontal stretching and compression, resulting in more fully developed cracks.
3.3.3. Analysis of Change Characteristics of Overlying Rock Plastic Zone
Figure 10 shows the development of the plastic zone in inclined direction when the working face is advanced at different distances. The closer to the working face, the greater the failure of coal and rock mass. H is the development height of the fracture zone, which is 101 m, 119 m, 126 m, and 130 m, respectively. When the upper part of coal reaches a certain height, it basically maintains a stable state. With the advance of the working face, the working face is in the state of nonuniform expansion and dynamic development and gradually intersects with other plastic areas of the goaf. The plastic zones of multiple working faces are interconnected, the rock mass above the roof is in the tensile and compressive stress zone, and the surrounding rock cracks in the goaf are fully developed.
(a)
(b)
(c)
(d)
3.3.4. Determination of the Height Range of the Two Zones of the Overlying Strata
The stress discrimination method is used to divide the height range of the caving zone, the fracture zone, and the bending subsidence zone formed by the movement and deformation of the overlying rock during the mining process. The relationship between the principal stress in the middle position of the alignment and the inclination directions and the distance of the roof of the goaf when the working face in the middle position of the overburden strata is mined at different distances is shown in Figure 11 and 12. As the advancing distance of the working face increases, the heights of both the caving zone and the fracture zone show a slight increase trend. In the alignment direction, the height of the caving zone ranges from 6.2 to 8.7 m, and the height of the fracture zone ranges from 33.5 to 38.6 m. In the inclination direction, the height of the caving is in the range of 6.5 to 9.2 m, and the height of the fracture zone is in the range of 34.8 to 40.5 m.
(a)
(b)
(c)
(d)
(a)
(b)
(c)
(d)
4. On-Site Measurement and Analysis
In order to further study the influence of the three-side goaf area in the mining process of the 120502 working face, analyze the behavior of the appearance of rock pressure in the overlapping and noncoinciding sections with the upper goaf area, and verify the accuracy of the numerical simulation stress and strain results, on-site monitoring of the leading support pressure and the displacement of the roadway roof at 120502 island working face was conducted.
4.1. Measuring Point Layout and Monitoring Method
The length of the 120502 working face is 400 m. After fully considering the impact of the open-off cut and the stop line to protect the coal pillars, three measurement points for both stress monitoring and displacement monitoring of the working face are, respectively, set up in the transportation roadway and the return air roadway at both ends of the working face. They are located at 100 m, 200 m, and 300 m away from the working face, respectively. The borehole stress gauge was used to measure the advance supporting pressure of the working face, and the roof separation gauge was used to measure the displacement of the surrounding rock of the overlying roof of the roadway.
4.2. Analysis of Monitoring Results of Two Roadways
Figure 13 shows the change curve of the leading support pressure at different advancing distances during the field measurement of the working face. The peaks of the leading support pressure of the working face in the coincident and noncoincident sections are located 10 m in front of the coal wall. In the overlap section, affected by the pressure relief of the upper goaf, the peak value of the leading support pressure is between 17.5 and 19.2 MPa, with an average of 18.2 MPa. The peak range of the leading support pressure in the noncoincident section is between 30.2 and 37.1, with an average of 33.3 MPa, which is larger than the coincident zone and consistent with the numerical simulation.
The measured vertical displacement curve of the roof of the roadway is shown in Figure 14. In the overlap section, the maximum displacement of the top plate is about 210 mm. The roof displacement of the noncoincident section is obviously larger than that of the coincident section. The displacement is about 300 mm, and it is not stable. It has a tendency to continue to increase as the working face advances, which is consistent with the simulation result.
5. Conclusions
(1)During the excavation of isolated island working face, the stress distribution presents the following characteristics. In the inclination direction, the pressure relief area expands with an arch upward movement as the mining distance increases and merges into a “saddle shape”. In the alignment direction, in the section coincident with the upper goaf, the stress concentration area is behind the open-off cut and in front of the working face. The vertical stress first decreases and then increases with the increase of the mining distance and, finally, presents a “horizontal three-zone” distribution. The support pressure in front of the working face is increasing, and the distance between the peak support pressure and the coal wall is stable at about 10 m. In the noncoincident zone, the stress concentration area remains unchanged, and the peak stress concentration is also equivalent. However, under the same advancing distance, the peak support pressure in front of the work is higher.(2)The roof displacement has the following characteristics. In the alignment direction, in the section that overlaps with the upper goaf, the vertical displacement direction of the overlying rock strata in the goaf changes from downward to upward, and the peak displacement remains stable. In the noncoincidence section, the vertical displacement of the overlying strata in the goaf is downward, and the peak value is increasing. The horizontal displacement increases continuously with the advancement of the working face in the alignment direction. The overlying rock is subjected to the combined action of horizontal stretching and compression, and thus, the cracks develop more fully. The failure range of the plastic zone is a “saddle shape” with a nonuniform expansion trend and an upward dynamic development. The plastic zone on both sides of the coal pillars has a large development height, while the roof plastic zone has a small development range. With the increase of the advancing distance of the isolated island working face, the height of the caving zone and the fracture zone both show a slight increase trend.(3)The field measurement shows that the leading support pressure and top plate displacement of the noncoincident section obviously exceed that of the coincident section, and the peaks are all 10 m in front of the coal wall of the working face. Similarly, the vertical displacement of the roof in the noncoincident section is also greater than that in the coincident section
Data Availability
The data used to support the findings of this study are available from the corresponding author upon request.
Conflicts of Interest
The author declares that there are no conflicts of interest regarding the publication of this paper.
Acknowledgments
This work was supported by the National Natural Science Foundation of China (51404010 and 51374012).