Study and Application of Horizontal Well Productivity Formula with Continuous Time-Varying Considered

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.


Introduction
Horizontal well productivity formula has been studied and widely applied from 1940s [1].It has been numerously researched by domestic and foreign scholars https://doi.org/10.4236/jpee.2018.69005[2] [3].The foreign researchers Giger [4] calculated the horizontal well productivity with hodograph curve method at the assumption of the well location being at the top of the infinite layer.Joshi [5] had the horizontal well 3-D flow field simply and held that the horizontal well real 3-D flow could be approximately made up of the inner radial vertical plane flow and external horizontal elliptical flow, after which the final productivity formula could be solving by principle of equivalent flow resistance.Babu [6] derived the horizontal well productivity formula through physical model.The domestic researcher Chen Yuanqian [7] got the other horizontal well productivity formula based on Joshi Formula.In the view of mathematical model, Chen Xiaofan [8] simplified the formula initial boundary condition and derived the analytical formula.According to Economides horizontal well productivity formula, Gao Jian [9] analyzed the main influencing factors.Gao Hongmei [10] researched the development of different well productivity formula.Zhang Feng [11] compared and analyzed differences of the horizontal well productivity formula at steady state time between Chen Formula and Joshi Formula.Ding Yiping [12] determined the productivity formula of fractured horizontal well multiple transverse cracks by the conception of equivalent radius and analyzed the geometric and physical parameters of reservoir model and fracture.Xu Guomin [13] studied remaining oil potential tapping in the maturing area combined with fine preliminary geological study and horizontal well development technology advantages.
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.

Methodology
Providing that at the mid of the horizontal homogeneous circular reservoir with constant thickness, there was a horizontal well producing, which showed as Fig- ure 1.The formation and well parameters were known which including that the formation radius is r 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.Journal of Power and Energy Engineering 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: During the pressure transmission, spread radius is investigation radius r i : 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 Figure 2.
From Figure 2, it can be seen that the flow of the fluid in the horizontal plane could be approximately recognized as elliptical flow which was showed as Figure 3.
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: The horizontal plane and vertical plane flow rate formulas of horizontal well were derived.

Horizontal Well Planar Flow Analysis
According to Zhukovsky conversion, the plan elliptical flow of horizontal well can be converted to circular flow, the formula is: , 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 Figure 4.
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 ( ) , the flow rate in the horizontal plane can be calculated with: ( ) 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:

Horizontal Well Vertical Flow Rate Analysis
According to the conformal transformation, the vertical flow of horizontal well was converted into circular flow.
. 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.

Application
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 Figure 5 and Figure 6.The error analysis was showed as Table 1 and Table 2.
In the view of Table 1 and Table 2, the production of some adjusting wells was designed with the horizontal well productivity formula.The actual production data showed that the errors were little enough to confirm the formula reliable, which provided the theoretical basis and be applied in determining the production in offshore oil field.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.
How to cite 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.

Figure 1 .
Figure 1.The figure of horizontal well production.

Figure 2 .
Figure 2. The figure of flow in the horizontal well production.

Figure 3 .
Figure 3.The figure of ellipse flow in horizontal well production.

Figure 4 .
Figure 4.The Figure of ellipse flow conformal transformation in horizontal well production.

Figure 5 .
Figure 5.The production calculation curve of the two horizontal wells in B oilfield.

Figure 6 .
Figure 6.The production calculation curve of the two horizontal wells in C oilfield.
et al.Journal of Power and Energy Engineering

Table 1 .
The comparison table of actual production data and formula calculation data in B oilfield.

Table 2 .
The comparison table of actual production data and formula calculation data in C oilfield.