2. РАСЧЕТ ПАРАМЕТРОВ РАБОТЫ МАГИСТРАЛЬНЫХ ГАЗОПРОВОДОВ
Пропускная способность и режим работы магистрального газопровода (МГ) определяется совместной работой КС и линейных участков, его составляющих. При этом режимы работы отдельных КС и участков, в связи с различием их геометрических размеров, давления и температуры газа, значительно отличаются, что диктует необходимость поочередного расчета всех элементов системы. Выходные параметры одного элемента являются входными параметрами следующего за ним. Таким образом, поочередно следуют расчеты работы участков и компрессорных станций. На каждом этапе ведется проверка соответствия полученных параметров условиям нормальной работы газопровода и оптимальности режима.
По аналогичной схеме будет работать математическая модель МГ при расчетах на ЭВМ. Для ее реализации требуется формализовать работу основных элементов системы: участок, газоперекачивающий агрегат (ГПА), пылеуловитель (ПУ), аппарат воздушного охлаждения (АВО).
Рассмотрим поочередно основные уравнения, позволяющие оценить изменения параметров газа при его перемещении по элементам МГ, и на их основе составим алгоритмы решения важнейших задач эксплуатации газопровода.
2.1. Физические свойства газа
Широкое использование ЭВМ диктует необходимость аналитического определения физических свойств газа. При гидравлических и тепловых расчетах МГ используются следующие физические величины: коэффициент сжимаемости, динамическая вязкость, удельная теплоемкость и коэффициент Джоуля - Томсона. Базовым параметром является относительная плотность газа , или плотность газа при стандартных условиях , между которыми существует следующая связь:
. (2.1)
Коэффициент сжимаемости (z) и динамическая вязкость газа () определяются через приведенные значения давления и температуры:
, (2.2)
где и - приведенные давление и температура;
и - давление, при котором определяются свойства и критическое давление газа;
и - температура, при которой определяются свойства и критическая температура газа:
(2.3)
. (2.4)
В настоящее время для определения коэффициента сжимаемости и динамической вязкости газа рекомендуется использовать следующие зависимости [2]:
, (2.5)
где ;
(2.6)
,
где - динамическая вязкость газа, Па c.
Удельная теплоемкость cp (кДж/(кгК)) и коэффициент Джоуля - Томсона Di (К/МПа) газа определяются из уравнений (2.7) и (2.8):
, (2.7)
. (2.8)
При решении задач следует постоянно следить за соответствием условий, при которых определяются физические свойства газа реальным, условиям рассчитываемого участка газопровода.
2.2. Расчет давления
Давление является основным параметром, по которому контролируется режим работы трубопроводов.
Газ поступает на КС с давлением и температурой в конце подводящего трубопровода (головная КС) или предшествующего участка (промежуточные КС) P2 и Т2. На выходе станции (в начале следующего участка) давление будет P1. Степень сжатия станции при этом составит
. (2.9)
Учитывая потери во входном и выходном коллекторах, степень сжатия нагнетателей должна быть более высокой:
, (2.10)
где - потери давления во входном и выходном коллекторах КС (приложение 3);
- потери давления в аппаратах воздушного охлаждения, = 0,0588 МПа.
Возможность реализации требуемой степени сжатия определяется располагаемой мощностью привода нагнетателя :
, (2.11)
где - потребляемая мощность двигателя и внутренняя мощность нагнетателя;
- механический кпд нагнетателя (приложение 4).
Для электроприводных ГПА .
Внутренняя мощность нагнетателя определяется с помощью приведенных характеристик:
, (2.12)
где - приведенная мощность нагнетателя,
- фактическая и номинальная частота вращения ротора нагнетателя;
- плотность газа при условиях входа в нагнетатель, кг/м3.
, (2.13)
где - давление (абсолютное) газа на входе нагнетателя и при стандартных условиях, Па;
- температура газа на входе нагнетателя и при стандартных условиях, Т;
- коэффициент сжимаемости газа при условиях входа в нагнетатель;
R - газовая постоянная, Дж/(кгК):
. (2.14)
Допускается определение внутренней мощности из (2.15)
, (2.15)
где - внутренняя мощность, кВт;
Q - производительность нагнетателя, млн.м3/сут;
- политропический кпд нагнетателя;
- показатель адиабаты сжатия, =1,31.
Располагаемая мощность ГТУ и электродвигателя зависит от условий их работы и для ГТУ определяется по формуле (4. 7).
Используя (2.11) и (2.15), можно определить максимальную производительность нагнетателя, а соответственно и КС, при требуемой степени сжатия и располагаемой мощности ГПА или максимальную степень сжатия при заданной производительности. Полученный таким образом результат не всегда может быть реализован при использовании конкретного типа нагнетателя с заданной проточной частью. Реализуемые значения степени сжатия определяются только по приведенным характеристикам ЦН.
Приведенная характеристика нагнетателя представляет собой графическое изображение трех функций:
.
Аналитически эти функции можно аппроксимировать полиномами
(2.16)
где QЛР - приведенная производительность, м3/мин:
, (2.17)
где - номинальная и фактическая частота вращения рабочего колеса нагнетателя;
QВ - производительность при условиях входа в нагнетатель, м3/мин
, (2.18)
где - производительность нагнетателя при стандартных условиях, млн. м3/сут.
Для определения коэффициентов a, b и с достаточно использовать значения при производительностях нагнетателя:
ПРmin- минимальная производительность;
ПРmax - максимальная производительность;
ПРmax1 - производительность при максимальном значении политропического кпд;
ПРmax2 - производительность при максимальном значении приведенной мощности.
Зависимость описывается для случая
, (2.19)
где - приведенные параметры нагнетателя (прило -
жение 6);
- приведенная частота вращения рабочего колеса ЦН.
Развиваемая ЦН степень сжатия при любой частоте вращения может быть пересчитана по уравнению
, (2.20)
где - степень сжатия ЦН при приведенной частоте вращения, равной 1,00.
Уравнение (2.20) позволяет определить частоту вращения рабочего колеса, при которой ЦН создает требуемую степень сжатия:
. (2.21)
При этом должны соблюдаться условия (2.11), (2.22), (2.23) и
, (2.22)
. (2.23)
Давление на выходе КС не должно превышать рабочего давления МГ :
. (2.24)
После компримирования газ, пройдя через АВО (на МГ большого диаметра), поступает в участок с параметрами P1 и T1.
Основным расчетным уравнением для участка МГ является уравнение пропускной способности. В общем случае пропускная способность участка зависит от его длины, внутреннего диаметра труб, перепада давления, физических свойств газа и рельефа трассы. Учитывая, что плотность газа мала, влиянием рельефа чаще всего можно пренебречь. Согласно ОНТП газопровод рассчитывается как горизонтальный при разности геодезических отметок менее 100 м. В этом случае уравнение пропускной способности используется в следующем виде:
, (2.25)
где - пропускная способность участка, млн. м3/сут;
- давление в начале и в конце участка, МПа;
- эквивалентный диаметр труб, м;
- коэффициент сжимаемости газа при среднем значении давления и температуры в участке;
- средняя температура газа в участке, К;
- длина участка, км;
- расчетное значение коэффициента гидравлического сопротивления.
Для определения теоретического значения при турбулентном течении газа рекомендуется формула ВНИИгаза:
, (2.26)
где - эквивалентная шероховатость труб;
- число Рейнольдса:
. (2.27)
При полной загрузке МГ чаще всего работают в квадратичной зоне, и в этом случае (2.26) при рекомендуемом значении k=0,03 мм принимает вид
, (2.28)
где - внутренний диаметр труб, мм.
Работа газопровода в квадратичной зоне наблюдается при производительностях, превышающих переходное значение :
. (2.29)
При определении гидравлического сопротивления учитывают возможность засорения трубопровода в процессе эксплуатации и наличие местных сопротивлений:
, (2.30)
где - коэффициент эффективности работы участка, принимаемый при регулярной очистке МГ равным 0, 95.
В связи с тем, что давление по длине участка меняется нелинейно, среднее давление P определяется как среднегеометрическое:
. (2.31)
Практически участок между КС представляет собой либо несколько параллельных трубопроводов, каждый из которых состоит из нескольких подучастков с различным внутренним диаметром, либо несколько подучастков, отличающихся числом ниток или внутренним диаметром. Расчет таких сложных трубопроводов заменяют расчетом простого, используя понятия эквивалентного диаметра или коэффициента расхода.
Эквивалентным диаметром называется диаметр простого трубопровода, имеющего пропускную способность, равную пропускной способности реального трубопровода при прочих равных условиях.
Коэффициентом расхода называют отношение пропускной способности реального трубопровода к пропускной способности эталонного трубопровода с произвольно выбранным эталонным диаметром при прочих равных условиях:
. (2.32)
Для случая простого трубопровода
, (2.33)
где и - диаметр и коэффициент гидравлического сопротивления простого трубопровода;
и - диаметр и коэффициент гидравлического сопротивления эталонного трубопровода.
При квадратичном режиме течения газа
. (2.34)
При параллельном соединении простых трубопроводов
, (2.35)
. (2.36)
При последовательном соединении трубопроводов
, (2.37)
, (2.38)
где - длина участка и подучастков.
Для сложного газопровода с последовательным и параллельным соединением участков коэффициент расхода или эквивалентный диаметр определяются последовательным использованием формул (2.35) и (2.37) или (2.36) и (2.38).
При расчете МГ с применением коэффициента расхода (2.25) имеет следующий вид:
, (2.39)
где - диаметр и коэффициент гидравлического сопротивления эталонного трубопровода.
Уравнения (2.25) и (2.39) позволяют определить давление газа в конце участка P2 (на входе в следующую КС).
2.3. Расчет температуры
Газ, поступающий на КС с температурой , при компримировании нагревается до температуры :
. (2.40)
На МГ малого диаметра температура на выходе станции равняется . Если станция оборудована АВО, то в этом случае
, (2.41)
где - теоретический теплосъем с одного АВО при двух работающих вентиляторах, Вт;
- коэффициенты тепловой эффективности АВО при,2 1и 0 работающих вентиляторах;
- количество АВО, работающих с 2,1 и 0 вентиляторов;
- массовый расход газа через все АВО, кг/с;
- теплоемкость газа при условиях АВО, Дж/(кг град) при
. (2.42)
Теплосъем Q0 удобно определять по номограммам теплового расчета АВО [4]. При расчетах на ЭВМ охлаждение газа можно представить следующим уравнением:
, (2.43)
где - постоянная для данного АВО величина, определяемая по данным эксплуатации или по номограмме;
- массовый расход одного АВО, кг/с.
В этом случае и .
Значения коэффициентов и также определяются по результатам эксплуатации АВО. В первом приближении можно принять =0,55-0,60, =0,18-0,20.
Количество работающих вентиляторов для реализации заданной температуры на выходе КС определяется из (2.41). Принимая во внимание, что при регулировании температуры сначала отключают поочередно по одному вентилятору на всех АВО и только после этого начинают отключение вторых, в сумме уравнения (2.41) никогда не будет больше двух слагаемых. Для определения схемы работы удобно воспользоваться величиной среднего коэффициента эффективности:
, 2.44)
где - количество работающих на КС АВО.
В зависимости от величины возможны следующие варианты:
;
;
;
.
Схема работы АВО определяется из условия минимума затрат электроэнергии для обеспечения оптимальной температуры газа за КС. Не рекомендуется принимать температуру на выходе КС более 500 С.
Газ в участке охлаждается вследствие теплообмена с окружающей средой и его расширения при снижении давления. С учетом этих факторов температуру в конце участка можно определить из уравнения ВНИИгаза:
, (2.45)
где - коэффициент, характеризующий интенсивность снижения температуры газа по длине участка:
; (2.46)
- полный коэффициент теплопередачи, =(1,5-2,0)Вт/(м2град);
- температура грунта на глубине заложения трубопровода.
Пренебрегая влиянием дроссельного эффекта получим уравнение Шухова В.Г.:
. (2.47)
В соответствии с (2.46) температура газа стремится к температуре окружающей среды. С учетом дроссельного эффекта температура газа в конце участка будет меньше температуры окружающей среды. При температуре грунта близкой к 00 С температура газа может быть отрицательной, что вызовет промерзание грунта вокруг труб и приведет к появлению дополнительных деформаций трубопровода. Рекомендуется ограничивать температуру газа в конце участка Т2 = 271-273 К, что приводит к ограничению температуры газа на выходе КС. Если МГ проложен в многолетнемерзлых грунтах, то температура газа в конце участка может быть равной температуре грунта. Минимальная температура на выходе КС определяется из (2.45).
Так как температура газа по длине участка меняется экспоненциально, то средняя температура определяется как среднегеометрическая:
. (2.48)
Приобететение доступа к файлам
ВНИМАНИЕ: Данная информация получена путем сканирования, цифровой обработки физических носителей или обмена с неравнодушными пользователями. Она не имеет отметок грифа секретности и тайны, если вы считаете, что эта информация нарушает Ваши авторские или другие права. Незамедлительно сообщите администратору для удаления ее из портала.