Research Article  Open Access
Haidong Wang, Yikun Liu, "Horizontal Well Completion with Multiple Artificial Bottom Holes Improves Production Performance in Bottom Water Reservoir", Mathematical Problems in Engineering, vol. 2020, Article ID 7247480, 11 pages, 2020. https://doi.org/10.1155/2020/7247480
Horizontal Well Completion with Multiple Artificial Bottom Holes Improves Production Performance in Bottom Water Reservoir
Abstract
The horizontal well completion with stinger is usually used to control the bottom water cone. Although the pressure profile and the inflow profile along the horizontal wellbore can be divided into two parts by the stinger, these profiles have not really flattened. In order to flatten the pressure distribution and inflow distribution further, it proposes a new technology. This new horizontal well has multiple artificial bottom holes (MABH) along the wellbore and it has application potential. In order to verify the effectiveness of MABH technology, a model of horizontal well completion with MABH was established, and the production performance of different water cone control technologies was analyzed: conventional horizontal well, stinger completion horizontal well, and MABH completion horizontal well. The results show that the MABH technology has more advantages than the stinger technology. The uniformity of pressure distribution of the 6MABH horizontal well is 55% higher than that of the horizontal well with string technology, and the uniformity of inflow distribution is increased by 65.25%. At the same time, although the operation of MABH technology is very simple, it should follow a rule of MABH installation: the position of the first MABH should be set at 242.5 m from the heel hole of the horizontal wellbore, and the other interval is 92.4 m.
1. Introduction
The water cone is formed due to the pressure loss from the toe to the heel of the horizontal well [1, 2]. It causes less inflow contribution from the upstream portion near the toe than from the downstream portion of the heel [3–5]. Once the bottom water breaks through into the wellbore, the water cut will rise rapidly, the production efficiency will drop significantly, and the life of the well will be greatly shortened. In order to solve this problem, some studies have been proposed as follows: the method of variable density perforation [2, 6] can delay the breakthrough time of water. Landman and Goldthorpe [4] established a predictive model that can calculate the inflow profile along the wellbore. In addition, the best variable density perforation technology [7] is also reported. Although this method has the effect of controlling the bottom water cone, it cannot solve the problem of pressure loss from the heel to the toe of the horizontal well.
In the early 1990s, Norwegian Hydropower engineers first proposed the concept of an inflow control device (ICD) and developed a tortuous flow control device to balance the influx of various parts along the horizontal well. Subsequently, the oil service company developed pipeline channel ICD, orifice plate ICD, automatic mixing ICD, and nozzle ICD [8–12]. Now ICD technology has been well developed. However, ICD still has some problems, such as the problem of being eroded and blocked by particles in the fluid.
On the other hand, technology of horizontal well completion with stinger was proposed [13]. It is pointed out that the stinger technology can delay the bottom water cone of long horizontal wells. The principle is to redistribute the pressure distribution of the hole by changing the actual position of the bottom hole. The optimal length of the stinger under various oilwater viscosity ratios and oil column thickness conditions were studied [5]. The stinger technology can make the wellbore pressure distribution flat, but the uniform effect is not ideal. Therefore, the unified problem of stinger technology should be studied in detail [14].
In this paper, a new method is proposed to further flatten the pressure profile and inflow profile along the horizontal wellbore. This method is horizontal well completion with multiple artificial bottom holes (MABH). At the same time, the pressure profile and the inflow profile along the new horizontal wellbore have been theoretically studied.
2. Horizontal Well Completion with MABH
In this article, we set up some independent annulus tubes through MABH to connect the base pipe (see Figure 1) and then use these MABH and annulus tubes to convert the unidirectional flow in the base pipe into multidirectional flow. Likewise, it is foreseeable that the pressure distribution along the base tube will exhibit periodic wavelets (red line) and that the red line will be flatter than the conventional pressure distribution (blue dotted line). Similarly, as can be seen from Figure 1, the pressure drop difference between the reservoir and the new wellbore (N_{1} ≈ N_{2} ≈ N_{3}) is much better than the pressure drop difference between the reservoir and the conventional wellbore (C_{1} < C_{2} < C_{3}). Then compare the new inflow profile (N_{11} ≈ N_{22} ≈ N_{33}) with the traditional inflow profile (C_{11} > C_{22} > C_{33}). Therefore, the heeltoe effect will be weakened.
3. Calculation Method
In order to explore the use of MABH technology to reduce the efficiency of the heeltoe effect, relevant theoretical calculations were carried out. The detailed calculation method is as follows.
Wellbore flow model: divide the horizontal wellbore into N segments, and remove one segment for force analysis (Figure 2).
According to the law of conservation of mass, the equation can be expressed as follows:where is the fluid density in the pipe, kg/m^{3}; is the crosssectional area of the tube, m^{2}; is the flow velocity, m/s; is the opening coefficient of the base pipe; t is the flow time, s; is the length of the infinitesimal segment, m; and is the inner diameter of the pipe, m.
When becomes 0, according to the definition of derivative, equation (1) can be acquired as follows:
Then, according to the principle of momentum, the equation can be expressed as follows:where is the gravitational acceleration, m/s^{2}; is the angle between wellbore axis and horizontal direction, degrees; and is the average frictional stress, Pa.
Assuming that the flow rate is stable, and when ∆x becomes 0, combining equations (2) and (3), the pressure drop in the base pipe can be expressed as follows:where is the flux of th segment in the wellbore, m^{3}/s; is the influx rate of the th wellbore segment, m^{3}/s; is the correction coefficient caused by wear and corrosion in the wellbore, dimensionless; and is the correction coefficient caused by variable mass pipe flow in the wellbore, dimensionless.
Similarly, pressure drop in the annulus tube can be represented as follows:
The average frictional stress can be represented as follows:
According to the research results of Ouyang [15], is shown in Table 1.

The following equation gives the friction factor:
Here,where is the correction coefficient caused by variable mass pipe flow in the wellbore, dimensionless; is the Reynolds number, dimensionless; and is the friction factor, dimensionless.
Reservoir inflow model: firstly, we make the following assumptions: ① The bottom water reservoir is homogeneous and infinite. ② There is a closed boundary at the reservoir top and a constant pressure boundary at the bottom. The distance between the upper and lower boundaries is h, and the distance between the horizontal wellbore and the lower boundary is . The reservoir temperature is constant. ③ The length of the horizontal wellbore is much larger than its diameter. ④ The model is assumed to be singlephase flow. If the pressure drops below the bubble point, the model is not available. ⑤ The flow in the reservoir follows Darcy’s law.
The above assumptions will affect the calculation results of the total output. However, the focus of this article is to study the performance of the inflow profile and the pressure profile of horizontal well. Therefore, these assumptions do not affect the results of this article.
First, the horizontal wellbore is divided into N sections, and the penetration velocity of the th section surface can be expressed as follows:where is the production rate of the th wellbore segment, m^{3}/s; is the perforated opening area on the pipe, m^{2}; is the seepage velocity, m/s; is the seepage area, m^{2}; is the opening coefficient of the base tube; is the length of th segment, m; and is the inner radius of the pipe, m.
In addition, the percolation rate can be represented as follows [16, 17]:where is the permeability in the reservoir, m^{2}; is the viscosity, ; is nonlinear permeability coefficient near wellbore, m^{2}; is pressure gradient, MPa/m; m is the experimental constant, and it is between 1.2∼1.5; represents the end point of linear pressure gradient, between 4∼10, MPa/m.
When , the potential of the th segment in the infinite space can be obtained as follows:
Through the mirror image method (see Figure 3) and the principle of potential superposition, the potential at a random point β (x, y, z) can be obtained as follows:
Here,
Furthermore, the potential of (x, y, z) at a random position caused by entire horizontal wellbore can be obtained as follows:
Here,
Therefore, the reservoir flow in th segment can be described as follows:where is the pressure at the constant pressure boundary, Pa; is the bottom hole pressure of the th segment, Pa; is the potential, Pa; is the xcoordinate value of the th segment, m; is the zcoordinate value of the th segment, m.
The vector representation is given by
Here,
When , combine equations (9) and (10), flowing near the wellbore of the th segment can be described as follows:
Therefore, the actual reservoir flow model is as follows:
In addition, there is oilwater flow interaction in production. Therefore, combining equations (4) and (20), the oilwater wellbore flow model and the reservoir flow model can be obtained as follows:where is the water density in the pipe, kg/m^{3}; is the oil density in the pipe, kg/m^{3}; is the water flux in the th section of the wellbore, m^{3}/s; is the water flux in the th section of the wellbore, m^{3}/s; is the water inflow rate in the th wellbore segment, m^{3}/s; is the th water production rate of the wellbore segment, m^{3}/s; is the nonlinear permeability coefficient of the water phase, m^{2}; is the nonlinear permeability coefficient of the water phase, m^{2}; is the th water production rate of the wellbore segment, m^{3}/s; is the oil production rate of the th wellbore segment, m^{3}/s.
4. Coupling Model Solution
The focus of this article is to study the performance of the inflow profile and pressure profile along the horizontal well. Therefore, singlephase flow or twophase flow does not affect the results of this article and then the coupling solution takes singlephase flow as an example. Step 1: if a set of bottom hole pressures along each ith section of the horizontal wellbore are given , the corresponding radial inflow can be calculated by equation (20) through the Gaussian elimination method. Step 2: divide the horizontal wellbore into segments (Figure 4). The part of the wellbore flow will be controlled by the heel hole. Assuming that the radial inflow along the horizontal wellbore is , the corresponding pressure difference of all pipe sections in the base pipe can be calculated by formula (4), . is known, and the value of can be represented as follows: The centre pressure of the th segment of wellbore is represented as follows: The end pressure of the th segment of wellbore is represented as where . A new group of bottom hole pressures can be obtained in section , . Step 3: the section is controlled by the first artificial bottom hole. In section , it is divided into two parts: is the left part, and is the right part.(1)For segment , firstly assume a group of radial influx, . Then the corresponding differential pressure of random segments in the wellbore M_{21} can be calculated by equation (4). From segment (M_{1} + M_{21}) to segment (M_{1}), the centre pressure of the th segment is represented as follows: From segment (M_{1} + M_{21}) to segment M_{1}, the end pressure of the th segment is represented as follows: where The pressure of can be calculated as follows: A new set of bottom hole pressures can be obtained in section , . Note. If the value of calculated in Step 3 is different from that of calculated in Step 2, there will be two kinds of discussions:①. Part of in Step 2 is not accurate and it should be recalculated. Calculation should be started from segment (M_{1} + M_{21}) to () until . Then, an accurate set of bottom hole pressures can be obtained in section , .②. Part of in Step 3 is not accurate and it should be recalculated. Calculation should be started from segment M_{1} to segment (M_{1} + i) until . Then, an accurate set of bottom hole pressures can be obtained in section , .(2)For segment , first assume a group of radial influx, . The corresponding differential pressure of random segments in the wellbore M_{22} can be calculated by equation (4). The centre pressure of the th segment of wellbore is represented as follows: The end pressure of the th segment of wellbore is represented as follows: where The pressure of can be calculated as A new set of bottom hole pressures can be obtained in section , . The main flow of segment is known and can be expressed as follows: Step 4: repeat Step 2 and Step 3 to calculate corresponding pressure distribution controlled by the next annulus bottom hole. Step 5: compare the pressure with , if the inaccuracy meets the engineering inaccuracy requirement, and stop calculating or replacing Step 1 until the desired value is reached. The pressure profile and the inflow profile will be obtained.
5. Results and Discussion
This article focuses on a new technology to improve the impact of horizontal wellbore pressure profile and production profile, so a singlephase flow model is used. Here, we have designed three schemes (Table 2) to predict production performance. The basic data of well XU2 is given in Table 3. The calculation method of scheme 2 is based on Permadi and Wibowo [5].


6. Optimal Water Control Method
Figure 5 shows three pressure distributions along the horizontal well. Here, (4), (5), and (6) represent the difference between the maximum pressure and the minimum pressure along the horizontal well (blue line). As shown in Scheme 1, the pressure distribution along the horizontal well increases monotonously, and the value of (4) is 0.0964 MPa. In Scheme 2, it can be seen that the pressure distribution along the horizontal well has been divided into two parts by the stinger and is funnelshaped. The value of (5) is 0.0889 MPa. In Scheme 3, the pressure distribution of MABH horizontal wells shows wavelets, indicating that the pressure distribution of horizontal wells is divided into several wavelets using MAHB. Here, the value of (6) is only 0.04 MPa. In summary, the pressure curve of Scheme 2 is better than the pressure curve of Scheme 1, and the uniformity of the former is 7.8% higher than that of the latter. The pressure curve in Scheme 3 is better than the pressure curve in Scheme 2, and the uniformity of the former is improved by 55% compared with the latter. Obviously, the new pressure profile of Scheme 3 is production desired.
Figure 5 also shows the distribution of inflow profiles along horizontal wells. Here, (1), (2), and (3) represent the difference between the maximum inflow and the minimum inflow along the horizontal wellbore, respectively. As can be seen from Scheme 1, the value of (1) is 0.0561 (m^{3}/d)/m. It can be seen from Scheme 2 that the difference can be reduced by using a stinger, but the inflow fluctuation range between the funnel point and the end of the wellbore is still very large, and the value of (2) is 0.0468 (m^{3}/d)/meter. It is worth noting that in Scheme 3, by using MABH, the value of (3) is only 0.01626 (m^{3}/d)/m. Compared with (3) and (2), the uniformity of the former is improved by 65.25% than the latter. Therefore, it can be concluded that MABH can further flatten the inflow profile.
7. Optimal Installation Interval of MABH
Figure 6 shows the pressure profile and inflow profile. The inflow profile gradually flattens and periodically flattens as the number of MABH increases. In addition, when the same amount of MABH is installed at the toes, the uniformity of the pressure profile and the inflow profile will be more uniform than the installation at the heel. When the length of the base pipe controlled by the wellbore hole is less than 242.5 m and the length of the MABH controlled base pipe is less than 92.4 m, the difference between the maximum inflow and the minimum inflow along the horizontal wellbore does not exceed 3% of the average value. Therefore, we set the first MABH to 242.5 m along the horizontal well, and other installation intervals to 92.4 m. But this rule only applies to wellbore with the above parameters. For other wellbores, we can recalculate the optimal length of the base pipe controlled by MABH according to the above method.
8. Operation and Challenges
The annulus tube can be oval and can be extruded through a seamless steel tube. It is locked between the base pipe and the screen pipe by the support ring and crossover. Here, you can set a MABH on the cross ring and connect the MABH to the annulus and the base pipe (see Figures 1 and 7(e)), respectively. Then, Figures 7(a)–7(e) show the ground installation process of the production unit with MABH. Install the screen spacers, wire wrap screen, crossover with MABH and support rings, annulus tubes, and screen shroud at the corresponding positions of the base pipe until the production unit shown in Step 6 is formed. Then, it is possible to form horizontal wells with MABH by connecting each production unit in accordance with the steps given in Figures 7(g), 7(h), and 7(i). In particular, after tightening the connection of the base pipe, annulus tubes are automatically further inserted upstream of annulus tubes on both sides, which has a good sealing and selflocking function. The completion structure of MABH horizontal well is complicated, but the operation is simple and convenient. Similarly, when the length of the annulus tube is designed to be 10 m and the cross section is designed to be 314 mm^{2}, the cost of an annulus tube is about $ 610. This is the challenge of developing lowyield fields. But it is suitable for highyield bottom water oil fields or gas fields.
(a)
(b)
(c)
(d)
(e)
(f)
(g)
(h)
(i)
9. Further Discussion
The reason for the multidirectional flow in the base pipe can be obtained from the following equation:where is the pressure drop along the wellbore, MPa; is the gravity pressure drop, MPa; is the friction pressure drop, MPa; and is the acceleration pressure drop, MPa.
The results show that, compared with the constant mass flow in the annulus tube, the various mass flows in the base pipe overcome the greater flow resistance. The reason is because there are perforated openings in the base pipe, but not in the annulus tube, so the additional in the base pipe will be formed by radial inflow acceleration. Therefore, the lower pressure in the base pipe at the toe can be controlled by using MABH, and then the fluid at the toe can easily bypass the heel (as shown in Figures 1 and 8). The influx of toes in the wellbore has increased, which is beneficial for reducing the heeltoe effect of horizontal wells. In short, MABH technology should be considered when controlling the water cone problem of horizontal wells for better reservoir management.
(a)
(b)
(c)
10. Conclusions
Here are some suggestions and conclusions:(1)A new MABH completion method is established to solve the water coning problem. The main mechanism of MABH technology is to reduce the original pressure in the conventional horizontal wellbore by using MABH. The original streamline from the toe to the heel in a conventional horizontal wellbore is changed by multidirectional flow. Then, the pressure distribution along the base pipe can be flattened, and the inflow profile along the horizontal wellbore can be made uniform.(2)A coupling model of horizontal well with MABH and reservoir is proposed, which can be used to calculate the pressure loss and inflow difference along the horizontal well. Then, the dynamics of water control technologies are compared based on the output. The influence of the type of wellbore on the development effect is also discussed. Both methods can effectively control the pressure profile and the inflow profile, but the system with 6MABH has significantly improved the pressure profile and the inflow profile than the system with stinger. The uniformity of the pressure distribution of the horizontal well using 6MABH is 55% higher than that of the horizontal well using the stinger, and the uniformity of the inflow profile is increased by 65.25%.(3)Although the MABH structure is complex, the operation is simple and convenient. The cost of a MABH and an annulus tube with a length of 10 m and a cross section of 314 mm^{2} is only $ 610, which is a challenge for the development of lowyield oil fields. At the same time, the first MABH should be set at 242.5 m away from the heel hole of the horizontal wellbore, and the other interval is 92.4 m.
Data Availability
The .xlsx data used to support the findings of this study are included within the article and can be downloaded from https://pan.baidu.com/s/1Bj6Ld1EHi4eH62N4200Pcg, and the supporting data extraction code is p7s1.
Conflicts of Interest
The authors declare no conflicts of interest.
Acknowledgments
This study was supported by Major Project of China National Offshore Oil Corporation 13th fiveyear plan (Grant no. CCL2018ZJFN0462).
References
 B. P. Marett and M. J. Landman, “Optimal perforation design for horizontal wells in reservoirs with boundaries,” in Proceedings of the SPE Asia Pacific Oil and Gas Conference, Singapore, February 1993. View at: Publisher Site  Google Scholar
 H. J. Yuan, CemSarica, and J. P. Brill, “Effect of perforation density on singlephase liquid flow behavior in horizontal well,” SPE Production & Facilities, vol. 14, 1999. View at: Publisher Site  Google Scholar
 I. Stanley Okafor, C. Ebenezer Ubani, and N. Anyanwu, “Development strategies for oil reservoirs with coning problems: a simulation case study,” in Proceedings of the SPE Nigeria Annual International Conference and Exhibition, Lagos, Nigeria, August 2018. View at: Publisher Site  Google Scholar
 M. J. Landman and W. H. Goldthorpe, “Optimization of perforation distribution for horizontal wells,” in Proceedings of the SPE AsiaPacific Conference, Perth, Australia, November 1991. View at: Publisher Site  Google Scholar
 P. Permadi and W. Wibowo, “Horizontal well completion with stinger for reducing water coning problems,” in Proceedings of the SPE Production Operations Symposium, Oklahoma City, OK, USA, March 1997. View at: Publisher Site  Google Scholar
 T. Yildiz, “Inflow performance relationship for perforated horizontal wells,” in Proceedings of the SPE Production and Operations Symposium, Oklahoma City, OK, USA, March 2001. View at: Publisher Site  Google Scholar
 T. Yildiz, “Productivity of selectively perforated horizontal wells,” in Proceedings of the SPE Production & Operations, Oklahoma City, OK, USA, February 2006. View at: Publisher Site  Google Scholar
 A. B. Elverhϕy and H. Aakre, “Autonomous inflow control for reduced water cut and/or gas oil ratio,” in Proceedings of the Offshore Technology Conference, Houston, TX, USA, May 2018. View at: Publisher Site  Google Scholar
 D. Zubarev, E. Nnebocha, and A. Vasper, “Feasibility study and selection process for icd applications in an offshore niger delta field,” in Proceedings of the SPE Nigeria Annual International Conference and Exhibition, Lagos, Nigeria, August 2018. View at: Publisher Site  Google Scholar
 D. C. Haeberle, A. Nwankwo, J. R. Mcdermott et al., “Application of flowcontrol devices for water injection in the erha field,” in Proceedings of the IADC/SPE Drilling Conference, Orlando, Florida, USA, March 2008. View at: Publisher Site  Google Scholar
 C. Jones, Q. Morgan, S. Beare et al., “Design, testing, qualification and application of orifice type inflow control devices,” in Proceedings of the International Petroleum Technology Conference, Doha, Qatar, December 2009. View at: Publisher Site  Google Scholar
 B. Least, S. Greci, and R. Huffer, “Steam flow tests for comparing performance of nozzle, tube, and fluidic diode autonomous ICDs in SAGD wells,” in Proceedings of the SPE Heavy Oil ConferenceCanada, Calgary, Canada, June 2014. View at: Publisher Site  Google Scholar
 K. Brekke and S. C. Lien, “New and simple completion methods for horizontal wells improve production performance in highpermeability, thin oil zones,” SPE Drilling & Completion, vol. 9, 1994. View at: Publisher Site  Google Scholar
 H. Li, B. Jiang, and Y. Wang, “A semianalytical model for horizontal wells with improved stinger completion in heterogeneous bottomwater reservoirs,” Journal of Canadian Petroleum Technology, vol. 53, 2014. View at: Publisher Site  Google Scholar
 L.B. Ouyang, “Single phase and multiphase fluid flow in horizontal wells,” Stanford University, Stanford, CA, USA, 1998, Ph.D. thesis. View at: Google Scholar
 S. Hansbo, “Consolidation equation valid for both Darcian and nonDarcian flow,” Géotechnique, vol. 51, no. 1, pp. 51–54, 2001. View at: Publisher Site  Google Scholar
 S. Hansbo, “Deviation from Darcy’s law observed in onedimensional consolidation,” Géotechnique, vol. 53, no. 6, pp. 601–605, 2003. View at: Publisher Site  Google Scholar
Copyright
Copyright © 2020 Haidong Wang and Yikun Liu. This is an open access article distributed under the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.