Investigation into Occurrence Mechanism of Rock Burst Induced by Water Drainage in Deep Mines

: Conﬁned aquifers widely exist in the strata of the Ordos mining area. Water drainage before mining is an effective measure to avoid water inrush disasters caused by the connectivity between mining-induced roof fractures and aquifers. However, rock burst disasters occur frequently in the mining process of many water drainage working faces. The statistics show that the surrounding rock of nine water drainage working faces in ﬁve mines has suffered different degrees of rock burst, which seriously restricts the safe, efﬁcient production and sustainable development of ten-million-ton modern mines in China. Based on this, taking the 22,106 longwall working face (LW22106) of the Shilawusu Mine, Ordos, as the engineering background, this paper investigates the occurrence mechanism of water drainage on rock burst in the roof water-rich area by using theoretical analysis, similar material simulation, microseismic measurement and other methods. The main conclusions are as follows: (1) After the drainage of the water-rich area, the pressure relief zone, pressurized zone and pressure stabilization zone are formed in sequence from the center to both sides. The width of the pressure relief zone is consistent with that of the water-rich area, and the width of the pressurized zone is about 35 m on one side. (2) When the mining is passing the pressurized zone, the coal rock mass is under the joint inﬂuence of gravity stress, mining disturbance stress and drainage transfer stress. The superimposed stress generally exceeds 2.5 times the uniaxial compressive strength (UCS) of the coal rock mass, and the maximum reaches 3.24 times, far exceeding the critical value of rock burst (1.5 times UCS), which is the main reason to induce rock burst. (3) The dynamic change in the rock burst risk areas before and after drainage in water-rich area no. 4 has been predicted, and the number of risk areas increases from 4 before drainage to 13 after drainage. Since the stress superposition effect differs due to different mining speeds, it is proposed that the dynamic regulation of mining speed in the front and at the back of the drainage area is an effective and efﬁcient method for rock burst prevention and control. The research results will provide a theoretical basis and technical support for the prevention and control in the roof water-rich area of deep mines.


Introduction
With the gradual depletion of shallow coal resources in China, coal resource exploitation presents a deep mining tendency, and deep mining is increasingly becoming normal [1,2].After coal mines step into the deep mining stage, geological conditions and stope structures become much more complex, the occurrence frequency and damage degree of dynamic disasters remarkably increase, and more and more mines will be faced with the threat of dynamic disasters such as rock burst [3][4][5][6].Since 2010, 25  have occurred in China, resulting in 119 deaths and 108 injuries.Rock burst has become a key scientific problem and a worldwide problem that needs to be solved urgently [7][8][9][10][11].Ordos coalfield, as an important coal base in China, stretches across Shaanxi, Gansu, Ningxia, Inner Mongolia and Shanxi provinces, as shown in Figure 1.Group Yan'an and Group Zhiluo artesian aquifers widely exist in the coal seam roof, and the large water-rich areas serve as the direct water source for coal seam mining [12].In the mining process on the working face below the roof water-rich areas, roof disasters such as water inrush and sand bursting may be faced [13].In order to avoid the roof water inrush and reduce the preparation time on the working face, the roof water drainage should be carried out during excavation on the working face.However, the drainage in the water-rich area causes the stress redistribution of the coal seam and roof [14], forming a stress concentration area and increasing the possibility of rock burst.In recent years, rock burst has appeared several times in the mining process of some roof water drainage working faces, which seriously restricts the safe and efficient production and coal supply.Therefore, the investigation into the influence mechanism of roof water drainage on the stress distribution of overlying rock is of great significance to reveal the occurrence mechanism of rock burst on the water drainage working face in the Ordos mining area, take measures for rock burst prevention and control and ensure the safety of the working face.
Sustainability 2023, 15, x FOR PEER REVIEW 2 of 21 degree of dynamic disasters remarkably increase, and more and more mines will be faced with the threat of dynamic disasters such as rock burst [3][4][5][6].Since 2010, 25 rock burst accidents have occurred in China, resulting in 119 deaths and 108 injuries.Rock burst has become a key scientific problem and a worldwide problem that needs to be solved urgently [7][8][9][10][11].Ordos coalfield, as an important coal base in China, stretches across Shaanxi, Gansu, Ningxia, Inner Mongolia and Shanxi provinces, as shown in Figure 1.Group Yan'an and Group Zhiluo artesian aquifers widely exist in the coal seam roof, and the large water-rich areas serve as the direct water source for coal seam mining [12].In the mining process on the working face below the roof water-rich areas, roof disasters such as water inrush and sand bursting may be faced [13].In order to avoid the roof water inrush and reduce the preparation time on the working face, the roof water drainage should be carried out during excavation on the working face.However, the drainage in the water-rich area causes the stress redistribution of the coal seam and roof [14], forming a stress concentration area and increasing the possibility of rock burst.In recent years, rock burst has appeared several times in the mining process of some roof water drainage working faces, which seriously restricts the safe and efficient production and coal supply.Therefore, the investigation into the influence mechanism of roof water drainage on the stress distribution of overlying rock is of great significance to reveal the occurrence mechanism of rock burst on the water drainage working face in the Ordos mining area, take measures for rock burst prevention and control and ensure the safety of the working face.Field practice shows that rock burst is affected by many factors, including geological and mining factors.Researchers have conducted extensive studies on the influencing factors of rock burst and made rich scientific research achievements.Jiang et al. [2,15] believed that rock burst, affected by various inducing factors, was a nonlinear dynamic process in the deformation of the coal rock mass system under different geological conditions and mining environments.During the process, energy stored in a stable state was transformed into energy released in an unstable state.They also pointed out that the essence of rock burst was the stress-induced sudden instability and failure of the coal rock mass.Dou et al. [16,17] and Cai et al. [18] held that rock burst was closely related to the physical and Field practice shows that rock burst is affected by many factors, including geological and mining factors.Researchers have conducted extensive studies on the influencing factors of rock burst and made rich scientific research achievements.Jiang et al. [2,15] believed that rock burst, affected by various inducing factors, was a nonlinear dynamic process in the deformation of the coal rock mass system under different geological conditions and mining environments.During the process, energy stored in a stable state was transformed into energy released in an unstable state.They also pointed out that the essence of rock burst was the stress-induced sudden instability and failure of the coal rock mass.Dou et al. [16,17] and Cai et al. [18] held that rock burst was closely related to the physical and mechanical properties of coal rock, as well as the static and dynamic loads on coal rock.Based on this, they further investigated the energy and stress conditions of rock burst induced by dynamic and static loads and clarified the mechanism and two types of rock burst induced by static and dynamic loads: high static load type and strong dynamic load type.Qi et al. [19] explained many determining factors for the occurrence of rock burst and put forward the idea of multi-scale and multi-source prevention and control of dynamic disasters in deep coals.Jiang et al. [20] and Zhu et al. [21,22] found that gravity and tectonic stress were the foundation stress of determining the rock burst in the stope, and the stress generated by the movement of overlying strata was the force source that induces the rock burst.They believed that the variable structure forms and movement laws of overlying strata were the main reasons for the complex mechanism and various forms of rock burst and finally raised the technical idea of classification evaluation on rock burst.Zhao et al. [23] studied the loading characteristics and instability mechanism of gob-side roadway in deep mining and revealed the "three-load" superimposition rock-burst-inducing mechanism of gob-side roadway under the influence of lateral static load, advanced static load and disturbance dynamic load in goaf.Cao et al. [24][25][26] conducted an investigation into the stress distribution law and multi-dimensional information of microseismic events near the residual coal pillar, holding that the residual coal pillar was the main influencing factor of rock burst.
In terms of the studies on the relationship between mining disturbance and rock burst, Gong et al. [27,28] found that the critical loading rate obviously existed in the coal-rock combined body, and the mutation rate of rock burst tendency was relatively obvious near the critical loading rate.Zhao et al. [29] illustrated the mechanism of roof energy accumulation and release under the influence of mining speed, believing that the increase in roof energy accumulation was essentially the increase in loading rate when the mining speed changed from slow to fast.Huang et al. [30] conducted uniaxial compression tests on the coal-rock combined body under different loading rates, explaining that the loading rate had a huge effect on failure forms and the rock burst tendency of coal rocks.Cui et al. [31] evaluated the rock burst risk of a coal body under different mining strengths and revealed the mining disturbance characteristics on the working face of deep coal seams under different mining strengths.Feng et al. [32] adopted theoretical analysis and microseismic monitoring to reveal how mining speed influenced the energy released by the movement of the hard roof, indicating that the energy from the hard roof breaking presented an obvious positive correlation with the mining speed.In the aspect of drainage-induced rock burst, scholars have conducted some research studies on the effect of water on the physical and mechanical properties of coal rock, but few studies focus on how drainage affects the stress field of coal and rock.Zhou et al. [33] explored the mechanical behavior of coal samples during water saturation, finding that when the coal around the mining space was subjected to water, the high degree of water saturation in the coal decreased the risks of coal bursts significantly.Zheng et al. [34] pointed out that there was a negative correlation between the compressive strength of coal and soaking time by studying the influence of soaking time on the macro and micro mechanical properties of coal samples.Zhang et al. [35][36][37] focused on the water weakening of rock strength and its influencing factors (water content, immersion time and wetting-drying cycles), and they found that the strength of the rock mass decreases to varying degrees with water content, immersion time and wetting-drying cycles depending on the rock mass type and mineral composition.
It can be found that the existing rock burst theories usually ignore water drainage as an inducing factor of rock burst, which leads to the lack of reliability and accuracy in the classification of rock burst risk areas on the working face with the roof water-rich areas and the absence of targeted prevention measures.For example, the serious rock burst accident in the drainage lane and no. 3 contact lane of LW1303 in Shandong Longyun Mining Company on 20 October 2018 caused 21 deaths [5] and damaged a 348 m roadway, including a 198 m severely damaged section, as shown in Figure 2. According to the accident investigation report, the influence of roof water drainage on the regional stress distribution was one of the indirect causes.The result showed that roof water drainage was one of the induced factors of rock burst in Ordos deep mining area.
accident investigation report, the influence of roof water drainage on the regional stress distribution was one of the indirect causes.The result showed that roof water drainage was one of the induced factors of rock burst in Ordos deep mining area.In view of this, taking LW22106 of the Shilawusu Mine in Ordos as the engineering background, this paper adopts theoretical analysis, similar material simulation and microseismic monitoring to investigate the occurrence mechanism of water drainage on rock burst in the roof water-rich area of deep mines and reclassified the risk areas on the working face, attempting to predict the rock burst risk more accurately.The research results can provide a theoretical basis and technical support for safe and efficient mining on the roof water-rich working face in the Ordos mining area.

Field Conditions of LW22106
The Shilawusu Mine, located in the Ordos mining area, is mainly mining the 2-2 coal seam, whose coal seam, roof and floor have a weak rock burst tendency.LW22106, with a width of 300 m and a length of 1120 m, is the first mining face of the south wing in the 221 mining area, and a strike longwall full-mechanized coal caving method is adopted for mining on this working face.The average surface elevation is 1340 m, the average mining depth is 660 m, and the average coal thickness is 9.02 m.Like a flat seam, the 2-2 coal seam has a simple geological structure.The immediate roof of the 2-2 coal seam is 0.77 m thick mudstone, the basic roof is 31.26m thick fine sandstone, the direct bottom is 6.10 m thick sandy mudstone, and the basic bottom is 9.25 m thick siltstone.According to the hydrogeological exploration results, there is a Group Zhiluo artesian aquifer with an average thickness of 110 m at 39.65~62.65 m and an average of 50.52 m above LW21106 of the 2-2 coal seam.The water-rich areas are unevenly distributed above the working face.The specific numbers and spatial distribution locations of the four water-rich areas are shown in Figure 3.In view of this, taking LW22106 of the Shilawusu Mine in Ordos as the engineering background, this paper adopts theoretical analysis, similar material simulation and microseismic monitoring to investigate the occurrence mechanism of water drainage on rock burst in the roof water-rich area of deep mines and reclassified the risk areas on the working face, attempting to predict the rock burst risk more accurately.The research results can provide a theoretical basis and technical support for safe and efficient mining on the roof water-rich working face in the Ordos mining area.

Field Conditions of LW22106
The Shilawusu Mine, located in the Ordos mining area, is mainly mining the 2-2 coal seam, whose coal seam, roof and floor have a weak rock burst tendency.LW22106, with a width of 300 m and a length of 1120 m, is the first mining face of the south wing in the 221 mining area, and a strike longwall full-mechanized coal caving method is adopted for mining on this working face.The average surface elevation is 1340 m, the average mining depth is 660 m, and the average coal thickness is 9.02 m.Like a flat seam, the 2-2 coal seam has a simple geological structure.The immediate roof of the 2-2 coal seam is 0.77 m thick mudstone, the basic roof is 31.26m thick fine sandstone, the direct bottom is 6.10 m thick sandy mudstone, and the basic bottom is 9.25 m thick siltstone.According to the hydrogeological exploration results, there is a Group Zhiluo artesian aquifer with an average thickness of 110 m at 39.65~62.65 m and an average of 50.52 m above LW21106 of the 2-2 coal seam.The water-rich areas are unevenly distributed above the working face.The specific numbers and spatial distribution locations of the four water-rich areas are shown in Figure 3.
During the mining process on the working face, the water drainage has been carried out through drainage holes.Specifically speaking, the drilling fields and the drainage holes are intensively arranged in the vicinity of the water-rich area, with one drilling field for every 200 m.The diameter of drilling holes for roof water drainage is 91 mm, and the final pitch of the holes is less than 80 m.After the drilling is completed, the gate valve and pressure gauge are installed at the orifice.The gate valve is connected with the drainage pipe and flowmeter for drainage water until the water pressure in the water-rich area is reduced up to the safety mining standard, which is that roof water disasters such as water inrush and sand bursting will not occur in the mining process.The microseismic monitoring system is SOS, developed in Poland [38].Eight microseismic sensors are arranged around LW22106, among which no. 5 and no.6 sensors are constantly moved and optimized along with the mining on the working face.These sensors are used to monitor microseismic events in the mining process and capture the precursor information of rock burst.During the mining process on the working face, the water drainage has been carried out through drainage holes.Specifically speaking, the drilling fields and the drainage holes are intensively arranged in the vicinity of the water-rich area, with one drilling field for every 200 m.The diameter of drilling holes for roof water drainage is 91 mm, and the final pitch of the holes is less than 80 m.After the drilling is completed, the gate valve and pressure gauge are installed at the orifice.The gate valve is connected with the drainage pipe and flowmeter for drainage water until the water pressure in the water-rich area is reduced up to the safety mining standard, which is that roof water disasters such as water inrush and sand bursting will not occur in the mining process.The microseismic monitoring system is SOS, developed in Poland [38].Eight microseismic sensors are arranged around LW22106, among which no. 5 and no.6 sensors are constantly moved and optimized along with the mining on the working face.These sensors are used to monitor microseismic events in the mining process and capture the precursor information of rock burst.

Description of Mine Seismicity on the Working Face
At 19:00 on 10 October 2019, the rock burst appeared in the normal mining process of LW21106, and the staff reported the big sound of coal blasting on the site.After further on-site investigation, obvious large rib splitting occurred on the coal wall of the working face, a large amount of slag fell off the roof in front of the working face, and part of the coal body in the roadway wall leaked out from the anchor net.Meanwhile, the SOS microseismic monitoring system monitored a large energy event of 4.0 × 10 5 J (source 1) and another large energy event of 1.13 × 10 5 J (source 2) at 20:06, with a time interval of 66 min.In addition to the sound of coal blasting, there was also a small amount of slag falling off the roof.By comparing the real-time online monitoring results of stress on the working face, it was found that the stress of the two gate roads before and after the incident did not change greatly, and the detection hole of pulverized coal presented a normal amount of drilling cuttings in the field construction.The plane projection and profile projection of the events were shown in Figure 4.

Description of Mine Seismicity on the Working Face
At 19:00 on 10 October 2019, the rock burst appeared in the normal mining process of LW21106, and the staff reported the big sound of coal blasting on the site.After further on-site investigation, obvious large rib splitting occurred on the coal wall of the working face, a large amount of slag fell off the roof in front of the working face, and part of the coal body in the roadway wall leaked out from the anchor net.Meanwhile, the SOS microseismic monitoring system monitored a large energy event of 4.0 × 10 5 J (source 1) and another large energy event of 1.13 × 10 5 J (source 2) at 20:06, with a time interval of 66 min.In addition to the sound of coal blasting, there was also a small amount of slag falling off the roof.By comparing the real-time online monitoring results of stress on the working face, it was found that the stress of the two gate roads before and after the incident did not change greatly, and the detection hole of pulverized coal presented a normal amount of drilling cuttings in the field construction.The plane projection and profile projection of the events were shown in Figure 4.As can be seen from Figure 4, Source 1 was horizontally located 12 m in front of the working face and at the edge of the no. 4 water-rich area and vertically located in the basic roof (fine sandstone) above the coal seam with 27.45 m away from the 2-2 coal roof.Source 2 was located in the coal seam.According to preliminary analysis, the water pressure in the no. 4 rich water area declined from the initial 3 MPa to 0.5 MPa after drainage.The pressure step interval of LW22106 was 20 m, and the mining speed reached 7.2 m/d on that day.Therefore, it was believed that the main reasons for this rock burst were as follows: the overlying basic roof failed to collapse in time due to the excessively fast mining As can be seen from Figure 4, Source 1 was horizontally located 12 m in front of the working face and at the edge of the no. 4 water-rich area and vertically located in the basic roof (fine sandstone) above the coal seam with 27.45 m away from the 2-2 coal roof.Source 2 was located in the coal seam.According to preliminary analysis, the water pressure in the no. 4 rich water area declined from the initial 3 MPa to 0.5 MPa after drainage.The pressure step interval of LW22106 was 20 m, and the mining speed reached 7.2 m/d on that day.Therefore, it was believed that the main reasons for this rock burst were as follows: the overlying basic roof failed to collapse in time due to the excessively fast mining speed on the working face, thus forming a suspended roof and causing energy accumulation; meanwhile, because of the no. 4 water-rich area on the overlying, a stress concentration appeared near the water-rich area after drainage, which increased the stress above the basic roof, accelerated the fracture of the basic roof and produced a large energy event.Moreover, the energy accumulated in the suspended roof could not be quickly released and stored in the coal body in front of the working face.The coal body above the working face was affected by the superposition of the advanced abatement pressure and the drainage transfer stress.The superimposed stress exceeded the own strength of the coal body, thus resulting in its destruction, cracks and energy release and inducing the rock burst on the working face, such as large rib splitting.

Principle of Effective Stress
There are irregular water-rich areas in the roof of the upper coal seam in dozens of newly built mines with a mining depth of more than 550 m in Ordos.In order to avoid water inrush and reduce preparation time on the working face during the mining process, the water drainage during mining is equal to a mine for a "quasi-liberating seam" in the water-rich area [39].The stress decreases in the drainage area and increases at the edge.Therefore, when the working face is advancing to the area affected by drainage, rock burst is more easily induced.It is of great necessity to explore how drainage affects the abutment pressure of the working face.
According to the column chart of the roof strata and the on-site pumping operation, the confined artesian aquifer of Group Zhiluo above the LW22106 coal seam mostly includes medium sandstones and fine sandstones, with the maximum value of confined water pressure up to 3 MPa.The artesian aquifer is saturated with high porosity and permeability, and the confined artesian aquifer of Group Zhiluo meets the principle of effective stress.The principle of effective stress means that the total stress carried by the saturated rock-soil mass is mainly borne by two parts [40].One is the effective stress, namely the skeleton formed by solid particles, which is transmitted through the contact surface between particles; the other is the pore water pressure, namely the water in the pores.The principle of effective stress can be expressed by Formula (1): where σ is the total stress of the saturated rock-soil mass (MPa); σ is the effective stress (MPa); and µ is the pore water pressure (MPa).
The bearing force of the artesian aquifer is borne by rock skeleton and pore water.After the drainage in the water-rich area, the pore water pressure and effective stress of the coal rock strata and the physical and mechanical properties of the coal rock mass will make changes within a certain range, resulting in the regional stress adjustment and redistribution of the coal seam and roof.

Stress Migration Law in the Process of Drainage
The above analysis demonstrates that the pore water pressure will gradually decrease and the effective stress of the artesian aquifer will change after the drainage in the waterrich area of the roof.The support of the artesian aquifer to the overlying strata will decrease within the drainage range, and the overlying strata will suffer from deformation and subsidence under the influence of its gravity.The rock mass on both sides of the water-rich area will provide the supporting reaction force.That is to say, the overlying strata load originally borne by the pore water pressure in the water-rich area is transferred to both sides of the water-rich area, which causes the stress redistribution of the coal seam and roof strata, further forming the pressure relief zone, pressurized zone and pressure stabilization zone.To further explain, the stress of the coal seam and roof reduces within the drainage area, the stress concentrates at the edge of the area, and the range not affected by drainage is the initial stress zone.
In order to investigate the influence of drainage on the stress distribution of the coal seam and roof before and after drainage in the water-rich area, taking LW22106 of the Shilawusu Mine in Ordos as an example, a section is made along the strike of the working face, and the corresponding estimation mechanical model of abutment pressure after drainage is established to analyze the stress state above the working face before and after drainage, as shown in Figure 5.
crease and the effective stress of the artesian aquifer will change after the drainage in the water-rich area of the roof.The support of the artesian aquifer to the overlying strata will decrease within the drainage range, and the overlying strata will suffer from deformation and subsidence under the influence of its gravity.The rock mass on both sides of the water-rich area will provide the supporting reaction force.That is to say, the overlying strata load originally borne by the pore water pressure in the water-rich area is transferred to both sides of the water-rich area, which causes the stress redistribution of the coal seam and roof strata, further forming the pressure relief zone, pressurized zone and pressure stabilization zone.To further explain, the stress of the coal seam and roof reduces within the drainage area, the stress concentrates at the edge of the area, and the range not affected by drainage is the initial stress zone.
In order to investigate the influence of drainage on the stress distribution of the coal seam and roof before and after drainage in the water-rich area, taking LW22106 of the Shilawusu Mine in Ordos as an example, a section is made along the strike of the working face, and the corresponding estimation mechanical model of abutment pressure after drainage is established to analyze the stress state above the working face before and after drainage, as shown in Figure 5.The following assumptions are made for drawing the general law: (1) for the model, the deformation of coal seam and roof strata meets three basic assumptions of material mechanics; (2) for the stress, the water level and water pressure in the area before and after drainage are uniform, all strata are uniformly stressed, and the differences in stress transfer caused by different rock joints are ignored; (3) for the boundary, only an ideal The following assumptions are made for drawing the general law: (1) for the model, the deformation of coal seam and roof strata meets three basic assumptions of material mechanics; (2) for the stress, the water level and water pressure in the area before and after drainage are uniform, all strata are uniformly stressed, and the differences in stress transfer caused by different rock joints are ignored; (3) for the boundary, only an ideal section of the working face is selected for the theoretical analysis because of the irregular shape in the actual water-rich area.In Figure 5a, before the implementation of drainage in the water-rich area, the strata of the artesian aquifer are intact, the overlying strata load q is uniformly transferred, and the load on the artesian aquifer is the original rock stress q.In Figure 5b, the origin O is the center of the water-rich area.After drainage in the water-rich area, the pore water pressure decreases, forming a pressure relief zone.The support capacity of the upper roof in the drainage range of the water-rich area decreases to q 1 (A A section), and the rock mass stress in A C and AC sections increases, forming a pressurized zone.The stress distribution in the pressurized zone can be simplified into a triangular linear distribution (the maximum stress is q 2 ).The original rock stress in the C D and CD sections remains unchanged, forming a pressure stabilization zone.Let the length of the OA (OA ) section be l, AC (A C ) section be a and the CD (C D ) section be b, and then the stress distribution in D D section is as follows: where q is the original rock stress (MPa); q 1 is the load in the pressure relief zone after drainage (MPa); q 2 is the maximum load in the pressurized zone (MPa); l is the half-width of the water-rich area (m); a is the range of the pressurized zone, expressed as a = h 2 / tan α (m); and b is the range of the pressure stabilization zone (m).
As can be seen from Figure 5 and Formula (2), the stress distribution after drainage is ideally symmetric with respect to the origin O, which means equal stress is transferred to both sides after drainage.Therefore, it just needs to make a detailed analysis of the stress distribution in the OD section on the right side of the origin O, as shown in Figure 5c.In Figure 5c, h 1 is the distance from the top of the artesian aquifer to the coal seam (m); h 2 is the distance from the top of the artesian aquifer to the surface (m); and α is the displacement angle of overlying strata in the water-rich area caused by drainage.(Compared with mining, the support of the artesian aquifer on its upper strata still exists in spite of a decline in support effect.Moreover, the overlying strata over the aquifer will suffer from fracture bands, but caving bands will not occur or the overlying strata will not stay in a completely overhanging state).
In this model, the water-rich area is simplified to be composed of confined water and a rock skeleton.The rock skeleton can be regarded as an equivalent spring, exerting a supporting reaction force on the overlying strata.The supporting reaction force of the water-rich area on the overlying strata before and after drainage is expressed as: where γ is the rock bulk density (kN/m 3 ) and k is the drainage pressure relief coefficient (0 < k < 1), negatively correlated with the drainage degree; that is, the higher degree of the drainage, the smaller the k.
After drainage, the gravity stress in the overlying OAFG region is transferred to the right of ∆ ACE , and the transferred weight Q OAFG can be expressed as: Assuming that stress transferred to the lower coal rock mass is approximately an isosceles triangular distribution [20], the load transferred to the right of the OAFG region σ OAFG is as follows: where x is the distance between the transfer stress and the center of the water-rich area.The peak stress transferred from the OAFG region to the pressurized zone q 2 can be expressed as: When the working face is under the water-rich area (pressure relief zone), at the edge of the water-rich area (pressurized zone) and in a more distant area (pressure stabilization zone), the gravity stress of the working face from the overlying strata after drainage differs greatly, which can be expressed as: Based on Formulas (2)-( 7), the distribution function of coal seam abutment pressure caused by drainage can be written: As can be seen from Formula (8), the stress concentration caused by drainage shows a negative correlation with the pressure relief coefficient in the water-rich area and is positively correlated with the drainage degree.That is to say, the better the drainage effect, the greater the stress peak in the pressurized zone.According to the research [12,14] and the field observation of water pressure, the water pressure in the rich-water area gradually decreases to the safe mining standard (that is, no water disasters occur, such as water inrush) and eventually tends to be stable with the progress of drainage work.The results show that the stress concentration degree caused by drainage tends to be stable during mining on the working face; the pressure relief zone, pressurized zone and pressure stabilization zone are formed in sequence from the center to both sides; and the pressurized zone applies additional stress to the coal seam and its roof, adding a static stress source for the occurrence of rock burst.

Physical Simulation Verification
Starting from the typical geological conditions of the Ordos mining area, the similarity ratios are determined based on the rock strata condition, physical and mechanical properties of rock [41], geometric size of the test bed and three basic similarity laws.Among them, the geometric ratio is 1/200, the bulk density ratio is 0.6, the strength ratio is 0.6/200, and the time ratio is 1/14.For the convenience of time simulation, the time ratio is set as 1/12; that is, 2 hours in the experiment equals the on-site 24 h.The model size is 5000 mm (length) × 250 mm (width) × 1680 mm (height).Sand is used as the skeleton material, and calcium carbonate and gypsum are used as cementing materials.Table 1 displays the proportions of similar simulation materials.A pressurized steel plate is arranged at the top of the model, and eight hydraulic jack devices are uniformly arranged above the steel plate to exert vertical forces on the model, which replaces the rock layer above that needs to be simulated.Water-rich areas exist in the coal seam roof.In order to analyze the stress evolution law before and after drainage, the no. 4 medium sandstone part in Table 1 is taken as an artesian aquifer.Several water pipes are placed in the aquifer during the building of the model.The water pipe is 2.1 m long, 1.4 m from the left boundary and 1.5 m from the right boundary of the model.The water pipe is provided with a cut-off valve, and the tail end is close to the rear steel plate for the convenience of drainage, as shown in Figure 6a.Ten UEI pressure boxes with high precision are staggered under the artesian aquifer and coal seam to monitor the stress transfer of the coal seam and roof during and after drainage, as shown in Figure 6b.The spacing between the test model and stress sensor is shown in Figure 6c.In order to study the influence of water drainage on coal seam mining, a 170 cm coal pillar is set in the middle to avoid the influence of excavation on both sides of the working face.The mining sequence is as follows: (1) first, excavate LW101 from the position 20 cm to 170 cm away from the right boundary of the model, with a total length of 150 cm; (2) secondly, open the stop valve of water pipe to release all of the water, and stand for a week after the completion of water drainage; (3) finally, excavate LW201 from the position 20 cm to 160 cm away from the left boundary of the model, with a total length of 140 cm.Before the initial excavation, the stress monitoring instrument is used to clear the sensor stress to zero for comparison.
Figure 7a shows the stress changing curve of no. 8 and no. 9 measuring points after drainage.It can be seen that the stress at the no. 9 measuring point under the water-rich area shows an upward trend after drainage, with a sharp increase from the initial 0 kPa to 13 kPa between 0.5 h and 1 h after drainage, and eventually stabilizes to 16 kPa after 5 h.Meanwhile, the stress at the no.8 measuring point at the edge of the water-rich area gradually declines from the initial 0 kPa to −7 kPa after 6 h.The drainage in the water-rich area will lead to the stress change in the coal seam, namely the stress decrease under the waterrich area and the stress increase at the edge of the water-rich area, thus forming the pressurized zone and pressure relief zone.Figure 7b is the stress changing curve of no.1~no.4 measuring points after drainage.It can be found that no apparent stress changes have taken place at no. 1 and no. 2 measuring points after drainage, which are far away from the water-rich area, and the change has been floating up and down within the allowable error range.The stress at no. 3 and no. 4 measuring points under the water-rich area decreases from 0 kPa to −1.75 kPa and from 0 kPa to −2 kPa, respectively.The stress changing curve of no. 8 and no. 9 measuring points displays that the influence range affected by drainage stress will not expand indefinitely, and the stress concentration simply appears in a certain range at the edge of the water-rich area.In order to study the influence of water drainage on coal seam mining, a 170 cm coal pillar is set in the middle to avoid the influence of excavation on both sides of the working face.The mining sequence is as follows: (1) first, excavate LW101 from the position 20 cm to 170 cm away from the right boundary of the model, with a total length of 150 cm; (2) secondly, open the stop valve of water pipe to release all of the water, and stand for a week after the completion of water drainage; (3) finally, excavate LW201 from the position 20 cm to 160 cm away from the left boundary of the model, with a total length of 140 cm.Before the initial excavation, the stress monitoring instrument is used to clear the sensor stress to zero for comparison.
Figure 7a shows the stress changing curve of no. 8 and no. 9 measuring points after drainage.It can be seen that the stress at the no. 9 measuring point under the water-rich area shows an upward trend after drainage, with a sharp increase from the initial 0 kPa to 13 kPa between 0.5 h and 1 h after drainage, and eventually stabilizes to 16 kPa after 5 h.Meanwhile, the stress at the no.8 measuring point at the edge of the water-rich area gradually declines from the initial 0 kPa to −7 kPa after 6 h.The drainage in the water-rich area will lead to the stress change in the coal seam, namely the stress decrease under the water-rich area and the stress increase at the edge of the water-rich area, thus forming the pressurized zone and pressure relief zone.Figure 7b is the stress changing curve of no.1~no.4 measuring points after drainage.It can be found that no apparent stress changes have taken place at no. 1 and no. 2 measuring points after drainage, which are far away from the water-rich area, and the change has been floating up and down within the allowable error range.The stress at no. 3 and no. 4 measuring points under the water-rich area decreases from 0 kPa to −1.75 kPa and from 0 kPa to −2 kPa, respectively.The stress changing curve of no. 8 and no. 9 measuring points displays that the influence range affected by drainage stress will not expand indefinitely, and the stress concentration simply appears in a certain range at the edge of the water-rich area.Figure 8a shows the stress changing curves at the no.6 measuring point with LW101 excavated for 90~95 cm and at the no. 2 measuring point with LW201 excavated for 60~65 cm.At this time, the no.6 measuring point is 100 cm away from the open cutting of LW101 and 5 cm away from the advancing position of LW101; the no. 2 measuring point is 70 cm away from the open cutting hole of LW201 and 5 cm away from the advancing position of LW201.It can be seen that the two measuring points basically keep a consistent stress growth trend, and the initial stress of the no. 2 measuring point is only 2.5 kPa lower than that of the no.6 measuring point, which is within the allowable error range of the instrument.Figure 8a shows the stress changing curves at the no.6 measuring point with LW101 excavated for 90~95 cm and at the no. 2 measuring point with LW201 excavated for 60~65 cm.At this time, the no.6 measuring point is 100 cm away from the open cutting of LW101 and 5 cm away from the advancing position of LW101; the no. 2 measuring point is 70 cm away from the open cutting hole of LW201 and 5 cm away from the advancing position of LW201.It can be seen that the two measuring points basically keep a consistent stress growth trend, and the initial stress of the no. 2 measuring point is only 2.5 kPa lower than that of the no.6 measuring point, which is within the allowable error range of the instrument.Figure 8a shows the stress changing curves at the no.6 measuring point with LW101 excavated for 90~95 cm and at the no. 2 measuring point with LW201 excavated for 60~65 cm.At this time, the no.6 measuring point is 100 cm away from the open cutting of LW101 and 5 cm away from the advancing position of LW101; the no. 2 measuring point is 70 cm away from the open cutting hole of LW201 and 5 cm away from the advancing position of LW201.It can be seen that the two measuring points basically keep a consistent stress growth trend, and the initial stress of the no. 2 measuring point is only 2.5 kPa lower than that of the no.6 measuring point, which is within the allowable error range of the instrument.Figure 8b shows the stress changing curves at the no. 10 measuring point with LW101 excavated for 110~115 cm and at the no.8 measuring point with LW201 excavated for 100~105 cm.At this time, the no. 10 measuring point is 120 cm distant from the open cutting of LW101 and 5 cm from the advancing position of LW101; the no.8 measuring point is 110 cm from the open cutting hole of LW201 and 5 cm from the advancing position of LW201.It can be obtained that the two measuring points almost present the same stress growth trend, and the initial stress of the no.8 measuring point is 19.5 kPa higher than that of the no. 10 measuring point.When the two working faces are advanced to 5 cm away from the measuring point, the stress of the no.8 measuring point is 15 kPa higher than that of the no. 10 measuring point.The stress difference is similar to the stress increment of 16 kPa at the no.8 measuring point after drainage in Figures 7a and 8c,d demonstrates the spatial relationship between the working face and the stress measuring point when LW101 is excavated for 110~115 cm and LW201 for 100~105 cm.
Above all, the influence range of drainage stress will not expand indefinitely, and stress concentration only appears in a certain range at the edge of the water-rich area.Compared with the working face not affected by water drainage, when the working face advances to the drainage-affected area, the mining disturbance stress and drainage transfer stress will be superimposed.If the superimposed stress reaches the critical stress, rock burst will be induced, which verifies the correctness of the theoretical calculation model.

Analysis of the Relationship between Drainage and Rock Burst
In order to evaluate the rock burst risk at point N, the rock burst risk evaluation index I c [42] proposed by our research team is introduced.I c stands for the ratio of vertical stress at point N to the uniaxial compressive strength of coal rock mass, expressed as: where σ N is the vertical stress at point N and [σ c ] is the uniaxial compressive strength of coal rock mass.The classification of rock burst risk is shown in Table 2.In the mining process of the roof water drainage working face, the superimposed stress on the coal body varies according to the position of point N (as shown in Figure 9).It can be obtained that the two measuring points almost present the same stress growth trend, and the initial stress of the no.8 measuring point is 19.5 kPa higher than that of the no. 10 measuring point.When the two working faces are advanced to 5 cm away from the measuring point, the stress of the no.8 measuring point is 15 kPa higher than that of the no. 10 measuring point.The stress difference is similar to the stress increment of 16 kPa at the no.8 measuring point after drainage in Figures 7a and 8c,d demonstrates the spatial relationship between the working face and the stress measuring point when LW101 is excavated for 110~115 cm and LW201 for 100~105 cm.Above all, the influence range of drainage stress will not expand indefinitely, and stress concentration only appears in a certain range at the edge of the water-rich area.Compared with the working face not affected by water drainage, when the working face advances to the drainage-affected area, the mining disturbance stress and drainage transfer stress will be superimposed.If the superimposed stress reaches the critical stress, rock burst will be induced, which verifies the correctness of the theoretical calculation model.

Analysis of the Relationship between Drainage and Rock Burst
In order to evaluate the rock burst risk at point N, the rock burst risk evaluation index Ic [42] proposed by our research team is introduced.Ic stands for the ratio of vertical stress at point N to the uniaxial compressive strength of coal rock mass, expressed as: where N σ is the vertical stress at point N and [ ] c σ is the uniaxial compressive strength of coal rock mass.The classification of rock burst risk is shown in Table 2.In the mining process of the roof water drainage working face, the superimposed stress on the coal body varies according to the position of point N (as shown in Figure 9).The following situations typically exist: (1) affected by gravity stress and mining disturbance stress (I, V); (2) affected by gravity stress and pressure relief from drainage and pressurization from mining disturbance under the water-rich area; and (3) affected by The following situations typically exist: (1) affected by gravity stress and mining disturbance stress (I, V); (2) affected by gravity stress and pressure relief from drainage and pressurization from mining disturbance under the water-rich area; and (3) affected by gravity stress and pressurization from drainage and mining disturbance at the edge of the water-rich area (II, IV).The stress state at point N can be approximated as: where σ N is the superimposed stress at point N, σ M is the mining disturbance stress at point N, and x is the horizontal distance between point N and the center of the water-rich area.
According to Formula (10), the vertical stress at point N is related to the mining speed and the drainage degree.When it is affected by both drainage pressurization and mining disturbance pressurization, the stress condition of rock burst induction can easily be reached.
The above analysis explains the rock burst mechanism induced by water drainage in the roof water-rich area of deep mines.That is to say, when the working face is passing the pressurization area by drainage, the coal rock mass is under the combined influence of gravity stress, mining disturbance stress and drainage transfer stress, and the superimposed stress easily exceeds the critical value of rock burst, which is the main reason to induce rock burst.

Field Application
Water drainage is taken before mining on LW21106.Water pressure in the water-rich area decreases from 3 MPa at the beginning to 0.5 MPa after drainage.After the drainage is completed, the pressure relief zone and pressurized zone are formed in front of the working face, which facilitates the changes in the rock burst risk areas on the working face.In this case, the rock burst risk areas can be re-predicted and divided according to the location of risk areas.
Due to the weak rock burst liability of seam 2-2 in the Shilawusu Mine, the rock burst risk assessment of LW22106 has been carried out before mining.Field practice and research show that the peak value of mining disturbance stress is generally 1.5~1.8times the original rock stress under the condition of a constant mining speed.It can be calculated that the peak value of mining disturbance stress is about 24.75~29.7 MPa with an average of 27.225 MPa, and the rock burst risk index I c = 27.225/17.6= 1.55.Thus, the working face is determined to have a weak rock burst risk.Taking the no. 4 water-rich area as an example, the multi-factor coupling evaluation method is adopted to divide the rock burst risk areas around the no. 4 water-rich area, including three weak rock burst risk areas and one medium rock burst risk area, as shown in Figure 10 and Table 3.
Sustainability 2023, 15, x FOR PEER REVIEW 15 of 21 gravity stress and pressurization from drainage and mining disturbance at the edge of the water-rich area (II, IV).The stress state at point N can be approximated as: [ -] [ ( ) ] ( ) where N σ is the superimposed stress at point N, M σ is the mining disturbance stress at point N, and x is the horizontal distance between point N and the center of the water- rich area.
According to Formula (10), the vertical stress at point N is related to the mining speed and the drainage degree.When it is affected by both drainage pressurization and mining disturbance pressurization, the stress condition of rock burst induction can easily be reached.
The above analysis explains the rock burst mechanism induced by water drainage in the roof water-rich area of deep mines.That is to say, when the working face is passing the pressurization area by drainage, the coal rock mass is under the combined influence of gravity stress, mining disturbance stress and drainage transfer stress, and the superimposed stress easily exceeds the critical value of rock burst, which is the main reason to induce rock burst.

Field Application
Water drainage is taken before mining on LW21106.Water pressure in the water-rich area decreases from 3 MPa at the beginning to 0.5 MPa after drainage.After the drainage is completed, the pressure relief zone and pressurized zone are formed in front of the working face, which facilitates the changes in the rock burst risk areas on the working face.In this case, the rock burst risk areas can be re-predicted and divided according to the location of risk areas.
Due to the weak rock burst liability of seam 2-2 in the Shilawusu Mine, the rock burst risk assessment of LW22106 has been carried out before mining.Field practice and research show that the peak value of mining disturbance stress is generally 1.5~1.8times the original rock stress under the condition of a constant mining speed.It can be calculated that the peak value of mining disturbance stress is about 24.75~29.7 MPa with an average of 27.225 MPa, and the rock burst risk index Ic = 27.225/17.6= 1.55.Thus, the working face is determined to have a weak rock burst risk.Taking the no. 4 water-rich area as an example, the multi-factor coupling evaluation method is adopted to divide the rock burst risk areas around the no. 4 water-rich area, including three weak rock burst risk areas and one medium rock burst risk area, as shown in Figure 10 and Table 3.However, the evaluation result does not take into account the influencing factors of drainage in the water-rich area.Therefore, it can be believed that the evaluation result is a classification method of rock burst risk areas simply considering the mining disturbance before drainage.
After the completion of water drainage, the pressure relief zone and pressurized zone are formed in front of the working face, which results in the change in the rock burst risk areas on the working face.According to the actual situation of LW22106, the values are as follows: h 1 = 160 m, where h 1 is defined as the distance from the top of the artesian aquifer to the coal seam; h 2 = 500 m, where h 2 represents the distance from the surface to the top of the artesian aquifer; 2l 1 = 40 m, where 2l 1 refers to the average width of the narrow area in the no. 4 water-rich area; 2l 1 = 260 m, where 2l 1 is the average width of the wide area; α = 86 • , where α means the fracture angle of the overlying strata in the water-rich area after drainage; k = 0.8, where k is taken as the pressure relief coefficient (the water pressure in the water-rich area is reduced from the initial 3 MPa to 0.5 MPa after drainage, which reports that the pressure relief value is about 20% of the original stress of water-rich artesian aquifer, so k = 1 − 0.2 = 0.8); [σ c ] = 17.6 MPa, where [σ c ] stands for the uniaxial compressive strength of coal; and γ = 25 kN/m 3 .By substituting each parameter into Equation (8), it can be obtained that the stress influence range after drainage is about 35 m, the peak value is about 17.5 m away from the edge of the water-rich area, the stress peak value is 33.4 MPa in the narrow zone of the no. 4 water-rich area, and the stress peak value is 46.3 MPa in the wide zone of the no. 4 water-rich area.
Given the mining disturbance stress and the stress transfer before and after drainage in the no. 4 water-rich area, it can be obtained that: When the working face is advanced to the pressurized zone in the no. 4 waterrich area, the stress peak value σ Nmax1 = 27.225+ 33.4 − 16.5 = 44.125Mpa, I c = 2.51, σ Nmax2 = 27.225+ 46.3 − 16.5 = 57.025Mpa and I c = 3.24.It can be comprehensively determined that the working face has a powerful rock burst risk at this time.
When the working face is advanced to the pressure relief zone in the no. 4 water-rich area, the stress peak value σ Mmax3 = 27.225+ 14 − 16.5 = 24.725MPa and I c = 1.40, which indicates there is no rock burst risk on the working face at this time.
It can be seen that it is easy to induce the rock burst when the working face is advanced to the pressurized zone, which is also consistent with the rock burst location on 10 October.On this basis, the rock burst risk areas after drainage are re-predicted and re-divided for the surrounding area of the no. 4 water-rich area, including two no rock burst risk areas, four weak rock burst risk areas, one medium rock burst risk area and six powerful rock burst risk areas, as shown in Figure 11 and Table 4.
By comparing Figures 10 and 11 and Tables 3 and 4, it can be obtained that the number of risk areas increases from 4 before drainage to 13 after drainage, one in-plane rock burst risk area is added, and the risk degrees change from weak and moderate rock burst before drainage to no, weak, medium and powerful rock burst after drainage.It can also be found that the pressurized zone and the pressure relief zone on the working face and two gate roads before and after drainage in the water-rich area lead to great changes in the rock burst risk zones and risk degrees on the working face.The rock burst risk degree is remarkably reduced within the pressure relief zone, whereas it significantly increases in the pressurized zone.All rock burst locations on 10 October 2019 are in the no.13 in-plane pressurized zone, which further demonstrates that the division of the rock burst risk areas on the working face after drainage is extremely consistent with the field practice.3 and 4, it can be obtained that the number of risk areas increases from 4 before drainage to 13 after drainage, one in-plane rock burst risk area is added, and the risk degrees change from weak and moderate rock burst before drainage to no, weak, medium and powerful rock burst after drainage.It can also be found that the pressurized zone and the pressure relief zone on the working face and two gate roads before and after drainage in the water-rich area lead to great changes in the rock burst risk zones and risk degrees on the working face.The rock burst risk degree is remarkably reduced within the pressure relief zone, whereas it significantly increases in the pressurized zone.All rock burst locations on 10 October 2019 are in the no.13 inplane pressurized zone, which further demonstrates that the division of the rock burst risk areas on the working face after drainage is extremely consistent with the field practice.

Field Verification
Each microseismic event has rich temporal and spatial information [43].Real-time monitoring and processing of microseismic events can identify the occurrence location and intensity of microseismic events, further determining the stress transition inside coal rock masses.Figure 12 shows the relationship between daily footage and daily microseismic energy of LW22106 from 1 September 2019 to 23 October 2019.Figure 13 is the plane projection of microseismic events above the 4th power of LW22106 during this period.It can be seen that mining Zone I from 1 September to 9 September was partly under the water-rich area and partly outside the influence range of drainage, and the microseismic energy just had a small change with the increase in mining speed.From 10 September to 18 September, mining Zone II was partly under the water-rich area and partly within the pressurized zone, and the microseismic energy increased significantly at a high mining speed compared with Zone I. Mining Zone III from 19 September to 4 October was almost under the water-rich area, and the microseismic energy was still very small in spite of a high mining speed (maximum 8.8 m/d).From 5 October to 14 October, most of mining Zone IV was within the pressurized zone, and the high mining speed resulted in the maximum daily energy and a significant increase in total microseismic energy.Most of mining Zone V from 15 October to 23 October was not affected by water drainage, and the microseismic energy was changed little with the changing mining speed.The red line in Figure 13 represents the footage of the working face.It can be found that the microseismic events mostly occurred within the pressurized zone, which proved that this zone had a high degree of concentration with the working face being advanced here.Therefore, rock burst appeared on the working face on 10 October.rock masses.Figure 12 shows the relationship between daily footage and daily microseismic energy of LW22106 from 1 September 2019 to 23 October 2019.Figure 13 is the plane projection of microseismic events above the 4th power of LW22106 during this period.It can be seen that mining Zone I from 1 September to 9 September was partly under the water-rich area and partly outside the influence range of drainage, and the microseismic energy just had a small change with the increase in mining speed.From 10 September to 18 September, mining Zone II was partly under the water-rich area and partly within the pressurized zone, and the microseismic energy increased significantly at a high mining speed compared with Zone I. Mining Zone III from 19 September to 4 October was almost under the water-rich area, and the microseismic energy was still very small in spite of a high mining speed (maximum 8.8 m/d).From 5 October to 14 October, most of mining Zone IV was within the pressurized zone, and the high mining speed resulted in the maximum daily energy and a significant increase in total microseismic energy.Most of mining Zone V from 15 October to 23 October was not affected by water drainage, and the microseismic energy was changed little with the changing mining speed.The red line in Figure 13 represents the footage of the working face.It can be found that the microseismic events mostly occurred within the pressurized zone, which proved that this zone had a high degree of concentration with the working face being advanced here.Therefore, rock burst appeared on the working face on 10 October.rock masses.Figure 12 shows the relationship between daily footage and daily microseismic energy of LW22106 from 1 September 2019 to 23 October 2019.Figure 13 is the plane projection of microseismic events above the 4th power of LW22106 during this period.It can be seen that mining Zone I from 1 September to 9 September was partly under the water-rich area and partly outside the influence range of drainage, and the microseismic energy just had a small change with the increase in mining speed.From 10 September to 18 September, mining Zone II was partly under the water-rich area and partly within the pressurized zone, and the microseismic energy increased significantly at a high mining speed compared with Zone I. Mining Zone III from 19 September to 4 October was almost under the water-rich area, and the microseismic energy was still very small in spite of a high mining speed (maximum 8.8 m/d).From 5 October to 14 October, most of mining Zone IV was within the pressurized zone, and the high mining speed resulted in the maximum daily energy and a significant increase in total microseismic energy.Most of mining Zone V from 15 October to 23 October was not affected by water drainage, and the microseismic energy was changed little with the changing mining speed.The red line in Figure 13 represents the footage of the working face.It can be found that the microseismic events mostly occurred within the pressurized zone, which proved that this zone had a high degree of concentration with the working face being advanced here.Therefore, rock burst appeared on the working face on 10 October.Reducing the mining speed can lower the peak of advanced abutment pressure and the overhanging distance of the basic roof, thus providing reserved time for rock energy release and decreasing the deformation rate of surrounding rock.However, a blind decrease in mining speed will cause huge economic losses.Taking into account the social resource utilization, corporate economic development and safety production, the dynamic control of the mining speed is suggested in the areas affected by water drainage.
In the mining process, according to the mining situation in the normal areas (no fault, no water-rich and other stress abnormal areas), the pressure step distance and preliminary reasonable mining speed can be determined through the microseismic statistical curve and theoretical calculation.When the working face is passing the pressurized zone, the mining speed should be reduced to ensure the stress peak cannot meet the condition of inducing

Figure 1 .
Figure 1.Location map of Ordos coal base in China.

Figure 1 .
Figure 1.Location map of Ordos coal base in China.

Figure 2 .
Figure 2. Scene damage of 20 October rock burst accident in Longyun Mine: (a) Anchor bolt bent and suspended; (b) Two sides convergence, roof subsidence.

Figure 2 .
Figure 2. Scene damage of 20 October rock burst accident in Longyun Mine: (a) Anchor bolt bent and suspended; (b) Two sides convergence, roof subsidence.

Figure 3 .
Figure 3. Distribution plan of sandstone water-rich area on roof of LW22106.

Figure 3 .
Figure 3. Distribution plan of sandstone water-rich area on roof of LW22106.

Figure 5 .
Figure 5. Stress distribution state before and after drainage: (a) before drainage; (b) after drainage; and (c) mechanical model for abutment pressure after drainage.

Figure 5 .
Figure 5. Stress distribution state before and after drainage: (a) before drainage; (b) after drainage; and (c) mechanical model for abutment pressure after drainage.

Figure 6 .
Figure 6.Physical simulation model: (a) water pipes laying; (b) pressure sensors; and (c) test model layout.The different numbers are stress measuring points.

Figure 6 .
Figure 6.Physical simulation model: (a) water pipes laying; (b) pressure sensors; and (c) test model layout.The different numbers are stress measuring points.

FigureFigure 7 .
Figure8bshows the stress changing curves at the no. 10 measuring point with LW101 excavated for 110~115 cm and at the no.8 measuring point with LW201 excavated for

FigureFigure 8 .
Figure8bshows the stress changing curves at the no. 10 measuring point with LW101 excavated for 110~115 cm and at the no.8 measuring point with LW201 excavated for

Sustainability 2023 ,
15,  x FOR PEER REVIEW 14 of 21 100~105 cm.At this time, the no. 10 measuring point is 120 cm distant from the open cutting of LW101 and 5 cm from the advancing position of LW101; the no.8 measuring point is 110 cm from the open cutting hole of LW201 and 5 cm from the advancing position of LW201.

Figure 9 .
Figure 9. Stress environment of coal body at a certain point in the mining process of roof water drainage working face.The I-V is the position of point N.

Figure 9 .
Figure 9. Stress environment of coal body at a certain point in the mining process of roof water drainage working face.The I-V is the position of point N.

Figure 10 .
Figure 10.Division diagram of rock burst risk areas around no. 4 water-rich area before drainage.Figure 10.Division diagram of rock burst risk areas around no. 4 water-rich area before drainage.

Figure 10 .
Figure 10.Division diagram of rock burst risk areas around no. 4 water-rich area before drainage.Figure 10.Division diagram of rock burst risk areas around no. 4 water-rich area before drainage.

Sustainability 2023 , 21 Figure 11 .
Figure 11.Division diagram of rock burst risk areas around no. 4 water-rich area after drainage.

Figure 11 .
Figure 11.Division diagram of rock burst risk areas around no. 4 water-rich area after drainage.

Figure 12 .
Figure 12.Curves of daily advance and microseismic energy.

Figure 13 .
Figure 13.Plane projection of microseismic events above the 4th power of LW22106.

Figure 12 .
Figure 12.Curves of daily advance and microseismic energy.

Figure 12 .
Figure 12.Curves of daily advance and microseismic energy.

Figure 13 .
Figure 13.Plane projection of microseismic events above the 4th power of LW22106.

Figure 13 .
Figure 13.Plane projection of microseismic events above the 4th power of LW22106.

Table 1 .
Proportion of similar simulation materials.

Table 2 .
Relationship between I c and rock burst risk.

Table 2 .
Relationship between Ic and rock burst risk.

Table 3 .
Division table of rock burst risk areas around no. 4 water-rich area before drainage.

Table 4 .
Division table of rock burst risk areas around no. 4 water-rich area after drainage.

Table 4 .
Division table of rock burst risk areas around no. 4 water-rich area after drainage.