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.
Pipeline flow simulator is addressing this problem. It should account for the varying pipeline trajectory, gravity effects, fluid friction with pipeline walls and varying heat exchange with surroundings.
Definition
Given
- space coordinates are
\{ x, \, y, \, z \} with
z-ccordinate facing down to the Earth Centre
- inflow pipeline coordinates
\{ x_s = 0, \, y_s = 0, \, z_s = 0 \}
- pipeline trajectory
\{ x_w(l), \, y_w(l), \, z_w(l) \}, where
l = \int_0^l \sqrt{dx^2 + dy^2 + dz^2} = \int_0^l \sqrt{\dot x^2 + \dot y^2 + \dot z^2} dl, is pipeline length from inflow point
\{ x_s = 0, \, y_s = 0, \, z_s = 0 \}
- pipeline cross-section area
A(l)
- earth gravity vector
{\bf g} = (0, \, 0, \, g) where
g = 9.81 \ \rm m/s^2
- inflow temperature
T_s, inflow pressure
p_s, inflow rate
q_s
- PVT properties of water
\rho(T, p),
\mu(T, p)
- surroundings initial temperature
T_g(l), thermal diffusivity
a_e(l), thermal conductivity
\lambda_e(l) of surrounding media
- heat exchange coefficient
U(l) based on pipeline schematics
Simulate
- along-pipe distribution of stabilized pressure
p(l), flow rate
q(l) and average flow velocity
u(l)
- along-pipe distribution of fluid flow temperature
T(t, l) after a flow period of time
t when the flow stabilization achieved
Как видно из вышеприведенных корреляций, коэффициент трения меняется в зависимости от скорости потока и соответствующего числа Рейнольдса.
Основным вкладом в вариабельность коэффициента трения вдоль трубы является диаметр трубы в данной точке траектории скважины, который может приводить к значительным изменениям скорости потока.
Тем не менее, зависимость от дебита является слабой. Из формулы
видно что изменение дебит в 10 раз приводит к изменению коэффициента трения в 10^{0.25} = 1.8 раз.
Еще более слабой является вариабельность коэффициента трения от давления вдоль ствола, что можно проиллюстрировать следующими соображениями.
Зависимость коэффициента трения от давления формируется только через число Рейнольдса: f = f(\rm Re(p)).
При этом число Рейнольдса
{\rm Re} = \frac{d \, \rho \, v}{\mu} с учетом
(1) | {\rm Re} = \frac{ d \, \rho_s \, q_s}{A \, \mu(p)} |
отсюда следует, что зависимость коэффициента трения от давления формируется вязкостью f = f(\mu(p)), которая для воды имеет слабую зависисмость от давления в широких практических пределах:
δμ/μ = 25 % при вариации μ = 2.4·10-5 Па · с для p = 1 атм до μ = 3.0·10-5 Па · с для 300 атм (cм. Свойства воды).
Это приводит к 25 % вариации коэффициента трения для ламинарного потока (в котором сила трения минимальна) и порядка 4.5 % для турбулентного потока (и максимальным вкладом трения).
Для оценки числа Рейнольдса для нагнетаемой по 2.5 " НКТ воды можно пользоваться формулой {\rm Re} = 230 \cdot \, q , где q дебит скважины на устье в м3/сут.
Отсюда видно, что при дебитах более 18 м3/сут число Рейнольдса становится больше 4,000 и режим течения является турбулентным и коэффициент трения можно считать практически постоянным вдоль ствола нагнетательной скважины.
А учитывая, что рост давления с глубиной сопровождается увеличением температуры, что компенсирует рост вязкости воды, то для большинства практических реализаций ППД можно полагать, что вариация коэффициента трения вдоль ствола не превышает 2-3 % и в оценках потери напора на трение принимать коэффициент трения постоянным f = f_s = \rm const.
References
https://en.wikipedia.org/wiki/Darcy_friction_factor_formulae
https://neutrium.net/fluid_flow/pressure-loss-in-pipe/
H. J. Ramey, Wellbore Heat Transmission - SPE-96-PA - 1992
R. Shankar, Pipe Flow Calculations, Clarkson University
solverbook.com – Коэффициент теплоотдачи