Research Article  Open Access
Yi Huang, Jin Yang, Lingyu Meng, Xuyue Chen, Ming Luo, Wentuo Li, "Numerical Investigation on Gas Accumulation and Gas Migration in the Wavy Horizontal Sections of Horizontal Gas Wells", Mathematical Problems in Engineering, vol. 2020, Article ID 7275209, 9 pages, 2020. https://doi.org/10.1155/2020/7275209
Numerical Investigation on Gas Accumulation and Gas Migration in the Wavy Horizontal Sections of Horizontal Gas Wells
Abstract
Wavy horizontal sections are typically encountered in horizontal gas wells, which will result in gas accumulation on top of the wavy horizontal sections. This gas accumulation can be a problem and may trigger gas kick or blowout accident while tripping and pulling this gas into the vertical section. In this paper, a numerical model for gas accumulation and gas migration in the wavy horizontal sections of the horizontal gas well is developed; meanwhile, the gas accumulation and gas migration process is numerically investigated. The results show that the gas exhausting time in the wavy horizontal section increases with the increase of the wellbore curvature and the critical drilling fluid flow velocity for gas exhausting increases with the increase of the wellbore curvature. When the drilling fluid flow velocity is higher than the critical drilling fluid flow velocity for gas exhausting, no gas accumulation will occur. With all other parameter values set constant, the number of the wavy horizontal sections has a great effect on the gasliquid flow pattern while it has little effect on the efficiency of the gas exhausting. This work provides drilling engineers with a practical tool for designing the drilling fluid flow velocity to avoid gas kick or blowout accident in horizontal gas well drilling.
1. Introduction
Horizontal wells are widely used in petroleum and natural gas development, and they have many advantages over traditional vertical wells, such as increased drainage area and high production [1–3]. However, since the wellbore trajectory in horizontal drilling is difficult to control, the horizontal section is always not completely horizontal, and sometimes wavy horizontal sections are formed. The gas in the wavy horizontal sections cannot migrate in the direction of flow owing to buoyancy; it results in pockets of gas accumulation (see Figure 1). The gas accumulation can be a problem and may trigger gas kick or blowout accident while tripping and pulling this gas into the vertical section [4–6].
In the past decades, several key studies have been conducted on horizontal well control and gas migration. Vefring et al. established new models for the gas slip and rise velocities in near horizontal wells and various models for different gas removal mechanisms [7]. Chexal et al. proposed a comprehensive drift flux model [8], but the model is complicated for field applications because the fluid distribution is formed by the correlation of multiple empirical curve fitting parameters through the distribution parameters. Hibiki and Ishii proposed a drift velocity equation; it is applied to slug flow [9]. Gao et al. simulated the storage and removal processes of the gas slug by experiments and analyzed the migration of the gas slug [10]. In recent years, the drift speed equation proposed by Woldesemayat and Ghajar considers the influences of surface tension and pipeline diameter on drift speed apart from the influences of pipeline direction and system pressure [11]. Wang et al. analyzed the variation of wellbore pressure along the depth of the well during the time of gas kick in a horizontal well [12]. However, previous studies are primarily for completely horizontal sections, there are few research studies on the gas accumulation and gas migration in the wavy horizontal sections of the horizontal gas well, and few people used numerical simulation to study them. In this paper, a numerical model for gas accumulation and gas migration in the wavy horizontal sections of the horizontal gas well is developed; meanwhile, the gas accumulation and gas migration process is numerically investigated.
2. Physical Model and Governing Equations
2.1. Physical Model
The diameter of the borehole in the wavy horizontal sections of the horizontal gas well is 0.2159 m. Gas is accumulated on the top of the wavy horizontal sections. The schematic diagram of the gas accumulation in the wavy horizontal sections of the horizontal gas well is shown in Figure 2.
Establishing a reasonable twodimensional physical model of the wavy horizontal sections of the horizontal gas well according to the actual situation is the key to the simulation of gas accumulation and gas migration. In order to reduce the calculation amount and improve the simulation efficiency, only the fluid domain is established. The borehole diameter is 0.2159 m, and the horizontal distance on both sides of the model is 1 m. A horizontal wellbore physical model with different wavy horizontal sections is established, as shown in Figure 2.
ICEM CFD is selected as the meshing software, and quadrilateral meshing is adopted [13, 14]. The ratio of the length and width of the mesh is not greater than 5. According to the characteristics of fluid motion, the wall grid is locally encrypted to ensure accurate simulation of gasliquid twophase characteristics, as shown in Figures 3 and 4.
(a)
(b)
(c)
2.2. Governing Equations
The gasliquid flow process in the wavy horizontal sections of a horizontal gas well is unstable. Even if the boundary conditions remain the same, the physical quantity of the fluid during the flow still has strong pulsation, so the simulation model should be run in a turbulent state [15–18]. By numerical analysis method and simulating it, the results can be compared with the actual situation. Considering the gravity factor, VOF model and RNG k turbulence model are selected to simulate the flow changes under different working conditions. The governing equation is established as follows:(1)Continuity equation:(2)Momentum equation:(3)Turbulence equation: Turbulent energy equation k is(4) Dissipation rate equation is
In the expression, is the fluid velocity, m/s; is gravitational acceleration, m/s^{2}; are the spatial coordinates; is the fluid density, kg/m^{3}; is static pressure, Pa; is the fluid viscosity, Pa·s; is volume force, N; is the effective fluid viscosity, Pa·s; is time, s; is turbulent kinetic energy, J; is the turbulent energy dissipation; are constants.
2.3. Physical Parameter
Gasliquid density and dynamic viscosity are shown in Table 1.

2.4. Equation Discretization and Solving Method
An unsteady, implicit separation, and solving algorithm is used. The governing equations to be solved are continuity equations that satisfy mass conservation, momentum conservation, energy conservation, momentum equations, energy equations, and turbulence equations that take turbulence properties into account [19]. The finite volume method is used to discretize the governing equations, and a suitable discretization format is selected. The pressure interpolation format selects the physical strength weighting format; the interpolation formats of density, momentum, turbulence energy, turbulence dissipation rate, and energy select the firstorder upwind style, high stability, and fast calculation speed, and volume fraction interpolation format selects the geometric reconstruction format; pressurespeed coupling algorithm selects PISO algorithm [20–22].
2.5. Setting of Definite Conditions
Definite solution conditions consist of a combination of boundary conditions and initial conditions.(1)Inlet conditions: velocity inlet is used, and the inlet speed is set to 1.2 m/s, 1.6 m/s, and 2 m/s. The inlet boundary is set to a gas inlet volume fraction of 0. The turbulence definition method selects the hydraulic diameter and turbulence intensity.(2)Outlet conditions: the pressure outlet and turbulence definition method are used to select the hydraulic diameter and turbulence intensity.(3)Initial conditions: considering the influence of gravity, set Y = −22129.81 m/s^{2} under the fluent environment panel. In the initialization panel, select inlet initialization and then perform partial repair on the initial volume fraction of the gas in the Patch panel to achieve the setting of the initial concentration of gas in the model. It is set to record a gasliquid composition cloud chart every 100 time steps. After the simulation is completed, an animation can be formed to observe the flow pattern.
2.6. Convergence Conditions
Fluent software uses residuals to reflect the convergence of the calculations and judges whether the iterative process converges through the final iterative residual output of each equation for each iteration step [23]. The calculation process ends when the residuals of each equation reach the set convergence criteria. In this model, all residuals are set to 1.0 × 10^{−4}, the time step is set to 0.01 s, and the maximum iteration step is set to 20 steps [24–28]. The total time steps are specifically set according to the actual flow of gasliquid twophase tube flow.
3. Simulation Results and Analysis
There are various combinations of gas accumulation simulations in the wavy horizontal sections of horizontal gas wells. Here, simulation studies of gas accumulation in single wavy horizontal section and complex wavy horizontal sections are performed.
3.1. Simulation Analysis of Single Wavy Horizontal Section
3.1.1. Influence of Curvature
According to the actual working conditions, wavy horizontal sections with different curvatures are selected, and models of gas accumulation of wavy horizontal sections with different curvatures are established to simulate gasliquid twophase flow. The diameter of the wavy horizontal section is 0.2159 m, the total horizontal length is 5 m, the drilling fluid flow velocity is 1.6 m/s, the initial gas accumulation at the top is 0.3 m^{3}, and the curvature of the wavy horizontal section is 0.2, 0.3, 0.4, and 0.5. As displayed, red is drilling fluid and blue is gas.
It can be seen from Figure 5 that, at a speed of 1.6 m/s and different curvatures of 0.2, 0.3, 0.4, and 0.5, the gas in the wavy horizontal section of the horizontal gas well is exhausted, and there are no phenomenon of gas accumulation and inability of exhaust. The gas pocket formed by the accumulated gas is carried forward by the drilling fluid. When the gas pocket enters the downward section, the gas pocket cannot be exhausted directly. The front end of the gas pocket is broken into small bubbles by the shear force, and wavy horizontal section of the horizontal gas well is exhausted in the form of small bubbles. With time, the gas pocket is gradually becoming smaller and the drilling fluid pushes the bladder out of the wavy horizontal section of the horizontal gas well. As the curvature increases, the time it takes to exhaust the gas increases. Therefore, the lower the curvature of the horizontal wellbore, the less likely the gas accumulation will occur.
(a)
(b)
(c)
(d)
3.1.2. Influence of Flow Rate
The borehole diameter is 0.2159 m, the total horizontal length is 5 m, the initial gas accumulation at the top is 0.3 m^{3}, the curvature is 0.3, and the drilling fluid flow velocity is 1.2 m/s, 1.6 m/s, and 2 m/s, as shown in Figure 5.
Observing Figure 6, it can be seen that when the drilling fluid flow velocity is 2 m/s, the drilling fluid directly carries the entire bladder to discharge the wavy horizontal section without shear fracture. When the drilling fluid flow velocity is 1.6 m/s, the gas formed by the accumulated gas is carried by the drilling fluid. Forward, when the gas pocket enters the downward section, the gas pocket cannot be exhausted directly, the front end of the gas pocket is broken into small bubbles by the shear force, the gas is exhausted in the form of small bubbles, and the drilling fluid is carried out by the remaining gas pocket. When the drilling fluid flow velocity is 1.2 m/s, the drilling fluid pushes the gas pocket to the downward section, the gas pocket is stationary, and the front section of the gas pocket is broken into small bubbles by the shear force to discharge the gas from the wavy horizontal section. Over time, part of the gas cannot be broken into small bubbles and stays in the wavy horizontal section, causing gas accumulation phenomenon.
(a)
(b)
(c)
From the above phenomenon, it can be seen that there is a critical flow rate to make the gas just exit the wavy horizontal section of the horizontal gas well. When the drilling fluid flow velocity is higher than the critical velocity, gas accumulation will not occur; when the drilling fluid flow velocity is lower than the critical velocity, gas retention will cause gas accumulation. At this time, the gas can only be exhausted by dissolving in the drilling fluid and increasing the speed of the drill pipe. In order to prevent the accumulation of gas in the wavy horizontal section, the critical flow rate needs to be determined. Therefore, this paper performed a series of gasliquid simulation of the wavy horizontal section of the horizontal gas well with different curvatures. The simulation results are shown in Table 2.
 
Y = exhaust; F = not exhausted. 
As known from Table 2, when the curvature is 0.2, the critical drilling fluid flow velocity is 0.8 m/s; when the curvature is 0.3, the critical drilling fluid flow velocity is 1.0 m/s; when the curvature is 0.4, the critical drilling fluid flow velocity is 1.2 m/s; when the curvature is 0.5, the critical drilling fluid flow velocity is 1.3 m/s. The critical drilling fluid flow velocity increases with the curvature.
3.1.3. Experimental Comparison and Verification
Gao et al. [10] analyzed the gas migration process in the undulating section of a horizontal well through experiments. In order to verify the correctness of the simulation results, the simulation results were compared with their experimental results. The gas traps in the elbows mainly bear the frictional resistance of the pipe wall in the initial state. When the bubbles enter the downdip sections, it is difficult for gas to discharge because of the buoyancy, i.e., the flow resistance. The front ends of the bubbles are crushed and separated under the liquid phase flow disturbance conditions, but only few bubbles are separated and migrate mostly in large bubble form through mainly surface tension and liquid phase friction.
The simulation results are exactly the same as the experimental results, verifying the accuracy of the simulation results.
3.2. Simulation Analysis of Complex Wavy Horizontal Sections
Horizontal gas wells may have multiple wavy horizontal sections, so the total horizontal length is 20 m, the borehole diameter is 0.2159 m, the curvature is 0.2, the initial top gas accumulation is 0.98 m^{3}, and the critical drilling fluid flow velocity is 1 m/s, 1.2 m/s, 1.6 m/s, 2 m/s, and 1.5 times wavy horizontal sections and 2 times wavy horizontal sections are simulated. The simulation results are shown in Figures 6 and 7.
(a)
(b)
(c)
(d)
Observing Figure 7, it can be seen that when the drilling fluid flow velocity is 2 m/s, the drilling fluid carries gas into the downward section, and the front end of the gas pocket is sheared and broken. The remaining part is still carried by the drilling fluid to the second upward section in the form of a large gas pocket and then exhausted; when the drilling fluid flow velocity is 1 m/s, 1.2 m/s, and 1.6 m/s, the drilling fluid carries gas into the downward section, and the front end of the gas pocket shears and breaks. When the gas pocket reaches the bottom of the wavy horizontal section, it will stay and can only be exhausted by breaking into small bubbles, and, as the speed decreases, the longer the crushing time, the smaller the broken bubbles.
Observing Figure 8 shows that when the drilling fluid flow velocity is 2 m/s, the gas pocket is carried by the drilling fluid to the bottom of the wavy horizontal section, and then it stays, and then it is sheared and broken into large bubbles and migrates to the second wavy horizontal section. There is no gas accumulation at the top and it is directly exhausted; when the drilling fluid flow velocity is 1 m/s, 1.2 m/s, and 1.6 m/s, the gas pocket is carried by the drilling fluid into the downward section, the front end of the gas pocket is sheared and broken, and the gas pocket is brought by the drilling fluid to the wavy horizontal sections. At the bottom of the section, it stays and breaks into medium and small bubbles and enters the second wavy section. At the top of the second wavy section, it begins to gather into an gas pocket, and then the front end of the gas pocket is sheared and exhausted.
(a)
(b)
(c)
(d)
Observing Figures 9 and 10, we can see that when the flow velocity is the same, the curvature and the total length of wavy horizontal sections are the same, and the same waveform; the exhaust time of 1.5 wavy horizontal sections is the same as that of two wavy horizontal sections, and the curve is almost the same; therefore, the number of wavy horizontal sections has little effect on the efficiency of the exhaust.
4. Conclusions
(1)In the case of the same length and the same flow velocity, the gas exhaust time increases with the increase of the curvature. Therefore, the lower the curvature of the wavy horizontal sections, the less likely it is to generate gas.(2)When the drilling fluid flow velocity is extremely large, the drilling fluid directly carries the entire gas pocket out of the wavy horizontal section without shear fracture; as the inlet flow velocity decreases, the gas pocket formed by the accumulated gas is carried forward by the drilling fluid. When the gas pocket enters the downward section, the front end of the gas pocket is broken into small bubbles by the shear force. The gas is exhausted in the form of small bubbles. When the drilling fluid flow velocity is reduced to a certain speed, the drilling fluid pushes the gas pocket downward section. Then, the gas pocket is stationary, and the front section of the gas pocket is broken into small bubbles by the shear force. Over time, some of the gas cannot be broken into small bubbles and stays in the wavy horizontal section, causing gas accumulation.(3)There is a critical flow velocity so that the gas is just completely exhausted from the wavy horizontal section. When the drilling fluid flow velocity is higher than the critical velocity, gas accumulation will not occur; when the drilling fluid flow velocity is lower than the critical velocity, gas retention will cause gas accumulation. And the critical drilling fluid flow velocity increases with the curvature.(4)When the total length and curvature of the wavy horizontal sections are the same, the number of wavy horizontal sections has a great effect on the gasliquid flow pattern but has little effect on the efficiency of the exhaust.
Data Availability
The data used to support the findings of this study are available from the corresponding author upon request.
Conflicts of Interest
The authors declare that they have no conflicts of interest.
Acknowledgments
The authors gratefully acknowledge the financial support by “Thirteenth FiveYear Plan” China National Offshore Oil Corporation (CNOOCKJ135ZDXM24LTDZJ01) and the National Science and Technology Major Project (2017ZX05009003).
References
 X. Zhang, J. Xie, and B. Yu, “Numerical simulation on the gasliquid twophase flow in the direct commissioning process of a hilly pipeline,” Journal of the University of Chinese Academy of Sciences, vol. 34, no. 2, pp. 265–272, 2017. View at: Google Scholar
 E. H. Vefring, “An advanced kick simulator for high angle and horizontal wellspart I,” in SPE/IADC Drilling Conference, Society of Petroleum Engineers, Amsterdam, Netherlands, March 1995. View at: Google Scholar
 N. Wei and Z. Cui, “The rule of carrying cuttings in horizontal well drilling of marine natural gas hydrate,” Energies, vol. 13, no. 5, p. 1129, 2020. View at: Publisher Site  Google Scholar
 N. Liu, W. Sun, Y. Meng et al., “Multiphase non equilibrium pipe flow behaviors in the solid fluidization exploitation of marine natural gas hydrate reservoir,” Energy Science & Engineering, vol. 6, no. 6, pp. 760–782, 2018. View at: Publisher Site  Google Scholar
 N. Wei, W.T. Sun, Y.F. Meng et al., “Change mechanism of transient gasliquid twophase flow in wellbore during marine natural gas hydrate reservoir drilling,” Thermal Science, vol. 23, no. 4, pp. 2179–2187, 2019. View at: Publisher Site  Google Scholar
 L. Yongwang, “Discussion on the trajectory control technology of horizontal section of horizontal well,” Science Technology and Engineering, vol. 11, no. 35, pp. 8872–8875+8881, 2011. View at: Google Scholar
 E. H. Vefring, Z. Wang, and R. Rommetveit, “An advanced kick simulator for high angle and horizontal wellspart II,” in Middle East Oil Show, Society of Petroleum Engineers, Manama, Bahrain, March 1995. View at: Google Scholar
 B. Chexal, G. Lellouche, J. Horowitz, and J. Healzer, “A void fraction correlation for generalized applications,” Progress in Nuclear Energy, vol. 27, no. 4, pp. 255–295, 1992. View at: Publisher Site  Google Scholar
 T. Hibiki and M. Ishii, “Onedimensional driftflux model and constitutive equations for relative motion between phases in various twophase flow regimes,” International Journal of Heat and Mass Transfer, vol. 46, no. 25, pp. 4935–4948, 2003. View at: Publisher Site  Google Scholar
 Y. Gao, X. Sun, T. Zhao, Z. Wang, X. Zhao, and B. Sun, “Study on the migration of gas kicks in undulating sections of horizontal wells,” International Journal of Heat and Mass Transfer, vol. 127, pp. 1161–1167, 2018. View at: Publisher Site  Google Scholar
 M. A. Woldesemayat and A. J. Ghajar, “Comparison of void fraction correlations for different flow patterns in horizontal and upward inclined pipes,” International Journal of Multiphase Flow, vol. 33, no. 4, pp. 347–370, 2007. View at: Publisher Site  Google Scholar
 B. Kang, H. Fan, P. Jiang et al., “Simulation and experiment on phase equilibrium of gas hydrate using the ttype pipe confluence model,” Mathematical Problems in Engineering, vol. 2020, Article ID 8932736, pp. 1–11, 2020. View at: Publisher Site  Google Scholar
 Y. Yu, Fluent Introduction and Advanced Course, Beijing Institute of Taechnology Press, Beijing, China, 2008.
 Q. Ma, “Simulation of slack line phenomena in the big drop pipeline under the different conditions,” Journal of Liaoning University of Petroleum & Chemical Technology, vol. 35, no. 1, pp. 37–40, 2015. View at: Google Scholar
 S. A. Morsi and A. J. Alexander, “An investigation of particle trajectories in twophase flow systems,” Journal of Fluid Mechanics, vol. 55, no. 02, pp. 193–208, 1972. View at: Publisher Site  Google Scholar
 Q. Li, “Numerical simulation of oil and water twophase dispersed flow in horizontal pipe based on fluent,” Oil and Gas Field Surface Engineering, vol. 09, pp. 6869, 2013. View at: Publisher Site  Google Scholar
 C. Chen, Threedimensional Simulation and Experimental Study of GasLiquid Twophase Flow in Vertical Riser. Diss, Tianjin University, Tianjin, China, 2009.
 Z. Shen, “Gasliquid twophase flow in pipelines,” Petrochemical Technology, vol. 25, no. 7, p. 175, 2018. View at: Google Scholar
 N. Wei, C. Xu, Y. Meng, G. Li, X. Ma, and A. Liu, “Numerical simulation of gasliquid twophase flow in wellbore based on drift flux model,” Applied Mathematics and Computation, vol. 338, pp. 175–191, 2018. View at: Publisher Site  Google Scholar
 S. Liu and Z. Zhu, “Application of composite deflecting model in horizontal well drilling,” Mathematical Problems in Engineering, vol. 2020, pp. 1–10, 2020. View at: Publisher Site  Google Scholar
 X. Li, Z. Han, S. Yang, and G. Chen, “Underwater gas release modeling and verification analysis,” Process Safety and Environmental Protection, vol. 137, pp. 8–14, 2020. View at: Publisher Site  Google Scholar
 S. J. Perkins and H. A. Li, “Elongated bubble centring in horizontal gasliquid slug flow,” International Journal of Multiphase Flow, vol. 123, Article ID 103158, 2020. View at: Publisher Site  Google Scholar
 B. Sun, X. Sun, Z. Wang, and Y. Chen, “Effects of phase transition on gas kick migration in deepwater horizontal drilling,” Journal of Natural Gas Science and Engineering, vol. 46, pp. 710–729, 2017. View at: Publisher Site  Google Scholar
 J.C. Feng, X.S. Li, G. Li, B. Li, Z.Y. Chen, and Y. Wang, “Numerical investigation of hydrate dissociation performance in the south China sea with different horizontal well configurations,” Energies, vol. 7, no. 8, pp. 4813–4834, 2014. View at: Publisher Site  Google Scholar
 P. V. Godbole, C. C. Tang, and A. J. Ghajar, “Comparison of void fraction correlations for different flow patterns in upward vertical twophase flow,” Heat Transfer Engineering, vol. 32, no. 10, pp. 843–860, 2011. View at: Publisher Site  Google Scholar
 A. J. Ghajar and C. C. Tang, “Void fraction and flow patterns of twophase flow in upward and downward vertical and horizontal pipes,” Advances in Multiphase Flow and Heat Transfer, vol. 4, pp. 175–201, 2012. View at: Google Scholar
 C. Y. Xu, “Experimental simulation and numerical modeling of dynamic variations in wellbore pressure during gaskicks,” Acta Petrolei Sinca, vol. 36, no. 1, pp. 120–126, 2015. View at: Google Scholar
 F. Ren, B. Wang, L. Zhao, and A. Zhu, “Experimental investigation and analysis of dynamic buckling of drill string in horizontal well,” Shock and Vibration, vol. 2017, Article ID 1658435, pp. 1–15, 2017. View at: Publisher Site  Google Scholar
Copyright
Copyright © 2020 Yi Huang et al. 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.