# Research on the Extraction Technology of Gas-Integrated Working-Face in Close Coal Seam Mining in Permafrost Area

Authors: Chao Wu; Rongbin Deng
Journal: Geofluids
 (2022)

---

## Abstract

In order to understand the mechanical deformation characteristics of the permafrost layer and the stress-fracture-gas seepage law of the superimposed mining of the proximity coal seam group, a coupled gas-solid model of the proximity coal seam in the permafrost area was constructed based on the theoretical model of small deflection of thin plate and the deformation control equation, gas transport equation, porosity evolution equation and permeability evolution equation of the coal rock body, and the analysis of the dependent project. The analysis was carried out at the NMT coal mine. The results show that the multiyear permafrost seam in the NMT coal mine has a tendency to move to the left and will produce extensive pull damage near the right end. The initial and periodic collapse steps are 192.87 m and 43.24 m, respectively; the pore and fissure structural characteristics of the NMT coal seam are not conducive to gas transport, and the stress and fissure fields have a nonlinear increasing effect. The study also proposed the “protected seam mining and bottom extraction lane directional drilling group comining technology” to increase the daily production of working face by 25.8%. The results of the study can be used as a reference for similar projects.


---

## Body

## 1. Introduction

The complex permafrost geological environment has caused a large number of damage problems in the permafrost layer due to freezing, which seriously threatens the safe and efficient mining of coal fields and makes the control work more difficult than ever [1, 2]. Xu et al. based on a large number of experimental studies and the summary of previous research results comprehensively analyzed and studied the influence of the main factors [3]. The study through the study of the weak interlayer in the Dajiang River of the Gezhouba Project puts forward some types of microstructures of the weak interlayer [4]. The relationship between ground stress and the physical and mechanical parameters of the weak layer zone was studied in the study, taking the weak interlayer at a hydroelectric dam site on the Jinsha River as an example, and it was concluded that ground stress is the most active factor controlling the engineering properties of the weak layer zone; the effect of freeze-thaw cycles on the shear strength properties of lime powder soils was discussed in the study [5, 6]. The relationship between the physical and mechanical parameters of the weak layer zone and the physical and mechanical parameters of the weak layer zone was studied in the study [7]. The mechanical properties of soils change before and after freeze-thawing, and the main reason for the change in soil mechanical properties due to repeated freeze-thawing is the change in soil structure, i.e., the destruction of the linkage between soil particles and the rearrangement of soil particles [8–12]. In summary, little research has been reported on the mechanical properties of coal mines, especially in high-altitude permafrost. Therefore, it is of practical importance to research the mechanical properties of strata in high-altitude permafrost areas. Close coal seam groups are widely distributed in China, with small seam spacing and large interaction between seams, and coal seam gas is always a threat to mine safety during mining. As the mining depth continues to increase, the gas content and pressure of the coal seam increases, making it more difficult to manage the gas. In the area of close coal seam group: the study explored the stress-fracture evolution characteristics of the close coal seam group under superposition mining conditions and realized the refined simulation of the distribution characteristics of the mining fracture zone [13, 14]. Qi et al. studied the mining influence range of the first mining layer under protective layer mining and analyzed the feasibility of the protective layer mining technology of the high gas coal seam group [15]. Ma et al. found that the overburden collapse characteristics under superposition mining conditions [16]. The scholars studied the spatial and temporal evolution characteristics of the three-dimensional stress field and the fracture mechanics behavior of the overburdened fracture zone under repeated mining of coal seam clusters [17, 18]. In the area of coal and gas comining: literature initially puts forward the new concept of “green mining” [19, 20]. Literature studied the fracture mechanism and fractal characteristics of the overlying rock layer and revealed the evolution law of mining fractures [21–24]. Papers changed the traditional “wind drainage” gas management measures and adopted pressure relief mining in the first mining layer to achieve a new model of efficient “extraction” of gas and formed a coal and gas comining technology under protective layer mining [25]. Papers established the theory of efficient coal and gas comining under coal pillarless mining and achieved a quantitative description of coal seam gas [26–29]. A large-aperture kilometer directional drilling rig was used to efficiently extract gas from the fissure zone of the roof, and the mechanism of coal and gas comining in kilometer directional drilling in high-gas coal seams was proposed [30, 31]. The nonlinear characteristics and failure characteristics of coal energy evolution under different gas pressures are studied [32].The above research results have effectively solved the technical problems that exist in the mining process of close coal seam clusters and constitute the basic content of coal and gas comining technology in China, but there is less research in the direction of the coupling mechanism of the stress field and gas pressure under the mining conditions of close coal seam clusters. Based on this, this paper takes the NMT mine as the research background and investigates the stress field and gas pressure crystallization under the superimposed mining conditions in the mine, to provide theoretical support for the coal and gas comining technology under the close coal seam group conditions.


## 2. Theoretical Model

Assumptions: coal is a porous medium consisting of pores-fissures with nonuniformity, anisotropy, and discontinuity, and also considering the characteristics of large depth of gas deposit, coal is taken as the object of study to analyze the gas extraction characteristics of proximity coal seam mining, and the following assumptions are made [33].
(1)
It is assumed that the coal rock mass is an isotropic homogeneous linear elastic material(2)
The coal seam gas adsorption model conforms to the Langmuir equation and the gas flow satisfies Darcy’s law(3)
The gas is at a constant temperature and the transport and diffusion of a gas is an isothermal processCoupling model: based on the above assumptions, a gas-solid coupling model is constructed. It is shown as Figure1.Figure 1
Schematic diagram of the coupling model.


### 2.1. Deformation and Mechanical Characteristics of Permafrost Coal (Rock) Seams

A permafrost environmental coal (rock) seam, subject to temperature loading, can be decomposed into two components, longitudinal load stress acting between the layers and a transverse load stress perpendicular to the midplane. For the longitudinal load stresses, they can be considered to be analyzed uniformly along the thickness direction, and therefore, the stresses, strains, and displacements caused by them can be calculated by solving for the plane stress problem. As shown in Figure2, the transverse loads will cause the sheet to bend, and the stresses, strains, and displacements caused can be calculated as a sheet bending problem. In the thin plate bending, the plane bent into the surface, known as the elastic surface of the thin plate, and the point in the surface perpendicular to the plane, along the z-direction of displacement, is called the deflection ω.Figure 2
Temperature loading in the permafrost zone.Because the thickness of the coal (rock) seam is much less than the planar dimension, the Kirchhoff hypothesis is invoked to establish a theoretical model for small deflections in thin slabs.(1)
There is no deformation in the midplane of the plate, i.e., the midplane bends to a neutral surface during bending. It follows thatεz, γzx, and γzy are negligible. This leads to the following:
(1)ω=ωx,y,where ε is a positive strain; γ is the shear strain; and ω is the deflection(2)
The deformation caused by the strain componentεz can be disregarded. By this assumption and in combination with assumption Equation (1), the physical Equation (2) for the small deflection of a thin plate can be obtained
(2)εx=1Eσx−υσy,εy=1Eσy−υσx,γxy=21+υEτxy,where E is the modulus of elasticity; υ is the Poisson’s ratio(3)
The points in the plane have no displacement parallel to the midplane. That is, whenz=0, u=0, and v=0, i.e., εx=0, εy=0, and γxy=0, then the differential equation for the elastic surface of a thin plate with a small deflection problem Equation (3), i.e., the basic equation of the problem, can be obtained
(3)Eh3121−υ2∇4ω=q,where q is the transverse load on the thin plate; h is the thickness of the thin plate; ∇ is the Laplace operator, ∇4ω=∂4ω/∂x4+2∂4ω/∂x2∂y4+∂4ω/∂y4.This gives the relationship between the cross-sectional internal force of the thin plate and the deflection of the thin plate, and the parametric elimination of the deflection gives Equations (4) and (5). From this, each stress extreme and its corresponding location can be found. In solving the multiyear permafrost deformation and damage problem, the load-induced bending and torsional stresses are numerically larger and are the main stresses; transverse shear stresses are numerically relatively small and are secondary stresses; extrusion stresses are numerically smaller and are more secondary stresses. Therefore, when solving thin-slab problems with multiyear permafrost layers, the focus is on solving for the major internal forces, i.e., bending moments and torsional moments.At the start of the coal seam excavation, before the initial break of the multiyear permafrost seam with the coal seam group, the mining mechanics analysis model is a thin plate surrounded by solid coal rock-solid support; the mechanical calculation model is shown in Figure3. In the figure, 2e is the length of the permafrost layer, and 2f is the width of the permafrost layer. The damage forces in the vertical direction are mainly caused by their own gravity G. The uniform load qt is the perennial permafrost capacity.
(4)Mx=−D∂2ω∂x2+v∂2ω∂y2,My=−D∂2ω∂y2+v∂2ω∂x2,Mxy=Mxy=−D1−v∂2ω∂x∂y,Qx=−D∂∂x∇2ω,Qy=−D∂∂y∇2ω,(5)σx=12Mxh3z,σy=12Myh3z,σz=−2q12−zh21+zh,τxy=τyx=12Mxyh3z,τxz=6Qxh3h24−z2,τyz=6Qyh3h24−z2.Figure 3
Mechanical model of initial permafrost break.The boundary conditions at this point are∂ω/∂xx±cωx=±e=0=0∂ω/∂xx+±ωx=±f=0=0.Using an approximate solution, assume that the deflection function of the perennial permafrost “slab cover” is(6)ω=Ax2−e22y2−f22.Taking the partial derivatives ofx and y, respectively, yields
(7)∂ω∂x=4Axx2−e22y2−f22,∂ω∂y=4Ayx2−e22y2−f22.Based on the principle of minimum potential energy for bending of thin plates, the Galyagin method was used to determine the coefficientA, A=49⋅qt/128⋅7e4+4e2f2+7f4.The equation for the deflection surface of a four-sided solidly supported perennial permafrost “slab cover” is(8)ω=49q11287e4+4e2f2+7f4x2−e22y2−f22.From this, the parametersMx, My, Mxy, Qx, Qy, etc. can be calculated, where Mx and My are
(9)Mx=−7qt32e4+4/7e2f2+f43x2−e2y2−f22+υx2−e23y2−f22,Mx=−7qt32e4+4/7e2f2+f43y2−e2x2−f22+υy2−e23x2−f22.From Equation (9), the perennial permafrost layer “plate cover” Mx,My bending moment maximum can be determined, when x=±e on, y=0, Mxmax for Equation (10); y=±f on, x=0, Mymax for Equation (11).
(10)Mxmax=−7qt16e4+4/7e2f2+f4e2f4,(11)Mymax=−7qt16e4+4/7e2f2+f4e4f2.Equation (10) and Equation (11) show that the bending moment reaches its maximum value at the midpoint of the long solid support and at the midpoint of the short solid support in the permafrost layer “slab cover.” As the stresses in the section are proportional to the bending moment, it can be concluded that the positive stresses at the midpoint of the long solid support and at the midpoint of the short solid support also reach their maximum values at the section z=±Ht/2, i.e., the positive stresses on the upper and lower sides of the plate reach their maximum values in
(12)σxmax=42qx16e4+4/7e2f2+f4h2e2f4,σymax=42q16e4+4/7e2f2+f4h2e4f2.The peak point of the positive stress caused by the bending moment on both the long and short sides is located at the midpoint of each side of the permafrost “slab cover,” and the bending moment at the midpoint of the long side is much greater than that at the midpoint of the short side, which means that the four-sided fixed permafrost “slab cover” should be at the middle of the upper surface of the long side. This means that the permafrost “cover” is most likely to break in the middle of the working face.The permafrost layer is a brittle material with low tensile strength, so the bending deformation before the permafrost layer breaks is small. When the longitudinal load reaches the critical load for the flexure of the thin plate, if it is only the effect of the longitudinal load, a slight increase in the load at this time will make the displacement and internal force increase a lot, while in fact, the presence of the transverse load makes the displacement and internal force of the plate increase a lot, which eventually leads to the breakage of the permafrost layer. The rupture will be a tensile rupture dominated by buckling instability. When the maximum tensile stress at a certain point in the permafrost layer reaches the permafrost tensile strength, the permafrost will fracture at that point, and the breaking conditions will be(13)σmax≥σt,where σmax is the maximum principal stress in the permafrost; σt is the tensile strength limit of the permafrost.As a result of the above analysis, the permafrost “slab cover” is prone to tensile damage on both sides of the inclined direction. In the vertical direction, the permafrost “slab cover” is subjected to its own gravity, and due to the movement deformation off-layer cracks between the protective rock beam and the permafrost “slab cover,” the permafrost “slab cover” will deform following plate mechanics. In addition, after the coal seam has been mined, the perennial permafrost layer will have a significant overhang in one direction, at which point the perennial permafrost layer “plate cover” is highly susceptible to significant plastic deformation on both sides. The deformation results in a state where the left and right sides become simply supported and the front and rear sides are solidly supported, the mechanical model is shown in Figure4.Figure 4
Mechanical model of suspended fracture.The boundary conditions at this point are shown as follows.(14)ωx=0x−2e=0∂2ω∂x2x=0x−2e=0,ωx=±f=0∂ω∂yx=±f=0.The deflection function is obtained by Levy’s solution as(15)ω=64qte4π5D∑m=1,3,5,⋯∞1m51−amcham+shamchmπy/achamamcham+sham−amsh2am+shamchamamcham+sham−amsh2ammπy2eshmπyasinmπx2e.Of which,am=mπf/2e.This level converges quickly, and a good approximate solution can be obtained by takingm=1. Now, let the permafrost “plate cover” be a square plate with side length 2e in the inclined direction and length 2e in the strike direction, and the maximum deflection can be found at the center of the plate. Substituting x=e, y=0, and a=π/2, into the above equation gives
(16)ωmax=0.0304q1e4D=0.0304121−U2qte4EHt3.It can be seen that the maximum tensile stress occurs atz=−Ht/2 at the lower end of the middle of the permafrost “slab cover,” so that after the deformation of the permafrost “slab cover” on both sides in the inclined direction, the limit of the collapse distance in the inclined direction is reached, and the mechanical deformation will occur in the middle of the coal seam in the inclined direction. The mechanical deformation occurs in the middle of the seam.


### 2.2. Deformation Control Equations for Coal Rock Masses

The deformation control equation of a gas-bearing coal body consists of three parts: the stress balance equation, the geometric equation, and the intrinsic equation. The stress balance equation for a gas-bearing coal body is(17)σij,j+Fi=0,where σij,j is the first-order partial derivative of the total stress component σij of the coal rock body for the coordinates (i,j in the principal direction), N/m3; Fi is the volume force component, N/m3.The geometric equation for the deformation of the coal rock body is:(18)εij=12ui,j+uj,i,where εij is the line strain component; ui,j,uj,i are the first-order partial derivatives of the displacement components ui,uj, respectively, with respect to the coordinates.From the generalized Hooke’s law, the intrinsic equation of the coal rock body is obtained as shown in Equation (19), where G is the shear modulus of the coal rock body, MPa; Ev is the volume strain of the coal rock body; δij is the Kronecker variable (1 for i=j, 0 for i≠j); αp is the Biot coefficient, and p is the gas pressure, MPa.
(19)σij=2Gεij+2Gv1−2vεVδij−αppδij.Combining Equations (17) to (19), the deformation control equation for coal rock masses considering pore pressure can be expressed in terms of displacement.
(20)Gui,jj+G1−2vuj,ji−αpp+Fi=0,where ui,jj,uj,ji are the second-order partial derivatives of the displacement components ui,uj for the coordinates, respectively.


### 2.3. Gas Transport Control Equations

The Langmuir isothermal sorption model is shown in(21)VEp=VLppL+p,where VEp is the adsorption volume, g/m3; VL is the Langmuir volume g/m3; pL is the Langmuir pressure, MPa, and the Langmuir pressure is the pressure at which the adsorption volume reaches 50% of the ultimate adsorption volume.The fracture system of the coal body provides a channel for the free gas to deposit and seep through, and the mass balance equation of the gas in the fracture system of the coal body per unit volume is(22)∂mf∂t+∇ρgq=Q,where mf is the mass of free gas per unit volume of coal, kg/m3; t is time, s; ρg is the gas density, kg/m3; q is the Darcy velocity ignoring the effect of gravity, m/s; q=−k▽p/μ, k is the coal permeability, m2; μ is the gas dynamic viscosity, Pa∙s; Q is the mass exchange rate between the matrix system and the fracture system per unit volume of coal, kg/(m3∙s).Under isothermal conditions, the free gas content is only related to the gas pressure. According to the gas equation of state, the mass of free gas in a unit volume of gas-bearing coal can be obtained as(23)mf=φρg=φρnpnp,whereφ is the porosity of the coal body, %; ρn is the density of gas at standard atmospheric pressure, kg/m3; pn is the standard atmospheric pressure, MPa.The diffusion of gas in the coal body is mainly in the gas phase. According to the law of conservation of mass, the rate of change of gas mass per unit volume of coal matrix system with time is equal to the rate of mass exchange per unit volume of coal matrix system with the fracture system.(24)∂mγ∂t=−Q,where mr is the mass of adsorbed gas per unit volume of coal, kg/m3. Considering the influence of moisture, ash, and combustible percentage in coal on gas adsorption, the mass of adsorbed gas per unit volume of coal can be obtained from the modified Langmuir equation, as shown in
(25)mγ=abp1+bpρgρc100−A−M1001+0.31M,where a is the ultimate gas adsorption capacity of the coal body, m 3/kg; b is the adsorption constant of the coal body, MPa-l; ρc is the density of the coal body, kg/m3; A is the ash content of the coal body, %: M is the moisture content of the coal body, %.LetR=ρgρc100−A−M/1001+0.31M, then Equation (26) can be obtained
(26)∂φρnp/pn∂t=∇kρnpμpn∇p=−∂abpR/1+bp∂t.


### 2.4. Porosity and Permeability Evolution Equations

Coal porosity and permeability are the key coupling terms in the coupling of deformation field and gas transport field. By the definition of porosity, the coal porosity can be obtained by deducing that(27)φ=VpVB=1−1−φ01+εv1+ΔVsΔVs0,where Vp is the pore volume, m3; VB is the coal body appearance volume, m3; φ0 is the initial porosity, %; ΔVs is the change in coal body skeleton volume, m3; ΔVs0 is the initial coal body skeleton volume change, m3.The volume deformationΔVs/ΔVs0 includes the deformation ΔVsp/ΔVs0 caused by changes in gas pressure, the deformation ΔVsf/ΔVΔVs0 caused by the expansion of gas adsorbed by coal particles, and the deformation ΔVst/ΔVΔVs0 caused by the thermal expansion of coal skeletal particles.
(28)ΔVsΔVs0=ΔVsp+ΔVsf+ΔVstΔVs0=31−2vEΔp+εLpp+pL1−φ0+αΔT,where E is the modulus of elasticity of the coal body, GPa; Δp is the change in gas pressure, MPa; εL is the Langmuir volume strain constant; PL is the Langmuir pressure constant, MPa; α is the volume expansion coefficient of the coal body, K-1; ∆T is the absolute temperature change, K.Considering the effects of gas pressure, coal rock volume strain, and temperature on porosity, the dynamic expression for the porosity of gas-bearing coal under isothermal conditions is derived as(29)φ=1−1−φ01+εV1−31−2vEΔp+εLpp+pL1−φ0,where φ0 is the porosity of the coal body under the initial conditions; εV is the volume strain of the coal body under the corresponding gas pressure conditions.Based on the Kozeny-Carman equation and defined by porosity and volumetric strain, this gives(30)kk0=ϕϕ03=11+εV1+β2Vp0+ΔVpVp03,where k0 is the initial coal body permeability, m2; β is the ratio of the change in pore internal surface area to the initial internal surface area of the pore, which is much smaller than 1 and therefore negligible; Vp0 is the initial pore volume of the coal body, m, and ΔVp is the change in the pore volume of the coal body, expressed by the following equation:
(31)ΔVp=εvVB0−εLp1−φ0PL+p−31−2vEΔpVs0,where VB0 is the initial appearance volume of the coal rock body, m3; VS0 is the initial skeletal volume of the coal rock body, VS0=VB0−VP0, m3.Letγ=εLp/1−φ0PL+p−31−2v/EΔp simplify the above equation to give
(32)kk0=11+εv1+γ+εV−γφ03.The resulting equation for the dynamic evolution of coal body permeability is(33)k=k01φ0−1−φ0φ01+εV1−31+2VEΔp+εLpp+PL1+φ03.


### 2.5. Gas-Solid Coupling Model

Gas movement with the gas flow in confined space is essentially a single-phase multicomponent diffusion problem without chemical reaction. Its motion law satisfies the basic equations of fluid mechanics described by general Navier-Stokes equations, including continuity equation, momentum conservation equation, energy conservation equation, and component mass conservation equation. The control equations of coal and rock deformation and gas migration are unified to form a gas-solid coupling model. The couplings (8), (15), (19), (25), (28), and (32) together form a mathematical model for gas-solid coupling in the permafrost region close to the coal seam, see
(34)Gui,jj+G1−2vuj,ji−αpp+Fi=0,∂φρnp/pn∂t−∇kρnpμpn∇p=−∂abpR/1+bp∂t,φ=1−1−φ01+εV1−31−2vEΔp+εLpp+PL1−φ0,k=k01φ0−1−φ01+εV1−31−2vEΔp+εLpp+PL1−φ03.


## 2.1. Deformation and Mechanical Characteristics of Permafrost Coal (Rock) Seams

A permafrost environmental coal (rock) seam, subject to temperature loading, can be decomposed into two components, longitudinal load stress acting between the layers and a transverse load stress perpendicular to the midplane. For the longitudinal load stresses, they can be considered to be analyzed uniformly along the thickness direction, and therefore, the stresses, strains, and displacements caused by them can be calculated by solving for the plane stress problem. As shown in Figure2, the transverse loads will cause the sheet to bend, and the stresses, strains, and displacements caused can be calculated as a sheet bending problem. In the thin plate bending, the plane bent into the surface, known as the elastic surface of the thin plate, and the point in the surface perpendicular to the plane, along the z-direction of displacement, is called the deflection ω.Figure 2
Temperature loading in the permafrost zone.Because the thickness of the coal (rock) seam is much less than the planar dimension, the Kirchhoff hypothesis is invoked to establish a theoretical model for small deflections in thin slabs.(1)
There is no deformation in the midplane of the plate, i.e., the midplane bends to a neutral surface during bending. It follows thatεz, γzx, and γzy are negligible. This leads to the following:
(1)ω=ωx,y,where ε is a positive strain; γ is the shear strain; and ω is the deflection(2)
The deformation caused by the strain componentεz can be disregarded. By this assumption and in combination with assumption Equation (1), the physical Equation (2) for the small deflection of a thin plate can be obtained
(2)εx=1Eσx−υσy,εy=1Eσy−υσx,γxy=21+υEτxy,where E is the modulus of elasticity; υ is the Poisson’s ratio(3)
The points in the plane have no displacement parallel to the midplane. That is, whenz=0, u=0, and v=0, i.e., εx=0, εy=0, and γxy=0, then the differential equation for the elastic surface of a thin plate with a small deflection problem Equation (3), i.e., the basic equation of the problem, can be obtained
(3)Eh3121−υ2∇4ω=q,where q is the transverse load on the thin plate; h is the thickness of the thin plate; ∇ is the Laplace operator, ∇4ω=∂4ω/∂x4+2∂4ω/∂x2∂y4+∂4ω/∂y4.This gives the relationship between the cross-sectional internal force of the thin plate and the deflection of the thin plate, and the parametric elimination of the deflection gives Equations (4) and (5). From this, each stress extreme and its corresponding location can be found. In solving the multiyear permafrost deformation and damage problem, the load-induced bending and torsional stresses are numerically larger and are the main stresses; transverse shear stresses are numerically relatively small and are secondary stresses; extrusion stresses are numerically smaller and are more secondary stresses. Therefore, when solving thin-slab problems with multiyear permafrost layers, the focus is on solving for the major internal forces, i.e., bending moments and torsional moments.At the start of the coal seam excavation, before the initial break of the multiyear permafrost seam with the coal seam group, the mining mechanics analysis model is a thin plate surrounded by solid coal rock-solid support; the mechanical calculation model is shown in Figure3. In the figure, 2e is the length of the permafrost layer, and 2f is the width of the permafrost layer. The damage forces in the vertical direction are mainly caused by their own gravity G. The uniform load qt is the perennial permafrost capacity.
(4)Mx=−D∂2ω∂x2+v∂2ω∂y2,My=−D∂2ω∂y2+v∂2ω∂x2,Mxy=Mxy=−D1−v∂2ω∂x∂y,Qx=−D∂∂x∇2ω,Qy=−D∂∂y∇2ω,(5)σx=12Mxh3z,σy=12Myh3z,σz=−2q12−zh21+zh,τxy=τyx=12Mxyh3z,τxz=6Qxh3h24−z2,τyz=6Qyh3h24−z2.Figure 3
Mechanical model of initial permafrost break.The boundary conditions at this point are∂ω/∂xx±cωx=±e=0=0∂ω/∂xx+±ωx=±f=0=0.Using an approximate solution, assume that the deflection function of the perennial permafrost “slab cover” is(6)ω=Ax2−e22y2−f22.Taking the partial derivatives ofx and y, respectively, yields
(7)∂ω∂x=4Axx2−e22y2−f22,∂ω∂y=4Ayx2−e22y2−f22.Based on the principle of minimum potential energy for bending of thin plates, the Galyagin method was used to determine the coefficientA, A=49⋅qt/128⋅7e4+4e2f2+7f4.The equation for the deflection surface of a four-sided solidly supported perennial permafrost “slab cover” is(8)ω=49q11287e4+4e2f2+7f4x2−e22y2−f22.From this, the parametersMx, My, Mxy, Qx, Qy, etc. can be calculated, where Mx and My are
(9)Mx=−7qt32e4+4/7e2f2+f43x2−e2y2−f22+υx2−e23y2−f22,Mx=−7qt32e4+4/7e2f2+f43y2−e2x2−f22+υy2−e23x2−f22.From Equation (9), the perennial permafrost layer “plate cover” Mx,My bending moment maximum can be determined, when x=±e on, y=0, Mxmax for Equation (10); y=±f on, x=0, Mymax for Equation (11).
(10)Mxmax=−7qt16e4+4/7e2f2+f4e2f4,(11)Mymax=−7qt16e4+4/7e2f2+f4e4f2.Equation (10) and Equation (11) show that the bending moment reaches its maximum value at the midpoint of the long solid support and at the midpoint of the short solid support in the permafrost layer “slab cover.” As the stresses in the section are proportional to the bending moment, it can be concluded that the positive stresses at the midpoint of the long solid support and at the midpoint of the short solid support also reach their maximum values at the section z=±Ht/2, i.e., the positive stresses on the upper and lower sides of the plate reach their maximum values in
(12)σxmax=42qx16e4+4/7e2f2+f4h2e2f4,σymax=42q16e4+4/7e2f2+f4h2e4f2.The peak point of the positive stress caused by the bending moment on both the long and short sides is located at the midpoint of each side of the permafrost “slab cover,” and the bending moment at the midpoint of the long side is much greater than that at the midpoint of the short side, which means that the four-sided fixed permafrost “slab cover” should be at the middle of the upper surface of the long side. This means that the permafrost “cover” is most likely to break in the middle of the working face.The permafrost layer is a brittle material with low tensile strength, so the bending deformation before the permafrost layer breaks is small. When the longitudinal load reaches the critical load for the flexure of the thin plate, if it is only the effect of the longitudinal load, a slight increase in the load at this time will make the displacement and internal force increase a lot, while in fact, the presence of the transverse load makes the displacement and internal force of the plate increase a lot, which eventually leads to the breakage of the permafrost layer. The rupture will be a tensile rupture dominated by buckling instability. When the maximum tensile stress at a certain point in the permafrost layer reaches the permafrost tensile strength, the permafrost will fracture at that point, and the breaking conditions will be(13)σmax≥σt,where σmax is the maximum principal stress in the permafrost; σt is the tensile strength limit of the permafrost.As a result of the above analysis, the permafrost “slab cover” is prone to tensile damage on both sides of the inclined direction. In the vertical direction, the permafrost “slab cover” is subjected to its own gravity, and due to the movement deformation off-layer cracks between the protective rock beam and the permafrost “slab cover,” the permafrost “slab cover” will deform following plate mechanics. In addition, after the coal seam has been mined, the perennial permafrost layer will have a significant overhang in one direction, at which point the perennial permafrost layer “plate cover” is highly susceptible to significant plastic deformation on both sides. The deformation results in a state where the left and right sides become simply supported and the front and rear sides are solidly supported, the mechanical model is shown in Figure4.Figure 4
Mechanical model of suspended fracture.The boundary conditions at this point are shown as follows.(14)ωx=0x−2e=0∂2ω∂x2x=0x−2e=0,ωx=±f=0∂ω∂yx=±f=0.The deflection function is obtained by Levy’s solution as(15)ω=64qte4π5D∑m=1,3,5,⋯∞1m51−amcham+shamchmπy/achamamcham+sham−amsh2am+shamchamamcham+sham−amsh2ammπy2eshmπyasinmπx2e.Of which,am=mπf/2e.This level converges quickly, and a good approximate solution can be obtained by takingm=1. Now, let the permafrost “plate cover” be a square plate with side length 2e in the inclined direction and length 2e in the strike direction, and the maximum deflection can be found at the center of the plate. Substituting x=e, y=0, and a=π/2, into the above equation gives
(16)ωmax=0.0304q1e4D=0.0304121−U2qte4EHt3.It can be seen that the maximum tensile stress occurs atz=−Ht/2 at the lower end of the middle of the permafrost “slab cover,” so that after the deformation of the permafrost “slab cover” on both sides in the inclined direction, the limit of the collapse distance in the inclined direction is reached, and the mechanical deformation will occur in the middle of the coal seam in the inclined direction. The mechanical deformation occurs in the middle of the seam.


## 2.2. Deformation Control Equations for Coal Rock Masses

The deformation control equation of a gas-bearing coal body consists of three parts: the stress balance equation, the geometric equation, and the intrinsic equation. The stress balance equation for a gas-bearing coal body is(17)σij,j+Fi=0,where σij,j is the first-order partial derivative of the total stress component σij of the coal rock body for the coordinates (i,j in the principal direction), N/m3; Fi is the volume force component, N/m3.The geometric equation for the deformation of the coal rock body is:(18)εij=12ui,j+uj,i,where εij is the line strain component; ui,j,uj,i are the first-order partial derivatives of the displacement components ui,uj, respectively, with respect to the coordinates.From the generalized Hooke’s law, the intrinsic equation of the coal rock body is obtained as shown in Equation (19), where G is the shear modulus of the coal rock body, MPa; Ev is the volume strain of the coal rock body; δij is the Kronecker variable (1 for i=j, 0 for i≠j); αp is the Biot coefficient, and p is the gas pressure, MPa.
(19)σij=2Gεij+2Gv1−2vεVδij−αppδij.Combining Equations (17) to (19), the deformation control equation for coal rock masses considering pore pressure can be expressed in terms of displacement.
(20)Gui,jj+G1−2vuj,ji−αpp+Fi=0,where ui,jj,uj,ji are the second-order partial derivatives of the displacement components ui,uj for the coordinates, respectively.


## 2.3. Gas Transport Control Equations

The Langmuir isothermal sorption model is shown in(21)VEp=VLppL+p,where VEp is the adsorption volume, g/m3; VL is the Langmuir volume g/m3; pL is the Langmuir pressure, MPa, and the Langmuir pressure is the pressure at which the adsorption volume reaches 50% of the ultimate adsorption volume.The fracture system of the coal body provides a channel for the free gas to deposit and seep through, and the mass balance equation of the gas in the fracture system of the coal body per unit volume is(22)∂mf∂t+∇ρgq=Q,where mf is the mass of free gas per unit volume of coal, kg/m3; t is time, s; ρg is the gas density, kg/m3; q is the Darcy velocity ignoring the effect of gravity, m/s; q=−k▽p/μ, k is the coal permeability, m2; μ is the gas dynamic viscosity, Pa∙s; Q is the mass exchange rate between the matrix system and the fracture system per unit volume of coal, kg/(m3∙s).Under isothermal conditions, the free gas content is only related to the gas pressure. According to the gas equation of state, the mass of free gas in a unit volume of gas-bearing coal can be obtained as(23)mf=φρg=φρnpnp,whereφ is the porosity of the coal body, %; ρn is the density of gas at standard atmospheric pressure, kg/m3; pn is the standard atmospheric pressure, MPa.The diffusion of gas in the coal body is mainly in the gas phase. According to the law of conservation of mass, the rate of change of gas mass per unit volume of coal matrix system with time is equal to the rate of mass exchange per unit volume of coal matrix system with the fracture system.(24)∂mγ∂t=−Q,where mr is the mass of adsorbed gas per unit volume of coal, kg/m3. Considering the influence of moisture, ash, and combustible percentage in coal on gas adsorption, the mass of adsorbed gas per unit volume of coal can be obtained from the modified Langmuir equation, as shown in
(25)mγ=abp1+bpρgρc100−A−M1001+0.31M,where a is the ultimate gas adsorption capacity of the coal body, m 3/kg; b is the adsorption constant of the coal body, MPa-l; ρc is the density of the coal body, kg/m3; A is the ash content of the coal body, %: M is the moisture content of the coal body, %.LetR=ρgρc100−A−M/1001+0.31M, then Equation (26) can be obtained
(26)∂φρnp/pn∂t=∇kρnpμpn∇p=−∂abpR/1+bp∂t.


## 2.4. Porosity and Permeability Evolution Equations

Coal porosity and permeability are the key coupling terms in the coupling of deformation field and gas transport field. By the definition of porosity, the coal porosity can be obtained by deducing that(27)φ=VpVB=1−1−φ01+εv1+ΔVsΔVs0,where Vp is the pore volume, m3; VB is the coal body appearance volume, m3; φ0 is the initial porosity, %; ΔVs is the change in coal body skeleton volume, m3; ΔVs0 is the initial coal body skeleton volume change, m3.The volume deformationΔVs/ΔVs0 includes the deformation ΔVsp/ΔVs0 caused by changes in gas pressure, the deformation ΔVsf/ΔVΔVs0 caused by the expansion of gas adsorbed by coal particles, and the deformation ΔVst/ΔVΔVs0 caused by the thermal expansion of coal skeletal particles.
(28)ΔVsΔVs0=ΔVsp+ΔVsf+ΔVstΔVs0=31−2vEΔp+εLpp+pL1−φ0+αΔT,where E is the modulus of elasticity of the coal body, GPa; Δp is the change in gas pressure, MPa; εL is the Langmuir volume strain constant; PL is the Langmuir pressure constant, MPa; α is the volume expansion coefficient of the coal body, K-1; ∆T is the absolute temperature change, K.Considering the effects of gas pressure, coal rock volume strain, and temperature on porosity, the dynamic expression for the porosity of gas-bearing coal under isothermal conditions is derived as(29)φ=1−1−φ01+εV1−31−2vEΔp+εLpp+pL1−φ0,where φ0 is the porosity of the coal body under the initial conditions; εV is the volume strain of the coal body under the corresponding gas pressure conditions.Based on the Kozeny-Carman equation and defined by porosity and volumetric strain, this gives(30)kk0=ϕϕ03=11+εV1+β2Vp0+ΔVpVp03,where k0 is the initial coal body permeability, m2; β is the ratio of the change in pore internal surface area to the initial internal surface area of the pore, which is much smaller than 1 and therefore negligible; Vp0 is the initial pore volume of the coal body, m, and ΔVp is the change in the pore volume of the coal body, expressed by the following equation:
(31)ΔVp=εvVB0−εLp1−φ0PL+p−31−2vEΔpVs0,where VB0 is the initial appearance volume of the coal rock body, m3; VS0 is the initial skeletal volume of the coal rock body, VS0=VB0−VP0, m3.Letγ=εLp/1−φ0PL+p−31−2v/EΔp simplify the above equation to give
(32)kk0=11+εv1+γ+εV−γφ03.The resulting equation for the dynamic evolution of coal body permeability is(33)k=k01φ0−1−φ0φ01+εV1−31+2VEΔp+εLpp+PL1+φ03.


## 2.5. Gas-Solid Coupling Model

Gas movement with the gas flow in confined space is essentially a single-phase multicomponent diffusion problem without chemical reaction. Its motion law satisfies the basic equations of fluid mechanics described by general Navier-Stokes equations, including continuity equation, momentum conservation equation, energy conservation equation, and component mass conservation equation. The control equations of coal and rock deformation and gas migration are unified to form a gas-solid coupling model. The couplings (8), (15), (19), (25), (28), and (32) together form a mathematical model for gas-solid coupling in the permafrost region close to the coal seam, see
(34)Gui,jj+G1−2vuj,ji−αpp+Fi=0,∂φρnp/pn∂t−∇kρnpμpn∇p=−∂abpR/1+bp∂t,φ=1−1−φ01+εV1−31−2vEΔp+εLpp+PL1−φ0,k=k01φ0−1−φ01+εV1−31−2vEΔp+εLpp+PL1−φ03.


## 3. Project Case

The Muli coalfield is located in the alpine region of the northern Qinghai-Tibet Plateau, in the upper basin of the Datong River at the junction of Haibei Tibetan Autonomous Prefecture and Haixi Mongolian Tibetan Autonomous Prefecture in Qinghai Province, in the adjoining area of Tianjun, Gangcha, Qilian, and Haiyan counties in northern Qinghai Province, starting from Haiyan Pass in the east to the watershed of the Datong and Shule rivers in the west, bounded by the Datong and Tolai mountains in the north and south, respectively, including Jiangcang, Juhuqian, Hot Water, Arc Mountain, Waihi Hada, Haider, and Merle 7 mining areas a total of more than 20 well fields. It is shown as Figure5.Figure 5
Illustration of the location of the Muli coalfield.The MNT coal mine is located on the south bank of the Datong River in the northeastern part of Qinghai Province and is the first large test shaft coal mine under construction in the Jiangcang mining area of Qinghai. The NMT coal mine has complex engineering conditions, and it is located in the permafrost area of Muli coalfield. The flake permafrost area of the NMT coal field accounts for 80% of the total area, and the permafrost layer is thick. The failure characteristics of the permafrost layer are different from ordinary strata, but continuous failure, with prominent creep mechanical characteristics. The minefield is approximately 3.47 km long from east to west and 1.81 km wide from north to south, covering an area of approximately 6.3 km2, with a resource reserve of 212 Mt, a design production capacity of 1.2 Mt/a, and a design service life of 67 years. The mine is located in the high cold, low mountainous, and hilly area of the plateau, the ground elevation reaches about +3820~4000 m, the mine is located at an elevation of about +3850 m, the four seasons are not very distinct, the winter is very cold, the daily temperature difference is large, the lowest temperature is -34°C, the highest is 19.8°C, the average annual temperature is -4.2°C~5.1°CC, and it belongs to the continental climate of the plateau, the surface of the fourth system The strata are perennial permafrost layers that do not melt all year round.The perennial permafrost in the Muli coalfield was formed during the Quaternary pale glacial period, and the neotectonic movement caused it to rise continuously, maintaining the alpine character and preserving the permafrost to this day. The permafrost in the coalfield is divided into a continuous permafrost distribution zone, a permafrost sheet distribution zone, and a permafrost island distribution zone. The perennial permafrost sheet distribution area is located between the snow line and 3,800 m below, and the MNT coal mine is located in this area. The perennial permafrost layer is found throughout the entire wellfield, and the Fourth Series strata are basically perennial permafrost layers.The thickness of the MNT permafrost layer varies greatly, thickening from south to north and thinning from west to east, and the distribution pattern of the thickness of the Quaternary strata is roughly the same, except for the melt zone under the riverbed of Jiangcang, there is generally a distribution of multiyear permafrost layers ranging from about 30 to 80 m thick, and the thickness of the permafrost layer increases with the thickness of the Quaternary surface layer, with the thickest point close to 100 m, and the permafrost layer penetrates deeply into the bedrock. The thickness of the permafrost layer measured in some boreholes of the geological exploration construction is shown in Table1, from which it can be seen that the maximum thickness of the multiyear permafrost layer is 89.15 m in hole 83 of exploration line 26, and the minimum thickness is 24.00 m in hole ZK28-1 of exploration line 28.Table 1
Statistical table of the thickness of permafrost layer drilled in the NMT minefield.
Exploration line-numberThe thickness of frozen soil (m)Floor elevation of frozen soil (m)29-9084.253895.826-8389.153780.025-8553.203819.129-8155.003829.5ZK30-136.293823.9ZK28-122.503851.9Combining the test results of mechanical parameters of permafrost specimens for many years and referring to the information about permafrost in the Qinghai area, physical composition, temperature, and water content of permafrost layer for many years, the mechanical parameters of the fourth series (Q) permafrost layer in NMT mine are: uniaxial compressive strength is 5.7 MPa, density 2000 kg/m3, modulus of elasticity 2.7 GPa, Poisson’s ratio 0.3, angle of internal friction 30%, cohesion 2 MPa, and tensile strength 1.63 MPa.A total of 8 coal seams are mined, divided into 2 groups, the upper coal group is nos. 2, 3, 4 (3+4), and 5 coal seams of the Shanxi Group, and the lower coal group is nos. 6, 8, and 9 coal seams. The main coal seams are nos. 3, 4, and 5, with average seam thicknesses of 1.15 m, 2.95 m, and 2.65 m, respectively. The seams contain a small number of intercalated gangue layers, all of which are stable and mineable in the whole area, with a minimum seam spacing of 5.6 m. Therefore, it is a typical close seam group mining, as shown in Figure 6.Figure 6
Coal seam distribution map.There is a correlation between the gas content of the coal seam and the depth of burial of the coal seam, the thickness of the overlying bedrock, and the lithology of the top and bottom slabs. Therefore, I took coal samples from the main mining seam of the mine area during the first period to study the relationship between the coal seam gas content and the thickness of the overlying bedrock, and the burial depth of the coal seam as shown in Figure7. As can be seen from Figure 7, the gas content of the coal seam has a linear relationship with the depth of burial of the coal seam and the thickness of the overlying bedrock, but as the depth of burial of the south wing coal seam in the mine area is generally greater than that of the north wing coal seam, the gas content of the south wing coal seam is generally higher than that of the north wing coal seam under the same coal seam base elevation conditions.Figure 7
Relationship between the gas content of the coal seam and the thickness of the overlying bedrock and depth of burial of the coal seam.The pore distribution characteristics of the coal seams in the mine area are shown in Table2, the largest proportions of medium, micropores, and small pores are in coal seams 3#, 4#, and 5#, accounting for 83.57%, 66.31%, and 69.57% of the total pores, while the main flow places of gas in the coal seams are visible pores, large and medium pores, and micropores, and small pores are the main storage spaces for adsorbed gas. The mechanical properties of the rock and coal seams are shown in Table 3.Table 2
Pore distribution characteristics.
Coal seamMicroporePinholeMesoporeMacroporeSsφ/%Vs/(mL/g)φs/%Vp/(mL/g)φ0/%Vp/(mL/g)φ0/%Vp/(mL/g)φ0/%3#0.016556.510.007927.060.00144.790.003411.643.7563.554#0.025943.550.013622.760.009515.900.010317.273.8932.685#0.022745.040.012324.530.009117.970.006911.612.4553.53Note:Vs: pore volume per unit mass, φs: specific pore volume, Ss: specific surface area.Table 3
Characteristic parameters of coal seam and rock stratum.
ParameterRock stratumCoal seamDirect roofWaste inclusionDirect base plate3#4#5#E/GPa23.115.221.94.74.13.3c/MPa3.13.21.53.33.02.8Φ/°33.231.737.824.123.923.8υ0.240.200.240.320.300.29ρ/(kg/m3)261026902690132013381350a/(m3/kg)///9.9×10−311.0×10−313.0×10−30///9.6×10−189.2×10−189.0×10−18p0/MPa///1.21.31.5μ/(Pa∙s)///1.08×10−61.08×10−61.08×10−6A/%///0.971.051.09M/%///15.9112.319.86PL/MPa///0.520.580.69L///0.010.010.01


## 4. Results and Discussion

According to the definition of the thin plate theory, when the ratio of the thickness of the plate to the minimum width is greater than 1100-1/80 and less than 1/8-1/5, it can be regarded as a thin plate. It is in a large area of overhang, which is consistent with the elastic thin plate theory. Therefore, the deformation characteristics and movement patterns of the tundra crust can be analyzed using the thin plate theory.Based on the established theoretical model, it can be seen that in the mining of the large inclined coal seam group under the permafrost layer in the NMT coal mine, the tilting and rotating composite roof structure of the roof slab has caused tensile damage at the two ends and the middle of the tilting direction of the permafrost layer “slab cover,” although it cannot form a macroscopic fracture phenomenon, with its continuous tilting. The left and right sides of the tilted direction and the central part of the tilted direction will be fractured by cracks. The permafrost “slab cover” in the direction of the strike forms a thick slab that is similar to the left and right ends of the break free, while the work is supported in the front and rear directions, at which point the permafrost “slab cover” can be seen as a thick strip beam structure in the direction of the strike. The collapse of the permafrost “slab cover” in the direction of travel is mainly caused by its own gravity.In the initial stage of mining of the large dip coal seam group under the permafrost layer in the NMT coal mine, the permafrost layer does not have a delamination structure with its lower part, or the delamination structure is small; the permafrost layer is more integral and has a higher tensile and shear strength. The perennial permafrost “slab cover” does not break down near the midpoint of the long and short supports. As the coal seam is mined and the protective rock beam is formed, the structure becomes more complex as the overhang distance increases and the forces on the permafrost “slab cover” become more complex. The slab cover is a composite structure with a perennial permafrost layer, creating a delaminated structure.According to the results of the theoretical analysis, there is a mechanical coupling relationship between the mechanical damage of the permafrost “slab cover” and the breakage of the lower coal rock body. The main relationship is: the rotational tilting movement of the protective rock beam will produce stress on the left support point of its tilting direction, and the supporting rock will have a tendency to move to the left, forming tensile stress on the upper part of the overlying ground surface of the permafrost layer “slab cover,” driving its overall tilting direction to the left. As the lower part of the permafrost layer “plate cover” is suspended, its own gravity sinks and bends to form tensile stresses on both sides of the permafrost layer, which is balanced by tensile stresses on the right side of the inclination. In the direction of the working face, the perennial permafrost layer “plate cover” due to its own gravity in the direction of both sides of the tensile stress for the perennial permafrost layer “plate cover” above the force is uniform. There is also a tendency to move to the left of the dependent works, resulting in widespread tensile damage near the fixed boundary at the right end of the mechanical model.In this paper, the mathematical model is established based on theoretical analysis, and the boundary conditions are given with the field test data. Then, the parameters are inversely calculated by the Matlab software. As the working face progresses, when the permafrost “slab cover” overhang distance 1 reaches the initial collapse stepLc, x=0, z=Ht/2, σx reaches the permafrost ultimate tensile stress subject to the maximum tensile stress σt, thus obtaining the initial permafrost “slab cover.” The first step across the break is calculated by
(35)Lc=Ht4σt+qt4+3vqt=Ht4σt+γtHt4+3vγtHt.After the initial span break of the permafrost layer, a cantilever beam structure is formed, and with the continuous advance of the working face direction, when the overhang distance 1 reaches the initial collapse stepLz, x=0, z=−Ht/2, σ reaches the permafrost ultimate tensile stress subject to the maximum tensile stress σt, thus obtaining the permafrost “plate cover” cycle span break step
(36)Ls=Htσt+1/5qt3qt=σt+1/5γtHtHt3γt.Substituting the mechanical parameters of the Quaternary (Q) permafrost layer at the NMT mine, based on an average permafrost thickness of 60 m, the initial collapse step in the strike direction and the periodic collapse step are
(37)Lc=60×4×1630000+20000×60×4+3×0.320000×60=192.87m,Ls=1630000+1/5×20000×60×603×20000=43.24m.In comparison to the initial and periodic collapse steps for general geological conditions, the perennial permafrost has significantly larger values for the initial and periodic collapse steps due to its low density and high tensile strength, which requires attention in coal mining practice. As the lower part of the permafrost, the crust is less detached, and the deformation breakage of the permafrost cover is also minimal and does not cause damage to water resources or surface ecology.The stress field-fissure field evolution and distribution characteristics of the overburden studied indicate that coal seam clusters are likely to produce large fissures or penetrating fissures under stacked mining conditions; however, the fissure penetration between coal seams has an important influence on the determination of gas transport and extraction parameters. Therefore, to investigate the gas transport and enrichment characteristics of the proximal coal seam group under stacked mining conditions, SF6 tracer gas identification was used to verify the penetration of interseam fissures and the degree of influence of penetrating fissures in the field.The test site was the underground drilling site of the return tunnel of the (3#+4#) coal seam, the target coal seam was the (3#+4#) coal seam, the sampling site was the extraction area of the working face of the no. 2 coal seam and the gas extraction borehole of the (3#+4#) coal seam; 5# measurement lines were arranged at the 2# ends of the working face and the floor of the extraction area, respectively, and gas samples were collected utilizing preburied pipes, as shown in Figure 8.Figure 8
Arrangement of gas sampling points in the mining area of the working face.As SF6 gas was found on the measurement lines in the mining area of the no. 2 coal face, it indicates that a penetration fissure had been created between the no. 2 coal face and the (3#+4#) coal face at this time; no SF6 gas was found in the first nine gas retrievals, but SF6 was found in the 10th gas retrieval, indicating that the face was mined back to the position at the time of the 10th gas retrieval, and a penetration fissure was created between the no. 2 coal face and the (3#+4#) coal face. SF6 gas was first detected at measurement points 2-3 and 4-3, respectively, indicating that a floor penetration-type fissure was generated at a location 30-40 m behind the protective seam workings. From the three-dimensional map of SF6 time concentration distribution in the mining area, it can be seen that according to the time when SF6 gas was first detected at measurement points 2-3 and 4-3 and the distance the working face advanced at this time, it indicates that the interseam penetration fissure was generated at about 40 m behind the working face of coal seam no. 2, and then, SF gas was also detected at other measurement points (measurement points 2-2, 2-4, 4-2, and 4-4) at the same time. This indicates that not only did the (3#+4#) seam produce a penetration fissure but also that secondary fissures were developed in its vicinity, and that SF66 gas originating from the (3#+4#) seam was transported through the penetration fissure to the area 20-60 m behind the working face of the no. 2 seam. The result are shown as Figure 9.Figure 9
Gradient diagram of SF6 concentration changes in the mining area of the working face.
(a)(b)(c)(d)(e)Therefore, according to the mining stress evolution-fissure development-gas transport characteristics under the overlapping mining conditions of the close coal seam group, it is found that the effect of overlapping mining on the overburden stress field and fissure field is not a simple superposition of effects, but a “1+1>2” effect, so that the overburden coal seam produces penetrating fissures under the overlapping mining, and a large number of secondary fissures are derived around it, which timely provides a superior channel for coal seam gas. This provides an advantageous channel for the transport of coal seam gas.Based on the requirement of “simultaneous management” of coal seam gas and protrusion hazards, the finite difference method was used for the two close coal seams 4# and 5# based on the established gas-solid coupling model. The gas pressure variation under different coal discharge (unloading) conditions at the same time were analyzed. The pressures at the boreholes in the 4# and 5# seams were calculated for coal volumes of 0, 3, 5, and 8 t, respectively. The gas pressure distribution at the calculation point after 60 d of extraction is shown in Figure10.Figure 10
Variation of gas pressure at different coal output after 60 d of extraction.
(a)(b)From the calculation results, it can be seen that under the condition of constant extraction time, the effective radius of influence of extraction is proportional to the amount of coal produced. At a certain extraction time, the coal output from the 4# coal seam is 0, 3, 5, and 8 t, respectively, and the gas pressure 0.74 MPa isopotential line is 1.65, 5.12, 7.30, and 8.06 m from the borehole wall, respectively, while the coal output from the 5# coal seam is 0, 3, 5, and 8 t, respectively, and the gas pressure 0.74 MPa isopotential line is 0.86, 2.55, 3.91, and 4.40 m from the borehole wall, respectively. This shows that the larger the coal output, the larger the effective radius of influence and the better the gas extraction effect. However, in the actual construction process, in order to ensure safety during the recovery period, the amount of coal exiting should not be too large.To break through the problem of gas extraction from close coal seams and realize the comining of coal and gas in a multiseam three-dimensional area, this paper proposes the “Protected Seam Mining + Directional Drilling Group in Bottom Extraction Lane Comining Technology” based on the spatial and temporal conditions and requirements of coal-gas comining at different stages, which can release high-quality coal and coal-bed methane production capacity while ensuring safe production in mines.In a mining area of the mine, for example, the no. 2 coal seam is used as a protective layer to mine a total of 5 working faces from 2203 to 2207 (Figure11), and the no. 2 bottom extraction road is arranged in the bottom plate of the no. 5 coal seam below the 2205 working face; the length of the roadway is 1284 m. 10 drilling fields (10 directional holes) are arranged in the north and south of the roadway, and a group of an upward directional hole is constructed along the coal seam, with a hole diameter of 113 mm, hole The final hole height is located in the no. 3+4 and no. 5 coal seams, respectively, each covering 5 working faces corresponding to the upper and lower coal of these 2 seams, and the no. 2 coal seam adopts 120 mm large hole diameter progressive presumping.Figure 11
Design of directional drilling group prepumping in protective layer mining + floor extraction roadway.Statistical analysis of the pure extraction volume of the 1,000 m drill hole and the common through seam drill hole in the 4203 working face of the second mining area shows that the pure volume of the 1,000 m drill hole is 5.05~10.88 m/min3, which is about 2.5 times of the common through the seam drill hole, and the pure volume of the common drill hole gradually decreases as the working face advances, while the extraction volume of the long directional drill hole is basically greater than 6 m/min. The daily production of the working face was increased from 4850 t 3to 6100 t, an increase of 25.8%, and both coal and CBM production were significantly increased. The gas drainage pure volume comparison of directional drillings and ordinary through layer drillings are shown as Figure 12.Figure 12
Gas drainage pure volume comparison of directional drillings and ordinary through layer drillings.On the one hand, the paper reveals the fracture forms, failure characteristics, and spatial structure relationship of coal roof and permafrost during the mining process of coal seam group with large dip angle under permafrost overburden. On the other hand, it is necessary to observe and study the permafrost surface deformation in the next production practice to verify the correctness of the research conclusions. The stability control of overlying rock in the mining of coal seam with large dip angle under the permafrost layer needs further study.


## 5. Conclusions

In this paper, a theoretical model of small deflection of thin plate is established based on Kirchhoff’s assumptions, and the deformation characteristics of the “cover” layer in the permafrost zone are analyzed; a gas-solid coupling model is derived and established for the close coal seam gas-integrated mining face, and the following conclusions are obtained for the NMT coal mine.(1)
There is a tendency for the multiyear permafrost layer at the NMT mine to move to the left, which will produce extensive pulling damage near the right end. The initial and periodic collapse steps are 192.87 m and 43.24 m, respectively, which are significantly larger than the initial and periodic collapse steps under normal geological conditions but will not cause damage to water resources or surface ecology(2)
The pore and fissure structural characteristics of the coal seams in the NMT mine are not conducive to gas transport. The effect of overburden mining on the overburden stress and fissure fields is not a simple superposition of effects, but a nonlinear incremental effect, with penetrating fissures generated under overburden mining and a large number of secondary fissures derived in the surrounding area, which provide an advantageous channel for coal seam gas transport in time(3)
Based on the established gas-containing gas-solid coupling model, the gas pressure evolution characteristics of the proximity coal seam are analyzed. The results show that under the condition of constant gas extraction time, the effective influence radius of extraction is directly proportional to the coal output. The larger the coal output, the larger the effective radius of influence and the better the gas extraction effect(4)
Given the difficulties of gas extraction, based on the spatial and temporal conditions and requirements of coal-gas comining in different stages, the “protected layer mining and bottom extraction lane directional drilling group comining technology” was proposed. Both coal and CBM productions were significantly increased



---

*Source: 2379040-2022-07-04.xml*