# Mathematical modeling of oil transportation by pipelines using anti-turbulent additives

## Marijonas Bogdevičius1, Jolanta Janutėnienė2, Kęstutis Jonikas3, Eleonora Guseinovienė4, Mindaugas Drakšas5

1, 4University of Klaipėda, Mechatronics Science Institute, Bijūnų str. 17, Klaipėda, LT-91225, Lithuania

1, 3Vilnius Gediminas Technical University, Plytinės str. 27, Vilnius, LT-10105, Lithuania

2, 5University of Klaipėda, Mechanical Engineering Department, Bijūnų str. 17, Klaipėda, LT-91225, Lithuania

Journal of Vibroengineering, Vol. 15, Issue 1, 2013, p. 419-427.
Received 15 January 2013; accepted 28 February 2013; published 31 March 2013

Copyright © 2013 Vibroengineering 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.
Views 144
Abstract.

The energy costs of oil transportation by pipeline are one of the most important questions for oil companies. Large amounts of energy are used in the oil pumping process. One of the ways to reduce energy costs is to use additives, which can change the rheological characteristics of oil. Anti-turbulent polymer additives reduce the hydraulic resistance and allow a significant reduction in pumping costs. The hydrodynamic processes in the pipeline and efficiency of oil transportation costs, dependent on the concentration of anti-turbulent additives, are investigated in this paper. The pressure and velocity of the flow in each point of the pipeline are obtained by solving a system of nonlinear algebraic equations, which were solved utilizing the Newton-Raphson method. The dependencies of the efficiency coefficient in the oil pumping process were deduced, depending on the concentration of the additives.

Keywords: oil transportation, additives, numerical method.

#### 1. Introduction

Oil transportation by pipeline requires a lot of energy. The load is constantly changing in the oil pump motors and consumes vast quantities of electricity. The optimization of oil pumping regimes increases the oil throughput and electric motor efficiency.

Sections of the oil pipeline in which inner sand, wax or tar aggregations have formed are a major contributor to increased energy consumption.

This problem is important when transporting high viscosity, heavy and extra-heavy crude oil. Innovative technological solutions concerned with viscosity and friction reduction to move such crude oils from the production site to the processing facilities are presented in the work . The stability and viscosity of water-oil emulsions and their application for heavy oil pipeline transportation were investigated in the work . Factors affecting the properties and stability of these emulsions were investigated and were limited to 60 % for crude oil content in the emulsions. Different methods for reducing the viscosity of heavy crude oil to enhance the flow properties were investigated in the work . The influence of shear rate, temperature and light oil concentration on the viscosity behavior was studied. It was observed that the blending of heavy crude oil with a limited amount of lighter crude oil provides better performance than the other alternatives.

One of ways to improve the process of oil transportation is to use a reagent, which modifies the paraffin wax properties in a way so that it does not adhere to the pipe walls. The usage of anti-turbulent polymer additives, which reduce the hydraulic resistance in oil during its flow, allows a significant reduction in pumping costs . The additives are substances which reduce the hydraulic resistance and increase the throughput of the oil pipeline . Some of the additives are high molecular weight polymers. The polymer additives increase the flow rate at a constant pressure . The most widely used additives are polymers such as polybutadiene, polyisoprene, polystyrene, its derivatives and others.

The object under investigation in the paper  was the reduction of the turbulent friction phenomenon in pipes when a slight amount of polymer additives was added to the pumped fluid. Experimental studies were carried out with oil transported by pipes of different diameters, using polymer additives . Investigation has shown the friction coefficient’s dependence on the Reynolds number and the influence of friction on the fluid flow rate. The authors  have found that the usage of high molecular weight α-olefin polymers increases the oil transportation capacity by 23 percent.

Cases were studied in the paper , where the introduction of 10 to 30 grams of additives to one metric ton of oil resulted in a 15-25 % increase in pipeline throughput, depending on the pipe diameter. The most common poly-α-olefins were used for the reduction in oil hydrodynamic resistance.

An aim of this work is to mathematically investigate the efficiency of oil transportation by pipeline, depending on the concentration of anti-turbulent additives.

#### 2. Mathematical model of the hydrodynamic system

Next we represent the oil pumping system by a pipeline (Fig. 1) and its dynamic model (Fig. 2). The hydrodynamic system being analyzed consists of centrifugal pumps and a long pipeline (Fig. 2).

Fig. 1. Altitude of oil transportation system Fig. 2. Hydrodynamics system The equation of the fluid continuity can be written in a differential form as follows:

(1)
$\frac{\partial }{\partial t}\left[S\left(x\right)\rho \right]+\frac{\partial }{\partial x}\left[S\left(x\right)\rho v\right]={F}_{1}\left(x\right)\text{,}$

where $\rho$, $v$ are the density and velocity of the fluid, $S\left(x\right)$ is the cross section area of the pipeline, ${F}_{1}\left(x\right)$ is the discharge of the fluid mass to the unit of length in the pipeline. The equation of the fluid flow impulse (momentum):

(2)
$\frac{\partial }{\partial t}\left[S\left(x\right)\rho v\right]+\frac{\partial }{\partial x}\left[S\left(x\right)\left(p+\rho {v}^{2}\right)\right]+\Pi \left(x\right)\tau +S\left(x\right)\rho {a}_{x}={F}_{2}\left(x\right)+p\frac{\partial S}{\partial x},$

where $p$ is the fluid pressure, $\mathrm{\Pi }\left(x\right)$ is the perimeter of the cross section of the pipeline, $\tau$ is the tangential fluid stress in the inner surface of the pipeline, ${a}_{x}$ is the acceleration along the $x$ axis, ${F}_{2}\left(x\right)$ is the kinetic energy of the fluid flow in the pipeline to the unit of area. The tangential fluid stress in the inner surface of the pipeline is determined:

(3)
$\tau =\rho \lambda \left(\mathrm{R}\mathrm{e},\epsilon \right)\frac{v|v|}{8}.$

Anti-turbulent fluid additives reduce the friction coefficient $\lambda$. This coefficient can be calculated from the equation :

(4)
$\frac{1}{\sqrt{\lambda }}=0.88\cdot \mathrm{l}\mathrm{n}\left(k\left(\theta \right)\mathrm{R}\mathrm{e}\sqrt{\lambda }\right)-3.745,$

where $\theta$ is the concentration of the turbulent additives.

The friction coefficient $\lambda$ is determined by solving the nonlinear algebraic equation:

(5)
$\Phi \left(\lambda \right)=\frac{1}{\sqrt{\lambda }}-0.88\cdot \mathrm{l}\mathrm{n}\left(k\left(\theta \right)\mathrm{R}\mathrm{e}\sqrt{\lambda }\right)+3.745=0.$

This equation is solved by the Newton-Raphson method:

(6)
$\Delta {\lambda }_{k}=-{\left(\frac{\partial \Phi \left({\lambda }_{k}\right)}{\partial \lambda }\right)}^{-1}\Phi \left({\lambda }_{k}\right),\text{\hspace{0.17em}\hspace{0.17em}\hspace{0.17em}\hspace{0.17em}}{\lambda }_{k+1}={\lambda }_{k}+\Delta {\lambda }_{k}.$

The dependence of hydraulic resistance coefficient $\lambda$ on the Reynolds number and the FLO-XL additive concentration is presented in Fig. 3.

Fig. 3. Dependence of hydraulic resistance coefficient $\lambda$ on the Reynolds number and the FLO-XL additive concentration Fig. 4. Dependence of coefficient $k$ on the FLO-XL additive concentration The results obtained from the calculation (Table 1) showed the dependence:

$k\left(\theta \right)=28.21+2.052\cdot \theta +0.246\cdot {\theta }^{2}-2.469\cdot {10}^{-3}\cdot {\theta }^{3},$ which is given in Fig. 4.

Table 1. Calculated values of coefficients $k\left(\theta \right)$ for additives FLO-XL

 $\theta$, ml/m³ 0 5.2 11.3 19 28.8 37.5 $k\left(\theta \right)$ 28 46 78 140 232 321

The system of equations (1) and (2) can be written as a system of second-order quasi-linear differential equations:

(7)
$\left[A\right]\left\{\frac{\partial u}{\partial t}\right\}+\left[B\right]\left\{\frac{\partial u}{\partial x}\right\}=\left\{f\right\},$

where $\left[A\right]$, $\left[B\right]$ are matrices and $\left\{f\right\}$ is a vector which depends on $t$, $x$ and the elements ${u}_{i}$ of vector ${\left\{u\right\}}^{\mathrm{T}}=\left[p,v\right]$:

Equating the determinant of matrix $\left[B\right]-\left[A\right]\frac{dx}{dt}$ to zero, the equation which allows determination of the $\frac{dx}{dt}$ derivative and characteristic direction is obtained. This equation has two various real roots $dx/dt={\lambda }_{i}$, $\left(i=1,2\right):$

(8)
${C}^{+}:\frac{dx}{dt}=v+a,\text{\hspace{0.17em}\hspace{0.17em}\hspace{0.17em}\hspace{0.17em}}{C}^{-}:\frac{dx}{dt}=v-a,$

$a$ – the sound velocity in a liquid with a certain amount of gas, which is stored in an elastic pipeline, is equal to:

(9)
$a=\sqrt{\frac{K\left(p\right)/\rho }{1+\frac{K\left(p\right)\cdot d}{E\cdot e}+\frac{\epsilon }{\gamma }\left[\frac{K\left(p\right)}{\gamma p}-1\right]}},$

where: $K\left(p\right)$ – bulk modulus of the elasticity of the liquid, $\rho$ – density of the liquid, $E$ – modulus of elasticity of the pipeline, – internal diameter of a pipeline, $e$ – thickness of the wall of a pipeline, $\gamma$ – index of adiabatic process, $\epsilon$ – ratio of the gas volume in the liquid and the total volume of the liquid (mixture).

Differential equations of liquid movement in the pipeline are solved by the characteristics method . The unknowns, variable velocity and pressure of the liquid at any specific moment in time $t+\Delta t$, are determined according to these parameters at the moment of time $t$ (Fig. 5). The length of the pipeline is divided into elements, the length of which is $\Delta x$. Pressure and velocity at a point $D$ at the moment of time $t$ are determined from the nonlinear algebraic equation system:

(10)
${C}^{+}:{\Phi }_{1}={v}_{D}-{v}_{L}+\frac{1}{2}\left({p}_{D}-{p}_{L}\right)\left[{\left(\frac{1}{\rho c}\right)}_{L}+{\left(\frac{1}{\rho c}\right)}_{D}\right]$
(11)
${C}^{-}:{\Phi }_{1}={v}_{D}-{v}_{R}+\frac{1}{2}\left({p}_{D}-{p}_{R}\right)\left[{\left(\frac{1}{\rho c}\right)}_{R}+{\left(\frac{1}{\rho c}\right)}_{D}\right]$

where: ${f}_{1}\left(v,p\right)=\frac{{c}_{2}}{S\left(x\right)}\left({F}_{1}-\rho v\frac{\partial S}{\partial x}\right)$; ${f}_{2}\left(v,p\right)=\frac{1}{\rho S\left(x\right)}\left({F}_{2}-v{F}_{1}\right)-\frac{\tau \Pi \left(x\right)}{\rho S\left(x\right)}-{a}_{x}$; $\tau =\frac{\rho f\left(\mathrm{R}\mathrm{e}\right)v|v|}{8}$; $\lambda \left(\mathrm{R}\mathrm{e}\right)$ is the coefficient of pressure losses along the pipe; $\text{Re}$ is the Reynolds number.

The equation system was solved:

(12)
$\left[B\left(u\right)\right]\left\{\frac{\partial u}{\partial x}\right\}=\left\{f\left(u\right)\right\},$

and the distributions of the initial pressure and velocity were determined.

This system of equations was solved by the finite element method and by using the Galerkin method.

Vector $\left\{u\right\}$ in the finite element is approximated by:

(13)
$\left\{u\right\}=\left[N\left(x\right)\right]\left\{{q}_{e}\right\}=\left[N\left(\xi \right)\right]\left\{{q}_{e}\right\},$

where: $\left[N\left(\xi \right)\right]=\left[1-\xi ,\xi \right]$, $\xi =x/{L}_{e}$; ${L}_{e}$ – length of the finite element;

${\left\{q\right\}}_{e}^{\mathrm{T}}=\left[\begin{array}{cccc}{p}_{e1}& {v}_{e1}& {p}_{e1}& {v}_{e1}\end{array}\right].$

Fig. 5. Scheme of determination of point parameters by characteristic method The equation system for the finite element was determined by applying the Galerkin method:

(14a)
$\underset{0}{\overset{1}{\int }}{\left[N\left(\xi \right)\right]}^{\mathrm{T}}\left(\left[B\left({q}_{e}\right)\right]\left[\frac{dN\left(\xi \right)}{{L}_{e}d\xi }\right]\left\{{q}_{e}\right\}-\left\{f\left({q}_{e}\right)\right\}\right){L}_{e}d\xi =0,$

or:

(14b)
$\left[{A}_{e}\left({q}_{e}\right)\right]\left\{{q}_{e}\right\}=\left\{{b}_{e}\left({q}_{e}\right)\right\},$
(15)
$\left[{A}_{e}\left({q}_{e}\right)\right]=\underset{0}{\overset{1}{\int }}{\left[N\left(\xi \right)\right]}^{\mathrm{T}}\left[B\left({q}_{e}\right)\right]\left[\frac{dN\left(\xi \right)}{d\xi }\right]d\xi ,\text{\hspace{0.17em}\hspace{0.17em}\hspace{0.17em}\hspace{0.17em}}\left\{{b}_{e}\left({q}_{e}\right)\right\}=\underset{0}{\overset{1}{\int }}{\left[N\left(\xi \right)\right]}^{\mathrm{T}}\left\{f\left({q}_{e}\right)\right\}{L}_{e}d\xi .$

The general equation system of the oil pipeline was generated from the finite elements equation system (14):

(16)
$\left[A\left(q\right)\right]\left\{q\right\}=\left\{b\left(q\right)\right\}.$

The nonlinear equation system (16) was solved using the Newton-Raphson method:

(17)
$\left[{J}_{k}\right]{\left\{\Delta q\right\}}_{k}=-\left\{\Phi \left({q}_{k}\right)\right\},\text{\hspace{0.17em}\hspace{0.17em}\hspace{0.17em}}{\left\{q\right\}}_{k+1}={\left\{q\right\}}_{k}+{\left\{\Delta q\right\}}_{k}.$

The power of the friction force in the inner surface of the pipeline is equal:

(18)
${N}_{frict}=\sum _{e=1}^{N}{N}_{e}=\sum _{e=1}^{N}\underset{xe,1}{\overset{xe,2}{\int }}v\tau dA={A}_{e}\sum _{e=1}^{N}\frac{1}{3}\left({v}_{e,2}-{v}_{e,1}\right)\left({\tau }_{e,2}-{\tau }_{e,1}\right)$

The coefficient of efficiency of the pipeline is equal:

(19)
$\eta =\frac{{N}_{input}-{N}_{output}-{N}_{frict}}{{N}_{input}-{N}_{output}},$

where ${N}_{input}$, ${N}_{output},$ are the power at the input and output of the pipeline, respectively.

#### 3. Theoretical and numerical analysis

The main parameters of the investigated pipeline are the length of the pipeline $l=$ 91.6 km, outer diameter of the pipe $D=$ 0.559 m, wall thickness of the pipe $\delta =$ 0.0095 m and inner diameter of the pipe $d=$ 0.540 m. Temperature of the oil being pumped: highest during the summer, +25°C and lowest in the winter, –6°C. Efficiency of the system – 1400 m³/h, maximum admissible pressure in the pipeline – ${P}_{max}=$ 10.2 MPa. Density of the transported oil – 860 kg/m³ (at 20°C).

Oil pumping is done by six pumps with a nominal shaft rotation speed of 2980 rpm, powered by three-phase electric motors of 900 and 2000 kW with an efficiency of 0.96.

The general characteristic of the pumps and piping at the first pump intake pressure ${H}_{p}=$ 13.2 m are shown in Fig. 6. The horizontal line in Fig. 6 shows the pressure limit of the section of investigated pipeline. The intersection point of this line with the characteristic curve of the pump shows the lowest capacity at which the relevant number of pumps can operate.

The intersection points of the characteristic curves of the pumps and piping show the capacity of the hydraulic system and the pressure developed by the pumps. The pipeline capacity at any intersection point of the curves can be calculated by the formula:

(20)
$Q=\sqrt{\frac{{H}_{p}+\sum _{i=1}^{n}{a}_{i}-\left({z}_{0}-{z}_{L}\right)-{H}_{L}}{\sum _{i=1}^{n}{b}_{i}+\frac{\lambda \cdot L}{2\cdot d\cdot g\cdot {S}^{2}\cdot {3600}^{2}}}},$

where: ${H}_{p}$ – the inlet pressure of the first pump; ${H}_{L}$ – the pressure at a distance $L$ from the beginning of the pipeline; ${z}_{L}$ – position of geometric height; $a$ and $b$ – coefficients of characteristics; $\lambda$ – coefficient of resistance. The pipeline is divided into 18320 elements, the length of each pipeline element is $\Delta x$ = 5 m.

The selection of an appropriate integration time step $\Delta t$ is very important to the accuracy of the analysis. In this case, a time step of $\Delta t=$ 0.00025 s was selected for use in the integration process. The pressure and velocity at each pipeline point was obtained by solving a system of nonlinear algebraic equations. The system was solved by using the Newton-Raphson method. The duration of the hydrodynamic processes for this simulation was chosen to be 160 s. Initial pressure distribution along the pipeline is shown in Fig 7. Oil pressure and discharge at the corresponding pipeline points $x$ are shown in Fig. 8 to Fig 13.

Oil pressure distributions after 160 s in the pipeline are shown in Fig. 14. The power of the friction forces in the inner surface of the pipeline with different concentrations of the additives are shown in Fig. 15. The efficiency coefficient of the pipeline is shown in Fig. 16 and the comparative energy dependency on oil pipeline productivity – in Fig. 17.

The relative energy costs of transporting one metric ton of oil decrease with increasing additive concentration.

Fig. 6. General characteristic of the pumps and piping Fig. 7. Initial pressure distribution in the pipeline Fig. 8. Dependencies of the pressure in the different pipeline points when concentration of additives in the oil is equal 0 ml/m3 Fig. 9. Dependencies of the discharge in the different pipeline points when concentration of the additives in oil is equal 0 ml/m3 Fig. 10. Dependencies of the pressure in the different pipeline points when the concentration of the additives in the oil is equal to 12 ml/m3 Fig. 11. Dependencies of the discharge in the different pipeline points when the concentration of the additives in the oil is equal to 12 ml/m3 Fig. 12. Dependencies of the pressure in the different pipeline points when the concentration of the additives in the oil is equal to 50 ml/m3 Fig. 13. Dependencies of the discharge in the different pipeline points when the concentration of the additives in the oil is equal to 50 ml/m3 Fig. 14. Distributions of pressure in the pipeline Fig. 15. Dependencies of the power of friction forces in the inner surface of the pipeline Fig. 16. Dependencies of the pipeline efficiency coefficient Fig. 17. Comparative energy dependence on oil pipeline productivity #### 4. Conclusions

The methodology developed allows an assessment of the technological methods of oil transportation by pipelines in consideration of the usage of additives for the hydrodynamic processes.

A numerical simulation of the hydrodynamic system has shown that increasing the additive concentration in oil increases pipeline capacity and decreases the relative energy consumption. However, the cost for the oil additive amount in transporting one metric ton of oil increases.

The mathematical model and numerical analysis created enable the investigation of fluid transportation systems, the selection of optimal operation regimes, determination of the necessary number of pumps, and an evaluation of the anti-turbulent additives’ effects on optimal system performance.

#### Acknowledgements

This work has been supported by the Research Council of Lithuania within the Project “Simulation Software and Investigation of Thermo-Hydrodynamic Processes in the Geothermal Loop”, Project No. MIP-090/2012.

1. Martínez-Palou M. et al. Transportation of heavy and extra-heavy crude oil by pipeline: A review. Journal of Petroleum Science and Engineering, Vol. 75, Issues 3-4, 2011, p. 274-282. [Search CrossRef]
2. Ashrafizadeh S. N., Kamran M. Emulsification of heavy crude oil in water for pipeline transportation. Journal of Petroleum Science and Engineering, Vol. 71, Issues 3-4, 2010, p. 205-211. [Search CrossRef]
3. Shadi W. Hasan, Mamdouh T. Ghannam,Nabil Esmail Heavy crude oil viscosity reduction and rheology for pipeline transportation. Fuel, Vol. 89, Issue 5, 2010, p. 1095-1100. [Search CrossRef]
4. Dmitrieva T. V. The main steps in the application of chemical additives and chemicals in the pipeline transportation of oil. Chemicals, Reagents and Small Volume Chemistry Processes. Ufa: Reaktiv, 2001, p. 152-153, (in Russian). [Search CrossRef]
5. Dujmovich T., Gallegos A. Drag reducers improve throughput, cut costs. Offshore, Vol. 65, Issue 12, 2005, p. 1-4. [Search CrossRef]
6. Chichkanov S. V., Myagchenkoyuv V. A. Toms effect – prospective applications. Herald of Kazan Technological University, Vol. 2, 2003, p. 314-329, (in Russian). [Search CrossRef]
7. Rakhmatullin Sh. I., Gareev M., Kim D. P. On a turbulent flow of a weakly concentrated polymer solutions in pipelines. Oil and Gas Business, Vol. 2, 2005, p. 1-6, (in Russian). [Search CrossRef]
8. Aguilar G., Gasljevic K., Matthys E. F. Reduction of friction in fluid transport: experimental investigation. Revista Mexicana de Fisica, Vol. 52, Issue 5, 2006, p. 444-452. [Search CrossRef]
9. Suleimanova J. V., Nesyn G. V. Technology and application of polysulfide α-olefin based elastomers. Materials of International Conference “Chemistry, Chemical Technology and Biotechnology”, 2006, p. 473-474, (in Russian). [Search CrossRef]
10. Nesyn G. V., Suleimanova J. V. et al. Anti-turbulen α-olefin based suspension type additives. Bulletin of the Tomsk Polytechnic University, Vol. 3, 2006, (in Russian). [Search CrossRef]
11. Ishmukhamedov I. T. et al. Pipeline Transportation of Petroleum Products. Oil and Gas. 1999, p. 299, (in Russian). [Search CrossRef]
12. Aladjev V., Bogdevičius M., Prendkovskis O. New Software for Mathematical Package Maple Releases 6, 7 and 8. Vilnius, Technika, 2002. [Search CrossRef]