Page tree

Versions Compared

Key

  • This line was added.
  • This line was removed.
  • Formatting was changed.


The base driving equations of a pipe flow are:

Steady-state 1D inviscid fluid flowPipe Flow Mass Conservation


LaTeX Math Block
anchor
9QRCZ
p
alignmentleft
\
bigg( 1 - \frac{c(p) \, \rho_s^2 \, q_s^2}{A^2} \bigg ) \frac{dp}{dl} =
frac{d p}{d l} =
-\rho \, u \, \frac{d u}{d l}  + \rho \, g \, \cos \theta + f_{\rm cnt, \, l}



LaTeX Math Block
anchorjm
alignmentleft
j_m(l) = j_m = \rho(l) \cdot u = \rm const


Equation of State (EOS)Darcy–Weisbach


LaTeX Math Block
anchorIY1I1
alignmentleft
\rho = \rho(p, T)



LaTeX Math Block
anchorf
alignmentleft
f_{\rm cnt, l} =  -  f \cdot \frac{  \rho \, u^2 \, }{2 d}


where

LaTeX Math Inline
bodyl

distance along the fluid flow streamline

LaTeX Math Inline
body\theta(l)

inclinational deviation,  

LaTeX Math Inline
body\displaystyle \cos \theta = dz/dl

LaTeX Math Inline
bodyz(l)

elevation along the 1D flow trajectory 

LaTeX Math Inline
bodyT(l)

fluid temperature

LaTeX Math Inline
bodyp(l)

fluid pressure

LaTeX Math Inline
body\rho(l)

fluid density

LaTeX Math Inline
body--uriencoded--%7B\bf u%7D(l)

fluid velocity vector 

LaTeX Math Inline
bodyu(l)

superficial velocity of the pipe flow

LaTeX Math Inline
body--uriencoded--%7B\bf f%7D_%7B\rm cnt%7D(l)

volumetric density of all contact forces exerted on fluid body

LaTeX Math Inline
body--uriencoded--f_%7B\rm cnt, l%7D(l) = %7B\bf e%7D_u \cdot %7B\bf f%7D_%7B\rm cnt%7D

projection of 

LaTeX Math Inline
body--uriencoded--%7B\bf f%7D_%7B\rm cnt%7D
 onto the unit fluid velocity vector 
LaTeX Math Inline
body--uriencoded-- %7B\bf e%7D_u = %7B %7C %7B\bf u%7D %7C%7D %5e%7B-1%7D \, %7B\bf u%7D

LaTeX Math Inline
bodyj_m = \rho(l) \cdot u(l)

fluid mass flux

LaTeX Math Inline
body\dot m

mass flowrate

LaTeX Math Inline
bodyg

standard gravity constant


Substituting 

LaTeX Math Block Reference
anchorjm
 and 
LaTeX Math Block Reference
anchorf
 into 
LaTeX Math Block Reference
anchorp
:

LaTeX Math Block
anchorXZ51K
alignmentleft
\frac{d p}{d l} =
-j_m \cdot \frac{d}{d l} \left( \frac{j_m}{\rho} \right)   + \rho \, g \, \

...

cos \theta -  f \cdot \frac{  \rho

...

  \,

...

 }{2 

...

d} \cdot \left( \frac{

...

j_m}{\rho} \right)^2


LaTeX Math Block
anchor

...

HERBX
alignmentleft

...

\frac{d p}{d l} =
j^2_m \cdot \frac{1}{\rho^2} \frac{d \rho}{dl}   + \rho \, g \, \cos \theta -  \frac{j_m^2}{2 d} \cdot \frac{f}{\rho}


LaTeX Math Block
anchor

...

ICKH6
alignmentleft

...

\frac{d p}{d l} =
j^2_m \cdot \frac{1}{\rho^2} \frac{d \rho}{dp} \cdot \frac{d p}{dl}   + \rho \, g \, \cos \theta -  \frac{j_m^2}{2 d} \cdot \frac{f}{\rho}


LaTeX Math Block
anchordrhoJ9867
alignmentleft
-\frac{d p}{d l} =
j^2_m \cdot \frac{1}{c}\rho} \cdot c \cdot \frac{d p}{dl}   + \rho \, g \, \cos \theta -  \frac{j_m^2}{2 d} \cdot \frac{f}{\rho}

and finally

LaTeX Math Block
anchor91WI7
alignmentleft
\left( 1 - j_m^2 \cdot \frac{1c}{\rho}   \right ) + \frac{dp}{dl} = \rho_s^2 q_s^2}{2A^2} \, g \, \cos \theta  - \frac{j_m^2 }{2  d} \cdot \frac{df}{dl} \left(\rho}


Alternative forms

...


LaTeX Math Block
anchorrho
alignmentleft
\left[ \rho -j_m^2 \, c \right] \cdot \frac{

...

d p}{dl} =
\rho^2

...

 \, g \, \cos \theta -  \frac{j_m^2 }{

...

2d} \cdot f(\rho)



LaTeX Math Block
anchorrho
alignmentleft
\left[ \frac{1}{c} - \frac{j_m^2}{\rho} \right] \cdot \frac{d \rho}{dl} =
\rho^2 \, g \, \cos \theta -  \frac{j_m^2 }{2d} \cdot f(\rho)




See Also

...

Petroleum Industry / Upstream / Pipe Flow Simulation / Water Pipe Flow @model / Stationary Isothermal Homogenous Pipe Flow Pressure Profile @model

Darcy friction factor ] [ Darcy friction factor @model ]

Euler equation ]



Show If
groupeditors


Panel
bgColor#FFDFDD


Expand
titleEditor

В процессе эксплуатации нагнетательной скважины движение флюида вдоль ствола

LaTeX Math Inline
bodyl(x,y,z)
происходит в стационарном режиме, при этом профиль скорости потока
LaTeX Math Inline
bodyv(l)
и давления
LaTeX Math Inline
bodyp(l)
удовлетворяют

условию баланса массы движущегося потока:

LaTeX Math Block
anchorMatBal2
alignmentleft
 A(l) \, \rho(l) \, v(l) = \rm const

и баланса сил действующих на единицу объема флюида в стволе скважины:

LaTeX Math Block
anchorgradP
alignmentleft
\frac{dp}{dl} = \rho \, g \, \cos \theta - \rho \, v \, \frac{dv}{dl} - \frac{ f \, \rho \, v^2 \, }{2 d}

где

LaTeX Math Inline
bodyl

длина ствола скважины, отсчитываемая вниз от поверхности


LaTeX Math Inline
body\rho(l)
 

профиль плотности воды

LaTeX Math Inline
body\theta(l)

Inclinational Deviation

LaTeX Math Inline
bodyd(l)

профиль диаметра скважины, вдоль которого идет поток

LaTeX Math Inline
bodyA(l)

профиль поперечного сечения ствола скважины

LaTeX Math Inline
bodyA(l) = 0.25 \, \pi \, d^2(l)

LaTeX Math Inline
bodyf(l)

профиль коэффициента трения Дарси

LaTeX Math Inline
bodyg

ускорение свободного падения ( = 9.87 м2/сек )




Эти замкнутая система уравнений для стационарного распределения давления и скорости потока вдоль трубы.


Для несжимаемой жидкости

LaTeX Math Inline
body\rho = \rm const
в отсутствии трения
LaTeX Math Inline
body f = 0
уравнение
LaTeX Math Block Reference
anchorgradP
принимает вид:

LaTeX Math Block
anchorgradP
alignmentleft
\frac{dp}{dl} = \rho \, g \, \cos \theta - \rho \, v \, \frac{dv}{dl}

и может быть явно проинтегрировано:

LaTeX Math Block
anchorgradP
alignmentleft
p(l) - \rho \, g \, l \, \cos \theta +  \frac{1}{2} \rho \, v^2 = \rm const

и называется уравнением Бернулли.




Если дебит скважины на устье составляет

LaTeX Math Inline
bodyq_s
, а плотность воды на устье
LaTeX Math Inline
body\rho_s
, то уравнение
LaTeX Math Block Reference
anchorMatBal2
можно записать в следующем виде:

LaTeX Math Block
anchorArhov
alignmentleft
A \, \rho \, v = \rho_s \, q_s

откуда можно выразить явно профиль скорости потока по стволу:

LaTeX Math Block
anchorv
alignmentleft
v(l) = \frac{\rho_s \, q_s}{\rho(p) \, A(l)}


Подставляя 

LaTeX Math Block Reference
anchorv
  в  
LaTeX Math Block Reference
anchorgradP
 получим уравнение на профиль давления вдоль ствола:

LaTeX Math Block
anchor6JNN1
alignmentleft
\frac{dp}{dl} = \rho \, g \, \cos \theta - \frac{\rho_s^2 \, q_s^2}{A}  \frac{d}{dl} \bigg( \frac{1}{A \, \rho} \bigg) - \frac{\rho_s^2 \, q_s^2 }{2 A^2 d} \frac{f}{\rho}


Далее учтем, что угол наклона к горизонту

LaTeX Math Inline
body\theta
может быть выражен через абсолютные отметки глубин  
LaTeX Math Inline
bodyz(l)
  вдоль траектории скважины 
LaTeX Math Inline
bodyl(x,y,z)
:

LaTeX Math Block
anchor9114L
alignmentleft
\cos \theta = \frac{dz}{dl}

и уравнение для давление примет вид:

LaTeX Math Block
anchor6JNN1
alignmentleft
\frac{dp}{dl} = \rho \, g \, \frac{dz}{dl} -  \frac{\rho_s^2 \, q_s^2}{A}  \frac{d}{dl} \bigg( \frac{1}{A \, \rho} \bigg) - \frac{\rho_s^2 \, q_s^2 }{2 A^2 d} \frac{f}{\rho}


Диаметр труб, вдоль которых идет движение воды, остается постоянным на долгом протяжении и меняется редко (например, километр НКТ и потов выход потока в колонну), и это позволяет решать задачу нахождения профиля давления на кусках постоянного диаметра 

LaTeX Math Inline
bodyd = {\rm const}, \quad A = {\rm const}
 и уравнение может быть переписано следующим образом:

LaTeX Math Block
anchordp_implicit
alignmentleft
\frac{dp}{dl} = \rho \, g \, \frac{dz}{dl} - \frac{\rho_s^2 \, q_s^2}{A^2}  \frac{d}{dl} \bigg( \frac{1}{\rho} \bigg) - \frac{\rho_s^2 \, q_s^2 }{2 A^2 d} \frac{f}{\rho}


Процесс движения воды вдоль трубы происходит в состоянии термодинамического равновесия и плотность воды является функцией только давления

LaTeX Math Inline
body\rho = \rho(p)
и, следовательно:

LaTeX Math Block
anchorAQKCW
alignmentleft
\frac{d}{dl} \bigg( \frac{1}{\rho} \bigg) = -\frac{1}{\rho^2} \frac{d \rho}{ dl} 
= - \frac{1}{\rho^2}\frac{d \rho}{dp} \frac{dp}{ dl}
=- \frac{c}{\rho} \frac{dp}{ dl}

где

LaTeX Math Inline
bodyc(p)= \frac{1}{\rho} \frac{d \rho}{dp}
– сжимаемость воды и уравнение профиля давления принимает вид:

LaTeX Math Block
anchordp_explicit
alignmentleft
\biggleft( 1 -  \frac{c(p) \, \rho_s^2 \, q_s^2}{\rho \, A^2}   \biggright )  \frac{dp}{dl} = \rho(p) \, g \, \frac{dz}{dl}  - \frac{\rho_s^2 \, q_s^2 }{2 A^2 d} \frac{f(p)}{\rho(p)}


Функция

LaTeX Math Inline
bodyz(l)
определяется траекторией скважины.

Cжимаемость
LaTeX Math Inline
bodyc(p)
и плотность
LaTeX Math Inline
body\rho(p)
воды слабо зависят от вариации давления вдоль ствола.

Как будет показано ниже коэффициент трения

LaTeX Math Inline
bodyf(p)
тоже слабо зависит от вариации давления и, следовательно, уравнение
LaTeX Math Block Reference
anchordp_explicit
представляет собой обыкновенное дифференциальное уравнение первого порядка на функцию
LaTeX Math Inline
bodyp(l)
со слабой нелинейностью.

See Also

Petroleum Industry / Upstream / Pipe Flow Simulation / Water Pipe Flow @model / Stationary Isothermal Homogenous Pipe Flow Pressure Profile @model

...