Abstract

The solid-gas coupling model for mining coal-rock mass deformation and pressure relief gas flow in protection layer mining is the key to determine deformation of coal-rock mass and migration law of pressure relief gas of protection layer mining in outburst coal seams. Based on the physical coupling process between coal-rock mass deformation and pressure-relief gas migration, the coupling variable of mining coal-rock mass, a part of governing equations of gas seepage field and deformation field in mining coal-rock mass, is introduced. Then, a new solid-gas coupling mathematical model reflecting the coupling effects of gas adsorption/desorption, gas pressure, and coal-rock mass deformation on the mining coal-rock mass deformation and pressure relief gas flow is established combined with the corresponding definite conditions. It lays a theoretical foundation for the numerical calculation of the deformation of mining coal-rock mass and the migration law of gas under pressure relief in the outburst coal seam group.

1. Basic Hypothesis

After mining protection seam, the initial rock stress balance is broken and the coal-rock mass deforms, which destroys the dynamic balance of original gas adsorption and desorption, and part of gas desorbs and diffuses to the mining fissure and gas pressure reduces, which makes the coal-rock mass pressure state changed; namely, deformation of mining coal-rock mass and flow of pressure relief gas are a multiphysical field coupling process. Therefore, the establishment of a solid-gas coupling mathematical model which can describe the deformation of mining coal-rock mass and pressure relief gas flow is the key to study the deformation of coal-rock mass and the migration law of gas under pressure relief in protection layer mining.

Due to the solid-gas coupling effect of coal-rock mass deformation and pressure relief gas flow are complicated, which involves seepage mechanics, rock mechanics, elastic-plastic mechanics and many other subjects [15]. The existing research results found that seepage pressure affects the stress field; moreover, the stress field in surrounding rock results in the change of permeability of rock masses [3, 5]. Based on the previous studies [15], the following hypotheses are introduced:(1)The solid-gas coupling system consists of both solid phase (coal-rock mass) and gas phase.(2)The solid skeleton deformation of coal-rock mass is small, which can be regarded as homogeneous and isotropic linear elastic medium; the coal seam gas can be regarded as an ideal gas, and its flow process can be regarded as approximately isothermal process, which obeys the Darcy seepage law.(3)The coal seam is saturated by single-phase gas, and the original gas pressure is uniformly distributed; gas content in coal seam follows the modified Langmuir equation.(4)The effective stress of coal-rock mass skeleton follows the modified Terzaghi effective stress principle.

2. Control Equations of Deformation Field of Mining Coal-Rock Mass

The control equations of displacement and stress fields of mining coal-rock mass deformation are mainly composed of constitutive equation, geometric equation, and equilibrium equation.

2.1. Constitutive Equation

For the single-phase saturated linear elastic porous media, the effective stress (tension is positive), the pore pressure p (compression is positive), and the total stress follow the modified Terzaghi effective stress principle [6, 7]:where is the Kronecker symbol (when i = j, , and when i ≠ j, ) and is the Biot coefficient, whose value depends on the compressibility of coal-rock mass and is calculated as follows:where is the bulk modulus of coal-rock mass skeleton (Pa) and is the matrix bulk modulus of coal-rock mass (Pa).

Previous studies have shown that [8, 9] adsorption/desorption will cause the bulk strain of coal matrix, and it can be calculated combined with the generalized Hooke’s law:where is the Poisson’s ratio of drainage, G is the shear modulus (Pa), and is the matrix bulk strain caused by coal adsorption/desorption; the experimental results show that the relationship between and p meets the Langmuir type curve [8, 9]:where is the bulk strain constant when the pore pressure is infinite and is the pore pressure when the bulk strain is (Pa).

2.2. Geometric Equation

Assuming that the deformation of coal-rock mass is small, the geometric equation is [10]where is the displacement in the i direction.

2.3. Equilibrium Equation

According to the virtual displacement principle of elastic mechanics, the momentum equation of coal-rock mass skeleton under quasi-static condition is [11]where and are the boundary surface force and the bulk force of the skeleton, and are the coordination of virtual displacement variation and virtual strain variation, and and are the regions occupied by solid skeleton and the known boundary of stress.

It can be obtained according to (5) and the symmetry of strain tensor [12]:

Substituting (1) into (6), we get

According to the Green-Gauss formula, the first item on the left of (8) is

The equilibrium equation of porous medium of coal-rock mass can be obtained by (8) and (9):

Connecting the constitutive equation (3), the geometric equation (5), and the equilibrium equation (10), the governing equations which take the solid-gas coupling deformation field under the effect of coal adsorption/desorption deformation, stress, and pore pressure into account can be obtained:

3. Governing Equation of Gas Flow in Mining Coal-Rock Mass

The governing equation of gas flow in coal-rock mass combines the continuity equation, state equation, gas content equation, and Darcy law of motion and eliminates some variables to deduct differential equations. Migration morphology of gas is closely related to the pore structure, but the pore structure of coal seam is complicated; when studying migration law of gas, people are often based on different types of gas diffusion or the primary and secondary roles of seepage to simplify the pore structure of coal seam into pure diffusion pore medium model, pure permeability fracture medium model, and pore-fissure dual-medium model with coexistence of diffusion and permeation. The pore-fissure dual-medium model has been accepted and used by many scholars, and it assumes gas migration in coal seam is a continuous process in two steps: diffusion and permeation are a cascade process that the gas diffuses from the matrix to the fissure in the form of diffusion and then flows through the fissures into the tunnel or borehole in the form of the Darcy seepage.

The continuity equation of gas diffusion motion is as follows [13, 14]:where mc is the mass diffusion flux (kg/(m2·t)), is the exchange mass source of the pore system flowing to the fracture system (kg/(m3·s)), and C is the quality of adsorbed gas in unit bulk coal:where is the gas density under the standard state (kg/m3), is density of coal-rock mass (kg/m3), and is the Langmuir bulk constant (kg/m3).

According to the law of mass conservation, the continuity equation of gas seepage in fracture can be obtained [1315]:where is the gas density (kg/m3), is the porosity of coal-rock mass, and is the velocity of gas seepage (m/s).

The gas flow is divided into two open systems: diffusion and seepage. In fact, there is a coupling effect between gas diffusion and seepage, and the adsorption and desorption of gas only change in the state of matter, and it does not cause the change of gas quality in representative bulk units. Therefore, in order to facilitate the analysis and solution, the mass conservation equation of gas flow can be obtained by adding (12) and (14):

For low permeability coal-rock mass, it can be assumed that the gas is desorbed and diffused into the fracture from the matrix instantaneously, and (15) can be further simplified:

In (16), , the first item on the right shows the free gas content in the coal-rock mass, and the second item represents the adsorbed gas content in the coal.

The coal gas, whose flow process could be assumed to be isothermal, is considered as the actual gas, and then the state equation of the gas can be obtained:where is the molar mass of gas, R is the universal gas constant, and T is the absolute temperature.

The outward coal-rock mass gas flow from the coal seam crack accord with the Darcy seepage law:where k is the permeability of the coal seam, is the dynamic viscosity of gas (Pa·s), and is the gas pressure gradient (Pa/m).

Substituting (17) and (18) into (16), we getwhere is the standard atmospheric pressure (Pa). In (19), the permeability k depends on the porosity of coal , and the porosity of coal is a function of the volumetric strain of coal-rock mass and the volumetric strain of matrix. Therefore, (11) and (19) are connected by the porosity of coal and physical parameters of permeability k.

4. The Structure of Solid-Gas Coupling Mathematical Model for Coal-Rock Mass Deformation and Pressure Relief Gas Flow in Protection Seam Mining

4.1. Solid-Gas Coupling Variables in Mining Coal-Rock Mass

Traditional classical percolation mechanics holds that in the pore pressure change process, the solid skeleton of porous media does not produce any elastic-plastic deformation, so the porosity of porous media is a constant. However, under the effect of surrounding rock stress and gas pressure, the solid skeleton of coal will deform in different degrees, which will cause the change of effective stress and lead to the dynamic change of porosity and permeability of coal-rock mass. When it comes to the deformation of coal-rock mass in protection layer mining and the solid-gas coupling of pressure relief gas flow, the dynamic changes of porosity, permeability, and other parameters of coal-rock mass must be taken into account. So, it is necessary to establish the theoretical calculation model for the dynamic changes of the basic physical parameters such as porosity and permeability of coal-rock mass.

4.1.1. Dynamic Variation Model of Coal-Rock Mass Porosity

The bulk V of porous medium in coal-rock mass is composed of the bulk of solid phase Vs and the bulk of pore volume Vp, namely, V = Vs + Vp. The definition of porosity is , and the bulk strain of coal-rock mass can be obtained by (11) [16, 17]:where , , and is the module of pore bulk.

It is assumed that the strain caused by coal adsorption is equal to the amount of pore variation. Without considering the adsorption effect of gas, the bulk change of porous media meets the equivalent law of Betti-Maxwell:and therefore,

Then, deducted from the definition of porosity

Connecting (20), (21), (24) and (25), we get

Substituting (20) and (23) into (26), we get

When the pore pressure of coal-rock mass is , the initial porosity is , and the initial bulk strain is 0; then, the porosity can be expressed as

In formula (28), and .

4.1.2. Dynamic Change Model of Coal-Rock Mass Permeability

Permeability, a physical parameter reflecting the difficult degree of gas flow in coal seam, is important for studying coupling process between deformation of mining coal-rock mass and gas flow. The relationship between the permeability of coal seam and the bulk strain of coal can be deduced by the Kozeny-Carman equation:

Neglecting the bulk deformation of coal, (29) can be simplified as

Connecting (28) and (30), we get

4.2. Solid-Gas Coupling Mathematical Model of Mining Coal Deformation and Pressure Relief Gas Flow

According to the antecedent, (11) is the control equation of mining coal deformation; substituting (4) into (11), we get

It can be obtained from (28):

Substituting (33) into (19), the control equation of the pressure relief gas of the mining coal can be obtained:

Formulas (32) and (34) constitute a gas-solid coupling mathematical model of coal seam deformation and pressure relief gas flow in protection seam mining.

4.3. Solid-Gas Coupling Mathematical Model of Rock Mass Deformation and Pressure Relief Gas Flow in Protective Seam Mining

When the rock mass is affected by mining, deformation and cracks occur, the change of rock porosity will lead to the change of permeability, and the pressure relief gas will flow under the effect of gas pressure gradient. The flow of gas in the mining rock mass is different from that of the mining coal, mainly in the absence of adsorption of the rock mass on the gas. Therefore, there is no need to consider the deformation of bulk expansion caused by the gas diffusion process and gas adsorption/desorption, so it is necessary to modify the solid-gas coupling mathematical model of mining coal deformation and pressure relief gas flow.

To sum up, the governing equation of mining rock mass deformation can be expressed as

Similarly, the control equation of pressure relief gas of mining coal can be expressed aswhere .

4.4. Definite Conditions

The control equations (32) and (3436) constitute the complicated solid-gas coupling mathematical model of the mining coal-rock mass deformation and the pressure relief gas flow. For a particular solution, additional conditions must be added, that is, initial conditions and boundary conditions. Definite conditions of the model include the boundary conditions of the deformation field of the mining coal-rock mass and the initial conditions and boundary conditions of the gas seepage [1722].

4.4.1. Boundary Conditions of Mining Coal-Rock Mass Deformation Field

The first kind of boundary condition: the surface force of the solid skeleton of the coal-rock mass is known as

The second kind of boundary condition: the surface displacement of the solid skeleton of the coal-rock mass is known aswhere is the unit normal vector of the surface solid skeleton of coal-rock mass and and are the known surface displacement functions and known surface stress distribution functions.

4.4.2. Initial Conditions and Boundary Conditions of Gas Seepage

(1) Initial Conditions of Gas Seepage

(2) Boundary Conditions of Gas Seepage. The first kind of boundary condition (the Dirichlet condition):

The second kind of boundary condition (the Neumann condition):

The third kind of boundary condition (equal flow at the interface):where and are the distribution function and the fluid flow function of the pressure p and time t at the point (x, y, z) of the outer boundary, respectively, and is the initial flow value at the outer boundary.

5. Conclusions

The solid-gas coupling effect of coal-rock mass deformation and pressure relief gas flow in protective layer mining is very complicated, involving rock mechanics, seepage mechanics, elastic mechanics, and other subjects. Based on the existing research results, the seepage control equation of mining coal-rock mass and the control equation of deformation field are established including solid-gas coupling variable of mining coal-rock mass, namely, the dynamic change equation of coal-rock mass porosity and permeability. According to the coupling equations, the corresponding definite conditions are put forward, and the solid-gas coupling mathematical model of the deformation of coal-rock mass and pressure relief gas flow in the protection layer mining are built. The details are as follows:(1)The control equations (32) and (35) of the coal-rock mass deformation in protection layer mining that contain the coupling term which reflects the pore gas pressure and gas adsorption/desorption effect is established, which can only be solved by connecting with the governing equations of gas seepage field.(2)The governing equations (34) and (36) of the flow field of the pressure relief gas flow in the mining coal-rock mass that include the coupling term reflecting the effect of coal-rock mass deformation is established, and it is necessary to connect with the control equation of mining coal-rock deformation field to be solved.(3)The combination of the solid-gas coupling control equation of mining coal-rock mass deformation and pressure relief gas flow with definite conditions is established, which constitute the solid-gas coupling mathematical model of coal-rock mass deformation and pressure relief gas flow, and it lays a foundation for numerical simulation research of coal-rock mass deformation and pressure relief gas flow law in following protection layer mining. Numerical simulations research will be showed in other papers.

Conflicts of Interest

The authors declare that there are no conflicts of interest regarding the publication of this paper.

Acknowledgments

This study was partially supported by the National Natural Science Foundation of China (51774131 and 51604111), the Natural Science Foundation of Hunan Province (2017JJ2082), and the Doctoral Scientific Research Foundation (E51778).