# Effects of Injection Pressure and Duration on Alternate High-Pressure Water-Gas Sequestration of Coalbed Methane

Authors: Ke Li; Chuanjie Zhu; Siyuan Liu; Deliang Chen; Guanghui Cai
Journal: Geofluids
 (2022)

---

## Abstract

To solve the problems of high ground stress, low permeability and low coalbed methane recovery, this paper proposes an optimization technique for improving coal seam permeability, i.e., water-gas alternating displacement enhances coalbed methane recovery (WGA-ECBM). The seepage field of gas-water two-phase flow is derived and a multiphase fluid-solid coupling model between physical parameters is established. This paper studies the effects of three displacement methods (water injection, air injection and water-air alternating injection) on methane recovery. The influence of injection pressure and injection duration on the WGA-ECBM technique is also studied. TheVCH4 (methane volume fraction) of the three methods decrease by 17.84%, 12.72% and 44.62% in the borehole, respectively. The influence ranges (diameter) are 24.68 m, 149.73 m and 135.21 m, respectively. Pore pressures and displacement area also increase as injection pressure increases in a 5 MPa gradient. And VCH4 in the borehole decrease by 35.57%, 43.82%, 46.78% and 51.72% from 10 MPa to 25 MPa, respectively. With the increase of injection duration, the influence range expands significantly but the pore pressure distribution changes little. And VCH4 decreases slightly in the displaced area and the displacement velocity gradually slows down. This technique promotes multiphase flow migration in coal seam, and increases the displacement range and methane recovery efficiency.


---

## Body

## 1. Introduction

The CBM production in China is significantly lower than that in other coal-rich countries, such as Australia and the United States[1]. Due to the different geological conditions, the permeability of some coal seams is low, which determines the difficulty of methane extraction and drainage[2]. Many researchers put forward a series of enhanced coalbed methane recovery (ECBM) techniques to solve the above problems, such as large-diameter borehole drainage [3], hydraulic punching [4], pulsating shock [5], hydraulic cutting [6] (or ultra-high pressure hydraulic slitting), hydraulic fracturing [7], CO2-ECBM[8] and pre-split blasting[9]. The basic mechanism of these techniques is to generate artificial fractures in coal seams, which can improve the permeability of coal seams and thus increase the methane recovery rate. However, methane is mainly adsorbed in nanopores of coal with high adsorption ability, which makes it difficult to effectively extract methane adsorbed[10]. Therefore, the problem of low methane recovery rate in coal seam cannot be fundamentally solved. Water displacement is considered an effective method to recover methane from small coal pores [11].Water flooding has encountered some problems in the field practice, such as poor methane recovery rate [12]. The water diffusion in coal seam is affected by viscous force, friction force, and normal stress of fracture wall, and water lock occurs in coal seam due to capillary effect and Jiamin phenomenon [13]. Recently, some researchers use acid-base salt solutions to improve the injected water [14], which still does not completely solve the water lock damage, and only to reduce the effect of bound water [15]. CO2 enhanced coalbed methane recovery (CO2-ECBM) is a good alternative method to obtain a high coalbed methane recovery rate [16], involving CO2 injected into deep coal seams [17, 18]. However, CO2 injection can cause coal matrix expansion, which has decreased gas permeability of coal [19]. The adsorption capacity of CO2 in coal is higher than that of methane, which may lead to high risk of coal and gas outburst or accidental emission of CO2 during coal mining. Considering the potential hazard of CO2-ECBM, nitrogen injection into coal seam was proposed [14]. N2 is mainly produced by separation from the atmosphere or decomposition of nitrogenous compounds, which requires a cost to produce, store and transport. The production of large amounts of N2 and CO2 is not economical. The adsorption of CH4 and N2 on coal matrix belongs to physical adsorption, and the adsorption capacity of CH4 on coal matrix is stronger than that of N2. Air containing 78% nitrogen is commonly used in field practice.Thus, we proposed a new technique named high-pressure water and gas alternating sequestration (H-P-WGAS) in our previous work [20]. In that work, we discussed the effectiveness of the H-P-WGAS technique and verified that the H-P-WGAS technique can enhance methane migration in coal seam. The field tests indicated that the average methane drainage flow rate of single borehole increased significantly and consequently the residual methane content in coal seam decreased dramatically as a result. In the present work, the effect of injection pressure and duration on water and gas alternating sequestration of coalbed methane will be discussed, which is important to the engineering practice.


## 2. Numerical Model

In our study, water and air are successively injected into coal seam to displace methane. Based on the theory of multiphase-flow in porous media, the model considered fluid seepage and capillary force in coal seam, and the basic assumptions were as follow:(1)
Coal is a double porous medium rich in pores and fractures. The pore system and fracture system are assumed to be continuum systems(2)
The coal seam is homogeneous, porous, and isotropic. In-situ methane pressure and methane content are uniform everywhere in the coal seam(3)
The initial volume fraction of methane in the coal seam is 85%, and the water and air saturation are 10% and 5%, respectively(4)
Temperature of the coal seam is constant (293.15 K), and there is no heat exchange between the boundary and the coal seam(5)
Methane in the coal seam follows the ideal gas equation


### 2.1. Model Parameters

The coal seam model is 300×200×10 m. The plan view is shown in Figure 1. A single borehole model is established. It is a 0.1 m borehole located in the center of calculated zone (X =150 m, Y =100 m). Eight detection points (A to H) are arranged at different distances away from injection borehole along line L. As shown in Figure 2, the boundary pressure of coal seam is set to 2 MPa. Coal seam parameters are shown in Table 1.Figure 1
Detection points arrangement of the computational model. Line L: a virtual straight line along which the eight detection points (A to H) are arranged.Figure 2
The 3D computational model.Table 1
Coal seam parameters.
SymbolParameterValueUnitpcCoal density1250Kg/m3∅Initial porosity0.031κInitial permeability3.8e-16m2ρmMethane standard density0.7Kg/m3ρaAir density1.293Kg/m3ρwWater density0.9814Kg/m3μmMethane viscosity1.84e-5Pa·sμaAir viscosity1.79e-5Pa·sμwWater viscosity0.96e-3Pa·spinInjection pressure15MPapoutBoundary pressure2MPapmInitial seam pressure2MPatReference temperature20°CswInitial water saturation0.11saInitial air saturation0.051smInitial methane volume fraction0.851


### 2.2. Multiphase Flow Model Equation

#### 2.2.1. Calculation Model of Fluid Seepage in Porous Media

The multiphase flow in coal seams follows Darcy’s law with the continuity equation [21, 22]:
(1)∂∂tml+∇∙ρlul=Qlswhere ρl is density of fluid, kg/m3; ul is seepage velocity of each phase, m/s, Qls is mass source of the fluid, kg/(m3·s). ml is liquid mass in the coal matrix per unit volume, which can be calculated by:
(2)ml=ρl∙ϕIn the present work, the swelling and deformation of the coal matrix due to water absorption are not considered, so the porosity can be expressed as:(3)ϕ=1+R0ϕ0+αR−R01+Rwhere R=εv+p/Ks;R0=εv0+p0/Ks;εv and εv0 are volumetric strain and initial volumetric strain of coal seam, respectively; Ks is bulk modulus of coal skeleton.The equation of permeability and porosity obtained by simultaneous formula is as follows:(4)klkl0=ϕϕ03=1+R0+αR−R0/ϕ01+R3where kf and kf0 are the permeability and initial permeability of coal seam, respectively; ϕ and ϕ0 are the porosity and initial porosity of coal seam, respectively.Darcy’s law for seepage in porous media is as follows:(5)qi=−kμi∇p+ρig∇Dwhere qi is seepage velocity of pressure water, m/s; μi is hydrodynamic viscosity, Pa·s; k is permeability, mD; ρiis density of pressure water, kg/m3,∇D is gravitational potential difference, m; g is gravitational acceleration, 9.8 m2/s.(1) and (3) can be derived as
(6)∂ϕ∂t=α−ϕ1+R∂εv∂t+1Ks·∂p∂tρl∂ϕ∂t+ϕ∂ρl∂t+∇ρlql=Qls


#### 2.2.2. Multiphase Fluid Migration Model

Coal seam be regarded as a continuous porous medium with fully developed pores and fractures. The multiphase fluid seepage follows Darcy’s seepage law[23]:
(7)u=−κriμiκ∇piwhere u is fluid velocity, m/s; μi is dynamic viscosity of the fluid, Pa·s; κ is effective permeability, mD; κri is relative permeability, mD; ∇pi is pressure gradient (Pa/m).The effective permeability (κ) is related to relative permeability (κri)and can be calculated by:
(8)κri=κiκTotal methane including adsorbed methane and free methane in the coal matrix can be calculated by:(9)q=VLqmpL+pm∗ρMVm+∅mMpmRTwhere q is mass of methane per unit volume in coal matrix, kg/m3; ρ is coal density, kg/m3; VL is Langmuir volume (maximum adsorption capacity of coal), m3/t;pL is Langmuir pressure, MPa; ∅m is porosity of the coal matrix, %;Vm is molar volume of methane under standard conditions, m3/mol [24].Multiphase fluid transport is based on macroscopic mass conservation equation of each phase. The average volume fraction (also known as saturation) of each phase is found. Describe single-phase fluid transport:(10)∂∂tεiρisi+∇∙ρiui=Qiwhere εi (dimensionless) is porosity; ρi is phase density, kg/m3; si is volume fraction (dimensionless); ui represents the velocity vector of phase i, m/s; Qi is mass source of phase i, kg/(m3·s).Assuming that the sum of the volume fractions of each phase is 1 (100%), the volume fraction of the remaining phaseSj is obtained by:
(11)Sj=1−∑i=1.i≠jNSiwhere Si is volume fraction of phase i; Sj is volume fraction of phase j.Combining (8) and (11), other N-1 phase fluid motion equations can be obtained:
(12)∂∂tεiρisi−∇·ρiκκriμi∇pi=Qi


## 2.1. Model Parameters

The coal seam model is 300×200×10 m. The plan view is shown in Figure 1. A single borehole model is established. It is a 0.1 m borehole located in the center of calculated zone (X =150 m, Y =100 m). Eight detection points (A to H) are arranged at different distances away from injection borehole along line L. As shown in Figure 2, the boundary pressure of coal seam is set to 2 MPa. Coal seam parameters are shown in Table 1.Figure 1
Detection points arrangement of the computational model. Line L: a virtual straight line along which the eight detection points (A to H) are arranged.Figure 2
The 3D computational model.Table 1
Coal seam parameters.
SymbolParameterValueUnitpcCoal density1250Kg/m3∅Initial porosity0.031κInitial permeability3.8e-16m2ρmMethane standard density0.7Kg/m3ρaAir density1.293Kg/m3ρwWater density0.9814Kg/m3μmMethane viscosity1.84e-5Pa·sμaAir viscosity1.79e-5Pa·sμwWater viscosity0.96e-3Pa·spinInjection pressure15MPapoutBoundary pressure2MPapmInitial seam pressure2MPatReference temperature20°CswInitial water saturation0.11saInitial air saturation0.051smInitial methane volume fraction0.851


## 2.2. Multiphase Flow Model Equation

### 2.2.1. Calculation Model of Fluid Seepage in Porous Media

The multiphase flow in coal seams follows Darcy’s law with the continuity equation [21, 22]:
(1)∂∂tml+∇∙ρlul=Qlswhere ρl is density of fluid, kg/m3; ul is seepage velocity of each phase, m/s, Qls is mass source of the fluid, kg/(m3·s). ml is liquid mass in the coal matrix per unit volume, which can be calculated by:
(2)ml=ρl∙ϕIn the present work, the swelling and deformation of the coal matrix due to water absorption are not considered, so the porosity can be expressed as:(3)ϕ=1+R0ϕ0+αR−R01+Rwhere R=εv+p/Ks;R0=εv0+p0/Ks;εv and εv0 are volumetric strain and initial volumetric strain of coal seam, respectively; Ks is bulk modulus of coal skeleton.The equation of permeability and porosity obtained by simultaneous formula is as follows:(4)klkl0=ϕϕ03=1+R0+αR−R0/ϕ01+R3where kf and kf0 are the permeability and initial permeability of coal seam, respectively; ϕ and ϕ0 are the porosity and initial porosity of coal seam, respectively.Darcy’s law for seepage in porous media is as follows:(5)qi=−kμi∇p+ρig∇Dwhere qi is seepage velocity of pressure water, m/s; μi is hydrodynamic viscosity, Pa·s; k is permeability, mD; ρiis density of pressure water, kg/m3,∇D is gravitational potential difference, m; g is gravitational acceleration, 9.8 m2/s.(1) and (3) can be derived as
(6)∂ϕ∂t=α−ϕ1+R∂εv∂t+1Ks·∂p∂tρl∂ϕ∂t+ϕ∂ρl∂t+∇ρlql=Qls


### 2.2.2. Multiphase Fluid Migration Model

Coal seam be regarded as a continuous porous medium with fully developed pores and fractures. The multiphase fluid seepage follows Darcy’s seepage law[23]:
(7)u=−κriμiκ∇piwhere u is fluid velocity, m/s; μi is dynamic viscosity of the fluid, Pa·s; κ is effective permeability, mD; κri is relative permeability, mD; ∇pi is pressure gradient (Pa/m).The effective permeability (κ) is related to relative permeability (κri)and can be calculated by:
(8)κri=κiκTotal methane including adsorbed methane and free methane in the coal matrix can be calculated by:(9)q=VLqmpL+pm∗ρMVm+∅mMpmRTwhere q is mass of methane per unit volume in coal matrix, kg/m3; ρ is coal density, kg/m3; VL is Langmuir volume (maximum adsorption capacity of coal), m3/t;pL is Langmuir pressure, MPa; ∅m is porosity of the coal matrix, %;Vm is molar volume of methane under standard conditions, m3/mol [24].Multiphase fluid transport is based on macroscopic mass conservation equation of each phase. The average volume fraction (also known as saturation) of each phase is found. Describe single-phase fluid transport:(10)∂∂tεiρisi+∇∙ρiui=Qiwhere εi (dimensionless) is porosity; ρi is phase density, kg/m3; si is volume fraction (dimensionless); ui represents the velocity vector of phase i, m/s; Qi is mass source of phase i, kg/(m3·s).Assuming that the sum of the volume fractions of each phase is 1 (100%), the volume fraction of the remaining phaseSj is obtained by:
(11)Sj=1−∑i=1.i≠jNSiwhere Si is volume fraction of phase i; Sj is volume fraction of phase j.Combining (8) and (11), other N-1 phase fluid motion equations can be obtained:
(12)∂∂tεiρisi−∇·ρiκκriμi∇pi=Qi


## 2.2.1. Calculation Model of Fluid Seepage in Porous Media

The multiphase flow in coal seams follows Darcy’s law with the continuity equation [21, 22]:
(1)∂∂tml+∇∙ρlul=Qlswhere ρl is density of fluid, kg/m3; ul is seepage velocity of each phase, m/s, Qls is mass source of the fluid, kg/(m3·s). ml is liquid mass in the coal matrix per unit volume, which can be calculated by:
(2)ml=ρl∙ϕIn the present work, the swelling and deformation of the coal matrix due to water absorption are not considered, so the porosity can be expressed as:(3)ϕ=1+R0ϕ0+αR−R01+Rwhere R=εv+p/Ks;R0=εv0+p0/Ks;εv and εv0 are volumetric strain and initial volumetric strain of coal seam, respectively; Ks is bulk modulus of coal skeleton.The equation of permeability and porosity obtained by simultaneous formula is as follows:(4)klkl0=ϕϕ03=1+R0+αR−R0/ϕ01+R3where kf and kf0 are the permeability and initial permeability of coal seam, respectively; ϕ and ϕ0 are the porosity and initial porosity of coal seam, respectively.Darcy’s law for seepage in porous media is as follows:(5)qi=−kμi∇p+ρig∇Dwhere qi is seepage velocity of pressure water, m/s; μi is hydrodynamic viscosity, Pa·s; k is permeability, mD; ρiis density of pressure water, kg/m3,∇D is gravitational potential difference, m; g is gravitational acceleration, 9.8 m2/s.(1) and (3) can be derived as
(6)∂ϕ∂t=α−ϕ1+R∂εv∂t+1Ks·∂p∂tρl∂ϕ∂t+ϕ∂ρl∂t+∇ρlql=Qls


## 2.2.2. Multiphase Fluid Migration Model

Coal seam be regarded as a continuous porous medium with fully developed pores and fractures. The multiphase fluid seepage follows Darcy’s seepage law[23]:
(7)u=−κriμiκ∇piwhere u is fluid velocity, m/s; μi is dynamic viscosity of the fluid, Pa·s; κ is effective permeability, mD; κri is relative permeability, mD; ∇pi is pressure gradient (Pa/m).The effective permeability (κ) is related to relative permeability (κri)and can be calculated by:
(8)κri=κiκTotal methane including adsorbed methane and free methane in the coal matrix can be calculated by:(9)q=VLqmpL+pm∗ρMVm+∅mMpmRTwhere q is mass of methane per unit volume in coal matrix, kg/m3; ρ is coal density, kg/m3; VL is Langmuir volume (maximum adsorption capacity of coal), m3/t;pL is Langmuir pressure, MPa; ∅m is porosity of the coal matrix, %;Vm is molar volume of methane under standard conditions, m3/mol [24].Multiphase fluid transport is based on macroscopic mass conservation equation of each phase. The average volume fraction (also known as saturation) of each phase is found. Describe single-phase fluid transport:(10)∂∂tεiρisi+∇∙ρiui=Qiwhere εi (dimensionless) is porosity; ρi is phase density, kg/m3; si is volume fraction (dimensionless); ui represents the velocity vector of phase i, m/s; Qi is mass source of phase i, kg/(m3·s).Assuming that the sum of the volume fractions of each phase is 1 (100%), the volume fraction of the remaining phaseSj is obtained by:
(11)Sj=1−∑i=1.i≠jNSiwhere Si is volume fraction of phase i; Sj is volume fraction of phase j.Combining (8) and (11), other N-1 phase fluid motion equations can be obtained:
(12)∂∂tεiρisi−∇·ρiκκriμi∇pi=Qi


## 3. Results and Discussion

### 3.1. Comparison of Injecting Different Fluids

#### 3.1.1. Water Injection Displacement

Figure3 shows the evolution of methane volume fraction recorded every 2 hours. The injected water expanded from the borehole to far-field area of coal seam. The influence range of water injection expanded with the increase of injection time. As shown in Figure 4, methane fraction volume (VCH4) near the injection borehole (point A, X =3 m) decreased significantly from 85% to 73.31% after 2 hours and decreased to 70.50% after 22 hours. After 24 hours water injection, VCH4 at point F to H still fluctuated around 85% and showed a slight upward trend, indicating a weak disturb by water injection. The influence effective radius of water injection is about 13 m.Figure 3
Contour ofVCH4 under water injection.Figure 4
Time-dependentVCH4 at each detection point. X =3 m: 3 m away horizontally from injection borehole. Avge: The average VCH4 of all detection points at the end of the displacement (24 h).VCH4 at all detection points increased at the initial stage (less than 5%) and then decreased sharply. The phenomenon also appears in the experimental studies of other scholars [11]. A reasonable explanation is: under water injection, affected by high injection pressure, pressure gradient forces the methane carried by water and the methane V1 displaced by water to passively flow to the outside, plus the external coalbed methane V2 that has not been affected by hydraulic permeability. Both V1 and V2 overlap, and the VCH4 increases briefly and in a small range.Figure5 gives time-dependent pore pressures at different detection points. Pore pressures increased rapidly in the first two hours under high-pressure water injection, which indicates that pore pressures quickly covered all detection points under water injection. Pore pressure is higher in areas closer to the injection hole, which will be beneficial to fracture the coal body and extend the primary fractures.Figure 5
Time-dependent pore pressures at each detection point. Min: Minimum pore pressure in all detection points at the end of displacement.


#### 3.1.2. Air Injection Displacement

As shown in Figure6, the reduction of VCH4 is lower than that of water displacement. The injected air expanded faster than water and had a wider displaced area. As shown in Figure 7 under air injection displacement, the decrease of VCH4 fluctuated little at all detection points. VCH4 decreased sharply within 2 h. The total VCH4 at detection point A decreased by 90% at the first 2 h, indicating a high displacement velocity. The air expansion in coal seam is relatively uniform.Figure 6
Contour ofVCH4 under air injection.Figure 7
Time-dependentVCH4 of each detection point.The pore pressures at different detection points finally reached a stable value when air was continuously injected. The average stable pore pressure of different detection points is 8.85 MPa (as shown in Figure8), which is slightly higher than 8.05 MPa under water injection. As shown in Figure 9, the pore pressures and influence range under air injection are both significantly larger than water injection at the same injection time (12 h) along line L. The velocity of pressure propagation is obviously faster than that of water injection.Figure 8
Time-dependent pore pressures at each detection point.Figure 9
Pore pressures under water injection and air injection for 12 hours, respectively, along line L (Orange dashed line shown in Figure1). Data is recorded and plotted every 2 hours.


#### 3.1.3. Characteristics of Water-Air Alternating Displacement

Water was injected into coal seam in the first 12 hours and air in the next 12 hours. As shown in Figure10, the VCH4 near the borehole decreased significantly under air injection. And the influence range was larger compared with water injection at the same time. As shown in Figure 11, VCH4 at detection points A, B, and C decreased twice and stabilized at about 43%, while the average stable VCH4 of other detection points was 58.5%. Under air injection (14-24 h), VCH4 at different detection points did not change much, indicating that the air injection completed the displacement within two hours. This phenomenon is consistent with air injection. Therefore, WGA-ECBM technique is combined with the advantages of air injection. Under the water-air altering displacement, the average VCH4 of all detection points decreased to 52.9%, while that of water injection and air injection under the same conditions decreased to 79.3% and 74.9%, respectively.Figure 10
Contour ofVCH4 under water-air alternating displacement.Figure 11
Time-dependentVCH4 at different detection points.Under the water-air alternating displacement, pore pressures first increased under water injection. Pore pressures firstly decreased and then increased under air injection, as shown in Figure12. At the time of air injection (t =12-13 h), pore pressures decreased because injected air drives the multiphase flow to migrate. According to Bernoulli equation: p+ρgz+1/2∗ρv^2=C, part of pressure potential energy is converted to kinetic energy. And the other part is lost when pressure gradient fractures coal body to expand the cracks and break through the water-blocked migration channel. Figure 13 shows the 3D distribution of pore pressures in coal seam. The decrease rate of pore pressure gradually slows down from borehole to the boundary.Figure 12
Time-dependent pore pressures at different detection points.Figure 13
Time-dependent pore pressures under water-air alternating displacement.


#### 3.1.4. Comprehensive Comparison of Three Displacement Methods

Figure14 shows VCH4 reduction at each detection point under three displacement methods. VCH4 decreases most significantly under water-air alternating displacement. VCH4 at detection points (F, G, H) do not decrease under water injection. However, VCH4 reduction at detection points (A, B, C) is higher under water injection than under air injection.Figure 14
VCH4 reduction at each detection point under three displacement methods.As shown in Figure15, VCH4 in the borehole decreased to 67.16% and 72.28% under water injection and air injection, respectively. The water-air alternating displacement decreased to 40.38%. VCH4 decreased again by 26.94% based on 17.84% reduction of water injection. After water injection, the aid of air injection on displacement is significant, and the second decrease in methane volume is greater. As shown in Figure 15, the influence ranges are 24.68 m, 149.73 m and 135.21 m (the diameter on L line) under water injection, air injection and water-air alternating displacement, respectively. As shown in Figure 16, pore pressure distribution under three displacement methods is similar at 24 h. This indicates that the injection pressure has completed a relatively uniform spread within 24 h.Figure 15
VCH4 under three displacement methods along line L shown in Figure 1.Figure 16
Pore pressures under three displacement methods at 24 h along line L.


### 3.2. Water-Air Alternating Displacement under Different Injection Pressures

The injection pressure affecting pore pressure and the permeability in coal seam is an important engineering parameter. As shown in Figure17, the pore pressure in the contour appears to be reduced at t =12 h compared to t =8 h. The injected air pushes the methane-water in the pore space to seep outward and overcomes the capillary pressure. Resolving water-lock damage causes the pore pressures to temporarily decrease and then increase.Figure 17
Contour of pore pressures under different injection pressures.As shown in Figure18, VCH4 decreases gradually as injection pressure increases in a pressure gradient of 5 MPa. However, VCH4 decreases slightly under water displacement. The influence range under air injection is greater than that under water injection.Figure 18
Contour ofVCH4 under different injection pressures.As shown in Figure19, at the end of water injection (t =8 h), the influence ranges from 10 MPa to 25 MPa are 10.50 m, 15.88 m, 19.82 m, and 22.27 m, respectively. VCH4 in the borehole are 75.06%, 66.79%, 64.01%, and 62.36%, respectively.Figure 19
VCH4 at the end of water injection (8 h) along line L.As shown in Figure20,with the foundation of the water injection, the displacement effect of the air injection is obviously better. VCH4 in the borehole are 49.43%, 41.18%, 38.22% and 33.28%, respectively. The pore pressure increases as the injection pressure increases. Pressure gradient contributes to seepage (Darcy’s law). The influence ranges from 10 MPa to 25 MPa are 80.92 m, 113.26 m, 135.18 m, and 154.54 m, respectively. Compared with the displacement characteristics in Figure 19, the influence range is obviously expanded.Figure 20
VCH4 at the end of air injection (16 h) following water injection.


### 3.3. Water-Air Alternating Displacement under Different Injection Durations

As shown in Figure21, the influence range increases continuously with the increasing injection duration. The area expansion of the contours is significantly larger under air injection than under water injection. The air injection is more dramatically influenced by injection duration.Figure 21
Contour ofVCH4 under different injection durations.As shown in Figure22, As injection duration increases in the WGA-ECBM, VCH4 in the displaced area is less and less affected by the injection duration. VCH4 in the borehole are 42.9%, 41.2%, 40.4%, 40.3% and 39.3%, respectively. Because the adsorption and desorption of coal matrix tend to be balanced after a period of time. The injection duration ranges from 8 to 96 h, and the influence ranges are 77.72 m, 110.43 m, 135.18 m, 188.46 m and 270.90 m, respectively. However, As the wetting area of coal seam becomes larger, the seepage channels increase, and the displacement velocity gradually slows down, as shown in Figure 23. The injection duration has little effect on pore pressures. The pore pressure can propagate to the model boundary at t =8 h in Figure 24. The pore pressure distribution characteristics are similar under different injection durations.Figure 22
VCH4 under water-air alternating displacement along line L.Figure 23
Average methane displacement velocity.Figure 24
Pore pressures under water-air alternating displacement along line L.


## 3.1. Comparison of Injecting Different Fluids

### 3.1.1. Water Injection Displacement

Figure3 shows the evolution of methane volume fraction recorded every 2 hours. The injected water expanded from the borehole to far-field area of coal seam. The influence range of water injection expanded with the increase of injection time. As shown in Figure 4, methane fraction volume (VCH4) near the injection borehole (point A, X =3 m) decreased significantly from 85% to 73.31% after 2 hours and decreased to 70.50% after 22 hours. After 24 hours water injection, VCH4 at point F to H still fluctuated around 85% and showed a slight upward trend, indicating a weak disturb by water injection. The influence effective radius of water injection is about 13 m.Figure 3
Contour ofVCH4 under water injection.Figure 4
Time-dependentVCH4 at each detection point. X =3 m: 3 m away horizontally from injection borehole. Avge: The average VCH4 of all detection points at the end of the displacement (24 h).VCH4 at all detection points increased at the initial stage (less than 5%) and then decreased sharply. The phenomenon also appears in the experimental studies of other scholars [11]. A reasonable explanation is: under water injection, affected by high injection pressure, pressure gradient forces the methane carried by water and the methane V1 displaced by water to passively flow to the outside, plus the external coalbed methane V2 that has not been affected by hydraulic permeability. Both V1 and V2 overlap, and the VCH4 increases briefly and in a small range.Figure5 gives time-dependent pore pressures at different detection points. Pore pressures increased rapidly in the first two hours under high-pressure water injection, which indicates that pore pressures quickly covered all detection points under water injection. Pore pressure is higher in areas closer to the injection hole, which will be beneficial to fracture the coal body and extend the primary fractures.Figure 5
Time-dependent pore pressures at each detection point. Min: Minimum pore pressure in all detection points at the end of displacement.


### 3.1.2. Air Injection Displacement

As shown in Figure6, the reduction of VCH4 is lower than that of water displacement. The injected air expanded faster than water and had a wider displaced area. As shown in Figure 7 under air injection displacement, the decrease of VCH4 fluctuated little at all detection points. VCH4 decreased sharply within 2 h. The total VCH4 at detection point A decreased by 90% at the first 2 h, indicating a high displacement velocity. The air expansion in coal seam is relatively uniform.Figure 6
Contour ofVCH4 under air injection.Figure 7
Time-dependentVCH4 of each detection point.The pore pressures at different detection points finally reached a stable value when air was continuously injected. The average stable pore pressure of different detection points is 8.85 MPa (as shown in Figure8), which is slightly higher than 8.05 MPa under water injection. As shown in Figure 9, the pore pressures and influence range under air injection are both significantly larger than water injection at the same injection time (12 h) along line L. The velocity of pressure propagation is obviously faster than that of water injection.Figure 8
Time-dependent pore pressures at each detection point.Figure 9
Pore pressures under water injection and air injection for 12 hours, respectively, along line L (Orange dashed line shown in Figure1). Data is recorded and plotted every 2 hours.


### 3.1.3. Characteristics of Water-Air Alternating Displacement

Water was injected into coal seam in the first 12 hours and air in the next 12 hours. As shown in Figure10, the VCH4 near the borehole decreased significantly under air injection. And the influence range was larger compared with water injection at the same time. As shown in Figure 11, VCH4 at detection points A, B, and C decreased twice and stabilized at about 43%, while the average stable VCH4 of other detection points was 58.5%. Under air injection (14-24 h), VCH4 at different detection points did not change much, indicating that the air injection completed the displacement within two hours. This phenomenon is consistent with air injection. Therefore, WGA-ECBM technique is combined with the advantages of air injection. Under the water-air altering displacement, the average VCH4 of all detection points decreased to 52.9%, while that of water injection and air injection under the same conditions decreased to 79.3% and 74.9%, respectively.Figure 10
Contour ofVCH4 under water-air alternating displacement.Figure 11
Time-dependentVCH4 at different detection points.Under the water-air alternating displacement, pore pressures first increased under water injection. Pore pressures firstly decreased and then increased under air injection, as shown in Figure12. At the time of air injection (t =12-13 h), pore pressures decreased because injected air drives the multiphase flow to migrate. According to Bernoulli equation: p+ρgz+1/2∗ρv^2=C, part of pressure potential energy is converted to kinetic energy. And the other part is lost when pressure gradient fractures coal body to expand the cracks and break through the water-blocked migration channel. Figure 13 shows the 3D distribution of pore pressures in coal seam. The decrease rate of pore pressure gradually slows down from borehole to the boundary.Figure 12
Time-dependent pore pressures at different detection points.Figure 13
Time-dependent pore pressures under water-air alternating displacement.


### 3.1.4. Comprehensive Comparison of Three Displacement Methods

Figure14 shows VCH4 reduction at each detection point under three displacement methods. VCH4 decreases most significantly under water-air alternating displacement. VCH4 at detection points (F, G, H) do not decrease under water injection. However, VCH4 reduction at detection points (A, B, C) is higher under water injection than under air injection.Figure 14
VCH4 reduction at each detection point under three displacement methods.As shown in Figure15, VCH4 in the borehole decreased to 67.16% and 72.28% under water injection and air injection, respectively. The water-air alternating displacement decreased to 40.38%. VCH4 decreased again by 26.94% based on 17.84% reduction of water injection. After water injection, the aid of air injection on displacement is significant, and the second decrease in methane volume is greater. As shown in Figure 15, the influence ranges are 24.68 m, 149.73 m and 135.21 m (the diameter on L line) under water injection, air injection and water-air alternating displacement, respectively. As shown in Figure 16, pore pressure distribution under three displacement methods is similar at 24 h. This indicates that the injection pressure has completed a relatively uniform spread within 24 h.Figure 15
VCH4 under three displacement methods along line L shown in Figure 1.Figure 16
Pore pressures under three displacement methods at 24 h along line L.


## 3.1.1. Water Injection Displacement

Figure3 shows the evolution of methane volume fraction recorded every 2 hours. The injected water expanded from the borehole to far-field area of coal seam. The influence range of water injection expanded with the increase of injection time. As shown in Figure 4, methane fraction volume (VCH4) near the injection borehole (point A, X =3 m) decreased significantly from 85% to 73.31% after 2 hours and decreased to 70.50% after 22 hours. After 24 hours water injection, VCH4 at point F to H still fluctuated around 85% and showed a slight upward trend, indicating a weak disturb by water injection. The influence effective radius of water injection is about 13 m.Figure 3
Contour ofVCH4 under water injection.Figure 4
Time-dependentVCH4 at each detection point. X =3 m: 3 m away horizontally from injection borehole. Avge: The average VCH4 of all detection points at the end of the displacement (24 h).VCH4 at all detection points increased at the initial stage (less than 5%) and then decreased sharply. The phenomenon also appears in the experimental studies of other scholars [11]. A reasonable explanation is: under water injection, affected by high injection pressure, pressure gradient forces the methane carried by water and the methane V1 displaced by water to passively flow to the outside, plus the external coalbed methane V2 that has not been affected by hydraulic permeability. Both V1 and V2 overlap, and the VCH4 increases briefly and in a small range.Figure5 gives time-dependent pore pressures at different detection points. Pore pressures increased rapidly in the first two hours under high-pressure water injection, which indicates that pore pressures quickly covered all detection points under water injection. Pore pressure is higher in areas closer to the injection hole, which will be beneficial to fracture the coal body and extend the primary fractures.Figure 5
Time-dependent pore pressures at each detection point. Min: Minimum pore pressure in all detection points at the end of displacement.


## 3.1.2. Air Injection Displacement

As shown in Figure6, the reduction of VCH4 is lower than that of water displacement. The injected air expanded faster than water and had a wider displaced area. As shown in Figure 7 under air injection displacement, the decrease of VCH4 fluctuated little at all detection points. VCH4 decreased sharply within 2 h. The total VCH4 at detection point A decreased by 90% at the first 2 h, indicating a high displacement velocity. The air expansion in coal seam is relatively uniform.Figure 6
Contour ofVCH4 under air injection.Figure 7
Time-dependentVCH4 of each detection point.The pore pressures at different detection points finally reached a stable value when air was continuously injected. The average stable pore pressure of different detection points is 8.85 MPa (as shown in Figure8), which is slightly higher than 8.05 MPa under water injection. As shown in Figure 9, the pore pressures and influence range under air injection are both significantly larger than water injection at the same injection time (12 h) along line L. The velocity of pressure propagation is obviously faster than that of water injection.Figure 8
Time-dependent pore pressures at each detection point.Figure 9
Pore pressures under water injection and air injection for 12 hours, respectively, along line L (Orange dashed line shown in Figure1). Data is recorded and plotted every 2 hours.


## 3.1.3. Characteristics of Water-Air Alternating Displacement

Water was injected into coal seam in the first 12 hours and air in the next 12 hours. As shown in Figure10, the VCH4 near the borehole decreased significantly under air injection. And the influence range was larger compared with water injection at the same time. As shown in Figure 11, VCH4 at detection points A, B, and C decreased twice and stabilized at about 43%, while the average stable VCH4 of other detection points was 58.5%. Under air injection (14-24 h), VCH4 at different detection points did not change much, indicating that the air injection completed the displacement within two hours. This phenomenon is consistent with air injection. Therefore, WGA-ECBM technique is combined with the advantages of air injection. Under the water-air altering displacement, the average VCH4 of all detection points decreased to 52.9%, while that of water injection and air injection under the same conditions decreased to 79.3% and 74.9%, respectively.Figure 10
Contour ofVCH4 under water-air alternating displacement.Figure 11
Time-dependentVCH4 at different detection points.Under the water-air alternating displacement, pore pressures first increased under water injection. Pore pressures firstly decreased and then increased under air injection, as shown in Figure12. At the time of air injection (t =12-13 h), pore pressures decreased because injected air drives the multiphase flow to migrate. According to Bernoulli equation: p+ρgz+1/2∗ρv^2=C, part of pressure potential energy is converted to kinetic energy. And the other part is lost when pressure gradient fractures coal body to expand the cracks and break through the water-blocked migration channel. Figure 13 shows the 3D distribution of pore pressures in coal seam. The decrease rate of pore pressure gradually slows down from borehole to the boundary.Figure 12
Time-dependent pore pressures at different detection points.Figure 13
Time-dependent pore pressures under water-air alternating displacement.


## 3.1.4. Comprehensive Comparison of Three Displacement Methods

Figure14 shows VCH4 reduction at each detection point under three displacement methods. VCH4 decreases most significantly under water-air alternating displacement. VCH4 at detection points (F, G, H) do not decrease under water injection. However, VCH4 reduction at detection points (A, B, C) is higher under water injection than under air injection.Figure 14
VCH4 reduction at each detection point under three displacement methods.As shown in Figure15, VCH4 in the borehole decreased to 67.16% and 72.28% under water injection and air injection, respectively. The water-air alternating displacement decreased to 40.38%. VCH4 decreased again by 26.94% based on 17.84% reduction of water injection. After water injection, the aid of air injection on displacement is significant, and the second decrease in methane volume is greater. As shown in Figure 15, the influence ranges are 24.68 m, 149.73 m and 135.21 m (the diameter on L line) under water injection, air injection and water-air alternating displacement, respectively. As shown in Figure 16, pore pressure distribution under three displacement methods is similar at 24 h. This indicates that the injection pressure has completed a relatively uniform spread within 24 h.Figure 15
VCH4 under three displacement methods along line L shown in Figure 1.Figure 16
Pore pressures under three displacement methods at 24 h along line L.


## 3.2. Water-Air Alternating Displacement under Different Injection Pressures

The injection pressure affecting pore pressure and the permeability in coal seam is an important engineering parameter. As shown in Figure17, the pore pressure in the contour appears to be reduced at t =12 h compared to t =8 h. The injected air pushes the methane-water in the pore space to seep outward and overcomes the capillary pressure. Resolving water-lock damage causes the pore pressures to temporarily decrease and then increase.Figure 17
Contour of pore pressures under different injection pressures.As shown in Figure18, VCH4 decreases gradually as injection pressure increases in a pressure gradient of 5 MPa. However, VCH4 decreases slightly under water displacement. The influence range under air injection is greater than that under water injection.Figure 18
Contour ofVCH4 under different injection pressures.As shown in Figure19, at the end of water injection (t =8 h), the influence ranges from 10 MPa to 25 MPa are 10.50 m, 15.88 m, 19.82 m, and 22.27 m, respectively. VCH4 in the borehole are 75.06%, 66.79%, 64.01%, and 62.36%, respectively.Figure 19
VCH4 at the end of water injection (8 h) along line L.As shown in Figure20,with the foundation of the water injection, the displacement effect of the air injection is obviously better. VCH4 in the borehole are 49.43%, 41.18%, 38.22% and 33.28%, respectively. The pore pressure increases as the injection pressure increases. Pressure gradient contributes to seepage (Darcy’s law). The influence ranges from 10 MPa to 25 MPa are 80.92 m, 113.26 m, 135.18 m, and 154.54 m, respectively. Compared with the displacement characteristics in Figure 19, the influence range is obviously expanded.Figure 20
VCH4 at the end of air injection (16 h) following water injection.


## 3.3. Water-Air Alternating Displacement under Different Injection Durations

As shown in Figure21, the influence range increases continuously with the increasing injection duration. The area expansion of the contours is significantly larger under air injection than under water injection. The air injection is more dramatically influenced by injection duration.Figure 21
Contour ofVCH4 under different injection durations.As shown in Figure22, As injection duration increases in the WGA-ECBM, VCH4 in the displaced area is less and less affected by the injection duration. VCH4 in the borehole are 42.9%, 41.2%, 40.4%, 40.3% and 39.3%, respectively. Because the adsorption and desorption of coal matrix tend to be balanced after a period of time. The injection duration ranges from 8 to 96 h, and the influence ranges are 77.72 m, 110.43 m, 135.18 m, 188.46 m and 270.90 m, respectively. However, As the wetting area of coal seam becomes larger, the seepage channels increase, and the displacement velocity gradually slows down, as shown in Figure 23. The injection duration has little effect on pore pressures. The pore pressure can propagate to the model boundary at t =8 h in Figure 24. The pore pressure distribution characteristics are similar under different injection durations.Figure 22
VCH4 under water-air alternating displacement along line L.Figure 23
Average methane displacement velocity.Figure 24
Pore pressures under water-air alternating displacement along line L.


## 4. Conclusions

The WGA-ECBM provide a new high-efficiency technique for the methane recovery. And in the multiphase fluid-solid coupling model, the influences of different parameters on fluid migration, pore pressure and methane content are discussed, which has important guidance value in field practice.(1)
The WGA-ECBM eliminated the high methane area generated by water injection while it also helped water injection to further expand the influence range. The influence ranges (diameter) are 24.68 m, 149.73 m and 135.21 m under water injection, air injection and water-air alternating injection, respectively. Under the same conditions, the pore pressure distribution characteristics of three methods are similar at 24 h(2)
Water-air alternating injection can quickly reduceVCH4 again based on methane reduction under water injection. VCH4 in the borehole are 67.16%, 72.28% and 40.38% under water injection, air injection and water-air alternating injection, respectively. The methane reduction effect of WGA-ECBM is the best among three methods(3)
The influence range of air injection is wider than water injection. The curve fluctuation ofVCH4 under air injection was smaller at each detection point. However, water injection is better than air injection in VCH4 reduction.(4)
The influence of injection pressure and injection duration under air injection is greater than under water injection. High pressure water injection for long time is not effective. In the WGA-ECBM technique,VCH4 decreases and influence range expands as the injection pressure increases in a pressure gradient of 5 MPa. With the increase of injection duration, the influence range expands significantly but the pore pressure distribution changes little. And VCH4 decreases slightly in the displacement area and the displacement velocity gradually slows down



---

*Source: 4563438-2022-11-30.xml*