Mathematical model of methane driven by hydraulic fracturing in gassy coal seams

Abstract

During hydraulic fracturing in gassy coal, methane is driven by hydraulic fracturing. However, its mathematical model has not been established yet. Based on the theory of ‘dual-porosity and dual-permeability’ fluid seepage, a mathematical model is established, with the cleat structure, main hydraulic fracture and methane driven by hydraulic fracturing considered simultaneously. With the help of the COMSOL Multiphysics software, the numerical solution of the mathematical model is obtained. In addition, the space–time rules of water and methane saturation, pore pressure and its gradient are obtained. It is concluded that (1) along the direction of the methane driven by hydraulic fracturing, the pore pressure at the cleat demonstrates a trend of first decreasing and later increasing. The pore pressure gradient exhibits certain regional characteristics along the direction of the methane driven by hydraulic fracturing. (2) Along the direction of the methane driven by hydraulic fracturing, the water saturation exhibits a decreasing trend; however, near the cleat or hydraulic fracture, the water saturation first increases and later decreases. The water saturation in the central region of the coal matrix block is smaller than that of its surrounding region, while the saturation of water in the entire matrix block is greater than that in the cleat or hydraulic fracture surrounding the matrix block. The water saturation at the same space point increases gradually with the time progression. The space–time distribution rules of methane saturation are contrary to those of the water saturation. (3) The free methane driven by hydraulic fracturing includes the original free methane and the free methane desorbed from the adsorption methane. The reduction rate of the adsorption methane is larger than that of free methane.

Keywords

Effect of methane driven by hydraulic fracturing water saturation methane saturation pore pressure pore pressure gradient

Introduction

During the process of hydraulic fracturing in gassy coal, the methane concentration of the airflow in the roadway increases considerably (Huang and Lu, 2018; Lu et al., 2019). This phenomenon is known as the effect of methane driven by hydraulic fracturing. The existence of this phenomenon has been demonstrated in laboratory experiments, and the effect has both advantages and disadvantages in production practice (Chen, 2015; Huang et al., 2013, 2016). Although the effect can provide novel approaches and technologies for the exploitation of coalbed gas and hydraulic measures used to eliminate methane outbursts, it can increase the methane concentration of the airflow in the roadway in addition to the methane content and methane pressure of the local region in the coal seam, which is detrimental to the prevention of coal and methane outbursts and methane explosions (Cheng, 2012; Huang et al., 2016). However, the mathematical model of methane driven by hydraulic fracturing, considering the effect of methane adsorption and desorption, has not been established yet. Furthermore, the tendencies of water and methane saturation, pore pressure and its gradient are yet unknown. Owing to the abovementioned reasons, it is desirable to establish a mathematical model of methane driven by hydraulic fracturing.

A coal seam is characterized by the complex and variable structure of the pore and fracture, which directly determines the transport mode of the fluid (methane, water, etc.) in the coal seam. At present, the two-phase seepage model of water and gas in the coal seam with double porosity is divided into the double porosity and single seepage model and the double porosity and dual seepage model. Darcy’s law is used to describe the process of fracturing fluid flow in the abovementioned two models (Hosking et al., 2018; Thararoop et al., 2012; Yun et al., 2018; Zhang, 2014; Zhang and Bian, 2015; Zhang et al., 2016; Zou et al., 2015). Only the fracture permeability is considered in the double porosity and single permeability model, while both the fracture and pore permeability are considered in the double porosity and dual permeability model. In other words, only one permeability exists in the double porosity and single permeability model; however, two different permeabilities exist in the double porosity and double permeability model (King et al., 1986; Remner et al., 1984a, 1984b). Most natural fractures in the coal seam are in a closed state, which has negligible effect on the permeability improvement of the coal seam. The connectivity of the pores in a coal seam, which is increased by the hydraulic fractures, is significantly different from that pertaining to the natural fractures (Chong et al., 2017; Jamison and Azad, 2017; Preisig et al., 2015). At present, the permeability of these two types of fractures has not been distinguished clearly in the mathematical models of water–methane two-phase seepage flow in dual–medium coal seams. A study on the mathematical model and rules, in which the permeability difference of the above two types of fractures is considered for the methane driven by hydraulic fracturing, has not been conducted yet.

Based on the theory of ‘dual-porosity and dual-permeability’ fluid seepage in a dual-porosity coal seam, the conceptual and mathematical model of methane driven by hydraulic fracturing is established, with the cleat structure of the coal seam, main hydraulic fracture and methane-driven process by hydraulic fracturing considered simultaneously. With the help of the COMSOL Multiphysics software, the numerical solution of the mathematical model is obtained. In addition, the space–time rules of water and methane saturation, pore pressure and its gradient are obtained.

Conceptual model of the methane driven by hydraulic fracturing

Pores and fractures exist in a coal seam, which make the coal reservoir capable of gas storage and allow the coalbed methane to participate in the diffusion-seepage-migration process. The pores are the main storage place of coal seam methane, and their structural characteristics have a significant effect on the accumulation and diffusion migration of the coal seam methane. Fractures are the main channels of methane migration; they are related to the permeability of the reservoir and determine the yield of the coal seam methane. The natural fractures in the coal seam are also known as cleats (Figure 1). The continuous cleats in the coal seam are called face cleats, and the discontinuous cleats, which end at the surface cleat or intersect with a face cleat, are called butt cleats. These two cleats are usually perpendicular or nearly orthogonal to each other.

Figure 1.

Diagram of the coal seam texture.

Both free and adsorbed methane are present in the coal seam, and they occur in an equilibrium state. With the change in the coal seam pressure, the adsorbed and free methane move and perform exchange between each other constantly. To fully utilize the methane released from the coal seam, it is necessary to ensure that the adsorbed methane is constantly changed into free methane. In other words, the essence of the methane output is a process of continuous release of free methane. Under the action of high-pressure water, hydraulic fractures are generated in the coal seam. In the context of natural fractures, hydraulic fractures can be considered artificial fractures. The main hydraulic fractures usually extend along the direction of the maximum principal stress. When high-pressure water penetrates the natural fractures, artificial fractures and matrix pores, it drives the free methane in the natural fractures, artificial fractures and matrix pores into free space. Based on the structural characteristics of the coal seam, the shape of the main hydraulic fracture and the process of methane driven by hydraulic fracturing, the conceptual model of the methane driven by hydraulic fracturing can be established, as shown in Figure 2.

Figure 2.

Conceptual model of methane driven by hydraulic fracturing in a coal seam.

Mathematical model of the methane driven by hydraulic fracturing

Basic assumptions

The basic assumptions for the mathematical model are as follows:

The solubility of the coal seam methane in water is ignored, and thus the process of methane driven by hydraulic fracturing is a water–gas two-phase flow in the coal seam.

The temperature of the seepage field is constant. A state of thermal balance exists among the water, methane and coal rock, and no heat transfer occurs among them. The flow of fluid in the artificial fractures, natural fractures and matrix pores pertains to laminar migration, which obeys Darcy’s law.

The solid skeleton and pores of the coal rock are compressible.

The porous medium of the coal is elastic and its deformation is small. The capillary pressure and gravity are considered.

Coal seam methane is an ideal gas, and the effect of methane adsorption occurs only in the matrix pores.

The coal seam is a double medium with pores and fractures, which are characterized by independent porosity and permeability.

Governing equations

Governing equations of stress

Based on the elastic theory of porous media, the constitutive equation, which considers the effect of methane adsorption and desorption, is as follows (Li et al., 2016; Liu et al., 2016).

\begin{array}{l} ε_{ij} = ︸ \underset{Strain c aused b y geostress}{\frac{1}{2 G} σ_{ij} - (\frac{1}{6 G} - \frac{1}{9 K}) σ_{kk} δ_{ij}} + ︸ \underset{Strain c aused b y pore pressure}{\frac{α_{m}}{3 K} p_{m} δ_{i j}} + ︸ \underset{Strain c aused b y frature pressure}{\frac{α_{f}}{3 K} p_{f} δ_{i j}} + ︸ \underset{Strain c aused b y adsorption swelling}{\frac{ε_{s}}{3} δ_{i j}} \end{array}

(1)

where

{\begin{array}{l} i, j, k = 1, 2, 3 \\ G = E / [2 (1 + ν)] \\ E = 1 / [1 / E_{s} + 1 / (a k_{n})] \\ K = 2 G (1 + ν) / 3 (1 - 2 ν) = E / 3 (1 - 2 ν) \\ α_{m} = 1 - K / K_{s} \\ α_{f} = 1 - K / (a K_{n}) \\ ε_{s} = β_{s g} V_{s g} = β_{s g} \frac{V_{L} p_{m}}{p_{m} + P_{L}} \end{array}

(2)

where G is the shear modulus; E is the elastic modulus;

ν

is Poisson’s ratio; K_s is the volume of the coal skeleton; E_s is the elastic modulus of the coal skeleton; a is the width of the coal matrix block; K_n is the stiffness of the fractures;

σ_{i j}

is the component of stress; p_m is the pore pressure; p_f is the fluid pressure in the fracture;

α_{m} (\leq 1)

is the Biot-effective stress coefficient corresponding to the pores;

α_{f} (\leq 1)

is the Biot-effective stress coefficient corresponding to the fractures; K is the volume modulus of coal;

δ_{i j}

is the Kronecker symbol (when i = j,

δ_{i j} = 1

; when i ≠ j,

δ_{i j} = 0

);

ε_{s}

is the adsorption strain of methane;

β_{s g}

is the adsorption strain coefficient of the coalbed methane;

V_{s g}

is the content of the adsorbed methane; P_L is the Langmuir pressure constant and V_L is the Langmuir volume constant.

According to the theory of elasticity, the relation between strain and displacement is as follows

ε_{i j} = \frac{1}{2} (u_{i, j} + u_{j, i})

(3)

It can be obtained from the static equilibrium condition that

σ_{i j, j} + F_{i} = 0

(4)

Substituting equations (1), (2) and (3) into equation (4), and the modified Navier balance equation (5) with the displacement as the fundamentally unknown quantity, the coupling term included can be obtained.

G u_{i, j}_{j} + (G + λ) u_{j, j i} - α_{m} p_{m}_{, i} - α_{f} p_{f, i} - K ε_{s, i} + F_{i} = 0

(5)

where u_i is the displacement along the i direction; and F_i is the volume force along the i direction

(i = x, y, z)

2. Mass conservation equations

The hydraulic fractures are denoted by subscript ‘f’, and the natural fractures are denoted by subscript ‘f₁’. The matrix micro-fractures or pores are denoted by subscript ‘m’, and the matrix is denoted by subscript ‘M’.

① Mass conservation equations of water and methane in the hydraulic fractures

The mass conservation equation of water in the hydraulic fractures is as follows

\nabla \cdot [\frac{ρ_{g} K_{f} K_{r g}}{μ_{g}} (\nabla p_{f g} - ρ_{g} g \nabla D)] + Q_{m f} + Q_{f}_{f 1} + Q_{g} = \frac{\partial}{\partial t} (ϕ_{f} ρ_{g} S_{f g})

(6)

The mass conservation equation of methane in the hydraulic fractures is as follows

\begin{array}{l} \nabla \cdot [\frac{ρ_{w} K_{f} K_{r w}}{μ_{w}} (\nabla p_{f w} - ρ_{w} g \nabla D)] - \frac{k_{m} k_{r}_{w} ρ_{w}}{μ_{w}} (p_{f w} - p_{m w}) \\ - \frac{k_{f 1 w} k_{r}_{w} ρ_{w}}{μ_{w}} (p_{f w} - p_{f 1 w}) + Q_{w} = \frac{\partial}{\partial t} (ϕ_{f} ρ_{w} S_{f w}) \end{array}

(7)

② Mass conservation equations of water and methane in the natural fractures

The mass conservation equation of water in the natural fracture is as follows

\nabla \cdot [\frac{ρ_{g} K_{f 1} K_{r g}}{μ_{g}} (\nabla p_{f 1 g} - ρ_{g} g \nabla D)] + Q_{M f 1} {+Q}_{m f 1} - Q_{f}_{f 1} = \frac{\partial}{\partial t} (ϕ_{f 1} ρ_{g} S_{f 1}_{g})

(8)

The mass conservation equation of methane in the natural fracture is as follows

\begin{array}{l} \nabla \cdot [\frac{ρ_{w} K_{f 1} K_{r w}}{μ_{w}} (\nabla p_{f 1 w} - ρ_{w} g \nabla D)] + \frac{k_{f 1 w} k_{r}_{w} ρ_{w}}{μ_{w}} (p_{f w} - p_{f 1 w}) \\ - \frac{k_{m w} k_{r}_{w} ρ_{w}}{μ_{w}} (p_{f 1 w} - p_{m w}) = \frac{\partial}{\partial t} (ϕ_{f 1} ρ_{w} S_{f 1 w}) \end{array}

(9)

③ Mass conservation equations of water and methane in the matrix

The mass conservation equations of methane in the matrix are as follows

\begin{array}{l} \nabla \cdot [\frac{ρ_{g} K_{m} K_{r g}}{μ_{g}} (\nabla p_{m g} - ρ_{g} g \nabla D)] + Q_{m M} - Q_{m f 1} \\ = \frac{\partial}{\partial t} [ϕ_{m} ρ_{g} s_{m g} + (1 - ϕ_{f} - ϕ_{f 1} - ϕ_{m}) \frac{V_{L} p_{mg}}{p_{L} + p_{mg}}] \end{array}

(10)

The mass conservation equations of water in the matrix are as follows

\begin{array}{l} \nabla \cdot [\frac{ρ_{w} K_{m} K_{mrw}}{μ_{w}} (\nabla p_{m w} - ρ_{w} g \nabla D)] + \frac{K_{m} K_{r}_{w} ρ_{w}}{μ_{w}} (p_{f w} - p_{m w}) \\ + \frac{K_{m} K_{r}_{w} ρ_{w}}{μ_{w}} (p_{f 1 w} - p_{m w}) = \frac{\partial}{\partial t} (ϕ_{m} ρ_{w} s_{m w}) \end{array}

(11)

④ Exchange quantity of methane between matrix system and fracture system

The methane diffuses into the fracture system from the coal matrix. The formula for calculating the methane exchange quantity between the matrix and fracture system is as follows (Feng et al., 2015; Li et al., 2016)

- Q_{M}_{f} - Q_{M} {_{f}}_{1} - Q_{M}_{m} = δ_{f} \frac{ρ_{g} K_{m} K_{r g}}{μ_{g}} (p_{m g} - p_{f g})

(12)

\frac{K_{m} K_{r}_{w} ρ_{w}}{μ_{w}} (p_{f w} - p_{m w}) + \frac{K_{m} K_{r}_{w} ρ_{w}}{μ_{w}} (p_{f 1 w} - p_{m w}) = δ_{f} \frac{ρ_{w} K_{m} K_{r}_{g}}{μ_{w}} (p_{f w} - p_{M w})

(13)

where Q_g is the exchange quantity of the methane diffusing from the matrix system to the fracture system; δ is the fracture shape factor of the matrix block, and

δ = π^{2} [(\frac{1}{L_{x}^{2}}) + (\frac{1}{L_{y}^{2}}) + (\frac{1}{L_{z}^{2}})]

L_x, L_y and L_z are the intervals of fractures along the x, y and z directions, respectively.

3. Relationship equation of capillary pressure with the pressure and saturation in the fracture and matrix

The expression of capillary pressure in the hydraulic fractures is as follows

p_{f c} = p_{f n}_{w} - p_{f w} = p_{f c} (S e_{f w})

(14)

The expression of capillary pressure in the natural fractures is as follows

p_{f 1 c} = p_{f 1 n}_{w} - p_{f 1}_{w} = p_{f 1 c} (S e_{f 1 w})

(15)

The expression of capillary pressure in the matrix pores is as follows

p_{c m} = p_{m n}_{w} - p_{m w} = p_{m}_{c} (S e_{m w})

(16)

The pores and fractures are always filled with fluid. The water and methane two-phase saturation in the hydraulic fracture is expressed as follows

S_{f g} + S_{f w} = 1

(17)

The expression of water and methane saturation in the natural fracture is as follows

S_{f 1}_{g} + S_{f 1}_{w} = 1

(18)

The water and methane two-phase saturation in the matrix pores can be expressed as follows

S_{mg} + S_{m w} = 1

(19)

4. Porosity and absolute permeability of coal matrix and fracture

The porosity of the coal matrix and fracture can be respectively expressed as follows (Fan et al., 2016; Liu et al., 2017; Zhang et al., 2017).

ϕ_{m} = \frac{1}{1 + S} [ϕ_{m 0} (1 + S_{0}) + α_{m} (S - S_{0})]

(20)

ϕ_{f} = ϕ_{f 0} - \frac{3 ϕ_{f 0} (Δ ε_{s} - Δ ε_{v})}{ϕ_{f 0} + 3 K_{f} / K}

(21)

where

φ_{m}

is the porosity of the coal matrix;

S = ε_{v} + p_{m} / K_{s} - ε_{s}

;

S_{0} = ε_{v 0} + p_{m 0} / K_{s} - ε_{s 0}

; K_s is the volume modulus of the coal skeleton;

α_{m}

is the Biot-effective stress coefficient corresponding to the pore;

φ_{f}

is the porosity of the fracture;

ε_{s}

is the skeleton strain caused by the adsorbed methane;

ε_{v}

is the volume strain of coal, and

ε_{v} = ε_{x} + ε_{y} + ε_{z}

;

ε_{x}

ε_{y}

ε_{z}

are the respective volume strains along the directions of x, y and z; K_f is the improved stiffness of the fracture, and K_f = bK_n; b is the width of the fracture and subscript ‘0’ represents the initial value.

A cubic relationship exists between the absolute permeability of the matrix and fracture and their porosity, and the corresponding expressions are as follows (Li et al., 2016; Palmer, 2009; Zhang, 2016).

\frac{k_{m}}{k_{m 0}} = {(\frac{ϕ_{m}}{ϕ_{m 0}})}^{3} {(\frac{1 - ϕ_{m 0}}{1 - ϕ_{m}})}^{2}

(22)

\frac{k_{f}}{k_{f 0}} = {(\frac{ϕ_{f}}{ϕ_{f 0}})}^{3} = {(1 - \frac{3 (Δ ε_{s} - Δ ε_{v})}{ϕ_{f 0} + 3 K_{f} / K})}^{3}

(23)

5. Relative permeability of water and methane

Based on the Brooks and Corey-Burdine (BCB) model, a new relative permeability model of water and methane, applicable to a coal reservoir with joints, was deduced by Chen Dong et al. The expression of the new model is as follows (Chen et al., 2013; Corey and Rathjens, 1956).

k_{r w} = S e_{w}^{η + 1 + 2 / λ}

(24)

k_{r g} = {(1 - S e_{w})}^{η} [1 - S e_{w}^{1 + 2 / λ}]

(25)

where λ is the index representing the scale distribution of coal joints; and η is the index representing the bend degree of coal joints.

The physical meaning of the parameters λ and η in equations (24) and (25) can be illustrated by Figure 3. Research has demonstrated that the parameters λ and η considerably influence the results of the numerical simulation; the range of parameter η is between 0 and 2, and the range of parameter λ is between 0.3 and 8.8 (Chen et al., 2013; Bertrand et al., 2017).

Figure 3.

Physical characteristics of fracture seepage presented by parameters λ and η.

Boundary conditions

Water enters into the interior of the coal body along the hydraulic fractures or cleats, and it causes the seepage of methane forward in the coal seam. The mathematical form of the boundary and initial conditions of the physical model is as follows

p_{w} = p_{w} (t) \partial Ω Borehole boundary

(26)

p_{g} = p_{g}_{0} \partial Ω Outlet boundary

(27)

n \cdot ρ_{w} [- λ_{w} (\nabla p_{w} + ρ_{w} g \nabla D)] = 0 \partial Ω Bottom and top boundaries

(28)

n \cdot ρ_{g} [- λ_{g} (\nabla p_{g} + ρ_{g} g \nabla D)] = 0 \partial Ω Bottom and top boundaries

(29)

{p = p}_{0} \partial Ω

(30)

Numerical simulation of the methane driven by hydraulic fracturing

Geometric models

The geometric model (Figure 4) of the coal seam used in this study has a length of 20 m and width of 6 m. The blue circle represents the fracturing borehole, whose diameter is 0.075 m. The natural fractures and the main hydraulic fractures are represented by the black and magenta lines, respectively. The roof of the coal seam, which bears the overburden gravity, is simulated by the upper boundary of the geometric model. The floor of the coal seam is simulated by the lower boundary of the geometric model, which is a fixed boundary. Neither water nor methane can penetrate the upper and lower boundaries. Both the right and left boundaries of the geometric model are sliding and permeable boundaries. The maximum principal stress is along the vertical direction, and the minimum principal stress direction is along the horizontal direction.

Figure 4.

Geometric model.

Parameter selection

To make the process of numerical simulation as factually accurate as possible, real physical parameters are used in the numerical simulation. The relevant parameters used in this numerical simulation are listed in Table 1.

Table 1.

Parameters used in the numerical simulation.

Symbol	Physical parameter	Value	Unit	Symbol	Physical parameter	Value	Unit
ρ_w	Density of water	1000	kg/m³	ρ_gsc	Density of methane under standard condition	0.716	kg/m³
E	Elastic modulus of coal	2713	MPa	K_n	Stiffness of fracture	4800	MPa/m
E_s	Elastic modulus of coal skeleton	8469	MPa	φ_m0	Initial porosity of matrix	0.04	1
$ν$	Poisson’s ratio	0.32	1	φ_f0	Initial porosity of fracture	0.08	1
μ_g	Dynamic viscosity of methane	1.84⁻⁵	Pa·s	S_r,w	Saturation of bound water	0.2	1
μ_w	Dynamic viscosity of water	1.01⁻³	Pa·s	S_r,g	Saturation of residual methane	0.1	1
ρ_s	Density of coal skeleton	1470	kg/m³	k_r,w0	Relative permeability of water at endpoint	1	1
k_r,g0	Relative permeability of methane at endpoint	1	1	K_m	Permeability of matrix	1e-18	m²
a	Langmuir limiting methane adsorption constant	0.036	m³/kg	p_w	Water pressure	10	MPa
p₀	Initial pore pressure	1	MPa	σ₃	Minimum principle stress	6	MPa
σ₁	Maximum principle stress	10	MPa	K_f	Permeability of fracture	1e-16	m²
b	Langmuir methane adsorption constant	3.304	MPa⁻¹

Analysis of numerical simulation results

Space–time distribution rules of pore pressure

The space distribution of the pore pressure at the times 0.1 h, 0.5 h, 1.0 h and 1.5 h is shown in Figure 5. It can be seen from Figure 5 that the (1) methane pressure at the boundary of the coal seam is lower than that inside the coal seam. This finding can be explained by the fact that the methane flows into the cleat or hydraulic fracture from the matrix block and migrates around the coal seam under the driving action of the high-pressure water. The free methane in the coal seam flows through the permeable boundary to the free space around the coal seam, resulting in the methane pressure at the boundary of the coal seam being lower than that inside the coal seam. (2) At a central time, the pore pressure in a matrix block is considerably higher than that in the adjacent cleat or hydraulic fracture. This phenomenon is due to the occurrence of the coal seam methane in the matrix block. (3) The methane pressure inside the hydraulic fracture and cleat connected to the hydraulic fracture is significantly higher than that in other parts. This aspect can be explained by the fact that the high-pressure water, which flows into the coal matrix through hydraulic fractures, has an effect of compression on the free methane in the coal matrix, resulting in the methane pressure of the internal coal matrix being higher than that within its surrounding area.

Figure 5.

Space–time distribution of pore pressure.

It can be seen from Figure 6 that (1) when the time is fixed, the pore pressure of the coal seam is continuously distributed along the direction of the methane driven by water. The initial pore pressure of the coal seam is 1.2 MPa, while the range of the pore pressure during the process of methane driven by water is 0.1–2.0 MPa. Among these values, the pore pressure near the borehole is higher than the initial pore pressure, while the pore pressure near the boundary of the coal seam is lower than the initial pore pressure. It is indicated that after the high-pressure water enters the coal seam, the pore pressure in the coal seam is redistributed, resulting in the appearance of zones involving an increase and decrease in the pore pressure. (2) When the time is fixed, the pore methane pressure in the same matrix block presents low values at both the ends and high values in the middle of the matrix block along the direction of the methane driven by water. The reason is that the methane in the matrix block flows from the centre of the matrix to the cleat or hydraulic fracture inside both the sides. The methane pressure in the middle of the matrix block is the highest, while the methane pressure at both ends of the same matrix block is relatively low. (3) When the time is fixed, along the direction of the methane driven by water, the pore pressure presents a notable abrupt change at the cleat. This phenomenon can be explained by the fact that the methane pressure in the cleat is lower than that in its adjacent upstream and downstream matrix blocks. (4) The pore pressure, in a certain range of the fracturing borehole, increases with time progression, while the pore pressure outside this range decreases with time progression.

Figure 6.

Space–time distribution curve of pore pressure.

Space–time distribution rules of pore pressure gradient

The pore pressure gradient on the monitoring line at times 0.1 h, 0.5 h, 1 h and 1.5 h is shown in Figure 7. The pore pressure gradients near the borehole, at the cleat or fracture and the boundary of the coal seam change significantly (Figure 7(a)). This phenomenon can be attributed to the regional characteristics of the pore pressure gradient in the direction of the methane driven by hydraulic fracturing.

Figure 7.

Curves of pore pressure gradient along the horizontal monitoring line.

Figure 7(b) shows that (1) when the time is fixed, the pore pressure gradient of the space point near the fracturing borehole gradually decreases along the direction of the methane driven by hydraulic fracturing. This phenomenon occurs because inside the coal matrix close to the fracturing borehole, a larger distance from the fracturing borehole corresponds to a weaker degree of methane compression; in addition, a lower methane pressure corresponds to a smaller methane pressure difference between two adjacent space points. (2) When the location of the space point is certain, the pore pressure gradient at the spatial point gradually increases with time progression. This finding is explained by the fact that with the continuous process of methane driven by hydraulic fracturing, the methane pressure on the space point increases under the influence of the water pressure, resulting in the increase in the methane pressure difference between the space point and its adjacent space point.

It can be seen from Figures 7(c), 7(d) and 7(e) that (1) when the time is fixed, the pore pressure gradient at the space point in the cleat or fracture is considerably lower than that inside its adjacent matrix block. This pore pressure gradient difference becomes the driving force of the methane flowing from the matrix block to the cleat or fracture. (2) When the position of the space point is fixed, the pore pressure gradient at this point decreases with time progression. The reason is that the free methane at this point constantly flows into the cleat or fracture of the coal seam, and the methane content in the coal matrix decreases as well. As a result, the methane pressure in the coal matrix decreases, eventually resulting in a decrease in the pore pressure gradient with time progression.

It can be seen from Figure 7(f) that (1) at times 0.1 h, 0.5 h, 1.0 h and 1.5 h, the corresponding pore pressure gradients on the 0.1 m space distance increase by 4.765 MPa/m, 4.311 MPa/m, 3.652 MPa/m and 3.264 MPa/m. This aspect indicates that when the time is fixed, the pore pressure gradient near the boundary of the coal seam exhibits a notable trend of increase along the direction of the methane driven by the hydraulic fracturing. (2) At times 0.1 h, 0.5 h, 1.0 h and 1.5 h, the pore pressure gradients at the 10 m space point are 5.814 MPa/m, 5.392 MPa/m, 4.741 MPa/m and 4.335 MPa/m, respectively. This aspect indicates that when the space point near the boundary of the coal seam is fixed, the pore pressure gradient at this point decreases with time progression.

Space–time distribution rules of water and methane saturation during the process of methane driven by hydraulic fracturing

(1) Visualization of space–time distribution of water-methane two-phase saturation

The space–time distributions of water and methane saturation at times 0.1 h, 0.5 h, 1.0 h and 1.5 h are shown in Figure 8. It can be observed from Figure 8(a) that (1) the borehole is located inside the matrix block, and the water saturation in the matrix block around the borehole is increasing. This finding is apparently a result of the continuous flow of water into the interior of the coal matrix. (2) Inside the coal matrix block, the water saturation in its central region is lower than the water saturation in its surrounding region. This finding occurs because when water penetrates the matrix block through cleats or hydraulic fractures, the water first infiltrates the surrounding matrix block and later flows into its central matrix block gradually. (3) The water saturation in the entire matrix block is larger than that in its surrounding fracture or cleat. The reason is that after water enters the inner part of the matrix block, the methane inside the matrix block is driven into the cleat or fracture surrounding the matrix block. As a result, the space occupied by water in the matrix space is greater than that in the cleat or fracture. Because the sum of the water and methane saturation is always 1, the space–time distribution rules of methane saturation (Figure 8(b)) are opposite to those of water saturation.

Figure 8.

Space–time distribution of water and methane saturation.

(2) Space–time distribution of water and methane saturation

The space–time distribution of water saturation on the monitoring line at times 0.1 h, 0.5 h, 1.0 h and 1.5 h is shown in Figure 9. It can be seen from Figure 9(a) that (1) when the time is constant, the water saturation demonstrates a decreasing trend along the direction of methane driven by hydraulic fracturing. This phenomenon occurs because when water is continuously injected into the coal seam from the fracturing borehole, the water flows from the closest point to a distant point, and the fracturing borehole is in the centre of these points. The closer it is to the fracturing borehole, the more space water occupies. The outward performance is that the water saturation gradually decreases along the direction of the methane driven by hydraulic fracturing in space. (2) When the time is fixed, the water saturation demonstrates a trend of first increasing and later decreasing near the cleat or hydraulic fracture along the direction of methane driven by hydraulic fracturing. This aspect indicates that when the space position transits from the upstream matrix block to the adjacent downstream matrix block through the cleat or hydraulic fracturing, the water saturation first increases and later decreases. The reason is that the water saturation of the upstream and downstream matrix blocks is different from that of the cleat or fracture along the direction of the methane driven by hydraulic fracturing. (3) To a certain position, water saturation gradually increases with time progression. This phenomenon is explained by the fact that with time progression, the amount of water flowing through the point increases, resulting in an increase in the water saturation at the space point. Similarly, because the sum of the methane and water saturation is always 1, the space–time distribution rules of methane saturation (Figure 9(b)) are contrary to those of water saturation.

Figure 9.

Space–time distribution curve of water and methane saturation.

Variation of methane content in coal seam over time with methane driven by hydraulic fracturing

In the process of methane driven by hydraulic fracturing, the curve of the total methane content, adsorbed methane content and free methane content of the coal seam over time is shown in Figure 10. It can be seen from Figure 10 that the free methane content, adsorbed methane content and total methane content decrease with the time progression, and the reduction rate of the adsorbed methane content is higher than that of the free methane content. This phenomenon occurs because during the process of methane driven by hydraulic fracturing, the free methane in the coal seam is first driven out by water. Subsequently, the dynamic equilibrium state of the free and adsorbed methane in the coal seam is disrupted, resulting in the adsorbed methane transforming into free methane. In addition, the hydraulic fracturing reduces the methane pressure in some areas of the coal seam, resulting in the transformation of the adsorbed methane into free methane, which is also driven out by the high-pressure water. The decrease in the free and adsorption methane finally leads to the decrease in the total methane content in the coal seam.

Figure 10.

Change in methane content over time during the process of methane driven by hydraulic fracturing.

Conclusions

During the process of methane driven by hydraulic fracturing, the pore pressure is redistributed. The methane pressure near the fracturing hole is higher than the original methane pressure; however, the methane pressure at the boundary of the coal seam is lower than the original methane pressure. The methane pressure in the hydraulic fracture and the cleat connected with it is higher than that in their surrounding area. The pore pressure exhibits low values at both ends and high values in the middle of the same matrix block along the direction of methane driven by hydraulic fracturing; furthermore, a notable abrupt change appears at the cleat. The pore pressure at the same space point near the fracturing borehole increases with time; however, the pore pressure at the same spatial point outside the range decreases with time progression.

Along the direction of methane driven by hydraulic fracturing, the pore pressure gradient in the coal seam exhibits a certain regional characteristic. The pore pressure gradient near the fracturing borehole demonstrates a notable decreasing trend along the direction of methane driven by hydraulic fracturing. The pore pressure gradient at the spatial point of the cleat is significantly lower than that in its adjacent matrix blocks. However, the pore pressure gradient near the boundary of the coal seam exhibits an obvious increasing trend. The pore pressure gradient at the same space point near the fracturing borehole and the boundary of the coal seam decreases with the time progression.

With the help of the COMSOL Multiphysics software, the space–time distribution of water and methane saturation in the coal seam matrix and cleat or fracture can be observed directly, and the visualization of the process of methane driven by hydraulic fracturing can be realized. From the cloud map of the water and methane saturation distribution, it can be noted that along the direction of the methane driven by hydraulic fracturing, the water saturation in the coal seam exhibits an overall decreasing trend; however, it first increases and later decreases near the cleat or hydraulic fracture. The water saturation at the inner space point of the coal seam increases gradually with the time progression. Inside the matrix block, the water saturation in the central region is less than that in the surrounding region, and the water saturation in the entire matrix block is greater than that in the cleat or hydraulic fracture surrounding the matrix block. The space–time distribution rules of methane saturation are contrary to those of water saturation.

The free methane driven by the hydraulic fracturing includes the original free methane and the free methane desorbed from the adsorbed methane. The reduction rate of the adsorbed methane is larger than that of the free methane.

Footnotes

Data Availability

The data in the manuscript can be available on request through Weiyong Lu, whose email address is 489698551@qq.com.

Declaration of Conflicting Interests

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

Funding

The author(s) declared the following potential conflicts of interest with respect to the research, authorship, and/or publication of this article: This study is funded by the National Natural Science Foundation of China (No. 51704285 and 51774272), Scientific and Technological Innovation Programs of Higher Education Institutions in Shanxi (STIP) (No. 2019L0938), Special project of 2019 plan for the introduction of high-level scientific and technological talents in development zone of Lvliang city (Development of automatic disassembly platform for hydraulic support pin shaft), Science and technology project of Lvliang city in 2019 (Pressure relief and permeability improvement technology by integrated hydraulic flushing and cutting for low permeability coal seam containing methane), Natural science foundation of Shanxi Province-Mechanical behaviour and breaking mechanism of three-axis fracture of weak and fractured coal body in coalbed methane mining (No. 201901D111330), the Key Innovation Team for Fault Diagnosis and Health Management Technology of Coal Mine Machinery Equipment of ‘1331 Project’ in Shanxi Province (No. TD201812), Research Start-up Funding of Doctors in Luliang University (1901010145) and Transformation of Scientific and Technological Achievements Programs of Higher Education Institutions in Shanxi (TSTAP), which is gratefully acknowledged.

ORCID iD

Weiyong Lu

References

Bertrand

Cerfontaine

Collin

(2017) A fully coupled hydro-mechanical model for the modeling of coalbed methane recovery. Journal of Natural Gas Science and Engineering, 46: 307–325.

Chen

Pan

Liu

, et al. (2013) An improved relative permeability model for coal reservoirs. International Journal of Coal Geology 109–110: 45–57.

Chen

SL.

Experimental Investigation on Methane Driving Caused by Hydraulic Fracturing in Gassy Coal Seams. Xuzhou: China University of Mining and Technology, 2015.

Cheng

QY.

Research on Permeability Improvement and Methane Driven Effect of Hydraulic Fracturing for Low Permeability Coal Seam. Xuzhou: China University of Mining and Technology, 2012.

Chong

Chen

, et al. (2017) Numerical investigation into the effect of natural fracture density on hydraulic fracture network propagation. Energies 10(7): 914.

Corey

Rathjens

(1956) Effect of stratification on relative permeability. Journal of Petroleum Technology 8(12): 69–71.

Fan

Luo

, et al. (2016) Deep CBM extraction numerical simulation based on hydraulic-mechanical-thermal coupled model. Journal of China Coal Society 41(12): 3076–3085.

Feng

Zhang

, et al. (2015) Research and application of fracture parameter optimization of fractured horizontal well in low permeability gas reservoir. Fault-Block Oil & Gas Field 21(05): 600–605.

Hosking

Thomas

Sedighi

(2018) A dual porosity model of high-pressure gas flow for geoenergy applications. Canadian Geotechnical Journal 55(6): 839–851.

10.

Huang

(2018) A fluid-solid coupling mathematical model of methane driven by water in porous coal. Geofluids 2018(9): 1–17.

11.

Huang

Cheng

Chen

(2016) Phenomenon of methane driven caused by hydraulic fracturing in methane-bearing coal seams. International Journal of Mining Science and Technology 26(5): 919–927.

12.

Huang

Chen

Cheng

(2016) Basic problems of hydraulic fracturing for mining and control zone gas in coal seams. Journal of China Coal Society 41(01): 128–137.

13.

Huang

Cheng

Chen

, et al. (2013) Study of coal seam outburst mitigation by deep hole hydro-fracturing and shallow hole methane drainage. Journal of China University of Mining & Technology 42(05): 701–711.

14.

Jamison

Azad

(2017) The hydraulic fracture–natural fracture network configuration in shale reservoirs: Geological limiting factors. Journal of Petroleum Science and Engineering 159: 205–229.

15.

King

Ertekin

Schwerer

(1986) Numerical simulation of the transient behavior of coal-seam degasification wells. SPE Formation Evaluation 1(2): 165–183.

16.

Fan

Han

, et al. (2016) A fully coupled thermal-hydraulic-mechanical model with two-phase flow for coalbed methane extraction. Journal of Natural Gas Science and Engineering 33: 324–336.

17.

Wang

(2016) Numerical simulation method of shale gas reservoirs after stimulated reservoir volume fracturing based on triple porous media model. Petroleum Geology and Recovery Efficiency 23(06): 105–111.

18.

Liu

Wang

Gao

(2016) Flow consistency between non-Darcy flow in fracture network and nonlinear diffusion in matrix to gas production rate in fractured shale gas reservoirs. Transport in Porous Media 111(1): 97–121.

19.

Liu

Lin

Yang

(2017) Impact of matrix–fracture interactions on coal permeability: Model development and analysis. Fuel 207: 522–532.

20.

Huang

Zhao

(2019) A review of recent research and development of the effect of hydraulic fracturing on gas adsorption and desorption in coal seams. Adsorption Science & Technology 37(5–6): 509–529.

21.

Palmer

(2009) Permeability changes in coal: Analytical modeling. International Journal of Coal Geology 77(1–2): 119–126.

22.

Preisig

Eberhardt

Gischig

, et al. (2015) Development of connected permeability in massive crystalline rocks through hydraulic fracture propagation and shearing accompanying fluid injection. Geofluids 15(1–2): 321–337.

23.

Remner

Ertekin

Sung

(1984a) A parametric study of the effects of coal seam properties on gas drainage efficiency. SPE Reservoir Engineering 1(6): 633–646.

24.

Remner

Ertekin

Sung

(1984b) A scalable massively parallel dual-porosity dual-permeability simulator for fractured reservoirs with super-k permeability. SPE Reservoir Engineering 1(6): 633–646.

25.

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.

26.

Yun

Rim

Han

(2018) A model for pseudo-steady and non-equilibrium sorption in coalbed methane reservoir simulation and its application. Journal of Natural Gas Science and Engineering 54: 342–348.

27.

Zhang

(2014) Numerical simulation of hydraulic fracturing coalbed methane reservoir. Fuel 136: 57–61.

28.

Zhang

Bian

(2015) Numerical simulation of hydraulic fracturing coalbed methane reservoir with independent fracture grid. Fuel 143: 543–546.

29.

Zhang

JX.

Research on Dynamic Coupling Model and Its Application Base on the Pore-Fracture Dual Media Characteristics of Coal Containing. Jiaozuo: Henan Polytechnical University, 2016.

30.

Zhang

Ranjith

Perera

MSA

, et al. (2016) Gas transportation and enhanced coalbed methane recovery processes in deep coal seams: A review. Energy & Fuels 30(11): 8832–8849.

31.

Zhang

Xie

, et al. (2017) An anisotropic coal permeability model that considers mining-induced stress evolution, microfracture propagation and gas sorption-desorption effects. Journal of Natural Gas Science and Engineering 46: 664–679.

32.

Zou

Wei

, et al. (2015) Modeling and application of coalbed methane recovery performance based on a triple porosity/dual permeability model. Journal of Natural Gas Science and Engineering 22: 679–688.