Advances in Chemical Engineering and Science, 2013, 3, 9-15

http://dx.doi.org/10.4236/aces.2013.31002 Published Online January 2013 (http://www.scirp.org/journal/aces)

Numerical Analysis of Heavy Oil-Water Flow and Leak
Detection in Vertical Pipeline
Joo Victor Nunes de Sousa
1*
, Cristiane Holanda Sodr
2
, Antonio Gilson Barbosa de Lima
1
,
Severino Rodrigues de Farias Neto
3

1
Department of Mechanical Engineering, Center of Sciences and Technology,
Federal University of Campina Grande, Campina Grande, Brazil
2
Department of Chemical Engineering, Technology Center, Federal University of Alagoas, Macei, Brazil
3
Department of Chemical Engineering, Center of Sciences and Technology,
Federal University of Campina Grande, Campina Grande, Brazil
Email:
*
joao.vns@hotmail.com

Received September 3, 2012; revised October 5, 2012; accepted October 14, 2012
ABSTRACT
Pipeline is a conventional, efficient and economic way for oil transportations. The use of a good system for detecting
and locating leaks in pipeline contribute significantly to operational safety and cost saving in petroleum industry. This
paper aims to study the heavy oil-water flow in vertical ducts including leakage. A transient numerical analysis, using
the ANSYS-CFX

11.0 commercial software is performed. The mathematical modeling considers the effect of drag and
gravitational forces between the phases and turbulent flow. Mass flow rate of the phases in the leaking orifice, the pres-
sure drop as a function of the time and the velocity distributions are presented and discussed. We can conclude that
volumetric fraction of phases and fluid mixture velocity affect pressure drop and mass flow rate at the leak hole.

Keywords: Heavy Oil; Multiphase Flow; Leakage; Numerical Simulation
1. Introduction
The activity of oil production is subject to high risks.
Even the petroleum industry running preventive meas-
ures, there is always the possibility of failure, making the
industrial plants susceptible to operational accidents with
loss of fluid to environment, causing great ecological,
social and economic damages, with delay in oil produc-
tion. A proper supervisory system must be capable of
detecting leaks in oil installations, enabling immediate
action to reduce the impacts of accidents and contribut-
ing significantly to operational safety. The simultaneous
flow of two immiscible liquids in vertical pipes is en-
countered in different industrials processes and particu-
larly in the petroleum industry [1].
Because of importance, many authors have focused
their researches in methods of leak detection in pipes on
oil production and transport [2-4].
However, in different applications, including oil trans-
portation, accurate locations of the leaks is still very dif-
ficult. In present day various leak detection techniques
based in the negative pressure wave, acoustic sensors,
satellite surveillance, mass and volume balance, analyti-
cal model-based method, among others, has been applied.
All these methods are based in process variables such as
pressure, mass and volumetric flow rates and temperature
[5].
According to Dong et al. [6], the negative pressure
method, which supplies high leak sensitivity and avail-
ability, is a relatively better method among them. Unfor-
tunately this method has a high possibility of false alarm
if there are some strong raises in the pressure measure-
ment records or if the leak is small (0.5% of nominal
flow) [4,5].
Thus, this paper aims to numerically study the hydro-
dynamic of heavy oil-water flow in a vertical pipe having
a small leak, which is much more difficult to detect by
conventional systems [7]. The interest in heavy oil is in
fact that recent studies indicate that in 2025 this kind of
oil will be the main source of fossil energy in the world
[8].
2. Methodology
2.1. The Geometry and Grid
The study domain (Figure 1) consists of a vertical pipe
with 800 cm (8 m) of length, with a constant circular
section 15 cm diameter. To simulate the leakage, the pipe
has a circular hole, with 0.6 cm diameter, located at the
*
Corresponding author.
Copyright 2013 SciRes. ACES
J. V. N. DE SOUSA ET AL. 10
midpoint of the length of pipe.
Figure 2 illustrates the mesh representing the study
domain, which was built with the support of ICEM-
CFD

11.0 software. This structured mesh was obtained

after various refinements, and it has 327,327 hexahedral
elements.
2.2. Mathematical Modeling
To investigate the multiphase flow of heavy oil-water in
vertical pipes leaking, it was considered three-dimen-
sional, transient and isothermal flow, non-homogeneous
model for the fluid mixture (particle model) and homo-
geneous model for turbulence (k- model). The homoge-
neous model considers a single field for both phases,
while in the non-homogeneous model is considered a
specific field for each phase [9].
On the mathematical modeling, the index represent
the continuous phase (oil) and the index represent the


Figure 1. Considered pipe layout.


Figure 2. Studied pipe, showing different regions of the
mesh.
dispersed phase (water). The dispersed water phase is
modeled as spherical particles.
The general equations used in this work are:
- Continuity Equations,
( )
( ) 0,
f
f
t
o o
o o o

c
+ V =
c
U (1)
( )
( )
0,
f
f
t
| |
| | |

c
+V =
c
U (2)
where f is the volume fraction, is the density and
( ) , , u v w = U is the velocity vector, each corresponding
to a given phase;
- Momentum Equations,
( )
( )
( )
{ }
T
,
,
eff M
f
f
t
f p f
o o o
o o o o
o o o o o o

c
( + V

c
(
' = V + V V + V + +

U
U U
U U S D
(3)
( )
( )
( )
{ }
T
,
,
eff M
f
f
t
f p f
| | |
| | | |
| | | | | |

c
(
+V

c
(
' = V +V V + V +
(

U
U U
U U S D
(4)
where p' is the modified pressure,
eff
, is the effective
viscosity, S
M
is the momentum sources due to external
body forces (when gravitational forces are includes) and
D is the drag force between the phases, that is modeled
by the equation
(
1
,
8
D
C A
o | o | o
= D U U U
)
U (5)
where C
D
is the drag coefficient and A is the interfacial
area density. For Re < 1000, the drag coefficient is mod-
eled by Schiller-Naumann model,
(
0.687
24
1 0.15 ,
D
C R
Re
= +
)
e (6)
and for Re 1000, the drag coefficient is considered 0.44,
where Re represents the particle Reynolds number, mod-
eled by
,
d
Re
o | o |
o

=
U U
(7)
where d

is the diameter of spherical particles (3 mm).

The interfacial area density, A, is modeled by the equa-
tion
6
;
f
A
d
|
|
= (8)
- Kinetic Energy Equation,
( )
( ) ,
m t
m m m k m
k
k
k k
t C

P
( c | |
+ V = V + V +
( |
c
( \ .
U c
(9)
Copyright 2013 SciRes. ACES
J. V. N. DE SOUSA ET AL. 11
where
m
is the density, k is the kinetic energy,
m
is
the velocity vector,
m
is the viscosity,
t
is the turbulent
viscosity, P
k
is the turbulence production and is the
turbulence eddy dissipation, each corresponding to a
mixture;
U
- Turbulence Eddy Dissipation Equation,
( )
( )
( )
2 3
1
,
m
m m
t
m k
t
C P C
C k
c c
c
c
c
c
m
c
c
+ V
c
( | |
= V + V +
( |
( \ .
U
c
(10)
where C
k
= 1.0, C
1
= 1.3, C
2
= 1.44 and C
3
= 1.92. In
the k- model the parameters as follows
2
,
3
din m
p p k ' = + (11)
,
,
eff t o o
= + (12)
,
,
eff t | |
= + (13)
,
m
f f
o o | |
= + (14)
( )
1
,
m
f f
o o o | | |

= + U U U (15)
,
m
f f
o o | |
= + (16)
2
,
t m
k
C

c
= (17)
( ) (
T
2
3
3
k t t
P = V V + V V V + U U U U U ), k (18)
where p
din
is the dynamic pressure,

is the viscosity
of phase and C

= 0.09. The dynamic pressure is calcu-

lated by the equation
( )
2
2 1
.
2
din
p f f
o o o | | |
= + U U (19)
In multiphase flow, the total pressure acting in the
phases, p
tot
, is modeled by
,
tot st din
p p p = + (20)
where p
st
is a term correspondent to the static pressure.
2.3. Initial and Boundary Conditions and Fluids
Properties
Initially (at time t = 0), the leak in the pipe does not exist,
and the multiphase flow occurs in steady state condition.
When t > 0 the leak appears abruptly in the pipeline.
Tables 1 and 2 shows the initial and boundary conditions
used in the simulations.
Simulations were performed with the following situa-
tions:
1) Oil volume fraction at the inlet, f
o
, ranging from
0.75 to 1.00 (step 0.05) and the water volume fraction at
Table 1. Boundary conditions for t = 0.
t = 0
Boundary
p [Pa] U [m/s] f
o
[-] f
w
[-]
Inlet - 0.75 - 2.00 0.75 - 1.00 0.00 - 0.25
Outlet 101,325 - - -
Leakage
a
- 0 - -
Wall - 0 - -
a
Impermeable.

Table 2. Boundary conditions for t > 0.
t > 0
Boundary
p [Pa] U [m/s] f
o
[-] f
w
[-]
Inlet - 0.75 - 2.00 0.75 - 1.00 0.00 - 0.25
Outlet 101,325 - - -
Leakage
a
101,325 - - -
Wall - 0 - -
a
Opening.

the inlet, f
w
, corresponds to the remaining fraction, with
the velocity mixture at the inlet, U, fixed in 1.00 m/s;
2) Mixture velocity, U, ranging from 0.75 to 2.00 m/s
(step 0.25), with the water volume fraction, f
w
, and oil
volume fraction, f
o
, fixed in 0.15 and 0.85, respectively.
The adopted properties of the fluids are showed in Ta-
ble 3.
Table 4 shows the considerations adopted for the nu-
merical solver.
3. Results and Discussion
The simulations were performed on a Intel

Core 2 Quad
2.4 GHz, 4 GB RAM and 1 Tb physical memory (HD)
computer. The solving time of the studied cases ranged
from 10 to 11 hours.
Figure 3 shows the total mass flow rate relationships
on the leakage, m
leak
/m
in
, as a function of oil volume frac-
tion at the inlet section, f
o,in
, to post-leak system. We can
see that the smaller the oil volume fraction, consequently,
greater the water volume fraction in the mixture, thus we
have a large amount of fluids exiting in the leak. This
behavior is explained by the reduction of the viscosity of
the fluid mixture, proportional to the water holdup con-
tained in it. Figure 4 shows the total and oil mass flow
rates in the leak, m
leak
, as a function of the oil volume
fraction at the inlet section of pipe, to post-leak system.
Note that the amount of water and oil mass flow rate in
leak is proportional to the holdup of phases present in the
mixture. Equations (21)-(23) were obtained by fitting the
bared on the numerical data.
Copyright 2013 SciRes. ACES
J. V. N. DE SOUSA ET AL. 12
Table 3. Fluids properties.
Fluid [kg/m
3
] [mPas] [N/m]
a

Water 997.0 0.8899
Oil 925.5 500
0.02905
a
is the oil-water surface tension.

Table 4. Numerical simulation characteristics.
Characteristic Consideration
Flow regime Transient
Total simulation time 0.100 s
Time step 0.001 s
Convergence criterion for mass
and momentum
10
7
(RMS)
Advection scheme High resolution
Transient scheme Second order
Interpolation scheme for pressure Trilinear
Interpolation scheme for velocity Trilinear


Figure 3. Total mass flow rate relationship in the leakage as
a function of the oil volume fraction at the pipe inlet (U =
1.00 m/s).

leak in ,in
2 3
,in ,in
1.602813447 5.122188379
6.011162647 2.480859856 .
o
o o
m m f
f f
=
+
(21)
leak, total ,in
2 3
,in ,in
262.9121767 834.7170159
974.6134762 401.0288889 .
o
o o
m f
f f
=
+
(22)
leak,oil ,in
2
,in ,in
236.8741911 828.5900693
1036.792032 443.2948148 .
o
o
m

Figure 4. Total and oil mass flow rates in the leakage as a
function of the oil volume fraction at the pipe inlet (U = 1.00
m/s).

Figure 5 shows the total mass flow rate relationships
in leakage, m
leak
/m
in
, as a function of the fluid mixture
velocity at the inlet, U
m,in
, to post-leak system. It is noted
that with increasing velocity, there is a reduction in the
total mass flow rate through in the leaking. However,
there is a reversal point (U
m,in
= 1.50 m/s), where by in-
creasing the velocity of flow, we have an increased total
mass flow rate in the leaking. This reversal point would
be the point where the inertial forces of the flow becomes
less than the forces caused by pressure differential be-
tween the inside and outside of the pipeline leaking
(pressure gradients). Figure 6 shows the mass flow rate
in the leak, m
leak
, as a function of the mixture velocity at
the inlet, to post-leak. It is visible that increasing the
mixture velocity at the inlet, the oil mass flow in the
leaking increase, and the water mass flow is practically
constant. Equations (24)-(26) were obtained by fitting the
bared on the numerical data.
leak in ,in
2 3
,in ,in
4
,in
0.2331399521 0.3712444597
0.2980465586 0.1059246175
0.01427994068 .
m
m m
m
m m U
U U
U
=
+
+

(24)
leak,total ,in
2 3
,in ,in
12.19420476 8.811055026
9.052387302 1.213925926 .
m
m m
m U
U U
=
+
(25)
3
o
f
f f
=
+
(23)
leak,oil ,in
2 3
,in ,in
10.14426714 7.349817778
7.577278095 1.010542222 .
m
m m
m U
U U
=
+
(26)
Table 5 shows the determination coefficients, R
2
, ob-
tained in the fittings correspondents to Equations (21) to
Copyright 2013 SciRes. ACES
J. V. N. DE SOUSA ET AL. 13

Figure 5. Total mass flow rate relationship in the leakage as
a function of the fluid mixture velocity at the pipe inlet (f
o,in

= 0.85).


Figure 6. Total and oil mass flow rates in the leakage as a
function of the fluid mixture velocity at the pipe inlet (f
o,in
=
0.85).

(26). These equations were obtained by the method of the
least squares.
Figure 7 shows the oil velocity vectors at the leakage
section, for some oil holdups in the mixture at the pipe
inlet (f
o,in
= 0.80, 0.90 and 1.00). Note that by reducing
the oil holdup, and thus, increasing the water holdup,
there is a greater spread in the leakage, with a larger an-
gle, taking as reference to the longitudinal direction of
the pipe.
Figure 8 shows the velocity vectors at the leakage
section, for some fluid mixture velocities at the pipe inlet
Table 5. Determination coefficients in the fittings.
Equation R
2

21 0.999837
22 0.999845
23 0.999652
24 0.999995
25 0.999974
26 0.999974


,
0.80
o in
f =
,
0.90
o in
f =
,
1.00
o in
f =

Figure 7. Velocity vectors at the leakage section, for some
oil volume fractions in the mixture at the pipe inlet (U
m,in
=
1.00 m/s).


,
1.00 m/s
m in
U =
,
1.50 m/s
m in
U =
,
2.00 m/s
m in
U =

Figure 8. Velocity vectors at the leakage section, for some
fluid mixture velocities at the pipe inlet (f
o,in
= 0.85).

(U
m,in
= 1.00, 1.50 and 2.00 m/s). It is visible the small
influence of the velocity in the leak spread and angle.
Figure 9 shows the total pressure drop, p
tot
, as a
function of the time, t, for different oil volume fractions
at the inlet, for a pipe section of 4 m length, where the
leak is in the midpoint of this section. In all cases, it is
visible reduction in the pressure drop at the initial in-
stants of the leakage. The transient period is short (less
than 0.03 s), as expected, due to low mass flow rate
through out the leak hole. Reduction in the pressure drop
is greater with the increase of the water holdup in the
mixture. After the flow reaches a new steady state (t >
0.025 s), the pressure drop reach the same level before
the leakage. This level is greater for the cases where the
Copyright 2013 SciRes. ACES
J. V. N. DE SOUSA ET AL. 14

Figure 9. Pressure drop as a function of the time, for dif-
ferent oil volume fractions in the mixture at the pipe inlet
(U
m,in
= 1.00 m/s).

water volume fraction in the mixture is greater, since the
water has a density higher than the oil, increasing the
static pressure drop portion, which is the most represen-
tative portion of the total pressure drop.
Figure 10 shows pressure drop as a function of the
time, for different mixture velocities at the inlet, for a
pipe section of 4 m length, where the leakage is located
in the midpoint of this section. Similarly, it is visible the
reduction in the pressure drop in the initial instants of the
leakage. The transient period is too short, being less than
0.025 s. The higher pressure drop in transient state is
obtained as the mixture velocity at the inlet is increased.
Figures 11 and 12 show, respectively, the total pres-
sure and the oil velocity fields near the leakage section
(in the transversal plan), for a case with the mixture ve-
locity in the inlet 1.00 m/s and the oil and water volume
fractions are 0.85 and 0.15, respectively. We can see that
the leakage region is one zone of low pressure and high
velocity. This behavior was found in all cases analyzed
in this paper.
4. Conclusions
In this paper the hydrodynamic of two-phase flow in ver-
tical pipe with a leakage is discussed. The study is related
to heavy oil-water flow in the turbulent regime by using
the ANSYS-CFX

11.0 commercial software.

The simulations revealed the difficulty of detecting
small leaks (less than 1% of the total mass flow rate

Figure 10. Pressure drop as a function of the time, for dif-
ferent mixture velocities at the pipe inlet (f
o,in
= 0.85).


Figure 11. Pressure field near the leakage.


Figure 12. Oil velocity field near the leakage.
Copyright 2013 SciRes. ACES
J. V. N. DE SOUSA ET AL.
Copyright 2013 SciRes. ACES
15
transported in the pipeline), since the time interval of
pressure drop observed is very short (less than 0.03 s).
This research revealed too that in heavy oil-water flow it
is more easy to detect leaks when the pipe operate with
high fluid velocities and greater water volume fraction in
the mixture at the pipe inlet. It was verified that the oil
volume fraction in the mixture affect strongly the spread
and the exit direction of the oil in the leakage hole.
5. Acknowledgements
The authors would like to express their thanks to Brazil-
ian researches agencies ANP/UFCG-PRH-25, CNPq,
CAPES, FINEP and PETROBRAS S/A for supporting this
work, and are also grateful to the authors of the refer-
ences in this paper, that helped in the improvement of
quality.
REFERENCES
[1] J. Y. Xu, D. H. Li, J. Guo and Y. X. Wu, Investigations
of Phase Inversion and Frictional Pressure Gradients in
Upward and Downward OilWater Flow in Vertical Pi-
pes, International Journal of Multiphase Flow, Vol. 36,
No. 11-12, 2010, pp. 930-939.
doi:j.ijmultiphaseflow.2010.08.007
[2] S. I. Kam, Mechanistic Modeling of Pipeline Leak De-
tection at Fixed Inlet Rate, Journal of Petroleum Science
and Engineering, Vol. 70, No. 3-4, 2010, pp. 145-156.
doi:j.petrol.2009.09.008
[3] X. H. Lu, Y. J. Sang, J. Z. Zhang and Y. Y. Fan, A Pipe-
line Leakage Detection Technology Based on Wavelet
Transform Theory, Proceedings of IEEE Annual Inter-
national Conference on Information Acquisition, Shandong,
20-23 August 2006, pp. 1432-1437.
doi:10.1109/ICIA.2006.305966
[4] R. Hu, H. Ye, G. Wang and C. Lu, Leak Detection in
Pipelines Based on PCA, Proceedings of 8th Interna-
tional Conference on Control, Automation, Robotics and
Vision, Kunming, 6-8 December 2004, pp. 1985-1989.
doi:10.1109/ICARCV.2004.1469466
[5] H. V. Silva, C. K. Morooka, I. R. Guilherme, T. C. Fon-
seca and J. R. P. Mendes, Leak Detection in Petroleum
Pipelines Using a Fuzzy System, Journal of Petroleum
Science and Engineering, Vol. 49, No. 3-4, 2005, pp.
223-238. doi:j.petrol.2005.05.004
[6] L. Dong, S. Chai and B. Zhang, Leak Detection and
Localization of Gas Pipeline System Based on Wavelet
Analysis, Proceedings of 2nd IEEE International Con-
ference on Intelligent, Control and Information Process-
ing, Harbin, 25-28 July 2011, pp. 478-483.
doi:10.1109/ICICIP.2011.6008290
[7] C. Verde, Multi-Leak Detection and Isolation in Fluid
Pipelines, Control Engineering Practice, Vol. 9, No. 6,
2001, pp. 673-682. doi:10.1016/S0967-0661(01)00026-0
[8] C. G. Moth and C. Silva, Heavy-oilReserves and
World Production, TN Petrleo, Vol. 57, 2007, pp. 76-81.
http://www.tnpetroleo.com.br/download.php/revista/dow
nload/i/33/nome/TN57_Artigos.pdf
[9] ANSYS, CFX-Solver Theory Guide (Release 11.0), AN-
SYS, Inc., Canonsburg, 2006.