# Optimal Inflow Performance Relationship Equation for Horizontal and Deviated Wells in Low-Permeability Reservoirs

Authors: Liqiang Wang; Zhengke Li; Mingji Shao; Yinghuai Cui; Wenbo Jing; Wei Zhang; Maoxian Wang
Journal: Geofluids
 (2021)

---

## Abstract

After Vogel proposed a dimensionless inflow performance equation, with the rise of the horizontal well production mode, a large number of inflow performance relationship (IPR) equations have emerged. In the productivity analysis of deviated and horizontal wells, the IPR equation proposed by Cheng is mainly used. However, it is still unclear whether these inflow performance models (such as the Cheng, Klins-Majcher, Bendakhlia-Aziz, and Wiggins-Russell-Jennings types) are suitable for productivity evaluations of horizontal and deviated wells in low-permeability reservoirs. In-depth comparisons and analyses have not been carried out, which hinders improvements in the accuracy of the productivity evaluations of horizontal wells in low-permeability reservoirs. In this study, exploratory work was conducted in two areas. First, the linear flow function relationship used in previous studies was improved. Based on the experimental pressure-volume-temperature results, a power exponential flow function model was established according to different intervals greater or less than the bubble point pressure, which was introduced into the subsequent derivation of the inflow performance equation. Second, given the particularity of low-permeability reservoir percolation, considering that the reservoir is a deformation medium, and because of the existence of a threshold pressure gradient in fluid flow, the relationship between permeability and pressure was changed. The starting pressure gradient was introduced into the subsequent establishment of the inflow performance equation. Based on the above two aspects of this work, the dimensionless IPR of single-phase and oil-gas two-phase horizontal wells in a deformed medium reservoir was established by using the equivalent seepage resistance method and complex potential superposition principle. Furthermore, through regression and error analyses of the standard inflow performance data, the correlation coefficients and error distributions of six types of IPR equations applicable to deviated and horizontal wells at different inclination angles were compared. The results show that the IPR equation established in this study features good stability and accuracy and that it can fully reflect the particularity of low-permeability reservoir seepage. It provides the best choice of the IPR between inclined wells and horizontal wells in low-permeability reservoirs. The other types of IPR equations are the Wiggins-Russell-Jennings, Klins-Majcher, Vogel, Fetkovich, Bendakhlia-Aziz, and Harrison equations, listed here in order from good to poor in accuracy.


---

## Body

## 1. Introduction

Research on horizontal well productivity began in the 1950s. First, people studied horizontal wells using electrical simulation in the laboratory. In 1958, Merkulov [1] first proposed a mathematical analysis method for calculating the production of horizontal wells. Then, in 1964, Borisov [2] systematically summarized the development process and production principle of horizontal wells and inclined wells and proposed an equivalent seepage resistance method. Since the 1980s, Giger et al. [3, 4] and Joshi [5], researchers in the United States [6–8], and Xu [9] and Li et al. [10], have derived various horizontal well productivity formulas.Many methods are used to study the productivity of horizontal wells. Two representative methods are the (semi)analytical method and the numerical simulation method. Analytical and semianalytical methods include applying integral transformation to solve the partial differential equation of seepage, removing the time variable through a Laplace transform, and obtaining a solution of a specific problem in the space variable using the classical Fourier method [11, 12]. For example, the equivalent seepage resistance method is used to establish the inflow performance relationship (IPR) of horizontal wells. The flow resistance of a horizontal well can be divided into two parts: the external resistance when flowing from the elliptic boundary to an imaginary vertical fracture of the horizontal well and the internal resistance of the oil well itself [5, 10]. The numerical simulation method generally uses a numerical simulator to calculate and predict the productivity of horizontal wells based on the established geological and numerical models. To verify the correctness of the numerical simulation results, it is necessary to maintain consistency with the analytical solution. The advantage is that the flexible grid system can be used to consider reservoir heterogeneity and overcome some limitations of the analytical solution [13–15].At present, the productivity analysis of deviated and horizontal wells, the determination of a reasonable working system, and the design of lifting technology [16–20] are mainly based on a series of Vogel IPR equations under different well inclination angles proposed by Cheng in 1990 [21] on the basis of reservoir numerical simulations. It is not clear whether other commonly used IPR equations, such as the Vogel equation (derived in 1968, 2012) [22–24], the Klins-Majcher equation (derived in 1989) [25], the Bendakhlia-Aziz equation (derived in 1989) [26], and the Wiggins-Russell-Jennings equation (derived in 1992) [27] (listed in Table 1), can be used for the productivity evaluations of inclined or horizontal wells.Table 1
Types and equations of inflow performance of oil wells.
TypeInflow performance equationComments on a formVogel type (type 1)qo/qomax=1−Vvpwf/pr−1−Vvpwf/pr2Vv=1+0.5bpr/a−1Klins-Majcher type (type 2)qo/qomax=1−Vkpwf/pr−1−Vkpwf/prnVk=1+bprm/1+ma−1Fetkovich type (type 3)qo/qomax=1−pwf/pr2nChange range ofn is 0.5,1Bendakhlia-Aziz type (type 4)qo/qomax=1−Vvpwf/pr−1−Vvpwf/pr2nWiggins-Russell-Jennings type (type 5)qo/qomax=1−Vw1pwf/pr−Vw2pwf/pr2−Vw3pwf/pr3−Vw4pwf/pr4V1=a0pr/aopr+a1pr2/2+a2pr3/3+a3pr4/4V2=a1pr2/2aopr+a1pr2+2a2pr3/3+a3pr4/2kro/μoBo=a0+a1p+a2p2+a3p3V3=a2pr3/3aopr+3a1pr2/2+a2pr3+3a3pr4/4V4=a3pr4/4aopr+2a1pr2+4a2pr3/3+a3pr4Harrison type (type 6)qo/qomax=1+Vh−VhexpChpwf/prVh=exp−bpe/1−exp−bprCh=bprThe IPR proposed (type 7)q/qmax=Ef′g1−M1pwf/pr−M2pwf/pr2−M3pwf/pr3−M4pwf/prm+1−M5pwf/prm+2−M6pwf/prm+3The purpose of this study was to determine the optimal IPR equations for horizontal and inclined wells and to provide a better evaluation tool for field productivity analysis. In this study, the productivity analysis solutions of single-phase and oil-gas two-phase horizontal wells were established by using equivalent seepage resistance and were transformed into the IPR equation in a dimensionless manner. Using the standardized inflow performance data of Cheng in 1990 [21], which is suitable for low-permeability reservoirs, the calculation errors of six IPR equations and the IPR equations established in this study under different well inclination angles were comprehensively compared.


## 2. Inflow Performance Relationship of Horizontal Wells in Low-Permeability Deformation Medium Reservoir

### 2.1. Establishment of Single-Phase Flow Stability Model

A schematic diagram of the horizontal well is shown in Figure1. The basic assumptions used in the derivation are as follows:
(1)
The reservoir is horizontal, has equal thickness, and has a top and bottom impermeable barrier(2)
There is a horizontal well parallel to the top and bottom in the middle of the reservoir(3)
The vertical plane flow is equivalent to a sink point in a parallel plate channel(4)
The formation is homogeneous and isotropic, ignoring the influence of gravity and capillary force(5)
The starting pressure gradient in the horizontal plane is the same as that in the vertical planeFigure 1
Schematic diagram of a horizontal well in top-bottom closed boundary reservoir.The flow resistance of a horizontal well can be divided into two parts: the external resistance when the flow is from the elliptic boundary to the imaginary vertical fracture of the horizontal well and the internal resistance of the oil well itself.


#### 2.1.1. External Resistance of the Horizontal Well

The conformal transformation diagram used to derive the external resistance of horizontal wells under elliptic boundary conditions is shown in Figure2. It is supposed that the semimajor axis of the drainage ellipse is a′, the semiminor axis is b′, the focal length is c, and the length of the horizontal well is L.Figure 2
Conformal mapping of horizontal well with elliptic boundary.The Zhukovsky transform is taken as follows:(1)z=c2w+1w.In this function, the circumference of an ellipse on planez is transformed into a circle of radius R=2a′+b′/L on plane w, and the length of a horizontal well in plane z is transformed into a circle of unit radius on plane w. Furthermore, the difference in the potential function is obtained as follows [10]:
(2)ϕR−ϕ1=q2πlna′+a′2−L2/4L/2,because
(3)dϕdr=−Kμodpdr−G0.If the permeability changes with pressure in accordance with the power exponent law [13]
(4)KK0=e−αkpr−p.The integral of equation (3) over rw,ab is as follows:
(5)qh2πlna′+a′2−L2/4L/2=K0μoαk1−exp−αkpr−pwf+αkGoa′b′−rw.Considering the influence of pollution, the deformation results are as follows:(6)qh=2πK0μo1−exp−αkpr−pwf−Goa′b′−rwαklna′+a′2−L2/4/L/2+sh.Therefore, the external resistance is(7)Rh=μo2πK0hαklna′+a′2−L2/4/L/2+shΔp1−exp−αkpr−pwf−Goa′b′−rw.


#### 2.1.2. Internal Resistance System

The conformal transformationw is taken, and its corresponding relationship is shown in Figure 3.
(8)w=eπ/hz=eπ/hxcosπyh+sinπyh=ζ+iη.Figure 3
Correspondence diagram of conformal transformation.In this transformation, the strip field−h/2<y<h/2 on the z-plane becomes the right half plane on the w-plane, and the origin on the z-plane becomes the point B1,0 on the w-plane. The complex potential of a point can be obtained by mirror image reflection in the w-plane:
(9)W=q2πlnw−1w+1+C1,W=q4πlnζ2+η2+12−4ζ2+iq2πarctanηζ−1+arctanηζ+1+C.2Therefore, the potential function is(10)ϕ=q4πlne2πx/h+12−4e2πx/hcos2πyh+C2.Whenx=0,y=h/2, the potential function at the boundary is
(11)ϕh/2=q2πln2+C2.The potential on the well wall of point sinkB can be considered the superposition of point sink A and point sink B, and the radius of point sink B is
(12)ρ=dw−1w+1dz0,0,rw=2πrwh.Then, the potential on the well wall of point sinkB is as follows:
(13)φwf=q2πln2+ln2πrwh+C2.If equation (3) is integrated within rw,h/2, the following results are obtained:
(14)Kμoαk1−exp−αkph/2−pwf+αkG0h2−rw=ϕh/2−ϕwf.By substituting equations (11) and (13) into equation (14) and considering the impact of pollution, it is concluded that
(15)qv=q=2πKμ1−exp−αkph/2−pwf−G0h/2−rwαklnh/2πrw+sv.According to equation (15), the internal resistance is as follows:
(16)Rv=μo2πK0LΔpαklnh/2πrw+sv1−exp−αkph/2−pwf−G0h/2−rw.Here,ph/2=pr is taken, because
(17)RV+RH=ΔpQ.Therefore, the horizontal well production formula is as follows:(18)Q=2πK0μoαk/Llnh/2πrw+sv/1−exp−αkpr−pwf−G0h/2−rw+μoαk/hlna′+a′2−L2/4/L/2+sh/1−exp−αkpr−pwf−Goa′b′−rw.


### 2.2. Establishment of Oil-Gas Two-Phase Flow Model

#### 2.2.1. Generalization of the Mobility Function Relation

The commonly used empirical expression of mobility functionKro/μoBo is Kro/μoBo=fp=a+bp, but the mobility function relationship in an actual oilfield is not approximately a straight line, but rather a convex or concave curve.In Figure4, the mobility curve of a well in a tight reservoir is shown. At p≤pb, comparing the curve and the linear fitting formula shows that the linear fitting relationship is not in accordance with the experimental data, and it cannot accurately reflect the relationship between Kro/μoBo and pressure. Therefore, the mobility function can be generalized as follows [27]:
(19)KroμoBo=fp=a+bpm,p≤pb,a,a+bp,p>pb,b.Figure 4
Mobility function curve.


#### 2.2.2. External Resistance System

We have the following hypothesis:(20)dH=kroμBodp.Equation (3) is arranged as follows:
(21)dϕdp=KdHdp−KGodHdL.By using a double integral and considering the influence of pollution, the horizontal production is obtained as follows:(22)qh=2πK0αklna′+a′2−L2/4/L/2+sh1−exp−αkpr−pwf+αkGoa′b′−rw×apr+bm+1prm+1−apwf−bm+1pwfm+1pr−pwf−1.Then, the external resistance is(23)Rh=Δp2αklna′+a′2−L2/4/L/2+sh2πK0h1−exp−αkpr−pwf+αkGoa′b′−rwapr+b/m+1prm+1−apwf−b/m+1pwfm+1.


#### 2.2.3. Internal Resistance System

By double integration of equation (21) and considering the influence of pollution, the yield of the vertical plane can be obtained as follows:
(24)qv=2πK0αklnh/2πrw+sv1−exp−αkph/2−pwf+αkGoh2−rw×aph/2+bm+1ph/2m+1−apwf−bm+1pwfm+1ph/2−pwf−1.Ifph/2=pr, the internal resistance is
(25)Rv=Δp2αklnh/2πrw+sv2πK0L1−exp−αkpr−pwf+αkGoh/2−rwapr+b/m+1prm+1−apwf−b/m+1pwfm+1,because
(26)RV+RH=ΔpQ.Therefore, the corresponding two-phase flow rate formula for the deformation medium reservoir is as follows:(27)1Q=lnh/2πrw+svpr−pwf−Goh/2−rw+1/2αkpr−pwf−Goh/2−rw2+1/6αk2pr−pwf−Goh/2−rw3×pr−pwf2πK0Lapr+b/m+1prm+1−apwf−b/m+1pwfm+1+pr−pwf2πK0hapr+b/m+1prm+1−apwf−b/m+1pwfm+1+lna′+a′2−L2/4/L/2+shpr−pwf−Goa′b′−rw+1/2αkpr−pwf−Goa′b′−rw2+1/6αk2pr−pwf−Goa′b′−rw3.


### 2.3. Inflow Performance of Oil and Gas Phases

Owing to the complexity of the two-phase flow production formula for horizontal wells in low-permeability and deformable media reservoirs, it is difficult to obtain the inflow performance equation with a simple form and strong generality by the usual mathematical means. In this study, the scaling method is used to deduce the inflow dynamic equation with a simple form and strong universality.


#### 2.3.1. The Reservoir Is Vertically Magnified

Ifh/2=ab, then equation (27) shows that
(28)1Q=pr−pwf2πK0Lapr+b/m+1prm+1−apwf−b/m+1pwfm+1lnh2πrw+sv+lna′+a′2−L2/4L/2+sh×pr−pwfpr−pwf−Goa′b′−rw+1/2αkpr−pwf−Goa′b′−rw2+1/6αk2pr−pwf−Goa′b′−rw3.The dimensionless equation of vertical enlargement is obtained as follows:(29)QQmax=Ef′g1−M1pwfpr−M2pwfpr2−M3pwfpr3−M4pwfprm+1−M5pwfprm+2−M6pwfprm+3.In that formula,(30)Ef′=1−hsv+Lshhlnh/2πrw+sv+Llna′+a′2−L2/4/L/2+sh,g=V+1−V11−pwf/pr,V=prpr−G0re−rw,K1=3αkpr+2αk2prpr−G0ab−rw6K3,K2=−αk2pr2K3,K3=1+12αkpr−Goa′b′−rw+16αk2pr−Goa′b′−rw2,U1=aa+b/m+1prm+1,M1=U1+K1,M2=−K1U1+K2,M3=−K2U1,M4=1−U1,M5=−K11−U1,M6=K21−U1.


#### 2.3.2. The Reservoir Shrinks Vertically

It is assumed that the oil layer is thin enough and that the vertical plane has no contribution to the production, i.e.,Rv=∞. The dimensionless equation of the vertical shrinkage of the reservoir can be obtained as follows:
(31)QQmax=Ef″g1−M1pwfpr−M2pwfpr2−M3pwfpr3−M4pwfprm+1−M5pwfprm+2−M6pwfprm+3.In the formula, we have(32)Ef″=lna′+a′2−L2/4/L/2lna′+a′2−L2/4/L/2+Sh.Comparing equations (29) and (31) reveals that the form of the inflow performance equation does not change when the liquid supply radius of the vertical plane changes from 0 to the maximum; however, the difference lies in the difference in flow efficiency. Because the flow efficiency in equation (31) is a special case result of equation (29), the IPR equation of the two-phase flow in horizontal wells of the deformation medium reservoir can be determined using equation (29).The first two parts of equation (29) are the flow efficiency terms considering the pollution degree and the additional term of seepage resistance considering the starting pressure gradient. When pollution and the starting pressure gradient are not considered, Ef=1, g=1. An exceptional case is when pwf=pr, because there is no pressure difference and no fluid flow, and the starting pressure gradient does not work, so G is equal to 0, then g=1. The third part, when m=0, can be transformed into the Wiggins-Russell-Jennings equation, and the deformation capacity of the medium is reflected in M1–M6. Therefore, the factors considered in the general formula are more comprehensive and maintain the continuity of theoretical research with the Wiggins-Russell-Jennings equation, which can be extended to more general cases.


## 2.1. Establishment of Single-Phase Flow Stability Model

A schematic diagram of the horizontal well is shown in Figure1. The basic assumptions used in the derivation are as follows:
(1)
The reservoir is horizontal, has equal thickness, and has a top and bottom impermeable barrier(2)
There is a horizontal well parallel to the top and bottom in the middle of the reservoir(3)
The vertical plane flow is equivalent to a sink point in a parallel plate channel(4)
The formation is homogeneous and isotropic, ignoring the influence of gravity and capillary force(5)
The starting pressure gradient in the horizontal plane is the same as that in the vertical planeFigure 1
Schematic diagram of a horizontal well in top-bottom closed boundary reservoir.The flow resistance of a horizontal well can be divided into two parts: the external resistance when the flow is from the elliptic boundary to the imaginary vertical fracture of the horizontal well and the internal resistance of the oil well itself.


### 2.1.1. External Resistance of the Horizontal Well

The conformal transformation diagram used to derive the external resistance of horizontal wells under elliptic boundary conditions is shown in Figure2. It is supposed that the semimajor axis of the drainage ellipse is a′, the semiminor axis is b′, the focal length is c, and the length of the horizontal well is L.Figure 2
Conformal mapping of horizontal well with elliptic boundary.The Zhukovsky transform is taken as follows:(1)z=c2w+1w.In this function, the circumference of an ellipse on planez is transformed into a circle of radius R=2a′+b′/L on plane w, and the length of a horizontal well in plane z is transformed into a circle of unit radius on plane w. Furthermore, the difference in the potential function is obtained as follows [10]:
(2)ϕR−ϕ1=q2πlna′+a′2−L2/4L/2,because
(3)dϕdr=−Kμodpdr−G0.If the permeability changes with pressure in accordance with the power exponent law [13]
(4)KK0=e−αkpr−p.The integral of equation (3) over rw,ab is as follows:
(5)qh2πlna′+a′2−L2/4L/2=K0μoαk1−exp−αkpr−pwf+αkGoa′b′−rw.Considering the influence of pollution, the deformation results are as follows:(6)qh=2πK0μo1−exp−αkpr−pwf−Goa′b′−rwαklna′+a′2−L2/4/L/2+sh.Therefore, the external resistance is(7)Rh=μo2πK0hαklna′+a′2−L2/4/L/2+shΔp1−exp−αkpr−pwf−Goa′b′−rw.


### 2.1.2. Internal Resistance System

The conformal transformationw is taken, and its corresponding relationship is shown in Figure 3.
(8)w=eπ/hz=eπ/hxcosπyh+sinπyh=ζ+iη.Figure 3
Correspondence diagram of conformal transformation.In this transformation, the strip field−h/2<y<h/2 on the z-plane becomes the right half plane on the w-plane, and the origin on the z-plane becomes the point B1,0 on the w-plane. The complex potential of a point can be obtained by mirror image reflection in the w-plane:
(9)W=q2πlnw−1w+1+C1,W=q4πlnζ2+η2+12−4ζ2+iq2πarctanηζ−1+arctanηζ+1+C.2Therefore, the potential function is(10)ϕ=q4πlne2πx/h+12−4e2πx/hcos2πyh+C2.Whenx=0,y=h/2, the potential function at the boundary is
(11)ϕh/2=q2πln2+C2.The potential on the well wall of point sinkB can be considered the superposition of point sink A and point sink B, and the radius of point sink B is
(12)ρ=dw−1w+1dz0,0,rw=2πrwh.Then, the potential on the well wall of point sinkB is as follows:
(13)φwf=q2πln2+ln2πrwh+C2.If equation (3) is integrated within rw,h/2, the following results are obtained:
(14)Kμoαk1−exp−αkph/2−pwf+αkG0h2−rw=ϕh/2−ϕwf.By substituting equations (11) and (13) into equation (14) and considering the impact of pollution, it is concluded that
(15)qv=q=2πKμ1−exp−αkph/2−pwf−G0h/2−rwαklnh/2πrw+sv.According to equation (15), the internal resistance is as follows:
(16)Rv=μo2πK0LΔpαklnh/2πrw+sv1−exp−αkph/2−pwf−G0h/2−rw.Here,ph/2=pr is taken, because
(17)RV+RH=ΔpQ.Therefore, the horizontal well production formula is as follows:(18)Q=2πK0μoαk/Llnh/2πrw+sv/1−exp−αkpr−pwf−G0h/2−rw+μoαk/hlna′+a′2−L2/4/L/2+sh/1−exp−αkpr−pwf−Goa′b′−rw.


## 2.1.1. External Resistance of the Horizontal Well

The conformal transformation diagram used to derive the external resistance of horizontal wells under elliptic boundary conditions is shown in Figure2. It is supposed that the semimajor axis of the drainage ellipse is a′, the semiminor axis is b′, the focal length is c, and the length of the horizontal well is L.Figure 2
Conformal mapping of horizontal well with elliptic boundary.The Zhukovsky transform is taken as follows:(1)z=c2w+1w.In this function, the circumference of an ellipse on planez is transformed into a circle of radius R=2a′+b′/L on plane w, and the length of a horizontal well in plane z is transformed into a circle of unit radius on plane w. Furthermore, the difference in the potential function is obtained as follows [10]:
(2)ϕR−ϕ1=q2πlna′+a′2−L2/4L/2,because
(3)dϕdr=−Kμodpdr−G0.If the permeability changes with pressure in accordance with the power exponent law [13]
(4)KK0=e−αkpr−p.The integral of equation (3) over rw,ab is as follows:
(5)qh2πlna′+a′2−L2/4L/2=K0μoαk1−exp−αkpr−pwf+αkGoa′b′−rw.Considering the influence of pollution, the deformation results are as follows:(6)qh=2πK0μo1−exp−αkpr−pwf−Goa′b′−rwαklna′+a′2−L2/4/L/2+sh.Therefore, the external resistance is(7)Rh=μo2πK0hαklna′+a′2−L2/4/L/2+shΔp1−exp−αkpr−pwf−Goa′b′−rw.


## 2.1.2. Internal Resistance System

The conformal transformationw is taken, and its corresponding relationship is shown in Figure 3.
(8)w=eπ/hz=eπ/hxcosπyh+sinπyh=ζ+iη.Figure 3
Correspondence diagram of conformal transformation.In this transformation, the strip field−h/2<y<h/2 on the z-plane becomes the right half plane on the w-plane, and the origin on the z-plane becomes the point B1,0 on the w-plane. The complex potential of a point can be obtained by mirror image reflection in the w-plane:
(9)W=q2πlnw−1w+1+C1,W=q4πlnζ2+η2+12−4ζ2+iq2πarctanηζ−1+arctanηζ+1+C.2Therefore, the potential function is(10)ϕ=q4πlne2πx/h+12−4e2πx/hcos2πyh+C2.Whenx=0,y=h/2, the potential function at the boundary is
(11)ϕh/2=q2πln2+C2.The potential on the well wall of point sinkB can be considered the superposition of point sink A and point sink B, and the radius of point sink B is
(12)ρ=dw−1w+1dz0,0,rw=2πrwh.Then, the potential on the well wall of point sinkB is as follows:
(13)φwf=q2πln2+ln2πrwh+C2.If equation (3) is integrated within rw,h/2, the following results are obtained:
(14)Kμoαk1−exp−αkph/2−pwf+αkG0h2−rw=ϕh/2−ϕwf.By substituting equations (11) and (13) into equation (14) and considering the impact of pollution, it is concluded that
(15)qv=q=2πKμ1−exp−αkph/2−pwf−G0h/2−rwαklnh/2πrw+sv.According to equation (15), the internal resistance is as follows:
(16)Rv=μo2πK0LΔpαklnh/2πrw+sv1−exp−αkph/2−pwf−G0h/2−rw.Here,ph/2=pr is taken, because
(17)RV+RH=ΔpQ.Therefore, the horizontal well production formula is as follows:(18)Q=2πK0μoαk/Llnh/2πrw+sv/1−exp−αkpr−pwf−G0h/2−rw+μoαk/hlna′+a′2−L2/4/L/2+sh/1−exp−αkpr−pwf−Goa′b′−rw.


## 2.2. Establishment of Oil-Gas Two-Phase Flow Model

### 2.2.1. Generalization of the Mobility Function Relation

The commonly used empirical expression of mobility functionKro/μoBo is Kro/μoBo=fp=a+bp, but the mobility function relationship in an actual oilfield is not approximately a straight line, but rather a convex or concave curve.In Figure4, the mobility curve of a well in a tight reservoir is shown. At p≤pb, comparing the curve and the linear fitting formula shows that the linear fitting relationship is not in accordance with the experimental data, and it cannot accurately reflect the relationship between Kro/μoBo and pressure. Therefore, the mobility function can be generalized as follows [27]:
(19)KroμoBo=fp=a+bpm,p≤pb,a,a+bp,p>pb,b.Figure 4
Mobility function curve.


### 2.2.2. External Resistance System

We have the following hypothesis:(20)dH=kroμBodp.Equation (3) is arranged as follows:
(21)dϕdp=KdHdp−KGodHdL.By using a double integral and considering the influence of pollution, the horizontal production is obtained as follows:(22)qh=2πK0αklna′+a′2−L2/4/L/2+sh1−exp−αkpr−pwf+αkGoa′b′−rw×apr+bm+1prm+1−apwf−bm+1pwfm+1pr−pwf−1.Then, the external resistance is(23)Rh=Δp2αklna′+a′2−L2/4/L/2+sh2πK0h1−exp−αkpr−pwf+αkGoa′b′−rwapr+b/m+1prm+1−apwf−b/m+1pwfm+1.


### 2.2.3. Internal Resistance System

By double integration of equation (21) and considering the influence of pollution, the yield of the vertical plane can be obtained as follows:
(24)qv=2πK0αklnh/2πrw+sv1−exp−αkph/2−pwf+αkGoh2−rw×aph/2+bm+1ph/2m+1−apwf−bm+1pwfm+1ph/2−pwf−1.Ifph/2=pr, the internal resistance is
(25)Rv=Δp2αklnh/2πrw+sv2πK0L1−exp−αkpr−pwf+αkGoh/2−rwapr+b/m+1prm+1−apwf−b/m+1pwfm+1,because
(26)RV+RH=ΔpQ.Therefore, the corresponding two-phase flow rate formula for the deformation medium reservoir is as follows:(27)1Q=lnh/2πrw+svpr−pwf−Goh/2−rw+1/2αkpr−pwf−Goh/2−rw2+1/6αk2pr−pwf−Goh/2−rw3×pr−pwf2πK0Lapr+b/m+1prm+1−apwf−b/m+1pwfm+1+pr−pwf2πK0hapr+b/m+1prm+1−apwf−b/m+1pwfm+1+lna′+a′2−L2/4/L/2+shpr−pwf−Goa′b′−rw+1/2αkpr−pwf−Goa′b′−rw2+1/6αk2pr−pwf−Goa′b′−rw3.


## 2.2.1. Generalization of the Mobility Function Relation

The commonly used empirical expression of mobility functionKro/μoBo is Kro/μoBo=fp=a+bp, but the mobility function relationship in an actual oilfield is not approximately a straight line, but rather a convex or concave curve.In Figure4, the mobility curve of a well in a tight reservoir is shown. At p≤pb, comparing the curve and the linear fitting formula shows that the linear fitting relationship is not in accordance with the experimental data, and it cannot accurately reflect the relationship between Kro/μoBo and pressure. Therefore, the mobility function can be generalized as follows [27]:
(19)KroμoBo=fp=a+bpm,p≤pb,a,a+bp,p>pb,b.Figure 4
Mobility function curve.


## 2.2.2. External Resistance System

We have the following hypothesis:(20)dH=kroμBodp.Equation (3) is arranged as follows:
(21)dϕdp=KdHdp−KGodHdL.By using a double integral and considering the influence of pollution, the horizontal production is obtained as follows:(22)qh=2πK0αklna′+a′2−L2/4/L/2+sh1−exp−αkpr−pwf+αkGoa′b′−rw×apr+bm+1prm+1−apwf−bm+1pwfm+1pr−pwf−1.Then, the external resistance is(23)Rh=Δp2αklna′+a′2−L2/4/L/2+sh2πK0h1−exp−αkpr−pwf+αkGoa′b′−rwapr+b/m+1prm+1−apwf−b/m+1pwfm+1.


## 2.2.3. Internal Resistance System

By double integration of equation (21) and considering the influence of pollution, the yield of the vertical plane can be obtained as follows:
(24)qv=2πK0αklnh/2πrw+sv1−exp−αkph/2−pwf+αkGoh2−rw×aph/2+bm+1ph/2m+1−apwf−bm+1pwfm+1ph/2−pwf−1.Ifph/2=pr, the internal resistance is
(25)Rv=Δp2αklnh/2πrw+sv2πK0L1−exp−αkpr−pwf+αkGoh/2−rwapr+b/m+1prm+1−apwf−b/m+1pwfm+1,because
(26)RV+RH=ΔpQ.Therefore, the corresponding two-phase flow rate formula for the deformation medium reservoir is as follows:(27)1Q=lnh/2πrw+svpr−pwf−Goh/2−rw+1/2αkpr−pwf−Goh/2−rw2+1/6αk2pr−pwf−Goh/2−rw3×pr−pwf2πK0Lapr+b/m+1prm+1−apwf−b/m+1pwfm+1+pr−pwf2πK0hapr+b/m+1prm+1−apwf−b/m+1pwfm+1+lna′+a′2−L2/4/L/2+shpr−pwf−Goa′b′−rw+1/2αkpr−pwf−Goa′b′−rw2+1/6αk2pr−pwf−Goa′b′−rw3.


## 2.3. Inflow Performance of Oil and Gas Phases

Owing to the complexity of the two-phase flow production formula for horizontal wells in low-permeability and deformable media reservoirs, it is difficult to obtain the inflow performance equation with a simple form and strong generality by the usual mathematical means. In this study, the scaling method is used to deduce the inflow dynamic equation with a simple form and strong universality.


### 2.3.1. The Reservoir Is Vertically Magnified

Ifh/2=ab, then equation (27) shows that
(28)1Q=pr−pwf2πK0Lapr+b/m+1prm+1−apwf−b/m+1pwfm+1lnh2πrw+sv+lna′+a′2−L2/4L/2+sh×pr−pwfpr−pwf−Goa′b′−rw+1/2αkpr−pwf−Goa′b′−rw2+1/6αk2pr−pwf−Goa′b′−rw3.The dimensionless equation of vertical enlargement is obtained as follows:(29)QQmax=Ef′g1−M1pwfpr−M2pwfpr2−M3pwfpr3−M4pwfprm+1−M5pwfprm+2−M6pwfprm+3.In that formula,(30)Ef′=1−hsv+Lshhlnh/2πrw+sv+Llna′+a′2−L2/4/L/2+sh,g=V+1−V11−pwf/pr,V=prpr−G0re−rw,K1=3αkpr+2αk2prpr−G0ab−rw6K3,K2=−αk2pr2K3,K3=1+12αkpr−Goa′b′−rw+16αk2pr−Goa′b′−rw2,U1=aa+b/m+1prm+1,M1=U1+K1,M2=−K1U1+K2,M3=−K2U1,M4=1−U1,M5=−K11−U1,M6=K21−U1.


### 2.3.2. The Reservoir Shrinks Vertically

It is assumed that the oil layer is thin enough and that the vertical plane has no contribution to the production, i.e.,Rv=∞. The dimensionless equation of the vertical shrinkage of the reservoir can be obtained as follows:
(31)QQmax=Ef″g1−M1pwfpr−M2pwfpr2−M3pwfpr3−M4pwfprm+1−M5pwfprm+2−M6pwfprm+3.In the formula, we have(32)Ef″=lna′+a′2−L2/4/L/2lna′+a′2−L2/4/L/2+Sh.Comparing equations (29) and (31) reveals that the form of the inflow performance equation does not change when the liquid supply radius of the vertical plane changes from 0 to the maximum; however, the difference lies in the difference in flow efficiency. Because the flow efficiency in equation (31) is a special case result of equation (29), the IPR equation of the two-phase flow in horizontal wells of the deformation medium reservoir can be determined using equation (29).The first two parts of equation (29) are the flow efficiency terms considering the pollution degree and the additional term of seepage resistance considering the starting pressure gradient. When pollution and the starting pressure gradient are not considered, Ef=1, g=1. An exceptional case is when pwf=pr, because there is no pressure difference and no fluid flow, and the starting pressure gradient does not work, so G is equal to 0, then g=1. The third part, when m=0, can be transformed into the Wiggins-Russell-Jennings equation, and the deformation capacity of the medium is reflected in M1–M6. Therefore, the factors considered in the general formula are more comprehensive and maintain the continuity of theoretical research with the Wiggins-Russell-Jennings equation, which can be extended to more general cases.


## 2.3.1. The Reservoir Is Vertically Magnified

Ifh/2=ab, then equation (27) shows that
(28)1Q=pr−pwf2πK0Lapr+b/m+1prm+1−apwf−b/m+1pwfm+1lnh2πrw+sv+lna′+a′2−L2/4L/2+sh×pr−pwfpr−pwf−Goa′b′−rw+1/2αkpr−pwf−Goa′b′−rw2+1/6αk2pr−pwf−Goa′b′−rw3.The dimensionless equation of vertical enlargement is obtained as follows:(29)QQmax=Ef′g1−M1pwfpr−M2pwfpr2−M3pwfpr3−M4pwfprm+1−M5pwfprm+2−M6pwfprm+3.In that formula,(30)Ef′=1−hsv+Lshhlnh/2πrw+sv+Llna′+a′2−L2/4/L/2+sh,g=V+1−V11−pwf/pr,V=prpr−G0re−rw,K1=3αkpr+2αk2prpr−G0ab−rw6K3,K2=−αk2pr2K3,K3=1+12αkpr−Goa′b′−rw+16αk2pr−Goa′b′−rw2,U1=aa+b/m+1prm+1,M1=U1+K1,M2=−K1U1+K2,M3=−K2U1,M4=1−U1,M5=−K11−U1,M6=K21−U1.


## 2.3.2. The Reservoir Shrinks Vertically

It is assumed that the oil layer is thin enough and that the vertical plane has no contribution to the production, i.e.,Rv=∞. The dimensionless equation of the vertical shrinkage of the reservoir can be obtained as follows:
(31)QQmax=Ef″g1−M1pwfpr−M2pwfpr2−M3pwfpr3−M4pwfprm+1−M5pwfprm+2−M6pwfprm+3.In the formula, we have(32)Ef″=lna′+a′2−L2/4/L/2lna′+a′2−L2/4/L/2+Sh.Comparing equations (29) and (31) reveals that the form of the inflow performance equation does not change when the liquid supply radius of the vertical plane changes from 0 to the maximum; however, the difference lies in the difference in flow efficiency. Because the flow efficiency in equation (31) is a special case result of equation (29), the IPR equation of the two-phase flow in horizontal wells of the deformation medium reservoir can be determined using equation (29).The first two parts of equation (29) are the flow efficiency terms considering the pollution degree and the additional term of seepage resistance considering the starting pressure gradient. When pollution and the starting pressure gradient are not considered, Ef=1, g=1. An exceptional case is when pwf=pr, because there is no pressure difference and no fluid flow, and the starting pressure gradient does not work, so G is equal to 0, then g=1. The third part, when m=0, can be transformed into the Wiggins-Russell-Jennings equation, and the deformation capacity of the medium is reflected in M1–M6. Therefore, the factors considered in the general formula are more comprehensive and maintain the continuity of theoretical research with the Wiggins-Russell-Jennings equation, which can be extended to more general cases.


## 3. Optimal IPR Equation for Horizontal and Deviated Wells

### 3.1. Regression Analysis

The types of inflow performance for the comparative analysis are shown in Table1, and standardized inflow performance data [10] are selected for the regression analysis (see Table 2). The data in Table 2 are the average values of the 16 reservoir conditions. The permeability value under the third reservoir condition does not meet the upper limit standard of a low-permeability reservoir (more than 200×10−3μm2). The permeability in the simulation of the other 15 reservoirs is 20, which belongs to the category of low-permeability reservoirs. Therefore, it has strong confidence and is ideal for comparative analyses of various IPR equations.Table 2
Average normalized dimensionless production of different well inclination angles [10].
pwf/pe0°15°30°45°60°75°85°88.56°90°0.01.0001.0001.0001.0001.0001.0001.0001.0001.0000.10.9720.9700.9750.9820.9860.9890.9900.9900.9950.20.9270.9240.9350.9470.9570.9620.9640.9640.9740.30.8680.8640.8790.8970.9100.9180.9210.9210.9350.40.7910.7860.8060.8280.8440.8540.8560.8570.8750.50.7000.6950.7170.7420.7600.7710.7740.7750.7960.60.5920.5870.6110.6360.6550.6660.6690.6690.6920.70.4680.4640.4860.5100.5270.5370.5400.5400.5620.80.3280.3250.3430.3620.3770.3850.3870.3870.4050.90.1720.1710.1820.1940.2020.2070.2090.2090.2191.00.0000.0000.0000.0000.0000.0000.0000.0000.000The seven inflow performance equations in Table1 are used for the regression analysis of dimensionless production under the different well inclination angles in Table 2. The undetermined parameter values are listed in Table 3. The overall fitting degrees of different IPR equations for well deviation data are quite different. Among them, the complex correlation coefficient of the type 2 and type 7 inflow performance equations is 1; the fitting degree is the highest, and the complex correlation coefficient of type 2 is close to 1, which also shows good applicability for inclined wells with different angles. The second is type 1, and the applicability of types 3, 4, and 6 is relatively poor.Table 3
Fitting values of undetermined parameters of inflow performance equations for different types of inclined wells.
TypeParameterInclination angle0°15°30°45°60°75°85°88.56°90°1Equationqo/qomax=V0−V1pwf/pr−V2pwf/pr2V00.999810.999800.996910.994550.992570.991520.991510.991410.98845V10.200800.220950.125360.02214-0.05487-0.100230-0.111990.11411-0.20545V20.798830.778320.868180.966321.039511.082871.094171.096391.18182R21.000001.000000.999960.999860.999760.999700.999690.999700.999432Equationqo/qomax=V0−V1pwf/pr−V2pwf/prnV00.999981.000061.000401.000791.000871.000841.000951.000861.00136V10.204470.226770.198350.147550.108070.081010.071640.069740.03726V20.795460.772960.801460.852210.891700.918670.927810.929830.96297n2.006122.010012.125112.206952.260932.283542.284532.284182.37119R21.000001.000001.000000.999990.999990.999990.999990.999990.999993Equationqo/qomax=V01−pwf/pr2nV00.973830.971250.979010.989180.997291.002521.004121.004401.01443n1.076231.082971.034680.986510.952910.934160.928770.928690.88890R20.998730.998450.999260.999830.999960.999890.999860.999840.999284Equationqo/qomax=V0−V1pwf/pr−V2pwf/pr2nV00.999831.000050.998930.998120.997220.996680.997190.997270.99537V10.201050.224710.154100.069210.00311-0.03811-0.04439-0.04432-0.12852V20.798640.775280.845360.929830.995471.036331.043711.044301.12610n1.000211.003111.024951.043221.055711.061411.067411.069741.08173R21.000001.000000.999990.999980.999960.999950.999950.999950.999905Equationqo/qomax=V0−V1pwf/pr−V2pwf/pr2−V3pwf/pr3−V4pwf/pr4V01.000100.999970.999941.000050.999860.999900.999820.999831.00001V10.206260.221250.171610.105910.055770.030110.013290.01690-0.02625V20.781060.787940.755570.762210.771880.753180.795080.764950.73091V30.01865-0.026810.058080.105480.134230.196000.141610.197160.26321V4-0.005830.017480.014570.026220.037880.020400.049530.020400.03205R21.000001.000001.000001.000001.000001.000001.000001.000001.000006Equationqo/qomax=V0−V1expC⋅pwf/prV01.287311.299551.253441.195961.171061.157121.153361.138511.13049V10.269470.282650.233180.175320.147960.132820.128680.118100.10358C1.572721.533671.688471.929812.079202.176302.205072.282522.40374R20.999110.999200.999090.998880.998640.998440.998380.998240.998067Equationqo/qomax=V0−V1pwf/pr−V2pwf/pr2−V3pwf/pr3−V4pwf/prm+1−V5pwf/prm+2−V6pwf/prm+3V01.000020.999990.999981.000000.999970.999980.999980.999980.99998V10.190090.276800.174090.070000.075350.022210.02088-0.00095-0.03362V21.05611-0.377910.741781.552530.445411.045040.806971.371560.86277V314.00211-58.991560.8278449.58385-16.1307723.189148.0411946.753397.29149V4-13.579557.24579-0.67021-47.5592915.81206-21.89913-7.28745-44.3826-6.80887V5-0.820333.54701-0.11652-3.400440.96540-1.79647-0.83814-3.61157-0.39675V60.15156-0.700180.042920.75325-0.167530.439150.256410.870070.08486R21.000001.000001.000001.000001.000001.000001.000001.000001.00000


### 3.2. Error Analysis

#### 3.2.1. Error Analysis of Inflow Performance Fitting in Inclined Wells

The normalized dimensionless data with inclination angles of 30°, 60°, and 85° were selected to analyze the relative errors of the seven types of IPR equations. The statistical analysis of the relative error is shown in Figures5–8. In the entire dimensionless pressure range, the third, fourth, and sixth IPR equations have the most obvious changes. For example, when the inclination angle of type 6 is 30°, the arch shape is large at both ends and small in the middle. When the deviation angle is 60°, the variation range is reduced. A well deviation angle of 85° results in the shape of the sinusoidal curve with large fluctuation amplitude and wave crest and trough.Figure 5
Relative error analysis of inflow performance equation fitting for seven types of 30° inclined wells.Figure 6
Relative error analysis of inflow performance equation fitting for seven types of 60° inclined wells.Figure 7
Relative error analysis of inflow performance equation fitting for seven types of 85° inclined wells.Figure 8
Relative error of normalized dimensionless yield fitting of type 7 IPR equation.


#### 3.2.2. Overall Evaluation of Fitting Error

According to the analysis of fitting errors of the seven types of IPR equations, the seventh type of IPR equation established in this study has the highest fitting accuracy with a maximum relative error of less than 0.25%. The fifth type of IPR equation has the second highest accuracy with a maximum relative error of less than 0.80%. The second type of IPR equation has a maximum relative error of less than 0.80%. The first type of IPR equation performs poorly, with a maximum relative error of no more than 3.00%. Types 3, 4, and 6 had the poorest performances, with maximum relative errors of no more than 5.00%.Using the inflow performance equation established in this study, the relative error of each well inclination angle fitted by the data in Table1 is shown in Figure 8. The fitting accuracy is very high, and it has good adaptability for wells with different well inclination angles; therefore, it can be applied to the performance analysis of oil wells on site.


## 3.1. Regression Analysis

The types of inflow performance for the comparative analysis are shown in Table1, and standardized inflow performance data [10] are selected for the regression analysis (see Table 2). The data in Table 2 are the average values of the 16 reservoir conditions. The permeability value under the third reservoir condition does not meet the upper limit standard of a low-permeability reservoir (more than 200×10−3μm2). The permeability in the simulation of the other 15 reservoirs is 20, which belongs to the category of low-permeability reservoirs. Therefore, it has strong confidence and is ideal for comparative analyses of various IPR equations.Table 2
Average normalized dimensionless production of different well inclination angles [10].
pwf/pe0°15°30°45°60°75°85°88.56°90°0.01.0001.0001.0001.0001.0001.0001.0001.0001.0000.10.9720.9700.9750.9820.9860.9890.9900.9900.9950.20.9270.9240.9350.9470.9570.9620.9640.9640.9740.30.8680.8640.8790.8970.9100.9180.9210.9210.9350.40.7910.7860.8060.8280.8440.8540.8560.8570.8750.50.7000.6950.7170.7420.7600.7710.7740.7750.7960.60.5920.5870.6110.6360.6550.6660.6690.6690.6920.70.4680.4640.4860.5100.5270.5370.5400.5400.5620.80.3280.3250.3430.3620.3770.3850.3870.3870.4050.90.1720.1710.1820.1940.2020.2070.2090.2090.2191.00.0000.0000.0000.0000.0000.0000.0000.0000.000The seven inflow performance equations in Table1 are used for the regression analysis of dimensionless production under the different well inclination angles in Table 2. The undetermined parameter values are listed in Table 3. The overall fitting degrees of different IPR equations for well deviation data are quite different. Among them, the complex correlation coefficient of the type 2 and type 7 inflow performance equations is 1; the fitting degree is the highest, and the complex correlation coefficient of type 2 is close to 1, which also shows good applicability for inclined wells with different angles. The second is type 1, and the applicability of types 3, 4, and 6 is relatively poor.Table 3
Fitting values of undetermined parameters of inflow performance equations for different types of inclined wells.
TypeParameterInclination angle0°15°30°45°60°75°85°88.56°90°1Equationqo/qomax=V0−V1pwf/pr−V2pwf/pr2V00.999810.999800.996910.994550.992570.991520.991510.991410.98845V10.200800.220950.125360.02214-0.05487-0.100230-0.111990.11411-0.20545V20.798830.778320.868180.966321.039511.082871.094171.096391.18182R21.000001.000000.999960.999860.999760.999700.999690.999700.999432Equationqo/qomax=V0−V1pwf/pr−V2pwf/prnV00.999981.000061.000401.000791.000871.000841.000951.000861.00136V10.204470.226770.198350.147550.108070.081010.071640.069740.03726V20.795460.772960.801460.852210.891700.918670.927810.929830.96297n2.006122.010012.125112.206952.260932.283542.284532.284182.37119R21.000001.000001.000000.999990.999990.999990.999990.999990.999993Equationqo/qomax=V01−pwf/pr2nV00.973830.971250.979010.989180.997291.002521.004121.004401.01443n1.076231.082971.034680.986510.952910.934160.928770.928690.88890R20.998730.998450.999260.999830.999960.999890.999860.999840.999284Equationqo/qomax=V0−V1pwf/pr−V2pwf/pr2nV00.999831.000050.998930.998120.997220.996680.997190.997270.99537V10.201050.224710.154100.069210.00311-0.03811-0.04439-0.04432-0.12852V20.798640.775280.845360.929830.995471.036331.043711.044301.12610n1.000211.003111.024951.043221.055711.061411.067411.069741.08173R21.000001.000000.999990.999980.999960.999950.999950.999950.999905Equationqo/qomax=V0−V1pwf/pr−V2pwf/pr2−V3pwf/pr3−V4pwf/pr4V01.000100.999970.999941.000050.999860.999900.999820.999831.00001V10.206260.221250.171610.105910.055770.030110.013290.01690-0.02625V20.781060.787940.755570.762210.771880.753180.795080.764950.73091V30.01865-0.026810.058080.105480.134230.196000.141610.197160.26321V4-0.005830.017480.014570.026220.037880.020400.049530.020400.03205R21.000001.000001.000001.000001.000001.000001.000001.000001.000006Equationqo/qomax=V0−V1expC⋅pwf/prV01.287311.299551.253441.195961.171061.157121.153361.138511.13049V10.269470.282650.233180.175320.147960.132820.128680.118100.10358C1.572721.533671.688471.929812.079202.176302.205072.282522.40374R20.999110.999200.999090.998880.998640.998440.998380.998240.998067Equationqo/qomax=V0−V1pwf/pr−V2pwf/pr2−V3pwf/pr3−V4pwf/prm+1−V5pwf/prm+2−V6pwf/prm+3V01.000020.999990.999981.000000.999970.999980.999980.999980.99998V10.190090.276800.174090.070000.075350.022210.02088-0.00095-0.03362V21.05611-0.377910.741781.552530.445411.045040.806971.371560.86277V314.00211-58.991560.8278449.58385-16.1307723.189148.0411946.753397.29149V4-13.579557.24579-0.67021-47.5592915.81206-21.89913-7.28745-44.3826-6.80887V5-0.820333.54701-0.11652-3.400440.96540-1.79647-0.83814-3.61157-0.39675V60.15156-0.700180.042920.75325-0.167530.439150.256410.870070.08486R21.000001.000001.000001.000001.000001.000001.000001.000001.00000


## 3.2. Error Analysis

### 3.2.1. Error Analysis of Inflow Performance Fitting in Inclined Wells

The normalized dimensionless data with inclination angles of 30°, 60°, and 85° were selected to analyze the relative errors of the seven types of IPR equations. The statistical analysis of the relative error is shown in Figures5–8. In the entire dimensionless pressure range, the third, fourth, and sixth IPR equations have the most obvious changes. For example, when the inclination angle of type 6 is 30°, the arch shape is large at both ends and small in the middle. When the deviation angle is 60°, the variation range is reduced. A well deviation angle of 85° results in the shape of the sinusoidal curve with large fluctuation amplitude and wave crest and trough.Figure 5
Relative error analysis of inflow performance equation fitting for seven types of 30° inclined wells.Figure 6
Relative error analysis of inflow performance equation fitting for seven types of 60° inclined wells.Figure 7
Relative error analysis of inflow performance equation fitting for seven types of 85° inclined wells.Figure 8
Relative error of normalized dimensionless yield fitting of type 7 IPR equation.


### 3.2.2. Overall Evaluation of Fitting Error

According to the analysis of fitting errors of the seven types of IPR equations, the seventh type of IPR equation established in this study has the highest fitting accuracy with a maximum relative error of less than 0.25%. The fifth type of IPR equation has the second highest accuracy with a maximum relative error of less than 0.80%. The second type of IPR equation has a maximum relative error of less than 0.80%. The first type of IPR equation performs poorly, with a maximum relative error of no more than 3.00%. Types 3, 4, and 6 had the poorest performances, with maximum relative errors of no more than 5.00%.Using the inflow performance equation established in this study, the relative error of each well inclination angle fitted by the data in Table1 is shown in Figure 8. The fitting accuracy is very high, and it has good adaptability for wells with different well inclination angles; therefore, it can be applied to the performance analysis of oil wells on site.


## 3.2.1. Error Analysis of Inflow Performance Fitting in Inclined Wells

The normalized dimensionless data with inclination angles of 30°, 60°, and 85° were selected to analyze the relative errors of the seven types of IPR equations. The statistical analysis of the relative error is shown in Figures5–8. In the entire dimensionless pressure range, the third, fourth, and sixth IPR equations have the most obvious changes. For example, when the inclination angle of type 6 is 30°, the arch shape is large at both ends and small in the middle. When the deviation angle is 60°, the variation range is reduced. A well deviation angle of 85° results in the shape of the sinusoidal curve with large fluctuation amplitude and wave crest and trough.Figure 5
Relative error analysis of inflow performance equation fitting for seven types of 30° inclined wells.Figure 6
Relative error analysis of inflow performance equation fitting for seven types of 60° inclined wells.Figure 7
Relative error analysis of inflow performance equation fitting for seven types of 85° inclined wells.Figure 8
Relative error of normalized dimensionless yield fitting of type 7 IPR equation.


## 3.2.2. Overall Evaluation of Fitting Error

According to the analysis of fitting errors of the seven types of IPR equations, the seventh type of IPR equation established in this study has the highest fitting accuracy with a maximum relative error of less than 0.25%. The fifth type of IPR equation has the second highest accuracy with a maximum relative error of less than 0.80%. The second type of IPR equation has a maximum relative error of less than 0.80%. The first type of IPR equation performs poorly, with a maximum relative error of no more than 3.00%. Types 3, 4, and 6 had the poorest performances, with maximum relative errors of no more than 5.00%.Using the inflow performance equation established in this study, the relative error of each well inclination angle fitted by the data in Table1 is shown in Figure 8. The fitting accuracy is very high, and it has good adaptability for wells with different well inclination angles; therefore, it can be applied to the performance analysis of oil wells on site.


## 4. Conclusions

(1)
In past studies, the linear flow function relationship has often been used to characterize the fluid pressure-volume-temperature model. Based on the results of the pressure-volume-temperature analysis of real reservoirs, it is considered that the improvement of the linear flow function to the power exponential flow function is helpful to describe the pressure-volume-temperature relationship of the fluid accurately under different pressure ranges and to improve the prediction accuracy of the IPR(2)
To reflect the special percolation law of a low-permeability reservoirs fully, the threshold pressure gradient effect and the pressure sensitivity of reservoir permeability were considered in the establishment of the IPR equation. This was different from previous research and resulted in some progress. In the dimensionless inflow dynamic relationship of the oil-gas two-phase model established by the analytical method, the influence of the starting pressure gradient is reflected inG, and the deformation ability of the medium is reflected in M1–M6. When m=0, it can be transformed into the Wiggins-Russell-Jennings equation. The factors considered in the equation are more comprehensive and maintain the continuity of theoretical research with the Wiggins-Russell-Jennings equation, which can be extended to more general cases(3)
Through regression and error analyses of standardized inflow dynamic data, the correlation coefficients and error distribution of six IPR equations for inclined and horizontal wells and the model established in this study were compared under different well inclination angles. The results show that the IPR equation established in this study is the best choice for the inflow dynamic relationship of inclined and horizontal wells in a low-permeability reservoir. The other types of IPR equations can be listed in order from good to poor as follows: Wiggins-Russell-Jennings, Klins-Majcher, Vogel, Fetkovich, Bendakhlia-Aziz, and Harrison.



---

*Source: 6640871-2021-03-18.xml*