Неизотермическое течение жидкости
В общем случае в неизотермическом трубопроводе может наблюдаться два режима течения: на начальном участке - турбулентный, а в конце - ламинарный. Чтобы решить, какие режимы существуют в трубопроводе, следует рассчитать критическую температуру, при которой происходит смена турбулентного режима на ламинарный:
(7.13)
где t* - температура в рабочем диапазоне, при которой известна вязкость n*; u – коэффициент крутизны вискограммы, 1/К; d – внутренний диаметр трубопровода, м; Q – объемный расход, м3/с; Rе КР= 2300.
Коэффициент крутизны вискограммы рассчитываем по формуле:
. (7.14)
При tКР ³ tH режим только ламинарный, а при tКР £ tКрежим только турбулентный. При tH>tКР>tК в трубопроводе имеют место два режима.
Средняя по сечению трубопровода температура на любом расстоянии от его начала определяется по формуле Шухова:
. (7.15)
где tH, tO – температура нефти в начале трубопровода и температура окружающей среды, соответственно, оС; k – коэффициент теплопередачи, Вт/(м2*К); d – внутренний диаметр трубопровода, м; x – расстояние от начала трубопровода до точки, для которой рассчитывается температура, м; G – массовый расход нефти, кг/с; СР – удельная теплоемкость нефти, Дж/(кг*К).
В конце трубопровода Х = L.
Коэффициенты теплопередачи различны для ламинарного и турбулентного режимов.
Для точки Х, находящейся в ламинарной зоне при двух режимах движения нефти в трубопроводе, формула имеет следующий вид:
, (7.16)
где tКР – критическая температура, оС; LT – длина турбулентного участка, м; L - длина всего трубопровода, м; ШуЛ:
. (7.17)
Длина турбулентного участка рассчитывается по формуле:
. (7.18)
Температура в конце трубопровода при двух режимах движения нефти рассчитывается по следующей формуле:
. (7.19)
ТИПОВАЯ ЗАДАЧА «НЕИЗОТЕРМИЧЕСКИЕ ТРУБОПРОВОДЫ»
По трубопроводу длиной 143 км, внутренним диаметром 359 мм перекачивается нефть в количестве 11409 т/сут. Плотность нефти 950 кг/м3, кинематическая вязкость при 30 и 70 оС, соответственно 26.5*10-4 и 1.61*10-4 м2/с, удельная теплоемкость 2093.5 Дж / (кг*К). Коэффициенты теплопередачи Кт = 12.76 и Кл = 12.38 Вт/(м2*К). Температура подогрева нефти 69 оС. Температура окружающей среды 0 оС. В конце трубопровода температура нефти должна быть 33 оС. Тепловая изоляция отсутствует.
Определить:
- режимы движения в трубопроводе,
- потерю напора на трение (по формуле Лейбензона, без учета поправки на неизотермичность, для турбулентного режима принять зону Блазиуса)
- температуру нефти в конце трубопровода,
- температуру нефти на расстоянии 12 км от начала трубопровода.
Алгоритм решения задачи 2.
1. Определяем коэффициент крутизны вискограммы.
2. Рассчитываем критическую температуру и устанавливаем режим движения нефти в трубопроводе.
3. Рассчитываем длину турбулентного участка (если установлено два режима движения нефти в трубопроводе).
4. Рассчитываем температуру нефти в конце трубопровода.
5. Рассчитываем температуру нефти по длине трубопровода – 6 точек.
Рекомендации к решению задачи 2.
· Для определения температуры по длине неизотермического трубопровода при отсутствии парафиновых отложений воспользуйтесь формулами (18, 22). Не забудьте учесть режим движения жидкости на отдельных участках (турбулентный или ламинарный) через коэффициент теплопередачиki.
· Таблица для построения графика по результатам решения задачи 2 имеет следующий вид:
Форма таблицы для построения графика по результатам расчетов
Таблица 7.4
Параметр | Номера точек | Примечание | |||||||
L, м | LK | - | |||||||
t, оС | Задача 2 |
· Для расчетов и построения графиков воспользуйтесь программой EXCEL.
Дано: l = 300 т/сут; Q = 11409 т/сут; d = 359 мм; υ1 = 26.5*10-4 м2/с; t1 = 30 оС; υ2 = 1.61*10-4 м2/с; t2 = 70 оС; Kt = 12.76 Вт/м2/К; tн = 69 оС; Kл = 12.38 Вт/м2/К; tо = 0 оC; Ср = 2093.5 Дж / кг / К; tк = 33 оС; αн = 0.000594 1/ К; X = 12 км; ρ293 = 950 кг/м3.
Найти: lт =?, lл =?, hн =?, tкон =?, tx =?
Решение:
Определим коэффициент крутизны вискограммы
u =[ln(υ1 / υ2)] / (t2 - t1) = [ln(26.5 / 1.61)] / (70 - 30) = 0.07 1/К
Рассчитываем критическую температуру
tкр = t2 + 1 / n * ln(υ2 * π * d * Reкр / 4 / Q) = 70 + 1 / 0.007 * ln(1.61 * 10 -4 * 3.14 * 359 * 10 -3 * 2320) / 4 / 0.139 = 66 оС
если
tk < tkp < Tн, 33 < 66 < 69, то в трубопроводе два режима течения жидкости
Рассчитываем среднюю температуру потока на турбулентном участке
tпт = (tн + tкр) / 2 = (69 + 66) / 2 = 67.5 оС
На ламинарном участке: tпл = (tkp + tk) / 2 = (66 + 33) / 2 = 49.5 оС
Рассчитываем плотность нефти при средних температурах потока:
а) турбулентный участок
ρпт = ρ293 / (1 + αн * (tп - 20)) = 950 / (1 + 0.000594 * (67.5 - 20)) = 950 / 1.0282 = 923.93 кг/м3
в) ламинарный участок ρпл = 950 / (1 + 0.000594 * (49.5 - 20)) = 933.64 кг/м3
Рассчитываем длины участков:
а) турбулентный lt = Q * ρпт * Cр *ln((tн - to) / (tкр - tо)) / Kt / π / d = 0.139 * 923.9 * 2093.5 / 12.76 / 3.14 / 359 / 0.001 * ln(69 / 66) = 830.9 м
в) ламинарный lл = Q * ρпл * Cp / Кл / π / d * ln((tкр - to) / (tk - to)) = 0.139 * 933.6 * 2093.5 / 12.38 / 3.14 / 0.359 * ln(66 / 33) = 13494 м
с) сравниваем расчетные длины участков с заданной длиной трубопровода:
lt + lл = 830.9 + 13494 = 14325 м (lt + lл) = 14.3 км
т.е. заданная температура нефти в конце трубопровода может быть обеспечена без тепловой изоляции.
Расчет потери напора на трение:
а) рассчитываем вязкость при начальной температуре потока:
υн = υ2 * exp[-u * (tн - t2)] = 1.61 * 10 - 4 * exp[-0.07 * (69 - 70)] = 1.727 * 10 - 4 м2/с
в) рассчитываем потери напора на турбулентном участке:
hnт = β * Q (2 - m) * υн m * lт / d (5 - m) = 0.241 / 9.81 * 0.139 1.75 * (0.0001727) 0.25 * 830.9 / 0.359 4.75 = 9.6 м
рассчитываем потери напора на ламинарном участке
hnл = β * Q (2 - m) * υн m * (l - lт) / d (5 - m) = 128 / 3.14 / 9.81 * 0.139 * 0.0001727 * (14300 - 830.9) / 0.359 4 = 80.9 м
Сумма потерь
hnт + hnл = 9.6 + 80.9 = 90.5 м
Рассчитываем температуру нефти в конце трубопровода:
tкон = to + [(tн - to) * exp(-Kл * φ * d * l / Q / ρпл / Cp)] * ((tкр - to) / (tн - to)) (1 – Kл / Кт) = [69 * exp(-12.38 * 3.14 * 0.359 * 14300 / 0.139 / 933.6 / 2093.5)] * (66 / 69) (1 – 12..38 / 12.76) = 69 * exp(-0.73457) * 0.9565 0.02978 = 33.06 оС
Рассчитываем температуру нефти на расстоянии 12 км от начала трубопровода:
точка X находится на ламинарном участке трубопровода, т.к. l > X > lт
Поэтому пользуемся следующей формулой:
tx = to + (tкр - to) * exp[-Кл * π * d * l / Q / ρпл / Cp * ( X / l – lт / l )] = 66 * exp[-12.38 * 3.14 * 0.359 * 14300 / 0.139 / 933.6 / 2093.5 * (12 / 14.3 - 0.831 / 14.3)] = 37.18
Ответ: В трубопроводе 2 режима течения,
lт = 0.831 км; lл = 13.469 км; hnт = 90.5 м; tкон = 33.1 оС; tx = 37.2 оС;
ИСХОДНЫЕ ДАННЫЕ К ЗАДАНИЮ «НЕИЗОТЕРМИЧЕСКИЕ ТРУБОПРОВОДЫ»
ЗАДАЧА 7.3.1.
Определить режимы движения нефти в трубопроводе длины L и внутреннего диаметра d при определенной его пропускной способности Q, температуре нефти в начале tН и необходимой температуре в конце трубы tК. Температура окружающей среды tO известна. Тепловая изоляция отсутствует. Рассчитать температуру нефти по длине трубопровода (минимум 6 точек) и температуру нефти в конце трубопровода.
ИСХОДНЫЕ ДАННЫЕ К ЗАДАЧЕ 7.3.1
Таблица 7.5
Исходные данные | Варианты | |||||||||
L, км | 6,4 | 14,3 | ||||||||
d, мм | ||||||||||
Q, м3/с | 0,035 | 0,069 | 0,05 | 0,04 | 0,15 | 0,05 | 0,139 | 0,083 | 0,08 | 0,139 |
r293, кг/м3 | ||||||||||
tН, оС | ||||||||||
tК, оС | ||||||||||
tо, оС | -8 | -10 | ||||||||
CP, Дж/(кг*К) | ||||||||||
t1, оС | ||||||||||
t2, оС | ||||||||||
n1 *104, м2/с | 0,65 | 0,65 | 0,65 | 13,2 | 0,312 | 26,5 | 0,339 | 0,312 | 13,2 | |
n2 *104, м2/с | 0,17 | 2,78 | 0,17 | 0,17 | 3,24 | 0,066 | 1,61 | 0,076 | 0,066 | 3,24 |
kТ, Вт/(м2*К) | 14,91 | 5,4 | 14,72 | 14,91 | - | 15,18 | 12,76 | 12,99 | 15,18 | - |
kЛ, Вт/(м2*К) | 12,97 | 5,35 | - | 12,97 | 12,49 | - | 12,38 | - | - | 12,49 |
Ответ: tК, оС | 25,9 | 35,1 | 24,0 | 25,6 | 14,9 | 40,7 | 33,5 | Требуется изоляция |
ГАЗОПРОВОДЫ
ЗАДАЧА 7.4.1.
Известно отношение давлений Р1/Р2 в сечениях 1 и 2 газопровода постоянного диаметра. Течение изотермическое, известна скорость газа v1 , м/с. Найти v2.
ЗАДАЧА 7.4.2.
Определить массовый суточный расход газа, который можно передать по газопроводу, уложенному из труб диаметром d мм, на расстояние L км. Абсолютное давление газа на выкиде компрессорной станции P1МПа, в конце участка P2МПа, плотность газа rг при атмосферном давлении (0,1 МПа) и температуре перекачки 20 ° С. Газ считать совершенным, течение изотермическим.
Указание. Для расчета коэффициента гидравлического сопротивления вспользоваться формулой Веймаута.
ИСХОДНЫЕ ДАННЫЕ К ЗАДАНИЮ ПО ТЕМЕ «ГАЗОПРОВОДЫ»
Таблица 7.6
Задача 7.4.1. | ||||||||||
Параметр | Варианты | |||||||||
Р1/Р2 | 3,5 | 4,2 | 2,8 | 3,2 | 1,9 | 2,4 | 1,5 | |||
V1, м/с | ||||||||||
Задача 7.4.2. | ||||||||||
Диаметр газопровода, мм | ||||||||||
Длина газопровода, км | ||||||||||
Р1, МПа | 6,5 | 5,2 | 6,6 | 5,6 | 8,0 | 7,5 | 7,0 | 7,2 | 6,8 | 5,0 |
Р2, МПа | 4,7 | 3,4 | 5,5 | 4,5 | 6,8 | 5,0 | 4,6 | 5,2 | 3,8 | 3,2 |
ρг, кг/м3 | 0,86 | 0,80 | 0,74 | 0,88 | 0,9 | 0,68 | 0,78 | 0,62 | 0,70 | 0,72 |
ОТСТОЙНИКИ
ЗАДАЧА 8.1.
Как будет влиять температура на скорость процесса разделения нефти и воды методом отстаивания в резервуарах?
Рекомендации. Сравните влияние температуры на изменение плотностей минерализованной воды и нефти и на разность плотностей нефти и воды. Результат представьте графически.
Диапазон изменения температуры и характеристику воды принять по условиям задачи 3.1 («Физико-химические свойства пластовых и сточных вод»). Для расчета плотности нефти при заданных температурах воспользуйтесь формулой из раздела 2 для вычисления величины коэффициента термического расширения в зависимости от плотности нефти.
Результат представьте графически.
ИСХОДНЫЕ ДАННЫЕ К ЗАДАНИЮ ПО ТЕМЕ «ОТСТОЙНИКИ»
Таблица 8.1
Вариант | Плотность нефти при 20 оС, кг/м3 |