Reader TTS

ABSTRACT Xinzhuang Coal Mine in Huaneng Qingyang Mining Area has geological characteristics such as large buried depth, thick topsoil and hard overburden, which is a typical kilometer-deep

rockburst mine.Currently, the mine remains in the fundamental construction phase. This provides a theoretical foundation for the mine’s subsequent production processes. , this paper studies

the stress distribution law of surrounding rock in the stope of Xinzhuang Coal Mine. First of all, select the Panel 1 west wing of Coal 8 as a typical representative, and establish a large

three-dimensional refined numerical model. Then, through numerical simulation research, the characteristics of surrounding rock abutment stress distribution zone in the mining process of the

working face under complex conditions are obtained, which provides support for the regional design of rockburst prevention and control. The main research conclusions are as follows: (1) The

research reveals the distribution characteristics of mining stress, elastic energy accumulation, and plastic zone at different widths of the first mining face, and determines the reasonable

width of the working face. It is considered that the width of the working face can be greater than 175 m and less than 250 m, and the design width 200 m of the working face of Xinzhuang

Coal Mine is reasonable. (2) The research reveals the evolution law and influence factors of mining stress, elastic energy accumulation and plastic zone distribution in the first mining

face. The research results can not only provide theoretical support for the subsequent production of Xinzhuang Coal Mine, but also provide reference for the rockburst prevention and control

work of mines with similar geological conditions. SIMILAR CONTENT BEING VIEWED BY OTHERS RESEARCH ON STRESS RELEASE AND PRESSURE RELIEF MECHANISM OF UNDERLYING COAL AND ROCK UNDER PROTECTIVE

LAYER MINING Article Open access 15 February 2025 STUDY ON FISSURE EVOLUTION OF OVERLYING ROCK IN LOWER PROTECTIVE MINING Article Open access 19 August 2024 EVOLUTION OF SURROUNDING ROCK

DISTORTION ENERGY DURING CLOSELY SPACED COAL SEAM MINING AND INNOVATIVE COAL PILLAR DESIGN METHOD Article Open access 20 May 2025 INTRODUCTION The safe and efficient mining of coal resources

is the eternal theme of the entire coal industry1.With the depletion of shallow resources, China’s mining depth extends to the deep at the rate of 10–15 m per year. According to statistics,

72% of mines buried more than 800 m face the threat of rockburst. Under the condition of thick topsoil and hard overburden, the stress field of surrounding rock shows strong nonlinear

characteristics after mining disturbance, and the traditional theory is difficult to accurately predict the mechanism of rockburst. In the mining process of mining coal mine, mining

activities will inevitably lead to the overburden movement and the redistribution of surrounding rock stress, so the prevention and control of mine pressure disasters is particularly

important2. Among all kinds of mine pressure disasters, rock burst disasters are extremely prompt and destructive, and have different mechanisms under different geological conditions, so it

is very difficult to efficiently and accurately monitor, warn and prevent3. In recent years, many scholars have conducted in-depth studies on the mechanical mechanism of rock burst

generation, providing important references for the early warning and prevention of such disasters4,5. Typical research achievements are as follows: Professor Dou Linming proposed the

mechanism of dynamic and static load superimposed to induce scourging, and explained the interrelationship between high static load condition of stope surrounding rock, dynamic load

generation caused by overburden migration, and rockburst disaster6. Researcher Pan Junfeng conducted an in-depth study on the scour prevention mechanism of artificial liberated layer, and

demonstrated the feasibility and effectiveness of active regional scour prevention measures7. Professor Wang Hongwei studied the abrupt effect of coupling instability of fault overlying rock

and identified the mechanical relationship between fault slip instability and rock burst8. Professor Song Dazhao studied the mechanism of impact in near-upright coal seam mining, and

proposed a static dynamic evaluation method for near-upright coal seam impact hazard. Huaneng Qingyang Mining area is currently planned to have two pairs of mines, Walnut Yu coal mine and

Xinzhuang Coal mine, both of which are kilometers deep and have typical geological characteristics such as thick topsoil and hard overlying rock. Walnut Yu coal mine has been put into

operation in October 2020, during the mining period, the mine earthquake appeared obviously, and there were many large energy vibration events. Xinzhuang Coal Mine, which is adjacent to

Walnut Yu Coal Mine, is expected to face more serious threat of rock burst disaster after putting into operation. In view of the fact that Xinzhuang Coal mine has not yet carried out

relevant research on the mechanism of rock burst disaster, in order to provide a theoretical basis for the subsequent coal mine operation, this paper studies the stress behavior law of the

surrounding rock of Xinzhuang Coal mine9. MINE ENGINEERING OVERVIEW GENERAL SITUATION OF QINGYANG MINING AREA Qingyang mining area has two coal mines, Walnut Yu Coal mine and Xinzhuang Coal

mine, with large geological reserves, high coal quality and good market prospects. It is an important support for Longdong energy base and occupies an important position in the coal industry

of Huaneng Group Company. However, the coal seams of the two mines are generally buried deep, and there are multiple disasters of Luohe Formation roof water damage, rock burst (mine

earthquake), oil-type gas, coal seam spontaneous ignition and other disasters, which are mutually induced, so it is extremely difficult to govern the disaster. Binchang Mining area adjacent

to similar conditions has now been defined by the State Administration of Mine Safety as “one of the mining areas with the heaviest coal disaster, the most difficult to govern, the highest

safety risk, and the most heavy task of effectively preventing and containing major and major accidents”10,11. Qingyang Mining area has complex geological conditions such as large buried

depth of coal seam, high stress of raw rock, complex hydrological conditions and the occurrence of thick hard roof. Its two coal mines are rock burst mines, among which the Huoyu coal mine

started construction in 2009 and put into trial production in 2020, and the Xinzhuang Coal mine started construction in 2013. At present, it is still in the construction period, and it is

planned to start production in 2024. The mine earthquake phenomena occurred in both mines during the period of roadway excavation and initial operation. In view of the above situation,

Qingyang mining area has planned a number of rock burst research projects in recent years, and carried out a certain amount of research work, but the overall research foundation is still

relatively weak. OVERVIEW OF XINZHUANG COAL MINE Xinzhuang Coal Mine is located about 15 km south of Ningxian County, Qingyang City, Gansu Province, the administrative division of Xinzhuang

Town and Zhongcun Township, Ningxian County, geographical coordinates: 107°43′59″–107°59′55″ east longitude; North latitude 35°14′41″–35°22′19″, east–west length of about 20.0 km,

north–south width of 7.6–12.5 km, an area of 206.2823 km2, the design capacity of 8Mt/a. The Xinzhuang minefield contains coal layers 5–1, 5–2 and 8, which are three recoverable coal layers

(see Table 1). Among them, the coal 5–1 layer and the coal 5–2 layer belong to relatively stable and mineable coal seam, which is the minor mineable coal seam in the minefield. The 8 layer

of coal belongs to relatively stable and recoverable coal seam, which is the main recoverable coal seam in the minefield.The spacing between coal layer 5–1 and coal layer 5–2 is 0.81–15.13 

m, and the spacing between coal layer 5–2 and coal layer 8 is 1.66–58.41 m. Coal layer 8 is the coal seam of the main mining seam and the coal seam to which the first mining face belongs and

the coal seam and the bottom floor have weak impact tendency. The roof of coal seam is mainly mudstone, followed by silty mudstone or argillaceous siltstone, and occasionally siltstone and

coarse sandstone roof. The coal seam floor is dominated by siltstone, sandy mudstone and carbonaceous mudstone, and partially by mudstone and argillaceous siltstone, and occasionally coarse

sandstone floor. The coal layer 8 is to adopt the full-mechanized mining method of caving top coal. The design mining height is 3.5 m, the width of working face is 200–250 m, and the

advancing length is 2500–5000 m. NUMERICAL MODEL ESTABLISHMENT AND SIMULATION SCHEME The first mining area of Xinzhuang Coal mine is a pan area, and the main recoverable coal seam is 8

layers of coal. Therefore,this simulation takes the mining of coal face in the 8-pan area as the main research object, in order to obtain the abutment stress, elastic energy and plastic zone

distribution characteristics of surrounding rock during the mining process under complex conditions such as large buried depth, thick topsoil and hard overlying rock, and provide support

for regional anti-impact design. MODELING SCOPE According to the supplementary exploration data of the previous mining area, the coal seams in Coal 8 Plate 1 area are relatively stable, and

the thickness of the coal seams is gradually thinning at the southern stripping boundary. Not affected by the stubble of the overlying coal seam, in order to put the mine into production as

soon as possible and reach production, the design of the first mining face of Coal 8 Plate Area is arranged in the western part of the plate area, and according to the mining and digging

succession plan of Coal 8 Plate Area, the mining range of the working face of the western wing of Coal 8 Plate Area is selected as the main modeling area. Geodetic coordinates X coordinate

range: 36484904.4203–36491404.4203 m, Y coordinate range: 3906759.8339–3910259.8339 m, and Z coordinate range: − 400 m to the surface of the ground, to carry out the large-size

three-dimensional refinement modeling, the model size of 6500 m × 3500 m × 1680 m, the modeling area is shown in Fig. 1. The complex conditions such as the change of surface soil thickness,

the change of buried depth of coal seam, the change of coal seam thickness, the change of coal seam inclination, the change of hard roof thickness, and the syncline structure are included in

the established numerical model, which is typical of Xinzhuang Coal mine. CONSTRUCTION OF NUMERICAL MODEL According to the connection plan of the stoping face in Coal8 No.1 block of

Xinzhuang Coal Mine, the layout of the working face in the west wing of Coal8 No.1 block of Xinzhuang Coal Mine and the layout of the central return air main roadway, the central auxiliary

transportation main roadway and the central belt conveyor main roadway, the layout of the working face in Coal8 No.1 block is shown in Fig. 2. The width of the first mining face is preset to

be 150–300 m (the design width of the mine is 200 m), and the width of the coal pillar in the section between the first mining face and the adjacent successive working face 1 is preset to

be 5 m, and the width of the coal pillar between the successive working face 2 and 3 is preset to be 5 m. The distance between the stop-mining line of each working face and the nearest

central return air main lane is 150 m, the section size of the central main lane is 6 m × 6 m, and the interval between the main lanes is 50 m. The 3D numerical simulation was established

based on the borehole exploration information and the continuous layout of the working face. The model is 6500 m long and 3500 m long. The bottom elevation of the model is − 400 m, the top

elevation of the model is + 1280 m and the minimum is + 880 m. There are 2,315,212 units with 1,021,362 nodes. In the model, 1802 face (X = 1350–5450 m, Y = 2165–2465 m), 1803 face, 1805

face and 1807 face are preset, and 5 m wide coal pillars are set between each adjacent face. Five central alleys are arranged on the east side of the preset stoping line of each face. From

west to east are the central return air main lane, the central auxiliary transport main lane, the central belt conveyor main lane. The 3D numerical model is shown in Fig. 3. According to the

ground stress conditions, normal displacement is fixed on both sides of the model in X and Y directions and horizontal stress is applied. Fixed normal displacement at the bottom of the

model; The top of the model is the surface, and the free boundary is set. G is set to 10 m/s2. The Moore-Coulomb criterion model is adopted for numerical calculation. The physical and

mechanical parameters of coal and rock mass in this study are shown in Table 212,13. NUMERICAL SIMULATION SCHEME EXCAVATION SCHEME FOR NUMERICAL SIMULATION OF WORKING FACE WIDTH Face width

plays a critical role in the mining process. The maximum height of the formation separation rupture is closely correlated with the size of the gob, extending from the cutting hole to the

full mining stage. Therefore, the width of the working face has an important influence on the appearance of rock burst disaster in the mine14,15,16,17,18. In order to study the influence of

different stoping widths on the adjacent working face to be worked, the distribution characteristics and peak position of lateral supporting stress, elastic strain energy and peak value of

surrounding rock, and the distribution of impact risk of surrounding rock, this paper takes the 1802 working face of Coal Mine No. 8 in Yipan District of Xinzhuang Coal Mine as an example.

When 1802 working face width is 150, 175, 200, 225, 250, 275, 300 m, the stress and energy field, plastic zone of surrounding rock are studied, as shown in Fig. 4. The excavation scheme of

the face width simulation study is as follows: * (1) Based on the stress balance results of the original rock, the central return air main lane, central auxiliary transportation main lane

and central belt conveyor main lane on the east side of the excavation model are calculated for balance; * (2) Based on the balance calculation results after the excavation of the central

main lane, the 1802 working face (4100 m long) is excavated in two times. The first excavation is 2000 m and the calculation is balanced; Then the 1802 working face was excavated 2100 m to

calculate the balance. The excavation diagram of the working face is shown in Fig. 5. * (3) The excavation scheme simulated by the “three fields” evolution law of surrounding rock during the

mining process of the working face. * (1) Numerical simulation excavation scheme of the “square” process of the first mining face In order to analyze and study the distribution and

evolution characteristics of stress and elastic strain energy of surrounding rock at the boundary of the 1802 face, the distribution law of rock burst risk of surrounding rock before and

after “square” of the 1802 face is studied. Taking the first mining face of 1802 as the research object, the 1802 working face is excavated to 100, 150, 200 (square), 250 and 300 m in turn

with the excavation step of 50 m. After each excavation, the balance is solved, and then the excavation is carried out according to the 400 m step, as shown in Figs. 6 and 7. * (2) Numerical

simulation excavation scheme for the mining process of the first mining face In order to study the distribution and evolution characteristics of supporting stress, peak stress, elastic

strain energy and peak elastic strain energy of surrounding rock during the mining process of working face, analyze the influence of complex geological conditions such as surface elevation,

yellow soil thickness, coal seam depth, hard roof strata, coal seam thickness and dip Angle changes, and delineate the concentrated areas of high stress and high elastic energy. It provides

theoretical basis for rock burst prevention and control in actual mining of working face. The following excavation scheme is designed. After the 1802 face is “squared”, the overall stoping

face 1802 with a excavation step of 400 m is adopted, and the equilibrium is solved after each excavation, as shown in Fig. 8. THE NUMERICAL SIMULATION ANALYSIS STUDY OF REASONABLE WIDTH OF

WORKING FACE DISTRIBUTION AND EVOLUTIONARY CHARACTERISTICS OF “MULTI-FIELD” OF SURROUNDING ROCKS WITH DIFFERENT WORKING FACE WIDTHS * (1) Evolutionary law of mining stress distribution

According to the excavation steps shown in Fig. 9, seven cases were simulated when the width of 1802 first mining face was 150, 175, 200, 225, 250, 275 and 300 m. The stress distribution of

the surrounding rock of the working face under the seven different widths of the 1802 first mining face was extracted respectively, as shown in Fig. 9. From the Fig. 9, it can be seen that

as the width of the 1802 working face increases from 150 to 300 m, the surrounding rock stress increases as a whole, and the peak stresses of the coal seam are 54.15, 62.31, 62.87, 62.49,

63.32, 67.76, 76.99 MPa in the following order, the peak stresses of the surrounding rock in the X = 3230 m section are 49.44, 52.26, 54.88, 57.28, 59.48, 61.53 and 63.37 MPa. Overall, the

surrounding rock stress gradually increased with the increase of the width of the working face, but the distribution characteristics of the surrounding rock stress in the working face did

not change significantly. * (2) Distribution and evolutionary law of plastic zone The distribution status of the plastic zone in the working face under 7 different widths of 1802 first

mining face was extracted respectively, as shown in Fig. 10. It can be seen from the plane distribution map of plastic zone, 1802 working face boundary plastic area is mainly distributed in

the range of 5–20 m from the goaf boundary, the plastic zone near the first cutting of the working face and the stopping line in the direction of the strike is smaller, and the plastic zone

in the central mining boundary of the working face is larger, and comparing the distribution of the plastic zone of the working face surrounding rock under different width conditions, it can

be seen that, with the increase of the width, the distribution characteristics of the plastic zone in the direction of the working face have not changed significantly. According to the

distribution of plastic zone profile in Fig. 10, it can be seen that the roof and floor of the goaf are affected by the mining activities and plastic failure occurs to a certain extent. The

plastic failure area of the floor under the goaf is within 15–20 m from the goaf boundary, mainly distributed in the mudstone bottom plate and thin sandstone bottom plate under the goaf. The

depth of the plastic zone under the goaf changes with the thickness of the mudstone and thin sandstone bottom plate, and the thicker the bottom plate, the deeper the plastic zone extends

downward. The plastic zone of the surrounding rock above the goaf is also mainly concentrated in the complex rock formation with the thickness of a single layer less than 10 m above Coal 8

and the thick sandstone with the thickness of the first layer above Coal 8 exceeding 10 m. Plastic failure areas in the rock layers above the thick sandstone are sparsely distributed, and

plastic failure rock masses that develop towards the interior of the goaf only appear above the boundaries on both sides of the goaf. Comparative analysis of the distribution of the plastic

zone above the 1802 goaf under different width conditions shows that with the increase of the width of the working face the plastic zone in the middle of the tendency of the goaf has not

changed significantly, and the height of the plastic zone development is 20–30 m above the coal seam; However, the plastic zone above the two sides of the goaf gradually increases and

develops towards the middle of the tendency of the goaf. When the width of the working face increases from 150 to 300 m, the development height of the plastic zone on both sides of the goaf

increases from 34 m above the goaf to 93 m, and develops towards the middle of the goaf by about 41 m. * (3) Evolutionary law of elastic energy distribution The elastic energy density

distribution states of the working face at seven different widths of the 1802 first mining face were extracted separately, as shown in Fig. 11. It can be seen from the planar distribution

diagram of elastic energy density of coal seam that the elastic energy of surrounding rock is mainly concentrated in the range of 20–200 m from the boundary of 1802 goaf, and the range of

elastic energy concentration along the strike direction of the working face near the open-off cut and the stopping line is smaller, and the range of elastic energy concentration in the

surrounding rock in the middle of the working face at the boundary of the goaf is larger. From the distribution diagram of elastic energy density profile, it can be seen that the elastic

energy density is mainly concentrated in the coal8 and its roof and floor rock layer, and the elastic energy density of the surrounding rock of other rock layer groups is small, mainly

because the elastic modulus of the coal mass is smaller than the elastic modulus of the rock mass. The peak elastic energy of coal is mainly concentrated at the goaf boundary on both sides

of the working face, and the peak elastic energy is about 15–20 m away from the goaf boundary. With the increase of the width of the working face, the concentration of the elastic energy of

the surrounding rock increases, but the characteristics of the distribution of the elastic energy have not changed significantly. DETERMINATION OF THE REASONABLE WIDTH OF THE WORKING FACE In

order to analyze the influence result of the working face width on the surrounding rock stress, the vertical stress curve distributed along the coal seam at the section where X = 3230 m is

extracted, as shown in Fig. 12. At the position where Y = 2165 m, it is the south goaf boundary of the 1802 working face, and at the position where Y = 2315–2465 m, it is the north goaf

boundary of the 1802 working face affected by geological factors such as the distance between the coal seam thickness and the hard roof, and the buried depth of the coal seam, as the width

of the 1802 working face increases from 150 to 300 m, the surrounding rock stress concentration area and the stress peak position continue to migrate and evolve. Therefore, there is no

obvious rule for the stress peak of the surrounding rock on the side of Y = 2315–2465 m, and the stress peak of the surrounding rock of the goaf boundary at Y = 2165 m increases from 43.43

to 55.75 MPa, indicating that the wider the width of the working face, the more intense the overlying roof rock stratum movement, and the greater the lateral bearing stress in the goaf.

Extract the evolution characteristics of the elastic strain energy distribution of the surrounding rock of the goaf boundary on the south side of the 1802 working face (Y = 2165 m, fixed

boundary) at the horizontal position of section X = 3230 m and Z = 115 m, as shown in Fig. 13. From Fig. 13a, it can be known that the stress distribution of the surrounding rock at the +

115 m level is basically consistent with the stress distribution characteristics along the coal seam, and the concentrated stress increases with the increase of the width of the working

face.When the width of the working face increases from 150 to 300 m, the vertical stress peak values are 47.7, 50.4, 53.0, 55.3, 57.5, 59.6, 61.5 MPa respectively, and the stress peak value

is about 20 m away from the goaf boundary.From Fig. 13b, it can be known that the elastic energy change of the surrounding rock of the 1802 goaf boundary is similar to the evolution

characteristics of the vertical stress distribution. As it extends from the goaf boundary to the solid coal, the elastic energy of the surrounding rock increases rapidly, and it reaches the

peak of elastic energy within the range of about 10–15 m away from the goaf boundary, with a peak value of about 400–640 kJ/m3. After that, the elastic energy of the surrounding rock rapidly

decreases within the range of 20–30 m away from the goaf boundary. When the distance between the surrounding rock and the goaf boundary is greater than 30 m, the rate of reduction of the

elastic strain energy slows down; when the distance between the surrounding rock and the goaf boundary is about 150 m, the increase in the elastic strain energy of the surrounding rock

compared with before mining has decreased to a relatively small value, indicating that the influence of the 1802 goaf on the stress and elastic strain energy distribution of the coal and

rock mass outside 150 m is relatively small. As the width of the working face increases from 150 to 300 m, the position of the peak elastic energy shifts from 10 m away from the goaf

boundary to 15 m away from the goaf boundary. The peak elastic strain energy of the goaf boundary is 400, 444, 486, 526, 566, 604, 640 kJ/m3 respectively, and the wider the working face, the

more concentrated the elastic energy of the surrounding rock. Extract the distribution and evolution characteristics of the elastic strain energy of the surrounding rock at the goaf

boundary on the north side (variation boundary) of the 1802 working face at the section of X = 3230 m and Z = 107 m, as shown in Fig. 14. From Fig. 14a, it can be known that as the width of

the working face increases from 150 to 300 m, the stress of the surrounding rock at the Z = 107 m level as a whole shows an increasing trend, and the stress peaks are 46.5, 48.4, 48.8, 50.7,

51.7, 53.3, and 55.4 MPa in turn. From Fig. 14b, it can be known that affected by geological factors such as the buried depth of the coal seam, the elevation of the coal seam, and the

distance from the thick and hard roof, the stress and elastic energy concentration area at the goaf boundary on the northside continuously migrates and evolves with the change of the width

of the working face. As the width of the working face increases from 150 to 300 m, the elastic energy density of the surrounding rock at the Z = 107 m level does not show a clearly

increasing or decreasing rule, and the elastic energy peaks are 420.3, 392.8, 372.6, 425.6, 419.9, 456.3, and 504.6 kJ/m3 in turn. By comparing the distribution of surrounding rock stress

and elastic energy at the goaf boundaries on both sides of the 1802 working face (Figs. 13 and 14), it can be known that the distance between the coal seam and the thick sandstone at the

goaf boundary on the southside of the working face is generally smaller than that on the northside of the working face, and the concentration degree of surrounding rock stress and elastic

energy at the goaf boundary on the south side (Y = 2165 m) of the working face is generally higher than that on the north side (Y = 2315–2465 m). In a word, as the width of the working face

increases, the concentration degree of the surrounding rock stress and elastic energy also goes up. The stress peak is about 15–20 m away from the goaf boundary, and the elastic energy peak

is about 20 m away from the goaf boundary. The lateral bearing stress influence range of the 1802 working face is approximately 200–250 m. Bas Based on the impact of the working face width

on the impact hazard in the impact ground pressure risk index evaluation form affected by mining technical factors, when the working face width is less than 100 m, the impact hazard

evaluation index is 3 (the highest); when the working face width is 100–150 m, the impact hazard evaluation index is 2; when the working face width is 150–300 m, the impact hazard evaluation

index is 1; when the working face width is greater than 300 m, the impact hazard evaluation index is 0. The higher the evaluation index, the higher the impact hazard. Therefore, when

arranging the working face, it is necessary to avoid the working face width causing a higher impact hazard evaluation index. From this perspective, the working face width should be as large

as possible. On the other hand, when considering the stress concentration factor of the surrounding rock and the development height of the plastic zone, it can be observed that an increase

in the working face width corresponds to a higher peak stress in the coal seam and a greater development height of the plastic zone. This indicates that the movement of the overlying rock

above the coal seam becomes more pronounced, thereby intensifying the degree of mining disturbance. From Fig. 15, it can be known that the peak stress of the dip section of the working face

approximately increases linearly with the increase of the working face width; the peak stress in the coal seam shows 3 step-like increases with the increase of the working face width, which

are located at the working face width of 175, 275, and 300 m respectively, and the difference in the peak stress of the coal seam when the working face width is 175–250 m is relatively

small; at the same time, the developing height of the plastic zone above the coal seam increases significantly with the increase of the working face width at 225–250 m; through comprehensive

analysis, it is considered that the reasonable width of the working face can betaken as greater than 175 m and less than 250 m. THE EVOLUTION LAW OF “MULTIPLE FIELDS” IN THE MINING OF FIRST

WORKING FACE Based on the mine working face connection plan and the layout of the first mining working face 1802, under the condition that the width of the first mining working face is 200 

m, based on the three-dimensional numerical model of the panel, the evolution laws of the surrounding rock stress field, energy field, etc. during the recovery process of the 1802 working

face are simulated and studied. THE STRESS DISTRIBUTION CHARACTERISTICS OF THE SURROUNDING ROCK DURING THE “SQUARING” PROCESS OF THE FIRST MINING FACE Based on the design of the 1802 working

face, the working face “squared” when it was excavated to 200 m. The 1802 working face was excavated sequentially to 100, 150, 200, 250, and 300 m to study the evolution characteristics of

the stress distribution of the surrounding rock before and after the working face “squared”. The stress distribution at the goaf boundary is shown in Fig. 16. From Fig. 16, it can be known

that during the “squaring” process of the 1802 working face, as the length of the working face excavation increases from 100 to 300 m, the degree of stress concentration at the goaf boundary

and the range of the peak zone gradually increase. By studying the evolution law of the stress distribution of the surrounding rock on the four sides of the working face boundary, it can be

known that as the recovery length increases, the stress peak and the high-stress zone gradually shift from the goaf boundaries on the east and west sides to the goaf boundaries on the north

and south sides, that is, from the advanced coal wall of the working face to the lateral coal body. In order to further analyze the evolution and transition process of the stress of the

surrounding rock on the four sides of the goaf boundary and the stress peak zone, the stress distribution curves of the goaf boundary on the side of the working face open-off cut and the

southside are extracted, as shown in Fig. 17. By comparing the evolution characteristics of the stress distribution of the surrounding rock on the side of the open-off cut and the southside

boundary of the goaf, it can be known before the working face “squared”, as the excavation length increases from 100 to 200 m, the stress peak of the surrounding rock on the side of the

open-off cut is 36.5, 39.0, and 40 MPa in turn, and the stress peak of the south side goaf boundary is 29.7, 34.4, and 38.8 MPa in turn. At this time, the excavation length is less than the

width of the working face, the side of the open-off cut is the long side of the goaf, and the southside is the short side, so the degree of stress concentration of the surrounding rock on

the side of the open-off cut and the mining field boundary is greater than that on the north and south sides of the goaf boundary. When the length of the working face excavation is 250 and

300 m, and the length of the excavation is greater than the width of the working face, the stress peak of the surrounding rock on the side boundary of the open-off cut is 40.6 and 41.0 MPa

in turn, and the stress peak of the surrounding rock on the south side of the goaf is 42.4 and 45 MPa in turn. The side of the open-off cut is the short side of the goaf, and the southside

is the long side. The surrounding rock stress and peak value shift from the side of the open-off cut to the north and southside boundaries. From Fig. 17a and b, it can be known that with the

increase of the excavation length, the distribution characteristics of the surrounding rock stress on the side of the open-off cut and the position of the stress peak do not change, but the

stress gradually increases and the increase rate gradually decreases, indicating that the movement of the overlying rock on the side of the open-off cut tends to be stable; the stress peak

of the surrounding rock on the goaf boundaries on the north and south sides of the working face maintains a large growth, and the position of the stress peak zone gradually shifts with the

increase of the mining length, and the whole is located in the middle of the goaf boundary. THE DISTRIBUTION AND EVOLUTION CHARACTERISTICS OF THE “MULTIPLE FIELDS” IN THE MINING PROCESS OF

THE FIRST WORKING FACE * (1) The evolution law of the distribution of mining-induced stress The tendency width of the 1802 working face is 200 m, and the strike length is as long as 4100 m.

The strike direction span of the working face is large, and the ground layers fluctuate and change, and the geological conditions are relatively complex. Therefore, it is necessary to study

the distribution and evolution characteristics of the surrounding rock stress, elastic energy, and plastic zone during the mining process of the 1802 working face. Figure 18 is the plane

distribution and profile distribution diagram of the surrounding rock stress during the mining process of the working face. Analyzing Fig. 18, it can be known that affected by the

distribution of stress in the panel area, during the mining process of the 1802 working face, the overall peak value of the advanced bearing stress of the working face (not higher than 45 

MPa) is smaller than the peak value of the lateral bearing stress of the working face. It can be seen from the figure that the elevation of Coal 8 at the section of Y = 2265 m, the buried

depth of the working face, and the ground elevation, etc. continuously change along the strike of the 1802 working face. When the working face is mined from 600 to 3900 m, the peak values of

the section stress are 42.82, 43.61, 43.60, 43.59, and 43.59 MPa in turn, indicating that although the length of the working face mining is constantly increasing, under the influence of

geological factors such as the buried depth of the coal seam, the distance between it and the thick and hard roof, and the thickness of the loesslayer, the peak value of the bearing stress

in front of the stope undergoes a certain degree of fluctuation but the change is not obvious. * (2) The evolution law of the distribution of elastic energy The evolution characteristics of

the distribution of elastic energy of the surrounding rock during the mining process of the 1802 working face are shown in Fig. 19. From the plan view of the distribution of elastic energy,

it can be known that as the length of the working face mining increases, the range of concentrated elastic energy of the surrounding rock at the goaf boundary of 1802 also gradually

increases. When the length of the working face mining increases to a certain extent, the concentration degree and range of elastic energy of the surrounding rock at the goaf boundary of 1802

tend to be stable. From the section view of the distribution of elastic energy, it can be known that the elastic energy in the model is mainly distributed in Coal 8 and the top and bottom

floor rock layers nearby, and the peak value of concentrated elastic energy is distributed within the range of both ends of the goafin 1802. With the increase of the length of the working

face mining, there is no obvious change in the distribution characteristics of the surrounding rock elastic energy. * (3) The evolution law of the distribution of the plastic zone The

extraction of the distribution of the plastic zone of the surrounding rock during the excavation process of the 1802 working face is shown in Fig. 20. From the plan e distribution map of the

plastic zone, it can be known that as the working face continues to be excavated, the range of the plastic zone on both sides of the 1802 goaf gradually increases and tends to be stable,

and the lateral plastic zone range of the goaf is about 5–20 m. From the sectional view of the distribution of the plastic zone, it can be known that the range of the plastic zone of the

surrounding rock above and below the 1802 goaf does not show the rule of increasing or decreasing with the excavation of the working face. When the working face increases from excavating 600

 m to excavating 1000 m, and from excavating 1800 m to excavating 2600 m, the development height of the plastic zone above the 1802 goaf increases significantly, increasing from 36 to 54 m,

and from 21 to 90 m respectively. However, when it increases from excavating 2600 m to excavating 3400 and 3900 m, the range of the plastic zone above the goaf is somewhat reduced. Combined

with the stratum distribution of the coal seam and itstop and bottom floor and the research results of the width of the working face, it can be known that the range of the plastic zone is

mainly concentrated and distributed in the complex stratum group composed of the rock layer with a single layer thickness less than 10 m above and below the coal seam. The thicker the

stratum group, the larger the plastic zone, and the thinner the stratum group, the smaller the plastic zone. THE CHARACTERISTIC OF THE ABUTMENT STRESS DISTRIBUTION ZONE IN THE MINING PROCESS

OF THE FIRST MINING FACE In conjunction with the geological conditions, a further analysis is conducted on the distribution and evolutionary characteristics of the stress field and energy

field of the surrounding rock during the mining of the 1802 face. Extract the distribution curves of the advanced abutment stress and elastic strain energy in the middle of the 1802 working

face (Y = 2265) when the mining length in the numerical study is sequentially 100, 150, 200, 250, 300, 600, 1000, 1400, 1800, 2200, 2600, 3000, 3600, 3850, 3900, 3950, 4000, 4050, 4100 m

(from the initial mining to the end of the working face mining), as shown in Fig. 21, the statistics of the peak stress value and its advanced distance in each mining stage are as shown in

Table 3. Analyzing Fig. 21 and Table 3, it can be known that in the early stage of mining in the 1802 working face, as the width of the working face mining increases, the disturbance of the

working face mining on the overlying rock formation becomes increasingly larger, the activity of the roof rock formation is more intense, and the advanced bearing stress of the surrounding

rock shows an overall increasing rule. After the 1802 working face is mined to “become square”, the distribution of the peak bearing stress of the working face is relatively stable, and the

influence of the increase in the mining length on the peak value of the advanced bearing stress gradually decreases. Statistical analysis shows that the peak values of the advanced bearing

stress and the elastic strain energy are mainly distributed at the position of 10–20 m in advance of the working face. The advanced position of the peak bearing stress increases gradually

from 10 to 20 m as the width of the working face mining increases. By comparing the size of the stress peak of the coal seam and the size of the original rock stress, the size of the peak

bearing stress in advance is between 31 and 40 MPa, and the concentration degree is about 1.45–1.7. It should be noted that in this numerical study, the size of the numerical simulation

model is large, and the overall size of the model grid is large. Therefore, there will be a certain deviation between the calculated results of the bearing stress advanced influence range

and the actual situation on site. By analyzing the curvature change of the advanced bearing stress curve, the advanced bearing stress influence range is approximately 150–200 m. CONCLUSION

AND PROSPECT The geological condition of the first pan area of Xinzhuang Coal mine is complicated, which is affected by many factors such as syncline structure, buried depth of coal seam,

thickness of surface soil layer, thickness of coal seam, inclination of coal seam and hard roof. In this study, the west wing of Coal 8 Pan area was selected as a typical representative, and

a large-scale 3D fine numerical model was established including the above complex factors, in order to obtain the characteristics of the surrounding rock supporting stress distribution zone

during the mining process of the working face under complex conditions through numerical simulation, and provide support for regional anti-scour design. The main conclusions are as follows:

* (1) The study reveals the distribution characteristics of mining stress, elastic energy accumulation and plastic zone distribution at different widths of the first mining face, and

determines the reasonable width of the working face The stress concentration degree and range, elastic energy concentration degree and range, plastic zone range of coal seam increase with

the increase of the width of the first mining face. When the width of the first mining face increases from 150 to 300 m, the peak stress of the coal seam increases from 54.15 to 76.99 MPa,

the elastic energy increases from 400 to 640 kJ/m3, the development height of the plastic zone increases from 34 to 93 m, and these parameters increase significantly when the width of the

working face is greater than 250 m. According to the requirement of impact risk assessment index on the width of the working face, it is considered that the width of the working face should

be greater than 175 m but less than 250 m, and the designed working face width of 200 m inXinzhuang Coal mine is more reasonable. * (2) The evolutional laws of mining stress, elastic energy

accumulation and plastic zone distribution and their influencing factors are revealed in the study The results show that the leading bearing stress is greater than the lateral bearing stress

before the working face is squared, while the leading bearing stress is smaller than the lateral bearing stress after the working face is squared. During the mining process of the first

mining face, there are three high mining dynamic concentration areas near the cut hole, near the strike middle syncline axis and in the second half of the working face, which are not

completely consistent with the location of the stress area of the high plain rock, mainly affected by the layer spacing between the hard overlying roof and the coal seam. The smaller the

layer spacing, the higher the mining stress, and the two are negatively correlated. In general, the concentration degree of mining-induced stress at the head mining face is 1.45–1.7, the

peak stress is located at the depth of the coal body 10–20 m, the influence range of mining-induced stress is 150–200 m, the concentration degree of elastic energy is 1.7–2.5, and the

concentration area of elastic energy is the coal body and its adjacent top and floor strata within 20–200 m of the goaf boundary. The peak elastic energy is located in the depth of coal body

about 20 m, the plastic zone range is 5–20 m, the development height of plastic zone is about 20–90 m, and the development depth of plastic zone under the coal seam is about 8–30 m. In the

future, the research in this field can be carried out from the following directions. In terms of theoretical research, the mining stress mechanics model of surrounding rock under

multi-factor coupling is further improved, and the nonlinear mechanical behavior of rock, the influence of geological structure and the interaction between different rock layers are

considered, so as to improve the accuracy and applicability of the theoretical model. In terms of numerical simulation, more advanced numerical simulation software and algorithms are

developed to improve the simulation accuracy of complex geological physical behavior, and fully consider the characteristics of rock mass heterogeneity, anisotropy and rheology, so that the

simulation results are closer to the actual situation. DATA AVAILABILITY The datasets used and/or analysed during the current study available from the corresponding author on reasonable

request. REFERENCES * He, S. Q. et al. Multi-parameter integrated early warning model and an intelligent identification cloud platform of rockburst. _J. China Univ. Min. Technol_ 51, 850–862

(2022). Google Scholar  * Xie, H. P. et al. Theoretical and technological conception of the fluidization mining for deep coal resources. _J. China Coal Soc._ 42, 547–556 (2017). * Qi, Q. X.

et al. Theory and technical framework of prevention and control with different sources in multi-scales for coal and rock dynamic disasters in deep mining of coal mines. _J. China Coal Soc._

7, 1801–1810 (2018). Google Scholar  * Dou L. M., He X. Q. _Impact Ore Pressure Prevention Theory and Technology_ (China Univ. Min. Univ., 2001). Google Scholar  * Qi, Q. X., Dou, L. M. 

_Intellectual Pressure Theory and technology_ (China Univ. Min. Univ., 2008). Google Scholar  * Cao, A. Y. et al. Dynamic manifestation mechanism of mining on the island coalface along fault

and dynamic pre-warning of seismic waves with seismic tomography. _J. Min. Saf. Eng._ 34, 411–417 (2017). Google Scholar  * Pan, J. F. et al. The top board “artificial liberation layer” to

prevent the impact ground pressure method, mechanism and application. _J. China Coal Soc._ 48, 636–648 (2023). * Wang, H. W. et al. Abrupt effect of coupling instability of fault overburden

induced by mining. _J. Coal Sci._ 48, 2961–2975 (2023). Google Scholar  * He, X. Q. et al. Typical nearly upright coal seam group impact ground pressing mechanism, monitoring and early

warning research. _Tunn. Undergr. Space Technol._ 49, 13–22 (2021). Google Scholar  * Pan, J. F. et al. Pressing the theory and technology of the static loading of static loads in the deep

lanes. _J. China Coal Soc. _45, 1607–1613 (2020). * Tan, Y. L. et al. Deep coal lane gangs to seduce the mechanism and “unload-solid” collaborative control research. _J. China Coal_ _Soc._ 

45, 66–81 (2020). Google Scholar  * Yang, R. S. et al. Model experiment on dynamic behavior of jointed rock mass under blasting at high-stress conditions. _Tunn. Undergr. Space Technol._ 74,

145–152 (2018). Article  Google Scholar  * Yang, L. Y. & Ding, C. X. Fracture mechanism due to blast-imposed loading under high static stress conditions. _Int. J. Rock Mech. Min. Sci._

107, 150–158 (2018). Article  Google Scholar  * Zhu, Q. J. et al. The design of the evaluation system of the coal mine impacting system and its implementation. _J. North China Univ. Sci.

Technol._ 19, 1–7 (2022). ADS  Google Scholar  * Wang, C. W. et al. Analysis on control action of geologic structure on rock burst and typical cases. _J. China Coal Soc._ 37, 263–268 (2012).

CAS  Google Scholar  * Wang, H. B. et al. Wudong coal mine urgent sanded coal sea anti-bringing technology. _Autom. Operating_ 48, 58–61 (2022). MathSciNet  Google Scholar  * Xu, W.

Research on monitoring technology and application of surrounding rock stress in mining space. China University of Mining and Technology (2012). * Si, X. Rock Burst simulation test and

mechanism study of deep high stress hard rock tunnel surrounding rock. Central South University (2022). Download references FUNDING This research was funded by the National Natural Science

Foundation of China (Grant No. 52064042), the China Postdoctoral Science Foundation (Grant No. 2023MD734191), the Program for Young Talents of Science and Technology in Universities of Inner

Mongolia Autonomous Region (Grant No. NJYT24040), the Fundamental Research Funds for Inner Mongolia University of Science & Technology (Grant No. 2023QNJS112), the Innovation Fund

Project of Inner Mongolia University of Science & Technology (Grant No. 2019QDL-B34), the Open Fund Project of the State Key Laboratory of Deep Geotechnical Mechanics and Underground

Engineering of China University of Mining and Technology (Beijing) (SKLGDUEK2020), the Huaneng Group Headquarters Science and Technology Project (HNKJ20-H32, HNKJ21-H56-10). AUTHOR

INFORMATION AUTHORS AND AFFILIATIONS * China University of Mining and Technology, Xuzhou, 221000, China Xingen Ma * School of Mining and Coal Engineering, Inner Mongolia University of

Science and Technology, Baotou, 014010, China Eryu Wang & Qiangdi Geng * Huaneng Coal Technology Research Co., Ltd., Beijing, 100070, China Xingen Ma, Bing Hu, Yongyuan Li, Fulong Sun, 

Bingshan Wang & Wei Wang * Inner Mongolia Yitai Group Co., Ltd., Ordos, 017000, China Eryu Wang * Huaneng Qingyang Coal Power Co., Ltd, Qingyang, 745000, China Zhiwen Liu * School of

Mechanics and Civil Engineering, China University of Mining and Technology (Beijing), Beijing, 100083, China Eryu Wang Authors * Xingen Ma View author publications You can also search for

this author inPubMed Google Scholar * Bing Hu View author publications You can also search for this author inPubMed Google Scholar * Eryu Wang View author publications You can also search

for this author inPubMed Google Scholar * Qiangdi Geng View author publications You can also search for this author inPubMed Google Scholar * Zhiwen Liu View author publications You can also

search for this author inPubMed Google Scholar * Yongyuan Li View author publications You can also search for this author inPubMed Google Scholar * Fulong Sun View author publications You

can also search for this author inPubMed Google Scholar * Bingshan Wang View author publications You can also search for this author inPubMed Google Scholar * Wei Wang View author

publications You can also search for this author inPubMed Google Scholar CONTRIBUTIONS Conceptualization, X.M.; methodology, E.W.; software, X.M. and E.W.; validation, W.W. and B.W.; formal

analysis, X.M. and E.W.; investigation, X.M. and Q.G.; resources, Z.L.; data curation, Q.G.; writing—original draft preparation, X.M. and E.W.; writing—review and editing, X.M. and E.W.;

visualization, Y.L.; supervision, B.H and F.S.; project administration, X.M. and B.H; funding acquisition, E.Y and X.M. All authors have read and agreed to the published version of the

manuscript. CORRESPONDING AUTHOR Correspondence to Eryu Wang. ETHICS DECLARATIONS COMPETING INTERESTS The authors declare no competing interests. ADDITIONAL INFORMATION PUBLISHER’S NOTE

Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations. RIGHTS AND PERMISSIONS OPEN ACCESS This article is licensed under a

Creative Commons Attribution-NonCommercial-NoDerivatives 4.0 International License, which permits any non-commercial use, sharing, distribution and reproduction in any medium or format, as

long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if you modified the licensed material. You do not

have permission under this licence to share adapted material derived from this article or parts of it. The images or other third party material in this article are included in the article’s

Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article’s Creative Commons licence and your intended use is not

permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this licence, visit

http://creativecommons.org/licenses/by-nc-nd/4.0/. Reprints and permissions ABOUT THIS ARTICLE CITE THIS ARTICLE Ma, X., Hu, B., Wang, E. _et al._ Study on mining-induced stress behavior on

surrounding rock of rockburst mine under large buried depth and hard overburden. _Sci Rep_ 15, 16819 (2025). https://doi.org/10.1038/s41598-025-00368-w Download citation * Received: 19

September 2024 * Accepted: 28 April 2025 * Published: 14 May 2025 * DOI: https://doi.org/10.1038/s41598-025-00368-w SHARE THIS ARTICLE Anyone you share the following link with will be able

to read this content: Get shareable link Sorry, a shareable link is not currently available for this article. Copy to clipboard Provided by the Springer Nature SharedIt content-sharing

initiative KEYWORDS * Large buried depth * Thick topsoil * Hard overburden * Rockburst mine * Mining-induced stress behavior