Abstract

Deep mining has become the normal state of coal mining; compared with the mine with shallow buried depth, the consequent high level of in situ stress and complex distribution have brought severe threats to the stability of the stope and the surrounding rock of the roadway. In this research, taking the 121304 working face of Kouzidong Mine as the engineering background, the characteristics of mining-induced stress distribution under complex in situ stress environment in deep mining are analyzed by using on-site measurement of the original rock stress and mining stress, establishing a theoretical model centered on the middle section of the working face, and establishing large-scale numerical calculation models for different advancing directions. It was found that under deep mining conditions, the maximum stress of the original rock is 25.12 MPa, and the direction is vertical. The advanced influence range of mining stress is about 150 m, and the abutment pressure presents a three-peak distribution characteristic in front of the working face. The research results provide important theoretical guiding value for guiding the mining of coal mines with similar geological conditions.

1. Introduction

As the level of modern industry increases year by year, the eastern part of China needs more energy supplies [16]. The shallow resource was exhausted, and the deep ones begin to develop. As a result, the level of ground stress has risen and the distribution law has become more complicated [711]. The stability maintenance of surrounding rock during mining and roadway excavation is more complicated [1216]. The critical factor which causes the fracture of the rock strata of the stope and roadway is in situ stress [1720].

Radwan et al. [21] have determined the ground stress field, pore pressure, stress regime, and fracture gradient through on-site measurement of ground stress. It is discovered that the stress field rotates affected by geologic structure and the influence of the stress field rotation on the effect of resource extraction, thereby providing a reliable guarantee for the efficient exploitation of resources. Saeidi et al. [22] analyzed the in situ stress distribution characteristics of Niobec Mine, Quebec (Canada) and the influence of in situ stress uncertainty on stability evaluation of open-pit slope. Whether a mining plan can be exploited successfully to a deeper mining level depends largely on whether the magnitude and direction of the in situ stress state can be fully considered. Many mines use in situ stress to evaluate the stability of structure of underground and, at the same time, the interpretation of various aspects of mine data and engineering judgments are generally based on mine depth functions and methods which are determined by the value of ground stress obtained from the results of stress measure. Guo et al. [23] statistically analyzed the distribution features and rules of ground stress in Chinese mines according to a massive data of measured results of ground stress and summarized the relationship between stress field type, stress magnitude, stress value, and buried depth. In order to study the splitting failure of high sidewall cavern under high in situ stress, Zhang et al. [24] carried out true three-dimensional (3D) geomechanical model test. In Zhang’s paper, the main powerhouse of Pubudou Hydropower Station is taken as an example. By using self-developed 3D loading system and analogue materials, Zhang reproduced the splitting failure of high sidewall cavern under high ground stress. The formation factor and collapse feature of the split failure of the high side rock cavern have been revealed in this model test. Furthermore, it summarized credible evidence for further researches on the formation mechanism of the split failure of the high sidewall cavern under the action of high ground stress. Yang et al. [25] put forward stress initialization methods. According to the results from simulations, the stress initialization results obtained by implicit-explicit sequence method and dynain file method had a close relationship to the material model of rock. Compared with other methods, the explicit DR method has an obvious advantage in solution time. According to the on-site measurement of a deep mine, Qin et al. [26] have found that the immediate cause of the roadway deformation is the ground stress. A kind of support form of high strength prestressed bolt and cable anchor is proposed, and the best prestress ratio of the bolt and cable anchor was proposed according to the corresponding geological conditions. Li et al. [27] proposed an intelligent monitoring and dynamic control method for coal gangue image recognition to provide the optimal illumination environment for coal gangue recognition.

In summary, the magnitude and direction of in situ stress play an important role in the control of surrounding rock. Especially in deep mines, in situ stress has a greater impact on mining work [2833]. When the working face is in the condition of large buried depth, due to the influence of high ground stress, the mine pressure behavior will be more severe under the mining disturbance, and the creep deformation of the roadway will occur for a long time after the support [3436]. The increase of working face length will lead to the increase of stress concentration factor at both ends of the working face, and the increase of the length will slow down the advancing speed, which makes it difficult to support the stope and roadway at the end of the working face and cause serious coal wall spalling. This paper analyzes the influence of the angle between the advancing direction and the initial minimum principal stress on the mining stress field. A better understanding of the influence of the angle on the mining stress field mentioned above can better strengthen the roadway in a targeted way and achieve safe and efficient production. The law of influence on the mining stress field is of great significance for guiding the study of the stability of the surrounding rock in deep mines.

2. Engineering Geological Background

The 121304 working face of Kouzidong Mine, AnHui, is taken as the engineering background. Due to the large buried depth of the working face, and the length of the working face is more than 350 m, it is a typical working face with complex ground stress and significant mine pressure. Rib spalling, roadway severe deformation, and floor heave are serious, as shown in Figure 1.

13-1 coal has been being mined at the working face, in which thickness is 5.2 m, and the comprehensive geological column chart is shown in Figure 2. The coal seam was mined by the process of mining all heights at one time. The on-site stress measurement of the working face found that the maximum principal stress is the vertical stress with a magnitude of 25.12 MPa. The intermediate principal stress is a horizontal stress with a magnitude of 21.84 MPa, and the direction is 35° west by north, and the minimum principal stress is the horizontal stress with a magnitude of 12.89 MPa, in which direction is perpendicular to the intermediate principal stress.

3. Distribution Characteristics of Coal Abutment Pressure

A simplified mechanical model of the working face section is established to analyze the distribution of abutment pressure in the coal seam area. As shown in Figure 3, the working face is regarded as a roadway with a circular cross section, surrounded by an equal-pressure homogeneous surrounding rock environment. According to the content and assumption of elastic mechanics, this is regarded as the stress distribution around the hole. In order to simplify the calculation, the weight of the surrounding rock within the range of 3–5r0 around the hole (working face) can be ignored in the calculation process. 4r0 is used in this paper.

Assuming the radius of the working face section is r0, the polar coordinate system is established with the center of the section as the polar coordinate origin. The coordinates of any point outside the working face are (r, θ); by analysing the stress and strain of the point, the equilibrium differential equation of the point is listed as the following equation:where σr is the radial stress, σt is the tangential stress, and τ represents the shear stress; in this model τθr = τ = 0, fr is the physical strength.

In order to simplify the calculation process, according to the characteristics of the engineering environment, equation (1) is simplified and further equation (2) is obtained,

For the convenience of calculation, the equation is approximated,

Analyze the differential volume at this point and get the geometric equation, as shown in the following equation:where εr is the radial strain, εt is the tangential strain, and u is the radial deformation of differential volume.

According to the generalized Hooke’s law, the strain expression of differential volume is obtained in the following equation,where E is the elastic modulus, μ is Poisson’s ratio, and σz is the axial stress of circular hole.

Taking the derivative of equation (4), one has

Substituting equation (5) into equation (6) leads to

From equation (6),

Combining equations (7) and (8),

Coordinate transformation in the rectangular coordinate system σr, σt means, as shown in the following equation,

According to Lame’s solution, let the internal pressure q1 = 0 and q2 = γH in the thick wall cylinder model.

By calculating the limit equilibrium area around the working face, the static equation is obtained as equation (3), according to the limit equilibrium condition.where C is the cohesion and φ is the internal friction angle.

Combining equations (3) and (12) and integral,where A is the integral constant.

Substituting the boundary conditions r = r1 and σr = 0,when the roof caving behind the working face fills the goaf, the abutment pressure value in this area gradually recovers from 0 to the original rock stress,where K is the stress recovery factor.

According to the comprehensive analysis for the abutment pressure at the coal seam position, from the front of the working face to the rear, it can be divided into elastic zone (a), limit equilibrium zone (b), complete pressure relief zone (c), and stress recovery zone (d), as shown in Figure 4.

4. Numerical Analysis of Characteristics of Coal Mining Stress Distribution

4.1. Model Establishment and Parameter Determination

Large numerical calculation models are established based on the engineering geological background of Kouzidong Mine, as shown in Figure 5. In order to study the impact of the advancing direction on the distribution of mining stress, six groups of models with the angle between the advancing direction and the minimum principal stress direction of 0°, 15°, 35°, 45°, 60°, and 90° were established, respectively. According to the results of on-site stress measurement, X axis is determined as the minimum principal stress direction, Y axis is the middle main stress direction, Z axis is the maximum main stress direction, and the top of the model is the stress boundary condition, where the vertical stress load is applied to simulate the weight of overburden which is not established, and the surrounding and bottom are fixed displacement boundary conditions. The length of each model working face is 350 m and the width of roadway is 6 m. In order to facilitate modeling and improve the accuracy of solution, the length, width, and height of each unit in the coal seam, roof, and floor located in the excavation range of the working face are 1 m, and the thickness of the coal seam is 6 m. After the working face is excavated, the excavated unit body is filled to simulate the goaf. Step by step excavation mode is adopted for the model with the same advancing direction as the actual site, which is calculated once for each excavation of 5 m to simulate the daily advance distance, and the goaf is filled with 20 m behind the coal wall.

In the process of simulation, the self-developed constitutive model is used to control the mechanical characteristics of coal and rock mass. In order to determine the model parameters, coal and rock samples are collected at 121304 working face to determine the rock mechanical parameters. According to the measured distribution characteristics of coal and rock fractures, combined with Hoek Brown criterion, the rock parameters of different layers are modified, and the rock mass parameters are obtained, as shown in Table 1. After excavation, the double yield model is used to simulate the compaction bearing characteristics of gangue in goaf. The model parameters are shown in Table 2. The parameters of the double yield model are determined by trial and error method: the model parameters are changed constantly to make the simulated stress-strain curve of gangue in goaf consistent with the prediction result of Salamon model.

4.2. Numerical Model Validation

In order to determine the reliability of the model, four borehole stress gauges were installed 200 m away from the front working face of the return airway, which will measure the distribution characteristics of air side leading abutment pressure of 121304 working face. The monitoring results are shown in Figure 6(a). When it is 150 m away from the working face, abutment pressure begins to be affected by mining, showing an upward trend, reaching a peak 10–12 m before the working face, and the coal body reached the limit equilibrium state; after that, the coal body entered the mining damage zone, the bearing capacity decreased, and the abutment pressure began to decrease. The coal wall near the working face is reduced to the minimum value, which is about the residual strength of the coal body. The abovementioned abutment pressure distribution characteristics indicate that the leading mining influence range of 121304 working face reaches 150 m, and the width of the coal body damage zone in front of working face reaches 10–12 m.

The on-site measurement and numerical calculation results are shown in Figure 6. Comparing the two results, it can be found that the trends of the two results are consistent on one side of the working face. In the length direction of the working face, there is a trend of first rising and then falling back from the end to the middle. Peaks appear on two sides of the working face. The measured results show that the peak value is close to 50 MPa. The numerical calculation results show that the peak position is 7 m away from the roadway. The abutment pressure reaches 57.8 MPa. Along the advancing direction, the change trend of the abutment pressure showed a first increase and then a decrease and finally approached the original rock stress. The peak points of numerical simulation and measured results are located 10–12 m in front of the working face, and the influence range was about 150 m. It can be seen that the numerical calculation model is basically consistent with the on-site measured results, and the calculation parameters are reliable.

4.3. Analysis of Numerical Results

In order to observe the variation characteristics of mining stress under different advancing directions, two measuring lines and a plane are arranged in the coal seam 5 m above the coal seam floor. The measuring plane is parallel to the horizontal plane. The position of measuring line is showed in Figure 5. As shown by the red line, one measuring line is arranged in the middle of the working face along the advancing direction, and the other measuring line is arranged at 5 m in front of the working face along the length direction, as shown by the pink line and the light blue line.

The influence of the advancing direction on the distribution characteristics of the maximum principal stress is shown in Figure 7. Under the influence of coal seam mining, the maximum principal stress of the working face is obviously concentrated in advance of the working face, and the maximum principal stress concentration area is mainly distributed in front of the working face, and at both ends, the leading maximum principal stress presents a three-peak shape on the middle of the working face and both sides.

In the length direction of working face, the three peaks at different included angles change significantly, as shown in Figure 8(a). When the angle between the advancing direction and the minimum principal stress is 0°, the maximum principal stress peaks on both sides of the working face are greater than the peak in the middle of the working face, and the peaks on the left and right sides are 51 MPa. As the angle between the advancing direction of the working face and the minimum principal stress increases, the leading maximum principal stress of the working face first decreases and then increases. When the angle is equal to 45°, the maximum principal stress is the smallest. At this time, the peak at the middle position is 49 MPa and the peak at both sides is 47 MPa. When the included angle increases to 90°, the peak at the middle position is 53 MPa, and the peak at both sides is 51 MPa. In the distribution form of the maximum principal stress of the tendency, with the increase of the included angle, the measuring line in the direction of the tendency gradually began to show an asymmetric distribution, and the peak value of the maximum principal stress in the middle part was gradually greater than the peaks at both sides.

In the advancing direction of working face, under different included angles, the maximum principal stress distribution of the working face is shown in Figure 8(b). The angle between the advancing direction and the minimum principal stress has a less impact to the maximum principal stress of the working face, and the leading influence range is about 40 m. When the included angle is 0°, the maximum principal stress peak position of the working face is 6 m in front of the working face, and the peak point is 51.9 MPa. As the included angle increases, the peak point position moves slightly to the front of the coal wall, and the maximum principal stress peak increases slightly. When the included angle increases to 90°, the maximum principal stress peak position of the working face is 4 m before the working face, and the peak point is 52.5 MPa. In general, the angle between the advancing direction and the minimum principal stress has little effect on the distribution of the maximum principal stress.

Figure 9 shows the maximum principal stress concentration factor of the working face under different angles between the advancing direction and the minimum principal stress. The stress concentration factor is expressed in η, where B is the peak stress of working face and A is the in situ stress without mining disturbance. As shown in Figure 9, the maximum stress concentration factor of principal stress appears to decrease first and then increase. The minimum value of the maximum principal stress concentration factor is at the angle of 15° with a value of 2.16. As the included angle continues to increase, the stress concentration factor gradually increases. When the included angle increases to 45°, the increase began to slow down. When the included angle increases to 90°, the stress concentration factor reaches the maximum value of 2.25.

In summary, the angle between the minimum principal stress and the advancing direction of the working face has a greater influence on the mining stress distribution in working face. In terms of stress, with the increase of the angle, the mining stress as a whole shows a trend of decreasing first and then increasing. In terms of morphological distribution, the spatial division characteristics of mining stress are also affected by the included angle. With the increase of the included angle, the peak values of the double peaks on both sides are gradually smaller than the peaks in the middle.

Under the condition of different included angles between the direction of the initial minimum ground stress (X-axis) and the advancing direction of the working face, the evolution characteristics of the maximum principal stress rotation trajectory are shown in Figure 10. The red, green, and blue measuring lines are arranged along the advancing direction of the working face. The red measuring line is arranged in the middle of the working face, and the green and blue measuring lines are arranged 20 m away from the roadway at both ends. The pink and light blue measuring lines are arranged from both ends of the tailgate and headgate to the middle of working face along the length direction of the working face, 5 m ahead of the working face. The maximum principal stress rotation trajectory on the advancing direction measuring line always deflects in the direction of the goaf in the vertical plane parallel to the advancing direction of the working surface. Therefore, the angle between the 90–270 axis direction and the maximum principal stress rotation trajectory on the strike measuring line gradually increases. It always maintains the angle between the advancing direction of the working face and the initial minimum ground stress direction. The maximum principal stress rotation trajectory on the length direction measuring line always maintains the arc of the symmetry axis with the advancing measuring line, roughly toward the roadway on both sides of the working face direction. In summary, the size of the rotation angle of the maximum principal stress is less affected by the direction of propulsion, but the direction of rotation is obviously affected by the advancing direction.

5. Conclusions

(1)The on-site measurement has obtained the in situ stress distribution of the deep part of the Kouzidong Mine. The maximum principal stress direction is the vertical direction with a magnitude of 25.12 MPa. The direction of the minimum principal stress is east by north 55° of the horizontal direction, and the magnitude is 12.89 MPa. The deep mine of Kouzidong Mine has an advanced mining influence range of 150 m, and the maximum stress concentration position is about 10 m in front of the coal wall. The results of numerical simulation agree well with this feature.(2)Under the influence of coal seam mining, the maximum principal stress of the working face is obviously concentrated in advance of the working face, and the maximum principal stress concentration area is mainly distributed in front of the working face, and at both ends, the leading maximum principal stress presents a three-peak shape on both sides of the working face and the middle of the working face.(3)Through the establishment of numerical models of different advancing directions, it is found that the stress concentration is located at 7–10 m in front of the working face, and in the length direction of the working face, the distribution characteristics of three-peak stress are all presented. With the increase of the angle between the advancing direction and the minimum principal stress, the mining stress concentration first decreases and then increases and tends to be stable. Therefore, in the safe and efficient production of super-long working face in deep mine, it is necessary to consider the influence between the advancing direction and the minimum principal stress angle, so as to minimize the stress concentration of the working face and strengthen the support at the end and the middle of the working face to prevent the occurrence of rib spalling and other accidents.

Data Availability

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

Conflicts of Interest

The authors declare that they have no conflicts of interest.

Acknowledgments

This work was financially supported by the National Natural Science Foundation of China (51934008), which is gratefully acknowledged.