Skip to main content

ORIGINAL RESEARCH article

Front. Energy Res., 01 December 2020
Sec. Advanced Clean Fuel Technologies
This article is part of the Research Topic Low-Carbon Technologies for the Petroleum Industry View all 19 articles

Transient Pressure Analysis of Inclined Well in Continuous Triple-Porosity Reservoirs With Dual-Permeability Behavior

Sheng-Zhi Qi,,
Sheng-Zhi Qi1,2,3*Xiao-Hua Tan,
Xiao-Hua Tan1,3*Xiao-Ping Li,Xiao-Ping Li1,3Zhan Meng,Zhan Meng1,3You-Jie Xu,You-Jie Xu1,3Daming TangDaming Tang4
  • 1State Key Laboratory of Oil and Gas Reservoir Geology and Exploitation, Southwest Petroleum University, Chengdu, China
  • 2Energy Engineering College, Longdong Uniersity, Qingyang, China
  • 3Key Laboratory of Seepage and Rock Mechanics of Oil and Gas Wells in Shanxi Province, Xian, China
  • 4Shuguang Oil Production Plant of Liaohe Oilfield, Panjin, China

Inclined wells has recently been adopted to develop fractured-vuggy carbonate hydrocarbon reservoirs (FVCHRs) composed by matrix, fracture and vug system. Therefore, it is significant for us to describe pressure transient of inclined well for FVCHRs. In this paper, it is assumed that vug and fracture system connect with wellbore, and inter-porosity flow from vug to fracture system appears. Therefore, the pressure transient responses model of inclined well with triple-porosity dual-permeability behavior was built. The model is solved by employing Laplace integral transform and finite cosine Fourier transform. Real-domain solution of the model is obtained by Stehfest inversion algorithm. On the basis model of the published paper, the solution of the simplified model of this paper was validated with horizontal and vertical well of FVCHRs with triple-porosity dual-permeability behavior, and results reach a good agreement. Type curve, according to pressure derivative curve characteristic, can be divided into eight flow regimes, which includes wellbore storage, skin reflect, early vertical radial flow, top and bottom boundary reflection, linear flow, inter-porosity flow from vug to fracture, inter-porosity from matrix to vug and fracture and pseudo-radial flow regimes. The influence of some vital parameters (inclination angle, inter-porosity flow coefficient, well length and permeability ratio etc.) on dimensionless pressure and its derivative curves are discussed in detail. The presented model can be used to understand pressure transient response characteristic of inclined wells in FVCHRs with dual-permeability behavior.

Introduction

Since the 1960s, numerous scholars mainly studied the flow characteristic of the fluid in naturally fractured reservoirs by the method of the pressure transient analysis (Warren and Root, 1963; Pruess and Narasimhan, 1985; Kuchuk and Biryukov, 2014; Kui et al., 2017; Tan et al., 2018; Yin, 2018; Song et al., 2019; Pang et al., 2020; Tontiwachwuthikul et al., 2020; Wang et al., 2020). Since a large number of the fractured-vuggy carbonate hydrocarbon reservoirs (FVCHRs) were found around the world, most of scholars have focused on studies of FVCHRs to understand flow characteristic of the fluid in the FVCHRs (Abadassah and Ershaghi, 1986; Liu et al., 2003; Sun et al., 2018). Compared with naturally fractured reservoirs, the complex structures of the FVCHRs posed significant research challenges towards scholars and scientists about pressure transient response of the FVCHRs. Therefore, the mathematic model descripting flow process of the fluid in FVCHRs is established to evaluate the properties and pressure transient characteristics of this type reservoirs.

During the last 20  years, pressure transient analysis is used widely to obtain some vital parameters of reservoirs. With the help of these parameters, reservoir engineers can evaluate oil and gas reservoirs effectively and develop a reasonable development plan for the FVCHRs. The source function method plays a vital role in the process of well test model solution. Gringarten and Ramey (1973) first put forward the real-domain point source analytical solution under the conditions of the different outer boundary. On the basis of the real-domain point source solution, Ozkan and Raghavan (1991; 2004) applied Green’s function theory to expand the point source to the Laplace domain and the corresponding calculation method is given.

Most of the world’s giant oil ans gas produce from naturally fractured and vuggy carbonate reserviors with complex pore system. Vug porosity is common pore system in FVCHRs, and the size of the vug have an imoprtant influence on pressure and production performance. The FVCHRs can be subdivided into triple-continuum medium model and discontinuous medium according to the size of the vug. For the research of the pressure and production of triple-continuum medium model, Abadassah and Ershaghi (1986) first presented triple-porosity single-permeability well-test mathematical model for vertical well, and semi-log and log-log curves were drawn to analysis flow characteristic of the triple-continuum medium. Camacho-Velázquez et al. (2005) established and solved mathematical model of triple-porosity dual-permeability of the vertical well, which assumed that fluid of the natural fracture and vuggy porosity flow into wellbore directly, and the characteristic of pressure transient and decline curve are analyzed. Wu et al. (2004; 2011) solved a triple-continuum pressure transient model of the multi-phase fluid by employing analytical and numerical method. Corbett et al. (2010) studied on fractured carbonate rocks by adopting the numerical well test, and summarized that numerical well testing has its limitations. Guo et al. (2012) established a triple-porosity single-permeability flow mathematical model of the horizontal well in a triple-continuum FVCHRs under different boundary condition, and standard pressure and its derivative log–log curves were drawn. Nie et al. (2012) presented pressure response mathematical model of the horizontal well in coal seam. Pressure and its derivative response curves are drawn in log-log plot, and influence of some vital parameters on type curves were analyzed. Jia et al. (2013) studied on flow issues of a porous–vuggy carbonate reservoir without fracture system. It is assumed that vug system is spherical and pressure distribute symmetrically. Corresponding mathematical model was established and solved by employing Laplace transform. Wang et al. (2018) developed a pressure transient mathematical model of acid fracturing wells in FVCHRs. Therefore, a composite system is considered, the inner region is made up with finite number of artificial fractures and wormholes, and the outer region is made up triple-porosity medium. Wang et al. (2014) derived well test mathematical model of hydraulic fractured well in FVCHRs by point sink integral method. This model assumed that fluid of matrix and natural fractures system flow into wormholes directly. The existence of large-scale cavities will have a noticeable influence on the transient pressure behavior. For the research of the pressure and production of discontinuous medium model, Gao et al. (2016) and Xing et al. (2018) presented a flow mathematical model, it is assumed that well is drilled through large-scale cavities and fluid flow into a filled cavity. Wellbore storage and skin factor were considered in their model. Du et al. (2019) presented an analytical well test model of fractured vuggy carbonate reservoirs. Pressure difference between the wellbore and large vug was considered in their model.

Because the reservoir geological condition becomes more and more complex, the technology of inclined wells has recently been adopted to solve the problem of onshore oil and gas reservoirs extraction. To estimate pressure transient behavior of inclined well, many scholars established pressure response mathematical model of inclined well. Cinco-Ley et al. (1975; 2007) presented an analytical solution of inclined well in real domain and shown the existence of the pseudo-radial flow period at large time values. Yildiz and Ozkan (Yildiz and Ozkan, 1999; Ozkan et al., 1999; Ozkan and Raghavan, 2000) studied transient pressure behavior and rate distributions of slanted well and perforated wells. Corresponding calculation method of slanted and perforated wells was shown to improve the computational efficiency and increase the accuracy. Feng and Liu (2015) established a well test mathematical model of inclined well with an impermeable fault by employing point source theory. Dong et al. (2018) presented a skin factor model of partially penetrated directionally-drilled wells in anisotropic reservoirs to study the stimulation effect of slanted.

This paper summarizes the published literature about well test model of the FVCHRs and slanted well, absorbs their highlights and studies on pressure transient behavior of the inclined well in FVCHRs. Therefore, the objective of this paper is to establish and solve well test mathematical model of inclined well in FVCHRs with triple-porosity dual-permeability behavior by employing Laplace transform and finite cosine Fourier transform. The fluid of natural fractures and vug system can flow into wellbore directly. Inter-porosity flows from matrix system to both natural fracture and vug system, and inter-porosity flow from vug system to a natural fracture system were also considered in this paper. Pressure and its derivative type curves under condition of constant production by adopting the Stehfest numerical inversion are drawn. Influence of some vital parameters on type curves are discussed in detail.

Physical Modeling and Assumptions

The fractured-vuggy carbonate reservoirs in Tahe oilfield of China is an example, reservoir is composed of vug, natural fracture and matrix systems (Figure 1). The fluid of the fractures and vug system flow into the inclined wellbore directly. There is an inclined well in top and bottom closed boundary in FVCHRs. The length of the well is denoted by hw, inclination angle of the well is denoted by θw, and height of reservoirs is denoted by h, the coordinates of point source in the reservoir is M(r,0,z), the coordinates of any point on the well is M(r,0,z) (Figure 2). The end view of an inclined well opened practically in FVCHRs is shown in Figure 3. The fluid of the fracture and vug systems flow into the inclined wellbore directly, Inter-porosity flows from matrix system to both natural fracture and vug systems, and inter-porosity flow from vug system to natural fracture system were considered. This process is named triple-porosity dual-permeable flow behavior. In order to better establish mathematical model of an inclined well in FVCHRs, the other assumptions are as follows:

(1) The production of the inclined well is denoted by qsc in FVCHRs of closed the top and bottom boundaries.

(2) The rock and fluid are slightly compressible, and the fluid flowing in the fractures and reservoirs is isothermal seepage.

(3) The inclined well is assumed to be penetrated fully or partially, and formation thickness and well length meet the relationships: |(zDzwD)cosy΄|hwD,ψ΄=tan1(kfv/kfhtanθw).

(4) The fluid in the natural fracture and vug system can flow into the inclined well and inter-porosity flow from vug system to a natural fracture system were considered. The fluid in the matrix system flow into inclined wellbore through vug and fracture system (see Figure 3A).

(5) The permeability of the matrix system is km, the permeability of the natural fracture is kf, the permeability of the vug is kf, and the permeability anisotropy of the horizontal and vertical directions is considered in FVCHRs.

(6) The original formation pressure of the matrix, fracture and vug system is equal to initial pressure pe.

FIGURE 1
www.frontiersin.org

FIGURE 1. Fractured-vuggy carbonate reservoirs scheme (Modified after Nie et al., 2012).

FIGURE 2
www.frontiersin.org

FIGURE 2. The scheme of an inclined well (Modified after Dong et al., 2018).

FIGURE 3
www.frontiersin.org

FIGURE 3. Flow scheme of FVCHRs with triple-porosity dual-permeability behavior.

Mathematical Modeling

Based on the scheme of an inclined well as shown in Figure 2, the mathematical model of the inclined well is established in FVCHRs. In order to help reader to understand coordinate system as shown in Figure 2, three different views of an inclined well in top and bottom closed boundary are given in Figure 4.

FIGURE 4
www.frontiersin.org

FIGURE 4. Three different view in Figure 2 (Modified after Nie et al., 2012). (A) Comparison with Guo et al.’ result. (B) Comparison with Camacho-Velázquez et al. result.

According to the state equation, the motion equation and the continuity equation, the differential equation of matrix, vug and natural fracture systems is established in FVCHRs (Camacho-Velázquez et al., 2005). Inclined well could be regarded as a line source in 3D space. Therefore, to obtain solutions for inclined well, we must establish a 3D point source model. Then we can get a surface source solution through point source integral method.

(1) The two-dimensional seepage differential equation of natural fracture system can be described as:

1rkfhμr(rpfr)+kfvμz(pfz)+αmkmhμ(pmpf)+αVkVhμ(pVpf)=(ϕCt)fpft(1)

(2) The two -dimensional seepage differential equation of vug system can be described as:

1rkVhμr(rpVr)+kVvμz(pVz)+αmkmhμ(pmpV)αVkVhμ(pfpV)=(ϕCt)VpVt(2)

(3) The seepage differential equation of matrix system can be described as:

αmkmhμ(pmpV)αmkmhμ(pmpf)=(ϕCt)mpmt(3)

(4) Initial condition can be described as:

pm(r,0)=pf(r,0)=pV(r,0)=pe(4)

(5) Inner condition boundary can be described as:

limε0zε/2z+ε/2[limδ02πrμ(kfhpfr+kVhpVr)r=ε]dz=q˜(t)(5)

(6) Outer condition boundary can be described as:

pf(r,z=0,t)z=pm(r,z=0,t)z=pV(r,z=0,t)z=0(6)
pm(r,z=h,t)z=pf(r,z=h,t)z=pV(r,z=h,t)z=0(7)
pf(r,z,t)=pm(r,z,t)=pV(r,z,t)=0(8)

where: pV is pressure of vug system, MPa; pf is pressure of natural fracture system, MPa; pm is pressure of matrix system, MPa; pe is initial reservoirs pressure, MPa; kfh is horizontal permeability of natural fracture system, mD; kfv is vertical permeability of natural fracture system, mD; kmh is horizontal permeability of matrix system, mD; kVh is horizontal permeability of vuggy system, mD; Ctf is total compressibility factor of natural fracture system, MPa−1; Cmf is total compressibility factor of matrix system, MPa−1; CVf is total compressibility factor of vuggy system, MPa−1; φf is porosity of natural fracture system, dimensionless; φm is porosity of matrix system, dimensionless; φV is porosity of vuggy system, dimensionless; t is production time, s; μ is fluid viscosity, cp; q˜ is strength of continuous point source of reservoirs, m2/s; α is shape factor; x,y,z is Descartes coordinates, m; ε,δ is infinitesimal distance, m; qsc is production rate, m3/s; rw is wellbore radial, m.

For convenience of derivation, the relevant dimensionless variables are defined in Table 1: Laplace transformation with respect to tD is then adopted to invert the above seepage model into Laplace domain (Van Everdingen and Hurst, 1949).

L[pjD(rD,tD)]=p¯jD(rD,s)=0pjD(rD,tD)estDdtD,(j=m,  f,  V)(9)

where: p¯jD is dimensionless pressure in Laplace domain;

TABLE 1
www.frontiersin.org

TABLE 1. Dimensionless variables definition.

According to the dimensionless definition, the dimensionless forms of the differential equation and boundary conditions can be given by:

{κ1rDrD(rDp¯fDrD)+zD(p¯fDzD)+λmf(p¯mDp¯fD)+λVf(p¯VDp¯fD)=ωfsp¯fD(1κ)1rDrD(rDp¯VDrD)+kDzD(p¯VDzD)+λmV(p¯mDp¯VD)λVf(p¯VDp¯fD)=ωVsp¯VDλmV(p¯mDp¯VD)λmf(p¯mDp¯fD)=ωmsp¯mDp¯mD(rD,0)=p¯fD(rD,0)=p¯VD(rD,0)=0limεD0zwDεD/2zwD+εD/2[limδD0(κrDp¯fDrD+(1κ)rDp¯VDrD)rD=εD]dzD=q˜¯Dp¯fD(rD,zD=0,tD)zD=p¯mD(rD,zD=0,tD)zD=p¯VD(rD,zD=0,tD)zD=0p¯fD(rD,zD=hD,tD)zD=p¯mD(rD,zD=hD,tD)zD=p¯VD(rD,zD=hD,tD)zD=0p¯fD(rD,zD,tD)=p¯mD(rD,zD,tD)=p¯VD(rD,zD,tD)=0(10)

where: kD is permeability ratio between vug and fracture, dimensionless; κ is permeability ratio, dimensionless; λmf is inter-porosity flow coefficient from matrix to fracture, dimensionless; λmv is inter-porosity flow coefficient from matrix to vuggy, dimensionless; λvf is inter-porosity flow coefficient from vuggy to fracture, dimensionless; ωf is fracture storativity ratio, dimensionless; ωm is matrix storativity ratio, dimensionless; ωv is vuggy storativity ratio, dimensionless; s is Laplace variables; subscripts D is dimensionless variables; superscripts – is variables in the Laplace domain;

By employing the finite cosine Fourier transform with respect to zD in Eq. 10, and expressing dimensionless pressures of natural fracture and vug systems by dimensionless matrix system pressure in Laplace domain after the finite cosine Fourier transform, Eq. 10 can be described as:

{κ1rDrD(rDp¯fDrD)+m1np¯fD+λp¯VD=0(1κ)1rDrD(rDp¯VDrD)+λp¯fD+m2np¯VD=0limεD0(κrDp¯fDrD+(1κ)rDp¯VDrD)rD=εD=hDcos(nπzDhD)q˜¯Dp¯fD(rD,n,s)=p¯VD(rD,n,s)=p¯mD(rD,n,s)=0(11)

where: m1n=λmf2λmf+λmV+sωmλmfλVf(sωf+n2π2hD2); λ=λmfλmVλmf+λmV+sωf+λVf; subscripted variable(n=0,1,2,,n).

By solving Eq. 10, the general point source solution for the top and bottom closed boundary, constant production and infinite lateral boundary FVCHRs with triple-porosity dual-permeability behavior can be obtained through Laplace transform and Fourier transform methods (Liu et al., 2019; Wang et al., 2019; Xue et al., 2019):

p¯D=q˜¯D[i=12Ci0K0(σi0rD)+2n=1i=12CinK0(σinrD)cos(nπzDhD)cos(nπzDhD)](12)

where: {C1n=(1a2)[κ+(1κ)a1](1a2)(1a1)[κ+(1κ)a2]C2n=(1a1)[κ+(1κ)a2](1a1)(1a2)[κ+(1κ)a1], (n=0,1,2,...,n); rD=(xDxD)2+(yDyD)2

It is useful to express the Eq. 12 in polar coordinates, let the polar coordinates of the points (xD,yD) and (xD,yD) be given by (rD,θ) and (rD,θ) respectively according to Raghavan and Ozkan (2017). We can do integration to the point source solution along inclined well in polar coordinates system as shown in Figure 4C. Therefore, the solution of the inclined well is obtain as following.

p¯wDincline=1shwDsinψ(hwD/2)sinψ+(hwD/2)sinψ[i=12Ci0K0(σi0RD)+2n=1i=12CinK0(σinRD)cos(nπzDhD)cos(nπzDhD)]drD(13)

where: K0(x) is the second kind modified Bessel function; RD=(rDcosθrD)2+(rDsinθ)2; zD=zwD+rDcotψ.

Although Gringarten et al. (1974) pointed out that the wellbore pressure of horizontal well and infinite conductivity hydraulic fractured well can be calculate by an equivalent-pressure point xD=0.732, there is difference between horizontal well and inclined well. Fortunately, Cinco-Ley et al. (1975) given the equivalent-pressure point of inclined well in the polar coordinates. As following:

rD=1+(0.3hwDsinψ)2;θ=cos1(0.3hwDsinψrD);zD=zwD±0.3hwDcosψ;

In order to enhance calculation speed of the inclined well, Ozkan and Raghavan (2000) provides an efficient algorithm to compute transient pressure responses of inclined wells in Laplace domain. The principle of convolution can be used to obtain solution with the effects of wellbore storage and skin factor, which could be given by

p¯wD=sp¯wDincline+Ss+CDs2(sp¯wDincline+S)(14)

where: pwD is dimensionless wellbore pressure with wellbore storage and skin; pwDincline is dimensionless wellbore pressure without wellbore storage and skin; CD is dimensionless wellbore storage; S is skin factor.

The solution of inclined well in real domain for FVCHRs with triple-porosity dual-permeability behavior is obtained by employing Stehfest (1970) numerical inversion algorithm. Dimensionless pressure and its derivative type curves are drawn analyzed.

Model Verification

Compared with vertical wells and horizontal wells, wellbore pressure solution of the inclined well is more universal, so our work expands analytical solution of vertical well and horizontal well in FVCHRs with triple-porosity dual-permeability behavior. However, in order to verify the correctness of our model and solution, the simplified model of this paper is compared with analytical solution of horizontal well for FVCHRs with triple-porosity dual-permeability behavior (Guo et al., 2012) and vertical well for FVCHRs with triple-porosity dual-permeability behavior (Camacho-Velázquez et al., 2005).

Guo et al. (2012) established the mathematical model of the horizontal well under condition of different outer boundary in FVCHRs with triple-porosity dual-permeability behavior. What is different with our model is that they assumed that fluid of the naturally fractured and the matrix system are directly flow into the wellbore, inter-porosity flow from matrix to fracture and inter-porosity flow from vug to matrix and fracture is considered. In order to make verification result reliable, if parameters λVf,λVm,λmf of the Guo’ model are replaced by parameters λmf,λmV,λVf of this paper correspondingly and corresponding value of these parameters is same. On the basis of replacement, we can prove that our model is same with Guo’ model. If θw is 90° in our paper, and hw=1,000 of our paper, inclined well of this paper will become horizontal well, which is similar to the Guo et al.’ model with horizontal well length is 500. Other common parameters are same and their values are: h=100, rw=0.1, khv=1, kD=1, κ=0.7, ωf=0.001, ωV=0.1, λmf=1012, λmV=1010, λVf=108. As shown in Figure 5A, result reaches a good agreement.

FIGURE 5
www.frontiersin.org

FIGURE 5. Comparison with Guo et al. (2012) and Camacho-Velázquez et al. (2005) models (h=100,hw=100.4,rw=0.1, khv=1,kD=1,κ=0.7,ωf=0.001, ωV=0.1,λmf=1012,λmV=1010,λVf=108)

Camacho-Velázquez et al. (2005) established the mathematical model of the vertical well in FVCHRs with triple-porosity dual-permeability behavior. It is assumed that fluid of the fracture and the vug system are directly flow into the wellbore. Inter-porosity flow from vug to fracture and inter-porosity flow from matrix to vug and fracture is considered, which is same with our assumption. If inclination angle of wellbore is less than 5° and well is penetrated fully, we can think that our model is equal to vertical well. Other common parameters are same and their values are: h=100, rw=0.1, khv=1,kD=1,κ=0.7,ωf=0.001,ωV=0.1,λmf=1012,λmV=1010,λVf=108. As shown in Figure 5A, result reaches a good agreement.

By comparing with horizontal and vertical well in FVCHRs with triple-porosity dual-permeability behavior, the results reached a good agreement, which verifies the correctness of our model and solution.

Type Curve and Sensitivity of Parameters

Type Curve

Figure 6 shows type curves of dimensionless pressure and its derivative of inclined well with dual-permeability behavior. According to the characteristic of the dimensionless pressure derivative curve, the type curve of inclined well can be divided into eight flow regimes for FVCHRs with triple-porosity dual-permeability behavior, which mainly includes well storage regime, skin effect regime, early vertical radial flow regime (Figure 7A), top and bottom boundary reflection (Figure 7B), linear flow regime (Figure 7C), inter-porosity flow and pseudo-radial flow regime (Figure 7D).

FIGURE 6
www.frontiersin.org

FIGURE 6. Pressure transient type curve of inclined well in FVCHRs with triple-porosity dual-permeability behavior (CD=0.1,S=0.1,h=100,hw=h/cos(ψ), rw=0.1,θw=75o,khv=1,kD=1,κ=0.7,ωf=0.01, ωV=0.1,λmf=1012,λmV=1010,λVf=108)

FIGURE 7
www.frontiersin.org

FIGURE 7. Main regime flow schemes for inclined well.

Regime①: Wellbore storage regime. During the pure wellbore storage, pressure and its derivative curves coincide and is a unit-slope straight line.

Regime②: Skin effect regime. Due to consideration of skin, pressure derivative curve looks like a “hump.”

Regime③: Early vertical radial flow regime. Pressure derivative curve is a constant value horizontal line in log-log plot (Figure 7A).

Regime④: Top and bottom boundary reflection regime. Pressure derivative curve take on a small step in log-log plot. However, this stage may disappear because of change of the length of wellbore and inclination angle of wellbore (Figure 7B).

Regime⑤: Linear flow regime. Pressure derivative curve is 1/2-slope straight line in log-log plot. Shorter well length lead to disappearance of 1/2-slope straight line (Figure 7C).

Regime⑥: Inter-porosity flow from vug to fracture. Because vug permeability is better than matrix permeability, inter-porosity flow from vug to fracture appear firstly. The curve of pressure derivative is V-shaped.

Regime⑦: Inter-porosity flow from matrix to vug and fracture. The curve of pressure derivative is V-shaped, which is controlled by parameters ωm and ωv.

Regime⑧: Pseudo-radial flow regime. Inter-porosity flows have already finished and the curve of pressure derivative is a 0.5-value horizontal line in log-log plot (Figure 7D).

Sensitivity of Parameters

Figure 8 shows type curve affected by length of inclined well in FVCHRs with triple-porosity dual-permeability behavior. Because short inclined well length can decrease contract area between reservoirs and wellbore, the short inclined well length can increase pressure loss that fluid flows into wellbore. Therefore, short inclined well length can lead to high position of pressure curve in the whole flow regime in log-log plot and make top and bottom boundary reflection regime disappear and spherical flow characteristic more obvious in pressure derivative curve.

FIGURE 8
www.frontiersin.org

FIGURE 8. Type curve affected by length of inclined well (CD=0.1,S=0.1,h=100,rw=0.1,θw=60o,khv=1,kD=1,κ=0.7,ωf=0.01,ωV=0.1,  λmf=1012,  λmV=1010,  λVf=108)

Figure 9 shows type curve affected by inclination angle in FVCHRs with triple-porosity dual-permeability behavior. It is assumed that length of inclined well keep constant when inclination angle changes. The large inclination angle change can only lead to disappearance of the top and bottom boundary reflection regime of pressure derivative curve without consideration of anisotropy (Figure 9A). However, if the inclined well is penetrated fully, inclination angle has a more effect on pressure and its derivative curves than inclined well penetrated partly. With the increasing of inclination angle length of inclined well become length, which leads to low position of pressure curve in whole regime and low position of pressure derivative curve during early vertical radial flow regime (Figure 9B).

FIGURE 9
www.frontiersin.org

FIGURE 9. Type curve affected by inclination angle of well (CD=0.1,S=0.1,h=100,hw=h/cos(ψ),rw=0.1,khv=1,kD=1,κ=0.7,ωf=0.01,ωv=0.1,  λmf=1012,  λmV=1010,  λVf=108)

Figure 10 shows type curve affected by vug storativity ratio in FVCHRs with triple-porosity dual-permeability behavior. According the flow model shown Figure 3A, fluid of the vug system do not flow into fracture, but also flow into wellbore directly. Therefore, as shown in Figure 10, small vug storativity ratio can make the second V-shaped of the pressure derivative curve become more and more deep while the first V-shaped become more and more shallow.

FIGURE 10
www.frontiersin.org

FIGURE 10. Type curve affected by parameter ωv(ωf=0.01).

Figure 11 shows type curve affected by fracture storativity ratio in FVCHRs with triple-porosity dual-permeability behavior. According the flow model shown Figure 3A, fluid of fracture system only flow into wellbore directly. Therefore, as is shown in Figure 11, fracture storativity ratio has an influence on the first V-shaped. The first V-shaped become more and more deeply with the decrease of fracture storativity ratio. The position of pressure curve become higher during vertical radial flow regime and top and bottom boundary reflection regime.

FIGURE 11
www.frontiersin.org

FIGURE 11. Type curve affected by parameter ωf(ωv=0.01). (CD=0.1,S=0.1,h=100,hw=200, rw=0.1,θw=60o,kD=1,κ=0.7, λmf=1012,λmV=1010,λVf=108)

Figure 12 shows type curve affected by inter-porosity flow coefficient λVf in FVCHRs with triple-porosity dual-permeability behavior. According the flow model shown Figure 3A, fluid of fracture system only flow into wellbore directly. Therefore, as shown in Figure 12, inter-porosity flow coefficient λVfλVf has an influence on the first V-shaped. As λVf increase, starting time of inter-porosity flow from vug system to fracture system become earlier, which makes the first V-shaped of pressure derivative move towards left in log-log plot.

FIGURE 12
www.frontiersin.org

FIGURE 12. Type curve affected by parameter λVf(λmV=1010)

Figure 13 shows type curve affected by inter-porosity flow coefficient λmV in FVCHRs with triple-porosity dual-permeability behavior. As shown in Figure 13, inter-porosity flow coefficient λmv has an influence on the second V-shaped. As λmV increase, starting time of inter-porosity flow from matrix system to vug system become earlier, which makes the second V-shaped of pressure derivative move towards left in log-log plot. Of course, with the flow coefficient λmV keep increasing, first V-shaped may be covered and first V-shaped become more deeply.

FIGURE 13
www.frontiersin.org

FIGURE 13. Type curve affected by parameter λmV(λVf=108)(CD=0.1,S=0.1,h=100,hw=200,rw=0.1,θw=65o,khv=1,kD=1,κ=0.7,ωf=0.01,ωV=0.1,λmf=1012)

Figure 14 shows type curve affected by inter-porosity flow coefficient λmf in FVCHRs with triple-porosity dual-permeability behavior. As shown in Figure 14, inter-porosity flow coefficient λmf has an influence on the second V-shaped. As λmf increase, starting time of inter-porosity flow from matrix system to fracture system become earlier, which makes the second V-shaped of pressure derivative move towards left in log-log plot. Of course, because some of fluid of fracture system come from supplement of vug system, but also others come from supplement of matrix system. Therefore, when inter-porosity flow coefficient λmf<λvf, the change of inter-porosity flow coefficient λmf has an influence on pressure derivative curve slightly.

FIGURE 14
www.frontiersin.org

FIGURE 14. Type curve affected by parameter λmf(κ=0.7)

Figure 15 shows type curve affected by permeability ratio κ in FVCHRs with triple-porosity dual-permeability behavior. The large permeability ratio indicates that permeability of fracture system is higher than vug system, and vug storativity capacity is lower than fracture storativity capacity. Therefore, the first V-shaped of pressure derivative curve becomes more and more deeply with the increase of permeability ratio.

FIGURE 15
www.frontiersin.org

FIGURE 15. Type curve affected by parameter κ(λmf=1012) (CD=0.1,S=0.1,h=100,hw=200,rw=0.1, θw=60o,khv=1,kD=1,ωf=0.01, ωV=0.1,  λmV=1012,λVf=108)

Figure 16 shows type curve affected by permeability ratio khv between horizontal permeability and vertical permeability in FVCHRs with triple-porosity dual-permeability behavior. The large khv indicates that permeability of horizontal permeability is higher than vertical permeability. Because inclination angle of wellbore is considered in vertical direction, difference between vertical and horizontal flow is shown obviously during early vertical radial flow regime. The greater the value of khv, the greater the loss of wellbore flow pressure in the vertical direction. Therefore, large khv will lead to higher position of pressure and its derivative curve in early flow regime.

FIGURE 16
www.frontiersin.org

FIGURE 16. Type curve affected by parameter khv(CD=0.1,S=0.1,h=100,hw=200,rw=0.1,θw=60o,kD=1,κ=0.7,ωf=0.01,ωV=0.1,  λmV=1012,λVf=108)

Conclusion

It is assumed that fluid of vug system and fracture system flows into incline wellbore directly, and inter-porosity flow from vug system to fracture system and inter-porosity flow from matrix system to vug and fracture system is considered. Pressure transient mathematical model of inclined well is established and solved by employing Laplace transform and finite cosine Fourier transform for FVCHRs with triple-porosity dual-permeability behavior. The following conclusions can be drawn:

(1) A new analytical model was established to analysis pressure transient behavior of inclined well in FVCHRs with triple-porosity dual-permeability behavior.

(2) The model was validated by comparing our simplified model with model of horizontal well and vertical well of FVCHRs with triple-porosity dual-permeability behavior. Result reach a good agreement.

(3) Eight flow regimes, according to the characteristic of pressure derivative curve, can be recognized, and flow characteristic of each flow regime was described in detail.

(4) The effect of some vital parameters (such as well length, well inclination angle and permeability ratio, etc.) on pressure and its derivative curves are analyzed.

Data Availability Statement

The raw data supporting the conclusions of this article will be made available by the authors, without undue reservation.

Author Contributions

XL and ZM designed the research; SQ and XT built the model; SQ and YX analyzed and simulated experimental data; ZM, TD and SQ were in charge of revision and writing.

Funding

This article was supported by the National Natural Science Foundation of China (51704246) and Sichuan Science and Technology Project (18SYXHZ0090).

Conflict of Interest

The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.

References

Abadassah, D., and Ershaghi, I. (1986). Triple-porosity systems for representing naturally fractured reservoirs. SPE Form. Eval. 1 (2), 113–127. doi:10.2118/13409-pa

CrossRef Full Text | Google Scholar

Camacho-Velázquez, R., Vásquez-Cruz, M, Castrejón-Aivar, R., and Arana-Ortiz, V. (2005). Pressure-transient and decline-curve behavior in naturally fractured vuggy carbonate reservoirs. SPE Reservoir Eval. Eng. 8 (2), 95–111. doi:10.2118/77689-pa

CrossRef Full Text | Google Scholar

Cinco-Ley, H., Miller, F. G., and Ramey, H. J. (1975). Unsteady-State pressure distribution created by a directionally drilled well. J. Petrol. Technol. 27, 1392–1400. doi:10.2118/5131-PA

CrossRef Full Text | Google Scholar

Cinco-Ley, H., Ramey, H. J., and Miller, F. (2007). “Pseudo-skin factors for partially-penetrating directionally-drilled wells,” in Fall meeting of the society of petroleum engineers of AIME, Dallas, TX, 28 September–1 October. doi:10.2523/5589-ms

CrossRef Full Text | Google Scholar

Corbett, P. W.M., Geiger, S., Borges, L., Garayev, M., Gonzalez, J. G., and Valdez, C. (2010). “Limitations in numerical well rest modelling of fractured carbonate rocks.” in 72nd European association of geoscientists and engineers conference and exhibition 2010: a new Spring for geoscience, Barcelona, Spain, 14–17 June, 1135–1145. doi:10.2523/130252-ms

CrossRef Full Text | Google Scholar

Dong, W., Wang, X., and Wang, J. (2018). A new skin factor model for partially penetrated directionally-drilled wells in anisotropic reservoirs. J. Petrol. Sci. Eng. 161, 334–348. doi:10.1016/j.petrol.2017.11.062

CrossRef Full Text | Google Scholar

Du, X., Lu, Z, Li, D, Xu, Y, Li, P, and Lu, D (2019). A novel analytical well test model for fractured vuggy carbonate reservoirs considering the coupling between oil flow and wave propagation. J. Petrol. Sci. Eng. 173, 447–461. doi:10.1016/j.petrol.2018.09.077

CrossRef Full Text | Google Scholar

Feng, G. Q., and Liu, Q.-G. (2015). Pressure transient behavior of a slanted well with an impermeable fault. J. Hydrodyn. 26 (6), 980–985. doi:10.1016/S1001-6058(14)60108-6

CrossRef Full Text | Google Scholar

Gao, B., Huang, X., Yao, J., Lv, X., and Wu, Y. (2016). Pressure transient analysis of a well penetrating a filled cavity in naturally fractured carbonate reservoirs. J. Petrol. Sci. Eng. 145, 392–403. doi:10.1016/j.petrol.2016.05.037

CrossRef Full Text | Google Scholar

Gringarten, A. C., and Ramey, H. J. (1973). Use of source and green’S functions in solving unsteady-flow problems in reservoirs. Soc. Pet. Eng. AIME J. 13 (5), 285–296. doi:10.2118/3818-PA

CrossRef Full Text | Google Scholar

Gringarten, A. C., Ramey, H. J., and Raghavan, R. (1974). Unsteady-state pressure distributions created by a well with a single infinite-conductivity vertical fracture. Soc. Pet. Eng. AIME J. 14 (4), 347–360. doi:10.2118/4051-pa

CrossRef Full Text | Google Scholar

Guo, J. C., Nie, R. S., and Jia, Y. L. (2012). Dual permeability flow behavior for modeling horizontal well production in fractured-vuggy carbonate reservoirs. J. Hydrol. 464–465, 281–293. doi:10.1016/j.jhydrol.2012.07.021

CrossRef Full Text | Google Scholar

Jia, Y. L., Fan, X. Y., Nie, R. S., Huang, Q. H., and Jia, Y. L. (2013). Flow modeling of well test analysis for porous-vuggy carbonate reservoirs. Transp. Porous Media 97 (2), 253–279. doi:10.1007/s11242-012-0121-y

CrossRef Full Text | Google Scholar

Kuchuk, F., and Biryukov, D. (2014). Pressure-transient behavior of continuously and discretely fractured reservoirs. SPE Reservoir Eval. Eng. 17 (1), 82–97. doi:10.2118/158096-PA

CrossRef Full Text | Google Scholar

Kui, M. Q., Tan, X. H., Li, X. P., and Cai, J. (2017). Analysis on unsteady flow for power-law fluids in dual fractal media. J. Porous Media 20 (12), 1071–1086. doi:10.1615/JPorMedia.v20.i12.20

CrossRef Full Text | Google Scholar

Liu, J., Bodvarsson, G. S., and Wu, Y. S. (2003). Analysis of flow behavior in fractured lithophysal reservoirs. J. Contam. Hydrol. 62–63, 189–211. doi:10.1016/S0169-7722(02)00169-9

CrossRef Full Text | Google Scholar

Liu, Q., Xu, Y., Peng, X., Liu, Y., and Qi, S. (2019). Pressure transient analysis for multi-wing fractured wells in dual-permeability hydrocarbon reservoirs. J. Petrol. Sci. Eng. 180, 278–288. doi:10.1016/j.petrol.2019.05.013

CrossRef Full Text | Google Scholar

Nie, R. S., Meng, Y. F., Guo, J. C., and Jia, Y. L. (2012). Modeling transient flow behavior of a horizontal well in a coal seam. Int. J. Coal Geol. 92, 54–68. doi:10.1016/j.coal.2011.12.005

CrossRef Full Text | Google Scholar

Ozkan, E., and Raghavan, R. (1991). New solutions for well-test-analysis problems. Part 2. Computational considerations and applications. SPE Form. Eval. 6 (3), 369–378. doi:10.2118/18616-PA

CrossRef Full Text | Google Scholar

Ozkan, E., and Raghavan, R. (2000). A computationally efficient, transient-pressure solution for inclined wells. SPE Reservoir Eval. Eng. 3 (5), 414–425. doi:10.2118/66206-PA

CrossRef Full Text | Google Scholar

Ozkan, E., and Raghavan, R. (2004). New solutions for well-test-analysis problems: Part 1-Analytical Considerations (includes associated papers 28666 and 29213). SPE Form. Evaluation. 6 (03), 359–368. doi:10.2118/18615-PA

CrossRef Full Text | Google Scholar

Ozkan, E., Yildiz, T., and Raghavan, R. (1999). “Pressure-transient analysis of perforated slant and horizontal wells,” in SPE Annual Technical Conference and Exhibition, Houston, TX, October 3–6, 1999.

CrossRef Full Text | Google Scholar

Pang, Y., Hu, X, Wang, S, Chen, S., Soliman, M. Y., and Deng, H. (2020). Characterization of adsorption isotherm and density profile in cylindrical nanopores: modeling and measurement. Chem. Eng. J. 396, 125212. doi:10.1016/j.cej.2020.125212

CrossRef Full Text | Google Scholar

Pruess, K., and Narasimhan, T. N. (1985). Practical method for modeling fluid and heat flow in fractured porous media. Soc. Petrol. Eng. J. 25 (1), 14–26. doi:10.2118/10509-pa

CrossRef Full Text | Google Scholar

Raghavan, R, and Ozkan, E. (2017). A method for computing unsteady flows in porous media. Oxon, UK: Routledge. doi:10.1201/9781315136387

CrossRef Full Text | Google Scholar

Song, F., Bo, L., Zhang, S., and Sun, Y. (2019). Nonlinear flow in low permeability reservoirs: modelling and experimental verification. Adv. Geo-Energy Res. 3, 76–81. doi:10.26804/ager.2019.01.06

CrossRef Full Text | Google Scholar

Stehfest, H. (1970). Algorithm 368: numerical inversion of Laplace transforms [D5]. Commun. ACM 13 (1), 47–49.

CrossRef Full Text | Google Scholar

Sun, Q., Zhang, N., Mohamed, F., and Wang, Y. (2018). Structural regeneration of fracture-vug network in naturally fractured vuggy reservoirs. J. Petrol. Sci. Eng. 165, 28–41. doi:10.1016/j.petrol.2017.11.030

CrossRef Full Text | Google Scholar

Tan, X. H., Jiang, L., Li, X. P., Zhang, B. J., and Li, X. C. (2018). Flow model of a multi-stage hydraulic fractured horizontal well based on tree-shaped fractal fracture networks. J. Petrol. Sci. Eng. 169, 494–503. doi:10.1016/j.petrol.2018.06.008

CrossRef Full Text | Google Scholar

Tontiwachwuthikul, P., Chan, C. W., Zeng, F., Liang, Z., Sema, T., and Chao, M. (Forthcoming 2020). Recent progress and new developments of applications of artificial intelligence (AI), knowledge-based systems (KBS), and machine learning (ML) in the petroleum industry. Petroleum. doi:10.1016/j.petlm.2020.08.001

CrossRef Full Text | Google Scholar

Van Everdingen, A. F., and Hurst, W. (1949). The application of the Laplace transformation to flow problems in reservoirs. J. Petrol. Technol. 1 (12), 305–324. doi:10.2118/949305-g

CrossRef Full Text | Google Scholar

Wang, L., Cheng, D., Xiang, L., Chen, X., and Xia, Z. (2019). Pressure transient analysis for asymmetrically fractured wells in dual-permeability organic compound reservoir of hydrogen and carbon. Int. J. Hydrog. Energy 44 (11), 5254–5261. doi:10.1016/j.ijhydene.2018.08.082

CrossRef Full Text | Google Scholar

Wang, L., Wang, X., Luo, E., and Wang, J. (2014). Analytical modeling of flow behavior for wormholes in naturally fractured–vuggy porous media. Transp. Porous Media 105 (3), 539–558. doi:10.1007/s11242-014-0383-7

CrossRef Full Text | Google Scholar

Wang, M., Fan, Z., Dong, X., Song, H., Zhao, W., and Xu, G (2018). Analysis of flow behavior for acid fracturing wells in fractured-vuggy carbonate reservoirs. Math. Probl. Eng. 2018 (PT. 3), 6431910. doi:10.1155/2018/6431910

CrossRef Full Text | Google Scholar

Wang, Q., Yang, S., Glover, P. W. J., Lorinczi, P., Qian, K., and Wang, L. (2020). Effect of pore-throat microstructures on formation damage during miscible CO2 flooding of tight sandstone reservoirs. Energy Fuels 34 (4), 4338–4352. doi:10.1021/acs.energyfuels.0c00158.

CrossRef Full Text | Google Scholar

Warren, J.E., and Root, P. J. (1963). The behavior of naturally fractured reservoirs. Soc. Petrol. Eng. J. 3 (3), 245–255. doi:10.2118/426-pa

CrossRef Full Text | Google Scholar

Wu, Y. S., Liu, H. H., and Bodvarsson, G. S. (2004). A triple-continuum approach for modeling flow and transport processes in fractured rock. J. Contam. Hydrol. 73 (1–4), 145–179. doi:10.1016/j.jconhyd.2004.01.002

CrossRef Full Text | Google Scholar

Wu, Y. S., Yuan, D., Kang, Z., and Fakcharoenphol, P. (2011). A multiple-continuum model for simulating single-phase and multiphase flow in naturally fractured vuggy reservoirs. J. Petrol. Sci. Eng. 78 (1), 13–22. doi:10.1016/j.petrol.2011.05.004

CrossRef Full Text | Google Scholar

Xing, C., Yin, H., Liu, K., Li, X, and Fu, J. (2018). Well test analysis for fractured and vuggy carbonate reservoirs of well drilling in large scale cave. Energies 11 (1), 1–15. doi:10.3390/en11010080

CrossRef Full Text | Google Scholar

Xue, L., Chen, X., and Wang, L. (2019). Pressure transient analysis for fluid flow through horizontal fractures in shallow organic compound reservoir of hydrogen and carbon. Int. J. Hydrog. Energy 44 (11), 5245–5253. doi:10.1016/j.ijhydene.2018.08.084

CrossRef Full Text | Google Scholar

Yildiz, T., and Ozkan, E. (1999). Pressure-transient analysis for perforated wells. SPE J. 4. doi:10.2118/56878-PA

CrossRef Full Text | Google Scholar

Yin, C. (2018). Test and analysis on the permeability of induced fractures in shale reservoirs. Nat. Gas. Ind. B 5 (5), 513–522. doi:10.1016/j.ngib.2018.03.006

CrossRef Full Text | Google Scholar

Keywords: inclined well, fractured-vuggy carbonate hydrocarbon reservoirs, triple-porosity dual-permeability, inter-porosity flow, pressure transient

Citation: Qi S-Z, Tan X-H, Li X-P, Meng Z, Xu Y-J and Tang D (2020) Transient Pressure Analysis of Inclined Well in Continuous Triple-Porosity Reservoirs With Dual-Permeability Behavior. Front. Energy Res. 8:601082. doi: 10.3389/fenrg.2020.601082

Received: 31 August 2020; Accepted: 02 October 2020;
Published: 01 December 2020.

Edited by:

Kaiqiang Zhang, Imperial College London, United Kingdom

Reviewed by:

Fanhua Zeng, University of Regina, Canada
Fuping Li, Curtin University, Australia

Copyright © 2020 Qi, Tan, Li, Meng, Xu and Tang. This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) and the copyright owner(s) are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.

*Correspondence: Sheng-Zhi Qi, ctfy742@163.com; Xiao-Hua Tan, xiaohua-tan@163.com

Disclaimer: All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article or claim that may be made by its manufacturer is not guaranteed or endorsed by the publisher.