Sage Journals: Discover world-class research

Abstract

Water in coal reservoirs serves as both the driving force of reservoir fluid and the resistance to coalbed methane (CBM) migration. However, the impact and control mechanism of water saturation on CBM development in the field have not been revealed yet. The evolution of reservoir dynamics and gas production rate under various initial water saturation conditions is simulated based on the COMSOL multiphysics numerical simulation software. The simulation results show that the reduction of initial water saturation can significantly enhance the effective permeability of the reservoir gas phase, which can be increased by up to 4764%. This increases the rate of gas production, which is the main mechanism for controlling production by water saturation in the reservoir. For high water saturation reservoirs, the positive effect of inter-well interference should be reasonably utilized, and high well network density should be used to rapidly increase gas-effective permeability. For low water saturation reservoirs, the negative effect of inter-well interference should be avoided, and low well network density should be used to increase the control of single well reserves.

Keywords

Coal reservoir water saturation coalbed methane development permeability evolution numerical simulation

Introduction

The water saturation of coal reservoirs is an important factor in the development of coalbed methane (CBM), which largely determines the flow and dynamic relative permeability of fluids in reservoirs. In coal reservoirs, formation water stores elastic potential energy (Zhao et al., 2022a), which drives the migration of reservoir fluids. Meanwhile, because gas and water share the flow channels in the reservoir, the formation of water also acts as a resistance to the migration of CBM (Bear, 1975). The saturation of reservoir water controls the migration and production of CBM (Akhondzadeh et al., 2018).

In the early stage of CBM exploration and development, it is generally recognized in the academic and engineering communities that the water in coal reservoirs is saturated (Chen et al., 2018). Only a few coal reservoirs (Dry Horseshoe Canon CBM Play) show low water saturation or are completely anhydrous, and good development results are achieved under low-permeability conditions (Clarkson, 2009; Xu et al., 2015b). It is also believed that most CBM reservoirs are undersaturated, meaning that the actual gas content is lower than the theoretical adsorption content under reservoir pressure (Xu et al., 2014). As a result, fluid flow in reservoirs begins with single-phase water flow, and as reservoir pressure decreases and reaches the critical desorption pressure, CBM begins to desorb and enter a two-phase gas water flow stage (Gerami et al., 2008). Therefore, in numerical simulations of CBM reservoir development, the water saturation degree is usually set to 100% (Tim, 2012). In recent years, with the exploration and development of deep and low-rank CBM, the free gas in coal reservoirs has gradually been discovered and recognized by the engineering community and academia (Chen et al., 2021). Water saturation in deep and low-rank coal reservoirs is often considered to be <1 (Sun et al., 2018), which challenges the previous assumption that coal reservoirs are always saturated with water (Guo and Kantzas, 2009).

Current research on the impact of water saturation in coal reservoirs on CBM development mainly involves its impact on coal matrix adsorption capacity and reservoir permeability, and they are all based on microscopic experimental studies (Xu et al., 2015a). Experimental studies are necessary for the study of microscopic mechanisms (Li et al., 2023a, 2023b). Regarding adsorption capacity, scholars generally believe that water inhibits CBM adsorption (Gao et al., 2018), and Joubert et al. (1973) and Joubert et al. (1974) suggested that the inhibitory effect disappears when the water saturation is higher than the equilibrium water saturation. Regarding reservoir permeability, Yin et al. (2011) pointed out that the lower the water saturation, the higher the effective methane permeability through experimental research. Liu et al. (2020a, 2020b) studied the evolution of CBM seepage velocity with effective stress in samples with different water saturation. These experimental-scale studies reveal the role of water saturation in controlling the single-factor physical mechanism during CBM production (Meng et al., 2018). This provides a theoretical basis for field studies. However, at the field scale, water and gas production from CBM wells is influenced by a combination of factors. With the development of CBM, the negative effect of effective stress (Zhang et al., 2023a) and the positive effect of matrix shrinkage (Xu et al., 2014) jointly control the absolute permeability of the reservoir, and the change of absolute permeability causes the change of effective permeability (Liu et al., 2023; Zhang et al., 2023b). Meanwhile, during the development process, the reservoir pressure keeps decreasing and the reservoir pressure is different at different locations from the wellbore (Zong et al., 2023). This results in different water saturation at different locations, which ultimately affects gas production. Currently, at the field scale, the comprehensive impact and control mechanism of water saturation in coal reservoirs on gas water production from CBM wells have not been fully studied (Liu et al., 2020a, 2020b; Tang et al., 2017; Thararoop et al., 2012; Xu et al., 2012).

This study constructs a multiphysics coupling numerical model for CBM reservoir development using COMSOL software, simulating the evolution of reservoir and gas production rates under various reservoir water saturation conditions. From the perspective of driving force for fluid flow and resistance for CBM migration, this study explores the control mechanism of reservoir water saturation on CBM production and provides corresponding development recommendations for reservoirs with different water saturation.

Multi-field coupling numerical simulation of development

Geological model

In this study, the geological model is constructed based on the coal reservoirs in the X-1 Block. The geological parameters can be divided into four categories (Table 1): (a) Based on on-site reports (such as well logging reports and well test reports), development design parameters were evaluated, including well radius, well control radius, and original reservoir pressure. (b) Basic coal properties, such as true density, initial permeability, and porosity, were obtained through experimental data of X-1 coal samples. (c) The physical properties of some basic components, such as water, air, and methane, including density and viscosity under standard conditions. (d) The mechanical properties of coal, including elastic modulus, Poisson's ratio, and adsorption strain-related parameters, were obtained from the literature.

Table 1.

Classification and origin of simulation parameters.

Parameter (unit)	Value	Source	Parameter (unit)	Value	Source
Well control radius r_a (m)	2000	Well-logging report, well-test report, development deployment map	True density $ρ_{s}$ (kg·m⁻³)	1417	Lab report
Well radius $r_{well}$ (m)	0.1		Initial permeability $k_{0}$ ( $μ m^{2}$ )	5 × 10⁻⁴
Original reservoir pressure p₀ (psi)	3000		Porosity $\emptyset_{0}$	0.02
Air density ρ_a (kg·m⁻³)	1.22	Physical properties of the basic ingredients (Liu et al., 2020a, 2020b)	Langmuir pressure $P_{L}$ (MPa)	2.22
Relative gravity of methane $γ_{g}$	0.6		Langmuir volume $V_{L}$ (m³·t⁻¹)	31.05
Water density $ρ_{wst}$ (kg·m⁻³)	1000		Adsorption strain-related parameters $ε_{max}$	0.0044	Meng et al. (2021)
Water volumetric coefficient B_w	1.0045		Pore compressibility coefficient $C_{pc}$ (MPa⁻¹)	0.04	Li et al. (2013)
Water compressibility factor C_w (MPa⁻¹)	1.9145 × 10⁻⁴		Elastic modulus E (MPa)	3400	Levine (1996)
The reference pressure for calculating the dynamic density of water p_rw (MPa)	6.8948		Biot coefficient $α$	0.8	Connell et al. (2010)
Water viscosity $μ_{w}$ (mPa·s)	0.73		Poisson's ratio $v$	0.22	Levine (1996)

Geometric model

This study is based on the multiphysics coupling software COMSOL and uses the finite-element method for numerical simulation. CBM development is simplified into a two-dimensional model, with a production well radius of 0.1 m and a well control radius of 2000 m. The reservoir thickness can be set in the two-dimensional model (referring to the thickness of the No. 8 coal seam in X-1 Block, the model is set to 5 m). In this case, the gas and water flow rates of the coal reservoir at different depths are the same, and the thickness is only used to calculate the final production data. The reservoir pressure in the X-1 Block is about 20 MPa, the initial formation pressure is 3000 psi, and the simulated burial depth of the coal seam is 2100 m. Considering symmetry, only one-quarter of the control field was simulated (Figure 1). Due to the strong nonlinearity of the mathematical model of CBM development and difficulty in convergence, an extremely fine linear triangular unit grid (a total of 13,339) was used.

Figure 1.

Numerical simulation of a geometric model of well control field (a quarter).

Mathematical model

This study couples the Darcy flow model, the porous media transfer model, and the geological mechanics model (Figure 2) to construct a mathematical model for CBM development. The detailed modeling process can be found in the author's previously published article (Zhao et al., 2022a). The specific mathematical models utilized in this study are improved in the geological mechanics and porous media transfer modules.

Figure 2.

Coupling diagram of numerical simulation of coalbed methane (CBM) development.

In the flow model, the material conservation equation is equation (1) when single-phase fluid flows in porous media and there is no source or sink (Zhao et al., 2022a).

\frac{\partial}{\partial t} (\emptyset ρ) + \nabla \cdot (ρ u) = 0

(1)

where

ρ

represents the fluid density (kg/m³) and u represents fluid velocity (m/s).

In the porous media phase transfer model, when there are two or more phases of fluid in a porous medium, there is a mass conservation relationship for each phase (Zhao et al., 2022a).

\partial \frac{\emptyset ρ_{i} S_{i}}{\partial t} + \nabla \cdot (ρ_{i} u_{i}) = q S_{i}

(2)

where S_i represents the fluid saturation, dimensionless, i represents the phase of the fluids, and q represents the source term (kg/s/m³).

The diffusion of methane is assumed to be instantaneous. Therefore, the source term can be considered as the derivative of the mass of desorbed gas with respect to time (Zhao et al., 2022b):

q S_{g} = \frac{\partial (M_{ad})}{\partial t}

(3)

where

M_{ad}

represents the adsorbed gas mass in the reservoir per unit volume (kg/m³).

The isothermal adsorption process of methane can be described by the Langmuir-like equation (4) (Zhao et al., 2022b).

M_{ad} = (1 - \emptyset_{0}) ρ_{s} γ_{g} V_{L} ρ_{a} (\frac{p}{P_{L} + p})

(4)

where

ρ_{s}

represents the true density of coal (kg/m³),

γ_{g}

represents the gravity of the gas,

V_{L}

represents the maximum adsorption volume (m³/kg),

ρ_{a}

represents the air density (kg/m³), and

P_{L}

represents the Langmuir pressure (MPa).

In the geological mechanics model, most of the absolute permeability evolution models are derived based on the porous media compression equation without considering adsorption effects, which does not adequately describe the characterization of matrix shrinkage effects (Zhao et al., 2015). This study uses an absolute permeability evolution model that considers the dual effects of matrix shrinkage in coal reservoirs to characterize the impact of geological mechanics during the development process (Xu et al., 2015b; Yuster, 1951; Zhang et al., 2020). This model is derived based on the compressibility equation of adsorptive porous media and can more accurately characterize the influence of matrix shrinkage on permeability evolution (Zhao et al., 2022a, 2022b).

\frac{k}{k_{0}} = \exp {- 3 C_{pc} [\frac{- v - 1}{3 (1 - v)} (p - p_{0}) + (\frac{2 E}{9 (1 - v)} + \frac{E}{3 α (1 - 2 v)}) (\frac{ε_{max} p}{p + p_{L}} - \frac{ε_{max} p_{0}}{p_{0} + p_{L}})]}

(5)

where v represents Poisson's ratio, E represents Young's modulus,

α

represents Biot's coefficient, and

ε_{max}

represents the maximum adsorption strain of coal rock.

In the phase transfer model, the gas water infiltration law cannot be effectively characterized by capillary tube theory in low permeability coal reservoirs. Therefore, this study adopts a plate relative permeability model based on the gas–water interface mechanism to describe the flow law of gas–water two-phase flow (equations (6) to (9)) (Zhao et al., 2021).

S_{w}^{'} = \frac{S_{w} - S_{iw}}{1 - S_{iw}}

(6)

S_{g}^{'} = 1 - S_{w}^{'}

(7)

K_{rw} = {(\frac{3 S_{w}^{' 2} [S_{w}^{' 2} (1 - μ) - 1] - 2 S_{w}^{' 3} [S_{w}^{'} (1 - μ) - 1]}{S_{w}^{'} (1 - μ) - 1})}^{θ}

(8)

K_{rg} {(S_{g}^{' 2} [\frac{S_{g}^{' 2} (1 - a) + 2 S_{g}^{'} - 3}{S_{g}^{'} (1 - a) - 1}])}^{θ}

(9)

where

S_{w}^{'}

represents the effective water saturation,

S_{w}

represents the water saturation,

S_{iw}

represents the bound water saturation,

S_{g}^{'}

represents the effective gas saturation,

K_{rw}

represents the water phase relative permeability,

μ

represents the viscosity ratio of gas phase and water phase, a represents the ratio of water phase viscosity and gas phase viscosity, and

θ

represents the strength of interfacial action, which is related to the reservoir pore structure.

When the water saturation $S_{w}$ is less than the bound water saturation $S_{iw}$ :

K_{rw} = 0, K_{rg} = 1

(10)

The bound water saturation and interfacial action strength of the reservoir are controlled by absolute permeability (equations (11) and (12)) (Zhao et al., 2021):

S_{iw} = - 0.104 \log (k) + 0.29

(11)

θ = - 0.19 \log (k) + 0.9762

(12)

Control effect of water saturation for reservoir exploitation

The pore and fracture networks within coal reservoirs can either contain water or be completely dry. However, when hydraulic fracturing is used for well completion in the development of CBM reservoirs, a large amount of water-based fracturing fluid enters the reservoir (Pan et al., 2022; Tang et al., 2016; Yan et al., 2020; Zheng et al., 2012). Although the water will go through back-draining after the completion of fracturing, there will still be bound water that can’t be discharged (Liu et al., 2014). Therefore, in CBM reservoirs being developed, the range of water saturation can vary from bound water saturation to 1. This paper sets the initial water saturation values to be 1.0 (the gas saturation is actually set to 0.001 so that the gas permeability is not 0), 0.8, 0.6 for free water, and 0.37 for bound water, representing different types of saturated reservoirs. According to the experimental study of coal samples, the bound water saturation in this area is 0.37. Additionally, the well control radius is set to 2000 m, simulating the production of gas and water from an unlimited reservoir.

Simulation results show that different types of saturated reservoirs exhibit significant differences in gas and water production rates (Figure 3). High water-saturation reservoirs are characterized by low gas production and high water production rates, while low water-saturation reservoirs are characterized by high gas production and low water production rates (Ihsan et al., 2013). After 1000 days of development, the gas production rate from a reservoir with an initial water saturation of 1.0 is only 85 m³/day, whereas that from a bound water saturation reservoir can reach up to 5625 m³/day. For water production rates, the reservoir with an initial water saturation of 1.0 has a production rate of 2.8 m³/day, while the reservoir with an initial water saturation of 0.6 has a production rate of <0.1 m³/day. In the X-1 Block, the production characteristics of deep CBM wells are typically characterized by high gas production and low water production rates, with gas production rates of 5000 m³/day and water production rates of 0.1 m³/day. Therefore, numerical simulation research shows that the X-1 Block conforms to the dynamic characteristics of low water saturation CBM development.

Figure 3.

Production of gas and water from coalbed methane (CBM) reservoirs with different water saturation during development.

This paper provides a detailed comparative analysis of pressure distribution, pressure gradient, water saturation, and gas-effective permeability evolution of CBM reservoirs with different water saturation during development. In terms of pressure distribution (Figure 4), numerical simulation results show that water saturation does not have a significant influence, with the bottom-hole pressure of a reservoir with bound water saturation being only 19% lower than that of a water-saturated reservoir after 1000 days of production (Table 2). A pressure gradient is also only slightly affected by water saturation, only 3% (Figure 5 and Table 2).

Figure 4.

Evolution of reservoir pressure for coalbed methane (CBM) reservoirs with different water saturation during development.

Figure 5.

Evolution of pressure gradient for coalbed methane (CBM) reservoirs with different water saturation during development.

Table 2.

Bottom-hole production performance of reservoirs with different water saturation after 1000 days of development.

Water saturation	Pressure (MPa)	Pressure gradient (MPa/m)	Water saturation	Gas effective permeability (mD)	Pressure gradient × gas effective permeability (N/m)	Gas production (m²)
1.0	4.96	32	0.83	1.71 × 10⁻³	5.46 × 10⁻¹¹	85
0.8	5.00	33	0.70	7.94 × 10⁻³	2.65 × 10⁻¹⁰	417
0.6	4.80	27	0.48	4.30 × 10⁻²	1.15 × 10⁻⁹	1860
Bound water saturation	4.00	31	0.37	8.30 × 10⁻²	2.58 × 10⁻⁹	5625
Total change	−19%	−3%	−55%	4764%	4624%	6539%

The water content of CBM reservoirs plays an important role in the evolution of reservoir water saturation. The higher the water saturation of the reservoir, the more water needs to be produced, and the higher the lowest achievable saturation limit (Figure 6) (Xu et al., 2015b; Zhang et al., 2020). After 1000 days of production, the bottom-hole water saturation of different water-saturated CBM reservoirs shows obvious step-like characteristics. The initial water-saturated reservoir saturation decreases to 0.83, but it is still much larger than the bound water saturation, with a 55% difference in proportion (Table 2).

Figure 6.

Evolution of water saturation in coalbed methane (CBM) reservoirs with different water saturation.

The difference in water saturation evolution has an important influence on the gas effective permeability of the gas phase. Numerical simulation results demonstrate that the evolution of gas-effective permeability in CBM reservoirs with different water saturation shows more obvious step-like characteristics (Figure 7). As the space occupied by reservoir water decreases, the number of channels available for gas phase flow increases. Additionally, the smaller the water content is, the smaller the dragging effect of water on the gas phase is. These two factors magnify the difference in water saturation evolution in the evolution of gas effective permeability by several orders of magnitude. After 1000 days of production, the gas effective permeability in the bottom reservoir of CBM reservoirs with bound water reached 4764% of the reservoir with water saturation of 1.

Figure 7.

Evolution of gas effective permeability in coalbed methane (CBM) reservoirs with different water saturation.

According to Darcy's law, the product of gas effective permeability in the bottom reservoir and pressure gradient is directly proportional to the gas production rate of the production well. The difference in initial water saturation of the reservoir affects the effective permeability of the gas phase in the bottom-hole reservoir up to 48 times, whereas it affects the pressure gradient by only 3%. Therefore, the impact of reservoir water saturation on the gas–water production in CBM reservoirs is mainly reflected in the evolution of water saturation and gas-phase effective permeability.

Adaptability analysis of well density to various reservoir water saturation

Gas production under different well densities

In this study, a well control radius of 100 m was set, and the gas production characteristics of CBM reservoirs with different water saturation were compared with those of a well control radius of 2000 m to investigate the impact of well network density on the development of different water saturation reservoirs. Numerical simulation results show that for high water saturation reservoirs (reservoirs with water saturation of 1 and 0.8), the gas production increases with the development after increasing the well network density, and the well interference caused by well network densification has a good production increasing effect. However, for low water saturation reservoirs (reservoirs with water saturation of 0.6 and under bound water conditions), increasing the well network density does not increase the maximum gas production, but rather, it decreases rapidly after reaching the production peak, and the increase in well network densification greatly reduces the single well production. Therefore, the well interference effect caused by high well network density is only applicable to the development of high-water saturation CBM reservoirs, not to low-water saturation CBM reservoirs (Figure 8).

Figure 8.

The effect of well control field on gas production in coal reservoirs with different water saturation.

Control mechanism of well density on gas production in reservoirs with various water saturation

The evolution laws of pressure distribution, pressure gradient, water saturation, and gas-phase effective permeability in different water-bearing reservoirs are analyzed in detail after increasing the well-network density, and compared with the simulation results of the low-well-network density. The control mechanism of wellbore density on the production of gas-bearing reservoirs with different water saturation is explored from two perspectives: the migration force of CBM (reservoir pressure gradient) and the permeability of CBM in the reservoir (gas effective permeability).

Numerical simulations show that under high wellbore density conditions, the reservoir pressure drops significantly after 1000 days of development, which is significantly different from the pressure evolution under low wellbore density. Under high wellbore density conditions, at 100 m from the wellbore, the reservoir pressure does not drop significantly due to continuous replenishment by distant formation water. The maximum decrease is about 3 MPa, and there is no significant difference in the variations of reservoirs with different water saturation.

Under high wellbore density conditions, as the initial water saturation of the reservoir increases, the pressure drop decreases. When the reservoir is initially in the state of bound water, the pressure drop can reach 13 MPa, while when the reservoir is saturated with water, the pressure drop of the reservoir decreases, and the drop is 5 MPa (Figure 9). This is because the water saturation of the reservoir to some extent reflects the size of the elastic potential energy of the fluid in the reservoir, and the lower the water saturation of the reservoir, the lower the energy of the fluid in the reservoir, and the attenuation of the transportation energy is more obvious. In terms of the evolution of the bottom-hole pressure gradient, after 1000 days of development, the pressure gradient of bound water reservoir decays by 68%, while the bottom-hole pressure gradient of saturated water reservoirs decays by only 21% (Figure 10 and Table 3).

Figure 9.

Evolution of pressure distribution in reservoirs with different water saturation.

Figure 10.

Evolution of pressure gradient in reservoirs with different water saturation.

Table 3.

Evolution of bottom-hole pressure gradient in reservoirs with different water saturation.

Time	1	0.8	0.6	Bound water saturation
200 days	37	39	37	28
400 days	35	36	29	18
600 days	33	32	23	13
800 days	32	29	19	11
1000 days	29	26	15	9
Total change	−21%	−33%	−58%	−68%

Under the condition of high well network density, as the initial water saturation of the reservoir increases, the decrease of reservoir water saturation significantly increases (Figure 11). When the reservoir is initially saturated with water, the reservoir water saturation decreases from 1 to 0.65 during the process of development, with a decrease of 0.35. Conversely, when the initial water saturation of the reservoir is low, especially when it is under the condition of bound water, the saturation does not change at all. The difference in saturation decrease rate causes variations in the increase of gas effective permeability in reservoirs with different water saturation (Figure 12). A high decrease rate in saturation leads to a great increase in gas-effective permeability. When the reservoir is initially saturated with water, the bottom-hole gas effective permeability increases by 526% during the development process, while in a reservoir of bound water condition, it only undergoes a slight increase of 2%, which is due to the evolution of absolute permeability (Table 4).

Figure 11.

Evolution of water saturation in reservoirs with different water saturation.

Figure 12.

Evolution of gas effective permeability in reservoirs with different water saturation.

Table 4.

Evolution of bottom-hole gas effective permeability in reservoirs with different water saturation.

Time	1.0	0.8	0.6	Bound water saturation
200 days	0.02	0.07	0.28	0.82
400 days	0.03	0.10	0.33	0.83
600 days	0.06	0.14	0.38	0.83
800 days	0.09	0.17	0.44	0.83
1000 days	0.12	0.20	0.51	0.83
Total change	526%	172%	82%	2%

In summary, the high well density results in a wider range of water saturation variation and a significant increase in gas phase effective permeability in high water saturation reservoirs. The gas phase's effective permeability can increase up to 526% when the reservoir is initially saturated with water. Moreover, due to the higher fluid elastic potential in high water saturation reservoirs, the decay of flow kinetic energy is slower, with only a 21% decrease in saturation level when the reservoir is initially saturated with water. The slow decay of kinetic energy and rapid increase in gas phase effective permeability result in the gas production of high water-bearing reservoirs increasing during development, so high well network density is required for development. In low water saturation reservoirs with high well density, the increase in gas phase effective permeability is limited, with a minimum increase of 2%. The fluid kinetic energy decay is significant, with a maximum decay of 68%. The rapid decay of kinetic energy and slow increase in gas phase effective permeability cause the gas production rate in low water saturation reservoirs to peak and then continuously decline. Interwell interference often has a negative effect, and therefore, low well density should be used to increase single well control of reserves.

Conclusions

Based on an absolute permeability evolution model that considers the dual impact of matrix shrinkage and a relative permeability model that considers gas–water interface interaction, a multiphysics coupling numerical model for CBM development that fully considers the permeability evolution characteristics of adsorptive low permeability coal reservoirs was constructed.

High-water saturation reservoirs show the characteristics of low gas production and high-water production, while low-water saturation reservoirs show high gas production and low water production. The simulation results show that the decrease in reservoir water saturation greatly enhances the gas phase effective permeability of the reservoir, which can be increased by up to 4764%, significantly increasing the gas production rate up to 6539%. This is the main mechanism by which reservoir water saturation controls production.

For high water saturation reservoirs, high well density should be used for development to make use of the positive effects of well interference, which can quickly increase the effective gas phase permeability by 526%. For low water saturation reservoirs, low well density should be used to increase the control reserves of individual wells and avoid the negative effects of well interference, such as a −68% drop in the bottom-hole pressure gradient.

Footnotes

Declaration of conflicting interests

The authors declared no potential conflicts of interest with respect to the research, authorship, and/or publication of this article.

Funding

The authors disclosed receipt of the following financial support for the research, authorship, and/or publication of this article: This work was supported by the National Natural Science Foundation Project, China (grant number 42172188).

ORCID iD

Hao Xu

References

Akhondzadeh

Keshavarz

Sayyafzadeh

, et al. (2018) Investigating the relative impact of key reservoir parameters on performance of coalbed methane reservoirs by an efficient statistical approach. Journal of Natural Gas Science and Engineering 53: 416–428.

Bear

(1975) Dynamics of fluids in porous media. Soil Science 120: 162–163.

Chen

Tang

Tao

, et al. (2018) Characteristics of in-situ stress distribution and its significance on the coalbed methane (CBM) development in Fanzhuang-Zhengzhuang Block, Southern Qinshui Basin, China. Journal of Petroleum Science and Engineering 161: 108–120.

Chen

Tao

Tian

, et al. (2021) Hydrogeological control on the accumulation and production of coalbed methane in the Anze block, southern Qinshui Basin, China. Journal of Petroleum Science and Engineering 198: 108138.

Clarkson

(2009) Case study: Production data and pressure transient analysis of horseshoe canyon CBM wells. Journal of Canadian Petroleum Technology 48(10): 27–38.

Connell

Pan

, et al. (2010) Coal permeability and its behaviour with gas desorption, pressure and stress. Paper presented at the SPE Asia Pacific Oil and Gas Conference and Exhibition, Brisbane, Queensland, Australia, 18–20 October. SPE-133915-MS.

Gao

Xie

Pan

, et al. (2018) Analysis of gas content and development geological model of deep coalbed methane in Linxing. Journal of China Coal Society 43: 1634–1640.

Gerami

Pooladi-Darvish

Morad

, et al. (2008) Type curves for dry CBM reservoirs with equilibrium desorption. Journal of Canadian Petroleum Technology 47(7): PETSOC-08-07-48.

Guo

Kantzas

(2009) Assessing the water uptake of Alberta coal and the impact of CO injection with low-field NMR. Journal of Canadian Petroleum Technology 48(7): 40–46.

10.

Ihsan

Talal

Sara

(2013) Coal seam gas and associated water: A review paper. Renewable and Sustainable Energy Reviews 22: 550–560.

11.

Joubert

Grein

Bienstock

(1973) Sorption of methane in moist coal. Fuel 52: 181–185.

12.

Joubert

Grein

Bienstock

(1974) Effect of moisture on the methane capacity of American coals. Fuel 53: 186–191.

13.

Levine

(1996) Model study of the influence of matrix shrinkage on absolute permeability of coal bed reservoirs. Geological Society, London, Special Publications 109: 197–212.

14.

, et al. (2023a) Experimental study on compressive behavior and failure characteristics of imitation steel fiber concrete under uniaxial load imitation steel fiber concrete under uniaxial load. Construction and Building Materials 399(8): 132599.

15.

, et al. (2023b) Research on acoustic emission multi-parameter characteristics in the failure process of imitation steel fiber reinforced concrete. Physics of Fluids 35(6): 107109.

16.

Tang

Pan

, et al. (2013) Characterization of the stress sensitivity of pores for different rank coals by nuclear magnetic resonance. Fuel 111: 746–754.

17.

Liu

, et al. (2023) Influence of hole diameter on mechanical properties and stability of granite rock surrounding tunnels. Physics of Fluids 35(6): 064121.

18.

Liu

Tang

, et al. (2020a) Analysis of hydraulic fracture behavior and well pattern optimization in anisotropic coal reservoirs. Energy Exploration and Exploitation 39(1): 299–317.

19.

Liu

Tang

, et al. (2020b) Coalbed methane production of a heterogeneous reservoir in the Ordos Basin, China. Journal of Natural Gas Science and Engineering 82: 103502.

20.

Liu

Zhang

, et al. (2014) Effect of moisture content on seepage of coalbed methane. Journal of China Coal Society 39: 1840–1844.

21.

Meng

Tang

, et al. (2021) Laboratory investigation on methane sorption-induced strain and permeability in middle and high rank coal samples. Journal of China Coal Society 46: 1915–1924.

22.

Meng

Wang

, et al. (2018) An improved productivity model in coal reservoir and its application during coalbed methane production. Journal of Natural Gas Science and Engineering 49: 342–351.

23.

Pan

Mou

, et al. (2022) Micro-nano-scale pore stimulation of coalbed methane reservoirs caused by hydraulic fracturing experiments. Journal of Petroleum Science and Engineering 214: 110512.

24.

Sun

Tian

Chen

, et al. (2018) Diversified accumulation characteristics and exploration directions of low-grade coalbed methane in China. Natural Gas Industry 38: 10–18.

25.

Tang

Zhao

, et al. (2016) Dynamic balance mechanism of material energy in medium–high coal rank coalbed methane system. Journal of China Coal Society 41: 40–48.

26.

Tang

, et al. (2017) Controlling factors of coalbed methane well productivity of multiple superposed coalbed methane systems: A case study on the Songhe mine field, Guizhou, China. Energy Exploration and Exploitation 35(6): 665–684.

27.

Thararoop

Karpyn

Ertekin

(2012) Development of a multi-mechanistic, dual-porosity, dual-permeability, numerical flow model for coalbed methane reservoirs. Journal of Natural Gas Science and Engineering 8: 121–131.

28.

Tim

(2012) Coalbed methane: A review. International Journal of Coal Geology 101: 36–81.

29.

Tang

Liu

, et al. (2012) Study on coalbed methane accumulation characteristics and favorable areas in the Binchang area, southwestern Ordos Basin, China. International Journal of Coal Geology 95: 1–11.

30.

Tang

, et al. (2014) A dynamic prediction model for gas–water effective permeability based on coalbed methane production data. International Journal of Coal Geology 121: 44–52.

31.

Tang

Zhao

, et al. (2015b) Geologic controls of the production of coalbed methane in the Hancheng area, southeastern Ordos Basin. Journal of Natural Gas Science and Engineering 26: 156–162.

32.

Tang

Zhao

, et al. (2015a) A new laboratory method for accurate measurement of the methane diffusion coefficient and its influencing factors in the coal matrix. Fuel 158: 239–247.

33.

Yan

Tang

Zhang

, et al. (2020) Analysis of productivity differences in vertical coalbed methane wells in the Shizhuangnan Block, Southern Qinshui Basin, and their influencing factors. Energy Exploration and Exploitation 38(5): 1428–1453.

34.

Yin

Jiang

, et al. (2011) Experimental study of influences for water content in coalbed gas reservoirs on methane seepage. Chinese Journal of Rock Mechanics and Engineering 30: 3401–3406.

35.

Yuster

(1951) Theoretical considerations of multiphase flow in idealized capillary systems. Proceedings of the Third World Petroleum Congress 2: 437–445.

36.

Zhang

Feng

Zhang

, et al. (2020) Multi-fractured horizontal well for improved coalbed methane production in eastern Ordos basin, China: Field observations and numerical simulations. Journal of Petroleum Science and Engineering 194: 107488.

37.

Zhang

Wang

, et al. (2023a) Research on macroscopic mechanical properties and microscopic evolution characteristic of sandstone in thermal environment. Construction and Building Materials 13(21): 12013.

38.

Zhang

Wang

Zhang

, et al. (2023b) Research on nonlinear variation of elastic wave velocity dispersion characteristic in limestone dynamic fracture process. Fractals 31(1): 2350008.

39.

Zhao

Tang

, et al. (2015) High production indexes and the key factors in coalbed methane production: A case in the Hancheng block, southeastern Ordos Basin, China. Journal of Petroleum Science and Engineering 130: 55–67.

40.

Zhao

Tang

, et al. (2021) A new relative permeability model of coal reservoir considering interface effect. Journal of Petroleum Science and Engineering 207: 109082.

41.

Zhao

Tang

, et al. (2022a) Dual mechanisms of matrix shrinkage affecting permeability evolution and gas production in coal reservoirs: Theoretical analysis and numerical simulation. Journal of Natural Gas Science and Engineering 108: 104844.

42.

Zhao

Tang

, et al. (2022b) The theoretical basis of model building for coal reservoir permeability: A review and improvement. Journal of Natural Gas Science and Engineering 106: 104744.

43.

Zheng

Pan

Chen

, et al. (2012) Laboratory study of gas permeability and cleat compressibility for CBM/ECBM in Chinese coals. Energy Exploration and Exploitation 30(3): 451–476.

44.

Zong

Tang

, et al. (2023) A dynamic prediction model of reservoir pressure considering stress sensitivity and variable production. Geoenergy Science and Engineering 225: 211688.

Influence of reservoir water saturation on coalbed methane development: Based on multi-field coupling numerical simulation

Abstract

Keywords

Introduction

Multi-field coupling numerical simulation of development

Geological model

Geometric model

Mathematical model

Control effect of water saturation for reservoir exploitation

Adaptability analysis of well density to various reservoir water saturation

Gas production under different well densities

Control mechanism of well density on gas production in reservoirs with various water saturation

Conclusions

Footnotes

Declaration of conflicting interests

Funding

ORCID iD

References