^{1}

^{1}

^{1}

^{1}

^{1}

The horizontal well technology has been widely applied to enhanced oil recovery for low permeability and heavy oil reservoir. It is the important basis for designing and optimizing horizontal well to determine the productivity. The productivity determination of horizontal wells in offshore oil fields is mainly based on the actual productivity data of producing directional wells in the similar reservoirs nearby. Considering pressure drop and oil layer thickness to calculate the productivity, this method lacks certain theoretical basis and requires rich working experience for reservoir engineers. The other method is Joshi Formula which needs the known horizontal well control radius to be known. But the control radius is man-made at certain degree. In order to address the shortcomings of existing methods, a new reservoir engineering method was proposed to determine the horizontal well productivity formula, horizontal flow pattern and control radius based on the principle of equivalent flow resistance and conformal transformation. This method has overcome the disadvantage of determining on person. It provided some theoretical basis for getting the horizontal well productivity and is of some guiding meaning for evaluating the productivity of adjustment wells and development wells.

Horizontal well productivity formula has been studied and widely applied from 1940s [

Previous studies have done similar to horizontal wells productivity formula. In the classical formula engineers usually using, there are patterns such as control radius that determines empirically. The productivity formula was improved mainly aimed at this question in this paper. Considering the travel time to the pseudo-steady state time and controlling radius after oil well producing, a new reservoir engineering method was proposed to determine the horizontal well productivity formula based on the principle of equivalent flow resistance, conformal transformation. It gets accuracy of the formula with the producing well data verifying and has the formula applied in the offshore reservoir adjusting wells proration production. It provides theoretical basis for determining the horizontal well productivity in offshore reservoir and has some directive significance.

Providing that at the mid of the horizontal homogeneous circular reservoir with constant thickness, there was a horizontal well producing, which showed as _{e}, the well radius is r_{w}, the horizontal interval is L, the energy at edge of the formation is too enough to keep the pressure P_{e} all the time, the bottom hole flowing pressure is P_{wf}, the thickness of oil layer is h, the permeability of formation is K, the viscosity of fluid is μ. Base on the assumption above the horizontal well flow characteristics and productivity formula can be studied.

The flow of horizontal well in the horizontal plane can be imaged as that countless vertical wells with length of differential units which radius is dL is producing in the horizontal interval L. As the time increasing, the flow of every well unit could be in pseudo-steady state. At that moment the formula of pseudo-steady state starting time was introduced:

The start time of pseudo-steady state:

t s = 0.0872 ϕ μ C t r e 2 K (1-1)

During the pressure transmission, spread radius is investigation radius r_{i}:

r i = 3.386 K t s ϕ μ C t (1-2)

According to the spread radius formula during pseudo-steady state, it could be known that the spread distance of the root and finger on horizontal interval was the longest. The pressure was spread gradually from horizontal interval to exterior of the layer. The flow characteristics of horizontal well were showed as

From

So that the flow of horizontal well can be equivalent to elliptical flow along the horizontal wellbore and plane radial flow along the vertical wellbore. The horizontal well productivity could be calculated with the principle of equivalent flow resistance. Providing the elliptical flow resistance in the horizontal plane is R_{h} and radial flow resistance in the vertical plane is R_{v}, the productivity formula is that:

Q h = P i − P wf R h + R v (1-3)

The horizontal plane and vertical plane flow rate formulas of horizontal well were derived.

According to Zhukovsky conversion, the plan elliptical flow of horizontal well can be converted to circular flow, the formula is:

z = x + i y , ξ = u + i v = r ξ e i θ (1-4)

Z L / 2 = 1 2 ( ξ + 1 ξ ) (1-5)

x L / 2 = 1 2 ( r ξ + 1 r ξ ) cos θ

y L / 2 = 1 2 ( r ξ + 1 r ξ ) sin θ (1-6)

r ξ = u 2 + v 2

cos θ = u u 2 + v 2

sin θ = v u 2 + v 2 (1-7)

According to conformal transformation Z L / 2 = 1 2 ( ξ + 1 ξ ) , the ellipse (long

axis is a, short axis is b) is converted to the circle (the radius is (a + b)/0.5L), the area from (−L/2, 0) to (L/2, 0) is mapped into the unit circle, showed as

Providing that there was a vertical well circular flow with radius 1 in circular area with radius (a + b)/0.5L of ξ plane. Considering the thickness of oil layer h and elliptical radius b = [ a 2 − ( 0.5 L ) 2 ] 0.5 , the flow rate in the horizontal plane can be calculated with:

q hh = 2 π K h ( P i − P wf ) μ B o ln ( a + a 2 − ( L / 2 ) 2 L / 2 ) (1-8)

In the formula (1-8) the long axis (a) of the ellipse was considered to equal the spreading radius (r_{i}) plus the half of horizontal interval (L/2) when the well was producing. So that the expression of value (a) in the formula (1-8) is:

a = r i + L / 2 = 3.386 K t s ϕ μ C t + L / 2

According to the conformal transformation, the vertical flow of horizontal well

was converted into circular flow. ξ = ( 1 − e − π Z h ) ( 1 + e − π Z h ) . The strip area in Z

plane can be converted into the circular area in ξ plane. r_{w} in the Z plane is just as ξ_{w} in the ξ plane.

ξ w = | d ξ d Z | ( 0 , 0 ) r w = 2 π r w h (1-9)

When the directional well which flow radius is ξ_{w} in ξ plane and horizontal interval is L, the vertical plan flow rate formula is:

q hv = 2 π K L ( P i − P wf ) μ B o ln ( h 2 π r w ) (1-10)

The flow resistance formula can be derived by formulas (1-8) and (1-10):

R h = μ B o 2 π K h ln ( a + a 2 − ( L / 2 ) 2 L / 2 ) (1-11)

R v = μ B o 2 π K L ln h 2 π r w (1-12)

The horizontal well productivity formula is:

Q h = P i − P wf R h + R v

Q h = 2 π K h ( P i − P wf ) μ B o [ ln ( a + a 2 − ( L / 2 ) 2 L / 2 ) + h L ln ( h 2 π r w ) ] (1-13)

Including the value a,

Q = 2 π K h ( P i − P wf ) μ B o [ ln 3.386 K t s ϕ μ C t + L / 2 + ( L / 2 + 3.386 K t s ϕ μ C t ) 2 − ( L / 2 ) 2 L / 2 + h L ln h 2 π r w + S ]

The horizontal well productivity formula is show as (1-14). The productivity formula is calculated with horizontal interval, oil layer properties, which solved the disadvantages that the existing formulas could not calculate well controlling radius and controlling area must be given.

According to the horizontal well productivity formula, the production of one horizontal well in Bohai A oilfield. And the formula calculation and actual production data were compared to verify the accuracy of the formula. It was known that according to the horizontal well productivity formula derived in this paper and the properties data of three wells, the theoretical curve of production changed with time. It had good results for the fitting which verified the accuracy of the formula.

The horizontal well technology has been widely applied in offshore reservoir development. It is the important index for producing to calculate the horizontal well productivity. The horizontal well productivity formula has been used to calculate the productivity of four adjusting wells. These wells have been drilled and put on production in Feb. 2014. The results showed that the errors between practical and the theoretical data were very small and the method has some practicality. The results were showed as

In the view of

1) In this paper, the new horizontal well productivity formula was derived according to the principle of equivalent flow resistance, conformal transformation and calculation of pseudo-steady state time. This method has overcome the

Oilfield | Well name | Actual production (m^{3}/d) | Calculate production (m^{3}/d) | Error (%) |
---|---|---|---|---|

B oilfield | B1H | 38.6 | 35.4 | 8.23 |

B oilfield | B2M | 43.9 | 39.5 | 9.96 |

Oilfield | Well name | Actual production (m^{3}/d) | Calculate production (m^{3}/d) | Error (%) |
---|---|---|---|---|

C oilfield | A4H | 145 | 148.4 | 2.36 |

C oilfield | A5H | 165 | 159.7 | 3.21 |

disadvantage of determining on person.

2) The theoretical formula was confirmed by comparison between the calculated and actual production data of three wells in Bohai A Oil field. The results were considered accurate enough to provide basis for adjusting wells production calculation.

3) The production of four adjusting wells in Bohai B oil field was designed with the horizontal well productivity formula. The actual production data showed that the errors were little enough to be applied in determining the production in offshore oil field.

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

Liu, C., Zhang, J.T., Wu, H.J., Zhang, B. and Wang, L. (2018) Study and Application of Horizontal Well Productivity Formula with Continuous Time-Varying Considered. Journal of Power and Energy Engineering, 6, 32-39. https://doi.org/10.4236/jpee.2018.69005