Research Article  Open Access
Xiangwei Kong, Yuanhua Lin, Yijie Qiu, Xing Qi, "A Novel Dynamic Model for Predicting Pressure Wave Velocity in FourPhase Fluid Flowing along the Drilling Annulus", Mathematical Problems in Engineering, vol. 2015, Article ID 134102, 17 pages, 2015. https://doi.org/10.1155/2015/134102
A Novel Dynamic Model for Predicting Pressure Wave Velocity in FourPhase Fluid Flowing along the Drilling Annulus
Abstract
A dynamic pressure wave velocity model is presented based on momentum equation, massbalance equation, equation of state, and small perturbation theory. Simultaneously, the drift model was used to analyze the flow characteristics of oil, gas, water, and drilling fluid multiphase flow. In addition, the dynamic model considers the gas dissolution, virtual mass force, drag force, and relative motion of the interphase as well. Finite difference and NewtonRaphson iterative are introduced to the numerical simulation of the dynamic model. The calculation results indicate that the wave velocity is more sensitive to the increase of gas influx rate than the increase of oil/water influx rate. Wave velocity decreases significantly with the increase of gas influx. Influenced by the pressure drop of fourphase fluid flowing along the annulus, wave velocity tends to increase with respect to well depth, contrary to the gradual reduction of gas void fraction at different depths with the increase of backpressure (BP). Analysis also found that the growth of angular frequency will lead to an increase of wave velocity at low range. Comparison with the calculation results without considering virtual mass force demonstrates that the calculated wave velocity is relatively bigger by using the presented model.
1. Introduction
In petroleum industry, managed pressure drilling (MPD) is considered to be one of the most important techniques, which allows accurate control of bottom hole pressure (BHP) by controlling the flow rate, drilling mud density, and back pressure (BP) at the wellhead [1]. As extensively used in the drilling of huge risk, uneconomical, or even abnormal formation, MPD has been of interest in the literature for at least the past decade, especially the related issues about pressure control inducing the dynamic well and kick detection [2]. Due to those efforts, research of the dynamic pressure wave velocity is of great significance to the detection of gas influx and effective control of the pressure at the bottom of well [3]. During the drilling operation of the socalled “microflux control,” an MPD technique developed by Santos et al. [4], the return flow is monitored and adjusted to control fluid loss or gain. In the light of control principal, simulation studies were performed to determine the most appropriate initial response to kicks arising due to MPD specific complications caused by BHP fluctuations [5]. During the MPD operations, all unsteady operating, such as changing of pumping rate, adjustment of choke, and controls of BP at wellhead, will generate a pressure wave and threaten the drilling equipment [6]. For the same reason, while tripping out of a drill string in the wellbore, bottom hole is submitted to a suddenly decrease in pressure, leading to fluid expansion and movement out of the annulus. The rapidly expanding fluids and dynamic pressure fluctuations can also lead to rock instability in a reservoir [7]. Particularly, the effects are more important in systems in which multiphase flows occur. New kickdetection tools are now available that are based on acoustic principles, which are of great benefit to potentially earlier and more sensitive detection of a gas influx than pitgain or paddle flow measurements [8]. The study of propagation of pressure wave is also relevant to control of downhole tool, such as intelligent well downhole control valves applied in different field for many purposes. With further development of oilfield, downhole tool technique for special casing wells is receiving much more attention [9]. At the meantime, an important problem in deep drilling is the propagation of measurementwhiledrilling (MWD) mud pulse, transmitting realtime various data from sensors located down hole near the drill bit. Hence, the propagation behavior of pressure wave is considered to provide reference for the MPD operations.
However, we also noticed that the conventional theories present are difficult to be employed in a systematic and accurate prediction when influxes generates. As the influxes fluid including gas, oil, and water will lead to variations of physic characteristics parameter of fluid in the annulus, the distribution of pressure wave velocity in gas, oil, water, and drilling mud fourphase flow along the annulus will be dynamic changing with respect to time and well depth. These effects may influence the safe operation of devices.
Pressure waves are disturbances that transmit energy and momentum along the wellbore through drilling fluid without significant displacement of matter. As migration of dispersed gas, water, and oil towards wellhead is quite complicated, the fundamental characteristics of the fourphase flow are still unknown and modeling results of pressure wave velocity are questionable. The proper treatment of propagation behavior of pressure wave in the fourphase flow in the annulus requires knowledge of description of influxes generation, development, and pressure wave propagation model in a twophase mixture.
In gas/liquid twophase flow, twophase medium interaction greatly changes the structure characteristic of flowing fluid, which results in greater compressibility of twocomponent flow than singlephase gas or liquid and further causes pressure waves propagating speed to be greatly reduced. This can be explained from the viewpoints of mixture density and compressibility of twophase fluid and the pressure drop along the flow direction in the wellbore. In the low void fraction range, the gas phase is dispersed in the liquid as bubble, so the wave velocity is influenced greatly by the added gas phase. According to the EOS, if gas invades into the wellbore with a small amount in the bottom hole, the density of the drilling mud has little variation while the compressibility increases obviously, which makes the wave velocity be decreased. This phenomenon was first proposed by Mallock [10] and attracts much attention for its important role in the development of science and technology applications. Extensive investigations involved in the issue of pressure wave propagation have been taken and some significant achievements associated with the theories have also been made in the researches. In early stage researching, Wood [11] extended the researches of Mallock and presented a succinct formula by assuming that the compressibility of twocomponent fluid is a function related to singlephase compressibility and elastic modulus E. Carstensen and Flody [12] measured the velocity of pressure wave under a lake. A dispersion relation for pressure waves propagating through a bubbly fluid was derived by using a linear scattering theory developed by Foldy. Under the hypothesis of homogeneous and adiabatic laminar flow, Thuraisingham [13] took the twophase media as a homogeneous fluid and derived the solution model concerning the problem of wave velocity model in twophase flow at low gas void fraction according to the analysis of state equation of mixture. According to the early stage investigations, Hsieh and Plesset [14] and Murray [15] researched the influence of thermal conductivity and viscosity coefficient on pressure wave propagation. Wallis [16] firstly studied the propagation mechanism of pressure wave and derived the propagation velocity in bubbly flow and separated flow using the homogeneous model. The proposed model is based on the respective compressibility of the vapor and the liquid. Also, the twophase mixture is treated as a compressible fluid with suitably averaged properties. In that expression, mass and heat transfers are neglected, so it can be applied to any gas/liquid mixture, in so far as no major effect due to vaporization or condensation must be taken into account. It may not be valid for such complex mixtures that never reach equilibrium. Assuming that no evaporation or condensation occurs when pressure wave is transporting, Moody [17] developed a simple acoustic wave model for bubbly flow and annular flow and established a relationship between sonic velocity and twophase critical flow. Similar model was developed by D’Arcy [18]. In the researches of Mcwilliam and Duggins [19], surface tension and compressibility of liquid phase were also considered. Henry et al. [20] calculated the velocity as a function of void fraction using a correlation to account for the change in bubble shape with void fraction. Martin and Padmanabhan [21] extended the simple model proposed by Henry by considering wave reflection and wave transmission at gasliquid interfaces. Researches of Mori et al. [22, 23] suggest that impact pipe elasticity on pressure wave propagation velocity is limited to the range of gas void of less than 1%. At high gas void range, the pressure wave velocity is in between the two velocities of singlephase. Nguyen et al. [24] proposed another type of model for bubbly flow (diluted gas phase in the liquid). The simple relations for prediction of the propagation of pressure disturbances in liquidgaseous twophase systems are presented. The model makes use of the wellknown physical behavior that the wave velocity in a singlephase fluid is influenced by the elasticity of the confining walls. The interface of the one phase is considered to act as the elastic wall of the other phase and vice versa. Mecredy and Hamilton [25] used the twofluid model to predict the pressure wave propagation in vaporliquid flow in detail. However, the analysis contained the important assumption that the evaporation or condensation was governed by kinetic theory. Michaelides and Zissis [26] developed a computational method which yields the sound velocity in terms of the thermodynamic coordinates of the substance without the use of diagrams. Corresponding velocities of sound for the four substances considered exhibit a certain similarity which is examined statistically. The relationship between the sound velocity and the critical mass flux is also investigated. Thuraisingham [13] studied the wave velocity in bubbly water at megahertz frequencies (1~10 MHZ). Numerical analytical results indicate that volume concentrations and the radius of the bubble relative to the incident wavelength of sound are the important parameters which determine the deviation of sound speed form that of bubblefree water.
Currently, the twofluid continuum model is the most common and reliable method to describe the gas/liquid twophase flow phenomenon [27]. In the model, the governing equation and phase interface relationship is established based on the assumption that each phase satisfies the continuum conditions. To obtain the practical flow equations, reasonable assumptions and constitutive equations should be introduced. In consequence, the predicted wave velocities were found to depend strongly on the introduced assumptions and equations. In recent years, the twofluid model was applied in determining the pressure wave propagation characteristics [28]. Ruggles et al. [29, 30] firstly performed the experimental investigation on the dispersion of pressure wave propagation in airwater bubbly flow and studied the propagation of pressure disturbance based on the twofluid model small perturbation analysis method. Through the comparison with experimental data, they found that the virtual mass force coefficient is a function of gas void. Chung et al. [31] calculated the sonic velocity versus angular frequency form the concept of bubble compressibility in a twocomponent bubbly flow regime. He also extended such a model to predict the sonic velocity of a vaporliquid system. Lee et al. [32] constructed the two fluids model to determine the pressure wave propagation speed for twophase bubbly flow, slug flow, and stratified flow by using pressure disturbance instead of virtual mass and other phase interfacial terms. The results fit well with the experimental in steam water and air water of Henry and theoretical analysis of Nguyen. Zhao and Li [33] derived the general formula of sonic velocity in gasliquid twophase flow linear analysis using the linear analysis of the closed fundamental equations of compressible gasliquid twophase flow. It is proposed that the appropriate formula for calculating sonic velocity in gasliquid twophase flows under usual conditions may be Wood adiabatic sonic velocity formula. By linearizing the conservation equations of twofluid model, Liu [34] derived a wave number equation of pressure wave for adiabatic gasliquid twophase flow. The effects of drag force and virtual mass force on propagation and dispersion of pressure wave were investigated. Xu and Chen [35] used the transient twofluid model to develop a general relation for acoustic waves with steamwater twophase mixture in onedimensional flowing system. Both the mechanical and thermal nonequilibrium are considered. Brennen [36] taken mass and heat exchanges into account and proposed more complete expressions of the speed of sound in twophase mixture. However, calibration of the mass and heat exchanges requires some further experimental investigations. Yeom and Chang [37] numerically investigated the wave propagation in the twophase flows. An assessment was made on the effect of interfacial friction terms. Zhang et al. [27] investigated the propagation of the pressure wave in the watergas twophase bubbly flow with a onedimensional twofluid model and employing small perturbation analysis. The governing equations are simplified and closured according to highspeed aerated flow characteristics in hydraulic engineering. The effects of aerated concentration, liquid pressure, perturbation frequency, and interfacial forces on the acoustic wave velocity and its attenuation in the aerated flow are also explored. With the application of thermal phase change model in computational fluid dynamics code CFX, Li et al. [38] proposed a pressure wave propagation model and investigated the pressure wave propagation characteristics in twophase fuel systems of liquidpropellant rocket. The propagation of pressure wave during the condensation of R404A and R134A refrigerants in pipe minichannels was given by Kuczyński [39]. Heat exchange between the phases in the condensation process was calculated by using the onedimensional form of Fourier’s equation.
In drilling industry, some scholars have been devoted to this aspect. In the late 1970s, the former Soviet AllUnion Drilling Technology Research Institute [40] began to study characteristics of pressure wave velocity in gasliquid twophase flow to detect early gas influx and achieved some important results. To study relationship between pressure wave velocity and gas void fraction, Li et al. [41] launched a gasdrilling mud twophase flow simulation experiment in vertical annulus. It is proved that the two phase flow of gasliquid patterns and the velocity of gas migration can be determined if the well depth, mud properties, and void fraction in bottom are given The method which is faster in detection time than the method of conventional kick detection was proposed. Starting with the analysis of transient flow, combining with the theory of transmission line, Wang [42] obtained the calculating model of frequency domain for the pulse velocity in drilling fluid built and the impendence and transmission operator of drilling fluid. Alternative initial responses to kicks for various well scenarios during MPD operations were also explored by Davoudi et al. [43]. Li et al. [44] proposed a mathematical model for predicting the attenuation and propagation velocity of measurement while drilling (MWD) pressure pulses in aerated drilling using the twophase flow model and considering the momentum and energy exchange at the phase interface, gravity of each phase, viscous pipe shear, and other closing conditions. According to the theory of unsteady flow, Xiushan [45] developed the formulas of transmission velocity for mud pulse signal. The formulas which cover all kinds of boundary conditions, including thin wall pipes and thick ones and interaction influence of gas content and solid content on transmission velocity are suitable for positive and negative mud pulse and accord with drilling practice. In a previous work, we proposed a united wave velocity model to predict the pressure wave velocity in gasdrilling mud twophase steady flow. The effect of well depth, back pressure, gas influx rate, virtual mass force, and angular frequency are all considered. However, under the effects of buoyancy and complicated turbulence interaction, the existing theoretical solutions are not involved in the dynamic model for predicting pressure wave velocity in fourphase fluid flowing along the drilling annulus when influx fluid migrate towards the top of wellbore.
In this paper, the drift model was used to analyze the flow characteristics of oil, gas, water, and drilling fluid multiphase flow. As the important characteristics of influx development, the relative motion of the interphase, such as slippage of gas phase and oil phase, is considered. Moreover, to predict pressure wave velocity in gasoilwaterdrilling mud fourphase flow in the annulus during MPD operations, a dynamic mathematical model is presented. By computing, the influence factors of pressure wave velocity, such as back pressure, gas void fraction, oil void fraction, influx time, influx rate, disturbance angular frequency, and virtual mass force, are analyzed.
2. Mathematical Model
Before introducing the new dynamic model and to make this point clear, this paper reviews the hydraulic system in MPD operations. The drilling system is a closed circulation with BP at wellhead. The key equipments include the rotating control device, dynamic well control system, conventional pressure control system, industrial personal computer, Coriolis meter, choke, and pressure sensor. First the drilling mud begins to circulate from mud tank, down the drill pipe, and the drill string and returns from the annulus travel back through mud pit where drilling solids are taken away and then to surface mud tank. An important function of the drilling fluid is to provide pressure support to the wellbore wall. The rock formation drilled through has some form of porosity filled with formation fluids. These fluids can be water, or in the case of a reservoir, hydrocarbons. The pressure in these fluids is referred to as the pore pressure. If the pore spaces are connected, these formations will also have permeability. Fluids can flow through them in response to a pressure gradient. The pressure in the annulus is controlled by varying BP to operate the fluid pressure in the wellbore. The aim in MPD is thus to maintain the pressure in the annulus between the two limits of pore pressure and fracture pressure [1].
2.1. The United Dynamic Model
When the bottom hole pressure is below the formation pressure, formation fluid will invade into the wellbore and the fourphase flow emerges in the annulus constituted by the drill string and wellbore. As seen in Figure 1, take any cross section of the wellbore as an infinitesimal control volume. In the infinitesimal control volume, the fourphase drilling fluid is consisted by drilling mud (considered as a pseudohomogeneous liquid), influx oil (considered as oil phase), influx natural gas (considered as gas phase), and influx water (considered as water phase). Appropriate assumptions and governing equations are critical to simulate realistic fourphase wellcontrol operations. The fourphase model was established based on the following assumptions:(1)it is unsteadystate fourphase flow;(2)the flow along the flow path is onedimensional;(3)the drilling mud is waterbased;(4)drilling mud is incompressible.
In the analysis of the multiphase flow characteristic, oil, water, and drilling mud which are all considered as liquid phase have great differences in physical and chemical properties. As waterbased medium, water and drilling mud have no substantive difference. In the flow processes, they blend quickly and have no clear phase boundary. Thus, the waterbased fluid phase is discussed as water phase, and the oil is considered as another phase.
As the presence of oil and gas interphase mass transfer in the multiphase flow system within the wellbore, the appropriate mass conservation equation can be listed according to oil, gas, and water three components. Given that is the mass fraction of components in phase, the mass conservation equations for fourphase mixture are
The momentum balance equation for the fourphase mixture is
Equation (2) is a general momentum balance equation including hydrostatic pressure gradient, frictional pressure loss gradient, and acceleration loss gradient.
Hence where
As propagation velocity is greatly affected by the gas void fraction and angular frequency of the pressure disturbance, the superficial velocity of flowing medium has almost no effect on the propagation velocity [46]; oil, water, and drilling can be considered as liquid phase for their similarity in mechanics. According to the twofluid model, the flow can be supposed to be gasliquid twophase flow from a macroscopic view.
To establish the wave velocity dispersion equation, the mass conservation equations for liquid and gas two phases can be written individually as follows: Hence, .
The gas momentum conservation equation is
The liquid momentum conservation equation is
The transfer of momentum and can be written by the following equations:
Virtual mass force is obtained by the equation in the following form: where and .
The momentum transfer term is described as [47]
The pressure difference between the liquid interface and liquid can be obtained by (11) where .
The gas interface pressure is defined as follows:
The pressure of the liquid is
The shear stress and the interphase shear stress can be described as
The Reynolds stress and interfacial average Reynolds stress are Hence, .
The wall shear stress of liquid phase is expressed as [48]
The pressure wave velocity of gas phase and that of liquid phase can be expressed in the following form:
Based on (17), the hydrodynamic equations of twofluid model (5)–(7) can be written in the following matrix form: Here, A is the matrix of parameters considered in relation to time; B is the matrix of parameters considered in relation to the spatial coordinate; C is the vector of extractions. By introducing the small disturbance theory, the disturbance of the state variable can be written as where is the wave number.
According to the solvable condition of the homogenous linear equations that the determinant of the equations is zero, the equation of pressure wave can be expressed in the following form:where , ,, and .
The real value of wave number is determined as the pressure wave velocity c, and pressure wave velocity in the fourphase flow is
2.2. Physical Models
To define the velocity of pressure wave propagation in the fourphase flow, the related physical models are required, such as equations of state, temperature distribution model, gas dissolution, and oil phase volume factor.
2.2.1. Equations of State for Gas
The equation of state (EOS) for gas can be expressed as
For MPa, the compression factor is obtained as follows: where , , .
For MPa, the compression factor under the condition of high pressure is [49] where Y is given by the follow equations:
2.2.2. Equations of State for Liquid
With °C, the density of drilling mud is expressed as follows: With °C, the density of drilling mud is expressed as follows:
2.2.3. Temperature Distribution Model
The temperature of the drilling mud at any depth of the wellbore is [50]
2.2.4. Gas Dissolution
Assuming gas goes into and comes out of solution instantaneously, gas solubility can be obtained by (29)
2.2.5. Oil Phase Volume Factor
The volume factor is calculated as follows:
2.3. Flow Pattern Prediction Models
The flow regime, the flow pattern, and structure of the flow are some of the important parameters to describe twophase gasquid flows, identify two phase gasquid flow regimes, and calculate the dynamic pressure wave propagation velocity. Thus, the transition among the three main flow regimes (bubbly, slug, and annular) is desirable to be known [51].
In the hydraulic calculation of annulus, the effective diameter of annulus [52] should be given as
The effective roughness of annulus can be calculated by
At low gas flow rate, the liquid is continuous phase, and the gas bubbles are dispersed in the liquid phase. Studies of Taitel [53] give the minimum diameter necessary to form bubbly flow as
The critical condition for forming bubbly flow is
For slug flow, the critical balance superficial flow rate of gas carrying droplets needs to meet the condition [54] that
For annular flow, the pattern transition criterions [55] is
2.3.1. Bubbly Flow
Gas void fraction of fourphase flow is described as
The value of the distribution factor can be determined by
Harmathy [56] established the calculation formula of gas slip velocity in bubbly flow based on the study of the migration velocity of the bubble in a stationary liquid as
The average density of fourphase mixture flow is
The oil void fraction for fourphase flow is
The value of the distribution factor is .
On the basis of the total liquid fluid, establish the oil phase velocity relationship as
According to crosssection flow rate, phase distribution, and slip mechanism of liquid phase, we can the draw the following relationship: where
Water void fraction is
Drilling mud void fraction is
Due to the similar physical properties of water and drilling fluids, , , and can be expressed by
The coefficient of virtual mass force for bubbly flow can be expressed as follows:
The coefficient of resistance coefficient C_{D} for bubbly flow can be expressed by
The friction pressure gradient for bubbly flow can be obtained from the following equation:
2.3.2. Slug Flow
Thevoid fraction of the four phases, , , , and , can be determined by (39), (43), (47), and (48) the same as bubbly flow.
The value of the distribution factor for slug flow can be described as
For slug flow, the slip velocity can be calculated as follows.
From experimental studies, Hasan and Kabir [57] established the calculation formula of drift velocity for slug flow on the basis of research on Taylor bubble migration rule of Davies and Taylor as
The coefficient of virtual mass force for slug flow can be expressed as follows:
The coefficient of resistance coefficient for slug flow can be expressed as
2.3.3. Annular Flow
As for annular flow, due to the miscible flow state of gas at center, the simplification can be .
The void fraction of gas can be determined by where is defined as
Oil void fraction is
Water void fraction is
The same as slug flow, and can be determined by (55) and (56).
3. Solution of the Dynamic Model
Now since obtaining the analytic solution of the aforementioned theoretical model directly is impossible, discretization of the model to a numerical model is required [58]. In this paper, the mathematical methods based on finite difference method provide a numerical solution approach for the dynamic model. As for the solution of the pressure wave velocity model, spatial domain includes the entire wellbore and the formation node; time domain is the time period, influx fluid flowing from the bottom hole to the wellhead along the wellbore. Discretizing the domain of determinacy, the entire spatial and time domain can be divided into discrete networked systems.
According to finite difference scheme, the four equations (5)–(7) are solved by using the finite difference method with computational cells shown in Figure 2; the difference equation systems of which described the basic principles of fourphase fluid motion in wellbore is presented as follows.
For the drilling mud phase,
For the water phase,
For the oil phase,
For the gas phase,
For momentum balance equation, where and
Besides finitedifference method, the characterized method and the NewtonRaphson iterative method are adopted to solve the united dynamic model in four phases along annulus. On the basis of discretization, the solution of models was realized by applying personally complied code on VB. NET and the solution procedure of that are shown in Figure 3.
The fourphase flow system is described completely by eight variables including pressure, temperature, gas void fraction and liquid void fraction, gas and liquid densities, and gas and liquid velocities. There are still five unknowns, such as gas and liquid velocities, gas fraction, pressure, and gas density based on the above assumptions. Therefore, five equations are required to compute the unknown variables with boundary conditions. At different annulus depth, we can obtain pressure, temperature, gas velocity, oil velocity, water velocity, drilling mud velocity and gas void fraction, oil void fraction, and water void fraction during MPD operations by application of the finitedifference. At initial time, the wellhead BP, wellbore structure, well depths, wellhead temperature, gasoilwaterdrilling mud properties, and so forth are known. On node , flow parameters, such as pressure, temperature, and void fraction of gasoilwaterdrilling mud multiphase, can be obtained by adopting finite difference. Then, the determinant equation (20) is calculated based on the calculated parameters. Omitting the two unreasonable roots, the pressure wave at different depth of annulus in MPD operations can be solved by (21). Upon substitution of the actual magnitudes, turned out to be the velocity of light. The process is repeated until the pressure wave velocity in the whole annulus has been obtained.
The published experimental data presented by Li et al. [41] is used to verify the new developed united dynamic model. In Figure 4, the points represent Li’s experimental data and the lines represent the calculation results. From the comparisons, it can be concluded that the developed united dynamic model agrees well with the experimental data. The averaged error between model and data in Figure 4 is 2.853%.
4. Analysis and Discussion
One of the potentially most useful capabilities of transient well control model is the ability to simulate various choke control procedures. Typically, the control of choke to maintain constant BHP depends on the experience and training of the hand on the choke. The gas and liquid flow rate measured by Coriolis meter and the BP measured by pressure senor is taken as the initial data for annulus pressure calculation. The well used for calculation is a gas well in Sichuan Region, China. The wellbore structure of this well is shown in Figure 5, and information about gasoilwaterdrilling mud properties, well design parameters, and operational conditions of calculation well are listed in Table 1.

In the following example, the drilling mud mixed with gas, oil, and water is taken as a fourphase flow medium when the influx of gas, oil, and water occurs at the bottom of the well. In the calculation, the BP is 0.2 MPa, and flow rate is m^{3}/h, m^{3}/h, m^{3}/h, and m^{3}/h, respectively. The propagation velocity of pressure wave in the fourphase flow is calculated and discussed by using the established model and well parameters.
4.1. Effect of Well Depth on the Pressure Wave Velocity
With bottom hole pressure being constant, the effect of well depth on pressure wave velocity is analyzed. When gasoilwater influx fluid invades into the wellbore at the bottom hole, gas, oil, water, and drilling mud fourphase flow appears in the annulus and migrates from the bottom of the well ( m) to the wellhead along the flow direction. It is assumed that the initial length of the gas pillar is 500 m under the bottom hole pressure and temperature condition in the analysis of the migration of influx gas. The migration and change of gas column including the solubility and length change lead to the variation of pressure wave velocity in the annulus. The formation and development of influx are a dynamic process, and the flow of fluid in the annulus is very complicated. The flow condition and flow parameter at any time are various.
Figure 6 illustrates the distribution of pressure wave velocity at different well depth after 5.16 minutes and 18.71 minutes’ development of influx. It is shown that the top of aerated mud column arrives at the position of m and m, respectively. Due to the circulation of drilling mud and slip of gas phase, the length of aerated mud column that represents the rising height of gas will increase with the gradual decrease of pressure along the annulus. Besides, in a segment of aerated mud column, the gas void fraction of the top is higher than that of the bottom for the slop of gas phase. So, in the declining stage, pressure wave velocity in aerated mud column presents an increasing tendency along the flow direction. As the gas migrates upwards, the casing pressure should be changed dynamically to keep the bottom hole pressure constant, which will lead to an increase of annular pressure and back pressure. As the influx develops for 18.71 minutes, the back pressure will be prompted to 2.39 MPa, increased by 0.29 MPa compared with that of 5.16 minutes.
4.2. Effect of BP on Pressure Wave Velocity
In Figures 7 and 8, the threedimensional figures illustrate the distributions of gas void fraction and variations of pressure wave velocity along the flow direction in the annulus at different time when BP is 0.1 MPa, 3.0 MPa, 6.0 MPa, and 9.0 MPa, respectively. When influx occurs, increase of applied BP is of the most commonly used approach for bottom hole pressure control. Control of back pressure has a great influence on the gas void fraction. In the two graphs, the two full lines give the shape of the void fraction and pressure wave velocity at different back pressure. In each graph, the two clusters of dotted lines at the twodimensional plane are the projection of the full lines, which, respectively, represent the influence of well depth and time on the distributions of gas void fraction or variations of pressure wave velocity at different back pressure. Increasing the casing pressure will make the compressibility decrease obviously and energy dissipation with the propagation of pressure wave decrease as well. Thus, it can be seen that, with the reduction of back pressure, void fraction shows an increasing trend and that the pressure wave velocity has a decrease at different depth of well. Meanwhile, at the high temperature and high pressure bottom hole, pressure wave velocity change in a comparatively smooth level in a consequence of compressibility of the fourphase fluid changing only in a small degree. As a result, as aerated drilling mud migrating from the bottom hole to the wellhead, gas void fraction significantly increases first slightly and then sharply along the flow direction a drop of annular pressure. Conversely, the pressure wave velocity shows a gradual fist and then remarkable decreasing tendency after the influx gas reaches the position close to the wellhead. Also, with the development of influx, the influence of BP on the gas void fraction and pressure wave velocity is larger over time.
4.3. Effect of Gas Influx Rate on Pressure Wave Velocity
Figures 9 and 10 graphically interpret the effect of gas influx rate on the distribution of gas void fraction and variations of pressure wave velocity along the flow direction in the annulus. The same as Figures 7 and 8, the two clusters of dotted lines at the twodimensional plane are the projection of the full lines. The pressure wave velocity at different well depth is the transient value with respect to time. It is noted that the changes of pressure wave velocity are obvious with the migration of gas, which is extremely significant at wellhead even at low gas influx rate. With the decreasing of pressure at the position near the wellhead, rapid expansion of gas volume appears, and as a result the gas void fraction increases sharply and the pressure wave velocity decrease obviously at the same time. However, the gas void fraction and pressure wave velocity have no obvious change at the bottom of well due to the high pressure at the bottom of well and insignificant compressibility of gas. When the mainly reflect fundamental factors, initial gas influx rate, is increased, pressure wave velocity shows an obvious descending trend with the increase of initial gas influx rate, the value of pressure wave velocity at the position close to the wellhead is much greater than that at bottom hole. At the gas void fraction range of about 5%~80%, as gas migrating along the annulus, the pressure wave velocity is decreased to a minimum at the wellhead with the gas void fraction increasing to a maximum.
4.4. Effect of Oil Influx Rate on Pressure Wave Velocity
Figure 11 shows the effect of oil influx on the gas void fraction in gasoilwaterdrilling mud fourphase flow. The increase of oil influx rate has little influence on the gas void fraction distribution, which is because gas void fraction is mainly impacted by pressure, temperature, and gas influx rate. When the oil influx rate increased from m^{3}/h to m^{3}/h, the change of mixture density is inconspicuous compared with drilling mud. Therefore, the gas void fraction keeps almost unchanged. Similarly, as shown in Figure 12, with the oil influx rate increase, the change of pressure wave velocity can be neglected due to the little influence of oil influx rate on the gas void fraction. Also, when the water influx invades the wellbore, the influx water will dissolve in drilling mud and reduce the density of drilling mud as the water phase and water base drilling mud have similar physical properties. But the influence of changing water influx rate on pressure wave velocity is too little to be neglected as the oil influx rate does.
4.5. Effect of Disturbance Angular Frequency on Pressure Wave Velocity
Figure 13 provides an illustration of the influence of disturbance angular frequency on the pressure wave velocity when gas, oil, and water influx generates. It is indicated that an increase in angular frequency results in a rise of pressure wave velocity. In addition, the influence on dynamic pressure wave velocity is enhanced with the continuation of time and the decrease of well depth as the gas migrates from the bottom hole to the wellhead. Accordingly, the disturbance angular frequency has a significantly distinct influence on the pressure wave velocity at the position near the wellhead due to the phenomenon that sharply increases of gas void fraction. For an increase in angular frequency from 1500 Hz to 10000 Hz, the influence on wave velocity is not as great as the increase form 50 Hz to 1500 Hz. This is consistent with the fact that the effect of angular frequency on the pressure wave appears at a low angular frequency.
5. Conclusions
Solution of dynamic pressure wave velocity in gasoilwaterdrilling mud fourphase flow is a lot more complicated than that in twophase flow. The dynamic model which takes consideration of gas solubility, oil compression coefficients, perturbation frequency, virtual mass force, and drag force is established. Owing to the solution of the migration of drilling mud which contains influx fluid (gas, water, and oil), the dynamic pressure wave varied with time and well depth is solved, and the main conclusions can be summarized as follows.
Agreeing well with the previous experimental data, the developed untied dynamic model can be used to calculate the wave velocity in the annulus. The model can accurately give the dynamic at different depth and time during the migration of influx, which will be beneficial to provide a reference for well control and MWD.
In the dynamic migration process, the pressure wave velocities at different well depth are changing constantly due to the circulation of drilling mud and slip of gas phase. The length of aerated mud column which represents the rising height of gas increases with the gradually decreasing pressure along the annulus. In the segment of aerated mud column, pressure wave velocity presents an increasing tendency along the flow direction.
The flow condition and flow parameter at any time are various. Both the static pressure and void fraction distribution along the annulus changed with the position in the annulus and time and have direct impact on the pressure wave velocity. The pressure wave velocity happens to decrease after an elapsed period of time in which influx fluid migrates along the annulus. The decreasing tendency will last until the aerated mud column pass through that point of annulus completely. Also, the time required is gradually shortened during process of influx migration and growth.
At the bottom hole, the void of gas, water, and oil is very low relative to drilling mud. With the migration of influx fluid, the variation of gas void fraction is much greater than other phases; in particular, a sharp increase of gas void fraction presents at the position near the wellhead for the rapid expansion of gas volume, while the void fractions of water, oil, and drilling mud are decreased for the nearly negligible variations of physical properties and drastic expansion of gas at the position near the wellhead.
When the BP at the wellhead increases, gas void fraction at any position of annulus decreases; instead, the pressure wave velocity increases. Moreover, the influence of BP on the gas void fraction is larger at position close to the wellhead than that at the bottom of well. With the development of influx, the influence of BP on the gas void fraction and pressure wave velocity is larger over time.
The pressure wave velocity is sensitive to the gas influx rate, in particular at the position close to the wellhead. Even at low range, the rise of gas influx rate also will result in a significantly downward trend of pressure wave velocity. On the contrary, the influence of changing water influx rate and oil influx rate on pressure wave velocity is too little to neglect.
Ignoring virtual mass force, the pressure wave velocity has an obvious dispersion characteristic at different position of well. The effect of angular frequency on the pressure wave appears at a low angular frequency, and the pressure wave velocity increases together with the growth of the angular frequency.
Nomenclature
Annulus effective cross area (m^{2})  
:  Oil phase volume factor 
:  Pressure wave velocity in gas (m/s) 
:  Pressure wave velocity in liquid phase (m/s) 
:  The coefficient of drag force 
:  Annulus effective diameter (m) 
:  Outer diameter of the inner pipe (m) 
:  Diameter of the wellbore (m) 
:  The minimum diameter to bubbly flow (m) 
:  Shear stresses coefficient of gas interface 
:  Shear stresses coefficient of liquid phase interface 
:  Coefficient of friction resistance 
:  Pressure applied on the twophase flow (N) 
:  Acceleration due to gravity (m^{2}/s) 
:  Conversion factor 
:  Geothermal temperature gradient (K/m) 
:  One discrete length of annulus (m) 
:  Annulus effective roughness (m) 
:  Roughness of outer pipe (m) 
:  Roughness of inner pipe (m) 
:  The depth of annulus (m) 
:  Air bubbles length of slug (m) 
:  Momentum transfer in gas interface (N/m^{3}) 
:  Momentum transfer in liquid interface (N/m^{3}) 
:  The nondrag force (virtual mass force) (N/m^{3}) 
:  The momentum transfer term (N/m^{3}) 
:  Parameter combining Thompson energy effects 
:  Gas interface pressure (MPa) 
:  Liquid interface pressure (MPa) 
:  Pressure (MPa) 
:  Initial pressure of annulus (MPa) 
:  Gas pressure (MPa) 
:  Critical pressure (kPa) 
:  Pressure of liquid phase (MPa) 
:  Reduced pressure 
:  Velocity of drilling mud (m^{3}/s) 
:  Velocity of gas phase (m^{3}/s) 
:  Velocity of liquid phase (m^{3}/s) 
:  Velocity of liquid phase within the annulus (m^{3}/h) 
:  Average diameter of the bubble (m) 
:  Standard oil density (kg/m^{3}) 
:  Standard gas density (kg/m^{3}) 
:  Constant of EOS (J/KgK) 
:  Air bubbles radius of slug (m) 
:  The bubble diameter (m) 
:  The distribution coefficient of gas velocity 
:  Temperature (K) 
:  Time (s) 
:  Critical temperature (K) 
:  Reduced temperature 
:  Undisturbed formation temperature at an depth (K) 
:  Undisturbed formation temperature of wellhead (K) 
:  Undisturbed temperature at the bottom hole (K) 
:  Gas and drilling mud flow velocity (m/s) 
:  Water mud flow velocity (m/s) 
:  Oil flow velocity (m/s) 
:  Gas flow velocity (m/s) 
:  The gas phase drift velocity 
:  Slip velocity (m/s) 
:  Critical flow rate of miscible phase (m/s) 
:  Superficial gas velocity (m/s) 
:  Liquid phase flow velocity (m/s) 
:  Superficial oil velocity (m/s) 
:  Superficial drilling mud velocity (m/s) 
:  Critical balance superficial flow rate of gas (m/s) 
:  Mass flow velocity (kg/m^{3}) 
:  Initial calculation height of annulus (m) 
:  Total well depth from surface (m). 
:  Gas density (kg/m^{3}) 
:  Reduced density 
:  Density under atmospheric pressure (kg/m^{3}) 
:  Oil density (kg/m^{3}) 
:  Drilling mud density (kg/m^{3}) 
:  Water density (kg/m^{3}) 
:  Liquid density (kg/m^{3}) 
:  Mixture density (kg/m^{3}) 
:  Shear stresses of noninterfacial mud (N/m^{2}) 
:  Shear stresses of liquid phase interface (N/m^{2}) 
:  Reynolds stress of liquid phase (N/m^{2}) 
Reynolds stress of liquid phase interface (N/m^{2})  
:  Reynolds stress of gas interface (N/m^{2}) 
:  Shear stresses of gas interface (N/m^{2}) 
:  Acceleration virtual force coefficient (m/s^{2}) 
:  Shear stresses of liquid phase along wall (N/m^{2}) 
:  Shear stresses of gas along wall (N/m^{2}) 
:  Reynolds stress of gas (N/m^{2}) 
:  Frictional pressure gradient (Pa/m) 
:  Frictional pressure coefficient 
:  Liquid phase density (kg/m^{3}) 
:  Surface tension (N/m^{2}) 
:  Surface tension of water phase and oil phase (N/m^{2}) 
:  Gas void fraction 
:  Drilling mud void fraction 
:  phase void fraction 
:  Oil void fraction 
:  Water void fraction 
:  Liquid void fraction 
:  Viscosity of liquid phase (Pa·s) 
:  Oil mass fraction 
:  Gas mass fraction 
:  Water mass fraction 
:  Drilling mass fraction 
:  The relative density of oil phase 
:  The relative density of gas phase. 
EOS:  Equations of state 
MPD:  Managed pressure drilling 
RCD:  The rotating control device 
DWC:  Dynamic well control system 
CPC:  Conventional pressure control system 
IPC:  Industrial personal computer. 
BP:  Back pressure (MPa) 
:  Pressure wave velocity (m/s) 
:  The coefficient of virtual mass force 
:  Depth of annulus (m) 
:  Gas influx rate at the bottom hole (m^{3}/h) 
Re:  The value of virtual mass force coefficient 
:  Angle frequency (Hz) 
:  Void fraction (%). 
Conflict of Interests
The authors declare that there is no conflict of interests regarding the publication of this paper.
Acknowledgments
This research work was cofinanced by the National Natural Science Foundation of China (no. 51274170) and Important National Science & Technology Specific Projects (2011ZX05022005005HZ). Without their support, this work would not have been possible.
References
 W. A. Bacon, Consideration of compressibility effects for appliedback pressure dynamic well control response to a gas kick in managed pressure drilling operations [M.S. thesis], University of Texas, Arlington, Tex, USA, 2011.
 S. Saeed and R. Lovorn, “Automated drilling systems for MPDthe reality,” in Proceedings of the IADC/SPE Drilling Conference and Exhibition, San Diego, Calif , USA, March 2012. View at: Google Scholar
 A. Malekpour and B. W. Karney, “Rapid filling analysis of pipelines with undulating profiles by the method of characteristics,” ISRN Applied Mathematics, vol. 2011, Article ID 930460, 16 pages, 2011. View at: Publisher Site  Google Scholar
 H. Santos, C. Leuchtenberg, and S. Shayegi, “Microflux control: the next generation in drilling process,” in Proceedings of the SPE Latin and Caribbean Petroleum Engineering Conference, Paper SPE 81183, PortofSpain, Trinidad, West Indies, 2003. View at: Google Scholar
 H. Guner, Simulation study of emerging well control methods for influxes caused by bottomhole pressure fluctuations during managed pressure drilling, Louisiana State University, Baton Rouge, La, USA, 2009.
 G. M. D. Oliveira, A. T. Franco, C. O. R. Negrão, A. L. Martins, and R. A. Silva, “Modeling and validation of pressure propagation in drilling fluids pumped into a closed well,” Journal of Petroleum Science and Engineering, vol. 103, pp. 61–71, 2013. View at: Publisher Site  Google Scholar
 I. Zubizarreta, “Pore pressure evolution, core damage and tripping out schedulesf: a computational fluid dynamics approach,” in Proceedings of the SPE/IADC Drilling Conference and Exhibition, Amsterdam, The Netherlands, 2013. View at: Google Scholar
 J. I. Hage and D. ter Avest, “Borehole acoustics applied to kick detection,” Journal of Petroleum Science and Engineering, vol. 12, no. 2, pp. 157–166, 1994. View at: Publisher Site  Google Scholar
 P. Thodi, M. McQueen, M. Paulin, G. Lanan, and INTECSEA, “Theory and Application of Probabilistic Method of Designing Customized Interval Control Valves Choke Trim for Multizone Intelligent Well Systems,” 2007. View at: Google Scholar
 A. Mallock, “The damping of sound by frothy liquids,” Proceedings of the Royal Society A: Mathematical, Physical and Engineering Sciences, vol. 84, no. 572, pp. 391–395, 1910. View at: Google Scholar
 A. B. Wood, A Textbook of Sound, G. Bill and Sons, London, UK, 1941.
 E. L. Carstensen and L. L. Foldy, “Propagation of sound through a liquid containing bubbles,” The Journal of the Acoustical Society of America, vol. 19, no. 3, pp. 481–501, 1941. View at: Google Scholar
 R. A. Thuraisingham, “Sound speed in bubbly water at megahertz frequencies,” Ultrasonics, vol. 36, no. 6, pp. 767–773, 1998. View at: Publisher Site  Google Scholar
 D. Hsieh and M. S. Plesset, “On the propagation of sound in a liquid containing gas bubbles,” Physics of Fluids, vol. 4, no. 8, pp. 970–975, 1961. View at: Publisher Site  Google Scholar  Zentralblatt MATH
 J. D. Murray, “Note on the propagation of disturbances in a liquid containing gas bubbles,” Applied Scientific Research, vol. 13, no. 1, pp. 281–290, 1964. View at: Publisher Site  Google Scholar
 G. B. Wallis, One Dimensional Two Phase Flow, vol. 7, McGrawHill, New York, NY, USA, 1969.
 F. J. Moody, “A pressure pulse model for twophase critical flow and sonic velocity,” Journal of Heat Transfer, vol. 91, no. 3, pp. 371–384, 1969. View at: Google Scholar
 D. F. D’Arcy, “On acoustic propagation and critical mass flux in twophase flow,” Journal of Heat Transfer, vol. 93, no. 4, pp. 413–421, 1971. View at: Publisher Site  Google Scholar
 D. McWilliam and R. K. Duggins, “Speed of sound in a bubbly liquids,” Proceedings of the Institution of Mechanical Engineers C, vol. 184, no. 3, pp. 102–107, 1969. View at: Google Scholar
 R. E Henry, M. A. Grolmes, and H. K. Fauske, “Pressurepulse propagation in twophase one and twocomponent mixtures,” Reactor Analysis and Safety Division, Argonne National Laboratory, 1971. View at: Google Scholar
 C. S. Martin and M. Padmanabhan, “Pressure pulse propagation in twocomponent slug flow,” Journal of Fluids Engineering, vol. 101, no. 1, pp. 44–52, 1979. View at: Publisher Site  Google Scholar
 Y. Mori, K. Hijikata, and A. Komine, “Propagation of pressure waves in twophase flow,” International Journal of Multiphase Flow, vol. 2, no. 2, pp. 139–152, 1975. View at: Publisher Site  Google Scholar
 Y. Mori, K. Huikata, and T. Ohmori, “Propagation of a pressure wave in twophase flow with very high void fraction,” International Journal of Multiphase Flow, vol. 2, no. 4, pp. 453–464, 1976. View at: Publisher Site  Google Scholar
 D. L. Nguyen, E. R. F. Winter, and M. Greiner, “Sonic velocity in twophase systems,” International Journal of Multiphase Flow, vol. 7, no. 3, pp. 311–320, 1981. View at: Publisher Site  Google Scholar
 R. C. Mecredy and L. J. Hamilton, “The effects of nonequilibrium heat, mass and momentum transfer on twophase sound speed,” International Journal of Heat and Mass Transfer, vol. 15, no. 1, pp. 61–72, 1972. View at: Publisher Site  Google Scholar  Zentralblatt MATH
 E. E. Michaelides and K. L. Zissis, “Velocity of sound in twophase mixtures,” International Journal of Heat and Fluid Flow, vol. 4, no. 2, pp. 79–84, 1983. View at: Publisher Site  Google Scholar
 H. W. Zhang, Z. P. Liu, D. Zhang, and Y. H. Wu, “Study on the sound velocity in an aerated flow,” Journal of Hydraulic Engineering, vol. 44, no. 9, pp. 1015–1022, 2013. View at: Google Scholar
 L. Y. Cheng, D. A. Drew, and R. T. Lahey Jr., “An analysis of wave propagation in bubbly twocomponent, twophase flow,” Journal of Heat Transfer, vol. 107, no. 2, pp. 402–408, 1985. View at: Publisher Site  Google Scholar
 A. E. Ruggles, R. T. Lahey Jr., D. A. Drew, and H. A. Scarton, “An investigation of the propagation of pressure perturbations in bubbly air/water flows,” Journal of Heat Transfer, vol. 110, no. 2, pp. 494–499, 1988. View at: Publisher Site  Google Scholar
 A. E. Ruggles, R. T. Lahey Jr., D. A. Drew, and H. A. Scarton, “The Relationship between standing waves, pressure pulse propagation, and critical flow rate in twophase mixtures,” Journal of Heat Transfer, vol. 111, no. 2, pp. 467–473, 1989. View at: Publisher Site  Google Scholar
 N. M. Chung, W. K. Lin, B. S. Pei, and Y. Y. Hsu, “A model for sound velocity in a twophase airwater bubbly flow,” Nuclear Technology, vol. 99, no. 1, pp. 80–89, 1992. View at: Google Scholar
 S. Lee, K. Chang, and K. Kim, “Pressure wave speeds from the characteristics of two fluids, twophase hyperbolic equation system,” International Journal of Multiphase Flow, vol. 24, no. 5, pp. 855–866, 1998. View at: Publisher Site  Google Scholar
 J. F. Zhao and W. Li, “Sonic velocity in gasliquid twophase flows,” Journal of Basic Science and Engineering, vol. 7, no. 3, pp. 321–325, 1999. View at: Google Scholar
 L. Liu, Y. Wang, and F. Zhou, “Propagation speed of pressure wave in gasliquid twophase flow,” Chinese Journal of Applied Mechanics, vol. 16, no. 3, pp. 22–27, 1999. View at: Google Scholar
 J. L. Xu and T. K. Chen, “Acoustic wave prediction in flowing steamwater twophase mixture,” International Journal of Heat and Mass Transfer, vol. 43, no. 7, pp. 1079–1088, 2000. View at: Publisher Site  Google Scholar
 C. E. Brennen, Fundamentals of Multiphase Flows, Cambridge University Press, 2005.
 G. Yeom and K. Chang, “The wave characteristics of twophase flows predicted by HLL scheme using interfacial friction terms,” Numerical Heat Transfer A: Applications, vol. 58, no. 5, pp. 356–384, 2010. View at: Publisher Site  Google Scholar
 Y. Li, C. Li, E. Chen, and Y. Ying, “Pressure wave propagation characteristics in a twophase flow pipeline for liquidpropellant rocket,” Aerospace Science and Technology, vol. 15, no. 6, pp. 453–464, 2011. View at: Publisher Site  Google Scholar
 W. Kuczyński, “Characterization of pressurewave propagation during the condensation of R404A and R134a refrigerants in pipe minichannels that undergo periodic hydrodynamic disturbances,” International Journal of Heat and Fluid Flow, vol. 40, pp. 135–150, 2013. View at: Publisher Site  Google Scholar
 Z. S. Fu, “The gas cut detect technology during drilling in Soviet Union,” Petroleum Drilling Techniques, vol. 18, no. 1, pp. 19–21, 1990. View at: Google Scholar
 X. F. Li, C. X. Guan, and X. X. Sui, “The theory of gas influx detection of pressure wave and its application,” Acta Petrolei Sinica, vol. 18, no. 3, pp. 128–133, 1997. View at: Google Scholar
 X. Wang, “Transmission velocity model in frequency domain of drillingfluid consecutive pulse signal,” Journal of Southwest Petroleum University, vol. 31, no. 3, pp. 138–141, 2009. View at: Publisher Site  Google Scholar
 M. Davoudi, J. R. Smith, B. Patel, and J. E. Chirinos, “Evaluation of alternative initial responses to kicks taken during managed pressure drilling,” in Proceedings of the IADC/SPE Drilling Conference and Exhibition, SPE128424MS, New Orleans, La, USA, February 2010. View at: Publisher Site  Google Scholar
 H. Li, G. Li, Y. Meng, G. Shu, K. Zhu, and X. Xu, “Attenuation law of MWD pulses in aerated drilling,” Petroleum Exploration and Development, vol. 39, no. 2, pp. 250–255, 2012. View at: Publisher Site  Google Scholar
 L. Xiushan, “Multiphase simulation technique of drilling fluid pulse transmission along well bore,” Acta Petrolei Sinica, vol. 27, no. 4, pp. 115–118, 2006. View at: Google Scholar
 F. Huang, M. Takahashi, and L. Guo, “Pressure wave propagation in airwater bubbly and slug flow,” Progress in Nuclear Energy, vol. 47, no. 1–4, pp. 648–655, 2005. View at: Publisher Site  Google Scholar
 J.W. Park, D. A. Drew, and R. T. Lahey Jr., “The analysis of void wave propagation in adiabatic monodispersed bubbly twophase flows using an ensembleaveraged twofluid model,” International Journal of Multiphase Flow, vol. 24, no. 7, pp. 1205–1244, 1998. View at: Publisher Site  Google Scholar
 G. B. Wallis, One Dimensional Two Phase Flow, vol. 7, McGraw Hill, New York, NY, USA, 1969.
 L. Yarborough and K. R. Hall, “How to solve equation of state for zfactors,” Oil and Gas Journal, vol. 72, no. 7, pp. 86–88, 1974. View at: Google Scholar
 A. R. Hasan and C. S. Kabir, “Wellbore heattransfer modeling and applications,” Journal of Petroleum Science and Engineering, vol. 86–87, pp. 127–136, 2012. View at: Publisher Site  Google Scholar
 M. AlOufi Fahd, An investigation of gas void fraction and transition conditions for twophase flow in an annular gap bubble column [Ph.D. thesis], Loughborough University, 2011.
 M. J. Sanchez, Comparison of correlations for predicting pressure losses in vertical multiphase annular flow [M.S. thesis], The University of Tulsa, 1972.
 Y. Taitel, D. Bornea, and A. E. Dukler, “Modelling flow pattern transitions for steady upward gasliquid flow in vertical tubes,” AIChE Journal, vol. 26, no. 3, pp. 345–354, 1980. View at: Publisher Site  Google Scholar
 G. Costigan and P. B. Whalley, “Slug flow regime identification from dynamic void fraction measurements in vertical airwater flows,” International Journal of Multiphase Flow, vol. 23, no. 2, pp. 263–282, 1997. View at: Publisher Site  Google Scholar  Zentralblatt MATH
 V. Casariego and A. T. Bourgoyne, “Generation, migration, and transportation of gascontaminated regions of drilling fluid,” in Proceedings of the SPE Annual Technical Conference and Exhibition, Houston, Tex, USA, October 1988. View at: Google Scholar
 T. Z. Harmathy, “Velocity of large drops and bubbles in media of restricted or infinite extent,” AIChE Journal, vol. 6, no. 2, pp. 281–288, 1960. View at: Google Scholar
 A. R. Hasan and C. S. Kabir, “Predicting multiphase flow behavior in a deviated well,” SPE Production Engineering, vol. 3, no. 4, pp. 474–482, 1988. View at: Publisher Site  Google Scholar
 X. W. Kong, Y. H. Lin, and Y. J. Qiu, “A new method for predicting the position of gas influx based on PRP in drilling operations,” Journal of Applied Mathematics, vol. 2014, Article ID 969465, 12 pages, 2014. View at: Publisher Site  Google Scholar
Copyright
Copyright © 2015 Xiangwei Kong 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.