Next Article in Journal
Effect of Municipal Solid Waste Compost on Yield, Plant Growth and Nutrient Elements in Strawberry Cultivation
Next Article in Special Issue
Dynamic Tensile Failure Characteristics and Energy Dissipation of Red Sandstone under Dry–Wet Cycles
Previous Article in Journal
A Sustainable Port-Hinterland Container Transport System: The Simulation-Based Scenarios for CO2 Emission Reduction
Previous Article in Special Issue
Prediction and Application of the Height of Water-Conducting Fracture Zone in the Composite Roof: A Case Study of Jinxinda Coal Mine
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Study on the Distribution Trend of Rockburst and Ground Stress in the Hegang Mining Area

1
School of Engineering for Safety and Emergency Management, Taiyuan University of Science and Technology, Taiyuan 030024, China
2
Intelligent Monitoring and Control of Coal Mine Dust Key Laboratory of Shanxi Province, Taiyuan University of Science and Technology, Taiyuan 030024, China
3
School of Civil Engineering, Hunan University of Science and Technology, Xiangtan 411201, China
*
Author to whom correspondence should be addressed.
Sustainability 2023, 15(12), 9445; https://doi.org/10.3390/su15129445
Submission received: 8 May 2023 / Revised: 1 June 2023 / Accepted: 7 June 2023 / Published: 12 June 2023
(This article belongs to the Special Issue Advances in Coal Mine Disasters Prevention)

Abstract

:
As the depth of coal mining increases, rockbursts become more severe, and multiple rockburst accidents have previously occurred, seriously affecting coal mine production safety. With the Hegang mining area as the research object and using field measurement and numerical simulation as the research methods, the geological and tectonic characteristics of each impact mine in the Hegang mining area are analyzed to obtain the tectonic stress field distribution characteristics of each mine and the tectonic stress field distribution trend in the Hegang mining area. The maximum horizontal principal stress and the minimum horizontal principal stress on the energy field in the overrun area of the retrieval working face is analyzed via numerical simulation, and the influence trend of the ground stress on the impact ground pressure is finally obtained. Results show that using the typical working face where rockburst occurs as the geological proto-type, the influence of the direction of the ground stress field, the maximum horizontal principal stress and the minimum horizontal principal stress on the energy field in the overrun area of the retrieval working face is analyzed via numerical simulation, and the influence trend of the ground stress on the impact ground pressure is finally obtained. When the angle is 70°~90°, an energy peak in the overrun area of the working face reaches its maximum. The ratio of the minimum horizontal stress to the vertical stress is positively correlated with the energy concentration in the overrun area of the working face, but its change has a minimal effect on the energy distribution in the overrun area of the working face. When this ratio is increased from 0.6 to 1.4, the peak energy of the simulated working face only increases by 8.22%, and the energy concentration area remains basically unchanged.

1. Introduction

With the increase in mining depth, the impact pressure situation in the Hegang mining area is becoming severe, with many rockburst accidents occurring that seriously affect the safe production of coal mines. The high ground stresses at deep mining levels and the high stress fields created by the superimposed stresses of mining are two of the root causes of rockburst hazards [1,2]. Rockburst is a typical dynamic disaster in coal mining [3]. Usually, when the mechanical system of coal and rock reaches its ultimate strength, it releases its elastic properties in a sudden, rapid and violent manner, causing instantaneous damage to the coal and rock layers; this is accompanied by the impact of coal powder and rock, resulting in damage to the shaft and personal injury accidents [4,5,6]. When rockburst occurs, it often causes great damage to the support of the roadway, leading to large-scale roof fall, wall spalling and a floor heave in the roadway, causing great damage to personnel and equipment and seriously affecting coal mine production safety. Moreover, the occurrence of rockburst is accompanied by the ejection of coal and rock, causing disasters such as coal and gas outbursts, coal and gas dust explosions, fires, floods and even ground vibration, which can cause the collapse of ground buildings and pose great harm [7,8]. Stress is the most fundamental cause of rockburst, and the stress of a coal–rock mass includes crustal stress and disturbance stress; thus, monitoring and controlling the stress in a coal–rock mass is very important for the research and prevention of rockburst.
The safe, efficient and green mining of deep coal seams are some of the major concerns with ensuring the healthy development of a country’s modern energy economy and an important guarantee of environmental sustainability. Deep coal seams not only have the essential property of high static load characteristics, but also have the additional property of strong disturbances resulting from large-scale mining activities, leading to the frequent occurrence of dynamic hazards, such as impact ground pressure, which can easily cause significant casualties and economic losses. Impact pressure control can overburden deformation, reduce surface subsidence and prevent impact pressure disasters, and is an effective way to achieve safe, efficient, sustainable and green mining in deep coal seams. Therefore, this study has far-reaching significance for the prevention and control of impact pressure and sustainable development.
Much research has been carried out by domestic and international scholars and field engineers on the relationship between ground stress and rockburst. Lianpeng Dai [9] focused on the geological and mining technical factors that influence the occurrence of impact ground pressure in a roadway. The impact ground pressure area of the roadway was quantitatively divided into several different risk level areas, and an impact ground pressure roadway zoning and grading management system was proposed. Masoud Ghorbani [10] investigated the behaviour of rockburst in high stress rock masses, described different types of energy absorbing bolts and other support components, and determined that a suitable support system was the main condition for preventing rockburst. Jian Zhou [11] used the stochastic gradient boosting (SGB) method to study impact ground pressure data, analyzed the effect of different indicators, such as the ground stress condition factor on the impact ground pressure and predicted the impact ground pressure damage dataset with accurate and reliable results. Charlie C. Li [12] analyzed a combination of seismic monitoring and numerical modeling as a promising approach to forecasting. Using existing knowledge and technology to predict impact ground pressure and large seismic events was extremely challenging, but prediction was possible. Jian Zhou [13] used a hybrid technique combining an artificial neural network (ANN) and an artificial bee colony (ABC) to establish a complex relationship between susceptibility to ground pressure and its influencing factors, which could be used as an effective tool for predicting susceptibility to ground pressure. Yigai Xiao [14] proposed a metal mine pressure monitoring and warning method based on deep learning data analysis. The theoretical basis of rockburst was elaborated, the inducing factors of deep well rockburst were analyzed and the classification of rockburst was examined. Hongpu Kang [15] conducted an on-site experiment involving large-scale hydraulic fracturing to weaken the hard and refractory rock layer above the long wall panels of underground coal mines. The strength of the long wall working face was significantly reduced, and the strong dynamic load pressure generated by the large-scale fracture of the roof was basically eliminated. Shengquan He [16] used the simultaneous compressional transform (SST) method and kinetic modeling to analyze the mechanism of remotely triggered rockbursts around coal mine roadways. Based on seismic monitoring data, the seismic waveform at the source was estimated using the empirical scaling method and calibrated in the model.
At present, with the further increase of mining depth and the further deterioration of the geological environment, the original research results are not suitable for deep underground engineering. The numerical simulation analysis methods of a geostress field are mainly applied in hydraulic and hydropower engineering, and the study area is small and the geological structure is simple. Reasonable boundary conditions can be obtained by using a linear regression analysis method, a nonlinear mathematical method, an optimization function and other methods to carry out numerical simulation of a geostress field in the study area. However, for the inversion of a geostress field in complex geological conditions and at large range, the existing numerical simulation analysis method is not suitable. With the advancement of technology, there have been great advances in the development of ground stress measurement methods and inversion analysis methods for ground stress fields [17]. Moreover, with the development of [18], the engineering problems encountered have become increasingly complex and serious, causing many issues for construction safety. With the increase in mining intensity and depth, rockburst is now a frequent occurrence in mining activities and a serious impediment to the safe development of coal mine production. As a fundamental source of force for the occurrence of impact ground pressure disasters in mines, ground stress is highly valued by scholars at home and abroad and has been thoroughly studied with effective results. However, with the further increase in mining depth, the geological environment further deteriorates, and the original research results are not applicable to deep underground studies.
The Hegang mine is located in the city of Hegang in the northeastern part of Heilongjiang Province at the southern foot of the Xiaoxinganling Mountains. The coal strata are 800 to 1200 m thick and contain 36 coal seams, of which thicker seams greater than 3.5 m account for 75.5%. The mines in this region are mostly exploited by vertical-slope mixed multilevel mining, with coal hardness coefficients ranging from 1.4 to 3.0 in mostly high gas mines; several mines in the Hegang mine area have already experienced impact ground pressure. At present, the mining depth of the mine has extended to 500~1100 m, and the situation of impact ground pressure is more severe. The high ground stress at depth and the high stress field formed by the superimposed stress of mining are the root causes of the impact ground pressure disasters and the large deformation damage of the roadway. The main objectives of our study are to test the ground stress in typical impact mines in the Hegang mining area to understand the distribution pattern of the deep ground stress field in the mines, to analyze the influence of ground stress on impact ground pressure and to evaluate the regional risk of impact ground pressure hazards in the Hegang mining area.

2. Geological Overview

2.1. Geological Formations

The mine is located in Hegang City, Heilongjiang Province. It is 100 km long from north to south and 28 km wide from east to west, with a total area of approximately 2800 km2. The average surface elevation of the mine area is 290 m, with elevations ranging from 250 to 340 m. The whole area has a topographic pattern of high elevation in the southwest and a low elevation in the northeast. A geological structure map of the Hegang mining area is shown in Figure 1.

2.2. Ground Stress Conditions

In situ crustal stress measurements are the most effective and direct means to obtain the state of crustal stress in a region. There are many measurement methods. The most widely used measurement methods are hydraulic fracturing [19,20] and stress relief [21,22]. Small-aperture hydrofracture measurement techniques and hollow envelope stress relief methods have been more widely used in ground stress testing at coal mine sites. The 56 mm diameter of the hydraulic fracture test borehole greatly reduces the need for deformation in the borehole, reduces the size of the test rig and increases the speed of testing. The measurement of ground stress is a precise measurement task, and factors such as geological formations, tunnel distribution characteristics, mining activities, construction space and the rock deposit characteristics of the test location can all have an impact on the measurement of ground stress. In addition, coal mine production is a complex production and operation system, and the operation of any link affects the operation of the entire production system.
The geostress measurements were mainly carried out for seven impact mines in the Hegang mining area. According to the principles of determining the ground stress measurement points, the mine area was arranged through the analysis of the geological conditions on site and combined with the actual production situation on site. The locations of the measurement points are shown in Figure 2 and Table 1.
The arrangement of the ground stress measurement points was determined in each typical impact mine in the Hegang mining area. The feasibility of the construction of the ground stress measurement points was examined by testing the plastic zone and the structure of the surrounding rock, and the depth of the ground stress test borehole was also determined. The strain data for each measurement point were then obtained according to the ground stress test procedure, and the optimal position for the ground stress measurement was obtained by screening the raw data and results. The results of the ground stress at each measurement point in the Hegang mine are shown in Table 2 and summarized in Figure 3.
The maximum horizontal principal stress at each measurement point in the Hegang mining area has a clear tendency to increase with depth, with the maximum horizontal principal stress ranging from 17.9 to 22.8 MPa at burial depths of 400 m to 500 m. The maximum horizontal principal stress ranges from 19 to 32.5 MPa at depths of 700 m to 800 m. The maximum horizontal principal stress ranges from 35.5 to 41.2 MPa at depths of 800 m to 900 m. The maximum horizontal principal stress ranges from 39.6 to 39.8 MPa at depths of 900 m to 1000 m. Although the maximum horizontal principal stress increases as the depth of burial increases, the difference between the maximum and minimum values of the maximum horizontal principal stress at each burial level is 10–20 MPa. There are evident regional characteristics that are influenced by the regional structure as shown in Figure 4.

3. Characteristics of Rockburst in the Hegang Mining Area

3.1. Characteristics of Rockburst during Driving

(1) Range of rockburst pressure: In the Hegang mining area, most rockburst occurs within the range of 15 m~110 m from the heading face, while most rockburst occurs within 50 m from the heading face, accounting for 86.7% of the total [23,24].
(2) Release energy: Use microseismic monitoring measurements of the energy released by rockbursts during roadway excavation in the Hegang mining area is calculated on the Richter scale, ranging from 1.25 to 2.86, with an average of 1.86. The energy ranges from 2.0 × 105 to 1.07 × 108 J, and the average energy is 2.15 × 106 J. By analyzing the impact energy and the impact energy rate during the mining process, solid coal mining and goaf mining, the impact energy curve as shown in Figure 5 is obtained.
(3) Failure characteristics: The failure characteristics of rockburst in the Hegang mining area are mainly manifested in four aspects: roadway deformation, supporting body failure, equipment damage and casualties. When it is subjected to artificial mining and other geological activities, it has a severe impact and destruction phenomenon, which can easily cause a coal mine rockburst accident, resulting in fully mechanized coal mine underground working face or roadway, resulting in sheet wall, floor bulge, roof collapse, support failure, roadway blockage, equipment damage, casualties and other destructive consequences.
Roadway deformation: The floor bulging ranges from 0.3 to 1.8 m, and the two sides move closer to the range from 0.1 to 1.7 m. The minimum section of the roadway after deformation is 1.3 m2. Roadway deformation is the most serious in the range from 15 m to 50 m outwards from the heading head. There is a patch phenomenon in the head, and the depth of the patch is 0.5 m.
Support body damage: The rigging equipment and anchor cable rope can be damaged, and the U-shaped steel can become deformed.
Equipment damage: The transport machine and boring machine can be sprung up or tilted; the rig can be overturned, plus many others.
Casualties: Partial rockburst can cause casualties.

3.2. Characteristics of Rockburst of the Working Face

(1) Impact range: The impact ground pressure in the Hegang mine area is most serious in the Junde and Nanshan coal mines. The location where the impact ground pressure causes serious damage occurs initially in the return air duct, followed by the working face, with relatively little impact damage in the machine duct. The impact damage in the return duct is generally distributed within 110 m from the exit position on the working face to the overhead working face, and the serious damage is generally within 20 m to 60 m of the overhead working face. A few shocks occur in the range of 100 m to 255 m of the overburden. The number of impacts on the working face is low, and they occur in the range of 20 to 100 groups of hydraulic supports, and the severe sections are mostly distributed in the range of 50 to 90 groups of hydraulic supports. The impact pressure in the machine path occurs less frequently, and the impact range is from the lower exit position of the working face to 110 m of the overrun working face.
(2) Release energy: The energy released from the impact ground pressure occurring during the workings retrieval process in the Hegang mining area ranges from 0.81 to 4 on the Richter scale, with an average of 2.23. The energy ranges from 3.57 × 104 J to 9.12 × 109 J, and the average energy is 9.12 × 106 J.
(3) Damage characteristics: The damage caused by the impact ground pressure that occurs during the workings retrieval process in the Hegang mining area can be studied in terms of both roadway damage and workings damage.
Roadway damage: Roadway damage is most serious in the return airway, the bottom drum range is 0.5~1.3 m, the top sink is 0.5~1.0 m and the two helpers shrink 0.6~2.0 m. Roadway deformation is serious, the section sharply shrinks, a part of the impact ground pressure causes the return airway to severely bubble and the number of impact ground pressure that causes bubbles accounts for 18.5% of the total statistics; when impact ground pressure occurs, it is accompanied by an upper helper, a lower helper and a broken bottom plate. The rocks are thrown out, equipment and pipelines near the upper and lower helpers are also ejected and overturned; the overhead support body is severely damaged, the bottom and top beams are bent and cracked, the single hydraulic pillar fails and tilts, the U-beam is bent and deformed, and the anti-impact bracket is damaged. Moreover, a part of the impact ground pressure leads to casualties.
Damage to the working face: the coal mining machine tilts towards the coal gang, the working face shelf shed tilts towards the soft gang, and some of the hydraulic brackets are damaged, accompanied by bottom drum, flaky gang and other phenomena.

4. Influence of Ground Stress on Rockburst

4.1. Model Establishment

The working section in the third level north of the 17 layers of the third four sections in the Junde coal mine is used as a simulation prototype; this area has a strike length of 920 m, an inclination length of 168 m, a mining height of 3.5 m and a mining depth of 451–551 m. The coal seam is 9.63 to 15.92 m thick and consists of bright and dark coals with an inclination of 25° to 35°, with an average of 30°. The direct top is fine sandstone with a thickness of 4–7 m, the old top is medium and fine sandstone with a thickness of 30–40 m, and the direct base is tuffaceous siltstone with a thickness of 4.5–5.5 m. The old bottom is coarse sandstone with a thickness of 48–50 m. The western windway side of the working face is the upper section of the mining void, the coal column of the section is 4–47 m wide and the eastern machine side is solid coal; in this section, the back of the northern working face is a stratified mining void, and the southern stop line is near the L1 and F7 large faults. The coal seams adjacent to this face are #3, #9, #11, #21 and #17, which are not fully mined, leaving a large number of irregular coal pillars. Due to the influence of various factors, such as ground stress, coal rock seam impact tendency, stage coal pillars, interval coal pillars left in the upper section, adjacent mining area and high intensity production, this working face has experienced several impacts of ground pressure.
The size of the model is 500 m × 300 m × 500 m, the width of the working face is 150 m, the dip angle of the rock seam is 30 degrees, the thickness of the coal seam is 4 m, the thickness of the direct top is 6 m, the thickness of the old top is 35 m and the thickness of the direct bottom is 5 m. This simulation adopts long-wall mining and roof caving management. According to the initial state close to the stratum, the initial stress caused by the dead weight of the overlying rock is simulated. The full thickness method was used for a simulation analysis, and the stress balance was calculated at each advance step before the next advance. The model is divided into 81,252 units and 93,340 nodes. According to the measured stress value at measurement point #3 of the Junde coal mine, the load value of the upper boundary was determined, and the horizontal load of the model was also determined according to the angle between the back mining face and the direction of ground stress. The bottom of the model is surface constrained. The model uses the Moore–Coulomb criterion (Figure 6), and the physical and mechanical parameters of each rock formation are shown in Table 3.
The superposition of high ground stresses and mining stresses is the root cause of many impact ground pressures, and the effect of the direction of ground stress on impact ground pressures during workface recovery is investigated. The influence of the angle between the direction of the maximum horizontal principal stress and the direction of the advance of the working face on impact ground pressure is analyzed, and the variation in energy in the area of the overrun working face is studied when the angles are 0, 20, 45, 70 and 90 degrees.

4.2. Results Analysis

According to the calculation results of each boundary load at different pinch angles, the load is applied to the model and calculated; the fish program is written to extract the energy value of the overrun area of the working face, and the energy value of the overrun area of the working face during the retrieval process is plotted as a graph. The calculated results and graphs are shown in Figure 7, Figure 8, Figure 9, Figure 10 and Figure 11.
As shown in Figure 7, Figure 8, Figure 9, Figure 10 and Figure 11, the area of energy concentration in the working face is biased towards the lower chute (machine path) due to the influence of the inclined rock formation. When the initial pressure is applied, the energy in the leading area of the working face is the highest. During the mining process, the energy of the solid coal seam on the lower side of the goaf along the dip direction is greater than that of the solid coal seam on the upper side. As the face advances, the energy in the solid coal seam below the mining area and in the area ahead of the face gradually increases and shifts away from the face, while the energy in the solid coal seam below the mining area is gradually greater than the energy in the area ahead of the face. This feature explains the concentration of severe impact pressure sections in the windrow (i.e., the upper chute) during face retrieval. Therefore, after the upper section is mined, the coal pillar is reserved and the mining of the next section begins; specifically, the mining of the next section is influenced by the empty area and the coal pillar of the upper section, and the mining activities of the working face and the return windway (i.e., the upper chute) are also carried out in the area of its influence and can definitely become an important cause of the occurrence of impact ground pressure.
As shown in Table 4, energy values extracted from the overrun area of the working face are as follows: (1) The peak energy of the coal seam at the initial incoming pressure is 2.7 × 105 J when the maximum principal stress direction is at an angle of 0 degrees to the working face advance direction. At an angle of 20 degrees, the peak initial incoming pressure energy is 3.3 × 105 J. At an angle of 45 degrees, the peak initial incoming pressure energy is 5.9 × 105 J. At an angle of 20 degrees, the peak initial incoming pressure energy is 3.3 × 105 J. At an angle of 45 degrees, the peak initial incoming pressure energy is 5.9 × 105 J. At an angle of 70 degrees, the peak initial incoming pressure energy is 9.5 × 105 J. At an angle of 90 degrees, the peak initial incoming pressure energy is 9.3 × 105 J. The angle between the direction of the maximum horizontal principal stress and the direction of the workface advance during the recovery process has a greater influence on the energy distribution in the overfront area of the workface; a larger angle correlates to more energy gathered in the overfront area of the workface. When the angle is greater than 70 degrees, the peak energy in the overfront area of the workface reaches its maximum. The gathered energy reaches the conditions for the occurrence of impact ground pressure, which also releases the most energy and causes the most damage. (2) When the angle between the direction of the maximum principal stress and the working face advance direction is 0 degrees and 20 degrees, the peak energy gathering position at the initial incoming pressure is 10 m from the working face. At 45 degrees, the peak energy gathering position is 15 m from the working face; at 70 degrees, the peak energy gathering position is 25 m from the working face; at 90 degrees, the peak energy gathering position is 20 m from the working face. As the angle between the direction of maximum horizontal principal stress and the working face advance direction increases, the energy in the overrun area of the working face is gradually transferred to the interior.

4.3. Horizontal Principal Stress and Weight Stress

The ratio of the maximum horizontal principal stress to the self-weight stress is K1, and the ratio of the minimum horizontal principal stress to the self-weight stress is K2. At K1 = 1.8, K2 has very little effect on the energy distribution in the overrun area of the working face, its peak energy changes minimally as K2 increases and the energy concentration area is relatively constant. As K2 increases, the peak energy in the overrun area of the working face increases, but to a lesser extent. Compared to the peak energy at K2 = 0.6, the peak energy increases by 1.59% at K2 = 0.8, 4.26% at K2 = 1.0, 6.29% at K2 = 1.2 and 8.22% at K2 = 1.4. The energy value of the coal seam in the overrun area of the working face is U1, the energy value of the coal seam in the stable area of the working face is U2, and the ratio of U1 to U2 is f. An area with f > 1.2 is considered the energy concentration area, and the distribution range of the energy concentration area under different K2 conditions is calculated as shown in Table 5. From the table, K2 has a minimal influence on the distribution range of the energy concentration area when K1 = 1.8 (see Figure 12).
As the value of K1 increases, the energy in the overrun area of the face and in the coal seam area below the mining area gradually increases, and the energy concentration area also increases. Figure 13 shows the graph of the energy distribution in the overrun area of the working face for different K1 conditions at K2 = 0.8. From the graph, the effect of the change in K1 value on the energy distribution in the overrun area of the working face has the following characteristics: (1) The peak energy is affected in the overrun area of the working face; with the increase in K1, the peak energy of the overrun area of the working face gradually increases. When K1 changes from 1.8 to 2.0, the greatest change in peak energy in the overrun area of the working face is observed. When K1 = 2.2, the peak energy in the overrun area of its working face doubles relative to the peak energy at K1 = 1.4; moreover, the change in K1 has a certain influence on the location where the peak energy gathers, and as K1 increases, the peak energy tends to shift towards the interior of the coal seam. When K1 = 2.2, the location where the peak energy appears in the overburden area of the working face shifts from 10 m to 15 m. (2) The energy concentration zone is also affected; the change in the K1 value has a greater impact on the distribution of the energy concentration zone of the coal seam. As shown in Table 6, with an increasing K1 value, the energy concentration area of the coal seam in front of the working face gradually increases.

4.4. Stress Control Technology

Controlling the occurrence of impact ground pressure hazards essentially avoids high stress concentrations due to the superposition of the two by controlling ground stress or mining stress to ensure that sudden impact destabilization damage does not occur in the coal–rock body.
(1) Mine ground stress measurements: Using a technically mature and highly accurate ground stress test method, the magnitude and direction of the ground stress are measured at three or more different levels in the mine by selecting measurement points.
(2) Three-dimensional geomechanical modeling: Combined with the engineering geological data of the mine, a three-dimensional geomechanical model including geological anomalies, such as faults and folds and different rock formations in the study area, is established, and the influence of these geological anomalies is fully considered when analyzing the ground stress field.
(3) Regional prevention of impact pressure: From the start of coal mining design, the requirements of ground stress are used, reasonable coal mining methods and techniques are selected, the mining sequence of coal seams is optimized and the coal seam is selected that is most conducive to impact pressure prevention and control as the first to be mined as a protective layer. After determining the order of coal seam mining, the orientation of the roadway layout for the identified coal seams is optimized such that the direction of the roadway axis is predominantly consistent with the direction of the maximum horizontal principal stress and the level of ground stress around the roadway is minimized, thereby decreasing the degree of stress concentration that can induce impact ground pressure and reducing the impact risk at the working face. For the identified working face, numerical simulation methods are combined to simulate the dynamic mining process in the coal mine, analyze the stress distribution state, evaluate the impact hazard at different locations of the working face and carry out impact ground pressure monitoring and prevention for areas with different risk levels.

5. Conclusions

(1) A numerical model was established using a typical impact ground pressure working face in the Hegang mining area as the geological proto-type; the influence of the direction of the ground stress field, the maximum horizontal principal stress and the minimum horizontal principal stress on the impact ground pressure was analyzed. As the angle between the maximum horizontal principal stress direction and the working face advance direction increases, the energy peak in the overtopping area of the working face shifts to the interior of the coal seam, the distribution range of the energy concentration area gradually increases, and the energy peak in the overtopping area of the working face reaches the maximum when the angle is 70°~90°.
(2) The ratio of the maximum horizontal principal stress to the vertical stress, K1, has a greater influence on the distribution of energy in the coal seam in the overrun area of the working face. As K1 increases, the peak energy in the overrun area of the working face increases, and the energy concentration area also increases. The ratio of minimum horizontal stress to vertical stress, K2, is positively correlated with the energy concentration in the overburden area of the working face; however, its change has less influence on the energy distribution in the overburden area of the working face. When K2 increases from 0.6 to 1.4, the peak energy of the simulated working face only increases by 8.22%, and the energy concentration area remains basically unchanged.
(3) The impact pressure area stress control technology includes four steps: ground stress measurement, 3D geomechanical modeling, 3D ground stress field inversion and impact pressure area prevention. The concept of the impact pressure area stress control technology is to cause the design and layout of the mining face to conform to the requirements of the ground stress to reduce the base stress level of the working face and to avoid the impact pressure occurrence caused by the superposition of the ground stress and the mining stress to form a high stress concentration.

Author Contributions

Conceptualization, J.X.; Methodology, Y.H. (Yongliang He); Software, Q.H.; Validation, Y.H. (Yongjiang Hao); Writing—original draft, J.P. All authors have read and agreed to the published version of the manuscript.

Funding

This research was funded by the National Natural Science Foundation of China (42202321), by the Scientific and Technological Innovation Programs of Higher Education Institutions in Shanxi (2019L0641, 2021L330), by the Fundamental Research Program of Shanxi Province (202203021222184), and by the Taiyuan University of Science and Technology Scientific Research Initial Funding (20222112).

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Data Availability Statement

The data used to support the findings of this study are available from the corresponding author upon request.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Li, Z.L.; He, X.Q.; Dou, L.M.; Song, D.Z.; Wang, G.F. Numerical investigation of load shedding and rockburst reduction effects of top-coal caving mining in thick coal seams. Int. J. Rock Mech. Min. Sci. 2018, 110, 266–278. [Google Scholar] [CrossRef]
  2. Wen, Z.; Wang, X.; Tan, Y.; Zhang, H.; Huang, W.; Li, Q. A Study of Rockburst Hazard Evaluation Method in Coal Mine. Shock Vib. 2016, 2016, 8740868. [Google Scholar] [CrossRef] [Green Version]
  3. Wu, Y.; Gao, F.; Chen, J.; He, J. Experimental Study on the Performance of Rock Bolts in Coal Burst-Prone Mines. Rock Mech. Rock Eng. 2019, 52, 3959–3970. [Google Scholar] [CrossRef]
  4. Wu, Y.; Gao, F.; Chen, J.; He, J. Mining-induced static and dynamic loading rate effect on rock damage and acoustic emission characteristic under uniaxial compression. Saf. Sci. 2019, 116, 86–96. [Google Scholar]
  5. Mottahedi, A.; Ataei, M. Fuzzy fault tree analysis for coal burst occurrence probability in underground coal mining. Tunn. Undergr. Space Technol. 2019, 83, 165–174. [Google Scholar] [CrossRef]
  6. Yang, Y.; Wei, S.; Li, K. Inverse analysis of dynamic failure characteristics of roadway surrounding rock under rockburst. Energy Sci. Eng. 2021, 9, 2298–2310. [Google Scholar] [CrossRef]
  7. Liu, D.; Ling, K.; Guo, C.; He, P.; He, M.; Sun, J.; Yan, X. Experimental simulation study of rockburst characteristics of Sichuan–Tibet granite: A case study of the Zheduoshan tunnel. Eng. Geol. 2022, 305, 106701. [Google Scholar] [CrossRef]
  8. Li, P.-X.; Feng, X.-T.; Feng, G.-L.; Xiao, Y.-X.; Chen, B.-R. Rockburst and microseismic characteristics around lithological interfaces under different excavation directions in deep tunnels. Eng. Geol. 2019, 260, 105209. [Google Scholar] [CrossRef]
  9. Dai, L.; Pan, Y.; Zhang, C.; Wang, A.; Canbulat, I.; Shi, T.; Wei, C.; Cai, R.; Liu, F.; Gao, X. New Criterion of Critical Mining Stress Index for Risk Evaluation of Roadway Rockburst. Rock Mech. Rock Eng. 2022, 55, 4783–4799. [Google Scholar] [CrossRef]
  10. Ghorbani, M.; Shahriar, K.; Sharifzadeh, M.; Masoudi, R. A critical review on the developments of rock support systems in high stress ground conditions. Int. J. Min. Sci. Technol. 2020, 30, 555–572. [Google Scholar] [CrossRef]
  11. Zhou, J.; Shi, X.-Z.; Huang, R.-D.; Qiu, X.-Y.; Chen, C. Feasibility of stochastic gradient boosting approach for predicting rockburst damage in burst-prone mines. Trans. Nonferrous Met. Soc. China 2016, 26, 1938–1945. [Google Scholar] [CrossRef]
  12. Li, C.C.; Mikula, P.; Simser, B.; Hebblewhite, B.; Joughin, W.; Feng, X.; Xu, N. Discussions on rockburst and dynamic ground support in deep mines. J. Rock Mech. Geotech. Eng. 2019, 11, 1110–1118. [Google Scholar] [CrossRef]
  13. Zhou, J.; Koopialipoor, M.; Li, E.; Armaghani, D.J. Prediction of rockburst risk in underground projects developing a neuro-bee intelligent system. Bull. Eng. Geol. Environ. 2020, 79, 4265–4279. [Google Scholar] [CrossRef]
  14. Xiao, Y.; Deng, H.; Xie, Z.; Lu, H. Mine Ground Pressure Monitoring and Early Warning Based on Deep Learning Data Analysis. Mob. Inf. Syst. 2022, 2022, 6255119. [Google Scholar] [CrossRef]
  15. Kang, H.; Jiang, P.; Feng, Y.; Gao, F.; Zhang, Z.; Liu, X. Application of Large-Scale Hydraulic Fracturing for Reducing Mining-Induced Stress and Microseismic Events: A Comprehensive Case Study. Rock Mech. Rock Eng. 2023, 56, 1399–1413. [Google Scholar] [CrossRef]
  16. He, S.; Chen, T.; Vennes, I.; He, X.; Song, D.; Chen, J.; Mitri, H. Dynamic Modelling of Seismic Wave Propagation due to a Remote Seismic Source: A Case Study. Rock Mech. Rock Eng. 2020, 53, 5177–5201. [Google Scholar] [CrossRef]
  17. Eremenko, A.A.; Mashukov, I.V.; Eremenko, V.A. Geodynamic and Seismic Events under Rockburst-Hazardous Block Caving in Gornaya Shoria. J. Min. Sci. 2017, 53, 65–70. [Google Scholar] [CrossRef]
  18. Li, Z.; Chen, Q.; Xue, Y.; Qiu, D.; Chen, H.; Kong, F.; Liu, Q. Numerical investigation of processes, features, and control of land subsidence caused by groundwater extraction and coal mining: A case study from eastern China. Environ. Earth Sci. 2023, 82, 82–93. [Google Scholar] [CrossRef]
  19. Cheng, Q.; Huang, B.; Shao, L.; Zhao, X.; Chen, S.; Li, H.; Wang, C. Combination of Pre-Pulse and Constant Pumping Rate Hydraulic Fracturing for Weakening Hard Coal and Rock Mass. Energies 2020, 13, 5534. [Google Scholar] [CrossRef]
  20. Huang, B.; Cheng, Q.; Zhao, X.; Kang, C. Hydraulic fracturing of hard top coal and roof for controlling gas during the initial mining stages in longwall top coal caving: A case study. J. Geophys. Eng. 2018, 15, 2492–2506. [Google Scholar] [CrossRef] [Green Version]
  21. Xue, Y.; Gao, F.; Gao, Y.; Cheng, H.; Liu, Y.; Hou, P.; Teng, T. Quantitative evaluation of stress-relief and permeability-increasing effects of overlying coal seams for coal mine methane drainage in Wulan coal mine. J. Nat. Gas Sci. Eng. 2016, 32, 122–137. [Google Scholar] [CrossRef]
  22. Liu, Z.; Cheng, Y.; Liu, Q.; Jiang, J.; Li, W.; Zhang, K. Numerical assessment of CMM drainage in the remote unloaded coal body: Insights of geostress-relief gas migration and coal permeability. J. Nat. Gas Sci. Eng. 2017, 45, 487–501. [Google Scholar] [CrossRef]
  23. Wang, Y.; Li, X.; Wu, Y.F.; Ben, Y.X.; Li, S.D.; He, J.M.; Zhang, B. Research on relationship between crack initiation stress level and brittleness indices for brittle rocks. Chin. J. Rock Mech. Eng. 2014, 33, 265–275. [Google Scholar]
  24. Zhang, Z.Z.; Gao, F. Confining pressure effect on rock energy. Chin. J. Rock Mech. Eng. 2015, 34, 1–11. [Google Scholar]
Figure 1. Geological structure map of the Hegang mining area.
Figure 1. Geological structure map of the Hegang mining area.
Sustainability 15 09445 g001
Figure 2. The location of in situ stress measurement points in the Hegang mining area.
Figure 2. The location of in situ stress measurement points in the Hegang mining area.
Sustainability 15 09445 g002
Figure 3. Results of the in situ stress measurements in the Hegang mining area.
Figure 3. Results of the in situ stress measurements in the Hegang mining area.
Sustainability 15 09445 g003
Figure 4. Relationship of the in situ stress measurements in the Hegang mining area.
Figure 4. Relationship of the in situ stress measurements in the Hegang mining area.
Sustainability 15 09445 g004
Figure 5. Rockburst energy. (a) Impact energy curve; (b) Impact energy ratio.
Figure 5. Rockburst energy. (a) Impact energy curve; (b) Impact energy ratio.
Sustainability 15 09445 g005
Figure 6. Numerical model.
Figure 6. Numerical model.
Sustainability 15 09445 g006
Figure 7. Energy nephograms of the coal seam when the angle is 0°.
Figure 7. Energy nephograms of the coal seam when the angle is 0°.
Sustainability 15 09445 g007
Figure 8. Energy nephograms of the coal seam when the angle is 20°.
Figure 8. Energy nephograms of the coal seam when the angle is 20°.
Sustainability 15 09445 g008
Figure 9. Energy nephograms of the coal seam when the angle is 45°.
Figure 9. Energy nephograms of the coal seam when the angle is 45°.
Sustainability 15 09445 g009
Figure 10. Energy nephograms of the coal seam when the angle is 70°.
Figure 10. Energy nephograms of the coal seam when the angle is 70°.
Sustainability 15 09445 g010
Figure 11. Energy nephograms of the coal seam when the angle is 90°.
Figure 11. Energy nephograms of the coal seam when the angle is 90°.
Sustainability 15 09445 g011
Figure 12. Energy distribution of the area in front of the work face under different K2 conditions when K1 is 1.8.
Figure 12. Energy distribution of the area in front of the work face under different K2 conditions when K1 is 1.8.
Sustainability 15 09445 g012
Figure 13. Energy distribution of the area in front of the work face under different K1 conditions when K2 is 0.8.
Figure 13. Energy distribution of the area in front of the work face under different K1 conditions when K2 is 0.8.
Sustainability 15 09445 g013
Table 1. Technical characteristics of the in situ measurement points in the Hegang mining area.
Table 1. Technical characteristics of the in situ measurement points in the Hegang mining area.
Measuring PointsBurial (m)LocationBorehole
Depth (m)Azimuth (°)Dip Angle (°)
Junde Coal Mine1#940Belt stone gate7.31813
2#470Return air stone gate7.41445
3#627Track stone gate7.1264
Xing ’an coal mine1#730Switchboard road7.21293
2#7542 #Switchboard Road6.93504
3#563.5Boundary stone gate7.41883
Fuli coal mine1#720−450 roadway7.52654
2#800−530 roadway7.21703
3#880−610 roadway7.31883
Xinlu coal mine1#940−650 Stone gates7.21733
2#990−700 water silos7.51435
3#830−540 Return air lane7.1824
Nanshan coal mine1#539−120 return roadway7.310304
2#539−120 return roadway6.9305
3#508−180 roadway6.81365
Yixin coal mine1#487Floor runway7.21213
2#635Oblique well7.3944
3#561Air inlet lane7.1794
Xingshan coal mine1#459Belt lane72254
2#706−340 bottom bend71754
3#458−90 rock lane754
Table 2. Results of the in situ stress measurement in the Hegang mining area.
Table 2. Results of the in situ stress measurement in the Hegang mining area.
Measuring PointsLocationBurial (m)Principal StressVertical Stress
Principal Stress (MPa)Azimuth (°)Dip Angle (°)
Junde Coal Mine1#Belt stone gate720σ133.4287−8.221.16
σ218.73177−23
σ310.81267−52
2#Return air stone gate470σ122.87798.713.2
σ210.87237−13
σ38.29152−62
3#Track stone gate627σ132.5111.5−2.517.49
σ216.42−62−14.4
σ314.41205−56.4
Xing ‘an coal mine1#Switchboard road730σ130.10113.0712.1223.61
σ216.8920.939.88
σ314.97252.7574.26
2#2 #Switchboard Road754σ132.7269.164.1624.03
σ216.85167.2362.63
σ314.75259.3526.99
3#Boundary stone gate563.5σ130.4886.61.7213.95
σ214.47193.684.14
σ314.09103.7−5.60
Fuli coal mine1#−450 roadway720σ135.990.91.621.67
σ222.318119.8
σ310.6−4.373.1
2#−530 roadway800σ139.2762.5918.69
σ221.226132.6
σ310.117057.2
3#−610 roadway880σ141.295.44.322.37
σ225.7184.3−5.9
σ318.978.864.2
Xinlu coal mine1#−650 Stone gates940σ139.64264.8910.2426.14
σ221.27−32.2778.43
σ318.81175.86−5.34
2#−700 water silos990σ139.79116.667.5727.93
σ224.0925.687.35
σ319.95252.0079.42
3#−700 water silos830σ135.4990.444.6120.98
σ217.57−62.5084.83.
σ315.40180.622.34
Nanshan coal mine1#−120 return roadway539σ127.878132.817−8.13814.223
σ223.383−80.934−78.472
σ319.972228.4267.949
2#−120 return roadway539σ125.729133.117−8.63814.438
σ223.267−80.93478.072
σ320.027228.4367.949
3#−180 roadway508σ123.892124.11618.00212.946
σ220.672226.34240.549
σ318.837177.91720.877
Yixin coal mine1#Floor runway487σ121.7107.37.7413.6
σ29.825.513.3
σ38.7236.774.5
2#Oblique well635σ119.089.30.2313.4
σ210.5179.424.6
σ39.9269.965.3
3#Floor runway561σ120.370.7−0.2212.0
σ211.5157.785.7
σ38.9247.74.22
Xingshan coal mine1#Belt lane459σ117.9116.53.0111.7
σ29.4−33.279.8
σ39.3219.39.65
2#−340 bottom bend708σ121.6107.211.0115.76
σ211.3−39.317.8
σ38.6162.872.0
3#−90 rock lane458σ118.6−84.3−1411.10
σ216.5184−6.6
σ313.2249.774
Table 3. Physical and mechanical parameters of the rock masses for the numerical model.
Table 3. Physical and mechanical parameters of the rock masses for the numerical model.
Rock MassDensity/kg·m−3Bulk Modulus/GPaShear Modulus/GPaCohesion/MPaFriction Angle/(°)Tensile Strength/MPa
Overlying strata25604.22.95.0341.5
Medium-fine sandstone27213.472.085.237.62.81
fine sandstone25582.011.452.441.982.16
Coal14200.460.190.8200.01
siltstone26305.03.86.0352.5
Coarse sandstone25604.22.95.0341.5
goaf20100.460.190.8200.01
Table 4. Energy distribution situation of the area in front of the work face at different angles.
Table 4. Energy distribution situation of the area in front of the work face at different angles.
AnglePeak Energy/JDistance from Working Face/mEnergy Concentration Areas/m
02.7 × 1051010~50
203.3 × 1051010~55
455.9 × 1051510~80
709.5 × 1052515~95
909.3 × 1052015~70
Table 5. Energy distribution situation of the area in front of the work face in different K2 conditions when K1 is 1.8.
Table 5. Energy distribution situation of the area in front of the work face in different K2 conditions when K1 is 1.8.
K2Peak Energy/JIncrease RateDistance from the Working Face/mEnergy Concentration Areas/m
0.6324,81001010~40
0.8329,9601.59%1010~40
1.0338,6604.26%1010~40
1.2345,2506.29%1010~40
1.4351,5008.22%1010~40
Table 6. Energy distribution situation of the area in front of the work face in different K1 conditions when K2 is 0.8.
Table 6. Energy distribution situation of the area in front of the work face in different K1 conditions when K2 is 0.8.
K1Peak Energy/JIncrease RateDistance from the Working Face/mEnergy Concentration Areas/m
1.4239,95001010~40
1.6270,17012.59%1010~45
1.8329,96037.51%1010~45
2.0442,58084.45%1010~50
2.2480,020100.05%1510~55
Disclaimer/Publisher’s Note: The statements, opinions and data contained in all publications are solely those of the individual author(s) and contributor(s) and not of MDPI and/or the editor(s). MDPI and/or the editor(s) disclaim responsibility for any injury to people or property resulting from any ideas, methods, instructions or products referred to in the content.

Share and Cite

MDPI and ACS Style

Pang, J.; Xie, J.; He, Y.; Han, Q.; Hao, Y. Study on the Distribution Trend of Rockburst and Ground Stress in the Hegang Mining Area. Sustainability 2023, 15, 9445. https://doi.org/10.3390/su15129445

AMA Style

Pang J, Xie J, He Y, Han Q, Hao Y. Study on the Distribution Trend of Rockburst and Ground Stress in the Hegang Mining Area. Sustainability. 2023; 15(12):9445. https://doi.org/10.3390/su15129445

Chicago/Turabian Style

Pang, Jiewen, Jianlin Xie, Yongliang He, Qiaoyun Han, and Yongjiang Hao. 2023. "Study on the Distribution Trend of Rockburst and Ground Stress in the Hegang Mining Area" Sustainability 15, no. 12: 9445. https://doi.org/10.3390/su15129445

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop