# Numerical Investigation of Complex Thermal Coal-Gas Interactions in Coal-Gas Migration

Authors: Xiaoyan Ni; Peng Gong; Yi Xue
Journal: Advances in Civil Engineering
 (2018)

---

## Abstract

Understanding the influence of temperature on the gas seepage of coal seams is helpful to achieve the efficient extraction of underground coal seam gas. Thermal coal-gas interactions involve a series of complex interactions between gas and solid coal. Although the interactions between coal and gas have been studied thoroughly, few studies have considered the temperature evolution characteristics of coal seam gas extraction under the condition of variable temperature because of the complexity of the temperature effect on gas drainage. In this study, the fully coupled transient model combines the relationship of gas flow, heat transfer, coal mass deformation, and gas migration under variable temperature conditions and represents an important nonlinear response to gas migration caused by the change of effective stress. Then, this complex model is implemented into a finite element (FE) model and solved through the numerical method. Its reliability was verified by comparing with historical data. Finally, the effect of temperature on coal permeability and gas pressure is studied. The results reveal that the gas pressure in coal fracture is generally higher than that in the matrix blocks. The higher temperature of the coal seam induces the faster increase of the gas pressure. Temperature has a great effect on the gas seepage behavior in the coal seams.


---

## Body

## 1. Introduction

Coal seam gas is an important natural source and an important part in the clean energy structure in China [1]. However, most Chinese coalbed methane cannot meet the requirements of effective recovery because of the low permeability characteristic [2, 3]. Therefore, many technological means are adopted to enhance the coal permeability and develop the gas extraction efficiency [4, 5], such as, hydraulic fracturing technology [6]. Besides, hot injection has been adopted to enhance the extraction of coalbed methane [7, 8]. Therefore, the thermal evolution characteristic evolution during the extraction of coalbed methane needs to be studied.The thermal stimulation of gas reservoir triggers the complex interaction between coal, gas, and temperature. It changes the deformation behavior of coal, as well as the heat and gas flow. Temperature has an obvious effect on the adsorption capacity of coal seam gas, which is negatively correlated with the increase of temperature [9–11]. At the same time, the thermal expansion coefficient is heterogeneous and anisotropic due to the different response of different components to temperature changes. Noorishad et al. [12] evaluated the coupled thermal-hydraulic-mechanical characteristics of rock and investigated the effects of thermal stresses on permeability through the deformation change of the fractures. Mctigue [13] developed a fluid-saturated, porous, thermoelastic model and calculated the heating of a half space at constant temperature conditions. Harpalani and Schraufnagel [14] found that desorption induced the shrinking of the coal matrix and increased permeability through seepage tests. Zhou et al. [15] consider thermodynamically coupled water and heat flow and demonstrate the influence of thermodynamic coupling through numerical and analytical solutions. Zhu et al. [16] analyzed the influence of temperature on coal permeability based on a complex thermal-hydrological-mechanical model and found that sorption-induced volumetric strain affected the seepage behavior of coal seam significantly. Cai et al. [17] found that thermal treatment can promote methane desorption and increase the permeability through the numerical simulation. Qu et al. [18] developed a new thermal model with gas flow and matrix adsorption and suggested that the increase in temperature leads to the decrease in coal swelling and larger cleat aperture and higher coal permeability.The interaction mechanism of heat injection has not yet fully understood for the design of engineering. It is indispensable to study the coal-gas interaction under variable temperature. In this study, a coupled thermal-hydrological-mechanical model is established, considering the sorption characteristic and permeability evolution. Then, this complex model is implemented into a finite element (FE) and solved through the numerical method. Finally, the effect of temperature on coal permeability and gas pressure is studied based on this coupled model.


## 2. Equations of Coupled Model

### 2.1. Gas Seepage Characteristic in Coal Fractures

For the coal seam reservoir, the non-Darcy flow can be expressed as(1)−∇pf=υkgμ⟶+βρgμ⟶μ⟶,where μ⟶ is velocity vector, ρg is density, β=1.75/150kgϕ3, ϕ is porosity, and kg is the permeability of coal.The above equation can be expressed in the following form:(2)μ⟶=−kg1+kg/υβρgμ⟶υ∇pf=−kgfqiυ∇pf,where fqi=1+kg/υβρgμ⟶ is the Forchheimer coefficient.The mass conservation equation of gas migration can be expressed as(3)∂∂tϕfρg+∇ρg⋅μ⟶=Qs1−ϕf,where μ⟶ is the velocity vector and Qs is the gas source, and the mass content can be expressed as [19](4)m=ρgϕf+ρgaρcVsg,where ϕf is the porosity, ρga is the gas density, ρc is the coal density, and Vsg is the absorbed content.In the coal seams, the gas absorption volume can be expressed as(5)Vsg=VLpfpf+PLexp−c21+c1pfTar+T−Tt.The Sorption Strainεs can be expressed as(6)εs=αsgVsg,where Vsg is the content of absorbed gas and αsg is the strain coefficient.The gas density can be obtained as(7)ρg=MgRTar+Tp.The mass conservation equation can be expressed as(8)ρgapa∂ϕfpf∂t+∇−kgμρgfqi∇pf=Qs1−ϕf.


### 2.2. Gas Diffusion in Coal Matrix

The gas migration process experiences three substeps: flow in fracture, gas diffusion, and sorption in matrix [20]. Figure 1 shows a conceptual model for gas transport. The source term can be expressed as(9)Qs=Dσccm−cf,∂mm∂t=−McτRTpm−pf,where Qs is the exchange from matrix to fractures, cm is the gas concentration in matrix, cf is the gas concentration, σc is the shape factor, and τ is the sorption time.Figure 1
A conceptual model for gas storage and transport in the coal seams. (a) Natural fracture networks. (b) Stage 1: desorption from internal surfaces. (c) Stage 2: flow/diffusion in the matrix. (d) Stage 3: flow in the natural fracture network.
(a)
(b)
(c)
(d)The relationship of gas concentration and gas pressure is written as(10)cm=McRTpm,cf=McRTpf,where Mc is the molar mass.The diffusion equation can be written as(11)dmbdt=−1τmb−me,where me is the equilibrium gas content.Then, the diffusion equation can be expressed as(12)∂mm∂t=−McτRTpm−pf.The permeabilitykg can be written as [21](13)kg=k∞1+bpf,where b is the Klinkenberg coefficient, b=αkk∞−0.36, where αk is the Klinkenberg effect coefficient, αk = 0.251.


### 2.3. Solid Mechanics Equation

The stress equation of coal seam can be written as(14)σeij=σij−αpfδij.The stress-strain equation can be written as(15)εij=12ui,j+uj,i.The overall stress balance equation can be written as(16)Gui,jj+G1−2μuj,ji−αpf,i−KαTTi−Kεs,i+fi=0.


### 2.4. Coal Permeability

The relationship between porosity and effective stress can be expressed as [22](17)Δϕf=1Kβf−ϕfσ¯+pf.Then, the porosity is expressed as(18)ϕf=α−α−ϕ0exp−1Kσ¯−σ¯0+pf−pf0.By substituting the porosity, it can be rewritten as(19)ϕf=α−α−ϕf0exp−εV+pfKs−εs−αTT−εV0+pf0Ks−εs0−αTT0,where S=εV+pf/Ks−εs and S0=εV0+pf0/Ks−εs0. p0 is the initial pressure, and ϕ0 is the initial porosity.The evolution characteristics of permeability can be expressed as(20)kgk∞0=k∞k∞01+bpf=ϕϕ031+bpf.


### 2.5. Energy Conservation

Total heat fluxqT is expressed as(21)qT=−λM∇T+ρgCgqgTar+T,where qT is the thermal flux and ρs is the mass density.The thermal balance can be written as [23](22)∂ρCMTar+T∂t+Tar+TKgαg∇⋅kgμ∇p+Tar+TKαT∂εV∂t=−∇⋅qT,where ρCM=ϕfρgCg+1−ϕfρsCs and ρs is the mass density.The conservation of mass yields(23)∂1−ϕfρs∂t=0,∂ϕfρg∂t=−∇⋅ρgqg.Considering1−ϕλs≫ϕλs and λM≈1−ϕfλs≈λs, it yields [16](24)ρCM∂T∂t−Tar+TKgαg∇⋅kgυ∇pf+Tar+TKαT∂εV∂t=λM∇2T+ρgapfTaCgpaTar+Tkgυ∇pf∇T.These equations (8), (12), (16), and (24) describe the fully coupling model of coal seam gas migration, including mechanical deformation of coal (gas desorption induces coal shrinkage and self-heat induces coal expansion), gas diffusion from the matrix, gas flow, and heat transfer in fractures. The coupled relationship is illustrated in Figure 2. The THM model indicates the nonlinear response of gas migration in underground coal seams. It is hard to get an analytic solution of these equations, so they are achieved through the numerical method with Comsol Multiphysics (a powerful partial differential equation solver) and appropriate boundary conditions.Figure 2
Couplings in the gas migration process.


## 2.1. Gas Seepage Characteristic in Coal Fractures

For the coal seam reservoir, the non-Darcy flow can be expressed as(1)−∇pf=υkgμ⟶+βρgμ⟶μ⟶,where μ⟶ is velocity vector, ρg is density, β=1.75/150kgϕ3, ϕ is porosity, and kg is the permeability of coal.The above equation can be expressed in the following form:(2)μ⟶=−kg1+kg/υβρgμ⟶υ∇pf=−kgfqiυ∇pf,where fqi=1+kg/υβρgμ⟶ is the Forchheimer coefficient.The mass conservation equation of gas migration can be expressed as(3)∂∂tϕfρg+∇ρg⋅μ⟶=Qs1−ϕf,where μ⟶ is the velocity vector and Qs is the gas source, and the mass content can be expressed as [19](4)m=ρgϕf+ρgaρcVsg,where ϕf is the porosity, ρga is the gas density, ρc is the coal density, and Vsg is the absorbed content.In the coal seams, the gas absorption volume can be expressed as(5)Vsg=VLpfpf+PLexp−c21+c1pfTar+T−Tt.The Sorption Strainεs can be expressed as(6)εs=αsgVsg,where Vsg is the content of absorbed gas and αsg is the strain coefficient.The gas density can be obtained as(7)ρg=MgRTar+Tp.The mass conservation equation can be expressed as(8)ρgapa∂ϕfpf∂t+∇−kgμρgfqi∇pf=Qs1−ϕf.


## 2.2. Gas Diffusion in Coal Matrix

The gas migration process experiences three substeps: flow in fracture, gas diffusion, and sorption in matrix [20]. Figure 1 shows a conceptual model for gas transport. The source term can be expressed as(9)Qs=Dσccm−cf,∂mm∂t=−McτRTpm−pf,where Qs is the exchange from matrix to fractures, cm is the gas concentration in matrix, cf is the gas concentration, σc is the shape factor, and τ is the sorption time.Figure 1
A conceptual model for gas storage and transport in the coal seams. (a) Natural fracture networks. (b) Stage 1: desorption from internal surfaces. (c) Stage 2: flow/diffusion in the matrix. (d) Stage 3: flow in the natural fracture network.
(a)
(b)
(c)
(d)The relationship of gas concentration and gas pressure is written as(10)cm=McRTpm,cf=McRTpf,where Mc is the molar mass.The diffusion equation can be written as(11)dmbdt=−1τmb−me,where me is the equilibrium gas content.Then, the diffusion equation can be expressed as(12)∂mm∂t=−McτRTpm−pf.The permeabilitykg can be written as [21](13)kg=k∞1+bpf,where b is the Klinkenberg coefficient, b=αkk∞−0.36, where αk is the Klinkenberg effect coefficient, αk = 0.251.


## 2.3. Solid Mechanics Equation

The stress equation of coal seam can be written as(14)σeij=σij−αpfδij.The stress-strain equation can be written as(15)εij=12ui,j+uj,i.The overall stress balance equation can be written as(16)Gui,jj+G1−2μuj,ji−αpf,i−KαTTi−Kεs,i+fi=0.


## 2.4. Coal Permeability

The relationship between porosity and effective stress can be expressed as [22](17)Δϕf=1Kβf−ϕfσ¯+pf.Then, the porosity is expressed as(18)ϕf=α−α−ϕ0exp−1Kσ¯−σ¯0+pf−pf0.By substituting the porosity, it can be rewritten as(19)ϕf=α−α−ϕf0exp−εV+pfKs−εs−αTT−εV0+pf0Ks−εs0−αTT0,where S=εV+pf/Ks−εs and S0=εV0+pf0/Ks−εs0. p0 is the initial pressure, and ϕ0 is the initial porosity.The evolution characteristics of permeability can be expressed as(20)kgk∞0=k∞k∞01+bpf=ϕϕ031+bpf.


## 2.5. Energy Conservation

Total heat fluxqT is expressed as(21)qT=−λM∇T+ρgCgqgTar+T,where qT is the thermal flux and ρs is the mass density.The thermal balance can be written as [23](22)∂ρCMTar+T∂t+Tar+TKgαg∇⋅kgμ∇p+Tar+TKαT∂εV∂t=−∇⋅qT,where ρCM=ϕfρgCg+1−ϕfρsCs and ρs is the mass density.The conservation of mass yields(23)∂1−ϕfρs∂t=0,∂ϕfρg∂t=−∇⋅ρgqg.Considering1−ϕλs≫ϕλs and λM≈1−ϕfλs≈λs, it yields [16](24)ρCM∂T∂t−Tar+TKgαg∇⋅kgυ∇pf+Tar+TKαT∂εV∂t=λM∇2T+ρgapfTaCgpaTar+Tkgυ∇pf∇T.These equations (8), (12), (16), and (24) describe the fully coupling model of coal seam gas migration, including mechanical deformation of coal (gas desorption induces coal shrinkage and self-heat induces coal expansion), gas diffusion from the matrix, gas flow, and heat transfer in fractures. The coupled relationship is illustrated in Figure 2. The THM model indicates the nonlinear response of gas migration in underground coal seams. It is hard to get an analytic solution of these equations, so they are achieved through the numerical method with Comsol Multiphysics (a powerful partial differential equation solver) and appropriate boundary conditions.Figure 2
Couplings in the gas migration process.


## 3. Model Establishment and Analysis

### 3.1. Verification with Field Test

To verify the validity of the model established in this paper in the calculation of simulated extraction, the numerical simulation is carried out according to the practical example of coal seam extraction. The actual size of the coal seam is 568 × 568 m2, and the parameters are obtained from Mora and Wattenbarger [24]. The numerical results and the field production data are compared in Figure 3. Good agreement was achieved between the results of numerical simulation and the field production data, which proves the validity of the model.Figure 3
Comparison between numerical result and field test.


### 3.2. Model Establishment

A model is established to analyze the effect of adsorption on the distribution of coal permeability and gas pressure. The size of the model is 0.1 m × 0.1 m, as shown in Figure4. The upper boundary is free. And the condition of the upper boundary is the constant gas pressure 2.5 MPa and the constant temperature 333 K. The initial pressure of the coal seam is 0.1 MPa and the parameters in the calculation are listed in Table 1.Figure 4
Calculation model.Table 1
Parameters of the model.
Parameter
Value
Sources
Parameter
Value
Sources
Young’s modulus of coal,E (MPa)
2713
[22]
CH4 Langmuir volume constant, VL (m3/kg)
0.043
[22]
Young’s modulus of the coal grains,Es (MPa)
4070
[22]
Specific heat capacity of gas,Cg (J/kg·K)
1.625 × 103
[16]
Initial porosity,ϕ0 (—)
0.01
[22]
Specific heat capacity of coal,Cs (J/kg·K)
1.25 × 103
[16]
Density of coal,ρc (kg/m3)
1250
[22]
Pressure coefficient,c1 (MPa−1)
0.07
[16]
Poisson’s ratio of coal,ν (—)
0.34
[22]
Sorption strain coefficient,αsg (kg/m3)
0.06
[18]
Initial gas permeability,k∞0 (m2)
1.09 × 10−18
[22]
Volumetric coefficient of matrix,αT (K−1)
2.4 × 10−5
[18]
Density of CH4 at standard condition, ρga (kg/m3)
0.717
[16]
Klinkenberg effect,b (Pa)
1.44 × 105
[18]
Gas dynamic viscosity,υ (N·s/m2)
1.84 × 10−5
[16]
Temperature coefficient,c2 (K−1)
0.02
[16]
CH4 Langmuir pressure constant, PL (MPa)
1.57
[16]
Thermal conductivity of coal,λs (J/(m·s·K))
0.2
[16]


### 3.3. Effect of Adsorption on Gas Pressure

There are two different pressures between matrix blocks and coal fractures. Because of the high permeability in coal fractures, gas can quickly flow through the fractures. The gas pressure in fractures is generally higher than that in the matrix blocks in the calculation model. The matrix also continues to adsorb gas, and the pressure rises gradually. As the gas is injected from the top of the model, the crack pressure reaches 2.5 MPa firstly at the upper boundary, and the matrix permeability is also increasing. Finally, pore pressure and matrix pressure reach the equilibrium state. Figure5 shows the gas pressure distribution at different times. As the time goes on, the difference of pressure becomes smaller and smaller. The initial pressure in the coal seam is 0.1 MPa. The gas pressure maximums of coal matrix are 0.3 MPa, 1.8 MPa, and 2.3 MPa at 1 h, 4 h, and 6 h, respectively. Finally, two pressures reach 2.5 MPa at about 24 hours.Figure 5
Gas pressure change law of the model.Figure6 shows the gas pressure distribution at the monitoring line. The difference between the matrix pressure and the fracture pressure can be seen more easily from the detection line. The gas pressure at the bottom of the model is low, while the gas pressure at the entrance is high. Within the equilibrium time of 1h, 4h, 8h, 12h, and 24h, the minimum pressure of the coal fracture is 0.1 MPa, 1.42 MPa, 2.16 MPa, 2.39 MPa, and 2.5 MPa, respectively; meanwhile, the minimum pressure of the coal matrix is 0.1 MPa, 0.64 MPa, 1.66 MPa, 2.19 MPa, and 2.5 MPa. With the passage of time, the pressure of internal coal mass increases gradually. The ratios of minimum matrix pressure and the minimum fracture pressure are 45%, 76%, and 92% at 4 hours, 8 hours, and 12 hours, respectively.Figure 6
Gas pressure distribution along the monitoring line.


### 3.4. Effect of Adsorption on Permeability Distribution

Based on the dual-porous medium model, the fracture pressure and the matrix pressure are calculated through different equations. Therefore, we can further understand the dual-porous medium model by analyzing the evolution rule of permeability. Figure7 is the variation diagram of coal permeability and gas content of coal mass at different time. The permeability of coal fractures decreases gradually. It is caused by the adsorption swelling of the matrix and the reduction of the crack opening in the coal mass.Figure 7
Permeability distribution along the monitoring line.Coal mass expends and the fracture aperture decreases when it absorbs gas. Meanwhile, the coal mass is limited by the boundary conditions, so coal deforms under the combined influence of the change of gas pressure, expansion stress, and effective stress [25]. Figure 7 shows that the permeability decreases with time, and the volumetric strain caused by adsorption is the main reason. There is an inverse relationship between the permeability and the content of gas; that is, the gas content of coal increases ceaselessly and permeability decreases gradually. The maximum of permeability ratio decreased from 0.901 to 0.876, while the minimum of gas content increased from 0.00239 m3/kg to 0.0481 m3/kg when the time increased from 1 h to 8 h. As time goes, the permeability of coal gradually stabilizes and tends to be at a uniform state.The changes of permeability and porosity of a point are analyzed. The point is at the center of the model, and its coordinates are (0.05 m, 0.05 m). Figure8 is the figure of permeability, porosity, gas content, and temperature change at the monitoring point A. Figure 8 shows that the evolution of permeability and gas content at point A are consistent with Figure 7. The temperature and gas content increase with the increase of time. The permeability and porosity of coal decrease rapidly before 150 h, which is mainly caused by the large pressure difference. Pore pressure and adsorption expansion affect together, which causes the change of permeability and porosity. With the decrease of pressure difference in coal mass, the change speed of permeability and porosity decreases gradually and tends to be stable finally.Figure 8
Permeability distribution of the monitoring point A.


### 3.5. Effect of Temperature on Permeability Distribution

Figure9 gives the temperature change law of the model. The temperature at the entrance is relatively high, and the temperature in the coal seam keeps rising. The temperature at point A increases from 303.8 K at 1 h to 308.6 K at 8 h, finally reaching 322.7 K at 40 h. The deformation characteristics of coal seam are given in Figure 10. After the coalbed is heated up, the adsorption increases and the volume expands. The bottom, left, and right boundaries of coal seam are restricted, so there is no normal displacement in this boundaries. When the gas pressure rises on the upper boundary of the coal mass, the coal mass absorbs gas constantly, which leads to expansive deformation of coal mass and expanding outwards of the upper boundary. These characteristics also show that the volumetric strain induced by gas adsorption plays a key role. If the effect of gas adsorption is not considered, the effective stress will decrease, and the permeability of coal may increase. If the influence of gas adsorption is ignored, pore pressure may be underestimated.Figure 9
Temperature change law of the model.Figure 10
Evolution of the coal volume with time under the effect of temperature.The distributions of gas pressure and permeability under different temperatures are given in Figures11 and 12. It indicates that temperature has a greater impact on the gas migration in the coal seam. The influence of temperature is reduced as the coal seam gas pressure approaches the setting pressure. High temperature can promote the desorption of gas and then change the fracture aperture and porosity of coal mass. The influence of temperature on the porosity of coal seam is a complicated process. The temperature evolution is obtained under the coupling effect of multiple factors, such as restriction on the boundary condition and adsorptive expansion deformation of coal mass. At the initial stage, the gas pressure and porosity of coal mass differ greatly at different temperatures. Gas pressure and permeability change constantly with the increase of time, and finally attain states of equilibrium. The gas pressure of point A is 1.05 MPa, 1.27 MPa, and 1.41 MPa at 1h when the initial boundary temperature is 333 K, 353 K, and 373 K, respectively.Figure 11
Gas pressure distribution with different temperatures.Figure 12
Permeability distribution with different temperatures.In order to quantify the influence of gas diffusion behavior, the fracture pressure distribution without Klinkenberg effect and different adsorption times was also plotted as shown in Figure13. The adsorption time significantly affects the distribution of fracture pressure. When the adsorption time is large, the gas pressure of coal seam will be underestimated. The influence of diffusion effect is the same as that of Klinkenberg effect at the same order of magnitude. In the sorption stage, the total volume expansion effect and Klinkenberg effect of coal increase the permeability. When the Klinkenberg effect is not considered, the adsorption pressure will be underestimated.Figure 13
Gas pressure distribution under different cases.


## 3.1. Verification with Field Test

To verify the validity of the model established in this paper in the calculation of simulated extraction, the numerical simulation is carried out according to the practical example of coal seam extraction. The actual size of the coal seam is 568 × 568 m2, and the parameters are obtained from Mora and Wattenbarger [24]. The numerical results and the field production data are compared in Figure 3. Good agreement was achieved between the results of numerical simulation and the field production data, which proves the validity of the model.Figure 3
Comparison between numerical result and field test.


## 3.2. Model Establishment

A model is established to analyze the effect of adsorption on the distribution of coal permeability and gas pressure. The size of the model is 0.1 m × 0.1 m, as shown in Figure4. The upper boundary is free. And the condition of the upper boundary is the constant gas pressure 2.5 MPa and the constant temperature 333 K. The initial pressure of the coal seam is 0.1 MPa and the parameters in the calculation are listed in Table 1.Figure 4
Calculation model.Table 1
Parameters of the model.
Parameter
Value
Sources
Parameter
Value
Sources
Young’s modulus of coal,E (MPa)
2713
[22]
CH4 Langmuir volume constant, VL (m3/kg)
0.043
[22]
Young’s modulus of the coal grains,Es (MPa)
4070
[22]
Specific heat capacity of gas,Cg (J/kg·K)
1.625 × 103
[16]
Initial porosity,ϕ0 (—)
0.01
[22]
Specific heat capacity of coal,Cs (J/kg·K)
1.25 × 103
[16]
Density of coal,ρc (kg/m3)
1250
[22]
Pressure coefficient,c1 (MPa−1)
0.07
[16]
Poisson’s ratio of coal,ν (—)
0.34
[22]
Sorption strain coefficient,αsg (kg/m3)
0.06
[18]
Initial gas permeability,k∞0 (m2)
1.09 × 10−18
[22]
Volumetric coefficient of matrix,αT (K−1)
2.4 × 10−5
[18]
Density of CH4 at standard condition, ρga (kg/m3)
0.717
[16]
Klinkenberg effect,b (Pa)
1.44 × 105
[18]
Gas dynamic viscosity,υ (N·s/m2)
1.84 × 10−5
[16]
Temperature coefficient,c2 (K−1)
0.02
[16]
CH4 Langmuir pressure constant, PL (MPa)
1.57
[16]
Thermal conductivity of coal,λs (J/(m·s·K))
0.2
[16]


## 3.3. Effect of Adsorption on Gas Pressure

There are two different pressures between matrix blocks and coal fractures. Because of the high permeability in coal fractures, gas can quickly flow through the fractures. The gas pressure in fractures is generally higher than that in the matrix blocks in the calculation model. The matrix also continues to adsorb gas, and the pressure rises gradually. As the gas is injected from the top of the model, the crack pressure reaches 2.5 MPa firstly at the upper boundary, and the matrix permeability is also increasing. Finally, pore pressure and matrix pressure reach the equilibrium state. Figure5 shows the gas pressure distribution at different times. As the time goes on, the difference of pressure becomes smaller and smaller. The initial pressure in the coal seam is 0.1 MPa. The gas pressure maximums of coal matrix are 0.3 MPa, 1.8 MPa, and 2.3 MPa at 1 h, 4 h, and 6 h, respectively. Finally, two pressures reach 2.5 MPa at about 24 hours.Figure 5
Gas pressure change law of the model.Figure6 shows the gas pressure distribution at the monitoring line. The difference between the matrix pressure and the fracture pressure can be seen more easily from the detection line. The gas pressure at the bottom of the model is low, while the gas pressure at the entrance is high. Within the equilibrium time of 1h, 4h, 8h, 12h, and 24h, the minimum pressure of the coal fracture is 0.1 MPa, 1.42 MPa, 2.16 MPa, 2.39 MPa, and 2.5 MPa, respectively; meanwhile, the minimum pressure of the coal matrix is 0.1 MPa, 0.64 MPa, 1.66 MPa, 2.19 MPa, and 2.5 MPa. With the passage of time, the pressure of internal coal mass increases gradually. The ratios of minimum matrix pressure and the minimum fracture pressure are 45%, 76%, and 92% at 4 hours, 8 hours, and 12 hours, respectively.Figure 6
Gas pressure distribution along the monitoring line.


## 3.4. Effect of Adsorption on Permeability Distribution

Based on the dual-porous medium model, the fracture pressure and the matrix pressure are calculated through different equations. Therefore, we can further understand the dual-porous medium model by analyzing the evolution rule of permeability. Figure7 is the variation diagram of coal permeability and gas content of coal mass at different time. The permeability of coal fractures decreases gradually. It is caused by the adsorption swelling of the matrix and the reduction of the crack opening in the coal mass.Figure 7
Permeability distribution along the monitoring line.Coal mass expends and the fracture aperture decreases when it absorbs gas. Meanwhile, the coal mass is limited by the boundary conditions, so coal deforms under the combined influence of the change of gas pressure, expansion stress, and effective stress [25]. Figure 7 shows that the permeability decreases with time, and the volumetric strain caused by adsorption is the main reason. There is an inverse relationship between the permeability and the content of gas; that is, the gas content of coal increases ceaselessly and permeability decreases gradually. The maximum of permeability ratio decreased from 0.901 to 0.876, while the minimum of gas content increased from 0.00239 m3/kg to 0.0481 m3/kg when the time increased from 1 h to 8 h. As time goes, the permeability of coal gradually stabilizes and tends to be at a uniform state.The changes of permeability and porosity of a point are analyzed. The point is at the center of the model, and its coordinates are (0.05 m, 0.05 m). Figure8 is the figure of permeability, porosity, gas content, and temperature change at the monitoring point A. Figure 8 shows that the evolution of permeability and gas content at point A are consistent with Figure 7. The temperature and gas content increase with the increase of time. The permeability and porosity of coal decrease rapidly before 150 h, which is mainly caused by the large pressure difference. Pore pressure and adsorption expansion affect together, which causes the change of permeability and porosity. With the decrease of pressure difference in coal mass, the change speed of permeability and porosity decreases gradually and tends to be stable finally.Figure 8
Permeability distribution of the monitoring point A.


## 3.5. Effect of Temperature on Permeability Distribution

Figure9 gives the temperature change law of the model. The temperature at the entrance is relatively high, and the temperature in the coal seam keeps rising. The temperature at point A increases from 303.8 K at 1 h to 308.6 K at 8 h, finally reaching 322.7 K at 40 h. The deformation characteristics of coal seam are given in Figure 10. After the coalbed is heated up, the adsorption increases and the volume expands. The bottom, left, and right boundaries of coal seam are restricted, so there is no normal displacement in this boundaries. When the gas pressure rises on the upper boundary of the coal mass, the coal mass absorbs gas constantly, which leads to expansive deformation of coal mass and expanding outwards of the upper boundary. These characteristics also show that the volumetric strain induced by gas adsorption plays a key role. If the effect of gas adsorption is not considered, the effective stress will decrease, and the permeability of coal may increase. If the influence of gas adsorption is ignored, pore pressure may be underestimated.Figure 9
Temperature change law of the model.Figure 10
Evolution of the coal volume with time under the effect of temperature.The distributions of gas pressure and permeability under different temperatures are given in Figures11 and 12. It indicates that temperature has a greater impact on the gas migration in the coal seam. The influence of temperature is reduced as the coal seam gas pressure approaches the setting pressure. High temperature can promote the desorption of gas and then change the fracture aperture and porosity of coal mass. The influence of temperature on the porosity of coal seam is a complicated process. The temperature evolution is obtained under the coupling effect of multiple factors, such as restriction on the boundary condition and adsorptive expansion deformation of coal mass. At the initial stage, the gas pressure and porosity of coal mass differ greatly at different temperatures. Gas pressure and permeability change constantly with the increase of time, and finally attain states of equilibrium. The gas pressure of point A is 1.05 MPa, 1.27 MPa, and 1.41 MPa at 1h when the initial boundary temperature is 333 K, 353 K, and 373 K, respectively.Figure 11
Gas pressure distribution with different temperatures.Figure 12
Permeability distribution with different temperatures.In order to quantify the influence of gas diffusion behavior, the fracture pressure distribution without Klinkenberg effect and different adsorption times was also plotted as shown in Figure13. The adsorption time significantly affects the distribution of fracture pressure. When the adsorption time is large, the gas pressure of coal seam will be underestimated. The influence of diffusion effect is the same as that of Klinkenberg effect at the same order of magnitude. In the sorption stage, the total volume expansion effect and Klinkenberg effect of coal increase the permeability. When the Klinkenberg effect is not considered, the adsorption pressure will be underestimated.Figure 13
Gas pressure distribution under different cases.


## 4. Conclusions

A new model about heat transfer was developed to study the dynamic problem of gas adsorption in coal seam under variable temperature conditions. In the model, the coal mass was regarded as a dual-porous medium, and each point of the coal mass has two different pressures, i.e., pore pressure and matrix pressure. The relationship between adsorption equilibrium and permeability evolution was analyzed through the numerical method. The main results are as follows:The gas pressure in fractures is generally higher than that in the matrix blocks in the calculation model. The matrix also continues to adsorb gas, and the pressure rises gradually. Finally, pore pressure and matrix pressure reach the equilibrium state.The adsorption of coal increases with the increase of gas pressure, and the gas content increases continuously. The porosity of coal mass decreases with the increase of time, while the temperature and gas content increase with the increase of time.The higher temperature induces the faster increase of the gas pressure. Temperature has a greater impact on the gas migration. As the coal seam gas pressure approaches the setting pressure, the influence of temperature is reduced.



---

*Source: 9020872-2018-12-05.xml*