# The Migration of Coalbed Methane under Mining Pressure and Air Injection: A Case Study in China

Authors: Liqiang Zhang; Yu Wu; Hai Pu; Xiaoping He; Pan Li
Journal: Geofluids
 (2018)

---

## Abstract

Gas outburst has always affected the safety of coal mining. To eliminate this risk by high-efficiency extraction of coalbed methane (CBM) in 4102 working face of number 3 coal seam in Hebi Number 3 coal mine, a model of CBM extraction in working face was established which was considering the mining impact of adjacent 4101 working face. In this model, the coupling relationships between stress, desorption, and migration of methane were analyzed. Moreover, we also studied the changes of methane pressure, plastic failure scope, and permeability of coal during the mining and then verified the results with the field data. And on this basis, a stimulation solution for methane extraction by injecting air into coal seam was presented, and the extraction effect was simulated. The simulation results show that the injection of air decreases the effective stress of coal which increases the permeability of coal and promotes the methane migration within the coal seam fractures. Besides, affected by the velocity of gas migration, the pressure drop between fractures and matrix will reduce with time while air injection can provide extra power for gas migration in fractures which causes the desorption and diffusion of methane in the matrix. So this stimulation solution can enhance the efficiency of gas extraction of coal seam and prevent gas outburst of the working face.


---

## Body

## 1. Introduction

China has abundant coal resources and mainly depends on coal as energy [1]. Coalbed methane (CBM) is a companion gas of coal, and it is also a kind of low-carbon energy. According to the statistics of China Energy Bureau, the total geological reserves of CBM are about 10 billion m3. In China, there are many problems on the extraction of CBM, such as slow methane pressure decreasing, rapid methane concentration reducing, and low permeability of coal seam. This makes it difficult to use the low concentration methane effectively which will be eventually discharged into the atmosphere. Efficient extraction of CBM and reduction of coal and gas outburst disasters have a profound impact on the simultaneous exploitation of coal and CBM [2].Due to the advantages of easy construction, eliminating gas outburst risks, the technologies of drilling holes along and across coal seam are widely used in CBM extraction. Many studies indicate that methane migration is influenced by gas-solid coupling effect between mining pressure and gas pressure of coal seam [3]. The deformation of coal induced by mining pressure would change the pore volume of matrix, which leads to the change of coal seam permeability. To quantitatively analyze the influence of matrix shrinkage and expansion on coal permeability, Liu et al. [4] proposed a coupled model of permeability based on the gas adsorption/desorption experiments and then analyzed the effect of matrix shrinkage on permeability. Since then, many scholars have done different researches on the coal seam permeability model with gas-solid coupling. Gilman and Beckie [5], Pillalamarry et al. [6], Chen [7], Sawyer et al. [8], and Seidle and Huitt [9] established relevant theoretical models of permeability based on different assumptions, respectively.In the mining process, Huy et al. [10] considered that with the increase of effective stress, the fracture opening of coal seam decreases gradually and some microfractures were completely closed, which made the permeability decrease sharply. Liang [11] introduced the geological storage conditions of Huainan mining area and studied the related technology of CBM extraction and drainage boreholes. Shi et al. [12] studied the relationship between loading rate and borehole destruction and proposed that the layout of methane drainage boreholes should take into account the mining pressure. Yin et al. [13] studied the methane extraction process under the mining pressure of the adjacent working face. The results show that the gas flow rate gradually decreases with the increase of mining pressure. With the help of KSE gas concentration recorder and SF6 gas conveyor, Wang et al. [14] studied the gas-solid coupling phenomenon between mining pressure and gas migration. The results showed that the change of mining pressure is the key factor of permeability change in the coal seam. Based on the research background of Zhangji coal mine, Xue et al. [15] analyzed the influence of different mining methods on the permeability of coal seam. Wei et al. [16] also pointed out that the coal permeability showed a u-shaped curve change in the coal sample full stress-strain experiment.In the previous studies, a series of theoretical models for methane extraction have been put forward on the safety of methane extraction. These models were applied to research the gas-solid coupling phenomenon in the process of methane extraction. However, most of the studies focused on the mechanism of gas migration and permeability change, and the influence of dynamic mining pressure change on the methane extraction in coal seam was seldom taken into account. In this study, we established a dual-porosity and dual-permeability methane drainage model with the periodic effect of mining pressure of the near working face to solve the specific problems in the methane drainage of 4102 working face in number 3 coal seam of Hebi Number 3 coal mine in China. The interaction relationship between gas and coal seam under the influence of mining was studied, and the methane concentration and flow rate were analyzed. On this basis, a solution to enhance the efficiency of methane extraction by injecting air was put forward. A simulation study was carried out to analyze the permeability change of coal seam and the gas migration rule after air injection.


## 2. Background of the Experimental Mine

### 2.1. Geological and Storage Conditions

Hebi Number 3 coal mine is located in Hebi City of Henan Province of China. Many coal and gas outburst accidents occurred since 1978. The geological conditions of number 3 coal seam are shown in Figure1(a), 4101 working face is the first mining working face, and 4102 working face is prepared for eliminating the danger of gas outburst before mining. The burial depth of two working faces and the average thickness of number 3 coal seam are 784 m and 5.2 m, respectively. And the methane content of number 3 coal seam is about 10.92–8.96 m3/t. The methane pressure measured in the field is between 1.58 and 2.36 MPa. Boreholes are drilled along the coal seam for methane extraction with a diameter of 94 mm, and negative extraction pressure is about 12–20 kPa.Figure 1
The location and extraction conditions of coal mining face.
(a)
Experimental coal mine location
(b)
Drilling and extraction pipeline of gas extraction in this coal seam


### 2.2. Main Problems in Gas Extraction

In fact, in the related study of mining pressure and methane emission, it is considered that under the influence of abutment pressure, due to the effect of “relieving pressure and increasing permeability of coal” [17], it will accelerate the migration of gas in the fracture. Finally, the amount of methane emission will increase, which will help to eliminate the danger of gas outburst. In view of the actual methane extraction of the 4102 working face, especially affected by the mining activities of 4101 working face, the mining pressure value and position of the 4102 working face are constantly changed. At this point, it is necessary to further study the coupling effect between the coal permeability change, free gas migration, and adsorption deformation of 4102 working face.As a gas outburst mine, Hebi Number 3 coal mine has taken a series of measures to eliminate coal and gas outburst, such as hydraulic fracturing and loose blasting of deep-hole explosives. At the same time, according to the regulations of coal mine safety, it is necessary to use the upward drilling of the roadway excavated in the floor of coal seam to extract methane technology [18]. But it is difficult to popularize because the effect of penetration and extraction cannot meet the requirement of mine safety production and the replacement of working face. Therefore, we have to take a more effective method to quickly reduce the methane pressure of coal seam.


## 2.1. Geological and Storage Conditions

Hebi Number 3 coal mine is located in Hebi City of Henan Province of China. Many coal and gas outburst accidents occurred since 1978. The geological conditions of number 3 coal seam are shown in Figure1(a), 4101 working face is the first mining working face, and 4102 working face is prepared for eliminating the danger of gas outburst before mining. The burial depth of two working faces and the average thickness of number 3 coal seam are 784 m and 5.2 m, respectively. And the methane content of number 3 coal seam is about 10.92–8.96 m3/t. The methane pressure measured in the field is between 1.58 and 2.36 MPa. Boreholes are drilled along the coal seam for methane extraction with a diameter of 94 mm, and negative extraction pressure is about 12–20 kPa.Figure 1
The location and extraction conditions of coal mining face.
(a)
Experimental coal mine location
(b)
Drilling and extraction pipeline of gas extraction in this coal seam


## 2.2. Main Problems in Gas Extraction

In fact, in the related study of mining pressure and methane emission, it is considered that under the influence of abutment pressure, due to the effect of “relieving pressure and increasing permeability of coal” [17], it will accelerate the migration of gas in the fracture. Finally, the amount of methane emission will increase, which will help to eliminate the danger of gas outburst. In view of the actual methane extraction of the 4102 working face, especially affected by the mining activities of 4101 working face, the mining pressure value and position of the 4102 working face are constantly changed. At this point, it is necessary to further study the coupling effect between the coal permeability change, free gas migration, and adsorption deformation of 4102 working face.As a gas outburst mine, Hebi Number 3 coal mine has taken a series of measures to eliminate coal and gas outburst, such as hydraulic fracturing and loose blasting of deep-hole explosives. At the same time, according to the regulations of coal mine safety, it is necessary to use the upward drilling of the roadway excavated in the floor of coal seam to extract methane technology [18]. But it is difficult to popularize because the effect of penetration and extraction cannot meet the requirement of mine safety production and the replacement of working face. Therefore, we have to take a more effective method to quickly reduce the methane pressure of coal seam.


## 3. Dual-Porosity Model for Methane Migration

### 3.1. Model Assumptions

The permeability of coal seam is the key influence factor for the methane extraction research. Compared with conventional natural gas reservoirs, the permeability of coal seam extracted by mines changes between 0.1 D and 0.001 mD (1 D = 10−12 m2) [19], which is completely different from the conventional natural gas extraction. In general, coal seam is a dual-porosity media containing matrix pores and fractures. Natural fractures of coal seam can be divided into horizontal bedding fracture and vertical cleat fracture. Methane is mainly stored in the pores by adsorbing on the micropore surface accounting for 70–95% [20, 21] of the total reserves. Free methane mainly migrates in coal seam fractures following the convection-diffusion equation. In the matrix system, there are many phenomena involved, but Fickian diffusion is the most dominant in the porous structure of the coal matrix [22, 23]. The main process of methane migration in the coal seam is shown in Figure 2.Figure 2
The process of methane migration in coal seam.In the following, a set of field equations for coal deformation, gas flow, and transport are defined. These field equations are coupled through new porosity and permeability models for matrix and fractures. The derivations are based on the following assumptions:(1)
Coal is a dual-poroelastic continuum.
(2)
The influence of coal seam moisture on gas migration is ignored.
(3)
Methane contained in the pores is ideal gas, and its viscosity is constant under isothermal conditions.
(4)
The air injected in coal seam is ideal saturated gas, and the adsorption effect of coal seam on air is ignored.


### 3.2. Coal Deformation Equation

As shown in Figure1, the ideal simplified model of coal seam is composed by fractures and matrices. The matrices are composed by many primary pores. The deformation of this dual-porosity model can be represented by a Navier-type equation [24]:
(1)Gui+G1−2νuk−αpmmix−βpfmix−KεLPLPL+Pm12+Fi=0,where G, ν, and K are the shear modulus, Poisson’s ratio, and bulk modulus of coal, respectively, α and β are the Biot coefficients of matrix and fracture, εs is gas sorption strain, and Pmmix and Pfmix are the pore pressure in fractures and matrix. The relationship of these parameters can be defined as follows:
(2)G=D21+ν,D=1E+1a0⋅Kf−1,α=1−KKs,β=1−Ka0⋅Kf,εs=εLpm1pm1+PL,where E is the elastic modulus of coal, Kf is the bulk modulus of fracture, Ks is the bulk modulus of matrix, Kf is the normal stiffness of fracture, Fi is the body force of element, usually is gravity, and εL and PL are the Langmuir adsorption strain and pressure constants.In (1), the first two terms represent the deformation of representative element volume (REV) under the influence of gas mixture pressure change and adsorption. The third and fourth terms represent the deformation of methane pressure in the matrix and fracture system, respectively. The fifth term represents the volume deformation due to the adsorption effect of matrix on gas.


### 3.3. Methane Migration Equation

According to Masoudian et al. [25], there are two types of coalbed methane reservoir models that may be used for methane flow in coal. The most primitive type is the diffusion which is assumed to be instantaneous, but this cannot fully describe the methane flow processes in coal seams. The second type, which is the most widely used for coalbed methane simulations, assumes coal is modeled with a dual-porosity system comprised of fractures and is microporous. Based on the methane migration speed, the model includes the Darcy convection equation for gas transfer in the fracture and Fick equation for gas diffusion in the matrix. The equation of gas mass transfers between them can be expressed as follows [26]:
(3)∂Mn∂t+∇vn⋅ρn−∇Dn⋅∇Mn=Qs,where M is the mass of gas in coal seam; the subscript “n” (n=1,2) represents methane and air; νn is the transfer velocity of gas in fracture; Dn is the diffusion coefficient of component gases; ρn is the gas density; t is the time; and Qs is the term of gas source, usually regarded as the gas initial pressure.Methane is mainly stored in matrix pores in adsorbed and free states. By ignoring the influence of matrix on the air, the mass of methane and air in matrix can be expressed as follows at any time [27]:
(4)Mm1=ρm1ϕm+ρgaρcVLpm1pm1+PL,Mm2=ρm2ϕm.Methane and air are mainly transferred in the fractures of coal, and the mass of mixed gas in the fractures can be expressed as follows at any time:(5)Mfn=ρgnϕf,where Mm1 is the mass (kg/m3) of gas. The subscripts of “m” and “f” indicate the matrix and fracture, respectively; ρga is the gas density at standard atmospheric pressure (kg/m3); ρc is the density of coal (kg/m3); PL is the Langmuir adsorption pressure constant (MPa); VL is the Langmuir adsorption volume constant (m3/t); ϕm is the matrix porosity (%); ρmn is the density of gas in the matrix, ρmn=mnPmn/RT; ρfn is the density of gas in the fracture, ρfn=mnPfn/RT.In regard to the term of gas migration in (4), the convection velocity of any gas can be expressed as follows (ignoring the effect of gravity):
(6)vfn=kfu∇pfmix,vmn=kmu∇pmmix,where kf and km are the fracture permeability and matrix permeability, respectively; u is the dynamic viscosity of gas mixture; and ∇pf and ∇pm are the pressure difference of gas flow in the fracture and matrix. Substituting (4), (5), and (6) into (3), the gas mixture flow control equation can be obtained.
(7)ϕm+ρcpgaVLPLpm1+PL2∂pm1∂t+pm1∂ϕm∂t=∇pmmixkmu∇pmmix+∇D1ϕm∇m1−ωpmmix−pfmix,(8)ϕf∂pf1∂t+pf1∂ϕf∂t−∇pfmixkfu∇pfmix−∇D1ϕf∇m1=ωpmmix−pfmix,(9)ϕm∂pm2∂t+pm2∂ϕm∂t−∇pmmixkmu∇pmmix−∇D2ϕm∇m2=−ωpmmix−pfmix,(10)ϕf∂pf2∂t+pf2∂ϕf∂t−∇pfmixkfu∇pfmix−∇D2ϕf∇m2=ωpmmix−pfmix,where ω=81+2/a02km/μ is the exchange coefficient between matrix and fracture [28] and a0 is the initial fracture width.


### 3.4. Dynamic Permeability Model of Coal Seam

For CBM extraction, the effective stress change and methane desorption cause deformation of the coal; Palmer and Mansoori, Shi and Durucan [29], and Cui and Bustin [30] had proposed a dynamic variation model of permeability based on different assumptions. On the basis of summarizing the development and application of coalbed methane, Pan and Connell [31] analyzed the existing reservoir permeability model and analyzed that the change of permeability is an important physical parameter for evaluating the extraction of coalbed methane. Mohsen et al. [32] also pointed out that the permeability model under vertical stress change should be considered and needs to find a balance between the complexity of the permeability model and the accuracy of the field permeability data. In this study, we use the dynamic permeability variation model proposed by Wu et al. [33] to study the CBM extraction under the mining pressure. However, unlike Wu’s model, the coupling effect of gas mixture pressure and different cycle mining pressures on gas migration was studied in this paper. The equation of coal seam matrix permeability can be expressed as
(11)ϕm=α+ϕm0−αexp1K1K+b0a0Kf−1εV−εV0−εS+εS0,where ϕm is the real-time porosity of the matrix and the subscript “0” indicates the initial state; εV=ε11+ε22+ε33 indicates the deformation of the representative element volume; and εS=εLpm1/pm1+PL is the gas adsorption strain in pores.In the same way as above, the porosity of coal seam fracture system can be expressed as(12)ϕfϕf0=1+1Kf1K+b0a0Kf−1εV−εV0−εS+εS0.During the methane extraction, with the continuous decrease of gas pressure, the pore volume of matrix changes continuously with time. The equation for porosity change within the matrix over time can be written as(13)∂ϕm∂t=ϕm−αexp1K1K+b0a0Kf−1∂εV∂t−∂εS∂t,(14)∂ϕf∂t=ϕf0Kf1K+b0a0Kf−1∂εV∂t−∂εS∂t.The permeability and porosity of coal seam matrix and fracture can be expressed in the cubic formula:(15)kfkf0=ϕfϕf03=1+1Kf1K+b0a0Kf−1εV−εV0−εS+εS03,kmkm0=ϕmϕm03=1−αKϕm01K+b0a0Kf−1εV−εV0−εS+εS03,where kf0 and km0 are the initial permeability of the fracture and matrix (in m2).


### 3.5. Coupled Field Equations

By substituting (13) and (14) into (7), (8), (9), and (10) of gas flow, respectively, the final equation of gas flow control equations can be expressed as
(16)ϕm+1−ϕm0ρcpgaPLVLpm1+PL2−αKfb0/a0K+KfPLεLpm1pm1+PL2∂pm∂t=∇kmupmmix∇pmmix+∇ϕmD1∇pmmix+81+2a02kmμpf1−pm1−αKfb0/a0K+Kfpm1∂εv∂t,(17)ϕf∂pf1∂t−∇kfupfmix∇pf−∇ϕfD1∇pfmix=−81+2a02kmμpf1−pm1−3ϕf0pf1ϕf0+3Kf/K∂εv∂t−εLPLpm1+PL2∂pm1∂t,(18)ϕm∂pm2∂t−∇kmupmmix∇pmmix−∇ϕmD2∇pmmix=81+2a02kmμpf2−pm2+pm2αK11/Kf+b0/a0KPLεLpm1+PL2∂pm1∂t−∂εV∂t,(19)ϕf∂pf2∂t−∇kfupfmix∇pfmix−∇ϕfD2∇pfmix=−81+2a02kmμpf2−pm2−pf2ϕf0Kf1K+b0a0Kf−1∂εV∂t−PLεLpm1+PL2∂pm1∂t.Equations (16), (17), (18), and (19) describe the flow of mixed gas within fracture of coal seam under mining pressure. In (16), the first term ∂pm/∂t at the left is composed of three subterms sequentially including (I) pressure variation of free gas in matrix pores; (II) deformation of the matrix caused by gas desorption; and (III) deformation of fracture volume. The third and fourth terms on the right of the equation indicate the exchange of free gas between the fracture and matrix, methane released by the deformation of coal. The first term on the left of (17) represents the pressure variation of methane in the fracture. The second term on the right of the equation contains two subterms, that is, the effect of coal seam volume deformation on gas and the effect of methane desorption strain. The physical meanings of (18) and (19) are the same as those of (16) and (17). The above coupling control equations contain the terms of time and space derivation, which need to be solved by Comsol Multiphysics. The interaction between the physical fields is shown in Figure 3.Figure 3
Coupling relationship between coal deformation and methane migration.


## 3.1. Model Assumptions

The permeability of coal seam is the key influence factor for the methane extraction research. Compared with conventional natural gas reservoirs, the permeability of coal seam extracted by mines changes between 0.1 D and 0.001 mD (1 D = 10−12 m2) [19], which is completely different from the conventional natural gas extraction. In general, coal seam is a dual-porosity media containing matrix pores and fractures. Natural fractures of coal seam can be divided into horizontal bedding fracture and vertical cleat fracture. Methane is mainly stored in the pores by adsorbing on the micropore surface accounting for 70–95% [20, 21] of the total reserves. Free methane mainly migrates in coal seam fractures following the convection-diffusion equation. In the matrix system, there are many phenomena involved, but Fickian diffusion is the most dominant in the porous structure of the coal matrix [22, 23]. The main process of methane migration in the coal seam is shown in Figure 2.Figure 2
The process of methane migration in coal seam.In the following, a set of field equations for coal deformation, gas flow, and transport are defined. These field equations are coupled through new porosity and permeability models for matrix and fractures. The derivations are based on the following assumptions:(1)
Coal is a dual-poroelastic continuum.
(2)
The influence of coal seam moisture on gas migration is ignored.
(3)
Methane contained in the pores is ideal gas, and its viscosity is constant under isothermal conditions.
(4)
The air injected in coal seam is ideal saturated gas, and the adsorption effect of coal seam on air is ignored.


## 3.2. Coal Deformation Equation

As shown in Figure1, the ideal simplified model of coal seam is composed by fractures and matrices. The matrices are composed by many primary pores. The deformation of this dual-porosity model can be represented by a Navier-type equation [24]:
(1)Gui+G1−2νuk−αpmmix−βpfmix−KεLPLPL+Pm12+Fi=0,where G, ν, and K are the shear modulus, Poisson’s ratio, and bulk modulus of coal, respectively, α and β are the Biot coefficients of matrix and fracture, εs is gas sorption strain, and Pmmix and Pfmix are the pore pressure in fractures and matrix. The relationship of these parameters can be defined as follows:
(2)G=D21+ν,D=1E+1a0⋅Kf−1,α=1−KKs,β=1−Ka0⋅Kf,εs=εLpm1pm1+PL,where E is the elastic modulus of coal, Kf is the bulk modulus of fracture, Ks is the bulk modulus of matrix, Kf is the normal stiffness of fracture, Fi is the body force of element, usually is gravity, and εL and PL are the Langmuir adsorption strain and pressure constants.In (1), the first two terms represent the deformation of representative element volume (REV) under the influence of gas mixture pressure change and adsorption. The third and fourth terms represent the deformation of methane pressure in the matrix and fracture system, respectively. The fifth term represents the volume deformation due to the adsorption effect of matrix on gas.


## 3.3. Methane Migration Equation

According to Masoudian et al. [25], there are two types of coalbed methane reservoir models that may be used for methane flow in coal. The most primitive type is the diffusion which is assumed to be instantaneous, but this cannot fully describe the methane flow processes in coal seams. The second type, which is the most widely used for coalbed methane simulations, assumes coal is modeled with a dual-porosity system comprised of fractures and is microporous. Based on the methane migration speed, the model includes the Darcy convection equation for gas transfer in the fracture and Fick equation for gas diffusion in the matrix. The equation of gas mass transfers between them can be expressed as follows [26]:
(3)∂Mn∂t+∇vn⋅ρn−∇Dn⋅∇Mn=Qs,where M is the mass of gas in coal seam; the subscript “n” (n=1,2) represents methane and air; νn is the transfer velocity of gas in fracture; Dn is the diffusion coefficient of component gases; ρn is the gas density; t is the time; and Qs is the term of gas source, usually regarded as the gas initial pressure.Methane is mainly stored in matrix pores in adsorbed and free states. By ignoring the influence of matrix on the air, the mass of methane and air in matrix can be expressed as follows at any time [27]:
(4)Mm1=ρm1ϕm+ρgaρcVLpm1pm1+PL,Mm2=ρm2ϕm.Methane and air are mainly transferred in the fractures of coal, and the mass of mixed gas in the fractures can be expressed as follows at any time:(5)Mfn=ρgnϕf,where Mm1 is the mass (kg/m3) of gas. The subscripts of “m” and “f” indicate the matrix and fracture, respectively; ρga is the gas density at standard atmospheric pressure (kg/m3); ρc is the density of coal (kg/m3); PL is the Langmuir adsorption pressure constant (MPa); VL is the Langmuir adsorption volume constant (m3/t); ϕm is the matrix porosity (%); ρmn is the density of gas in the matrix, ρmn=mnPmn/RT; ρfn is the density of gas in the fracture, ρfn=mnPfn/RT.In regard to the term of gas migration in (4), the convection velocity of any gas can be expressed as follows (ignoring the effect of gravity):
(6)vfn=kfu∇pfmix,vmn=kmu∇pmmix,where kf and km are the fracture permeability and matrix permeability, respectively; u is the dynamic viscosity of gas mixture; and ∇pf and ∇pm are the pressure difference of gas flow in the fracture and matrix. Substituting (4), (5), and (6) into (3), the gas mixture flow control equation can be obtained.
(7)ϕm+ρcpgaVLPLpm1+PL2∂pm1∂t+pm1∂ϕm∂t=∇pmmixkmu∇pmmix+∇D1ϕm∇m1−ωpmmix−pfmix,(8)ϕf∂pf1∂t+pf1∂ϕf∂t−∇pfmixkfu∇pfmix−∇D1ϕf∇m1=ωpmmix−pfmix,(9)ϕm∂pm2∂t+pm2∂ϕm∂t−∇pmmixkmu∇pmmix−∇D2ϕm∇m2=−ωpmmix−pfmix,(10)ϕf∂pf2∂t+pf2∂ϕf∂t−∇pfmixkfu∇pfmix−∇D2ϕf∇m2=ωpmmix−pfmix,where ω=81+2/a02km/μ is the exchange coefficient between matrix and fracture [28] and a0 is the initial fracture width.


## 3.4. Dynamic Permeability Model of Coal Seam

For CBM extraction, the effective stress change and methane desorption cause deformation of the coal; Palmer and Mansoori, Shi and Durucan [29], and Cui and Bustin [30] had proposed a dynamic variation model of permeability based on different assumptions. On the basis of summarizing the development and application of coalbed methane, Pan and Connell [31] analyzed the existing reservoir permeability model and analyzed that the change of permeability is an important physical parameter for evaluating the extraction of coalbed methane. Mohsen et al. [32] also pointed out that the permeability model under vertical stress change should be considered and needs to find a balance between the complexity of the permeability model and the accuracy of the field permeability data. In this study, we use the dynamic permeability variation model proposed by Wu et al. [33] to study the CBM extraction under the mining pressure. However, unlike Wu’s model, the coupling effect of gas mixture pressure and different cycle mining pressures on gas migration was studied in this paper. The equation of coal seam matrix permeability can be expressed as
(11)ϕm=α+ϕm0−αexp1K1K+b0a0Kf−1εV−εV0−εS+εS0,where ϕm is the real-time porosity of the matrix and the subscript “0” indicates the initial state; εV=ε11+ε22+ε33 indicates the deformation of the representative element volume; and εS=εLpm1/pm1+PL is the gas adsorption strain in pores.In the same way as above, the porosity of coal seam fracture system can be expressed as(12)ϕfϕf0=1+1Kf1K+b0a0Kf−1εV−εV0−εS+εS0.During the methane extraction, with the continuous decrease of gas pressure, the pore volume of matrix changes continuously with time. The equation for porosity change within the matrix over time can be written as(13)∂ϕm∂t=ϕm−αexp1K1K+b0a0Kf−1∂εV∂t−∂εS∂t,(14)∂ϕf∂t=ϕf0Kf1K+b0a0Kf−1∂εV∂t−∂εS∂t.The permeability and porosity of coal seam matrix and fracture can be expressed in the cubic formula:(15)kfkf0=ϕfϕf03=1+1Kf1K+b0a0Kf−1εV−εV0−εS+εS03,kmkm0=ϕmϕm03=1−αKϕm01K+b0a0Kf−1εV−εV0−εS+εS03,where kf0 and km0 are the initial permeability of the fracture and matrix (in m2).


## 3.5. Coupled Field Equations

By substituting (13) and (14) into (7), (8), (9), and (10) of gas flow, respectively, the final equation of gas flow control equations can be expressed as
(16)ϕm+1−ϕm0ρcpgaPLVLpm1+PL2−αKfb0/a0K+KfPLεLpm1pm1+PL2∂pm∂t=∇kmupmmix∇pmmix+∇ϕmD1∇pmmix+81+2a02kmμpf1−pm1−αKfb0/a0K+Kfpm1∂εv∂t,(17)ϕf∂pf1∂t−∇kfupfmix∇pf−∇ϕfD1∇pfmix=−81+2a02kmμpf1−pm1−3ϕf0pf1ϕf0+3Kf/K∂εv∂t−εLPLpm1+PL2∂pm1∂t,(18)ϕm∂pm2∂t−∇kmupmmix∇pmmix−∇ϕmD2∇pmmix=81+2a02kmμpf2−pm2+pm2αK11/Kf+b0/a0KPLεLpm1+PL2∂pm1∂t−∂εV∂t,(19)ϕf∂pf2∂t−∇kfupfmix∇pfmix−∇ϕfD2∇pfmix=−81+2a02kmμpf2−pm2−pf2ϕf0Kf1K+b0a0Kf−1∂εV∂t−PLεLpm1+PL2∂pm1∂t.Equations (16), (17), (18), and (19) describe the flow of mixed gas within fracture of coal seam under mining pressure. In (16), the first term ∂pm/∂t at the left is composed of three subterms sequentially including (I) pressure variation of free gas in matrix pores; (II) deformation of the matrix caused by gas desorption; and (III) deformation of fracture volume. The third and fourth terms on the right of the equation indicate the exchange of free gas between the fracture and matrix, methane released by the deformation of coal. The first term on the left of (17) represents the pressure variation of methane in the fracture. The second term on the right of the equation contains two subterms, that is, the effect of coal seam volume deformation on gas and the effect of methane desorption strain. The physical meanings of (18) and (19) are the same as those of (16) and (17). The above coupling control equations contain the terms of time and space derivation, which need to be solved by Comsol Multiphysics. The interaction between the physical fields is shown in Figure 3.Figure 3
Coupling relationship between coal deformation and methane migration.


## 4. Gas Migration Characteristics under the Mining Pressure

### 4.1. Numerical Model and Input Parameters

The model is established based on a section of 4102 working face of number 3 coal seam with dimensions of 40 m (length) × 5 m (height). Negative-pressure drainage boreholes were drilled in advance and arranged in the horizontal center of number 3 coal seam. The diameter of boreholes is 100 mm, and the horizontal spacing is 3 m. The numerical model of the 4102 working face is shown in Figure4 and the input parameters in Tables 1 and 2 as required. In this part, we only simulated the methane extraction under mining pressure.Figure 4
Numerical model of methane extraction in 4102 working face.Table 1
Numerical model physical parameters.
Parameter
Value
Parameter
Value
Coal elastic modulusE (MPa)
2813
Adsorption pressure constantPL (MPa)
2.07
Matrix elastic modulusEs (MPa)
8439
Adsorption volume constantVL (m3/kg)
0.0256
Poisson’s ratioν
0.34
Adsorption pressure strainεL
0.015
Coal densityρc (kg/m3)
1250
Initial matrix porosity.φm0
0.05
Air density (kg/m3)
0.717
Initial matrix permeabilitykm0
2.8·10−18
Facture stiffnessKf (MPa)
4800
Methane diffusion rateDn (m2)
3.6·10−12
Dynamic viscosityu (Pa·s)
1.84·10−5
Initial fracture porosityφf0
0.002
Bulk modulusK (MPa)
5400
Initial fracture permeabilitykf0
2.08·10−17
Matrix widtha0 (m)
5·10−6
Initial temperatureT (K)
300
Frack widthb0 (m)
0.005
Air diffusion rateDn (m2)
5.8·10−12Table 2
Numerical model boundary conditions.
Coal seam deformation
Value
Methane migration
Value
Upper boundary (MPa)
14 transition to 21
Pressure of methane extraction borehole (kPa)
20
Lower boundary
Fixed boundary
Left boundary
Roller support
Right boundary
Methane pressure (MPa)
2
Borehole boundary
Free deformation
Initial displacement
0
Others (MPa)
0As the advancing of 4101 working face, the mining pressure applied in 4102 working face would change as shown in Figure5. In order to present the mining pressure, the pressure applied on the top of model can be divided into three iteration cycles. When the 4101 working face advanced by 10 meters totally, the mining pressure boundary was set again and the advance repeated in sequence until the end of the methane extraction. The initial pressure of methane in the coal (fracture and matrix) was set to 2 MPa, and the air pressure was set to 0 MPa. The migration of methane into the boreholes was only considered. Three points of A, B, and C of the model were selected as the data source for postprocessing with coordinates of (0.2, 2.5), (22.5, 2.5), and (39.8, 2.5).Figure 5
4102 preparation work face.


### 4.2. Simulation Results

Under the influence of the mining in the 4101 working face, the initial stress equilibrium state of the coal seam is broken, and the different stress distribution regions are formed along the direction of the coal seam, which further affected the permeability of the coal seam. In the vicinity of gas drainage holes, there were three zones, that is, plastic failure zone, stress concentration zone, and original rock stress zone [34]. The permeability variations at the plastic failure zone and data points around the boreholes on both sides of the model are shown in Figure 6.Figure 6
Different calculation steps, permeability, and plastic zone distribution.
(a)
The first calculation step, permeability, and plastic zone change
(b)
The last calculation step, permeability, and plastic zone changeAs shown in Figure6(a), the permeability in the abutment pressure zone (right side of the model) was significantly lower than that in the original rock stress zone. During the whole extraction time, the variation of permeability was affected by the effective stress (Fy−Pfmix) and the strain of methane desorption. The permeability of coal seam (Kf) depends on the competition between them. In the early stage of extraction, the effective stress near the abutment pressure zone was obviously higher than the original rock stress zone, so the permeability was lower than that of other zones inhibiting the gas migration. With the increase in extraction time, the gas in adsorbed state was desorbed gradually, resulting in matrix shrinkage and gradual increase in the permeability. At the end of extraction, the permeability of coal seam in the original rock stress zone changed by about 5% (0.8–1.2). In the vicinity of the drainage boreholes, the drilling and excavation improved the expansion and extension of the original fractures leading to a sharp increase in the gas permeability of the coal seam. The permeability of the coal seam was significantly higher than the external permeability, and the fracture permeability increased by 5 times (0.8–4). The fractures in plastic failure zone became the dominant factor of permeability.Figure7 shows the methane pressure distribution at different extraction times. The gas migration during the extraction represented a typical gas-solid coupling. In the effect of extraction negative pressure, the gas pressure showed an exponentially decreasing trend. The trend line of methane pressure at point C went down significantly faster than that in the matrix. The decrease rate of gas pressure at point C was significantly slower than that at point B as the effective pressure of the coal seam was higher and the opening of the original fractures was smaller in the abutment pressure zone. The methane pressure difference (Pm−Pf) significantly increased at the early stage of the extraction (10 days ago) as there was pressure difference of gas transferred between the matrix and fracture. At the late stage of the extraction, the gas pressure difference gradually decreased until it was zero.Figure 7
Coal seam methane pressure changes with time.


### 4.3. Comparisons with Field Data

The field data of methane pressure at number 1 and number 3 drainage boreholes were collected as shown in Figure1. The methane pressure data obtained from the model were compared with the field data as shown in Figure 8. In the figure, although the methane pressure showed an overall decreasing trend, its downward trend was affected by the geological and storage conditions of the coal seam and the negative pressure of the drainage drilled holes resulting in fluctuations in pressure. The methane extraction numerical model was established with ideal conditions, ignoring the ash content and moisture in coal seam or other geological factors. Although their data were different, their curves showed the same change trend. This model can be used to predict and evaluate the long-term effect of borehole gas extraction.Figure 8
Comparison of field methane pressure and numerical model results.
(a)
Methane pressure varies with extraction time
(b)
Methane concentration changes with the time


## 4.1. Numerical Model and Input Parameters

The model is established based on a section of 4102 working face of number 3 coal seam with dimensions of 40 m (length) × 5 m (height). Negative-pressure drainage boreholes were drilled in advance and arranged in the horizontal center of number 3 coal seam. The diameter of boreholes is 100 mm, and the horizontal spacing is 3 m. The numerical model of the 4102 working face is shown in Figure4 and the input parameters in Tables 1 and 2 as required. In this part, we only simulated the methane extraction under mining pressure.Figure 4
Numerical model of methane extraction in 4102 working face.Table 1
Numerical model physical parameters.
Parameter
Value
Parameter
Value
Coal elastic modulusE (MPa)
2813
Adsorption pressure constantPL (MPa)
2.07
Matrix elastic modulusEs (MPa)
8439
Adsorption volume constantVL (m3/kg)
0.0256
Poisson’s ratioν
0.34
Adsorption pressure strainεL
0.015
Coal densityρc (kg/m3)
1250
Initial matrix porosity.φm0
0.05
Air density (kg/m3)
0.717
Initial matrix permeabilitykm0
2.8·10−18
Facture stiffnessKf (MPa)
4800
Methane diffusion rateDn (m2)
3.6·10−12
Dynamic viscosityu (Pa·s)
1.84·10−5
Initial fracture porosityφf0
0.002
Bulk modulusK (MPa)
5400
Initial fracture permeabilitykf0
2.08·10−17
Matrix widtha0 (m)
5·10−6
Initial temperatureT (K)
300
Frack widthb0 (m)
0.005
Air diffusion rateDn (m2)
5.8·10−12Table 2
Numerical model boundary conditions.
Coal seam deformation
Value
Methane migration
Value
Upper boundary (MPa)
14 transition to 21
Pressure of methane extraction borehole (kPa)
20
Lower boundary
Fixed boundary
Left boundary
Roller support
Right boundary
Methane pressure (MPa)
2
Borehole boundary
Free deformation
Initial displacement
0
Others (MPa)
0As the advancing of 4101 working face, the mining pressure applied in 4102 working face would change as shown in Figure5. In order to present the mining pressure, the pressure applied on the top of model can be divided into three iteration cycles. When the 4101 working face advanced by 10 meters totally, the mining pressure boundary was set again and the advance repeated in sequence until the end of the methane extraction. The initial pressure of methane in the coal (fracture and matrix) was set to 2 MPa, and the air pressure was set to 0 MPa. The migration of methane into the boreholes was only considered. Three points of A, B, and C of the model were selected as the data source for postprocessing with coordinates of (0.2, 2.5), (22.5, 2.5), and (39.8, 2.5).Figure 5
4102 preparation work face.


## 4.2. Simulation Results

Under the influence of the mining in the 4101 working face, the initial stress equilibrium state of the coal seam is broken, and the different stress distribution regions are formed along the direction of the coal seam, which further affected the permeability of the coal seam. In the vicinity of gas drainage holes, there were three zones, that is, plastic failure zone, stress concentration zone, and original rock stress zone [34]. The permeability variations at the plastic failure zone and data points around the boreholes on both sides of the model are shown in Figure 6.Figure 6
Different calculation steps, permeability, and plastic zone distribution.
(a)
The first calculation step, permeability, and plastic zone change
(b)
The last calculation step, permeability, and plastic zone changeAs shown in Figure6(a), the permeability in the abutment pressure zone (right side of the model) was significantly lower than that in the original rock stress zone. During the whole extraction time, the variation of permeability was affected by the effective stress (Fy−Pfmix) and the strain of methane desorption. The permeability of coal seam (Kf) depends on the competition between them. In the early stage of extraction, the effective stress near the abutment pressure zone was obviously higher than the original rock stress zone, so the permeability was lower than that of other zones inhibiting the gas migration. With the increase in extraction time, the gas in adsorbed state was desorbed gradually, resulting in matrix shrinkage and gradual increase in the permeability. At the end of extraction, the permeability of coal seam in the original rock stress zone changed by about 5% (0.8–1.2). In the vicinity of the drainage boreholes, the drilling and excavation improved the expansion and extension of the original fractures leading to a sharp increase in the gas permeability of the coal seam. The permeability of the coal seam was significantly higher than the external permeability, and the fracture permeability increased by 5 times (0.8–4). The fractures in plastic failure zone became the dominant factor of permeability.Figure7 shows the methane pressure distribution at different extraction times. The gas migration during the extraction represented a typical gas-solid coupling. In the effect of extraction negative pressure, the gas pressure showed an exponentially decreasing trend. The trend line of methane pressure at point C went down significantly faster than that in the matrix. The decrease rate of gas pressure at point C was significantly slower than that at point B as the effective pressure of the coal seam was higher and the opening of the original fractures was smaller in the abutment pressure zone. The methane pressure difference (Pm−Pf) significantly increased at the early stage of the extraction (10 days ago) as there was pressure difference of gas transferred between the matrix and fracture. At the late stage of the extraction, the gas pressure difference gradually decreased until it was zero.Figure 7
Coal seam methane pressure changes with time.


## 4.3. Comparisons with Field Data

The field data of methane pressure at number 1 and number 3 drainage boreholes were collected as shown in Figure1. The methane pressure data obtained from the model were compared with the field data as shown in Figure 8. In the figure, although the methane pressure showed an overall decreasing trend, its downward trend was affected by the geological and storage conditions of the coal seam and the negative pressure of the drainage drilled holes resulting in fluctuations in pressure. The methane extraction numerical model was established with ideal conditions, ignoring the ash content and moisture in coal seam or other geological factors. Although their data were different, their curves showed the same change trend. This model can be used to predict and evaluate the long-term effect of borehole gas extraction.Figure 8
Comparison of field methane pressure and numerical model results.
(a)
Methane pressure varies with extraction time
(b)
Methane concentration changes with the time


## 5. Enhanced Gas Extraction Test by Injection Air

### 5.1. Air Injection Parameters

To eliminate gas outburst risks and reduce methane pressure as soon as possible, by injecting air into the existing boreholes of the 4102 working face, the gas migration was improved in the coal seam to increase the methane extraction effect. Figure9 shows the injection boundary conditions and drilled boreholes in the model. The deformation field of the coal seam was established with the same settings as that of Section 4.1. The calculation parameters required by the model were the same as those in Table 1, and the air transport equation is shown in Table 3.Figure 9
Setting of numerical model for methane extraction under air injection.Table 3
Air migration equation setting.
Air migration
Value
Injection borehole position
Value
Initial pressure (MPa)
0
Extraction borehole number
1, 3, 5,7, 9, 11
Injection pressure (MPa)
0.3, 0.6, 1
Borehole boundary
Free deformation
Air injection borehole number
2, 4, 6,8, 10
Initial displacement
0


### 5.2. Gas Pressure Distribution

Figure10 shows the effect of injection air pressure on methane pressure drop. The dotted line in the figure shows the change of methane pressure without air injection, and the solid line shows the change of methane pressure under different injected air pressures. When there was no air injected into the coal seam, the methane pressure at point C decreased to the safety pressure of 0.75 MPa, which took about 12 days. When the air pressure was increased to 1 MPa, the methane pressure decreased to 0.75 MPa, which only took about 8 days. The injection of air could significantly reduce the methane pressure.Figure 10
Methane pressure changes with time during air injection.


### 5.3. Evaluation of Coal Seam Permeability

Figure11 shows the change trend of fracture permeability of coal seam at points B and C. During the whole extraction period, the permeability at the two points showed an increasing trend. When the air pressure in the fracture continued to increase, the effective stress of the coal seam continuously decreased, so the permeability at point B in the initial stage of injection was significantly higher than that at point C without injection. After 30 days of air injection, the air pressure in the model reached an equilibrium; the coal seam permeability at points B and C increased by 10% and 8%, respectively, compared to that in the initial stage of the extraction.Figure 11
Coal seam permeability change under air injection.


### 5.4. Gas Flow Characters in Fracture

Figure12 shows the variation of methane concentration. In the figure, the solid black line represents the methane concentration which decreased slowly when there was no air injection and showed an approximately linear decreasing trend in the extraction. After 15 days of extraction, the concentration only decreased by about 23%. The methane concentration decreased by 50% and 80%, respectively, compared with that under the air injection pressure of 0.6 MPa and 1 MPa in the condition that the extraction time was the same. By comparing the drop curve of methane pressure at point B in Figure 10, the methane pressure decreased to less than half of the initial value.Figure 12
Methane concentration varies with time.Figure13 shows the variation of air concentration at point B and the change curves of methane flow rate under different air pressures. Figure 13(b) shows that the methane flow rate was significantly affected by the air pressure, that is, the higher the air pressure, the faster the methane flow rate at point B decreased. This was because the increasing of air pressure provided an additional power for methane to migrate in the fracture. Meanwhile, the fracture permeability was affected by the deformation of adsorbed methane desorption and the volume deformation of the coal. The desorption of methane caused continuous shrinkage of matrix eventually improving the permeability by 10% (Figure 11). The additional air pressure and permeability improvement accelerated the flow of methane in the fracture.Figure 13
Methane flow rate changes with time.
(a)
Air concentration varies with time
(b)
Methane flow rate changes with the time


## 5.1. Air Injection Parameters

To eliminate gas outburst risks and reduce methane pressure as soon as possible, by injecting air into the existing boreholes of the 4102 working face, the gas migration was improved in the coal seam to increase the methane extraction effect. Figure9 shows the injection boundary conditions and drilled boreholes in the model. The deformation field of the coal seam was established with the same settings as that of Section 4.1. The calculation parameters required by the model were the same as those in Table 1, and the air transport equation is shown in Table 3.Figure 9
Setting of numerical model for methane extraction under air injection.Table 3
Air migration equation setting.
Air migration
Value
Injection borehole position
Value
Initial pressure (MPa)
0
Extraction borehole number
1, 3, 5,7, 9, 11
Injection pressure (MPa)
0.3, 0.6, 1
Borehole boundary
Free deformation
Air injection borehole number
2, 4, 6,8, 10
Initial displacement
0


## 5.2. Gas Pressure Distribution

Figure10 shows the effect of injection air pressure on methane pressure drop. The dotted line in the figure shows the change of methane pressure without air injection, and the solid line shows the change of methane pressure under different injected air pressures. When there was no air injected into the coal seam, the methane pressure at point C decreased to the safety pressure of 0.75 MPa, which took about 12 days. When the air pressure was increased to 1 MPa, the methane pressure decreased to 0.75 MPa, which only took about 8 days. The injection of air could significantly reduce the methane pressure.Figure 10
Methane pressure changes with time during air injection.


## 5.3. Evaluation of Coal Seam Permeability

Figure11 shows the change trend of fracture permeability of coal seam at points B and C. During the whole extraction period, the permeability at the two points showed an increasing trend. When the air pressure in the fracture continued to increase, the effective stress of the coal seam continuously decreased, so the permeability at point B in the initial stage of injection was significantly higher than that at point C without injection. After 30 days of air injection, the air pressure in the model reached an equilibrium; the coal seam permeability at points B and C increased by 10% and 8%, respectively, compared to that in the initial stage of the extraction.Figure 11
Coal seam permeability change under air injection.


## 5.4. Gas Flow Characters in Fracture

Figure12 shows the variation of methane concentration. In the figure, the solid black line represents the methane concentration which decreased slowly when there was no air injection and showed an approximately linear decreasing trend in the extraction. After 15 days of extraction, the concentration only decreased by about 23%. The methane concentration decreased by 50% and 80%, respectively, compared with that under the air injection pressure of 0.6 MPa and 1 MPa in the condition that the extraction time was the same. By comparing the drop curve of methane pressure at point B in Figure 10, the methane pressure decreased to less than half of the initial value.Figure 12
Methane concentration varies with time.Figure13 shows the variation of air concentration at point B and the change curves of methane flow rate under different air pressures. Figure 13(b) shows that the methane flow rate was significantly affected by the air pressure, that is, the higher the air pressure, the faster the methane flow rate at point B decreased. This was because the increasing of air pressure provided an additional power for methane to migrate in the fracture. Meanwhile, the fracture permeability was affected by the deformation of adsorbed methane desorption and the volume deformation of the coal. The desorption of methane caused continuous shrinkage of matrix eventually improving the permeability by 10% (Figure 11). The additional air pressure and permeability improvement accelerated the flow of methane in the fracture.Figure 13
Methane flow rate changes with time.
(a)
Air concentration varies with time
(b)
Methane flow rate changes with the time


## 6. Conclusions

To study the impact of mining pressure on the methane extraction in coal, a dual-porosity and dual-permeability model was established with gas seepage, effective stress, and adsorption deformation. The gas-solid coupling between methane migration and borehole damage under the influence of abutment mining pressure was studied. On this basis, a stimulation treatment for methane extraction by injecting air was presented. Based on this study, the following conclusions can be obtained:(1)
When the boreholes are affected by the abutment mining pressure, the plastic failure zone around the boreholes increases the permeability of the coal and promotes the methane pressure drops rapidly, eliminating the danger of gas outbursts. On the contrary, radial extension along the boreholes, especially in the abutment pressure zone, with the increase of effective stress, not only the permeability decreases continuously but also limits the methane migration within the fracture. In the original rock stress zone, the change trend of coal permeability is between plastic failure zone and abutment pressure zone.
(2)
The simulation results show that the methane in the fracture first flows into the drainage boreholes, then the methane in the matrix gradually diffuses. This is the reason why the methane pressure within coal is not on average always decreasing. Compared with the field data, the correctness of the methane drainage model under the abutment pressure near the working face was verified.
(3)
When air is injected into the coal, it first flows into the fractures and then reduces the effective stress of the coal, caused the fracture opening to increase continuously, and finally improves the permeability of coal. Extra air provides a new source of power for gas flow and promotes the flow of gas into the boreholes at the same time. This is the reason why the methane pressure within the fracture is always decreasing.



---

*Source: 4034296-2018-08-05.xml*