WJETWorld Journal of Engineering and Technology2331-4222Scientific Research Publishing10.4236/wjet.2019.72021WJET-92279ArticlesChemistry&Materials Science Engineering Study on Productivity Model of Herringbone-Like Laterals Wells and Optimization of Morphological Parameters Considering Threshold Pressure Gradient in Heavy Oil Reservoirs EnhuiSun1*JieTan1DongZhang1WeiWang1SongruMu1Bohai Oilfield Research Institute of CNOOC Ltd.-Tianjin Branch, Tianjin, China0103201907023023131, April 20195, May 2019 8, May 2019© Copyright 2014 by authors and Scientific Research Publishing Inc. 2014This work is licensed under the Creative Commons Attribution International License (CC BY). http://creativecommons.org/licenses/by/4.0/

Compared with conventional well, herringbone-like laterals wells can increase the area of oil release, and can reduce the number of wellhead slots of platforms, and also can greatly improve the development efficiency. Based on threshold pressure gradient in heavy oil reservoir, and the applied principle of mirror reflection and superposition, the pressure distribution equation of herringbone-like laterals wells is obtained in heavy oil reservoir. Productivity model of herringbone-like laterals wells is proposed by reservoir-wellbore steady seepage. The example shows that the productivity model is great accuracy to predict the productivity of herringbone-like laterals wells. The model is used to analyze the branching length, branching angle, branching symmetry, branching position and spacing and their effects on productivity of herringbone-like laterals wells. The principle of optimizing the well shape of herringbone-like laterals wells is proposed.

Threshold Pressure Gradient Herringbone-Like Laterals Wells Heavy Oil Reservoirs Productivity Model Optimization of Morphological Parameters1. Introduction

Compared with the onshore oilfields, the number of wellhead slots is limited in offshore oilfields. The herringbone-like laterals wells not only increase the drainage area and single well-controlled reserves, but also increase the production of oil wells; it has certain advantages in offshore oilfield development. Therefore, it is necessary to study the productivity prediction of herringbone-like laterals wells in reservoirs. In 1996, Salas [1] assumed each branch is divided into several segments and established analytical model of single-phase seepage in fishbone multi-branch horizontal wells. In 2004, Han Guoqin [2] and He Haifeng [3] established a steady seepage mathematical model of fishbone multi-branch wells considering the flow in the wellbore. Liu Xiangping, Yang Xiaosong, Zhao Guang, et al. [4] [5] [6] [7] [8] deduced the productivity calculation model of fishbone horizontal well by potential superposition principle. Based on the principle of equivalent seepage resistance and productivity formula of horizontal wells, a simple productivity formula for fishbone horizontal wells is derived by Li Chunlan [9] . Fan Yuping, Ye Shuangjiang, Zhang Shiming, et al. [10] [11] [12] used numerical simulation method to study the productivity of herringbone-like laterals wells. Because heavy oil is affected by threshold pressure gradient, the above research is not suitable for predicting the productivity of wells in heavy oil reservoirs. Considering the influence of threshold pressure gradient on heavy oil reservoir, herringbone-like laterals wells coupling productivity model for reservoir-wellbore steady seepage is proposed; the example shows that the productivity model is high accuracy in predicting the productivity of herringbone-like laterals wells, and the principle of well shape optimization for herringbone-like laterals wells is proposed.

2. Pressure Distribution Considering Threshold Pressure Gradient

In reference [13] , porous media conditions, heavy oil reservoir seepage law is non-Darcy seepage with threshold pressure gradient. When the driving pressure gradient exceeds its initial pressure gradient, heavy oil begins to flow and its seepage characteristics are shown in Figure 1. Among them, A is the minimum threshold pressure gradient, B is the average threshold pressure gradient; C is the maximum threshold pressure gradient.

The law of non-Darcy seepage can be described by the following formula:

v = − K μ ( d p d r − λ ) ,         d p d r > λ (1)

where v is seepage velocity, m3/s; K is reservoir permeability, mD; d p d r is pressure gradient, MPa/m; λ is threshold pressure gradient, MPa/m.

when Φ = K μ p , have:

v = − K μ d Φ d r + λ K μ (2)

Assuming that the formation is infinitely homogeneous and isotropic, there is a horizontal well in which the length of the horizontal well is L, the coordinates of the two ends of the well are (x1, 0, zw), (x1, 0, zw), setting the horizontal well is homogeneous line sink, the productivity of the well is Q, the productivity of per unit length L is q, selecting micro-element dx0 at x0 of horizontal well. It can be regarded as the seepage velocity of M(x0, y0, z0) at any point:

v = q 4 π r 2 (3)

A new velocity potential function considering threshold pressure gradient at point M is obtained:

d Φ = − Q 4 π L d x 0 ( x 0 − x ) 2 + ( y 0 − y ) 2 + ( z 0 − z ) 2 + λ K μ d R D + C (4)

where R D is the shortest distance between the micro-element and the moving boundary.

According to the superposition principle of potential, the velocity potential caused by the whole horizontal well is as follows:

Φ = ∫ x 1 x 2 − Q 4 π L d x 0 ( x 0 − x ) 2 + ( y 0 − y ) 2 + ( z 0 − z ) 2 + λ K μ R D + c (5)

In this paper, a horizontal well is divided into N segments by mathematical discretization method. As the productivity of the well varies little in a small section, it can be assumed that the productivity of this section is a fixed value, and the production of each section is different:

Φ = ∑ i = 0 N − 1 ( − Q 4 π L i ∫ x i x i + 1 d x 0 ( x 0 − x ) 2 + ( y 0 − y ) 2 + ( z 0 − z ) 2 ) + λ K μ R D + c ′ (6)

As the horizontal well is parallel to the X-axis and on the XOY plane, so y 0 = y , z 0 = z w , the Formula (6) is changed to:

Φ ( M ) = ∑ i = 0 N − 1 − Q 4 π L i ln ( x i + 1 − x ) + r i + 1 ( x i − x ) + r i + λ K μ R D + c ′ (7)

where r i = ( x i − x ) 2 + y 2 + ( z w − z ) 2 .

The herringbone-like laterals wells have M branching (including main branching). The herringbone-like laterals wells are divided into several segments. And it is setted the flow rate of section s of the branching t. According to the superposition principle of mirror reflection and potential, potential produced by herringbone-like laterals wells system at any point M in reservoir:

Φ t s ( M ) = ∑ k = − ∞ + ∞ ∑ t = 0 M − 1 ∑ s = 1 N [ φ ( ( x [ t ] [ s ] , y [ t ] [ s ] , 4 k h + z w ) , ( x [ t ] [ s + 1 ] , y [ t ] [ s + 1 ] , 4 k h + z w ) , α t )     + φ ( ( x [ t ] [ s ] , y [ t ] [ s ] , 4 k h + 2 h − z w ) , ( x [ t ] [ s + 1 ] , y [ t ] [ s + 1 ] , 4 k h + 2 h − z w ) , α t )     − φ ( ( x [ t ] [ s ] , y [ t ] [ s ] , 4 k h − z w ) , ( x [ t ] [ s + 1 ] , y [ t ] [ s + 1 ] , 4 k h − z w ) , α t )     − φ ( ( x [ t ] [ s ] , y [ t ] [ s ] , 4 k h − 2 h + z w ) , ( x [ t ] [ s + 1 ] , y [ t ] [ s + 1 ] , 4 k h − 2 h + z w ) , α t ) ]     + λ K μ ( R e − x ) 2 + y 2 + z 2 (8)

where:

φ ( ( x [ t ] [ s ] , y [ t ] [ s ] , 4 k h + 2 h − z w ) , ( x [ t ] [ s + 1 ] , y [ t ] [ s + 1 ] , 4 k h + 2 h − z w ) , α t ) = − q t s 4 π L t s ln r t s + r t s + 1 + L t s r t s + r t s + 1 − L t s + c t s

r t s = ( x t s − x ) 2 + ( y t s − y ) 2 + ( 4 k h + 2 h − z w − z ) 2 , m; L t s is the length of the s segment of the branching t, m; k is infinite well row reflected by mirror image of branching in Z direction.

The wellbore pressure distribution equation of herringbone-like laterals wells wells under bottom water reservoir is obtained according to Equation (8):

p w t s = p e + μ K ( Φ e t s − Φ t s ) (9)

For bottom water reservoirs, there is:

Φ e t s = 0 (10)

The FORMULA (9) is changed to:

p w t s = p e − μ k Φ t s (11)

where μ is oil viscosity, mPa・s; K is reservoir permeability, mD; p w t s is the well bottom hole flow pressure of the s segment of the branching t, MPa; p e is initial reservoir pressure, MPa.

3. Pressure Drop Equation of Wellbore Flow3.1. Computational Model of Flow Pressure Drop in Main and Branch Wellbore

The main wellbore and branching wellbore are divided into many units. Considering that the length of each unit is enough small, the pressure drop of the second section of the branching t of herringbone-like laterals wells is calculated in reference [14] :

Δ p w t s = 2 f h w ρ π 2 D 5 ( 2 Q t s + q t s ) 2 Δ x + 16 ρ q t s π 2 D 4 ( 2 Q t s + q t s ) (12)

where q t s is radial inflow of micro-element section in the s segment of the branching t, m3/s; Q t s is the upstream flow of the mainstream for the s segment of branching t, m3/s; f h w is friction resistance coefficient of the well tube wall with radial inflow, f; ρ is fluid density, kg/m3; D is the wellbore diameter, m.

3.2. Computational Model of Convergent Flow Pressure Drop of Main and Branching Wellbore

While the fluid of the main wellbore and branching wellbore is confluencing, a mixed loss will occur, and generated local pressure drop. Based on the principle of fluid mechanics, a calculation model of local pressure drop at the confluence point of main wellbore and branching wellbore is established.

Assuming that the fluid in the wellbore flows steadily, adiabatically and isothermally, without considering the friction between the fluid and the pipe wall, and ignoring the influence of gravity, the confluence flow diagram of the branching wellbore is shown in Figure 2.

The momentum equation of the main wellbore direction of the fluid at the confluence point is as follows:

p 1 π 4 D 2 − p 2 π 4 D 2 + F x = ρ Q 2 V 2 − ρ Q 1 V 1 (13)

Continuity equation:

V 1 π 4 D 2 + q t = V 2 π 4 D 2 (14)

Energy equation:

p 1 ρ g + V 1 2 2 g = p 2 ρ g + V 2 2 2 g + h 12 (15)

where p 1 , p 2 are pressure at the inflow and outflow ends along the direction of the main wellbore at the confluence point, MPa; D is main wellbore and branching wellbore diameter, m; q t is from the branching wellbore t to main wellbore, m3/s.

The force F x of the wall acting on the gas at the junction point can be derived from the momentum equation:

F x = ρ q t V 3 cos φ (16)

Combination the Formula (13) and (16):

p 1 − p 2 = 4 ρ π D 2 ( Q 2 V 2 − Q 1 V 1 − q t V 3 cos φ ) (17)

Combination the Formula (14), (15) and (17), pressure drop equation at the confluence point of main and branching wellbore is:

Δ p w t = p 1 − p 2 = 16 ρ q t 2 π 2 D 4 + 8 V 1 ρ q t π D 2 − 4 V 3 ρ q t π D 2 cos φ (18)

Substitute V 1 = 4 Q t π D 2 , V 3 = 4 q t π D 2 into the Formula (18):

Δ p w t = p 1 − p 2 = 16 ρ q t 2 π 2 D 4 ( 1 − cos φ ) + 32 ρ Q t q t π 2 D 4 (19)

where Q t is flow from upstream end of the point of main wellbore t, m3/s; q t is flow from branching wellbore t to the main wellbore, m3/s; φ is the angle of branching, degrees; D is diameter of the main wellbore and branching wellbore, m.

4. Coupling Model of Seepage and Wellbore Flow in Herringbone-Like Laterals Wells

In addition to flowing along the length of horizontal wellbore, reservoir fluid also flows into wellbore along the horizontal wellbore direction. There is a coupling relationship between seepage flow in reservoir and in the wellbore.

4.1. Establishment of Coupling Model

The pressure distribution in wellbore can be calculated by the pressure drop calculation model as:

p w t s = p w t ( s − 1 ) + 0.5 ( Δ p w t ( s − 1 ) + Δ p w t s )       ( 0 ≤ t ≤ M − 1 , 1 ≤ s ≤ N ) (20)

Δ p w t 0 = 0 , p w t 0 = p w f t

where p w f t is the flow pressure at the heel of the branching wellbore t, MPa; p w t s is the flow pressure of s segment of the branching t, MPa; p w t ( s − 1 ) is the flow pressure of s−1 segment of the branching t, MPa; Δ p w t ( s − 1 ) is the pressure drop of s−1 segment of the t branch, MPa; Δ p w t s is the pressure drop of s segment of the branching t, MPa.

According to the principle of material balance, the inflow of each branching wellbore equals the sum of the inflow of each small section at the upstream:

Q t s = ∑ t = 0 M − 1 ∑ s = 1 N q t s (21)

As can be seen from the above, the reservoir seepage model has m × n equation, the wellbore pressure drop model has m × n equation, there are 2 m × n equations. The variables to be solved are q t s and p w t s ( 0 ≤ t ≤ M − 1 , 1 ≤ s ≤ N ) , which are also 2 m × n , so the equations are closed.

4.2. Solution of Coupled Model

The coupling model is solved by iteration method. The specific steps are as follows: 1) Assuming that the initial value of p w t s is p w t s 0 , in actual calculation it can be assumed p w t s 0 = p w f t ; 2) Substitute p w t s into the Formula (11), used gauss elimination method to find q t s ; 3) Substitute q t s into the Formula (21), find Q t s ; 4) Substitute q t s and Q t s into the Formula (12) and (19), find Δ p w t s ; 5) Substitute Δ p w t s into the Formula (20), to update p w t s . This value is taken as the initial value of the next iteration; 6) repeat (2) - (5), comparing p w t s n + 1 , q t s n + 1 after n iterations with p w t s n and q t s n after n iterations, when both of them satisfied certain accuracy, the iteration stops, otherwise repeat (2) - (5) steps until the accuracy is satisfied; 7) Last, the Formula (21) can be used to calculate the total production of herringbone-like laterals wells.

5. Case Study and Optimization of Morphological Parameters5.1. Example Analysis

The herringbone-like laterals wells in a heavy oil reservoir in Bohai Oilfield as an example, the well pattern is shown in Figure 3. There are four branching along the main wellbore. The distance between each branching is 50 m, the angle between the main wellbore and the branching wellbore is 45 degrees, the length of the main wellbore is 400 m, the length of each branching wellbore is 200 m, and the radius of the main wellbore and the branching wellbore is 0.119 m, reservoir thickness h = 6.9   m , distance from the well to bottom of reservoir Z w = 3.45   m , reservoir permeability K = 3159 × 10 − 3   μ m 2 , initial reservoir pressure p e = 11.47   MPa , bottom hole flow pressure p w = 10.47   MPa , production pressure drop Δ p = 1   MPa , branching wellbore radius r w = 0.119   m , wall roughness e = 0.001   m , oil viscosity μ = 50   mPa ⋅ s , volume coefficient of oil B = 1.07 , oil density ρ = 0.969   g / cm 3 , the threshold pressure gradient of the reservoir is 0.02 MPa−1/m.

The productivity coupling model deduced by the author is used to predict the productivity of herringbone-like laterals wells. Compared with the actual production data, as shown in Table 1, it can be seen from Table 1 that the relative error between the calculated results and the actual production is less than that calculated by Liu Xiangping formula, which is 7.2%. The main reason is that the effect of threshold pressure gradient on productivity is considered in this paper. This formula has high practicability for predict the productivity of herringbone-like laterals wells.

Comparison of actual data and productivity of herringbone-like laterals wells
WellA50H
The methodProductivity (m3/d)Relative error
Formula of Liu Xiangping2677.2
Formula in this paper2542.4
Actual data250-
5.2. Study on Optimization of Morphological Parameters

The reservoir-wellbore steady seepage coupling model is used to optimize the shape of herringbone-like laterals wells, for give full play to the advantages of herringbone-like laterals wells.

1) Branching length optimization

The optimization of branching length mainly studies during the total wellbore length is equal, the productivity difference between equal and unequal branching length. In the study, the structure of two branching shown in Figure 4. The total length of the main wellbore and branching wellbore is 800 m, and the angle of each branching is 45 degrees.

Figure 5 shows that when the total length of the wellbore is equal, during the length of the main wellbore and branching wellbore is increased; the productivity of the well is increased. When the branching length is 0 m, the well will be transformed into a single horizontal well, that is, the productivity of a single horizontal well is greater than an equal length herringbone-like laterals wells.

2) Branching angle optimization

Using the three branching structures as an example, the effect of branching angle on productivity is studied. As can be seen from Figure 6, the total angle of the three branching structures is 135 degrees. The productivity variation law is studied by changing the angle of each branching. The main wellbore and branching wellbore length are 400 m and 200 m respectively.

From Figure 7, it can be seen that the productivity of the well is the smallest when the branching angle is equal and the productivity is increased with the increase of the angle difference. Generally, the change of branching angle has little effect on the total productivity of the well, less than the influence of branching length.

3) Branching symmetry optimization

In order to study the effect of branching symmetry on productivity of the well, four branching structures are designed as shown in Figure 8. This paper mainly analyses whether there are common convergence points between the branching and the influence of the branching on the productivity of the well. The main wellbore length is 400 m, the branching wellbore length is 200 m, and the angle of each branching is 45 degrees.

As can be seen from Figure 9, that the ipsilateral branching structure will contribution more productivity than opposite side branching structure. During the number of branching on one side of the main wellbore is increasing, the interference on the one side of the main wellbore is increasing, and the driving area of each branching is increasing. The result of comprehensive action increases the total productivity of the well.

4) Branching location and spacing optimization

In order to study the influence of branching location and spacing on the productivity of the well, four branching structures are designed as shown in Figure 10. The main wellbore length is 400 m, the branching wellbore length is 200 m, the angle of branching is 45 degrees, and the distance between the main wellbore and the branching heel is shown in Figure 10.

As can be seen from Figure 11, that during the branching is closer to the heel of the main wellbore, the productivity of the well is larger. When branching spacing is increased the interference of branching is decreased, but the interference range of main wellbore is increased, the result of comprehensive action decreases the total productivity of the well.

6. Conclusion

Based on the threshold pressure gradient, the productivity coupling model of herringbone-like laterals wells is established in heavy oil reservoir-wellbore steady seepage. The productivity coupling model is suitable for predicting the productivity of herringbone-like laterals wells in heavy oil reservoir. Using the productivity coupling model in this paper, the well shape parameters of the well are optimized, and the principle of optimizing the well shape of herringbone-like laterals wells is proposed.

Conflicts of Interest

The authors declare no conflicts of interest regarding the publication of this paper.

Cite this paper

Sun, E.H., Tan, J., Zhang, D., Wang, W. and Mu, S.R. (2019) Study on Productivity Model of Herringbone-Like Laterals Wells and Optimization of Morphological Parameters Considering Threshold Pressure Gradient in Heavy Oil Reservoirs. World Journal of Engineering and Technology, 7, 302-313. https://doi.org/10.4236/wjet.2019.72021

ReferencesSalas, J.R., Clifford, P.J. and Jenkins, D.P. (1996) Multilateral Well Performance Prediction. SPE 35711, SPE Western Regional Meeting, Anchorage, 22-24 May 1996, 1-9. https://doi.org/10.2118/35711-MSHan, G.Q., Wu, X.D. and Chen, H. (2004) Analysis of Factors Affecting the Productivity of Dual-Branch Wells in Multi-Layered Heterogeneous Reservoirs. Journal of Petroleum University (Natural Science Edition), 28, 81-85.He, H.f., Zhang, C., Fu, X., et al. (2004) Calculating the Productivity of Fishbone Branch Wells by Nodal Method. China Offshore Oil and Gas, 16, 263-265.Liu, X.P., Yu, G.D. and Li, Z.P. (2006) Study on Productivity of Complex Branch Horizontal Wells. Petroleum Exploration and Development, 33, 729-733.Yang, X.S., Liu, C.X., et al. (2008) Study on Productivity Law of Fishbone Multi-Branch Horizontal Gas Wells. Journal of Petroleum, 29, 727-733.Zhao, G.Y. (2014) Flowing Feature and Optimization of Multi-Branch Horizontal Well in Offshore Low Permeability Reservoirs. China University of Petroleum (East China), Dongying.Huang, Y., Cheng, S.L., He, Y.W., et al. (2016) Transient Pressure Analysis of Fishbone Multi-Lateral Horizontal Well with Non-Uniform Flux Density. Journal of Shenzhen University Science and Engineering, 33, 202-208. https://doi.org/10.3724/SP.J.1249.2016.02202Huang C. ,et al. (2017)Research on Application Technology of Multi-Branch Horizontal Wells Inner Mongolia Petrochemical Industry 34, 73-75.Li, C.L. and Zhang, S.C. (2010) Steady-State Productivity Formula for Fishbone Branch Wells. Journal of Daqing Petroleum Institute, 34, 56-59.Fan, Y.P., Qing, K. and Yang, C.C. (2006) Productivity Prediction of Yugu Well and Shape Optimization of Branch Wells. Journal of Petroleum, 27, 101-104.Ye, S.J., Jiang, H.Q., Zhu, G.J., et al. (2010) Productivity Prediction and Influencing Factors Analysis of Fishbone Wells. Fault Block Oil and Gas Fields, 17, 341-344.Zhang, S.M., Zhou, Y.J., Song, Y., et al. (2013) Shape Design Optimization of Fishbone Branch Horizontal Wells. Petroleum Exploration and Development, 38, 606-612.Zhang, D.Y., Peng, J., Gu, Y.L., et al. (2012) Start-Up Pressure Gradient Experiment for Heavy Oil Reservoirs. Xinjiang Petroleum Geology, 33, 201-204.Liu, X.P., Zhang, Z.S., Cui, G.X., et al. (2000) Inflow Performance Relationshio of a Herringbone Multilateral Well. Journal of Petroleum, 21, 57-60.