Fluid–Solid Coupling Model and Simulation of Gas-Bearing Coal for Energy Security and Sustainability

The optimum design of gas drainage boreholes is crucial for energy security and sustainability in coal mining. Therefore, the construction of fluid–solid coupling models and numerical simulation analyses are key problems for gas drainage boreholes. This work is based on the basic theory of fluid–solid coupling, the correlation definition between coal porosity and permeability, and previous studies on the influence of adsorption expansion, change in pore free gas pressure, and the Klinkenberg effect on gas flow in coal. A mathematical model of the dynamic evolution of coal permeability and porosity is derived. A fluid–solid coupling model of gas-bearing coal and the related partial differential equation for gas migration in coal are established. Combined with an example of the measurement of the drilling radius of the bedding layer in a coal mine, a coupled numerical solution under negative pressure extraction conditions is derived by using COMSOL Multiphysics simulation software. Numerical simulation results show that the solution can effectively guide gas extraction and discharge during mining. This study provides theoretical and methodological guidance for energy security and coal mining sustainability.


Introduction
Since the 21st century, global resource shortage and environmental pollution have become difficult problems in human sustainability development [1,2]. Policy makers in various countries have focused on sustainable energy and low-carbon development by proposing sustainable strategies and methods [3]. With the advent of low-carbon economy, coalbed methane, as a clean, efficient, and safe energy source, has been eliciting considerable research attention. This unconventional natural gas is mainly present in coal seams in free and adsorption states [4]. Coal is a typical dual-porosity/permeability system containing a porous matrix surrounded by fractures. The coal matrix is separated by a natural fracture network composed of butt cleats and face cleats. The cleat system provides an effective flow channel for gas. The developed pore structure is the main space for coalbed methane. During gas percolation of porous media, the effective stress of the porous medium skeleton changes because of the change in pore pressure; the porosity and permeability of porous media also change to a certain extent [5]. These changes can affect the gas flow in pores and the redistribution of gas pressure within a certain range. Therefore, in studying the migration rule and deformation characteristics of gas in a porous medium, such as coal, the interaction between the gas flow in a porous medium and the deformation of the porous medium body should be considered. The mutual coupling between the gas seepage flow field and the stress field in the porous medium should also be considered.
Through relevant laboratory tests and field practice, people have gradually realized that the seepage property of gas in coal is related to the mechanical properties of coal, gas pressure, in situ stress, temperature, and other factors. Sommerton et al. studied the effect of stress on the gas seepage property of coal [6]. Brace investigated the law of permeability change of rock mass under stress [7]. McKee et al. analyzed the relationship between stress and coal porosity and permeability and studied the phenomenon in which the depth of coal seams and effective stress increase, the width of cleats in coal seam decreases, and permeability decreases exponentially [8]. Enever et al. obtained the influence rule between the effective stress and permeability of coal [9]. In accordance with the generalized form of the power law, Sun established a mathematical model for the flow of compressible gas in coal seams. This model, which is called the nonlinear gas flow model, was based on the measured gas parameters in the Zhongmacun Coal Mine of Jiaozuo Mining Bureau. A numerical simulation of the pressure distribution of the homogeneous gas seepage flow field was carried out using various models [10][11][12]. Li et al. studied the relationship among coal adsorption expansion, deformation, porosity, and permeability in consideration of the adsorption deformation characteristics of the coal skeleton and obtained the relationship among porosity, permeability, and swelling deformation [13,14]. Tao et al. analyzed the problems existing in the theory of nonlinear gas flow and theory of fluid-solid coupling of coal-bed gas. The research achievements in the fields of nonlinear gas flow theory, gas flow theory of the geothermal field effect, and coal-bed gas fluid-solid coupling theory have been scrutinized extensively [15,16]. Zhu et al. considered the Klinkenberg effect and proposed a coupled mathematical model of solid deformation and gas flow [17][18][19][20].
However, the fluid-solid coupling model that scholars have established still has certain limitations. For example, the Klinkenberg and adsorption expansion effects on gas migration in coal were not considered simultaneously. Existing research shows that strain on adsorption expansion occurs after the coal body adsorbs gas. When the strain is limited by certain factors, adsorption expansion stress follows, which causes a certain degree of primary deformation of the coal skeleton and affects the development of coal pores. Thus, the effect of adsorption expansion on gas migration should not be ignored [21]. On the basis of previous studies, the author comprehensively considers the influence of the two factors of migration of coalbed methane in coal and establishes a mathematical model of fluid-solid coupling in a low-permeability coal seam. Thereafter, relevant partial differential equations are derived. The establishment of the model expands the theory of fluid-solid coupling of gas-bearing coals under multi-field conditions and clarifies the law of gas occurrence and seepage in coal. The establishment of the fluid-solid coupling equation can characterize gas flow in coal seams from the perspective of time and improve the research method of gas dynamic flow law. The rationality of the established mathematical model is verified by using a specific example of the effective extraction radius of coal mine gas. This model provides a theoretical basis for the design and layout of gas drainage boreholes in coal mining and a reasonable reference for decision makers to control coal mine gas effectively [22,23].

Basic Assumptions
(1) Coal containing gas can be regarded as an isotropic elastic medium.
(2) The coal seam is considered homogeneous, that is, the physical properties of each part of the coal seam are similar everywhere and do not change with a change in position. (3) The coal seam temperature is constant. (4) The gas contained in the coal seam is regarded as ideal gas and obeys the ideal gas state equation.
The gas desorption obeys the Langmuir equation. (5) The seepage characteristic of coalbed methane in coal meets the Klinkenberg effect. (6) The deformation of coal is small. (7) The overall deformation of coal rock consists of pore and fracture deformation.
(8) Single-phase saturated gas fluid exists in the coal seam, and only free and adsorbed states are available. (9) The model is isolated from the outside world and does not undergo any form of energy and material exchange.

Physical Property Parameter Model of Coal
The compression and adsorption expansion of deformation result in different degrees of deformation of the coal skeleton because of the changes in conditions of the coal seam, such as crustal stress and gas pressure. Owing to the different depths of the coal seam, the gas pressure and crustal stress also change to varying degrees, resulting in changes in coal seam porosity and permeability.

Deformation Mechanism of Coal Containing Gas
According to previous research results [24], two kinds of deformation mechanisms exist under the joint action of internal and external stresses of coal containing gas.
(1) Structural deformation: Relative dislocation occurs between coal particles because of external stress. The particle arrangement becomes increasingly compact, which causes the compression of the coal particle skeleton. The deformation caused by the two is called structural deformation, as shown in Figure 1. (2) Body deformation: The deformation of coal particles is mainly caused by the adsorption expansion and desorption shrinkage of coalbed methane, the compression of the coal skeleton by gas pressure, and the thermal expansion and contraction effect of temperature.  (7) The overall deformation of coal rock consists of pore and fracture deformation. (8) Single-phase saturated gas fluid exists in the coal seam, and only free and adsorbed states are available. (9) The model is isolated from the outside world and does not undergo any form of energy and material exchange.

Physical Property Parameter Model of Coal
The compression and adsorption expansion of deformation result in different degrees of deformation of the coal skeleton because of the changes in conditions of the coal seam, such as crustal stress and gas pressure. Owing to the different depths of the coal seam, the gas pressure and crustal stress also change to varying degrees, resulting in changes in coal seam porosity and permeability.

Deformation Mechanism of Coal Containing Gas
According to previous research results [24], two kinds of deformation mechanisms exist under the joint action of internal and external stresses of coal containing gas.
(1) Structural deformation: Relative dislocation occurs between coal particles because of external stress. The particle arrangement becomes increasingly compact, which causes the compression of the coal particle skeleton. The deformation caused by the two is called structural deformation, as shown in Figure 1. This study assumes that the coal seam is thermostatic. Thus, the deformation of the particles under the influence of temperature changes is ignored.

Porosity Mathematical Model
According to the relevant definition of porosity, the porosity change of coal can be expressed as follows: where is coal porosity, is the initial porosity of coal, is the pore volume of coal, is the initial pore volume of coal, is the total apparent volume of coal, is the initial total apparent volume of coal, ∆ is the variation in the pore volume of coal, ∆ is the total apparent volume change in coal, is the volume of the coal skeleton, ∆ is the volume variation of the coal skeleton, and is the volumetric strain of coal. This study assumes that the coal seam is thermostatic. Thus, the deformation of the particles under the influence of temperature changes is ignored.

Porosity Mathematical Model
According to the relevant definition of porosity, the porosity change of coal can be expressed as follows: where ϕ is coal porosity, ϕ 0 is the initial porosity of coal, V p is the pore volume of coal, V p0 is the initial pore volume of coal, V b is the total apparent volume of coal, V b0 is the initial total apparent volume of coal, ∆V p is the variation in the pore volume of coal, ∆V b is the total apparent volume change in coal, V s is the volume of the coal skeleton, ∆V s is the volume variation of the coal skeleton, and e is the volumetric strain of coal.
In the actual state, the volume strain increment of coal particles caused by the bulk deformation of coal particles is mainly composed of three parts. The first part is the strain increment caused by pore gas pressure compressing the coal particles. The second part is the increment in expansion strain caused by coal particles adsorbing the coalbed methane. The third part is the strain increment caused by thermoelastic expansion. The strain increment caused by thermoelastic expansion is zero because we assume that the temperature of the coal seam is constant. In reference to Figure 2, the relationship among the three parts can be expressed as follows: where ∆V sp V s0 is the strain increment caused by pore gas pressure compressing the coal particles, is the swelling strain increment caused by coal particles adsorbing the coal bed gas, and ∆V st V s0 is the strain increment caused by thermal elastic expansion. ∆V st V s0 is zero because the temperature of the coal seam is assumed to be constant. The total volume strain of coal particle deformation can be expressed as follows: where ε p is the expansion strain generated by the adsorption of gas per unit volume.
where T is the thermodynamic temperature of the coal seam (K), a is the limit adsorption capacity per unit mass of combustibles under a reference pressure (m 3 /Mg), b is the adsorption constant (MPa −1 ), R is the universal gas constant (R = 8.3143 J/(mol·K)), ρ is the coal density (kg/m 3 ), V m is the molar volume of gas (22.4 × 10 −3 m 3 /mol), and Ks is the volume modulus of the coal skeleton (Pa).
In the actual state, the volume strain increment of coal particles caused by the bulk deformation of coal particles is mainly composed of three parts. The first part is the strain increment caused by pore gas pressure compressing the coal particles. The second part is the increment in expansion strain caused by coal particles adsorbing the coalbed methane. The third part is the strain increment caused by thermoelastic expansion. The strain increment caused by thermoelastic expansion is zero because we assume that the temperature of the coal seam is constant. In reference to Figure 2, the relationship among the three parts can be expressed as follows: where ∆ is the strain increment caused by pore gas pressure compressing the coal particles, ∆ is the swelling strain increment caused by coal particles adsorbing the coal bed gas, and ∆ is the strain increment caused by thermal elastic expansion. ∆ is zero because the temperature of the coal seam is assumed to be constant. The total volume strain of coal particle deformation can be expressed as follows: where is the expansion strain generated by the adsorption of gas per unit volume.
where T is the thermodynamic temperature of the coal seam (K), is the limit adsorption capacity per unit mass of combustibles under a reference pressure (m 3 /Mg), is the adsorption constant (MPa −1 ), R is the universal gas constant (R = 8.3143 J/(mol•K)), is the coal density (kg/m 3 ), is the molar volume of gas (22.4 10 −3 m 3 /mol), and Ks is the volume modulus of the coal skeleton (Pa). Equations (3) and (4) are substituted into Equation (1) to obtain a mathematical model of the dynamic evolution of porosity.

Permeability Evolution Mathematical Model
Permeability is an important indicator that describes the difficulty of gas migration in gas-bearing coal seams, the permeability of coal seams, and gas drainage difficulty. Therefore, a  (1) to obtain a mathematical model of the dynamic evolution of porosity.

Permeability Evolution Mathematical Model
Permeability is an important indicator that describes the difficulty of gas migration in gas-bearing coal seams, the permeability of coal seams, and gas drainage difficulty. Therefore, a correct mathematical model of permeability evolution must be established for gas control in coal mines.
The permeability of gas-containing coal is closely related to the stress state of coal. Different stress states cause changes in coal-rock skeleton deformation and pore volume. When porosity changes, permeability also changes. To establish the relationship between coal permeability and porosity, we can refer to the Kozeny-Carman equation [25][26][27].
where k is permeability (md), k z is a dimensionless constant (k z = 5), S p is the pore surface area per unit pore volume of coal (cm 2 ), and A s is the total surface area of coal pore (cm 2 ). The permeability in the initial state is assumed to be where k 0 is the initial permeability (md), A s0 is the total surface area of the coal pore in the initial state (cm 2 ). The total volume of coal and the volume change of a single coal particle are ∆V b and ∆V s , respectively, when the initial state changes to a new one. Based on the definition of porosity, the new porosity is The new pore surface area can be expressed as follows: where ∂ is the increasing coefficient of the pore surface area of coal (%). For Equations (6) and (7), the ratio of the new permeability to the original permeability is computed by the following: The total surface area of the unit volume of coal particles is almost unchanged in the stress and strain process of coal [28]. This occurrence can be ignored. Thus, ∂ is approximately zero. According to [24], Equation (10) can be simplified as follows: The combination of Equations (3), (4) and (11) can be obtained: The above equation is a mathematical model for the evolution of the permeability of coal containing gas.

Effective Stress of Gas-Bearing Coal
Gas-containing coal is a complex deformable pore-fracture dual media. Coal has a strong adsorption capacity for gas and produces a certain adsorption expansion stress, which changes the stress distribution of coal.
Therefore, when studying the problem of fluid-solid coupling of coal containing gas and rock, the relationship between the effective stress of the coal seam and its adsorption-expansion stress should be considered simultaneously. In view of this problem, Wu et al. [29,30] derived the following formula for calculating the expansive stress of coal. The expansion stress formula of coal can be expressed as follows: where σ p is the expansion stress (Pa), v is the Poisson ratio, and E is the elastic modulus of coal (Pa).
In accordance with the effective stress law of Terzaghi and in consideration of the expansion stress absorbed by coal, the effective stress equation of gas-bearing coal can be expressed as follows: where σ ij is the effective stress of gas-bearing coal (MPa), σ ij is the overall stress of gas-containing coal (MPa), and α is the Biot coefficient.

Gas Content Equation
The existing mine gas in the coal bed occurs mostly in free, adsorption, and dissolved states. The dissolved gas content is not considered in this study because the amount is notably small. The gas content is the sum of the free-state and adsorbed-state gas contents. The adsorption gas content accounts for more than 90% of the total, and the content of adsorption gas is related to moisture, coal ash, and gas pressure. The free gas content mainly depends on coal porosity and the magnitude of gas pressure. According to previous research and the modified Langmuir adsorption equilibrium equation, the gas content of a coal seam can be obtained as follows [31]: where P n is the gas pressure under standard conditions (P n = 0.10325 MPa), and ρ n is the coalbed methane density under standard conditions (kg/m 3 ).
In the above equation, where Q is the gas content (kg/m 3 ), c is the coal quality correction parameter (kg/m 3 ), A is the ash content of coal (%), and M is the coal moisture (%).

Stress Field Equation of Gassy Coal
Assuming that the gas-bearing coal is an isotropic linear elastic medium, the stress field changes obey the following equation.
(1) Balance equation where F i is the bulk stress (N/m 3 ). According to the effective stress Equation (14) of gas-bearing coal, the equilibrium differential equation expressed by effective stress is obtained by the introduction of Equation (16). Then the equilibrium differential equation can be expressed as follows: (2) Geometric equation In the spatial distribution of gas-bearing coals, let u (x, y, z), v (x, y, z), and w (x, y, z) be the displacement components in directions x, y, and z, respectively, and continuous single-valued functions of coordinates. Then, the strain and displacement components satisfy the following geometric equations, which can be expressed as tensor symbols. (

3) Constitutive equation of gas-bearing coal
The constitutive equation of gas-bearing coal describes the relationship between the stress and strain of coal. The constitutive relationship in this study is based on the strains caused by the adsorption expansion of gas-bearing coal, compression of the coal particle body, and rock stress.
The linear strain caused by gas adsorption by coal particles is as follows: The linear compression strain of coal particles caused by the change in pore gas pressure is as follows: According to Hooke's law, the strain due to crustal stress is computed using the following: According to the above expressions, the total strain of gas-bearing coal is as follows By using the above formula as a reference, the following equation can be derived.
where G is the shear modulus (MPa), and Θ is the effective volume stress.
The following equation can be obtained by arranging Equation (23) after introducing the Lame constant.
where λ is the Lame constant.
Assuming that the coal is a linear elastic medium, the constitutive equation of gas-bearing coal-rock deformation conforms to Hooke's law, as follows: According to the stress-strain relationship of coal in Equation (32) and combined with the above formula, the effective stress constitutive equation of gas-bearing coal expressed in tensor form is derived as follows: (4) Stress field equation of gas-bearing coal Substituting Equation (26) into Equation (17) yields the following: This formula is the fluid-solid coupling stress field equation of gas-containing coal, and the deformation field of gas-bearing coal is represented by displacement. The variation of the strain field of gas-containing coal caused by crustal stress, gas adsorption by coal particles, and gas pressure change are considered in the equation.

Fluid-Solid Coupling Gas Seepage Equation of Gas-Bearing Coal
(1) Gas flow equation The previous experimental results reveal that when gas migrates in low-permeability gas-bearing coal seams, the gas molecules near the surface of the coal wall show the phenomenon of non-zero velocity, which does not conform to Darcy's law [32]. This occurrence is called the slippage effect or Klinkenberg effect in seepage mechanics. Its permeability can be expressed as follows: where k g is the Klinkenberg permeability, ω is the scale factor, λ 1 is mean free path of gas molecules, r is the average pore radius, m is the Klinkenberg coefficient (MPa), and ∇P is the gas pressure gradient in the coal seam (Pa/m). Therefore, when the Klinkenberg effect is considered, the equation of gas flow in the coal seam can be expressed by the following: where q is the velocity vector of gas flow (m/s), and µ is the gas dynamic viscosity (1.08 × 10 −5 Pa·s).
(2) Continuity equation According to the hypothesis, if the model is isolated from the outside and no exchange of substance and energy in any form occurs, the gas flow in the coal seam will conform to the law of conservation of mass, expressed in the form of a differential equation: where Q is the gas content in coal, ρ g is the gas density when the gas pressure is P (kg/m 3 ), and I is the source sink term.
(3) Gas seepage field equation of gas-bearing coal According to the state equation of coalbed methane in [33], the gas content equation (15), and the gas flow equation (29), the equation of gas seepage flow field can be obtained by combining them with Equation (30). The result is as follows: where M g is the molar mass of gas (kg/mol). Z is the gas compressibility factor, and its value is approximately 1 in the case of a small temperature difference. In summary, Equations (5), (12), (27), and (30) constitute a fluid-solid coupling model of gas-bearing coal.

Numerical Simulation of the Model and Analysis of Its Results
Deduction and analysis show that the fluid-solid coupling mathematical model of gas-bearing coal is a complex nonlinear equation group. The coupled numerical solution requires the use of a numerical method and COMSOL Multiphysics simulation software.

Geometric Model Establishment
The  Figure 3 shows the geometric diagram of the model. The model is divided into three layers, in which the coal seam is located in the two rock layers, the thickness is 4.22 m, and the upper part is loaded with a stable load of 11.7 MPa. The borehole is located in the middle of the model and has a diameter of 0.089 m and a negative pressure of 13 kPa. The initial gas pressure in the coal seam is 0.9 MPa, and the gas only migrates in the coal seam. The initial time is t = 0 day, and the simulation time is 100 days.  Figure 3 shows the geometric diagram of the model. The model is divided into three layers, in which the coal seam is located in the two rock layers, the thickness is 4.22 m, and the upper part is loaded with a stable load of 11.7 MPa. The borehole is located in the middle of the model and has a diameter of 0.089 m and a negative pressure of 13 kPa. The initial gas pressure in the coal seam is 0.9 MPa, and the gas only migrates in the coal seam. The initial time is t = 0 day, and the simulation time is 100 days.
According to the field measurement, the drainage radius of the drilling borehole down the seam is the effective extraction radius of 1.5-1.8 m when the gas is extracted for 30 days. The measured drainage radius of the mine is the effective extraction radius of 3 m when the gas is extracted for 80 days. These data are used to verify the practical application value of the model.  Table 1 shows the relevant parameters of the fluid-solid coupling model of gas-bearing coal.  According to the field measurement, the drainage radius of the drilling borehole down the seam is the effective extraction radius of 1.5-1.8 m when the gas is extracted for 30 days. The measured drainage radius of the mine is the effective extraction radius of 3 m when the gas is extracted for 80 days. These data are used to verify the practical application value of the model. Table 1 shows the relevant parameters of the fluid-solid coupling model of gas-bearing coal.

Input of the Mathematical Model
The mathematical model in Chapter 2 is inputted into the COMSOL Multiphysics numerical simulation software (COMSOL Multiphysics 5.4.0.388, COMSOL Inc., Stockholm, Sweden, 1986) to verify the rationality of the mathematical model. Figure 4 shows the specific input.

Input of the Mathematical Model
The mathematical model in Chapter 2 is inputted into the COMSOL Multiphysics numerical simulation software (COMSOL Multiphysics 5.4.0.388, COMSOL Inc., Stockholm, Sweden, 1986) to verify the rationality of the mathematical model. Figure 4 shows the specific input.

Initial Conditions and Boundary Conditions
(1) Initial conditions of gas seepage flow field

Initial Conditions and Boundary Conditions
(1) Initial conditions of gas seepage flow field where P 0 is the initial gas pressure in the coal seam.
(2) Boundary conditions When t = 0, the displacement boundary conditions of the model are: where u 0 is the initial displacement on the boundary.
(3) Stress boundary conditions where f i is the surface force on the boundary.
where q m is the gas flow on the boundary. The roof and floor of the coal seam in the model are rock layers with poor gas permeability. Thus, both are assumed to be flow boundaries, and the flux is zero.

(4) Pressure boundary conditions
where P m is the gas pressure on the boundary. Figure 5 shows the distribution of instantaneous gas pressure during the fluid-solid coupling deformation of coal. Figure 6 shows the contour figure of gas pressure for 30 and 80 days of extraction. The graph indicates that the maximum gas pressure in the coal seam is distributed on both sides of the geometric model. When t = 1 day, the gas pressure is 0.6 MPa at 0.32 m away from the drilling hole. With the increase in extraction time, the gas pressure around the borehole decreases gradually. When t = 30 days, the gas pressure is 0.6 MPa at 1.73 m from the drilling hole. When t = 80 days, the gas pressure is 0.6 MPa at 2.76 m away from the drilling hole.

Analysis of the Numerical Simulation Results of the Model
According to Henan's regulations on the prevention and control of coal and gas, the critical value of gas pressure should not be more than 0.6 MPa. This value can be used as a reference for the effective extraction radius of extraction boreholes. The numerical simulation results show that when t = 30 days, the effective extraction radius of the extraction borehole is 1.73 m. When t = 80 days, the effective extraction radius of the extraction borehole is 2.76 m. The drilling drainage radius measured by Henan University of Science and Technology is 1.5-1.8 m for 30 days and 3 m for 80 days [34,35]. The numerical simulation results are basically consistent with the measured results, indicating their strong practical application value. Figure 7 shows the instantaneous gas pressure evolution curve of the fluid-solid coupling of coal during deformation. The graph shows that the closer the distance from the drainage hole, the faster the gas pressure drops and the more noticeable the pressure relief effect is. The pressure relief of coalbed methane is a nonlinear process. Within a certain limit, the change in gas pressure gradient decreases with the increase in time. Figure 5 shows the distribution of instantaneous gas pressure during the fluid-solid coupling deformation of coal. Figure 6 shows the contour figure of gas pressure for 30 and 80 days of extraction. The graph indicates that the maximum gas pressure in the coal seam is distributed on both sides of the geometric model. When t = 1 day, the gas pressure is 0.6 MPa at 0.32 m away from the drilling hole. With the increase in extraction time, the gas pressure around the borehole decreases gradually. When t = 30 days, the gas pressure is 0.6 MPa at 1.73 m from the drilling hole. When t = 80 days, the gas pressure is 0.6 MPa at 2.76 m away from the drilling hole.  According to Henan's regulations on the prevention and control of coal and gas, the critical value of gas pressure should not be more than 0.6 MPa. This value can be used as a reference for the effective extraction radius of extraction boreholes. The numerical simulation results show that when t = 30 days, the effective extraction radius of the extraction borehole is 1.73 m. When t = 80 days, the effective extraction radius of the extraction borehole is 2.76 m. The drilling drainage radius measured by Henan University of Science and Technology is 1.5-1.8 m for 30 days and 3 m for 80 days [34,35]. The numerical simulation results are basically consistent with the measured results, indicating their strong practical application value. Figure 7 shows the instantaneous gas pressure evolution curve of the fluid-solid coupling of coal during deformation. The graph shows that the closer the distance from the drainage hole, the faster the gas pressure drops and the more noticeable the pressure relief effect is. The pressure relief of coalbed methane is a nonlinear process. Within a certain limit, the change in gas pressure gradient decreases with the increase in time. Figures 8 and 9 show the evolution curves of the porosity and permeability of coal around the borehole after drilling along the seam, respectively. The graph shows that the minimum values of porosity and permeability is lower than the initial ones at the initial state on day 0 because of the stress concentration area around the drill hole. In this area, the pores are compressed, the pore channel of gas migration and production is small, and the permeability is reduced. With the increase in extraction time, the closer the distance from the drilling hole, and the larger the values of porosity and permeability of coal are. Coal porosity and permeability increase with the increase in time, but the rate of increase declines gradually. The changing trend of permeability is almost similar to that of porosity because the closer the distance from the borehole, the more evident the coal disruption by artificial disturbance is. The coal rock breaks and forms a new pore, and its permeability increases. After gas drainage, the coal gas pressure near the extraction borehole, the gas content, and the gas adsorbed by coal particles are reduced. The coal body shrinks, coal pores and fissures are developed, and the permeability increases. The shrinkage and deformation of the coal-rock matrix are also the main factors determining coal adsorption characteristics. Although the regional pressure of the coal seam cannot be reduced, the gas dissolution in coal induces the shrinkage of the coal matrix, which enlarges the fissures and creates internal ones in the coal seam. In this way, coal porosity, the channel of gas migration and output, and permeability increase.   Figures 8 and 9 show the evolution curves of the porosity and permeability of coal around the borehole after drilling along the seam, respectively. The graph shows that the minimum values of porosity and permeability is lower than the initial ones at the initial state on day 0 because of the stress concentration area around the drill hole. In this area, the pores are compressed, the pore channel of gas migration and production is small, and the permeability is reduced. With the increase in extraction time, the closer the distance from the drilling hole, and the larger the values of porosity and permeability of coal are. Coal porosity and permeability increase with the increase in time, but the rate of increase declines gradually. The changing trend of permeability is almost similar to that of porosity because the closer the distance from the borehole, the more evident the coal disruption by artificial disturbance is. The coal rock breaks and forms a new pore, and its permeability increases. After gas drainage, the coal gas pressure near the extraction borehole, the gas content, and the gas adsorbed by coal particles are reduced. The coal body shrinks, coal pores and fissures are developed, and the permeability increases. The shrinkage and deformation of the coal-rock matrix are also the main factors determining coal adsorption characteristics. Although the regional pressure of the coal seam cannot be reduced, the gas dissolution in coal induces the shrinkage of the coal matrix, which enlarges the fissures and creates internal ones in the coal seam. In this way, coal porosity, the channel of gas migration and output, and permeability increase.    Figures 8 and 9 show that the porosity and permeability of coal also increase from both sides of the coal seam to the direction of the extraction boreholes, respectively. The closer to the extraction boreholes, the greater the increase in range. With the passage of time, the increase rate of permeability and porosity decreases significantly. This situation is mainly due to the artificial disturbance around the drilling hole. The coal becomes unstable and pressure-relieved, and pores and cracks increase. With the decrease in gas content in the coal seam pore, the effective stress of coal increases, the coal compresses, and porosity and permeability decrease. At the same time, the absorbed gas is constantly used to supply pore gas, and the volume shrinkage of coal particles increases the porosity, which affects the increase in coal permeability. With the continuation of extraction time, the increase in porosity and permeability of coal decreases under the two effects.  Figures 8 and 9 show that the porosity and permeability of coal also increase from both sides of the coal seam to the direction of the extraction boreholes, respectively. The closer to the extraction boreholes, the greater the increase in range. With the passage of time, the increase rate of permeability and porosity decreases significantly. This situation is mainly due to the artificial disturbance around the drilling hole. The coal becomes unstable and pressure-relieved, and pores and cracks increase. With the decrease in gas content in the coal seam pore, the effective stress of coal increases, the coal compresses, and porosity and permeability decrease. At the same time, the absorbed gas is constantly used to supply pore gas, and the volume shrinkage of coal particles increases the porosity, which affects the increase in coal permeability. With the continuation of extraction time, the increase in porosity and permeability of coal decreases under the two effects.

Discussion
After establishing the fluid-solid coupling model of gas-bearing coal, the model's field practicability needs to be verified. A good model should be applied to practical applications to reflect its importance and guiding significance for practical design and production work. The SF6 gas tracing method is mainly used to measure the effective extraction radius of the bedding borehole in this coal mine. This method overcomes many of the limitations of traditional measuring methods, such as numerous drilling holes, complicated working procedures, high requirement for sealing quality, long time, and large deviation, leading to a highly accurate, reliable, and simple measurement [36]. The measured effective extraction radius is 1.5-1.8 m at 30 days of extraction. The relevant parameters of the measured coal seams are inputted into the mathematical model, and the coupled numerical solution is calculated with COMSOL Multiphysics simulation software. The calculation results show that the effective extraction radius is 1.73 m. Afterward, the effective extraction radius is measured independently in the mine. The measured effective extraction radius is 3 m at 80 days, and the effective extraction radius is 2.76 m by numerical simulation.
A comparison with the results of field measurement shows that the numerical simulation findings of this model are close to actual conditions. The results provide reference values for the rational layout of extraction boreholes in coal mines. Based on the fluid-solid coupling equation, the

Discussion
After establishing the fluid-solid coupling model of gas-bearing coal, the model's field practicability needs to be verified. A good model should be applied to practical applications to reflect its importance and guiding significance for practical design and production work. The SF6 gas tracing method is mainly used to measure the effective extraction radius of the bedding borehole in this coal mine. This method overcomes many of the limitations of traditional measuring methods, such as numerous drilling holes, complicated working procedures, high requirement for sealing quality, long time, and large deviation, leading to a highly accurate, reliable, and simple measurement [36]. The measured effective extraction radius is 1.5-1.8 m at 30 days of extraction. The relevant parameters of the measured coal seams are inputted into the mathematical model, and the coupled numerical solution is calculated with COMSOL Multiphysics simulation software. The calculation results show that the effective extraction radius is 1.73 m. Afterward, the effective extraction radius is measured independently in the mine. The measured effective extraction radius is 3 m at 80 days, and the effective extraction radius is 2.76 m by numerical simulation.
A comparison with the results of field measurement shows that the numerical simulation findings of this model are close to actual conditions. The results provide reference values for the rational layout of extraction boreholes in coal mines. Based on the fluid-solid coupling equation, the effective extraction radius simulation of drilling boreholes along the seam reveals the gas pressure distribution state at different extraction times. This condition eliminates the blank belt of gas drainage with the time effect, changes the current disadvantage of the restricted effective radius test of gas extraction by many factors, and ensures the safe production of mines to a great extent. Moreover, the process of establishing the model reveals the mechanism of gas migration in the coal seam and describes factors one by one from the perspectives of coal pore, coal particle adsorption expansion deformation, compression deformation, stress field, and gas seepage flow field. Then, a relatively complete practical mathematical model is obtained. This model considers the effects of adsorption expansion and the Klinkenberg effect on gas migration in coal seams. The model improves the application of fluid-solid coupling theory in coal. The numerical simulation results show the dynamic changes in gas pressure in the coal seam and the dynamic evolution of porosity and permeability. The model can simulate the dynamic evolution law of coalbed methane in the corresponding coal seam of the coal mine based on different coal seam parameters, such as the influence of different types of water and ash on gas migration and various extraction negative pressures on the drainage effect. The numerical simulation shows the extraction effect and scope of gas drainage boreholes and provides important theoretical support and basis for the rational optimization and layout of gas drainage boreholes in mines.

Conclusions
The mathematical model combines the findings of previous studies on the effects of adsorption expansion and the Klinkenberg effect on gas migration in coal seams. On the basis of the basic definitions of porosity and permeability, a mathematical model of dynamic evolution of porosity and permeability is derived. On the basis of fluid-solid coupling theory, a mathematical model of fluid-solid coupling with gas-bearing coal is established.
According to the results of the numerical simulation, when t = 30 days and t = 80 days, the effective extraction radius of the bedding borehole reaches 1.73 m and 2.76 m, respectively. The simulation results are consistent with the actual measured extraction radius values. Figure 5 shows the effect of gas pressure reduction around the borehole with the continuous change in simulation time. The gas pressure decreases, and the porosity and permeability of coal increase with the increase in gas extraction time. In addition, the growth rate of permeability and porosity decreases with the increase in gas extraction time. These results are consistent with the field permeability test law and can be used as reference to further understand the mechanism of gas extraction and mine gas control. The research results also have theoretical significance and practical application value.
The dynamic evolution mathematical model of fluid-solid coupling for gas bearing coal can reflect the coalbed methane migration in a mining area of the mine through the coal seam parameters measured by coal miners. According to the simulation results of the model, the effective extraction radius of the borehole can be predicted. Thus, the extraction borehole can be optimized and reasonably arranged for safety reasons and scientific purposes, to effectively control mine gas, and to provide strong support for decision-makers as they formulate efficient coal mining schemes.
Although the model is an extension of the theory of fluid-solid coupling, its simulation results can be suitable for field applications. However, several problems should be considered in the multi-field coupling model, thus indicating the need for further research and improvement. For example, with the increase in mining depth, the influence of temperature on gas adsorption and migration in the coal seam considerably influences gas extraction. Another example is the influence of movable and residual water in the coal seam on expansion stress in the gas-bearing coal seam [37]. Studying these problems is crucial for the further improvement of the model.