...
...
...
...
Motivation
...
One of the key problems in designing the pipelines and controlling the pipeline fluid transport is to predict the temperature and pressure losses during the stationary fluid transport.challenges in Pipe Flow Dynamics is to predict the pressure distribution along the pipe during the steady-state fluid transport.
In many practical cases the stationary pressure distribution can be approximated by Isothermal or Quasi-isothermal homogenous fluid flow model.
Pipeline Flow Pressure Model is addressing this problem with account of Pipeline flow simulator is addressing this problem. It should account for the varying pipeline trajectory, gravity effects , and fluid friction with pipeline walls and varying heat exchange with surroundings.
Definition
Given
...
Outputs
...
...
...
...
...
LaTeX Math Inline |
---|
body | \{ x_w(l), \, y_w(l), \, z_w(l) \} |
---|
|
...
| Flow velocity distribution along the pipe |
Inputs
...
...
| Fluid temperature at inlet point ( |
...
...
...
...
...
| Fluid pressure at inlet point ( |
...
...
...
...
...
...
...
...
...
...
...
...
...
...
...
--uriencoded--\displaystyle \cos \theta (l) |
|
|
...
Simulate
...
...
...
...
Assumptions
...
Steady-State flow | Quasi-isothermal flow |
LaTeX Math Inline |
---|
body | --uriencoded--\displaystyle \frac%7B\partial p%7D%7B\partial t%7D = 0 |
---|
|
| LaTeX Math Inline |
---|
body | --uriencoded--\displaystyle \frac%7B\partial T%7D%7B\partial t%7D =0 \rightarrow |
---|
|
|
...
...
...
user | ama@naftacollege.com |
---|
group | sofoil |
---|
Профиль давления
В процессе эксплуатации нагнетательной скважины движение флюида вдоль ствола
происходит в стационарном режиме, при этом профиль скорости потока и давления удовлетворяютусловию баланса массы движущегося потока:
LaTeX Math Block |
---|
anchor | MatBal2 |
---|
alignment | left |
---|
|
A(l) \, \rho(l) \, v(l) = \rm const |
и баланса сил действующих на единицу объема флюида в стволе скважины:
LaTeX Math Block |
---|
|
\frac{dp}{dl} = \rho \, g \, \sin \theta - \rho \, v \, \frac{dv}{dl} - \frac{ f \, \rho \, v^2 \, }{2 d} |
где
...
...
Image Removed
...
...
...
...
...
профиль поперечного сечения ствола скважины
LaTeX Math Inline |
---|
body | A(l) = 0.25 \, \pi \, d^2(l) |
---|
|
...
...
...
Эти замкнутая система уравнений для стационарного распределения давления и скорости потока вдоль трубы.
Уравнение
LaTeX Math Block Reference |
---|
|
часто в литературе записывают как разложение изменения давление вдоль ствола скважины на компоненты: LaTeX Math Block |
---|
anchor | gradP_General |
---|
alignment | left |
---|
|
\frac{dp}{dl} = \bigg( \frac{dp}{dl} \bigg)_g + \bigg( \frac{dp}{dl} \bigg)_v + \bigg( \frac{dp}{dl} \bigg)_f |
где
...
anchor | gradP_G |
---|
alignment | left |
---|
...
|
Homogenous flow | |
LaTeX Math Inline |
---|
body | --uriencoded--\displaystyle \frac%7B\partial p%7D%7B\partial \tau_x%7D =\frac%7B\partial p%7D%7B\partial \tau_y%7D =0 \rightarrow p(t, \tau_x,\tau_y,l) = p(l) |
---|
|
| |
Equations
...
LaTeX Math Block |
---|
| \left( \rho(p) - j_m^2 \cdot c(p) \right) \cdot \frac{dp}{dl} |
|
...
...
...
...
LaTeX Math Block |
---|
anchor | gradP_v |
---|
alignment | left |
---|
|
\bigg( \frac{dp}{dl} \bigg)_v = - \rho \, v \, \frac{dv}{dl}, |
кинетическая компонента вариация давления, формируемая вариацией скорости потока
вдоль ствола скважины, которая вызвана сжатием-расжатием флюида и изменением диаметра труб...
(l) - \frac{ j_m^2 }{2 d} \cdot f(p) |
| |
LaTeX Math Block |
---|
anchor | gradP_f |
---|
alignment | left |
---|
|
\bigg( \frac{dp}{dl} \bigg)_f = - \frac{ f \, \rho v^2}{2 d}, |
фрикционная компонента вариации давления, формируемая трением флюида со стенкой скважины
она всегда имеет отрциательный знак и приводит к потере давления вдоль направления движения потока
Для несжимаемой жидкости
в отсутствии трения уравнение LaTeX Math Block Reference |
---|
|
принимает вид: LaTeX Math Block |
---|
|
\frac{dp}{dl} = \rho \, g \, \sin \theta - \rho \, v \, \frac{dv}{dl} |
и может быть явно проинтегрировано:
...
и называется уравнением Бернулли.
...
title | Вывод уравнений движения флюида в стволе |
---|
...
Уравнение неразрывности одномерного потока с линейной плотностью
массы:...
...
...
\frac{\partial (\rho A)}{\partial t} + \frac{\partial}{\partial l} \big( A \, \rho \, v \big) = 0
для стационарного режима течения принимает вид:
...
u(l) = \frac{j_m}{\rho(l)} |
| |
...
...
\frac{\partial (\rho A)}{\partial t} = - \frac{\partial}{\partial l} \big( A \, \rho \, v \big) = 0
where
откуда и следует формула
LaTeX Math Block Reference |
---|
|
....
LaTeX Math Block Reference |
---|
|
...
...
LaTeX Math Inline |
---|
body | dm = \rho \, A \, dl |
---|
|
...
--uriencoded--\displaystyle j_m =\frac%7B \rho_0 \, q_0%7D%7BA%7D= \rm const |
|
| mass flux |
|
...
– сила гравитации, – сила трения со стенками трубы, – номральная реакция опоры стенок трубы.Рассмотрим стационарное (то есть установившееся во времени) течение потока по трубе.
Движение поперек трубы отсуствует и, следовательно, сумма проекций всех сил на трансверсальное направление
к трубе должно равняться нулю: LaTeX Math Block |
---|
|
dF_p \bigg |_{l_{\perp}} + dF_g \bigg |_{l_{\perp}} + dF_f \bigg |_{l_{\perp}}+ dF_N \bigg |_{l_{\perp}} =0 |
и выполняется автоматически, при наличии достаточного запаса прочности трубы
LaTeX Math Inline |
---|
body | dF_N \bigg |_{l_{\perp}} |
---|
|
.Уравнение движения флюида вдоль оси трубы
имеет вид: LaTeX Math Block |
---|
|
dF_p \bigg |_l + dF_g \bigg |_l + dF_f \bigg |_l+ dF_N \bigg |_l = \frac{d I}{dt}\bigg |_l |
где
представляет собой изменение импульса LaTeX Math Inline |
---|
body | I = \delta m \, v = \rho \, A \, \delta l \, v |
---|
|
элементарного объема флюида под действием внешних сил.Изменение импульса c учетом стационарности скорости потока
и сохранения массы LaTeX Math Inline |
---|
body | \frac{d (\delta m)}{dt}=0 |
---|
|
имеет вид: LaTeX Math Inline |
---|
body | \frac{dI}{dt} = \frac{d}{dt} (\delta m \, v) = \frac{d (\delta m \, v)}{\delta l} \frac{dl}{dt} = v \frac{ d (\delta m) \, v + \delta m \, dv}{\delta l} =v \frac{\delta m}{ \delta l} dv = \rho \, A \, v \, dv |
---|
|
.Сила, формируемая гидравлическим напором
LaTeX Math Inline |
---|
body | dF_p \bigg |_l = A (p - dp) - A p = - A \, dp |
---|
|
.Проекция гравитационной силы
LaTeX Math Inline |
---|
body | dF_g \bigg |_l = \delta m \, g \, \sin \theta = \rho \, A\, \delta l \, g \, \sin \theta |
---|
|
.Сила трения со стенками трубы дается феноменологическим уравнением Дарси-Вейсбаха:
LaTeX Math Inline |
---|
body | dF_f \bigg |_l = - \frac{f}{d} \frac{dm \, v^2}{2} = - \frac{f \, \rho \, v^2}{2 d} \, A \, \delta l |
---|
|
.Аксиальная компонента реакции опоры труб по определению отсутствует
.Подставляя вышеприведенные выражения в уравнение
LaTeX Math Block Reference |
---|
|
получим: LaTeX Math Block |
---|
|
- A dp + \rho \, A\, \delta l \, g \, \sin \theta - \frac{f \, \rho \, v^2}{2 d} \, A \, \delta l = \rho \, A \, v \, dv
|
Разделив уравнение на бесконечно малый объем элемента
получим LaTeX Math Block Reference |
---|
|
.Если дебит скважины на устье составляет
, а плотность воды на устье , то уравнение LaTeX Math Block Reference |
---|
|
можно записать в следующем виде: LaTeX Math Block |
---|
|
A \, \rho \, v = \rho_s \, q_s |
откуда можно выразить явно профиль скорости потока по стволу:
LaTeX Math Block |
---|
|
v(l) = \frac{\rho_s \, q_s}{\rho(p) \, A(l)} |
Подставляя
LaTeX Math Block Reference |
---|
|
в LaTeX Math Block Reference |
---|
|
получим уравнение на профиль давления вдоль ствола:...
| Fluid flowrate at inlet point () |
LaTeX Math Inline |
---|
body | \rho_0 = \rho(T_0, p_0) |
---|
|
| Fluid density at inlet point () |
LaTeX Math Inline |
---|
body | \rho(l) = \rho(T(l), p(l)) |
---|
|
| Fluid density at any point |
LaTeX Math Inline |
---|
body | --uriencoded--\displaystyle с(p) = \frac%7B1%7D%7B\rho%7D \left( \frac%7B\partial \rho%7D%7B\partial p%7D \right)_T |
---|
|
| Fluid Compressibility |
LaTeX Math Inline |
---|
body | --uriencoded--f(T, \rho) = f(%7B\rm Re%7D(T, \rho), \, \epsilon) |
---|
|
| Darcy friction factor |
LaTeX Math Inline |
---|
body | --uriencoded--\displaystyle %7B\rm Re%7D(T,\rho) = \frac%7Bj_m \cdot d%7D%7B\mu(T, \rho)%7D |
---|
|
| Reynolds number in Pipe Flow |
| dynamic viscosity as function of fluid temperature and density |
LaTeX Math Inline |
---|
body | --uriencoded--\displaystyle d = \sqrt%7B \frac%7B4 A%7D%7B\pi%7D%7D= \rm const |
---|
|
| Characteristic linear dimension of the pipe (or exactly a pipe diameter in case of a circular pipe) |
Alternative forms
...
LaTeX Math Block |
---|
| \frac{dp}{dl} = \ |
|
...
Далее учтем, что угол наклона к горизонту
может быть выражен через абсолютные отметки глубин вдоль траектории скважины : LaTeX Math Block |
---|
|
\sin \theta = \frac{dz}{dl} |
и уравнение для давление примет вид:
...
...
left( \frac{dp}{dl} \right)_G + \left( \frac{dp}{dl} |
|
...
\right)_K + \left( \frac{ |
|
...
...
where
...
...
LaTeX Math Block |
---|
anchor | dp_implicit |
---|
alignment | left |
---|
|
\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 Block |
---|
|
\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 |
---|
body | c(p)= \frac{1}{\rho} \frac{d \rho}{dp} |
---|
|
– сжимаемость воды и уравнение профиля давления принимает вид: LaTeX Math Block |
---|
anchor | dp_explicit |
---|
alignment | left |
---|
|
\bigg( 1 - \frac{c(p) \, \rho_s^2 \, q_s^2}{A^2} \bigg ) \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 Block Reference |
---|
|
представляет собой обыкновенное дифференциальное уравнение первого порядка на функцию со слабой нелинейностью.Если предположить постоянство коэффициента трения
и несжимаемость флюида LaTeX Math Inline |
---|
body | \rho(p) = \rho_s = \rm const |
---|
|
, то уравнение LaTeX Math Block Reference |
---|
|
можно явно проинтегрировать: LaTeX Math Block |
---|
|
p(l) = p_s + \rho \, g \, z(l) - \frac{\rho_s \, q_s^2 }{2 A^2 d} \, f_s \, l |
Pressure gradient will be:
LaTeX Math Block |
---|
|
\frac{dp}{dl} = \cos \theta(l) - \frac{\rho_s \, q_s^2 }{2 A^2 d} \, f_s |
where
LaTeX Math Inline |
---|
body | \cos \theta(l) = \frac{dz(l)}{dl} |
---|
|
The first term defines the hydrostatic column of static fluid while the last term defines the friction losses under fluid movement:
LaTeX Math Block |
---|
|
\frac{dp}{dl} \Bigg|_{loss} = \frac{\rho_s \, q_s^2 }{2 A^2 d} \, f_s |
В калькуляторе Well Flow Performance Calculator можно оценить величину потерь на трения для различных сценариев диаметров труб и дебитов скважин.
Как видно из вышеприведенных корреляций, коэффициент трения меняется в зависимости от скорости потока и соответствующего числа Рейнольдса.
Основным вкладом в вариабельность коэффициента трения вдоль трубы является диаметр трубы в данной точке траектории скважины, который может приводить к значительным изменениям скорости потока.
Тем не менее, зависимость от дебита является слабой. Из формулы
LaTeX Math Block Reference |
---|
|
видно что изменение дебит в 10 раз приводит к изменению коэффициента трения в раз....
--uriencoded--\displaystyle \left( \frac%7Bdp%7D%7Bdl%7D \right)_G = \rho \cdot g \cdot \cos \theta |
|
| gravity losses which represent pressure losses for upward flow and pressure gain for downward flow |
LaTeX Math Inline |
---|
body | --uriencoded--\displaystyle \left( \frac%7Bdp%7D%7Bdl%7D \right)_K = u%5e2 \cdot \frac%7Bd \rho%7D%7Bdl%7D |
---|
|
| kinematic losses, which grow contribution at high velocities and high fluid compressibility (like turbulent gas flow) |
LaTeX Math Inline |
---|
body | --uriencoded--\displaystyle \left( \frac%7Bdp%7D%7Bdl%7D \right)_f = - \frac%7B j_m%5e2%7D%7B2 d%7D \cdot \frac%7Bf%7D%7B\rho%7D |
---|
|
| friction losses which are always negative along the flow direction |
Approximations
See also
...
Зависимость коэффициента трения от давления формируется только через число Рейнольдса:
....
LaTeX Math Inline |
---|
body | {\rm Re} = \frac{d \, \rho \, v}{\mu} |
---|
|
...
LaTeX Math Block Reference |
---|
|
...
LaTeX Math Block |
---|
|
{\rm Re} = \frac{ d \, \rho_s \, q_s}{A \, \mu(p)} |
отсюда следует, что зависимость коэффициента трения от давления формируется вязкостью
, которая для воды имеет слабую зависисмость от давления в широких практических пределах:δμ/μ = 25 % при вариации μ = 2.4·10-5 Па · с для p = 1 атм до μ = 3.0·10-5 Па · с для 300 атм (cм. Свойства воды).
...
Для оценки числа Рейнольдса для нагнетаемой по 2.5 " НКТ воды можно пользоваться формулой
LaTeX Math Inline |
---|
body | {\rm Re} = 230 \cdot \, q |
---|
|
, где дебит скважины на устье в м3/сут.Отсюда видно, что при дебитах более 18 м3/сут число Рейнольдса становится больше 4,000 и режим течения является турбулентным и коэффициент трения можно считать практически постоянным вдоль ствола нагнетательной скважины.
А учитывая, что рост давления с глубиной сопровождается увеличением температуры, что компенсирует рост вязкости воды, то для большинства практических реализаций ППД можно полагать, что вариация коэффициента трения вдоль ствола не превышает 2-3 % и в оценках потери напора на трение принимать коэффициент трения постоянным
.See Also
Petroleum Industry / Upstream / Pipe Flow Simulation
...
solverbook.com – Коэффициент теплоотдачи
Show If |
---|
|
Panel |
---|
bgColor | papayawhip |
---|
title | ARAX |
---|
| PipeFlow |
|
References
...