Слайд 1Магистральные трубопроводы
Лекция №8
Слайд 2ТЕХНОЛОГИЧЕСКИЙ РАСЧЕТ МАГИСТРАЛЬНОГО ГАЗОПРОВОДА
Магистральные трубопроводы.
Лекция №8
Слайд 3Рекомендуемая литература:
Общесоюзные нормы технологического проектирования на магистральные трубопроводы. Часть 1.
Газопроводы. ОНТП 51-1-85. – М.: ВНИИЭгазпром, 1985. – 220 с.
Альбом
характеристик центробежных нагнетателей природного газа. – М.: ВНИИГАЗ, 1985. – 87 с.
Гольянов А.И., Муфтахов Е.М. Технологический расчет газопровода. Часть 1. Уфа . – 1995 . – 35 с.
Слайд 4ВЫБОР РАБОЧЕГО ДАВЛЕНИЯ, ОПРЕДЕЛЕНИЕ КОЛИЧЕСТВА КОМПРЕССОРНЫХ СТАНЦИЙ И РАССТОЯНИЯ МЕЖДУ
НИМИ
Слайд 5Расчет выполняется в соответствии с требованиями норм технологического проектирования [1].
Для
определения числа компрессорных станций необходимо уточнить рабочее давление в газопроводе
(55 или 76 кг/см2) и давление на входе компрессорной станции. Выбранные давления должны соответствовать нормативным давлениям на входе и выходе центробежных нагнетателей в соответствии с их характеристиками [2].
Слайд 61.1. Выбор рабочего давления.
Современные газопроводы работают с рабочим (избыточным) давлением
75 кг/см2.
При этом абсолютное давление на нагнетании РНАГ центробежного нагнетателя
(ЦН) не должно превышать 76 кг/см2 (76 · 0,0981 = 7,456 МПа).
1.2. Выбор давления на входе в компрессорный цех.
Анализ характеристик ЦН [2] показывает, что давление на всасывании ЦН лежат в пределах 50,7…52,8 кг/см2.
Принимаем РВС = 52 · 0,0981 = 5,10 МПа.
Слайд 71.3. Расчет характеристик транспортируемого газа.
Плотность газа при стандартных условиях (20°С
и 0,101325 МПа) определяется по формуле аддитивности
(1)
где а1,…, аn –
доля каждого компонента в смеси для данного состава газа;
ρ1,…, ρn – плотность компонента при стандартных условиях.
Слайд 9Псевдокритические температура и давление:
(4)
(5)
Относительная плотность газа:
(6)
Слайд 101.4. Суточная производительность газопровода, млн м3/сут
(7)
где kН – оценочный коэффициент
пропускной способности газопровода, который ориентировочно можно принять kН = 0,88.
Слайд 111.5. Определение расстояния между компрессорными станциями
Воспользуемся формулой расхода
(8)
откуда
(9)
где D –
внутренний диаметр газопровода, м;
PН и РК – соответственно абсолютные значения
давления в начале и в конце участка газопровода, МПа;
λ – коэффициент гидравлического сопротивления;
ZСР – средний по длине коэффициент сжимаемости газа ZСР=f(PСР, TСР);
Δ – относительная плотность газа.
Слайд 12Условный диаметр газопровода в зависимости от принятого рабочего давления ориентировочно
можно определить по таблице 1.
Таблица 1 – Ориентировочные значения диаметра
газопровода
Слайд 13Для расчета расстояния между КС можно принять ориентировочное значение средней
температуры, например
(10)
Температуру газа на входе в линейный участок TН можно
принять равной ТН = 303 К.
Слайд 14Давление в начале участка газопровода определяется по формуле [1]
(11)
где δPВЫХ
– потери давления в трубопроводе между компрессорным цехом и узлом
подключения к линейной части магистрального газопровода (без учета потерь давления в системе охлаждения транспортируемого газа);
δPОХЛ – потери давления в системе охлаждения газа, включая его обвязку. Для АВО следует принимать δPОХЛ = 0,0588 МПа. При отсутствии охлаждения газа δPОХЛ = 0. Потери давления могут быть приняты по таблице 2.
Слайд 15Таблица 2 – Потери давления газа на КС
Слайд 16Давление в конце участка газопровода
(12)
где δPВС – потери давления газа
на входе КС с учетом потерь давления в подводящих трубопроводах
и на узле очистки газа (принимается по таблице 2).
Слайд 17Коэффициент гидравлического сопротивления λ определяется по формуле
(13)
где Е – коэффициент
гидравлической эффективности, принимается по результатам расчетов диспетчерской службы в соответствии
с отраслевой методикой. При отсутствии этих данных коэффициент гидравлической эффективности принимается равным 0,95, если на газопроводе имеются устройства для периодической очистки внутренней полости трубопроводов, а при отсутствии указанных устройств принимается равным 0,92.
Слайд 18Коэффициент сопротивления трению для всех режимов течения газа в газопроводе
определяется по формуле
(14)
где kЭ – эквивалентная шероховатость труб: для монолитных
труб без внутреннего антикоррозионного покрытия принимается равной 0,03 мм;
Re – число Рейнольдса, которое определяется по формуле
(15)
где пропускная способность участка газопровода Q принимается в млн м3/сут, внутренний диаметр D – в м, коэффициент динамической вязкости μ – в Па · с.
Слайд 19В первом приближении можно принять квадратичный режим течения газа и
(16)
Коэффициент
сжимаемости газа определяется по формуле
(17)
где значения приведенных давления и температуры
определяются как
(18) (19)
(20)
Слайд 20Среднее давление в газопроводе
(21)
После определения расстояния между КС по формуле
(10) определяем число компрессорных станций
(22)
После округления найденного числа КС до
nКС (в большую сторону), уточняем значения расстояния между КС, км
(23)
Слайд 21УТОЧНЕННЫЙ ТЕПЛОВОЙ И ГИДРАВЛИЧЕСКИЙ РАСЧЕТ УЧАСТКА ГАЗОПРОВОДА МЕЖДУ ДВУМЯ КОМПРЕССОРНЫМИ
СТАНЦИЯМИ
Слайд 22Уточненный тепловой и гидравлический расчет участка газопровода между двумя компрессорными
станциями производится в соответствии с порядком расчета, изложенным в примере
1.11 методических указаний [3].
Абсолютное давление в конце участка газопровода определяем из формулы расхода (9)
(24)
В этом уравнении величина λ определяется с учетом коэффициента динамической вязкости μ при средних значениях температуры и давления.
Слайд 23Порядок дальнейшего расчета будет следующий:
1. Принимаем в качестве первого приближения
значения
λ = 0,01 ; ZСР = 0,9 .
Значение ТСР определим
по формуле (10).
2. По формуле (24) определяем в первом приближении значение РК .
3. Определяется среднее давление по формуле (21).
4. По формулам (18) и (19) с учетом средних значений давления и температуры определяем средние приведенные давление и температуру.
Для расчета конечного давления во втором приближении вычисляются уточненные значения ТСР, λ и ZСР. Для этого при определении ТСР будем использовать величины средней удельной теплоемкости СР, эффекта Джоуля-Томсона Di и коэффициента а, вычисленные для РСР и ТСР первого приближения.
Слайд 245. Удельная теплоемкость газа (кДж/(кг · К)) определяется по формуле
(25)
6.
Коэффициент Джоуля – Томсона, К/МПа
(26)
7. Рассчитываем коэффициент a по формуле
(27)
где
КСР – средний на участке общий коэффициент теплопередачи от газа в окружающую среду, Вт/(м2 · К).
Слайд 25Средняя температура по формуле
(28)
8. Коэффициент сжимаемости ZСР по формуле (17).
9.
Коэффициент динамической вязкости по формуле
(29)
10. Число Рейнольдса по формуле (15).
11.
Коэффициент сопротивления трению по формуле (14) и коэффициент гидравлического сопротивления по формуле (13).
12. Определяем конечное давление во втором приближении по формуле (24).
Слайд 2613. Если полученный результат отличается от предыдущего приближения более чем
на 1 %, имеет смысл уточнить расчеты, выполняя третье приближение,
начиная с п.3. Если результат удовлетворяет требованиям точности расчетов, переходим к следующему пункту.
14. Уточняется среднее давление по формуле (21).
15. Определяется конечная температура газа
(30)
На этом тепловой и гидравлический расчет участка газопровода заканчивается.
Слайд 27ВЫБОР ТИПА ГПА И РАСЧЕТ РЕЖИМА РАБОТЫ КС
Слайд 28Для расчетов режимов работы КС применяются характеристики ЦН, представляющие зависимость
степени повышения давления ε, политропического к. п. д. ηПОЛ и
приведенной относительной внутренней мощности
(31)
от приведенной объемной производительности
(32)
Слайд 29при различных значениях приведенных относительных оборотах
(33)
где ρВС, zВС, TВС, QВС
– соответственно плотность газа, коэффициент сжимаемости, температура газа и объемная
производительность ЦН, приведенные к условиям всасывания;
R – газовая постоянная;
zПР, RПР, TПР – условия приведения, для которых построены характеристики;
Ni – внутренняя (индикаторная) мощность;
n, nН – соответственно рабочая частота вращения вала ЦН и номинальная частота вращения.
Слайд 30Одним из универсальных видов характеристик ЦН является приведенная характеристика ВНИИГаза.
Слайд 31Порядок определения рабочих параметров следующий:
1. По известному составу газа, температуре
и давлению на входе в ЦН, используя формулы (17)…(20), определяется
коэффициент сжимаемости zВС;
2. Определяется плотность газа ρВС и производительность нагнетателя при условиях всасывания
(34) (35)
где QКС = Q – производительность КС при стандартных условиях; mН – число параллельно работающих ЦН (групп ЦН), определяемое из соотношения
(36)
Qном – номинальная производительность ЦН при стандартных условиях. Значение mН округляется вверх до целого числа.
Слайд 323. Задаваясь несколькими (не менее трех) значениями оборотов ротора в
диапазоне возможных частот вращения ГПА, определяются QПР и [n/ nН]ПР.
Полученные точки наносятся на характеристику и соединяются линией (плавная кривая abc на рисунке 1).
4. Определяется требуемая степень повышения давления
(37)
где РНАГ – номинальное давление на выходе ЦН; РВС = РК – δРВС – уточненное давление на входе (всасывании) ЦН.
Слайд 33Проведя горизонтальную линию из ε до кривой abc найдем точку
пересечения. Восстанавливая перпендикуляр до пересечения с горизонтальной осью, находим QПР.
Аналогично определяются ηПОЛ и [Ni / ρВС]ПР. Значение QПР должно удовлетворять условию QПР ≥ QПРmin, где QПРmin – приведенная объемная производительность на границе зоны помпажа (расход, соответствующий левой границе характеристик ЦН).
Слайд 345. Определяется внутренняя мощность, потребляемая ЦН
(38)
6. Определяется мощность на муфте
привода
(39)
где NМЕХ – механические потери мощности в редукторе и подшипниках
ЦН при номинальной загрузке (таблица 3).
Слайд 35Таблица 3 – Механические потери мощности для некоторых типов ЦН
Слайд 367. Вычисляется располагаемая мощность ГТУ
(40)
где NeН – номинальная мощность ГТУ;
kН – коэффициент технического состояния по мощности; kОБЛ – коэффициент,
учитывающий влияние системы противообледенения (при отключенной системе kОБЛ = 1); kУ – коэффициент, учитывающий влияние системы утилизации тепла; kt – коэффициент, учитывающий влияние атмосферного воздуха на мощность ГТУ; TВОЗД, TВОЗДН – соответственно фактическая и номинальная температура воздуха, К.
Значения NeН, kН, kОБЛ, kУ, kt, TВОЗДН принимаются по справочным данным ГТУ.
Слайд 378. Производится сравнение Ne и NeP. Должно выполняться условие Ne
≤ NeP. При невыполнении этого условия следует увеличить число mН
и повторить расчет начиная с пункта 2.
9. Определяется температура газа на выходе ЦН
(41)
где k – показатель адиабаты природного газа, k = 1,31.