# Transient Coupled Flow Model for Matrix-Fracture in High-Pressure Gas Reservoirs Opened by Gas Drilling

Authors: Jiajie Liu; Gao Li; Zhi Ye; Wang He; Zhihao Zeng
Journal: Geofluids
 (2023)

---

## Abstract

Gas drilling technology is an effective method to protect low permeability tight sandstone reservoirs from damage. However, using this technology to drill high-pressure fractured reservoirs, the fluid flow characteristics are extremely complex, leading to drilling safety challenges. Therefore, based on the conservation equations of mass, momentum, and energy, a new transient coupled flow model in matrix-fracture is established for high-pressure gas reservoirs opened by gas drilling, considering convection term and energy equation, to research the pressure and flow velocity transient behaviors. This model is solved by Shaher Momani’s algorithm and CyabJiBe format, and the type curves of the pressure and flow velocity are plotted. The three flow stages are developed, and the impact of parameters and equation terms on the transient flow behavior is analyzed and discussed, such as the energy equation, the matrix permeability (k), the fracture aperture (Lf), and the fracture outlet pressure. It is found that the energy equation term of the model has a great influence on the transient characteristics such as velocity and pressure, which cannot be ignored. In the first flow stage, the flow velocity characteristics are not affected by the matrix parameters being extremely unstable, but in the second and third stages, the influence of the gas supplied to fracture cannot be ignored. The proposed transient coupled flow model provides a new understanding of the transient flow behavior in matrix-fracture for high-pressure reservoirs during gas drilling, which can be used to interpret the characteristics of pressure and velocity more realistically, to provide support for safe drilling of gas drilling.


---

## Body

## 1. Introduction

The development of high-pressure tight gas reservoirs is an unsolved problem in the petroleum industry due to the low productivity and high pressure of gas wells. Producing gas from tight reservoirs presents a unique challenge due to the contamination of water from the fracturing and drilling fluids [1–3]. As an advanced drilling technology, using air as a drilling fluid began in the 1950s [4–6]. Gas drilling technology has shown to be promising to solve the problem, which is characterized by the easier discovery of reservoirs, low reservoir damage, and increased rate of penetration [7, 8]. However, there are many problems during practice, especially the safety of gas drilling, which is very uncontrollable in many areas due to the lack of understanding of the transient flow law and pressure characteristics of high-pressure gas, when the high-pressure reservoir is opened by gas drilling [9, 10]. Therefore, it is very important for gas drilling engineers to understand the flow characteristics of gas in fractured high-pressure reservoirs.A matrix-fracture medium model is commonly used for natural fracture gas reservoirs. Warren and Root have developed a dual-porosity medium model, which is based on the assumption that the matrix is the main storage space, the natural fractures are the main seepage channels, and the matrix provides a gas source to fracture [11–13]. Kazemi et al. have developed a double-porosity dual-permeability model; the gas of the matrix system and small-scale fractures not only can provide a gas source to fracture but also can provide a gas source to the wellbore [14]. Although there are abundant studies of the coupled flow behavior for matrix-fracture reservoirs, few models research the characteristics of unstable complex flow for high-pressure reservoirs during gas drilling.For naturally fractured gas reservoirs, most researchers have focused on transient flow behaviors. Slimani and Tiab [15] established a theoretical model, to investigate the influence of transient flow behaviors on type curves for partially penetrating vertical wells in naturally fractured reservoirs. Wu et al. researched a triple-porosity medium model to analyze the fluid flow behavior in matrix-fracture-cave reservoirs [16]. Du et al. conducted a multiregion composite model, to investigate the transient flow behavior for fractured-cave reservoirs [17]. For composite fractured gas reservoirs, Nie et al. [18] studied a transient flow model with the high-velocity non-Darcy effect, to simulate the non-Darcy flow for the inner region as a dual-porosity medium, and the Darcy flow for the outer region as a homogeneous medium.With the development of hydraulic fracturing technology, which is one of the effective ways to develop tight sandstone and shale reservoirs with low permeability characteristics, many researchers attach importance to the study of complex fluid behaviors in such reservoirs [19]. For a hydraulic fractured sandstone tight reservoir, a series of new fluid flow models are presented to help engineers understand the transient pressure characteristics of gas flow. A new mathematical model has been proposed to study the transient pressure characteristics for hydraulic fractured wells with nonuniform fracture geometry and permeability [20]. The semianalytical model has been established and solved to analyze the transient pressure behavior for a fractured vertical well with hydraulic/natural fracture networks by considering stress-sensitive effects [21]. The model for multistage fractured horizontal wells with discretely distributed natural fractures has been researched to investigate pressure performance [22]. The new fractally fractional diffusion model of composite dual-porosity has been developed to evaluate the performance of multiple fractured horizontal wells with stimulated reservoir volume in tight gas reservoirs [23]. The discrete fracture-matrix method based on the numerical well testing model has been proposed to study the pressure transient responses of discretely distributed natural fractures [24]. Due to the characteristics of ultralow porosity and ultralow permeability in shale gas fields, which contain both adsorbed gas and free gas, hydraulic fracturing technology and horizontal well technology are commonly used in shale gas development. In recent years, discrete fracture models have been vigorously developed to study gas seepage in gas reservoirs in shale [25, 26], and many scholars believe that the changes in matrix permeability and fracture permeability with increasing effective stress cannot be ignored [27, 28]. A series of multiple-porosity model transient flow models are proposed, including the dual-porosity dual-permeability model and triple-porosity model, to study the transient pressure behaviors [29, 30]. Although there are abundant studies of transient flow behavior for naturally and hydraulically fractured reservoirs, few models mentioned above have touched on the convection term and considered the influence of the energy equation.In fact, owing to the high pressure of the reservoirs and high conductivity of the fractures, the transient flow of gas cannot ignore the influence of the gas convection term and energy equation. Based on the law of conservation of mass, momentum, and energy, a transient coupled flow in matrix-fracture for the high-pressure reservoirs opened by gas drilling model is established. The flow in fracture is considered as a radial flow to simulate the convection and diffusion process of gas, and the flow in matrix reservoir is considered a parallel seepage into fracture to simulate the non-Darcy flow, as shown in Figure1. By using Shaher Momani’s algorithm and CyabJiBe format, the analytical solution of the coupled transient flow in matrix-fracture model can be obtained [31–33]. Then, the type curves with the influence of energy conversion and physical characteristics of matrix-fracture are plotted.Figure 1
Physical model of gas flow process in matrix-fracture reservoirs opened by gas drilling.


## 2. Physical Model and Assumptions

Figure1 is the schematic of the physical model, which represents the gas transient flow process in fractured reservoirs opened by gas drilling with no damage. The gas flow in reservoirs is divided into two regions, which have different flow characteristics and different reservoir characteristics. The first gas flow area is the gas transient radial flow at high velocity from fracture to wellbore, considering the convective and diffusion terms. The second gas flow area is the gas transient parallel flow perpendicular to the fracture from the matrix. The main assumptions of this model are the following: (1) The fracture has the same geometric characteristics. (2) In the second flow area, the matrix with homogeneous characteristics provides a gas source to fracture vertically. (3) The permeability of the matrix is independent of pressure. (4) The influence of the secondary gradient on the flow in the matrix is ignored.


## 3. Mathematical Model

### 3.1. Mathematical Model for the Gas Transient Parallel Flow Perpendicular to Fracture from Matrix

Owing to the high pressure of the reservoirs, the influence of non-Darcy flow cannot be ignored for the gas flow in the matrix; the equation is as follows [18]:
(1)−dpdzm=Kμv+βρvn.In Equation (1), K represents the selected reservoir matrix permeability. v and μ represent gas velocity and viscosity. ρ represents gas density. n represents the value related to the properties of the porous medium. β represents the non-Darcy flow factor.For plane flow in homogeneous formation, Equation (1) can be rewritten as follows:
(2)v=−δKμ∂p∂zm.In Equation (2), δ represents the inertial-turbulence correction coefficient factor and zm represents vertical distance.The continuity equation of gas transient parallel flow perpendicular to fracture from the matrix is as follows:(3)∂ρϕ∂t+∂ρv∂zm=0.Binding Equation (2), Equation (3) can be written as follows (the detailed derivation of Equation (4) is presented in the appendix:
(4)∂p2∂zm2=ϕμCgδK∂p∂t.


### 3.2. Mathematical Model for the Gas Transient Radial Flow from Fracture to Wellbore

In the study of transient radial flow in fracture, the influence of the convection term and energy equation is mostly ignored; the transient flow characteristics cannot be truly reflected. Therefore, the convection term and energy equation are introduced to describe the transient flow characteristics of natural gas in fractures, during the instantaneous opening of fractured reservoirs by gas drilling. The transient mathematical model of gas transient radial flow in fractures is established.The governing equations of the model are as follows (the detailed derivation of Equations (5)–(7) is presented in the appendix) [34–36].The continuity equation is as follows:(5)∂ρ∂t+∂rρvrr∂r=0.The momentum conservation equation is as follows:(6)∂ρvr∂t+∂rρvrvr+p−τrrr∂r=pr.The energy conservation equation is as follows:(7)∂ρE∂t+∂rρHvr−vrτrr+qrr∂r=ρq.In the equation, energyE=CvT+1/2vr2, enthalpy H=CpT+1/2vr2=E+p/ρ, radial flow qr=−λg∂T/∂r, and tangential stress τrr=2/3μ2∂rvr/r∂r−2μvr/r.


### 3.3. Mathematical Model for Transient Coupled Flow Model in Matrix-Fracture

The transient coupled flow model in matrix-fracture is applied to describe and analyze the coupling flow characteristics of gas in matrix and fracture, when the fractured reservoir is opened by gas drilling, including the radial flow characteristics of gas from fracture to wellbore and the flow characteristics of gas flowing vertically from matrix to fracture. The gas source is continuously supplied from matrix to fracture with the decrease of fracture pressure. And the convection term and energy equation are introduced to the transient coupled flow model. The physical model is shown in Figure1.According to Equations (5)–(7), governing equations of the model are shown in Equations (8)–(10).Continuity equation for transient coupled flow in the matrix-fracture system [34–36]:
(8)∂ρ∂t+∂rρvrr∂r=q̇.In Equation (8), q̇ represents the gas source which is continuously provided from matrix to fracture.Momentum conservation equation for transient coupled flow in the matrix-fracture system:(9)∂ρvr∂t+∂rρvrvr+p−τrrr∂r=pr.The energy conservation equation for transient coupled flow in the matrix-fracture system:(10)∂ρE∂t+∂rρHvr−vrτrr+qrr∂r=ρq.


### 3.4. Mathematical Model Solutions

In this section, Equations (7)–(9) are rewritten as the vector form of the convection diffusion equation; the mathematical model is as follows:
(11)∂u∂t=−a∂u∂r+b∂2u∂r2+fr,t.The vector parameteru,a,b,c,fr,t in Equation (11) are as follows:
(12)u=pvrT,a=vrρ01pvr−23μrρ00p+23μvrrpvrCv,b=000043μρ0000,fx,t=q̇pρ−pvrr0pq−pvrρr.Shaher Momani’s algorithm and CyabJiBe format for solving the convection-diffusion equation with a nonlinear source term is used to address the mathematical model. The basic definitions of the algorithm are as follows.Definition 1.
(13)Jfr=1Γ∫ozftdt.
In the equation,J can be found in [37].Definition 2.
(14)Dfx=JmDmfx=1Γm∫0xx−tm−2fmtdt.
Differential operatorD is described in detail in Reference [38].Definition 3.
(15)Lx=∂∂x,Lxx=∂2∂x2,Dt=∂∂t.Shaher Momani’s algorithm for solving the convection-diffusion equation consists in finding the solution in the form:(16)ur,t=∑n=0∞unr,t.The operatorJ and the inverse of the operator D are applied to solve Equation (11):
(17)ur,t=∑k=0m−1∂ku∂tkr,0+tkk!+Jfr,t+bJLrrur,t−aJLrur,t.The componentsun recursively is defined as follows:
(18)u0r,t=∑k=0m−1∂ku∂tkr,0+tkk!+Jfr,t,(19)u1r,t=bJLrru0−aJLru0,(20)un+1r,t=bJLrrun−aJLrun.The exact series solution Equation (16) is defined by Equations (18)–(20), which give an approximate solution to Equation (11). Finally, the standard characteristic curves of fluid velocity with time and pressure with space in fracture can be plotted by Mathcad as computing software, when the fractured reservoir is opened by gas drilling.


## 3.1. Mathematical Model for the Gas Transient Parallel Flow Perpendicular to Fracture from Matrix

Owing to the high pressure of the reservoirs, the influence of non-Darcy flow cannot be ignored for the gas flow in the matrix; the equation is as follows [18]:
(1)−dpdzm=Kμv+βρvn.In Equation (1), K represents the selected reservoir matrix permeability. v and μ represent gas velocity and viscosity. ρ represents gas density. n represents the value related to the properties of the porous medium. β represents the non-Darcy flow factor.For plane flow in homogeneous formation, Equation (1) can be rewritten as follows:
(2)v=−δKμ∂p∂zm.In Equation (2), δ represents the inertial-turbulence correction coefficient factor and zm represents vertical distance.The continuity equation of gas transient parallel flow perpendicular to fracture from the matrix is as follows:(3)∂ρϕ∂t+∂ρv∂zm=0.Binding Equation (2), Equation (3) can be written as follows (the detailed derivation of Equation (4) is presented in the appendix:
(4)∂p2∂zm2=ϕμCgδK∂p∂t.


## 3.2. Mathematical Model for the Gas Transient Radial Flow from Fracture to Wellbore

In the study of transient radial flow in fracture, the influence of the convection term and energy equation is mostly ignored; the transient flow characteristics cannot be truly reflected. Therefore, the convection term and energy equation are introduced to describe the transient flow characteristics of natural gas in fractures, during the instantaneous opening of fractured reservoirs by gas drilling. The transient mathematical model of gas transient radial flow in fractures is established.The governing equations of the model are as follows (the detailed derivation of Equations (5)–(7) is presented in the appendix) [34–36].The continuity equation is as follows:(5)∂ρ∂t+∂rρvrr∂r=0.The momentum conservation equation is as follows:(6)∂ρvr∂t+∂rρvrvr+p−τrrr∂r=pr.The energy conservation equation is as follows:(7)∂ρE∂t+∂rρHvr−vrτrr+qrr∂r=ρq.In the equation, energyE=CvT+1/2vr2, enthalpy H=CpT+1/2vr2=E+p/ρ, radial flow qr=−λg∂T/∂r, and tangential stress τrr=2/3μ2∂rvr/r∂r−2μvr/r.


## 3.3. Mathematical Model for Transient Coupled Flow Model in Matrix-Fracture

The transient coupled flow model in matrix-fracture is applied to describe and analyze the coupling flow characteristics of gas in matrix and fracture, when the fractured reservoir is opened by gas drilling, including the radial flow characteristics of gas from fracture to wellbore and the flow characteristics of gas flowing vertically from matrix to fracture. The gas source is continuously supplied from matrix to fracture with the decrease of fracture pressure. And the convection term and energy equation are introduced to the transient coupled flow model. The physical model is shown in Figure1.According to Equations (5)–(7), governing equations of the model are shown in Equations (8)–(10).Continuity equation for transient coupled flow in the matrix-fracture system [34–36]:
(8)∂ρ∂t+∂rρvrr∂r=q̇.In Equation (8), q̇ represents the gas source which is continuously provided from matrix to fracture.Momentum conservation equation for transient coupled flow in the matrix-fracture system:(9)∂ρvr∂t+∂rρvrvr+p−τrrr∂r=pr.The energy conservation equation for transient coupled flow in the matrix-fracture system:(10)∂ρE∂t+∂rρHvr−vrτrr+qrr∂r=ρq.


## 3.4. Mathematical Model Solutions

In this section, Equations (7)–(9) are rewritten as the vector form of the convection diffusion equation; the mathematical model is as follows:
(11)∂u∂t=−a∂u∂r+b∂2u∂r2+fr,t.The vector parameteru,a,b,c,fr,t in Equation (11) are as follows:
(12)u=pvrT,a=vrρ01pvr−23μrρ00p+23μvrrpvrCv,b=000043μρ0000,fx,t=q̇pρ−pvrr0pq−pvrρr.Shaher Momani’s algorithm and CyabJiBe format for solving the convection-diffusion equation with a nonlinear source term is used to address the mathematical model. The basic definitions of the algorithm are as follows.Definition 1.
(13)Jfr=1Γ∫ozftdt.
In the equation,J can be found in [37].Definition 2.
(14)Dfx=JmDmfx=1Γm∫0xx−tm−2fmtdt.
Differential operatorD is described in detail in Reference [38].Definition 3.
(15)Lx=∂∂x,Lxx=∂2∂x2,Dt=∂∂t.Shaher Momani’s algorithm for solving the convection-diffusion equation consists in finding the solution in the form:(16)ur,t=∑n=0∞unr,t.The operatorJ and the inverse of the operator D are applied to solve Equation (11):
(17)ur,t=∑k=0m−1∂ku∂tkr,0+tkk!+Jfr,t+bJLrrur,t−aJLrur,t.The componentsun recursively is defined as follows:
(18)u0r,t=∑k=0m−1∂ku∂tkr,0+tkk!+Jfr,t,(19)u1r,t=bJLrru0−aJLru0,(20)un+1r,t=bJLrrun−aJLrun.The exact series solution Equation (16) is defined by Equations (18)–(20), which give an approximate solution to Equation (11). Finally, the standard characteristic curves of fluid velocity with time and pressure with space in fracture can be plotted by Mathcad as computing software, when the fractured reservoir is opened by gas drilling.


## 4. Results and Discussion

### 4.1. Flow Stage Analysis

The conditions of the model are as follows.The outer boundary pressure of the matrix reservoir is constant pressurepi=60 MPa. The inner boundary pressure of the matrix pmlz=0=p is the dynamic boundary pressure, which is the outer boundary pressure of the fracture. The inner boundary pressure of the fracture is the constant pressure pmlr=0=58 MPa. The matrix permeability is 0.01 mD.Figure2 shows the velocity behavior characteristics of coupled flow for matrix-fracture in a high-pressure reservoir for gas drilling. The three flow stages can be represented by curves. Figure 3 shows the amplification characteristic curve of the first flow stage.Figure 2
Three-stage characteristics of transient coupled flow in matrix-fracture.Figure 3
Magnification curve of the first and second stage characteristics of transient coupled flow in matrix-fracture.The first stage describes gas velocity characteristic in the fractured reservoir, which is opened by gas drilling instantaneously. The magnification curve is shown in Figure3. In this flow stage, the gas flow velocity in fracture is much higher than that of the second and third stages, it is extremely unstable, and it may have a huge impact on the safety of the bottom hole. The flowing gas mainly comes from the gas stored in fracture, and the matrix will not provide gas to the fracture. Because the total amount of gas stored in the fracture of reservoirs is much less than that in the matrix, the duration of the flow process is extremely short, which is less than 1 second.In the second stage, the flowing gas in the fracture mainly includes the gas from the reservoir in the fracture and the gas supplied by the matrix to the fracture. The gas flowing velocity drops rapidly, and it is quasisteady flow. The third stage reflects the flow characteristics of the stable flow stage. The flowing gas in the fracture is supplied mainly from the matrix. In this process, the pressure and the flow velocity drop very slowly.Figure4 shows the characteristics of the pressure funnel in fracture from 0.05 seconds to 1 second after being opened by gas drilling. The pressure attenuation mainly exists in the fracture near the wellbore, that is, within 10 m from the wellbore. The closer to the wellbore, the more serious the pressure attenuation. As time goes on, the pressure in the fracture becomes lower and lower. The pressure drop is mainly consumed near the wellbore, during the fluid flowing from the reservoir edge to the bottom of the well.Figure 4
Radial pressure drop funnel in fracture.


### 4.2. Impact of the Parameters

In this section, the characteristic parameters of the transient coupled flow model in the matrix-fracture reservoir are analyzed in detail, including the influence of the energy conversion, the matrix permeabilityk, the fracture aperture Lf, and the fracture outlet pressure p. These parameters have different influences on the characteristics of gas flow velocity and pressure in the fracture, as shown in Figures 5–9.Figure 5
Influence of energy conversion equation.Figure 6
Influence of matrix permeability (k).Figure 7
Influence of matrix permeability (k) on pressure drop funnel in fractures.Figure 8
Influence of fracture apertureLf.Figure 9
Influence of fracture outlet pressure.


#### 4.2.1. Influence of the Energy Conversion

Figure5 compares the influence of considering energy equation and without energy equation on curves of the transient coupled flow in matrix-fracture model. As a compressible fluid, the gas flow process has the conversion of kinetic energy and internal energy under the condition of high pressure difference and high flow velocity, as shown in Figure 5, by comparing the flow velocity in fracture without considering the energy conversion and considering the energy conversion. In the first stage, it is found that if the energy conversion process is ignored, the flow velocity cannot peak rapidly when the fracture is opened by gas drilling but rises slowly, and the characteristics of the curve cannot reflect the real flow situation. In the second and third flow stages, the characteristics of flow velocity without considering energy conversion are higher than those considering energy conversion. In the condition of high velocity of gas flow in the fracture, the part of kinetic energy can be converted into internal energy. According to the law of energy conservation, the velocity of gas flow decreases due to the reduction of kinetic energy. Therefore, the effect of energy conversion cannot be ignored in the study of the characteristics of gas flow at high velocity for high-pressure reservoirs.


#### 4.2.2. Impact of the Matrix Permeability (k)

Figure6 shows the influence of different matrix permeability (k) on the curves of flow velocity with time in fracture simulated by the transient coupled flow model. The values of the matrix permeability (k) are shown in Figure 6. Due to the fracture supplied with a source of gas from the matrix, it can be observed that the matrix permeability mainly influences the second and third flow stages of the type curves for the transient coupled flow model. In these two flow stages, when the matrix permeability (k) is increased, the flow velocity curves rises. On the contrary, if the matrix permeability (k) is decreased, the flow velocity curves goes down. In the first flow stage, the matrix permeability cannot affect the flow characteristic of the transient coupled flow model; the main reason is that the flowing gas mainly comes from the gas source stored in the fracture, while the gas in the matrix does not flow.Figure7 illustrates the impact of matrix permeability (k) on the curves of the pressure drop funnel in fracture for the transient coupled flow model. The values of matrix permeability (k) are shown in Figure 7. The characteristic curves of the pressure drop funnel also shows that matrix permeability (k) affects the attenuation of fracture pressure. The larger the matrix permeability (k), the slower the fracture pressure falls. If the matrix permeability (k) is increased, the supply of gas from the matrix to the fracture can increase, so the rate of pressure drop decreases.


#### 4.2.3. Impact of the Fracture ApertureLf

Figure8 shows the influence of different fracture aperture (Lf) on the curves of flow velocity with time in fracture simulated by the transient coupled flow model. The values of the fracture aperture (Lf) are shown in Figure 8. The curves also illustrate that the fracture aperture (Lf) influences the overall three flow stages. The larger the fracture aperture (Lf), the higher the instantaneous peak for gas production in the first stage and the greater the gas velocity in the second and third fracture flow stages. Therefore, the larger the fracture aperture (Lf) encountered during gas drilling, the greater the impact on the wellbore safety. The reason is that the smaller the fracture aperture, the smaller the flow channel and the greater the flow resistance.


#### 4.2.4. Influence of the Fracture Outlet Pressure

Figure9 illustrates the influence of the fracture outlet pressure on the curves of flow velocity with time in fracture for the transient coupled flow model. The parameters of the fracture outlet pressure are shown in Figure 9. It can be seen in the curves of flow velocity with time; the fracture outlet pressure influences the overall three flow stages during gas drilling. Undoubtedly, the fracture outlet pressure determines the driving force of the flow system, including gas flow system from matrix to fracture and gas flow system from fracture to wellbore. The lower the fracture outlet pressure, the higher the driving force of the fluid system, the more gas supplied from matrix to fracture, and the greater the gas flow velocity in the fracture.


## 4.1. Flow Stage Analysis

The conditions of the model are as follows.The outer boundary pressure of the matrix reservoir is constant pressurepi=60 MPa. The inner boundary pressure of the matrix pmlz=0=p is the dynamic boundary pressure, which is the outer boundary pressure of the fracture. The inner boundary pressure of the fracture is the constant pressure pmlr=0=58 MPa. The matrix permeability is 0.01 mD.Figure2 shows the velocity behavior characteristics of coupled flow for matrix-fracture in a high-pressure reservoir for gas drilling. The three flow stages can be represented by curves. Figure 3 shows the amplification characteristic curve of the first flow stage.Figure 2
Three-stage characteristics of transient coupled flow in matrix-fracture.Figure 3
Magnification curve of the first and second stage characteristics of transient coupled flow in matrix-fracture.The first stage describes gas velocity characteristic in the fractured reservoir, which is opened by gas drilling instantaneously. The magnification curve is shown in Figure3. In this flow stage, the gas flow velocity in fracture is much higher than that of the second and third stages, it is extremely unstable, and it may have a huge impact on the safety of the bottom hole. The flowing gas mainly comes from the gas stored in fracture, and the matrix will not provide gas to the fracture. Because the total amount of gas stored in the fracture of reservoirs is much less than that in the matrix, the duration of the flow process is extremely short, which is less than 1 second.In the second stage, the flowing gas in the fracture mainly includes the gas from the reservoir in the fracture and the gas supplied by the matrix to the fracture. The gas flowing velocity drops rapidly, and it is quasisteady flow. The third stage reflects the flow characteristics of the stable flow stage. The flowing gas in the fracture is supplied mainly from the matrix. In this process, the pressure and the flow velocity drop very slowly.Figure4 shows the characteristics of the pressure funnel in fracture from 0.05 seconds to 1 second after being opened by gas drilling. The pressure attenuation mainly exists in the fracture near the wellbore, that is, within 10 m from the wellbore. The closer to the wellbore, the more serious the pressure attenuation. As time goes on, the pressure in the fracture becomes lower and lower. The pressure drop is mainly consumed near the wellbore, during the fluid flowing from the reservoir edge to the bottom of the well.Figure 4
Radial pressure drop funnel in fracture.


## 4.2. Impact of the Parameters

In this section, the characteristic parameters of the transient coupled flow model in the matrix-fracture reservoir are analyzed in detail, including the influence of the energy conversion, the matrix permeabilityk, the fracture aperture Lf, and the fracture outlet pressure p. These parameters have different influences on the characteristics of gas flow velocity and pressure in the fracture, as shown in Figures 5–9.Figure 5
Influence of energy conversion equation.Figure 6
Influence of matrix permeability (k).Figure 7
Influence of matrix permeability (k) on pressure drop funnel in fractures.Figure 8
Influence of fracture apertureLf.Figure 9
Influence of fracture outlet pressure.


### 4.2.1. Influence of the Energy Conversion

Figure5 compares the influence of considering energy equation and without energy equation on curves of the transient coupled flow in matrix-fracture model. As a compressible fluid, the gas flow process has the conversion of kinetic energy and internal energy under the condition of high pressure difference and high flow velocity, as shown in Figure 5, by comparing the flow velocity in fracture without considering the energy conversion and considering the energy conversion. In the first stage, it is found that if the energy conversion process is ignored, the flow velocity cannot peak rapidly when the fracture is opened by gas drilling but rises slowly, and the characteristics of the curve cannot reflect the real flow situation. In the second and third flow stages, the characteristics of flow velocity without considering energy conversion are higher than those considering energy conversion. In the condition of high velocity of gas flow in the fracture, the part of kinetic energy can be converted into internal energy. According to the law of energy conservation, the velocity of gas flow decreases due to the reduction of kinetic energy. Therefore, the effect of energy conversion cannot be ignored in the study of the characteristics of gas flow at high velocity for high-pressure reservoirs.


### 4.2.2. Impact of the Matrix Permeability (k)

Figure6 shows the influence of different matrix permeability (k) on the curves of flow velocity with time in fracture simulated by the transient coupled flow model. The values of the matrix permeability (k) are shown in Figure 6. Due to the fracture supplied with a source of gas from the matrix, it can be observed that the matrix permeability mainly influences the second and third flow stages of the type curves for the transient coupled flow model. In these two flow stages, when the matrix permeability (k) is increased, the flow velocity curves rises. On the contrary, if the matrix permeability (k) is decreased, the flow velocity curves goes down. In the first flow stage, the matrix permeability cannot affect the flow characteristic of the transient coupled flow model; the main reason is that the flowing gas mainly comes from the gas source stored in the fracture, while the gas in the matrix does not flow.Figure7 illustrates the impact of matrix permeability (k) on the curves of the pressure drop funnel in fracture for the transient coupled flow model. The values of matrix permeability (k) are shown in Figure 7. The characteristic curves of the pressure drop funnel also shows that matrix permeability (k) affects the attenuation of fracture pressure. The larger the matrix permeability (k), the slower the fracture pressure falls. If the matrix permeability (k) is increased, the supply of gas from the matrix to the fracture can increase, so the rate of pressure drop decreases.


### 4.2.3. Impact of the Fracture ApertureLf

Figure8 shows the influence of different fracture aperture (Lf) on the curves of flow velocity with time in fracture simulated by the transient coupled flow model. The values of the fracture aperture (Lf) are shown in Figure 8. The curves also illustrate that the fracture aperture (Lf) influences the overall three flow stages. The larger the fracture aperture (Lf), the higher the instantaneous peak for gas production in the first stage and the greater the gas velocity in the second and third fracture flow stages. Therefore, the larger the fracture aperture (Lf) encountered during gas drilling, the greater the impact on the wellbore safety. The reason is that the smaller the fracture aperture, the smaller the flow channel and the greater the flow resistance.


### 4.2.4. Influence of the Fracture Outlet Pressure

Figure9 illustrates the influence of the fracture outlet pressure on the curves of flow velocity with time in fracture for the transient coupled flow model. The parameters of the fracture outlet pressure are shown in Figure 9. It can be seen in the curves of flow velocity with time; the fracture outlet pressure influences the overall three flow stages during gas drilling. Undoubtedly, the fracture outlet pressure determines the driving force of the flow system, including gas flow system from matrix to fracture and gas flow system from fracture to wellbore. The lower the fracture outlet pressure, the higher the driving force of the fluid system, the more gas supplied from matrix to fracture, and the greater the gas flow velocity in the fracture.


## 4.2.1. Influence of the Energy Conversion

Figure5 compares the influence of considering energy equation and without energy equation on curves of the transient coupled flow in matrix-fracture model. As a compressible fluid, the gas flow process has the conversion of kinetic energy and internal energy under the condition of high pressure difference and high flow velocity, as shown in Figure 5, by comparing the flow velocity in fracture without considering the energy conversion and considering the energy conversion. In the first stage, it is found that if the energy conversion process is ignored, the flow velocity cannot peak rapidly when the fracture is opened by gas drilling but rises slowly, and the characteristics of the curve cannot reflect the real flow situation. In the second and third flow stages, the characteristics of flow velocity without considering energy conversion are higher than those considering energy conversion. In the condition of high velocity of gas flow in the fracture, the part of kinetic energy can be converted into internal energy. According to the law of energy conservation, the velocity of gas flow decreases due to the reduction of kinetic energy. Therefore, the effect of energy conversion cannot be ignored in the study of the characteristics of gas flow at high velocity for high-pressure reservoirs.


## 4.2.2. Impact of the Matrix Permeability (k)

Figure6 shows the influence of different matrix permeability (k) on the curves of flow velocity with time in fracture simulated by the transient coupled flow model. The values of the matrix permeability (k) are shown in Figure 6. Due to the fracture supplied with a source of gas from the matrix, it can be observed that the matrix permeability mainly influences the second and third flow stages of the type curves for the transient coupled flow model. In these two flow stages, when the matrix permeability (k) is increased, the flow velocity curves rises. On the contrary, if the matrix permeability (k) is decreased, the flow velocity curves goes down. In the first flow stage, the matrix permeability cannot affect the flow characteristic of the transient coupled flow model; the main reason is that the flowing gas mainly comes from the gas source stored in the fracture, while the gas in the matrix does not flow.Figure7 illustrates the impact of matrix permeability (k) on the curves of the pressure drop funnel in fracture for the transient coupled flow model. The values of matrix permeability (k) are shown in Figure 7. The characteristic curves of the pressure drop funnel also shows that matrix permeability (k) affects the attenuation of fracture pressure. The larger the matrix permeability (k), the slower the fracture pressure falls. If the matrix permeability (k) is increased, the supply of gas from the matrix to the fracture can increase, so the rate of pressure drop decreases.


## 4.2.3. Impact of the Fracture ApertureLf

Figure8 shows the influence of different fracture aperture (Lf) on the curves of flow velocity with time in fracture simulated by the transient coupled flow model. The values of the fracture aperture (Lf) are shown in Figure 8. The curves also illustrate that the fracture aperture (Lf) influences the overall three flow stages. The larger the fracture aperture (Lf), the higher the instantaneous peak for gas production in the first stage and the greater the gas velocity in the second and third fracture flow stages. Therefore, the larger the fracture aperture (Lf) encountered during gas drilling, the greater the impact on the wellbore safety. The reason is that the smaller the fracture aperture, the smaller the flow channel and the greater the flow resistance.


## 4.2.4. Influence of the Fracture Outlet Pressure

Figure9 illustrates the influence of the fracture outlet pressure on the curves of flow velocity with time in fracture for the transient coupled flow model. The parameters of the fracture outlet pressure are shown in Figure 9. It can be seen in the curves of flow velocity with time; the fracture outlet pressure influences the overall three flow stages during gas drilling. Undoubtedly, the fracture outlet pressure determines the driving force of the flow system, including gas flow system from matrix to fracture and gas flow system from fracture to wellbore. The lower the fracture outlet pressure, the higher the driving force of the fluid system, the more gas supplied from matrix to fracture, and the greater the gas flow velocity in the fracture.


## 5. Conclusions

(1)
A new transient coupled flow model is established for matrix-fracture in high-pressure reservoir during gas drilling, based on the laws of conservation of mass, momentum, and energy(2)
The model considers more realistic assumptions about the convection term and the energy equation and reflects the more real transient flow characteristics(3)
The model is solved by using Shaher Momani’s algorithm and CyabJiBe format. The standard pressure and velocity characteristic curves of the proposed flow model are drawn into three flow regimes. In the first stage, the characteristics of pressure and flow velocity in fracture are extremely unstable, with a rapid emergence of peak, and the gas flowing from the fracture to the wellbore mainly comes from the gas stored in the fracture, without matrix gas. In the second and third flow stages, the flow characteristics are relatively stable(4)
The characteristic parameters of the transient coupled flow model are also analyzed for high-pressure fractured reservoirs opened by drilling. It is found that the pressure and velocity characteristics are mainly affected by the matrix permeability (k), the fracture aperture (Lf), and the fracture outlet pressure. Different characteristic parameters have different influences on the type curves



---

*Source: 9706407-2023-02-10.xml*