UNIVERSITY OF OKLAHOMA GRADUATE COLLEGE MODELING OF MULTIPHASE FLOW IN WELLS UNDER NONISOTHERMAL AND NONEQUILIBRIUM CONDITIONS A THESIS SUBMITTED TO THE GRADUATE FACULTY In partial fulfillment of the requirements for the Degree of MASTER OF SCIENCE By GUILLERMO GERMAN MICHEL VILLAZÓN Norman, Oklahoma 2007
110
Embed
Gillermo Michel. Modeling of multiphase flow in wells
This document is posted to help you gain knowledge. Please leave a comment to let me know what you think about it! Share it to your friends and learn new things together.
Transcript
UNIVERSITY OF OKLAHOMA
GRADUATE COLLEGE
MODELING OF MULTIPHASE FLOW IN WELLS UNDER
NONISOTHERMAL AND NONEQUILIBRIUM CONDITIONS
A THESIS
SUBMITTED TO THE GRADUATE FACULTY
In partial fulfillment of the requirements for the
Degree of
MASTER OF SCIENCE
By
GUILLERMO GERMAN MICHEL VILLAZÓN Norman, Oklahoma
2007
MODELING OF MULTIPHASE FLOW IN WELLS UNDER NONISOTHERMAL AND NONEQUILIBRIUM CONDITIONS
A THESIS APPROVED FOR THE MEWBOURNE SCHOOL OF PETROLEUM AND GEOLOGICAL
“Science without religion is lame, religion without science is blind”
Albert Einstein (1879 - 1955)
"Science, Philosophy and Religion: a Symposium", 1941
To my parents, Ricardo and Stina, for their unconditional and unlimited support
and faith.
To my lovely wife, Alejandra, for her infinite love, kindness and care.
ACKNOWLEDGMENTS
I wish to acknowledge and thank many people for their cooperation during the
course of my studies at the University of Oklahoma.
In particular, I would like to express my most sincere gratitude to Dr. Faruk
Civan, chairman of my committee, for his advice and assistance in completing the
present work, for his patience and guidance, and for the trust he put in my work.
I acknowledge the time and dedication given by the members of my committee
Dr. Roy Knapp and Mr. Robert Hubbard.
I would like to thank to the ConocoPhillips Company for providing a fellowship
during my graduate studies.
I am grateful to our Creator, for all the blessings received in the path that he has
chosen for me.
iv
TABLE OF CONTENTS ACKNOWLEDGMENTS ..................................................................................... iv TABLE OF CONTENTS........................................................................................ v LIST OF FIGURES .............................................................................................. vii LIST OF TABLES................................................................................................ vii ABSTRACT......................................................................................................... viii 1. INTRODUCCION .......................................................................................... 1
1.1. OVERVIEW ................................................................................................ 1 1.2. DESCRIPTION OF THE PROBLEM......................................................... 2 1.3. PRESENT STUDY...................................................................................... 7 1.4. ORGANIZATION OF THE THESIS.......................................................... 9
2. LITERATURE REVIEW ............................................................................. 12 2.1. OVERVIEW .............................................................................................. 12 2.2. THE ANSARI ET AL. APPROACH ....................................................... 13 2.3. THE ASHEIM APPROACH ..................................................................... 14 2.4. THE AYALA AND ADEWUMI APPROACH ........................................ 15 2.5. THE DOWNAR-ZAPOLSKI ET AL. APPROACH................................. 16 2.6. THE BADUR AND BANASZKIEWICZ APPROACH........................... 17 2.7. THE FEBURIE ET AL. APPROACH....................................................... 18 2.8. THE CIVAN APPROACH........................................................................ 19 2.9. SUMMARY............................................................................................... 20
3. DETERMINATION AND CONSTITUTIVE EQUATIONS OF THE PHYSICAL PROPERTIES .................................................................................. 21
3.1. OVERVIEW .............................................................................................. 21 3.2. PHYSICAL PROPERTIES OF A MULTIPHASE FLUID ...................... 22 3.3. STANDARD CONSTITUTIVE EQUATIONS........................................ 26 3.4. PROPOSED MODEL FOR LIQUID HOLDUP ....................................... 28 3.5. RELAXATION TIME FOR PRODUCING WELLS................................ 32
4. DESCRIPTION OF THE MULTIPHASE FLOW OF A RESERVOIR FLUID IN WELLS ............................................................................................... 34
4.1. OVERVIEW .............................................................................................. 34 4.2. MODELING MULTIPHASE FLOW IN WELLS.................................... 35 4.3. STEADY-STATE MODEL UNDER NONISOTHERMAL AND NONEQUILIBRIUM CONDITIONS.............................................................. 42
5. SIMULATION OF THE MULTI-PHASE FLOW OF A RESERVOIR FLUID IN WELLS ............................................................................................... 45
5.1. OVERVIEW .............................................................................................. 45 5.2. SELECTION OF THE NUMERICAL DIFFERENTATION SCHEME.. 46 5.3. SCHEME OF THE NUMERICAL DIFFERENTATION......................... 47 5.4. COMPUTING THE CHANGE OF STATE.............................................. 48 5.5. COMPUTATIONAL PROCEDURE ........................................................ 52
6. VALIDATION AND APPLICATION......................................................... 56 6.1. OVERVIEW .............................................................................................. 56
v
6.2. DATA SELECTION.................................................................................. 57 6.3. SIMULATOR VALIDATION .................................................................. 58 6.4. UPWARD MOTION OF BLACK-OIL IN VERTICAL WELLS ............ 62 6.4.1 CASE 1: TWO-PHASE FLOW OF HEAVY CRUDE OIL ................... 63 6.4.2 CASE 2: TWO-PHASE FLOW OF LIGHT CRUDE OIL ..................... 68 6.4.3 CASE 3: THREE-PHASE FLOW OF A HEAVY CRUDE OIL............ 73
REFERENCES ..................................................................................................... 83 APPENDIX A: A DISCUSSION OF THE EQUILIBRIUM CONDITION FOR A MULTIPHASE FLUID SYSTEM........................................................................ 86 APPENDIX B: CORRELATIONS AND BASIC RELATIONSHIPS ................ 90 APPENDIX C: NOMENCLATURE.................................................................... 98
vi
LIST OF FIGURES
Figure 1-1 : Schematic of the motion in wells........................................................ 3 Figure 1-2 : Phases distribution for a cross-sectional area ..................................... 4 Figure 4-1 : Schematic for local properties in a conduit....................................... 35 Figure 4-2 : Example of local velocity, pressure and temperature distributions .. 36 Figure 5-1 : Segmentation of the production pipe in a well ................................. 48 Figure 5-2 : Flowchart for simulating multiphase flow in wells .......................... 54 Figure 6-1: Correlation for Case 1 ........................................................................ 60 Figure 6-2 : Correlation for Case 2 ....................................................................... 61 Figure 6-3 : Correlation for Case 3 ....................................................................... 61 Figure 6-4 : Pressure drop for Case 1 ................................................................... 64 Figure 6-5 : Void fraction for Case 1.................................................................... 65 Figure 6-6 : Temperature drop for Case 1 ............................................................ 66 Figure 6-7 : Temperature difference for Case 1.................................................... 66 Figure 6-8 : Dryness gradient for Case 1 .............................................................. 67 Figure 6-9 : Relaxation time for Case 1................................................................ 68 Figure 6-10 : Pressure drop for Case 2 ................................................................. 69 Figure 6-11 : Void fraction for Case 2.................................................................. 69 Figure 6-12 : Temperature drop for Case 2 .......................................................... 70 Figure 6-13 : Temperature difference for Case 2.................................................. 71 Figure 6-14 : Dryness gradient for Case 2 ............................................................ 71 Figure 6-15 : Relaxation time for Case 2.............................................................. 72 Figure 6-16 : Pressure drop for Case 3 ................................................................. 73 Figure 6-17 : Void fraction for Case 3.................................................................. 74 Figure 6-18 : Temperature drop for Case 3 .......................................................... 75 Figure 6-19 : Temperature difference for Case 3.................................................. 75 Figure 6-20 : Dryness gradient for Case 3 ............................................................ 76 Figure 6-21 : Relaxation time for Case 3.............................................................. 77
LIST OF TABLES Table 2-1 : Literature Review............................................................................... 20 Table 6-1 : Data considered for application.......................................................... 57 Table 6-2 : Adjustable parameters and Coefficient of Determination.................. 62
vii
ABSTRACT
The multiphase flow of reservoir fluids in producing wells has been a subject of
investigation in various previous studies. In general, the motion of reservoir fluids
undergoing a gas separation along the well has been modeled by using empirical
correlations. Recently, however, the emphasis has shifted to theoretical modeling.
The present study provides a rigorous theoretical approach for modeling of the
upward motion of reservoir fluids considering the gas separation phenomenon in
the production wells.
The reservoir fluid is represented as a mixture of three phases, consisting of the
gas, oil, and water phases. A homogenous fluid model is formulated for general
purposes for describing the upward motion of a multiphase fluid system in pipes.
But, its application is demonstrated for well operations under the steady-state
conditions. The upward motion is considered under the non-isothermal and non-
equilibrium conditions by taking into account the irreversible loss in energy. The
loss in energy is mainly due to the interaction of the system with the
surroundings. The homogeneous model is simplified for the steady-state motion
in pipes having constant and circular cross-sectional areas.
The separation of the gas phase is considered to cause a non-equilibrium effect in
the upward motion. The non-equilibrium effect occurs when the phase velocities
are not equal. Two approaches are presented for describing the non-equilibrium
viii
effect on the bases of the prediction of the liquid holdup and the estimation of the
relaxation time occurring in the gas phase separation.
A new improved model for prediction of the liquid holdup is formulated. The
liquid holdup is predicted by the means of a constitutive equation. The
constitutive equation is based on the mixture density and the slip ratio. The
proposed holdup model provides a closure for the developed homogenous model
and it is employed for the application in the present study.
A practical means for solving the resulting differential equations is developed. A
series of simulated case studies are performed using the selected data. The data
was acquired from producing vertical wells and published in a previous study.
After validating the output data of the simulations, the motion of the studied cases
is described and characterized. The characterization includes the behavior of the
relaxation time occurring in the gas phase separation. The model developed here
provides important improvements over the existing models, which do not take
into account accurately the effects of the relaxation phenomenon and the liquid
holdup.
ix
1. INTRODUCCION
CHAPTER 1
INTRODUCTION
1.1. OVERVIEW
The particular phenomenon of concern of this thesis is the upward motion of
reservoir fluids in producing wells. In this chapter, the motivation and the scope
of the present study are established. A description of the fluid flow in petroleum
wells in terms of the governing physical phenomena is addressed. Then, the
specific objectives of the present study are defined. The specific objectives are
considered to accomplish the solution of the main problem. At the end, the
organization of the study towards the fulfillment of objectives is presented.
1
1.2. DESCRIPTION OF THE PROBLEM
In general, hydrocarbon fluids present in reservoirs contain a large number of
various substances. Each of these substances has different physical properties and
behavior affecting in specific ways the properties of the fluid phases. Moreover,
the interfaces or surface borders between the fluid phases have physical properties
and behavior on their own. Consequently, large amounts of measurements have to
be done in order to determine the required properties by means of a detailed
model. For that reason, theoretical models of fluid dynamics for reservoir fluids in
producing wells have been proposed in various types and successes.
Typically, the reservoir fluid consists of three distinct phases1,22. These are the
gas, oil, and water phases. Thus, the flow of the reservoir fluid in wells can be
modeled as the flow of a multiphase-fluid system of several phases.
For a producing well, the motion of the reservoir fluid is depicted in figure 1-1.
By considering the well fluid as a single multiphase-fluid system containing gas,
oil, and water phases, the flow in the production pipe can be described by the
fundamental equations governing the flow of fluids in conduits.
2
Oil WaterGas
Figure 1-1 : Schematic of the motion in wells
In the present modeling approach, it is assumed that the three fluid phases (gas,
oil, and water) are homogeneous and uniformly distributed over a cross-sectional
area (figure 1-2a). As the multi-phase fluid flows upward along the pipe from the
well-bore to the wellhead, an interface mass transfer is considered to occur across
the gas and liquid (oil and water) interphases14. The mass transfer may be
bidirectional. However, only the separation of the gaseous phase (gas) from the
liquid phases (oil and water) is considered in this study. Because the pressure
continuously decreases in the upward motion of the fluid, there is no dissolution
of the gas phase into the liquid phases occurring during flow. Within a particular
cross-sectional area, the multiphase fluid has a distribution of the mass fraction
3
for the various phases set by the local state of properties. While moving upward,
the multiphase fluid of various phases undergo a change in mass fraction
distribution along the well (figure 1-2b and figure 1-2c).
Multiphase Water Oil Gas
(a) (b) (c)
Figure 1-2 : Phases distribution for a cross-sectional area
Usually, depending on the prevailing conditions in a pipe, the interface mass
transfer between the liquid and gaseous phases occurs without reaching an
equilibrium state when the flow is sufficiently fast. Hence, it is reasonable to
consider that the mass transfer between the various phases occurs at a non-
equilibrium state17 (flashing) process. This means that the mass transfer occurs
dynamically backward and forward between the various phases. Unfortunately,
there is no well-proven and satisfactory model available for such cases involving
the flashing hydrocarbons.
A generalized model for flashing fluids has been developed in a limited number
of previous studies6,17. This flashing model considers a relaxation in time for gas
4
separation from the liquid phases due to the slow mass transfer between the gas
and liquid phases. Consequently, a unidirectional and cumulative mass transfer
from the liquid phases to the gaseous phase is assumed for the present study.
The mass transfer from the liquid phases to the gaseous phase begins when the
multiphase-fluid system pressure drops to below the bubble-point pressure. As the
multiphase-fluid flows along the pipe length, the pressure and temperature of the
fluid system decrease. The motion of the multiphase-fluid causes a pressure drop.
Simultaneously, the heat transfer by conduction and convection, the effect of the
fluid expansion and the effect of friction cause a temperature change. The
temperature change by expansion is referred to as the Joule-Thompson effect. The
heat transfer can be computed knowing the temperature of the surroundings. The
surrounding temperature is set mainly by the insulation technique of the conduit
and the geothermal gradient of the surrounding rock formation.
Another approach to modeling the motion of the multiphase fluid system is to
estimate the volumetric fraction of the liquid phases, referred to as liquid holdup.
Several studies have been performed for predicting the liquid holdup in wells.
These studies model the deviation from equilibrium in terms of a slippage
occurring between the gas phase and the liquid phases (oil and water) rather than
as a flashing process.
5
However, the slippage phenomenon have been proven to be complex enough to
be modeled by a single correlation for the liquid holdup9,18. All developed models
use a set of these correlations for predicting the liquid holdup. Usually, different
correlations are employed depending on the local conditions along the producing
pipe.
In the field facilities, the hydrocarbon fluid can be separated into three
components1,22 (gas, oil, and water). They are called the pseudo-components.
These components are at atmospheric conditions and behave differently than the
phases flowing through the conduit. Therefore, the gas, oil, and water pseudo-
components are different substances than the gas, oil, and water phases.
Because both the pressure and temperature are changing along the pipe, it is
impractical to measure directly all the physical properties of the reservoir fluid
phases during flow. Hence, several correlations have been developed for
estimating the properties of these phases. In general, the properties of the pseudo-
components and the conditions of the local state are required for these
correlations. Thus, by knowing the afore-mentioned properties, the physical
properties of the various phases can be estimated.
6
1.3. PRESENT STUDY
The scope of the present study is to develop an improved model for the flow of a
reservoir fluid as a multiphase fluid system in wells producing under steady-state
conditions. The flow is assumed to be non-adiabatic and non-isothermal
considering the convective and conductive heat transfer as the energy losses and
the effect of the friction. The reservoir fluid is flowing along the production pipe
with a non-equilibrium mass transfer across the interface between the liquid
phases and the gas phase. The study cases are considered based on the published
data for producing vertical wells.
Then, the fundamental laws of mass, momentum and energy conservation are
applied to describe the change in velocity, pressure, and temperature of reservoir
fluids flowing through the wells. However, the change in density of the fluid
cannot be obtained by predicting the previously mentioned changes alone because
the gas mass transfer from the liquid phases to the gas phase is not at equilibrium
during flow. Therefore, the flashing process occurring inside the production pipe
has to be modeled by other means.
The velocities of the phases are equal when the system has reached an
equilibrium9,18 as shown in Appendix A. For this reason, the deviation from
equilibrium is predicted by estimating the phase velocities. In this study, a new
7
equation for obtaining the ratio of the gas phase velocity to the liquid phase
velocity is presented. With this velocity ratio, the liquid holdup can be obtained
accurately as well as the flowing density.
By using the fundamental laws of conservation and the proposed method for
liquid holdup prediction, a series of simulations are then performed to accurately
predict the drop in pressure and temperature along the wells of each study case.
Both the relaxation time and the liquid holdup models describe the same
phenomenon satisfactorily which is the deviation from the equilibrium. Thus, the
behavior of the relaxation along the pipe is estimated with the data yielded by the
simulations.
The fundamental laws of conservation are formulated in their differential forms.
Thus, all the properties of the multiphase fluid system are either spatially
averaged in nature or homogenous. A numerical method is developed for solving
the differential equations given by the conservation laws. This numerical method
is extensively described.
The main objective of the present study is to model and characterize the flow of a
reservoir fluid in producing wells. The main objective is accomplished by the
following specific objectives:
8
• Develop a technique for estimating the properties for a multiphase fluid
system.
• Introduce a new improved model for estimating the liquid holdup as a
deviation from equilibrium.
• Develop a homogenous model applicable to the flow of reservoir fluids in
wells under non-isothermal and non-adiabatic conditions.
• Prove the relaxation time as a property that characterizes the deviation
from equilibrium for flowing reservoir fluids.
• Solve the developed homogenous model for simulating the flow with a
numerical scheme.
• Validate the results of the simulations by using a correlation developed for
experimental measurement of the void fraction.
1.4. ORGANIZATION OF THE THESIS
The contents of this thesis are organized and reported in seven chapters and two
appendixes as described in the following.
The current chapter, Chapter One, provides an overview of the problem of interest
and presents the scope of the present study. Chapter Two presents a
comprehensive review of the relevant literature.
9
Chapter Three provides a technique to approximate the properties of a multiphase
fluid system under non-equilibrium conditions. The reservoir fluid is represented
as a multiphase fluid system.
Chapter Four describes the modeling of a reservoir fluid in motion in wells. The
laws of mass, momentum, and energy conservation are expressed in differential
forms. A homogenous model for pipes with circular cross-sectional area is
developed. The cross-sectional area can be either constant or variable.
Chapter Five presents the numerical method developed in order to perform the
simulations using the technique specified in Chapter Three and the homogenous
model developed in Chapter Four.
Chapter Six shows the relevant results obtained by the simulations. The results are
validated with a model developed in the literature for correlating experimental
data. The application is illustrated by means of three study cases for the upward
motion of the gas/oil/water mixtures in wells.
Chapter Seven contains the discussion and conclusions after analyzing the results
obtained for the application.
10
Appendix A shows that a multiphase fluid system is at equilibrium condition if
the phase velocities are equal.
Appendix B presents a collection of correlations required for estimating the
properties of the gas, oil, and water phases as well as the wall surface properties
of the pipe.
Appendix C illustrates the adopted nomenclature for the various properties,
parameters and variables employed by the formulations in the present study.
11
2. LITERATURE REVIEW
CHAPTER 2
LITERATURE REVIEW
2.1. OVERVIEW
In this chapter, a review of the relevant studies about the flow and behavior of the
flashing fluids at steady-state is presented. These studies describe the flashing
phenomenon with different approaches. The description of each approach is
properly addressed towards detailing the features of interest for the present study.
Usually, some simplifications were made in order to enable the measurement of
the pertinent properties. At the end, a table summarizes the key features covered
by the current and the reviewed studies involved in modeling the flashing fluids.
12
2.2. THE ANSARI ET AL. APPROACH
The mechanistic approach proposed by Ansari et al. 2 modeled the upward flow of
reservoir fluids in pipes. The model compiled and systematized the use of several
correlations for predicting the liquid holdup and flow pattern distribution along
the well. Separate models and correlations were proposed for each flow pattern.
Consequently, the estimation of the flowing density is not continuous when a
change in flow pattern is predicted.
The chosen correlations were selected to minimize the error as it was
demonstrated in the error analysis section of the previous studies. The validation
of the model was executed with data measured in producing vertical wells
although the formulations can be applied for all angles of inclination.
In this study, the multiphase fluid system is defined as a mixture of phases
flowing within a pipe having a constant and circular cross-sectional area.
However, this approach is not a homogeneous model because the velocity of the
mixture is set equal to the volumetric flux even though the system is not at
equilibrium.
The pressure drop is estimated mainly by the prediction of the liquid holdup.
However, there is no specification on how to incorporate the effect of a
13
simultaneous drop in pressure and temperature. Thus, it is implied that the system
is isothermal having taken the average between the inlet and outlet temperatures.
2.3. THE ASHEIM APPROACH
The mathematical approach proposed by Asheim3 modeled the slippage occurring
in an upward motion of reservoir fluids in producing wells with constant diameter
and variable inclination. The deviation from the equilibrium is modeled by the
liquid holdup prediction.
The prediction of the liquid hold up is achieved by assuming a linear relationship
between the velocity of the gas phase and the velocity of the liquid phase. The
linear parameters have to be assumed a priori for the phase velocity relationship.
This assumption resulted in a quadratic relationship between the velocity of the
various phases and the liquid holdup. Therefore, there is no assurance for a
continuous estimation of the flowing density when the mixture is undergoing a
change from the saturated to the unsaturated fluid conditions. In the error
analysis, it was proven that this approach minimizes the error in history matching.
Although the multiphase fluid system is defined as a mixture of phases, this
approach is not a homogeneous model because the velocity of the mixture is set
equal to the volumetric flux.
14
In this study, the liquid holdup prediction is mainly set by the pressure drop. The
flow is assumed to be isothermal having taken the average between the inlet and
outlet temperatures.
2.4. THE AYALA AND ADEWUMI APPROACH
The multi-fluid approach proposed by Ayala and Adewumi4 modeled the flow of
gas and condensates along a transmission pipeline with constant diameter. The
multiphase fluid system is defined as a collection of two completely separated
phases undergoing mass transfer across the interface.
The pressure drop is mainly set by the liquid hold up and the mechanical loss of
momentum occurring at the interface. The mechanical loss of momentum is the
free term in the modeling that provides closure in the pressure, temperature, gas
velocity, and liquid velocity formulations. The gas density and liquid density are
obtained by the equation of state for each phase. Furthermore, the mass transfer is
estimated by a numerical scheme based on the gas density equation of state.
The mechanical loss of momentum is estimated by the correlations describing
several flow pattern distributions. For that reason, there is no continuous
15
estimation of the liquid phase velocity when the mixture is undergoing a change
from a saturated to an unsaturated fluid.
2.5. THE DOWNAR-ZAPOLSKI ET AL. APPROACH
The homogenous approach proposed by Downar-Zapolski et al.17 modeled the
flow of water and steam along a horizontal conduit with variable cross-sectional
area. However, the effect caused by the change in the cross-sectional area was
only considered in the velocity formulation. The cross-sectional area effect was
omitted in the pressure and temperature formulations.
The multiphase fluid system is defined as a mixture of water and steam phases
flowing at the non-equilibrium and adiabatic conditions. The critical flow of the
mixture in pipes with small diameter is the key testing condition.
The deviation from the equilibrium is described by the means of a relaxation time
occurring in the steam separation. The pressure drop and the void fraction are
known a priori by experimental measurement. Then, a correlation is developed to
estimate the relaxation time by using the experimental data at various flow rates.
16
2.6. THE BADUR AND BANASZKIEWICZ APPROACH
The homogenous approach proposed by Badur and Banaszkiewicz5 modeled the
flow of water and steam along a conduit with small cross-sectional area. The main
feature is testing and describing the mass transfer of the gas phase by the means
of a constitutive equation. The multiphase fluid system is defined as a mixture of
the water and steam phases where the flowing conditions cause a deviation from
the equilibrium.
The pipe is horizontal with a variable cross-sectional area. The effect of a variable
area is omitted in the pressure and temperature formulations but it is considered in
the velocity formulation.
The homogenous model is closed by a constitutive equation for the flowing fluid
quality. This constitutive equation includes the relaxation time as a coefficient by
assuming that the flow is adiabatic. The remaining constant parameters of this
equation can be correlated by analyzing the experimental data.
An adequate correlation was developed for two different flow rates. The predicted
and experimental pressure drops are compared by the means of a plot as well as
the predicted and experimental void fractions.
17
2.7. THE FEBURIE ET AL. APPROACH
The homogenous approach proposed by Feburie et al.19 modeled the flow of
steam and water derived from a multi-fluid model. The model was applied to the
flow of steam/water mixtures along horizontal conduits with variable and small
cross-sectional area.
The multiphase fluid system is defined as a mixture of superheated water,
saturated water and saturated steam. Although the phases are assumed to flow at
equilibrium conditions, the deviation from the equilibrium was addressed by
partitioning the water phase into the superheated water phase and the saturated
water phase.
The homogenous model is closed by a constitutive equation for the relaxation in
the mass transfer occurring at the interface between the superheated water and the
saturated steam/water mixture. The temperature change is formulated by the
change in entropy considering irreversible heat transfer towards the surroundings.
However, the effect caused by the variable cross-sectional area is omitted in the
pressure and entropy formulations.
The validity of the constitutive equation was tested by comparing the predicted
pressure drop with the experimental pressure drop at various flowing conditions.
18
2.8. THE CIVAN APPROACH
The mechanistic approach proposed by Civan14 modeled the upward flow of
reservoir fluids in wells at non-equilibrium conditions. The flow is assumed to be
isothermal. It was implied that the constant temperature considered in this model
is the average between the inlet and outlet temperatures.
The key feature of the study is to demonstrate that the deviation from the
equilibrium in producing wells can be modeled by means of the relaxation time
concept even though this property was originally developed for tubes with small
diameter. It was shown that the law of conservation for the gas phase can be used
to give closure in a homogenous model for producing wells. Nonetheless, this
approach is not a homogeneous model because the velocity of the mixture is set
equal to volumetric flux.
The multiphase fluid system is defined as a mixture of gas, oil and water flowing
within a pipe having a constant and circular cross-sectional area. Although the
model was formulated for all angles of inclination, the application only
considered a vertical well.
The relaxation time is estimated by a correlation developed for the steam/water
mixture flowing in small tubes. The mass transfer of the gas phase and the density
19
of the mixture are set mainly by the relaxation time. Consequently, the pressure
drop and the quality gradient are set by this property as well.
2.9. SUMMARY
The table 2-1 summarizes the main attributes of the present study and all the
The motion of a reservoir fluid along the production pipe from the well-bore to
the wellhead was modeled by means of a homogenous model with liquid holdup
and applied to three study cases. The study cases considered the flow of gas/oil
and gas/water/oil mixtures in vertical wells. Consequently, the reservoir fluid was
considered as a multiphase fluid system. For this reason, all the properties of this
fluid were predicted by estimating the properties of its phases.
78
A homogenous model applicable for reservoir fluids flowing along a pipe with
constant cross-sectional area was developed. This model was simplified for
flowing fluids across a circular pipe at steady state. The differential conservation
laws of mass, momentum, and energy were adopted.
By estimating the multiphase fluid properties, the differential laws of
conservation were solved to compute the change in pressure, the change in
temperature and the flowing velocity. The change in pressure was computed
considering the gravitational force, the friction loss and the fluid compressibility.
The change in temperature was computed considering the change in pressure, the
friction effect, and the energy dissipation towards the surroundings. The velocity
is computed by knowing the fluid density.
Because the flow might not be at equilibrium, the density of the multiphase fluid
system was not obtained from the conservation laws. Two approaches were
presented to determine the density. One is the differential conservation law of
gaseous mass that can be solved by estimating the relaxation time, and the other is
predicting the liquid holdup by estimating the slip ratio. The later was used for the
application.
A new approach for predicting the liquid hold up was introduced. The liquid
holdup was associated to the deviation from equilibrium. It was related to the slip
79
ratio which is the ratio of the phase velocities. The slip ratio was computed by
interpolation. The interpolation is based on a Lagrange’s polynomial of the
second order relating the slip ratio with the equilibrium density. The slip ratio for
saturated oil and saturated gas is set to be the unit which represents the
equilibrium. An apparent slip ratio is assumed at the surface. The equilibrium
density is computed by knowing the system pressure and temperature.
The simulation of the reservoir fluid motion was executed by applying both the
homogenous model and liquid holdup model in conjunction. The homogenous
model was solved by using the classical fourth-order of Runge-Kutta method.
Because the slip ratio and the external temperature are not known at the wellhead,
the pressure and the temperature of the multiphase fluid at the well-bore were
used to give a closure to the system. The shooting method was applied for this
purpose. Thus, the slip ratio and the external temperature at the wellhead were
guessed until the pressure and temperature of the multiphase fluid at the well-bore
are predicted by the simulation. The simulation results were validated by a
generalized model for void fraction prediction.
Having validated the simulation results, the relaxation time was computed. The
relaxation time is the main property that delineates the conservation law for the
gas phase. This law completes the homogenous model in order to compute the
flowing density.
80
The assumption made of a constant geothermal gradient and a quadratic
relationship between the slip ratio and the equilibrium density did not introduce a
significant error for modeling the present phenomenon. The deviation of the
results from the Butterworth’s model8 was negligible for all samples.
The pressure drop tended to be linear even for the samples with considerable
deviation from equilibrium.
The temperature drop tended to be non-linear for all samples. The motion of the
fluid proved to be quick enough to delay the heat dissipation towards the
surroundings. This became evident when showing the difference between the
fluid temperature and the apparent external temperature.
It was proven that the relaxation time characterize the deviation from the
equilibrium for flowing fluids in wells. It was suggested by Downar-Zapolski et
al.17 that the relaxation time is a fluid property. In their study, they developed a
single correlation for the relaxation time of flashing water. However, the
relaxation time presented a unique curve for each sample of the study cases. Thus,
the present formulation of this property stills depends on the prevailing conditions
of the flowing reservoir fluid. Nevertheless, the relaxation time curves proved to
be a family of curves. This suggests that the present formulation of the relaxation
time can be adjusted towards becoming a fluid property for the reservoir fluids.
81
7.2. CONCLUSIONS
Having analyzed the results presented in the applications, the following
conclusions concerned with the present study cases have been reached:
1. Assuming a constant geothermal gradient does not introduce a significant
error.
2. The proposed approach predicts a continuously varying liquid-holdup by
means of interpolating the slip ratio.
3. The heat dissipation to the surroundings and the fluid expansion and the
friction effect, cause a non-linear temperature drop.
4. The upward motion of reservoir fluids in producing wells can be
successfully modeled by the developed homogenous model in conjunction
with the proposed model for liquid holdup prediction.
5. The relaxation time of gas separation proved to be an adequate property
for characterizing the deviation from the equilibrium for reservoir fluids.
6. The conservation law for the gas phase and the relaxation time of gas
separation from the liquid phases can be applied in order to achieve a
closure in the area-averaged homogenous model.
82
REFERENCES
[1] Amyx, J., Bass, D., Whiting, R., “Petroleum Reservoir Engineering: physical properties”, McGraw-Hill, USA, pp. 211-472, 1960
[2] Ansari, A.M., Sylvester, N.D., Sarica, C., Shoham, O., Brill, J.P., “A Comprehensive Mechanistic Model for Upward Two-Phase Flow in Wellbores”, SPE Production & Facilities, pp. 143-152, May, 1994
[3] Asheim, H., “MONA, An Accurate Two-Phase Well Flow Model Based on Phase Slippage”, SPE Production Engineering, pp 221-230, May 1986
[4] Ayala, L. F., Adewumi, M. A., “Low-Liquid Loading Multiphase Flow in Natural Gas Pipelines”, J. of the Energy Resources Technology, Vol. 125, pp. 284-293, 2003.
[5] Badur, J., Banaszkiewicz, M., “A Model of two-phase flow with relaxation-gradient microstructure”, Third International Conference on Multiphase Flow, held in Lyon, France, June 8-12, 1998.
[6] Bilicki, Z., Kestin, J., “Physical Aspects of the Relaxation Model in Two-Phase Flow”, Proceedings of the Royal Society of London, Series A, Mathematical and Physical Sciences, Vol.428, No. 1875, pp 379-397, Apr. 9, 1990.
[7] Bird, R., Stewart, W., Lightfoot, E., “Transport phenomena”, Wiley, USA, pp. 71-110, pp. 310-342, 1960
[8] Butterworth, D., “A Comparison of Some Void-Fraction Relationships for Cocurrent Gas-Liquid Flow”, International Journal of Multiphase Flow, Vol. 1, pp 845-850, 1975
[9] Brill, J. P., Mukherjee, H., “Multiphase Flow in Wells”, SPE, Richardson, p.16, pp. 102-122, 1999.
83
[10] Cazaraez-Candia, O., Vásquez-Cruz, M., “Prediction of Pressure, Temperature and Velocity Distribution of Two-Phase Flow in Oil Wells”, Journal of Petroleum Science and Engineering, Vol. 46, pp. 195-208, 2005.
[11] Cengel, Y., Boles, M., “Thermodynamics: an engineering approach”, McGraw-Hill, USA, pp. 150-155, pp. 603-626, 2002.
[12] Chapra, S., Canale, R., “Numerical Methods for Engineers”, McGrw-Hill, USA, pp. 675-718, 1998
[13] Chierici, G.L., Ciucci, G.M., Sclocchi, G., “Two-Phase Vertical Flow in Oil Wells – Prediction of Pressure Drop”, SPE Journal of Petroleum Technology, pp. 927-938, August 1974.
[14] Civan, F., “Including Non-equilibrium Effects in Models for Rapid Multiphase Flow in Wells”, SPE Paper 90583, the 2004 SPE Annual Technical Conference and Exhibition, held in Houston, Texas, 26-29 September 2004.
[15] Civan, F., “Including Non-equilibrium Relaxation in Models for Rapid Multiphase Flow in Wells”, SPE Production&Operations Journal, pp. 98-106, February 2006.
[16] Crowe, C., Elger, D., Roberson, J., “Engineering fluid mechanics”, Wiley, USA, pp. 368-434,2005
[17] Downar-Zapolski, P., Bilicki, Z., Bolle, L. and Franco, J., “The Non-equilibrium Relaxation Model for One-Dimensional Flashings Liquid Flow”, International J. Multiphase Flow, Vol. 22, No. 3, pp. 473-483, 1996.
[18] Faghri, A., Zhang, Y., “Transport phenomena in multiphase systems”, Elsevier Academic Presss, pp. 238-320, pp. 853-945, 2006
[19] Feburie, V., Goit, M., Granger, S., Seyhaeve, J. M., “A Model for Chocked Flow through Cracks with Inlet Subcooling”, International J. Multiphase Flow, Vol. 19, No. 4, pp 541-562, 1993
84
[20] Hagoort, J., “Prediction of wellbore temperatures in gas production wells”, J. of Petroleum Science and Engineering, Vol. 49, pp. 22-36, 2005.
[21] Himmelblau, D., Bischoff, K., “Process analysis and simulation: deterministic systems”, USA, Wiley, 1967, pp 9-37
[22] Lee, J. and Wattenbarger, R. A., “Gas Reservoir Engineering”, SPE, Richardson, TX, pp. 1-28, 2004.
[23] Pattillo, P.D., Bellarby, J.E., Ross, G.R., Gosch, S.W., McLaren, G.D., “Thermal and Mechanical Considerations for Design of Insulated Tubing”, paper SPE 79870 presented at IADC/SPE Drilling Conference, Amsterdam, 19-21 February 2003.
[24] Ros, N. C. J., “Simultaneous Flow of Gas and Liquid as Encountered in Well Tubing”, J. of Petroleum Technology, pp. 1037-1049, October 1961.
[25] Yoshioka, K., Zhu, D., Hill, A.D, Dawkrajai, P., Lake, L.W., “A Comprehensive Model of Temperature Behavior in a Horizontal Well”, paper SPE 95656 presented at the 2005 SPE Annual Technical Conference and Exhibition, Dallas, 9-12 October 2005.
85
APPENDIX A: A DISCUSSION OF THE EQUILIBRIUM CONDITION FOR A MULTIPHASE FLUID SYSTEM
This section shows that a multiphase fluid system is at equilibrium condition if
the velocities of the phases are equal. A multiphase fluid system is defined to be
at equilibrium conditions when18:
• A liquid holdup equal to the liquid fractional flow.
• A void fraction equal to the gaseous fractional flow.
• A flowing density equal to the equilibrium density.
• A flowing velocity equal to the volumetric flux.
Recall the definitions for the mixture density, velocity, volumetric flux and
quality:
LLgg HH ρρρ += …………………………………………….……..(A-1)
LLgg
LLLggg
HHvHvH
vρρρρ
+
+= ………………………………………………(A-2)
AVu&
= ………………………………………………………………....(A-3)
gg Hxρρ
= ……………………………………………………………(A-4)
The mass flow rate can be written as:
LLLggg vAHvAHm ρρ +=& …………………………………………...(A-5)
86
By rearranging the equation A-5 and combining the equations A-1 and A-2, the
next relationships are obtained:
LLLggg vHvHv ρρρ += ………….……………………..………..…..(A-6)
Amv&
=ρ …………………………………………………………...…..(A-7)
The volumetric flow rate of the multiphase fluid system, the gas phase and the
liquid phases can be defined as:
…………………………………………………………...(A-8) Lg VVV &&& +=
ggg vAHV =& ..………………………………….……………………...(A-9)
…..……………………………….………………….....(A-10) LLL vAHV =&
By replacing the equations A-3, A-9 and A-10, the equation A-8 becomes:
…………………………….………...……....……(A-11) LLgg vHvHu +=
Assuming that vL=v* and vg= λv* where λ is the slip ratio, the equations A-6 and
A-11 take the form:
*)( vHHv LLgg ρλρρ += ……………..……………………………(A-11)
*)( vHHu Lg += λ …………………….……..……………………..(A-12)
87
By introducing the equations A-1 and A-4 into the equation A-11 and the
expression Hg+HL=1 into the equation A-12, the following relationships are