The Spatial Evolution Law and Water Inrush Mechanism of Mining-Induced Overburden in Shallow and Short Coal Seam Group

: In order to grasp the overburden evolution law of the shallow and short coal seam group, based on the key bed theory, a mechanical analysis of the longitudinal expansion of mining-induced overburden fractures in the coal seam group was carried out, and the spatial evolution characteristics of mining-induced overburden fractures were simulated by the numerical simulation method. The results show that in the process of working face advancement, with the continuous instability and failure of the overburden, the size and shape of the fracture network are also changed. The repeated mining of the lower coal seam further causes the secondary activation of the upper overburden, which makes the roof fractures of the partially compacted goaf violently move again. The “channel source” and “space source” continue to carry out the process of “generation-expansion-compression-generation-expansion”, in combination with pore fracture elastic theory. The water inrush characteristics of the whole coal seam are divided into three “solid-liquid” coupling stages: the original gap seepage stage, the initial water discharge stage of mining ﬁssures and the water inrush stage of fractured rock mass. The steady value of water inﬂow and its variation characteristics with time are predicted by using the formula of deep well ﬂow in a conﬁned aquifer.


Introduction
The coal seams in the mining areas of central and western China generally have the characteristics of shallow buried depth, thin bedrocks, multi-layer coal and small seam spacing. These characteristics can be collectively referred to as shallow and short coal seam group extract. First, the underlying bedrock thickness of coal seam is small, and the burial depth is no more than 150 m, so it is easy to break through under the influence of coal mining. Secondly, due to the small interval between layers in the short coal seam group, the coal seam mining seams will have a great influence on each other, and the repeated mining will continue to affect the development characteristics of overburdening. Whether there will be a water inrush phenomenon in coal seam mining is mainly affected by the following main control factors: water inrush source, abundance of water source and water pressure, water inrush channel and stagnant water space, which can be referred to as the "four source elements", including the water inrush source, dynamic source, channel source and spatial source. When the "four source elements" constitute a complete relationship, water inrush will occur in the working face. For the problems of overburden failure and water inrush in coal seam mining under shallow buried conditions, Huang Qingxiang [1][2][3] studied the ore pressure characteristics of shallow buried coal seams and proposed the theory of "double key layers" in the roof of a near-shallow buried coal seam with large mining height. Jia Mingkui [4] studied the failure height of the overburdening after mining of thin bedrock working face and made clear the conditions of water and sand gushing. For shallow coal seam mining, Yang Dengfeng [5] studied the instability mechanism of overburdened rock under load using catastrophe theory and analyzed the main influencing factors of overburdened rock instability. Fan Limin et al. [6,7] studied the strata control and water conservation mining technology of coal seam mining under the condition of shallow burial in the ecologically fragile region of western China. Liu Quanming [8] studied the buried depth effect of the overburdened migration rule in fully mechanized caving mining under the condition of shallow buried depth and thin bedrock by using the method of a similar simulation experiment. More research focused on single coal seam mining under the condition of impact, for the shallow buried close distance coal seam group mining roof damage rule and the study of water inrush mechanisms is not very perfect, so that under the condition of close distance coal seam group occurrence, a shallow buried coal seam repeated mining group will cause the overburden rock fracture channels and severe changes of the space, causing a water inrush mechanism different from those of single coal seam mining conditions. Therefore, it is of great significance to study the failure law of roof overburdening and the water inrush mechanisms of shallow and short-distance coal seam group mining ( Figure 1 research methodology flowchart) and grasp its evolution law to prevent the occurrence of water disasters.
Sustainability 2022, 14, x FOR PEER REVIEW 2 of 16 relationship, water inrush will occur in the working face. For the problems of overburden failure and water inrush in coal seam mining under shallow buried conditions, Huang Qingxiang [1][2][3] studied the ore pressure characteristics of shallow buried coal seams and proposed the theory of "double key layers" in the roof of a near-shallow buried coal seam with large mining height. Jia Mingkui [4] studied the failure height of the overburdening after mining of thin bedrock working face and made clear the conditions of water and sand gushing. For shallow coal seam mining, Yang Dengfeng [5] studied the instability mechanism of overburdened rock under load using catastrophe theory and analyzed the main influencing factors of overburdened rock instability. Fan Limin et al. [6,7] studied the strata control and water conservation mining technology of coal seam mining under the condition of shallow burial in the ecologically fragile region of western China. Liu Quanming [8] studied the buried depth effect of the overburdened migration rule in fully mechanized caving mining under the condition of shallow buried depth and thin bedrock by using the method of a similar simulation experiment. More research focused on single coal seam mining under the condition of impact, for the shallow buried close distance coal seam group mining roof damage rule and the study of water inrush mechanisms is not very perfect, so that under the condition of close distance coal seam group occurrence, a shallow buried coal seam repeated mining group will cause the overburden rock fracture channels and severe changes of the space, causing a water inrush mechanism different from those of single coal seam mining conditions. Therefore, it is of great significance to study the failure law of roof overburdening and the water inrush mechanisms of shallow and short-distance coal seam group mining ( Figure 1 research methodology flowchart) and grasp its evolution law to prevent the occurrence of water disasters.

Overview of the Study Area
Xialiyuan Coal Mine belongs to the Pingshuo Mining Area. Coal 4 and 9 of the Taiyuan formation are mainly mined in the mine field. The second mine area's first working face of 4 coal is a 40,201 working face. The extension length is 765 m, the dip length is 150 m and the average thickness of the coal seam is 11.2 m with a dip angle of 2°~5°. The structure of the coal seam is simple and stable, the roof contains K3 and K4 sandstone aquifers, the water yield is weak, the average thickness of the bedrock is 100 m, the

Overview of the Study Area
Xialiyuan Coal Mine belongs to the Pingshuo Mining Area. Coal 4 and 9 of the Taiyuan formation are mainly mined in the mine field. The second mine area's first working face of 4 coal is a 40,201 working face. The extension length is 765 m, the dip length is 150 m and the average thickness of the coal seam is 11.2 m with a dip angle of 2 •~5• . The structure of the coal seam is simple and stable, the roof contains K3 and K4 sandstone aquifers, the water yield is weak, the average thickness of the bedrock is 100 m, the average thickness of the Quaternary is 30 m, the water yield is medium and the spring water is exposed in Sustainability 2022, 14, 5320 3 of 16 the low-lying part of the ground. The average thickness of a 9 coal seam is 8 m, which, under the 4 coal and the coal seam spacing, is 30.50 m. The fully mechanized top coal caving mining method is adopted in the working face. According to the construction drilling results and the laboratory physical and mechanical property test results of the construction rock samples during the mine field exploration (the test is based on DZ/T 0276 rock physical and mechanical property test specification), the statistical regional borehole columnar characteristics and the physical and mechanical parameters of each coal and rock stratum are shown in Table 1. According to the above "four source elements" of water inrush. K3, K4 and the Quaternary aquifer in Xialiyuan coal mine are all water inrush sources. The Quaternary aquifer has a medium water yield and is a stable drastic source. Once qualified channel and spatial sources are produced in coal seam mining, water inrush may occur in the working face. Therefore, by studying the spatial evolution characteristics of mining overburdening in the coal seam, the occurrence process of channel source and spatial source can be obtained, and then the law of water inrush can be mastered and prevented in the future.

"The Channel Source"-Mechanical Analysis of Longitudinal Propagation of Mining-Induced Overburden Fractures in Short Coal Seam Group
The strata above the direct roof of coal seams are the key strata that control the activities of the whole overburden in the stope or to the surface, so the structural position of the key strata can affect the overall development height of the mining-induced overburden. When the main key strata is located under the overburden and close to the coal seam, and is less than a critical distance, because the compressible rotary space resistance of the lower part of the main key strata is relatively large, the subsidence and rotation of the structure block are larger when the key strata is broken, and the fracture of the key strata can lead to the larger opening of the rock fracture of the strata, which extends to the top of the overburden.

Key Strata Structure of Overburden in Coal Seam Group
According to the key strata theory [9][10][11], when in the rock strata above the coal seam (q n+1 ) 1 < (q n ) 1 , the first hard rock layer and the N + 1 rock layer can both be the key strata, controlling the synchronous fracture and subsidence of other soft rock strata. When L n+1 < L 1 , above the coal seam, only the first hard strata is the key strata, which is a single key stratum. When L n+1 > L 1 , there are at least two key layers in the overburden: the upper one controlling the upper layer to the surface of the hard rock as the main key layer and the lower hard rock as the sub-key layer. According to the judgment criterion of key strata, and referring to the columnar structure of boreholes and the physical and mechanical parameters of strata in Table 1, quantitative analysis and judgment are made on all strata in the working faces of Xialiyuan Coal Mine. The judgment results are as follows: the rock No. 10 (medium-grained sandstone, 9.63 m thick) is the main key layer, and the rock No. 14 (medium-grained sandstone, 6.75 m thick) is the sub-key layer.
In the coal seam No. 4 mining, the main key layer (rock No. 10, referred to as key layer 1) plays a dominant role in controlling the activities of the overburden ( Figure 2). When the coal seam No. 4 is exhausted, the coal seam No. 9 is mining, and the sub-key layer (rock layer No. 14, referred to as key layer 2) becomes the main key strata of the overburden of the coal seam No. 9. The key strata structure has changed a single key layer of the upper structure ( Figure 3), which can present typical rock pressure characteristics of strata behavior.

"The Channel Source"-Mechanical Analysis of Longitudinal Propaga ing-Induced Overburden Fractures in Short Coal Seam Group
The strata above the direct roof of coal seams are the key strata that con ities of the whole overburden in the stope or to the surface, so the structur the key strata can affect the overall development height of the mining-indu den. When the main key strata is located under the overburden and close to and is less than a critical distance, because the compressible rotary space res lower part of the main key strata is relatively large, the subsidence and r structure block are larger when the key strata is broken, and the fracture of can lead to the larger opening of the rock fracture of the strata, which exte of the overburden.

Key Strata Structure of Overburden in Coal Seam Group
According to the key strata theory [9][10][11], when in the rock strata above , the first hard rock layer and the N + 1 rock layer can both be controlling the synchronous fracture and subsidence of other soft rock strat L1, above the coal seam, only the first hard strata is the key strata, which i stratum. When Ln+1 > L1, there are at least two key layers in the overburd one controlling the upper layer to the surface of the hard rock as the main the lower hard rock as the sub-key layer. According to the judgment criterion and referring to the columnar structure of boreholes and the physical and m rameters of strata in Table 1, quantitative analysis and judgment are made o the working faces of Xialiyuan Coal Mine. The judgment results are as foll No. 10 (medium-grained sandstone, 9.63 m thick) is the main key layer, and 14 (medium-grained sandstone, 6.75 m thick) is the sub-key layer.
In the coal seam No.4 mining, the main key layer (rock No.10, referred t 1) plays a dominant role in controlling the activities of the overburden (Fig the coal seam No.4 is exhausted, the coal seam No.9 is mining, and the sub-k layer No.14, referred to as key layer 2) becomes the main key strata of the o the coal seam No.9. The key strata structure has changed a single key laye structure ( Figure 3), which can present typical rock pressure characteristic havior.

Failure Mechanism of Key Strata in Coal Seam Group
There are two types of overburden strata failure instability [12]: Structure bility (S instability, shear failure) and structure deformation instability (R insta sion bending failure).

Shear Failure Condition
The shear failure condition of overburden is when the working face is pu cracking end of overburden, the middle part of the strata has not cracked and stress of the remaining section has exceeded the ultimate strength. In this case, is a small amount of space under the strata, the overburden will be cut off a resulting in the overburden being cut off as a whole in the goaf [13,14]. Accor rock beam stress analysis and material mechanics [15], the shear failure expr mula is as follows:

Tensile Bending Failure Condition
The bending failure condition of the overburden on the roof is as follow continuous advance of the stope, the overburden (assumed to be plastic strata) and under the action of overlying weak strata and its own gravity, the strata ar subsided. The original fixed beam is transformed into a simply supported beam maximum bending moment transferred from the end of the beam to the mid

Failure Mechanism of Key Strata in Coal Seam Group
There are two types of overburden strata failure instability [12]: Structure slide instability (S instability, shear failure) and structure deformation instability (R instability, tension bending failure).

Shear Failure Condition
The shear failure condition of overburden is when the working face is pushed to the cracking end of overburden, the middle part of the strata has not cracked and the shear stress of the remaining section has exceeded the ultimate strength. In this case, when there is a small amount of space under the strata, the overburden will be cut off as a whole, resulting in the overburden being cut off as a whole in the goaf [13,14]. According to the rock beam stress analysis and material mechanics [15], the shear failure expression formula is as follows: where L c -shear distance of rock strata, m: M c -residual thickness of rock beam end after cracking, m; [τ]-Shear strength of rock, MPa; M-actual thickness of rock beam, m; M i -thickness of the following layer above the rock beam, m; γ-average bulk density of rock strata, MPa.

Tensile Bending Failure Condition
The bending failure condition of the overburden on the roof is as follows: With the continuous advance of the stope, the overburden (assumed to be plastic strata) is exposed, and under the action of overlying weak strata and its own gravity, the strata are bent and subsided. The original fixed beam is transformed into a simply supported beam. With the maximum bending moment transferred from the end of the beam to the middle of the beam, the middle part of the beam begins to crack and fall, and finally forms a "pseudo plastic rock beam"; if the height of the lower free space of the rock beam exceeds its allowable settlement value, the rock beam will collapse and the tensile bending failure and instability will occur. The expression formula is: where S 0 -Allowable settlement value of "pseudoplastic beam", m; S n -The free space height that allows the movement of the rock, m; K A -Dilatancy coefficient of caved rock; M-Thickness of caving rock, m; h-Mining height, m.

Theoretical Criterion of Overburden Failure form
When L c > L 0 , the rock strata are pulled away from the middle. When L c = L 0 , the strata are sheared.
Where L c -Rock shear distance, m; L 0 -Limit span at first crack of beam end.

Mechanical Analysis of Longitudinal Propagation of Overburden Fracture in Coal Seam Group Mining
According to the analysis of the strata structure in Xialiyuan coal mine, there are two key strata that play a key role in controlling the whole overburden in the coal seam group working face. The fracture instability characteristics of two key layers, key layer 1 and key layer 2, are analyzed.
In coal seam mining, when the overlying key layer reaches the limit span, it will not collapse immediately, but will form a three hinged arch structure. During the continuous advancement of the mining face, the stability of the three hinged arch structure is mainly controlled by two key blocks. The interaction of the two key blocks presents two failure forms: sliding instability (s) and deformation instability (R).

Failure Judgment of Key Layer "S-R" Judgment of Key Stratum Sliding Instability (S)
After the initial fracture of the key strata, the two rock blocks form an arch equilibrium state ( Figure 4). The horizontal extrusion force T can be considered as the action point at a/2 ( Figure 5), the load borne by the key strata is p and the length and height of the rock block are l and h, respectively. The moment balance is obtained for point O to obtain the balance curve formula without sliding instability: where i-Fracture degree of rock block, ratio of height H to length L of fracture block; tanϕ-The friction coefficient between rocks is 0.3; α-Block rotation angle, M-Mining height; H-Height of rock strata falling; K P -Breaking expansion coefficient, take 1.5; L-Length of fractured rock block.
key strata that play a key role in controlling the whole ove working face. The fracture instability characteristics of tw layer 2, are analyzed. In coal seam mining, when the overlying key layer r collapse immediately, but will form a three hinged arch s advancement of the mining face, the stability of the three controlled by two key blocks. The interaction of the two forms: sliding instability (s) and deformation instability (     According to the relationship between T/P and the rotation angle α, when the fracture degree is constant, the horizontal thrust T between the key blocks increases with the increase of the rotation angle, which will lead to the crushing of the corner blocks and deformation and rotation instability. In order to avoid the deformation and instability of the key blocks, the following conditions should be met: where σ p = T/a -The extrusion stress formed at the bite; η-The extrusion coefficient at the end angle of the key block, ησ c is the extrusion strength, η is 0.3. According to the above formula, the judgment curve formula of rock block without rotary instability is: where n-Ratio of compressive strength to tensile strength of rock; k-Calculated as simply supported beam, k = 0.3; η-Ratio of extruding strength to compressive strength of rock block, η is 0.3. According to the above analysis, the instability form of rock block mainly depends on the fracture degree and rotation angle of the block. Calculated i = 0.793, tanϕ + sinα = 0.646, i > tanϕ + sinα, Therefore, key strata 1 will slip and lose stability.
Therefore, there will be no rotary instability in key strata 1. 2 Judge the instability form of key strata 2, and substitute the relevant parameters of the above judgment formula: tanϕ is 0.3; H is 6.75 m; M: The mining height of coal seam No. 9 is 8 m; h is 11.5 m; k p is 1.5; l is 10.3 m.
Therefore, there will be no rotary instability in key strata 2.

Longitudinal Development Characteristics of Mining-Induced Overburden Fractures in Close Coal Seam Groups
Through the instability of key strata fracture mechanism analysis, the mining face in the working face end composite beams above will produce open fractures, and gradually, from the bottom up, when the rotary rock fracture gradually increases, the main roof structure sliding instability occurs, being exposed for the first time to pressure, formatting a water bursting fracture source "channel", pushing the continuous mining face forward. As a result, the key rocks in the overburden roof break periodically quickly, and the longitudinal fractures in the overburden rock continue to develop upward through the wave base rock. During the mining of coal seam No. 9 in the lower group, the fracture instability of the key strata will directly lead to the fracture of the roof rock in the goaf of coal seam No. 4, resulting in the secondary activation of the upper overburden, resulting in the violent activity of the roof cracks in the partially compacted goaf again, as well as the expansion of the crack channel and the secondary development of overburden cracks. Longitudinal development model of overlying rock fissure in close coal seam mining (see Figure 6). inability 2022, 14, x FOR PEER REVIEW Through the instability of key strata fracture mechanism ana the working face end composite beams above will produce open from the bottom up, when the rotary rock fracture gradually structure sliding instability occurs, being exposed for the first tim a water bursting fracture source "channel", pushing the continuo As a result, the key rocks in the overburden roof break periodica gitudinal fractures in the overburden rock continue to develop up base rock. During the mining of coal seam No.9 in the lower grou of the key strata will directly lead to the fracture of the roof rock No.4, resulting in the secondary activation of the upper overbur lent activity of the roof cracks in the partially compacted goaf aga sion of the crack channel and the secondary development of over dinal development model of overlying rock fissure in close coal s 6).

"Spatial Source"-Spatial Evolution Characteristics of Minin den Fractures in Short Distance Coal Seam
The structure of goaf and overburdening after mining can b complete layered structure, massive structure and broken struct gree of fragmentation and form, and the overall combination m the size of the water accumulation space. In the original strata of mulation space mainly exists in rock micro fractures. After coal m opment of overburden fractures leads to the rapid expansion of space, and rocks in the overburden fracture zone are massive str 7), so the water accumulation space is mainly distributed in the p rock mass in the caving zone, which has been seriously damaged, ture, and the overburden space is fully developed, as is the main w The spatial evolution characteristics of mining-induced overbu

"Spatial Source"-Spatial Evolution Characteristics of Mining-Induced Overburden Fractures in Short Distance Coal Seam
The structure of goaf and overburdening after mining can be divided into relatively complete layered structure, massive structure and broken structure, according to the degree of fragmentation and form, and the overall combination mode is closely related to the size of the water accumulation space. In the original strata of mining, the water accumulation space mainly exists in rock micro fractures. After coal mining, the violent development of overburden fractures leads to the rapid expansion of the water accumulation space, and rocks in the overburden fracture zone are massive structure rocks (see Figure 7), so the water accumulation space is mainly distributed in the pore fracture channel. The rock mass in the caving zone, which has been seriously damaged, presents a broken structure, and the overburden space is fully developed, as is the main water accumulation area. The spatial evolution characteristics of mining-induced overburden fractures are simulated by numerical simulation. The following 40,201 working face area of Xialiyuan Coal Mine is taken as an example. According to the strata distribution and physical and mechanical parameters of coal and rock in Table 1, UDEC numerical simulation software is used to establish a numerical analysis model to simulate the spatial variation characteristics of overburdened fractures in the process of coal seam group mining (see Figures 8 and 9) [16]. The following 40201 working face area of Xialiyuan Coal Mine is taken as an example. According to the strata distribution and physical and mechanical parameters of coal and rock in Table 1, UDEC numerical simulation software is used to establish a numerical analysis model to simulate the spatial variation characteristics of overburdened fractures in the process of coal seam group mining (see Figures 8 and 9) [16].  (c) (d)  In the mining of single coal seams, when the working face advances 20 m, the direct roof collapses and begins to form a water storage space. The overhanging area of the main roof is large, local cracks appear in the lower part of the main roof and fracture instability begins to appear. At this time, after the direct roof collapses, the main roof does not completely collapse under the support of coal pillars on both sides of the goaf, so the goaf is not fully filled, and there is a large space on its left and right. When the aquifer is connected with the goaf, this area becomes the main water accumulation space. When it is excavated to 40 m, initial breaking of main roof occurs. Therefore, it can be inferred that the initial pressure step distance is between 35 and 40 m. In terms of fracture form, due to the large lower space after coal caving, the whole stratification under the main roof first cuts off, the cracks continue to develop in the main roof top space and the water accumulation space keeps increasing. When the working face reaches 160 m, spatial fractures in overburdened rock are basically fully developed, and overburdened rock fractures have spread to the surface, At this time, the middle zone of the goaf fully sinks, and the rock in the caving fracture zone is compressed due to the vertical compressive stress. The water source cannot gather in this area, mainly in the incompletely compacted cracks.
In lower coal seam No. 9, repeated mining and the advancement of the working face to 40 m have caused the overburden fracture to spread to the 4 coal floor; when the working face advances to 50 m, the overburden fracture breaks completely through the 4 coal floor, which occurs as a result of repeated 9 coal mining and has carried on the secondary development of overburdening, which increases in a crack in the channel, increasing the water space, making the risk of water inrush bigger.
In the process of the working face advancing, with the continuous instability and failure of overburdening, the space, size and shape of the fracture network changes accordingly. The general evolution process of the fracture network space is compaction and generation at the same time.

"Solid-Liquid" Coupling Water Inrush Mechanism in Coal Seam Mining
As mentioned above, water inrush in coal seam mining is the interaction of water inrush source, dynamic source, channel source and spatial source. The sudden water source and dynamic source are the original conditions that existed before the mining of the coal seam, that is, there is a strong aquifer. The mining of the coal seam leads to the change of the structure of the upper overburden, and then the channel source and spatial source are generated. Overburden fracture is the channel of water seepage, migration and water storage space. The seepage disaster process of rock strata is the result of the coupling effect of a water-bearing body and consolidated rock mass [17]. According to the deformation and instability process of roof strata, the elastic theory of pore-fissure is referred to [18][19][20][21], and the water inrush process of coal seam mining is divided into three "solid-liquid" coupling action stages:

Original Gap Seepage Stage
Before the coal seam is mined, there is an original fissure in the overburden mass. At this time, the water channel and water accumulation space are insufficient, and the aquifer carries out low-speed seepage under the action of water flow. The solid-liquid coupling equation is as follows: where m:m = 1 represents rock foundation, m = 2 represents rock gap; G is the shear modulus, Pa; λ is the lame constant; u i,jj , u k,ki is displacement, m; α m is Biot constant; p m,i is the fluid pressure, Pa; k is the equivalent permeability value or the average permeability of the overall system, m 2 ; µ is the dynamic viscosity of fluid, m 2 /s; c * is lumped compressibility; ΓIs the fluid exchange rate caused by differential pressure ∆P.

Initial Water Discharge Stage of Mining-Induced Fracture
During coal seam mining, under the influence of mining, the cracks of overburden rock mass expand and gradually form connectivity; under the action of the water pressure of the aquifer, the water from an aquifer can move down into the goaf through the cracks in the overburden, which is manifested as seepage in the roof of the working face. The solid-liquid coupling equation satisfied by the seepage of mining fractured rock mass in this process is: where k m is permeability of m phase, m 2 .

Water Inrush Disaster Stage of Fractured Rock Mass
With the advance of coal seam mining, the overburden fractures will further expand until the rock mass structure is damaged and collapses. At this time, the aquifer water has a sufficient migration channel and water storage space, and breaks into the goaf in a large area in the form of water inrush at the working face. The solid-liquid coupling equation satisfied by the seepage of broken rock mass in this process is: where m = 1, 2 and 3 represent pores, fissures and fractures, respectively; k 1 , k 2 and k 3 are pores, fractures and fracture permeability, m 2 ; k 13 is the average permeability of rock foundation and fracture, m 2 ; Γ ij is the fluid exchange rate between phase I and phase J, and it is assumed that there is interstitial flow caused by internal pressure difference between the two phases; others conform to the same meaning.

Forecast of Water Inrush and Inflow in Working Face
The accurate prediction of water inflow plays an important role in preventing and controlling the occurrence and control of water disasters and can guide mines to formulate prevention and control measures. According to the characteristics of water inrush, the inflow formula of a confined aquifer with a constant depth is used to predict the mine water inflow [22][23][24], and the stable value of inflow and its variation characteristics with time are predicted.

Mathematical Model
In the early stage of coal seam mining, because the drainage time was not long enough, the falling funnel kept expanding and did not reach the stable flow state. However, the water level has been reduced to the roof of the coal seam, so the mathematical model of the water flow movement should be a definite solution of the variable flow with constant depth. Its mathematical model is This mathematical model is obtained by Laplace transform and Darcy's law, and is: where s w is the drawdown in the well; λ = Tt r 2 w S is dimensionless time; G (λ) is the flow function of constant drawdown well flow in confined aquifer without overflow recharge.

Forecast of Water Inflow
In order to obtain the process Q(t) of water inflow with time, different time points t 1 , t 2 , . . . and t n take in (6.17); the corresponding values at different time points can be obtained as G (λ1), G (λ2), . . . and G (λn) and Q 1 , Q 2 , . . . and Q n . The value of G (λ) is shown in Table 2. The calculation parameters of water inflow of coal seam No. 4 are shown in Table 3, and the value of time t is shown in Table 4. Bring the parameter values in Table 3 into formula 11 and take values according to Table 2 to obtain (Table 4 and Figure 10).  Table 4. Bring the parameter values in Table 3 into formula 11 and take values accord Table 2 to obtain (Table 4 and Figure 10).  It can be seen from Figure 10 that after about 300 days of drainage, the water i of coal 4 is basically stable, and the flow is 27.8 m 3 /h. Water Inflow of Coal Seam No.9 The calculation parameters of the water inflow of coal seam No.9 are shown in 2, and the value of time t is shown in Table 5.  It can be seen from Figure 10 that after about 300 days of drainage, the water inflow of coal 4 is basically stable, and the flow is 27.8 m 3 /h. Water Inflow of Coal Seam No. 9 The calculation parameters of the water inflow of coal seam No. 9 are shown in Table 2, and the value of time t is shown in Table 5. Bring the parameter values in Table 5 into formula 11 and take values according to  Table 2 to obtain (Table 6 and Figure 11). Bring the parameter values in Table 5 into formula 11 and take values accord  Table 2 to obtain (Table 6 and Figure 11).  It can be seen from Figure 11 that after about 500 days of drainage, the water of coal 9 is basically stable, and the flow is 105.10 m 3 /h.

Water Inrush Prevention and Control Measures
(1) The water-rich area of the working face should be appropriately accelerated.
one hand, the processing speed can shorten the overburden dynamic deform process time; the better the overburden compaction effect is, the more the com sion of "spatial source". On the other hand, the accelerating speed can short development time of the overlying rock fissures and can quickly close, prev the development of a "channel source". (2) Drainage of water-rich aquifers in the upper working face before mining, as w drilling and drainage in advance and drainage and pressure reduction to redu "water inrush source" and eliminate the "dynamic source". It can be seen from Figure 11 that after about 500 days of drainage, the water inflow of coal 9 is basically stable, and the flow is 105.10 m 3 /h.

Water Inrush Prevention and Control Measures
(1) The water-rich area of the working face should be appropriately accelerated. On the one hand, the processing speed can shorten the overburden dynamic deformation process time; the better the overburden compaction effect is, the more the compression of "spatial source". On the other hand, the accelerating speed can shorten the development time of the overlying rock fissures and can quickly close, preventing the development of a "channel source". (2) Drainage of water-rich aquifers in the upper working face before mining, as well as drilling and drainage in advance and drainage and pressure reduction to reduce the "water inrush source" and eliminate the "dynamic source". (3) Thickness limited mining or local filling in the water-rich area of the working face.
This method can restrain the development of longitudinal fractures in working face.

Conclusions
(1) Water inrush from coal seam mining consists of "four source elements", which are water inrush source, dynamic source, channel source and spatial source, The water inrush source and dynamic source are the original conditions existing before coal seam mining; the mining of the coal seam changes the upper overburden structure and then produces the channel source and spatial source; when the "four source elements" constitute a complete relationship, water inrush will occur in the working face. (2) Using the key strata theory, this paper carries out a mechanical analysis of the longitudinal expansion of overburden fractures during mining of short distance coal seam groups and grasps the development and evolution characteristics of a "channel source". Through mechanical model analysis, the two key strata that play a key role in controlling the whole overburden in the mining of coal seam groups in the study area have the form of shear failure (sliding instability) during mining. (3) The UDEC numerical simulation method is used to simulate the spatial evolution characteristics of mining-induced overburden fractures. The results show that in the process of the working face advancing, with the continuous instability and failure of overburden, the spatial size and shape of fracture network are also changing, and the general evolution process is as follows: generation-compaction and generationsynchronous action. (4) The seepage catastrophe process of rock stratum is the result of the intercoupling of a water bearing body and a consolidated rock mass. According to the deformation and instability process of roof strata and the theory of pore-fracture elasticity, the water inrush process of coal seam mining is divided into three "solid-liquid" coupling action stages: the original gap seepage stage, the initial water discharge stage of a mining-induced fracture and the water inrush disaster stage of a fractured rock mass.