Abstract
Producing an oilfield in a cost-effective way depends on how long water production could be delayed in the reservoir. Many flow mechanisms, correlations, and methods to calculate maximum water-free oil production rate have been published, However, those methods have generally failed to not consider the skin effect which affects the flow into the wellbore. In this paper, the semi-analytical perforation skin model as presented by Karakas and Tariq is incorporated into the Meyer and Garder correlation for critical oil rate from a perforated vertical well interval to obtain the maximum water-free oil production rate and optimal perforation parameters. The resulting coupled computational model is used to determine the sensitivity of the maximum water-free oil production rate to wellbore perforation parameters. Whilst an increase in perforation length and decrease in spacing between perforation increase the critical flow rate, an increase in perforation radius did not translate to higher productivity. The optimal perforation angles are 45° and 60°, however, for the data used in this work the maximum water-free oil rate of 23.2 std/d was obtained at 45° of phasing angle, 1 in of spacing between perforation, 0.36 in of perforation radius and 48 in of perforation length. Thus, the perforation strategy can be optimized prior to drilling and completion operations to improve productivity using the computational model presented in this work.
Similar content being viewed by others
Introduction
An increase in the cost of production operations, environmental problems, reduction in depletion mechanism efficiency and processing of the produced water are major challenges to production and reservoir engineers. Oil or gas reservoirs will produce with increasing water cuts especially from zones with active underlying aquifers that tend to elevate the water–oil contact in the vicinity of the completion interval. This necessitates the determination of the maximum water-free oil production rate that a completed well can deliver. The flow mechanisms, correlations and methods to calculate this critical rate have been developed for wells completed in a continuous pay zone with oil–water contact or gas-oil contact or both (Meyer and Garder 1954; Chierici et al. 1964; Høyland et al. 1989; Menouar and Hakim 1995; Zhang et al. 2011; Tabatabaei et al. 2012). Correlations have also been developed to predict the water breakthrough time and time for a well to produce above its critical rate (Papatzacos et al. 1989; Yang and Wattenbarger 1991). These common analytical formulas have been developed to calculate the critical rate without considering the effect of perforation or using the perforation interval of zero, which is practically impossible. Meanwhile, perforation has been reported to be one of the factors that affect both coning (Li et al. 2015) and productivity (Locke 1981; Karakas and Tariq 1991; Zhang et al. 2012; Luo et al. 2015) of a well. For this reason, Abass and Bass (1988) studied the critical production rate in a water-coning system using numerical and physical models. Abass and Bass (1988) considered the radius of the cone, the maximum height of the water cone and the thickness of the oil zone as the control parameters for the calculation of the critical oil rate. Although their analysis was carried out for a restricted case which can only be approximated to the real conditions, they concluded that the optimum perforation location is at the top of the formation whilst the length should be half of the oil zone thickness without considering the effect of perforation on skin effects. Pasztor and Toth (2017) investigated the effect of perforation parameters (shot density, perforation length, perforation radius and phase angle) on good productivity; however, they did not consider how these parameters could be used to determine the critical flow rate.
Three fundamental elements have been identified to form the basis for an optimized perforation strategy critical to maximizing hydrocarbon recovery. Firstly, perforations must be properly oriented via optimal perforation-gun orientation; second, debris from the perforation tunnels must be effectively removed via dynamic underbalanced perforating; and third, formation damage must be minimized through proper wellbore fluid selection during perforating (Behrmann et al. 2007; Bellarby 2009; Onuh et al. 2019). In addition, perforation done closer to the water or gas layer has the possibility of producing water within a short time and thereby reducing hydrocarbon productivity.
In view of the fact that coning can have an important influence on operations, recovery and economics, it is of great value to optimize perforation length, position and interval to minimize both the skin factor and coning of water into the completion interval. This work investigates all the key perforation parameters as a function of the skin effect in critical flow rate determination. It also presents a workflow/protocol to predict the maximum water-free oil production rate. The results obtained from a computational model developed to obtain the maximum water-free production rate in terms of perforation length, radius, spacing between perforations and phasing angle are presented in this paper. The mathematical details of the equations used and the presented results are in appendices.
Theoretical basis
Meyer and Garder (1954) developed a correlation for estimating the critical oil rate in the water coning system as shown schematically in Fig. 1.
A homogeneous system with a uniform permeability throughout the reservoir was assumed. The correlation shows that the critical oil rate is directly proportional to the difference between the water and oil densities, effective oil permeability and inversely proportional to oil viscosity (Eq. 1). If the total skin factor, Sdp, is added to the equation then critical flow rate becomes
The skin factor as modelled by Karakas and Tariq (1991) captured perforation parameters and how they contribute to the skin damage (this is presented in “Appendix” section) and the critical flow rate. To determine either overall productivity and/or skin from the perforated interval, the performance of a single perforation must be considered. The overall perforation performance can then be evaluated by combining the phasing and shots per foot of multiple adjacent perforations. A schematic of a perforated wellbore showing the altered near wellbore region is shown in Fig. 2.
Computational model
The Meyer and Garder critical coning correlation for the vertical well is employed in the computational model because it has been found to be the most conservative correlation from the literature (Joshi 1991; Ahmed 2019). Firstly, the computational model computes the overall skin from perforation within and beyond damage zones as a function of near wellbore parameters which include perforation length and radius, phasing angle and shot density. Secondly, the skin value obtained is incorporated into the maximum water-free oil production rate correlation. This implies that the optimum length, radius of perforation, phasing angle and shot density required to obtain maximum water-free oil production rate can be estimated prior to the drilling and perforating of the interval to be completed. The computational model flow sheet is presented in Fig. 3. The computational model was programmed in C++.
The related well and reservoir data used for the computation are given in Table 3 in “Appendix”; however, for the purpose of analysis and discussion, few results were selected and plotted in Figs. 4, 5, 6, 7 and 8 at oil viscosity, μo = 0.73 cp; formation volume factor, Bo = 1.0 bbl/STB and water density, ρw = 63.76 lb/ft3.
For the effect of skin, we considered a vertical well drilled into an oil reservoir with perforation parameters of spacing between perforations, 2 in (6 shots/ft); perforation radius, rp = 0.16 in; phasing = 45 degrees; oil viscosity, μo = 1 cp; formation volume factor, Bo = 1.1 bbl/STB and water density, ρw = 62.4 lb/ft3. Considering the perforations that terminate within and beyond the damaged zones, the overall skin factor varies with perforation length as shown in Table 1.
There was a large skin factor when the perforation was done within the damaged zone; however, when the perforation length was increased by 2 in, a drastic reduction of 6.2 from skin factor was recorded. As perforation length was increased further and crossed a little bit beyond damaged zone, there was a further reduction of 0.9 in skin factor. Generally, the combined effect of perforation done beyond the damaged zone and increase in perforation length was evident as the overall skin decreases with increasing perforation length. Hence, the effect of phasing, perforation radius and other perforation parameters on overall skin can be evaluated prior to drilling and completion operations. The optimal combination of perforation parameters that minimizes skin can be obtained from the computational model.
Using the modified Mayer and Garder correlation (Eq. 1), the influence of skin factor on the maximum water-free oil production rate is shown in Table 1. The loss in productivity due to overall skin is evident in Table 1. This demonstrates how perforation parameters can add to the effect of skin factors on flow rate. Thus, the completion operation could be carried out with optimal perforation length as well as other perforation parameters to minimize the skin and maximize the water-free oil production rate.
Results and discussion
Effect of phasing angles
Critical rate increases with phasing angle in the order of 0°, 180°, 90°, 120°, 45° and 60° for a perforation length of 6 in and for the spacing between perforation of 1 to 3 in (Fig. 4). The same pattern was recorded for a perforation length of 18 in and 30 in but for the spacing between perforation above 1.2 in and above 2.4 in, respectively (Figs. 5, 6). It was observed that both the spacing between perforations and perforation length have influence on the phasing angle yielding an optimum critical rate between 45° and 60°. Hence, the above phasing angle order changed to 0°, 180°, 90°, 120°, 60° and 45° from perforation length of 18 in except for perforation length in favour of 60° (Figs. 5, 6, 7, 8). Hence, the only condition that will make a 60° phasing angle a better option is when the perforation length is shorter with wider perforation spacing.
Effect of spacing between perforations
The impact of a decrease in spacing between perforations (which is an increase in perforation density) on the critical flow rate is small. For instance, from Fig. 6 at a perforation length of 30 in, perforation radius of 0.36 in and at a phasing angle of 45°, the critical rate recorded is 21.29 bbl/day, 21. 16 bbl/day, 20. 80 bbl/day, 20.29 bbl/day and 19.87 bbl/day at spacing between perforation of 1 in, 1.2 in, 1.7 in, 2.4, and 3 in respectively. There was just 1.42 bbl/day of oil that will be added for decreasing the spacing between perforations from 3 to 1 in. Hence, the need to consider the cost of perforation to determine optimal perforation density to aid cost-effective completion.
Effect of perforation radius
Increasing the perforation radius from 0.12 in to 0.36 in has a minimum and maximum increment on critical flow rate to be 0.3 bbl/day and 4.25 bbl/day. These minimum and maximum are recorded at a perforation length of 48 in, the spacing between perforations of 1 in and a phasing angle of 0° and the maximum effect at a perforation length of 6 in, the spacing between perforations of 3 in and a phasing angle of 60°, respectively. However, as the perforation length increases the effect of perforation radius reduces as the gradient reduces/flattens out (Figs. 4, 5, 6, 7, 8). It implies that the effect of an increase in perforation length overshadows the effect of an increase in perforation radius. For example, at perforation length of 6 in, the spacing between perforations of 1 in and at phasing angle of 180° there is an increment of 2.33 bbl/day critical flow rate by increasing the perforation radius from 0.12 in to 0.36 in whilst the increment has reduced to 0.45 bbl/day at perforation length of 48 in (Figs. 4i and 8i). Similarly, the increment reduces from 3.97 bbl/day to 1.48 bbl/day when the perforation length was increased from 6 to 48 in at spacing between perforations of 3 in and at a phasing angle of 180° (Figs. 4v and 8v).
Effect of perforation length
An increase in perforation length increases the critical flow rate, however, the effect reduces drastically as the perforation length reached 30 in (Figs. 4, 5, 6, 7, 8). For instance, Figs. 4i to 8i show that at a phasing angle of 120°, the spacing between perforations of 1 in and perforation length of 6 in, 18 in, 30 in, 42 in and 48 there are critical flow rate increment percentage of 35%, 9.1%, 5.95% and 2.4%. Also, Figs. 4v to 8v show that at phasing angle of 120° but the spacing between perforations of 3 in and perforation length of 6 in, 18 in, 30 in, 42 in and 48 there is a critical flow rate increment percentage of 59%, 12.4%, 7.3% and 2.8%.
Similar observations were made by Pasztor and Toth (2017), when the authors investigated sensitivity tests of perforation parameters on the production rate.
Conclusion
A computation model that predicts the maximum water-free oil production rate has been presented. This model incorporates the skin factor as a function of the perforation parameters. The model will aid engineers in obtaining the optimal combination of perforation parameters that minimize the total skin factor. The critical flow rate generally increases with an increase in perforation radius for all phasing angles considered, although the gradient reduces/flatten out as the perforation length increases. The effect of perforation radius reduces with an increase in the perforation length. The model shows that the optimal perforation angles are 45° and 60° with a marginal difference in rates, and 0° is the worst perforation angle. Considering the input data used for the computation in this work, the perforation parameters that produced maximum water-free oil rate of 23.2 std/d are 45°, 1 in, 0.36 in and 48 in of phasing angle, spacing between perforation, perforation radius and perforation length, respectively.
References
Abass HH, Bass DM (1988) The critical production rate in water-coning system. Society of Petroleum Engineers. https://doi.org/10.2118/17311-MS
Ahmed T (2019) Gas and Water Coning. In: Ahmed T (ed) Reservoir engineering handbook, 3 and 5th edn. Gulf Professional Publishing, pp. 605–661. https://doi.org/10.1016/B978-0-12-813649-2.00009-8.
Behrmann L, Walton IC, Chang FF, Fayard A, Khong CK, Langseth B, Mason S, Mathisen A, Pizzolante I, Xiang T, Svanes G (2007) Optimal fluid systems for perforating. Oilfield Rev 19:14–25
Bellarby J (2009) Well completion design. Elsevier Science, Amsterdam
Chierici GL, Ciucci GM, Pizzi G (1964) A systematic study of gas and water coning by potentiometric models. J Pet Technol 16:923–929. https://doi.org/10.2118/871-PA
Høyland LA, Papatzacos P, Skjaeveland SM (1989) Critical rate for water coning: correlation and analytical solution. Society of Petroleum Engineers, London. https://doi.org/10.2118/15855-PA
Joshi S (1991) Horizontal well technology. Pennwell Publishing Company, Tulsa
Karakas M, Tariq SM (1991) Semianalytical productivity models for perforated completions. Society of Petroleum Engineers, London. https://doi.org/10.2118/18247-PA
Locke S (1981) An advanced method for predicting the productivity ratio of a perforated well. Society of Petroleum Engineers, London. https://doi.org/10.2118/8804-PA
Li Y, Li H, Li Y (2015) Prediction method of bottom water coning profile and water breakthrough time in bottom water reservoir without barrier. Math Prob Eng. https://doi.org/10.1155/2015/149490
Luo X, Jiang L, Su Y, Huang K (2015) The productivity calculation model of perforated horizontal well and optimization of inflow profile. Petroleum 1:154–157. https://doi.org/10.1016/j.petlm.2015.04.002
Menouar HK, Hakim AA (1995) Water coning and critical rates in vertical and horizontal wells. Society of Petroleum Engineers, London. https://doi.org/10.2118/29877-MS
Meyer HI, Garder AO (1954) mechanics of two immiscible fluids in porous media. J Appl Phys 25:1400–1406. https://doi.org/10.1063/1.1721576
Onuh CY, Ogunkunle TF, Oluwatosin JR, Efeoghene E (2019) The impact of perforation geometry on oil well productivity. IJMET 10:1012–1021
Papatzacos P, Herring TR, Martinsen R, Skjaeveland SM (1989) Cone breakthrough time for horizontal wells. Society of Petroleum Engineers, London. https://doi.org/10.2118/19822-PA
Pasztor AV, Toth A (2017) Effect of perforation parameters on the productivity of geothermal wells. In: The publications of the multiscience—XXXI. MicroCAD international scientific conference. https://doi.org/10.26649/musci.2017.002
Tabatabaei M, Ghalambor A, Guo B (2012) An analytical solution for water coning in vertical wells. Society of Petroleum Engineers, London. https://doi.org/10.2118/113106-PA
Yang W, Wattenbarger RA (1991) Water coning calculations for vertical and horizontal wells. Society of Petroleum Engineers, London. https://doi.org/10.2118/22931-MS
Zhang L, Zhao Y, Liu Z (2011) A novel steady-state productivity equation for horizontal wells in bottom water drive gas reservoirs. Petrol Sci 8:63–69. https://doi.org/10.1007/s12182-011-0116-2
Zhang HJ, Yang ZZ, Guan PJ, Xu GP, Fu L, Li HY, Xu L (2012) Sectional optimization study on perforation parameters of horizontal wells in Daqing oilfield. Adv Mater Res 594–597:2586–2589. https://doi.org/10.4028/www.scientific.net/amr.594-597.2586
Funding
This research received no specific grant from any funding agency in the public, commercial, or not-for-profit sectors.
Author information
Authors and Affiliations
Corresponding author
Additional information
Publisher's Note
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
Appendix: Perforation skin (S p)
Appendix: Perforation skin (S p)
A semi-analytical perforation skin model for calculating the overall skin due to drilling damage and perforation is expressed as Eq. 3 (Karakas and Tariq 1991). The perforation skin (Sp) excluding the drilling damage skin is the sum of the horizontal skin (Sh), wellbore skin (Swb), vertical skin (Sv) and the crushed zone skin (Sc):
where
The constants α, α1, α2, b1, b2, C1, C2, are obtainable based on gun-phasing angle in Table 2.
The calculation of overall skin varies depending on whether or not the perforation terminates within the drilling damaged zone. For perforations terminating within the damaged zone (i.e. lp < ld) the calculation of the overall skin which is the combination of drilling damage and perforation (Sdp) is given as follows:
The parameter \(S_{\aleph }\) is used for boundary effect correction and is often ignored. Computation of the overall skin in this paper is based on data in Table 3.
The overall skin (Sdp) for the more relevant case of perforations that extend beyond the damage zone is expressed as follows:
Meanwhile, the perforation length and wellbore radius are modified as follows:
\(l_{\rm p}^{{\prime}}\) and \(r_{\rm w}^{{\prime}}\) are used instead of lp and rw in Eqs. (6) and (9). The modified perforation length and wellbore radius also affect Swb, Sv, and Sc to become modified (\(S_{\rm wb}^{{\prime}} , S_{\rm v}^{{\prime}} , S_{\rm c}^{{\prime}} )\) thus:
Meanwhile, Sh is partly calculated using rw, such that \(S_{\rm h }^{{\prime}}\) is given by
Rights and permissions
Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if changes were made. The images or other third party material in this article are included in the article's Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article's Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this licence, visit http://creativecommons.org/licenses/by/4.0/.
About this article
Cite this article
Adeleye, J.O., Olamigoke, O. & Mumuni, O.T. A predictive protocol to obtain maximum water-free oil production rate for perforated vertical wells. J Petrol Explor Prod Technol 11, 347–357 (2021). https://doi.org/10.1007/s13202-020-01014-z
Received:
Accepted:
Published:
Issue Date:
DOI: https://doi.org/10.1007/s13202-020-01014-z