Next Article in Journal
Measuring Energy Efficiency and Environmental Performance: A Case of South Asia
Next Article in Special Issue
The Influence and Optimization of Geometrical Parameters on Coast-Down Characteristics of Nuclear Reactor Coolant Pumps
Previous Article in Journal
Optimization of Distillation Sequences with Nonsharp Separation Columns
Previous Article in Special Issue
Off-Grid Solar PV Power Generation System in Sindh, Pakistan: A Techno-Economic Feasibility Analysis
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Productivity Models of Infill Complex Structural Wells in Mixed Well Patterns

Research Institute of Petroleum Exploration & Development, PetroChina, Beijing 100083, China
*
Author to whom correspondence should be addressed.
Processes 2019, 7(6), 324; https://doi.org/10.3390/pr7060324
Submission received: 14 May 2019 / Revised: 24 May 2019 / Accepted: 27 May 2019 / Published: 31 May 2019

Abstract

:
The mathematical models of productivity calculation for complex structural wells mainly focus on the single well or the regular well pattern. Previous research on the seepage theory of complex structural wells and vertical wells in mixed well pattern is greatly insufficient. Accordingly, this article presents a methodology of evaluating the productivity of infill complex structural wells in mixed well patterns. On the basis of the mirror-image method and source–sink theory, two semi-analytical models are established. These models are applied to the productivity prediction of an infill horizontal well inhorizontal-vertical well pattern and an infill multilateral well inmultilateral-vertical well pattern, respectively, in which the interference of other wells, the randomicity of well patterns, and the pressure drawdown along the horizontal laterals are taken into account. The semi-analytical models’ results are consistent with those calculated by the Eclipse reservoir simulator with the relative error of less than 15%. Results indicate that the bottom hole flowing pressure decreases logarithmically while the wellbore flow rate increases monotonically from the toe to the heel of the horizontal well. Due to the pseudo-hemispherical flow at each endpoint and the pseudo-linear flow at the center of the horizontal well, the drainage area at each endpoint is relatively larger than that at the center. The radial inflow at each endpoint of the horizontal segment is considerably greater than that at the center, which presents the U-shape distribution. The proposed methodology enhances and promotes the theory of productivity evaluation for complex structural wells in mixed well patterns.

1. Introduction

Complex structural wells including horizontal wells and multilateral wells have become a popular alternative for the development of oil and gas fields around the world because of their high flow efficiency due to larger contact area made with the reservoir and lower pressure drawdown at the same liquid volume [1]. As the process of oilfield development enters into the intermediary and later phase, the mixed well patterns of complex structural wells and conventional vertical wells have been widely applied to the implementation of adjustment plans [2]. The infilling horizontal wells or multilateral wells in mixed well patterns are of great significance to optimization of development strategies. Due to the complexity of seepage mechanism near the horizontal wellbore, the coupling between reservoir flow and wellbore conduit flow, and the interference of other wells in mixed well patterns [3], the original mathematical models based on time invariant flow are no longer applicable as a result of the change of flow regimes. Therefore, it is necessary to establish new productivity models of infill complex structural wells in mixed well patterns.
The productivity evaluation of complex structural wells under different reservoir conditions is an important topic in the field of complex structural wells. At present, studies have been conducted on the methods of production calculation for horizontal wells [4,5,6,7,8,9,10,11,12,13,14,15,16,17,18] and multilateral wells [19,20,21,22,23], which are mainly divided into analytical methods and semi-analytical methods. The analytical model aims to directly build a calculation formula based on ideal assumptions. In the semi-analytical model, each branch of the complex structural well is divided into several infinitesimal sections, thus the productivity can be obtained by solving the system of linear equations combined with fluid flow rate and pressure drawdown for each infinitesimal in the wellbore. However, these methods mainly focus on building models of single well or regular well patterns, and little research has been conducted on the seepage theory of mixed well patterns of vertical wells and complex structural wells [24,25,26,27,28,29,30]. In view of the field problems involving the productivity prediction of infill wells in irregular well patterns, there are limitations in current mathematical models [31,32,33,34,35]. On the one hand, they are only suitable for the productivity calculation of the entire well pattern not for single infill wells; on the other hand, they are merely applied to the regular five-spot pattern, seven-spot pattern, and nine-spot pattern not for the irregularly mixed well pattern. More importantly, the pressure drop caused by wall friction and fluid acceleration along the horizontal lateral is not comprehensively taken into account during the coupling of reservoir seepage and wellbore conduit flow. In order to overcome the deficiencies of the existing models, Ye et al. [36] presented a productivity evaluation model for infill horizontal wells considering the interference of other wells and the wellbore friction, which is suitable for mixed horizontal injection and production patterns. Although this method did not take into consideration more complicated conditions, such as the mixed well pattern including multilateral-vertical wells, it provides us with an effective approach to solve these problems.
The objective of this article here is to present two semi-analytical models for the productivity evaluation of infill complex structural wells in mixed well patterns on the basis of the mirror-image method and source–sink theory. The first model is suitable for the infill horizontal well in horizontal-vertical well pattern, and the other for the infill multilateral well in multilateral-vertical well pattern. Then, the two models are applied to the study the seepage mechanism in terms of the bottom hole flowing pressure and the distribution of wellbore flow and radial flow along the horizontal segment. The main feature of this methodology is that the models take into account the interference of other wells, the randomicity of well patterns, and the pressure drawdown along the horizontal laterals. The application in the optimization of infill well location also indicates the significant practical value of the proposed models.

2. Productivity Model

2.1. Productivity Model of Mixed Horizontal-Vertical Well Pattern

2.1.1. Reservoir Flow Model

As is shown in Figure 1, in the deployment of mixed horizontal-vertical well pattern, we assume that the number of vertical producers including P 1 z , P 2 z , , P m z is mz, the number of horizontal producer including P 1 s , P 2 s , , P m s is ms, the number of vertical injectors including I 1 z , I 2 z , , I n z is nz, the number of horizontal injectors including I 1 s , I 2 s , , I n s is ns, and P n e w is the infill horizontal producer. The top and bottom boundaries of the reservoir are both closed, and the surrounding area is infinite. In order to analyze the well performance, the vertical interval and horizontal interval are both divided into N segments (Figure 2) [2,36].
The coordinates of an arbitrary point M in the vertical segment k (1 ≤ kN) of the uth vertical producer can be expressed as follows:
{ x z p ( u , k , t ) = x z p , y z p ( u , k , t ) = y z p , z z p ( u , k , t ) = z z p + ( k + t 1 ) L / N ,         0 t 1 .
The coordinates of an arbitrary point M in the horizontal segment k (1 ≤ kN) of the uth horizontal producer can be expressed as follows:
{ x s p ( u , k , t ) = x s p + ( k + t 1 ) L / N , y s p ( u , k , t ) = y s p , z s p ( u , k , t ) = z s p ,        0 t 1 .
The coordinates of an arbitrary point M in the vertical segment k (1 ≤ kN) of the vth vertical injector can be expressed as follows:
{ x z i ( v , k , t ) = x z i , y z i ( v , k , t ) = y z i , z z i ( v , k , t ) = z z i + ( k + t 1 ) L / N ,      0 t 1 .
The coordinates of an arbitrary point M in the horizontal segment k (1 ≤ kN) of the vth horizontal injector can be expressed as follows:
{ x s i ( v , k , t ) = x s i + ( k + t 1 ) L / N , y s i ( v , k , t ) = y s i , z s i ( v , k , t ) = z s i ,           0 t 1 .
The coordinates of an arbitrary point M in the horizontal segment k (1 ≤ kN) of the infill horizontal producer can be expressed as follows [37]:
{ x ( n e w , k , t ) = x n e w p + ( k + t 1 ) L / N , y ( n e w , k , t ) = y n e w p , z ( n e w , k , t ) = z n e w p ,          0 t 1 ,
where x z p ( u , k , t ) , y z p ( u , k , t ) , z z p ( u , k , t ) are the coordinates of an arbitrary point in the segment k of the uth vertical producer; x z p , y z p , z z p are the coordinates of the left end (well heel) of the uth vertical producer; x s p ( u , k , t ) , y s p ( u , k , t ) , z s p ( u , k , t ) are the coordinates of an arbitrary point in the segment k of the uth horizontal producer; x s p , y s p , z s p are the coordinates of the left end (well heel) of the uth horizontal producer; x z i ( v , k , t ) , y z i ( v , k , t ) , z z i ( v , k , t ) are the coordinates of an arbitrary point in the segment k of the vth vertical injector; x z i , y z i , z z i are the coordinates of the left end (well heel) of the vth vertical injector; x s i ( v , k , t ) , y s i ( v , k , t ) , z s i ( v , k , t ) are the coordinates of an arbitrary point in the segment k of the vth horizontal injector; x s i , y s i , z s i are the coordinates of the left end (well heel) of the vth horizontal injector; x ( n e w , k , t ) , y ( n e w , k , t ) , z ( n e w , k , t ) are the coordinates of an arbitrary point in the segment k of the infill horizontal producer; x n e w p , y n e w p , z n e w p are the coordinates of the left end (well heel) of the infill horizontal producer. L is the length of the entire horizontal section, m; N is the total number of the segments, dimensionless.
The mirror-image method and source–sink theory have been widely applied to the boundary effect reduction. As is shown in Figure 3a, assuming the constant pressure boundary is a mirror, we can project the producer (we call sink in the method) into an image injector (we call source in the method) at the symmetric coordinate to counteract the constant pressure boundary effect. As is shown in Figure 3b, in the same way, we can project the producer into an image producer at the symmetric coordinate to counteract the closed boundary effect [2].
On the basis of the method of images and the principle of superposition, we take the closed boundaries of the top and the bottom of the reservoir as mirrors. Thus, it is transformed into the problem of infinite well rows for horizontal producers, which is easy to solve. The potential of the producing segments of the infill horizontal well at any point M (x, y, z) in the infinite formation is:
Φ ( x , y , z ) = k = 1 N Φ k ( x , y , z ) = 1 4 π Δ L k k = 1 N [ q r ( k ) φ k ( x , y , z ) ] + C .
With:
φ i ( x , y , z ) = n = + { ξ k [ x ( n e w , k , t ) , y ( n e w , k , t ) , 2 n h + z ( n e w , k , t ) ] + ξ k [ x ( n e w , k , t ) , y ( n e w , k , t ) , 2 n h z ( n e w , k , t ) ] + u = 1 m z q u , z p ( k ) q r ( k ) ξ k [ x z p ( u , k , t ) , y z p ( u , k , t ) , 2 n h + z z p ( u , k , t ) ] + u = 1 m z q u , z p ( k ) q r ( k ) ξ k [ x z p ( u , k , t ) , y z p ( u , k , t ) , 2 n h z z p ( u , k , t ) ] + u = 1 m s q u , s p ( k ) q r ( k ) ξ k [ x s p ( u , k , t ) , y s p ( u , k , t ) , 2 n h + z s p ( u , k , t ) ] + u = 1 m s q u , s p ( k ) q r ( k ) ξ k [ x s p ( u , k , t ) , y s p ( u , k , t ) , 2 n h z s p ( u , k , t ) ] v = 1 n z q v , z i ( k ) q r ( k ) ξ k [ x z i ( v , k , t ) , y z i ( v , k , t ) , 2 n h + z z i ( v , k , t ) ] v = 1 n z q v , z i ( k ) q r ( k ) ξ k [ x z i ( v , k , t ) , y z i ( v , k , t ) , 2 n h z z i ( v , k , t ) ] v = 1 n s q v , s i ( k ) q r ( k ) ξ k [ x s i ( v , k , t ) , y s i ( v , k , t ) , 2 n h + z s i ( v , k , t ) ] v = 1 n s q v , s i ( k ) q r ( k ) ξ k [ x s i ( v , k , t ) , y s i ( v , k , t ) , 2 n h z s i ( v , k , t ) ] } ,      n = 0 , ± 1 , ± 2 ,
where q u , z p ( k ) Q u , z p N ,   q v , z i ( k ) Q v , z i N ,   q u , s p ( k ) Q u , s p N ,   q v , s i ( k ) Q v , s i N ; Q u , z p is the production of the uth vertical producer, m3/d; Q v , z i is the injection rate of the vth vertical injector; Q u , s p is the production of the uth horizontal producer, m3/d; Q v , s i is the injection rate of the vth horizontal injector, m3/d.
For:
ξ k [ x ( n e w , k , t ) , y ( n e w , k , t ) , z ( n e w , k , t ) ] = ln r 1 ( n e w , k ) + r 2 ( n e w , k ) + L N r 1 ( n e w , k ) + r 2 ( n e w , k ) L N ,
r 1 ( n e w , k ) = [ x ( n e w , k , t = 0 ) x ] 2 + [ y ( n e w , k , t = 0 ) y ] 2 + [ z ( n e w , k , t = 0 ) z ] 2 , r 2 ( n e w , k ) = [ x ( n e w , k , t = 1 ) x ] 2 + [ y ( n e w , k , t = 1 ) y ] 2 + [ z ( n e w , k , t = 1 ) z ] 2 .
In the same way, we can obtain the expressions of ξ k [ x ( u , k , t ) , y ( u , k , t ) , z ( u , k , t ) ] and ξ k [ x ( v , k , t ) , y ( v , k , t ) , z ( v , k , t ) ] .
Once the distributions of the flow rates of q r ( k ) ,   q u , z p ( k ) ,   q v , z i ( k ) ,   q u , s p ( k ) ,   q v , s i ( k ) are known, all terms of flow pressure p w f ( k ) of the infill horizontal producer can be calculated. Hence Equation (6) is simplified to a linear equation system, and can be arranged in the form as follows:
[ φ 11 φ 12 φ 13 φ 1 N φ 21 φ 22 φ 23 φ 2 N φ 31 φ 32 φ 33 φ 3 N φ N 1 φ N 2 φ N 3 φ N N ] [ q r ( 1 ) q r ( 2 ) q r ( 3 ) q r ( N ) ] = 4 π K Δ L μ o [ p e p w f ( 1 ) p e p w f ( 2 ) p e p w f ( 3 ) p e p w f ( N ) ] ,
where φij is the value of φj at the midpoint of the segment i of the infill horizontal producer, dimensionless; qr(i) is the radial flow rate of the segment i of the infill horizontal producer, m3/d; pe is the reservoir pressure, MPa; pwf(i) is the flowing pressure of the segment i of the infill horizontal producer, MPa; K is the reservoir permeability, 10−3 µm2; µo is the oil viscosity, mPa·s; h is the net pay thickness, m.
The wellbore flow rate along the horizontal well can be expressed as follows:
q l ( k ) = j = k N q r ( j ) ,
where ql(k) is the wellbore flow rate of the segment k of the infill horizontal producer, m3/d.

2.1.2. Wellbore Flow Model

The wellbore pressures of each segment along the horizontal segment are not independent of each other, instead, they are related to each other via wellbore hydraulics. More specifically, the pressure difference between two adjacent segment midpoints is dependent on the radial flow into the two segments, the local pressure, and the fluid property [2,36].
The pressure drop along the horizontal well is caused by the wall friction and the fluid acceleration [3,11,12,38], and can be calculated by the correlation of pressure drop in the horizontal interval under the condition of open-hole completion proposed by Liu et al. [12]:
Δ p w f ( i ) = L N { 2 f ρ π 2 D 5 [ 2 q l ( i ) q r ( i ) L N ] 2 + 16 ρ q r ( i ) π 2 D 4 [ 2 q l ( i ) q r ( i ) L N ] } , p w f ( i ) = p w f ( i 1 ) + 0.5 ( Δ p w f ( i 1 ) + Δ p w f ( i ) ) ,       ( 1 i N + 1 )
where pwf is the flowing pressure of the well heel, pwf(0) = pwf, Δpwf(0) = Δpwf(N + 1) = 0, MPa; ρ is the fluid density, g/cm3; f is the wall friction factor, and can be calculated as follows:
f = 64 N Re ,        N Re 2100 1 f = 1.14 2 lg ( τ D + 21.25 N Re 0.9 ) ,     N Re 2100
By combing Equations (9) and (10), the following system of equations can be obtained:
f [ q r ( i ) , p w f ( i ) ] = 0 .
As expected, the total number of Equations (8) and (11) is equal to 2N, which is the same as that of unknowns. Thus, the model has a unique solution.

2.1.3. Solution Procedure

The iterative algorithm is an efficient way to solve this problem, and the solution procedure is shown in Figure 4.

2.2. Productivity Model of Mixed Multilateral-Vertical Well Pattern

2.2.1. Reservoir Flow Model

As is shown in Figure 5, in the deployment of mixed multilateral-vertical well pattern, we assume that the number of vertical producers including P 1 z , P 2 z , , P m z is mz, the number of vertical injectors including I 1 z , I 2 z , , I n z is nz, and P n e w is the infill multilateral well. The top and bottom boundaries of the reservoir are both closed, and the surrounding area is infinite. In order to analyze the well performance, the vertical interval is divided into N segments.
The coordinates of an arbitrary point M in the vertical segment k (1 ≤ kN) of the uth vertical producer can be expressed as follows:
{ x z p ( u , k , t ) = x z p , y z p ( u , k , t ) = y z p , z z p ( u , k , t ) = z z p + ( k + t 1 ) L / N ,        0 t 1 .
The coordinates of an arbitrary point M in the vertical segment k (1 ≤ kN) of the vth vertical injector can be expressed as follows:
{ x z i ( v , k , t ) = x z i , y z i ( v , k , t ) = y z i , z z i ( v , k , t ) = z z i + ( k + t 1 ) L / N ,      0 t 1
where x z p ( u , k , t ) , y z p ( u , k , t ) , z z p ( u , k , t ) are the coordinates of an arbitrary point in the segment k of the uth vertical producer; x z p , y z p , z z p are the coordinates of the left end (well heel) of the uth vertical producer; x z i ( v , k , t ) , y z i ( v , k , t ) , z z i ( v , k , t ) are the coordinates of an arbitrary point in the segment k of the vth vertical injector; x z i , y z i , z z i are the coordinates of the left end (well heel) of the vth vertical injector.
Figure 6 is the three-dimensional (3D) schematic of the infill multilateral well, we assume that each branch is symmetrically distributed on the same horizontal plane, the length of which is identical. The xoy coordinate system is established by taking the subpoint of the main hole in the xoy plane as the origin of coordinates, and the direction of the main hole as the z axis. Each branch is distributed counterclockwise around the z axis with the x axis as the starting direction.
M0(x, y, z0) is the bottom-hole coordinate of the main wellbore, and Mi(xi,0, yi.0, zi,0) is the starting coordinate of the ith branch with the length of Li (1 ≤ iM). M is the number of branches. The ith branch of the infill multilateral well is divided into Ni (1 ≤ iM) segments, thus the infinitesimal length of the ith branch is ΔLi = Li/Ni.
Accordingly, the coordinates of an arbitrary point M in the horizontal segment i (1 ≤ iM) of the infill multilateral well can be expressed as follows:
{ x ( i , j ) = x i , 0 + Δ L i [ k = 1 j 1 ( sin θ i , k cos α i , k ) + t sin θ i , j cos α i , k ] , y ( i , j ) = y i , 0 + Δ L i [ k = 1 j 1 ( sin θ i , k sin α i , k ) + t sin θ i , j sin α i , k ] , z ( i , j ) = z i , 0 + Δ L i [ k = 1 j 1 cos θ i , k + t cos θ i , j ] ,     ( 1 i M , 0 j N i , 0 t 1 )
where x(i, j), y(i, j), z(i, j) are the coordinates of an arbitrary point in the ith branch of the infill multilateral producer; xi,0, yi,0, zi,0 are the coordinates of the left end (well heel) in the ith branch of the infill multilateral producer; ɵ is the deviation angle of the ith segment of the infill multilateral producer; α is the azimuth angle of the ith segment of the infill multilateral producer.
As mentioned above, we can also transform this problem into the one of infinite well rows by using the method of images and the principle of superposition. Thus the potential of the producing segments of the infill multilateral well at any point M(x, y, z) in the infinite formation is:
f ( x , y , z ) = i = 1 M f i ( x , y , z ) = 1 4 π i = 1 M j = 1 N j [ q r ( i , j ) φ i , j ( x , y , z ) ] + C .
With:
φ i , j ( x , y , z ) = n = + [ 1 Δ L i ξ i , j ( x i , j , y i , j , 2 n h + z i , j ) + 1 Δ L i ξ i , j ( x i , j , y i , j , 2 n h z i , j ) + 1 Δ L u = 1 m z q u , z p ( k ) q r ( i , j ) ξ k [ x z p ( u , k , t ) , y z p ( u , k , t ) , 2 n h + z z p ( u , k , t ) ] + 1 Δ L u = 1 m z q u , z p ( k ) q r ( i , j ) ξ k [ x z p ( u , k , t ) , y z p ( u , k , t ) , 2 n h z z p ( u , k , t ) ] 1 Δ L v = 1 n z q v , z i ( k ) q r ( i , j ) ξ k [ x z i ( v , k , t ) , y z i ( v , k , t ) , 2 n h + z z i ( v , k , t ) ] 1 Δ L v = 1 n z q v , z i ( k ) q r ( i , j ) ξ k [ x z i ( v , k , t ) , y z i ( v , k , t ) , 2 n h z z i ( v , k , t ) ]
where q u , z p ( k ) Q u , z p N , q v , z i ( k ) Q v , z i N ; Q u , z p is the production of the uth vertical producer, m3/d; Q v , z i is the injection rate of the vth vertical injector, m3/d.
Once the distributions of the flow rates of q r ( i , j ) ,   q u , z p ( k ) ,   q v , z i ( k ) are known, all terms of flow pressure p w f ( i , j ) of the infill multilateral producer can be calculated. Hence Equation (15) is simplified to a linear equation system, and can be arranged in the form as follows:
A [ q r ( 1 , 1 ) q r ( 1 , 2 ) q r ( 1 , N 1 ) q r ( M , N M ) ] = 4 π k h k v μ o [ P e P w f ( 1 , 1 ) P e P w f ( 1 , 2 ) P e P w f ( 1 , N 1 ) P e P w f ( 1 , N M ) ] .
With:
A = [ ϕ ( 1 , 1 ) ( 1 , 1 ) ϕ e ( 1 , 1 ) ϕ ( 1 , 1 ) ( 1 , N 1 ) ϕ e ( 1 , N ) ϕ ( 1 , 1 ) ( 2 , 1 ) ϕ e ( 2 , 1 ) ϕ ( 1 , 1 ) ( M , N M ) ϕ e ( M , N M ) ϕ ( 1 , N 1 ) ( 1 , 1 ) ϕ e ( 1 , 1 ) ϕ ( 1 , N 1 ) ( 1 , N 1 ) ϕ e ( 1 , N ) ϕ ( 1 , N 1 ) ( 2 , 1 ) ϕ e ( 2 , 1 ) ϕ ( 1 , N 1 ) ( M , N M ) ϕ e ( M , N M ) ϕ ( 2 , 1 ) ( 1 , 1 ) ϕ e ( 1 , 1 ) ϕ ( 2 , 1 ) ( 1 , N 1 ) ϕ e ( 1 , N ) ϕ ( 2 , 1 ) ( 2 , 1 ) ϕ e ( 2 , 1 ) ϕ ( 2 , 1 ) ( M , N M ) ϕ e ( M , N M ) ϕ ( M , N M ) ( 1 , 1 ) ϕ e ( 1 , 1 ) ϕ ( M , N M ) ( 1 , N 1 ) ϕ e ( 1 , N ) ϕ ( M , N M ) ( 2 , 1 ) ϕ e ( 2 , 1 ) ϕ ( M , N M ) ( M , N M ) ϕ e ( M , N M ) ]
where qr(i, j) is the radial flow rate of the ith branch of the infill multilateral producer, 1 ≤ iM, 1 ≤ jNi, m3/d; Pe is the reservoir pressure, MPa; Pwf(i, j) is the flow pressure of the ith branch of the infill multilateral producer, 1 ≤ iM, 1 ≤ jNi, MPa; kh is the reservoir horizontal permeability, 10−3 µm2; kv is the reservoir vertical permeability, 10−3 µm2; µo is the oil viscosity, mPa·s.

2.2.2. Wellbore Flow Model

1. Pressure drop model of the horizontal lateral
The pressure drop along the horizontal interval is caused by the wall friction and the fluid acceleration. According to Equation (11), the pressure drop of the ith branch of the multilateral well along the horizontal interval can be calculated by:
F i , j [ q r ( i , j ) , p w f ( i , j ) ] = 0 .
2. Pressure drop model of the deviated segment
The branch holes enter the reservoir through the deviated segment. Although each deviated segment does not contribute to the production, the pressure drop in the deviated segments cannot be ignored. According to fluid seepage theory in the elbow, the pressure drop can be calculated by the following formula:
Δ p s = C w 4 × 10 6 ρ λ c R c Q W 2 π d w 5 + ρ g R c 10 6 ,
where ΔPs the pressure drop along the deviated segment, MPa; Rc is curvature radius of the deviated segment, m; ρ is the fluid density, kg/m3; dw is the wellbore diameter of the deviated segment, m; Cw is the similarity factor of the wellbore and deviated segment, dimensionless; λc is the friction factor of the deviated segment, dimensionless.
3. Pressure drop model of the vertical segment
According to the law of conservation of energy, the equation of pressure gradient in the pipe with inclination can be expressed as follows:
d P d z = ρ g sin θ + ρ v d v d z + f ρ v 2 2 d ,
where dP/dz is the pressure loss per unit length, MPa/m; ρgsinθ is the pressure drop caused by the fluid gravity, MPa/m; ρv × dv/dz is the pressure drop caused by the fluid acceleration, MPa/m; f × ρv2/2d is the pressure drop caused by the friction, MPa/m; ρ is the fluid density, kg/m3; d is the pipe diameter, m. The friction factor f can be calculated by:
f = 64 N Re , N Re 2100 1 f = 1.14 2 lg ( τ d + 21.25 N Re 0.9 ) , N Re 2100
where NRe is the Reynolds number, dimensionless.
4. Internal pressure drop model of the branch hole
The flowing pressure at the heel of the ith branch can be calculated by simultaneous Equations (19) and (20):
P w f ( i ) = P m b Δ P v i Δ P w i ( 1 i M ) ,
where Pmb is the flowing pressure of the main wellbore, MPa; ΔPvi is the pressure drop of the vertical segment of the ith branch, MPa; ΔPwi is the pressure drop of the deviated segment of the ith branch, MPa.
By combining Equations (20) and (24), the following system of equations can be obtained:
F i , j [ q r ( i , j ) , p w f ( i , j ) ] = 0 , ( 1 i M , 1 j N i ) .
As expected, the total number of equations of Equations (17) and (22) is equal to 2 ( M × N i ) , which is the same as that of unknowns. Therefore, the solution of the coupling model is unique.

2.2.3. Solution Procedure

The coupling model can also be solved by using the iterative algorithm, and the solution procedure is shown in Figure 7.

3. Results and Discussion

3.1. Model Validation

The semi-analytical models can be applied to the productivity prediction of the infill horizontal well and the multilateral well in mixed well patterns. In order to verify the reliability of the models, two cases are considered, and the results calculated by the semi-analytical models are compared with those calculated by the Eclipse software (Schlumberger). In the first case, we selected a typical mixed well pattern of horizontal-vertical wells from a thin carbonate reservoir in the Middle East. There are three vertical producers (P11, P12, and P13), one horizontal producer (P14), one vertical injector (I11), one horizontal injector (I12), and P1new is the infill horizontal producer. In the other case, the reservoir is a thick, multi-layer reservoir in the Middle East, multilateral wells are adopted to improve the degree of the reservoir development. The mixed well pattern of multilateral-vertical well includes two vertical producers (P21 and P22), one vertical injector (I21), and P2new is the infill multilateral producer with two branch holes.
The productivity models can be solved by VB modular programming of the object-oriented technology. The reservoir properties and wellbore geometry are listed in Table 1 and Table 2. By inputting the same model parameters, the productivity evaluation of the two wells were conducted respectively with the semi-analytical models and the Eclipse simulator. The calculation results are shown in Table 3.
Table 3 indicates that the semi-analytical models have a high accuracy with the relative error of less than 15%. The results calculated by the models are consistent with those by the Eclipse numerical simulator, which verify the reliability of the productivity models. By taking into account the actual conditions, namely the interference of other wells in the well pattern, the coupling between the reservoir flow and wellbore flow, and the pressure drop along the horizontal lateral, the accuracy of the productivity prediction by the semi-analytical models is greatly improved.

3.2. Model Application

3.2.1. Study on Seepage Mechanism of Horizontal Well

In the first case of model validation, we can further analyze the bottom hole flowing pressure, the distribution of wellbore flow rate and radial flow rate along the horizontal segment. Figure 8 and Figure 9 show that the bottom hole flowing pressure decreases logarithmically while the wellbore flow rate increases monotonically from the toe to the heel of the horizontal well. The reservoir flow coupled with variable mass wellbore flow results in this near wellbore dynamics. On one hand, the mass flow rate of fluid from the toe to the heel increases gradually. In this case, the fluid velocity along the main flow direction also increases with an accelerating pressure drop. The radial inflow of reservoir fluid along the horizontal wellbore disturbs the boundary layer of the main stream tube and affects its velocity profile, thus altering the wall friction resistance determined by the velocity distribution. These factors lead to the bottom hole flowing pressure decreasing in the form of logarithmic function. On the other hand, the radial inflow performance affects the distribution of pressure and pressure drop, and vice versa. Due to the pseudo-hemispherical flow at each endpoint (well heel and well toe) and the pseudo-linear flow at the center of the horizontal segment, the drainage area at each endpoint is relatively larger than that at the center. The radial flow rate at each endpoint of the horizontal segment is considerably greater than that at the center, and decreases quickly toward the intermediate section, which generally presents a U-shape distribution.

3.2.2. Optimization of Infilling Well Location

In the intermediary and later phase of oilfield development, infill wells are usually needed to improve the well pattern and the effect of tapping potential of remaining oil. Therefore, the key is the optimization of infill well location. The productivity models provide a shortcut for this target. Based on the understanding of geological conditions and development status, the productivity of infill well can be accurately predicted, and in combination with the findings of remaining oil distribution, the infill well location could be optimized.
Taking the mixed five-spot well pattern of vertical injectors and horizontal producer as an example, we conducted a study on the productivity variation of horizontal well at different locations in the well pattern (Figure 10). The basic reservoir parameters refer to the first case of model validation. Figure 11 shows that the lateral shifting of horizontal well has some effects on the productivity, and the productivity at the center of the well pattern is relatively higher than that at other locations. However, this does not mean that the central location is the optimum selection for the horizontal well. In the adjustment of well pattern, geological condition and development factors, such as the distribution of flow field and remaining oil, should also be considered.
We made further study on the flow field distribution in the mixed well pattern so as to provide references for the optimization of infilling well location. Figure 12 and Figure 13 indicate that the equipotential lines and stream lines are symmetrically distributed along the x-axis as the horizontal well moves laterally. When the horizontal well is located at the center of the well pattern, both the equipotential lines and stream lines present a uniform distribution; when the horizontal well moves toward the right side, both of the two types of the lines on the right side become dense, while those on the left side are sparsely distributed, and vice versa. The more the horizontal well location changes, the more obvious this phenomenon is. The flow pattern almost presents a linear flow as the horizontal well locates at the line of the two vertical wells on the same side, which indicates a relatively higher swept efficiency. Therefore, in combination with the findings of productivity variation and flow fluid distribution, the infilling well location could be optimized. This improves the effect of well pattern adjustment considerably.

3.3. Discussion

In this study, the semi-analytical models for productivity evaluation of infill complex structural wells are built with the method of images and source–sink theory. Compared with the existing models, the advantages of the semi-analytical model lie in the following three aspects. Firstly, the models are suitable for mixed well patterns of vertical wells and complex structure wells; Secondly, both the interference of other wells in the well pattern and the pressure drawdown along the horizontal segment are taken into account, which makes the models more reliable. In addition, the models provide some insight into the seepage mechanism of complex structural wells such as the distribution of wellbore flow and radial flow along the horizontal segment. However, the semi-analytical models also have some limitations, the basic assumptions of the models are steady-state flow and closed top and bottom boundaries, and the productivity evaluation of different reservoir boundaries under the condition of unsteady-state flow still needs to be further studied.

4. Conclusions

In this paper, on the basis of the mirror-image method and the source–sink theory, two semi-analytical models are established to predict the productivity of infill horizontal wells and multilateral wells in mixed well patterns respectively, in which the interference of other wells, the randomicity of well pattern, and the pressure drawdown along the horizontal lateral are taken into account.
The application of the productivity models verifies the reliability and practicability. The results calculated by the models are consistent with those by the Eclipse numerical simulator, and the semi-analytical models have a high accuracy with the relative error of less than 15%.
The results indicate that the bottom hole flowing pressure decreases logarithmically while the wellbore flow rate increases monotonically from the toe to the heel of the horizontal segment. Due to the pseudo-hemispherical flow at each endpoint and the pseudo-linear flow at the center of the horizontal segment, the drainage area at each endpoint is relatively larger than that at the center. The radial flow rate at each endpoint of the horizontal segment is considerably greater than that at the center, which generally presents a U-shape distribution.
This study also proposes a practical and efficient approach for the study on the horizontal seepage mechanism, and the optimization of infilling well locations.

Author Contributions

Conceptualization, L.S.; Data curation, L.S.; Formal analysis, L.S.; Funding acquisition, Y.L.; Methodology, L.S.; Project administration, B.L.; Supervision, Y.L.; Validation, L.S.; Writing – review & editing, L.S., Y.L.

Funding

This research was founded by the National Science and Technology Major Project of China, grant number 2016ZX05015-002, and the National Natural Science Foundation of China, grant number 2017ZX05030-001.

Acknowledgments

The authors are grateful for the technical and financial support from the Research Institute of Exploration and Development, PetroChina, Beijing 100083, China.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Adesina, F.; Paul, A.; Oyinkepreye, O.; Adebowale, O. An improved model for estimating productivity of horizontal drain hole. In Proceedings of the SPE Nigeria Annual International Conference and Exhibition, Lagos, Nigeria, 2 August 2016. [Google Scholar] [CrossRef]
  2. Ye, S.J. A productivity evaluation model and its application for an infill horizontal well in different types of reservoirs. Pet. Sci. Technol. 2012, 30, 1677–1691. [Google Scholar]
  3. Liangbiao, O.; Arbabi, S.; Aziz, K. General wellbore flow model for horizontal, vertical and slanted well completions. SPE J. 1998, 3, 124–133. [Google Scholar]
  4. Merkulov, V.P. The flow of slanted and horizontal well. Neft. Khoz 1958, 6, 51–56. [Google Scholar]
  5. Borisov, J.P. Oil Production Using Horizontal and Multiple Deviation Wells; Nedra: Moscow, Russia, 1964. [Google Scholar]
  6. Giger, F.M. Horizontal wells production techniques in heterogeneous reservoirs. In Proceedings of the Middle East Oil Technical Conference and Exhibition, Manama, Bahrain, 14–17 March 1983. [Google Scholar] [CrossRef]
  7. Joshi, S.D. Augmentation of well productivity using slant and horizontal wells. In Proceedings of the SPE Annual Technical Conference and Exhibition, New Orleans, Louisiana, 5 October 1986. [Google Scholar] [CrossRef]
  8. Raghavan, R.; Joshi, S.D. Productivity of multiple drainholes or fractured horizontal wells. SPE Format. Eval. 1993, 8, 11–16. [Google Scholar] [CrossRef]
  9. Thomas, L.K.; Todd, B.J.; Evans, C.E.; Pierson, R.G. Horizontal well IPR calculations. SPE Reserv. Eval. Eng. 1998, 1, 392–399. [Google Scholar] [CrossRef]
  10. Zhang, W.; Han, D. 3D potential distribution and precise productivity equation of horizontal well. Pet. Explor. Dev. 1999, 26, 49–52. [Google Scholar]
  11. Dikken, B.J. Pressure drop in horizontal wells and its effect on production performance. J. Pet. Technol. 1990, 42, 1426–1433. [Google Scholar] [CrossRef]
  12. Liu, X.; Zhang, Z.S.; Liu, X.E.; Guo, S. A model to calculate pressure drops of horizontal wellbore variable mass flow coupled with flow in a reservoir. J. Southwest Pet. Inst. 2000, 22, 36–39. [Google Scholar]
  13. Wang, R.H.; Zhang, Y.Z. A segmentally numerical calculation method for estimating the productivity of perforated horizontal wells. Pet. Explor. Dev. 2006, 33, 630. [Google Scholar]
  14. Li, X.P.; Guo, C.Z.; Jiang, Z.X.; Liu, X.E.; Guo, S.P. The model coupling fluid flow in the reservoir with flow in the horizontal wellbore. Acta Pet. Sin. 1999, 3, 82–86. [Google Scholar]
  15. Butler, R.M. Discussion of augmentation of well productivity with slant and horizontal wells. Author’s reply. J. Pet. Technol. 1992, 44, 942–943. [Google Scholar]
  16. Dang, L. Analysis of productivity formulae of horizontal well. Pet. Explor. Dev. 1997, 5, 21. [Google Scholar]
  17. Jiang, H.; Ye, S.; Lei, Z.; Wang, X.; Zhu, G.; Chen, M. The productivity evaluation model and its application for finite conductivity horizontal wells in fault block reservoirs. Pet. Sci. 2010, 7, 530–535. [Google Scholar] [CrossRef] [Green Version]
  18. Huang, S.J.; Cheng, L.S.; Zhao, F.L. The productivity evaluation model of the stepped horizontal well in thin interbeded reservoirs. J. Southwest Pet. Univ. 2007, 3, 16. [Google Scholar]
  19. Борисов, Ю.П.; Табаков, В.П. Расчет взаимoдействия батарей наклoнных и мнoгoзабoйных скважин в слoистoм пласте; НТС пo дoбыче нефти: Оренбургская oбл, г. Оренбург, ул., Рoссия, 1961. [Google Scholar]
  20. Wang, W.H.; Li, D. Productivity study on branch horizontal wells. Oil Drill. Prod. Technol. 1997, 4, 12. [Google Scholar]
  21. Li, C.L. Derivation of productivity formulae of a fishbone well. J. Southwest Pet. Inst. 2005, 27, 36. [Google Scholar]
  22. Cheng, L.S.; Li, C.L.; Lang, Z.X.; Zhang, L.H. The productivity study of branch a horizontal well with multiple branched wells. Acta Pet. Sin. 1995, 16, 49–55. [Google Scholar]
  23. Zhao, L.X.; Jiang, M.H.; Zhao, X.F. Research on deliverability relationship of complicated horizontal well. J. Univ. Pet. China (Ed. Nat. Sci.) 2006, 30, 77–80. [Google Scholar]
  24. Johansen, T.E.; James, L.; Cao, J. Analytical coupled axial and radial productivity model for steady-state flow in horizontal wells. IJPE 2015, 1, 290. [Google Scholar] [CrossRef] [Green Version]
  25. Wang, H.; Guo, J.; Zhang, L. A semi-analytical model for multilateral horizontal wells in low-permeability naturally fractured reservoirs. J. Pet. Sci. Eng. 2017, 149, 564–578. [Google Scholar] [CrossRef]
  26. Simonov, M.V.; Akhmetov, A.V.; Roshchektaev, A.P. Semi-analytical model of transient fluid flow to multilateral well. In Proceedings of the SPE Annual Caspian Technical Conference and Exhibition, Baku, Azerbaijan, 1–3 November 2017. [Google Scholar] [CrossRef]
  27. Hassan, A.; Abdulraheem, A.; Elkatatny, S.; Ahmed, M. New approach to quantify productivity of fishbone multilateral well. In Proceedings of the SPE Annual Technical Conference and Exhibition, San Antonio, TX, USA, 9–11 October 2017. [Google Scholar] [CrossRef]
  28. Liu, G.; Meng, Z.; Cui, Y.; Wang, L.; Liang, C.; Yang, S. A semi-analytical methodology for multiwell productivity index of well-industry-production-scheme in tight oil reservoirs. Energies 2018, 11, 1054. [Google Scholar] [CrossRef]
  29. Vodorezov, D.D. Estimation of horizontal-well productivity loss caused by formation damage on the basis of numerical modeling and laboratory-testing data. SPE J. 2018. [Google Scholar] [CrossRef]
  30. Al-Rbeawi, S.; Artun, E. Fishbone type horizontal wellbore completion: A study for pressure behavior, flow regimes, and productivity index. J. Pet. Sci. Eng. 2019. [Google Scholar] [CrossRef]
  31. Suprunowicz, R.; Butler, R.M. The productivity and optimum pattern shape for horizontal wells arranged in staggered rectangular arrays. J. Can. Pet. Technol. 1992, 31. [Google Scholar] [CrossRef]
  32. Chunlan, L.; Lingsong, C.; Lihua, Z.; Zhaoxin, L. The Study on Productivity of Horizontal Well 9-spot Patterns. J. Southwest Pet. Inst. 1998, 20, 56–58. [Google Scholar]
  33. Cheng, L.S.; Zheng, J.Q.; Li, C.L. Productivity study of horizontal wells pattern. Oil Drill. Prod. Technol. 2002, 24, 39–41. [Google Scholar]
  34. Liu, Y.T.; Zhang, J.C. Stable permeating flow and productivity analysis for anisotropic reservoir in horizontal well networks. Pet. Explor. Dev. 2004, 31, 94–96. [Google Scholar]
  35. Yin, G.F.; Xu, H.M.; Ye, S.J.; Li, Y.R.; Qiu, J.P. Productivity evaluation model for infill horizontal well in horizontal injection and production pattern. J. China Univ. Pet. (Ed. Nat. Sci.) 2011, 4, 18. [Google Scholar]
  36. Ye, S.J.; Jiang, H.Q.; Li, J.J. Productivity calculation of infill horizontal wells in mixed well pattern. Chin. J. Comput. Phys. 2011, 28, 693–697. [Google Scholar]
  37. Huang, S.J.; Cheng, L.S.; Zhao, F.L.; Li, C.L. The flow model coupling reservoir percolation and variable mass pipe flow in production section of the stepped horizontal well. J. Hydrodyn. 2005, 20, 463–471. [Google Scholar]
  38. Tabatabaei, M.; Ghalambor, A. A new method to predict performance of horizontal and multilateral wells. SPE Prod. Oper. 2011, 26, 75–87. [Google Scholar] [CrossRef]
Figure 1. The deployment of mixed horizontal-vertical well pattern.
Figure 1. The deployment of mixed horizontal-vertical well pattern.
Processes 07 00324 g001
Figure 2. Schematic of the segmented horizontal well.
Figure 2. Schematic of the segmented horizontal well.
Processes 07 00324 g002
Figure 3. The mirror-image method and source–sink theory applied in different boundaries: (a) constant pressure boundary; (b) closed boundary.
Figure 3. The mirror-image method and source–sink theory applied in different boundaries: (a) constant pressure boundary; (b) closed boundary.
Processes 07 00324 g003
Figure 4. The flow chart of solution procedure.
Figure 4. The flow chart of solution procedure.
Processes 07 00324 g004
Figure 5. The mixed well pattern of multilateral-vertical wells.
Figure 5. The mixed well pattern of multilateral-vertical wells.
Processes 07 00324 g005
Figure 6. 3D schematic of infill multilateral well.
Figure 6. 3D schematic of infill multilateral well.
Processes 07 00324 g006
Figure 7. The flow chart of solution procedure.
Figure 7. The flow chart of solution procedure.
Processes 07 00324 g007
Figure 8. Distribution of bottom hole flowing pressure along the horizontal lateral.
Figure 8. Distribution of bottom hole flowing pressure along the horizontal lateral.
Processes 07 00324 g008
Figure 9. Distribution of wellbore flow and radial flow along the horizontal lateral.
Figure 9. Distribution of wellbore flow and radial flow along the horizontal lateral.
Processes 07 00324 g009
Figure 10. Mixed five-spot well pattern of vertical-injectors and horizontal producer (①–⑦ show the code of different horizontal locations).
Figure 10. Mixed five-spot well pattern of vertical-injectors and horizontal producer (①–⑦ show the code of different horizontal locations).
Processes 07 00324 g010
Figure 11. Productivity variation of horizontal well in mixed well pattern.
Figure 11. Productivity variation of horizontal well in mixed well pattern.
Processes 07 00324 g011
Figure 12. Equipotential line pattern for lateral shifting of the horizontal well. ((ac) show the equipotential line pattern of different lateral location.)
Figure 12. Equipotential line pattern for lateral shifting of the horizontal well. ((ac) show the equipotential line pattern of different lateral location.)
Processes 07 00324 g012
Figure 13. Stream line pattern for lateral shifting of the horizontal well. ((ac) show the stream line pattern of different lateral location.)
Figure 13. Stream line pattern for lateral shifting of the horizontal well. ((ac) show the stream line pattern of different lateral location.)
Processes 07 00324 g013
Table 1. Productivity calculation parameters of P1new.
Table 1. Productivity calculation parameters of P1new.
Parameters (Unit)Value
Net thickness (m)6.0
Porosity (%)15.3
Permeability (10−3 µm2)88.5
Reservoir pressure (MPa)38.6
Oil viscosity (mPa·s)8.1
Horizontal length of P1new (m)335
Wellbore radius of P1new (m)0.1
Bottom hole flowing pressure (MPa)36.6
Production rate of P11 (m3/d)18.5
Production rate of P12 (m3/d)21.6
Production rate of P13 (m3/d)24.3
Production rate of P14 (m3/d)56.8
Injection rate of I11 (m3/d)68.5
Injection rate of I12 (m3/d)79.2
Table 2. Productivity calculation parameters of P2new.
Table 2. Productivity calculation parameters of P2new.
Parameters (Unit)Value
Net thickness-the 1st branch of P2new (m)5.1
Net thickness-the 2nd branch of P2new (m)9.4
Porosity (%)22.1
Ratio of vertical to horizontal permeability (-)0.1
Reservoir pressure (MPa)20.3
Oil viscosity (mPa·s)6.2
Length of the 1st branch of P2new (m)200
Length of the 2nd branch of P2new (m)140
Bottom hole flowing pressure (MPa)10.0
Production rate of P21 (m3/d)52.9
Production rate of P22 (m3/d)67.8
Injection rate of I21 (m3/d)125.5
Table 3. Results comparison of two methods.
Table 3. Results comparison of two methods.
Well NameActual Productivity
(m3/d)
The Semi-Analytical ModelThe Eclipse Simulator
Productivity
(m3/d)
Relative Error (%)Productivity
(m3/d)
Relative Error (%)
P1new25.027.39.127.28.8
P2new1st branch197.7220.711.7221.311.9
2nd branch147.1166.513.2168.014.2
Total344.7387.212.3389.312.9

Share and Cite

MDPI and ACS Style

Sun, L.; Li, B.; Li, Y. Productivity Models of Infill Complex Structural Wells in Mixed Well Patterns. Processes 2019, 7, 324. https://doi.org/10.3390/pr7060324

AMA Style

Sun L, Li B, Li Y. Productivity Models of Infill Complex Structural Wells in Mixed Well Patterns. Processes. 2019; 7(6):324. https://doi.org/10.3390/pr7060324

Chicago/Turabian Style

Sun, Liang, Baozhu Li, and Yong Li. 2019. "Productivity Models of Infill Complex Structural Wells in Mixed Well Patterns" Processes 7, no. 6: 324. https://doi.org/10.3390/pr7060324

APA Style

Sun, L., Li, B., & Li, Y. (2019). Productivity Models of Infill Complex Structural Wells in Mixed Well Patterns. Processes, 7(6), 324. https://doi.org/10.3390/pr7060324

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop