3.1 Характеристики топливного газа
По заданному составу топлива определяются внутренняя теплота сгорания, молекулярная масса и теплоемкость на входе в КС.
В качестве примера рассмотрим природный газ, имеющий следующий состав [18]:
метан (СН4) - 98%; этан (С2Н6) - 0,45%; пропан (С3Н8) - 0,1%; бутан (С4Н10) - 0,02%; азот (N2) - 0,63%; кислород (О2) - 0,78%; диоксид углерода (СО2) - 0,02%.
Плотность топливного газа сПГ = 0,7231 кг/м3..
Температура топливного газа tПГ = 16 0С.
Молекулярная масса топливного газа (смеси газов) для данного состава природного газа, кг/кмоль
µПГ = (3.1)
Примечания.
1. Молекулярные массы газов, входящих в с10остав газообразного топлива, см. в табл. П 1.3. приложения 1.
2. Учитывая, что в смеси газов преобладает метан (98%), можно принять µПГ =16,04 (погрешность ).
Низшая теплота сгорания сухой массы топлива (1 м3 при НФУ (нормальных физических условиях) в данном случае, кДж/нм3,
(3.2)
Соотношение между рабочей массой топлива и сухой [1] пересчитывается при помощи множителя
, (3.3)
где WР - содержание влаги в составе топлива.
В составе топлива, приведенном выше, влага отсутствует (WР = 0), поэтому .
Примечания.
1. В тепловых расчетах необходимо учитывать, в каких единицах, объемных (кДж/нм3) или массовых (кДж/кг), применяется теплота сгорания топлива. Удобнее в формулах применять в кДж/кг. Это упрощает вид формулы, однако при этом целесообразно проверять размерности величин, входящих в формулу и размерность получаемого результата вычислений.
2. Перевод теплоты сгорания рабочей массы топлива из (кДж/нм3) в (кДж/кг) выполняется по соотношению
.
= R
. (3.4)
Здесь - давление в камере сгорания (принято равным давлению воздуха в камере сгорания), бар; - газовая постоянная природного (топливного) газа, кДж/(кг · К); - плотность природного газа, кг/м3.
Удельная массовая изобарная теплоемкость топливного газа, отнесенная к 1 м3 сухого газа, согласно [1] определяется по формуле, кДж/(кг·К),
, (3.5)
где - удельная массовая изобарная теплоемкость компонентов топливного газа, кДж/(кг · град);
, … - объемное содержание газа в газовой смеси.
Удельную массовую изобарную теплоемкость топливного газа можно оценить по приближенной формуле (расчет по метану)
.
Энтальпия топливного газа, кДж/кг,
. (3.6)
3.2 Продукты сгорания топлива
Определяются продукты сгорания топлива, теоретически необходимый объем воздуха, состав продуктов сгорания при заданном избытке воздуха. Окончательный избыток воздуха перед 1-й ступенью ГТ находится после расчета температуры газов за последней ступенью ГТ. Избыток воздуха в КС принимается в пределах 1,2 ч 1,6 с таким условием, чтобы температура ядра факела не превышала 1600 оС.
Теоретический объем воздуха V0, необходимый для сгорания 1 м3 газа, согласно [1] определяется по формуле, нм3/нм3,
(3.7)
Количество продуктов сгорания:
а) объем сухих трехатомных газов (), м3/м3,
), (3.8)
б) объем азота, м3/м3,
, (3.9)
в) объем водяных паров, м3/м3,
(3.10)
где объемная доля горючих газов в топливе, %, - влажность газа согласно нормам расчета [1].
Коэффициент избытка воздуха в камере сгорания [4]
(3.11)
Избыток воздуха в КС, м3/м3,
= () (3/12)
Суммарный объем продуктов полного сгорания, м3/м3,
(3.13)
Примечания:
1. Обычно расчет ведут в объемных долях:
;
;
;
;
2. Состав продуктов сгорания задают или представляют по результатам расчетов в процентах:
3. Размерность (м3/м3) читается как объем продуктов сгорания (м3), получающийся при сжигании 1 м3 топливного (природного) газа.
Избыток воздуха за НТД
. (3.14)
Здесь - коэффициент, учитывающий содержание кислорода в уходящих газах, принят по данным эксплуатации (по показаниям приборов) ГТД - 110.
3.3 Давление и температура воздуха в камере сгорания
При попадании воздуха после компрессора в камеру сгорания его давление () уменьшается вследствие дросселирования в жаровых трубах (). Следовательно, давление воздуха в камере сгорания (см. рис. 9)
, (3.15)
где бар.
Примечание. В зависимости от степени открытия ВНА (%) расход воздуха через компрессор будет изменяться, следовательно в общем случае . Учитывая, что положение ВНА может изменяться в пределах от 00 до - 350, согласно инструкции по эксплуатации [11]. Можно считать величину постоянной.
Относительное давление с учетом дросселирования воздуха в жаровых трубах
/ . (3.16)
Стандартное отношение относительных давлений
. (3.17)
Стандартная энтропия воздуха в камере сгорания
= . (3.18)
Действительная энтропия воздуха в камере сгорания
= / .
, кДж/(кг•К) (3.19)
3.4 Расчет теоретической температуры горения
Суммарная энтальпия продуктов сгорания (hТ), соответствующая теоретической температуре горения (tТ), равна теплу, подведенному в камеру сгорания (QKС) и складывается из тепла реакции горения топлива (QРН), тепла, внесенного с воздухом (hB) и физического тепла топливного газа (hПГ). Следовательно, теоретической температуре горения отвечает равенство:
hТ(tТ) = QKС = QРН + hB(t2 В) + hПГ(tПГ). (3.20)
Учитывая, что
hТ(tТ) = cР Г(tТ) • tТ, (3.21)
зависимость для теоретической температуры горения имеет вид
tТ = t2 В + tТ = t2 В + hТ(tТ) / cР Г(tТ). (3.22)
Расчет этой температуры проводится с учетом зависимости теплоемкости продуктов сгорания от температуры. Учитывая, что процесс сгорания протекает при постоянном давлении, расчеты ведутся по изобарной теплоемкости (массовой или молярной). Учитывая тот факт, что расчеты теоретической и действительной температур горения можно проводить через молярные или массовые величины, алгоритм вычислений выглядит следующим образом.
1. Молекулярная масса газообразного продукта сгорания топлива (смеси газов), кг/кмоль:
Г = ?(i • ri) = мH2O • rH2O + м RO2 • rRO2 + м N2 • rN2+ В • rB =
= 18,016 • rH2O + 44,01 • rRO2 + 28,15 • rN2+ 28,97 • rВ. (3.23)
2. Газовая постоянная газообразного продукта сгорания топливной смеси определится из соотношения, кДж/(кг•К):
RСМ = мR /мСМ > RГ = 8,3145 /мГ. (3.24)
3. Вычисляется молярная энтальпия газообразного продукта сгорания топлива, соответствующая теоретической температуре горения, кДж/кмоль:
HТ(tТ) = QРН • VН + HB(t2 В) + HПГ(tПГ) = QРН • VН + h2В КС • В + hПГ • ПГ, (3.25)
где: Q - кДж/м3; H - кДж/кмоль; VН = 22,414 м3/кмоль - объем 1 кмоля любого газа при НФУ; h - кДж/кг; - кг/кмоль.
4. Задаемся рядом возможных теоретических температур горения (например, tТ = 1000; 1050; 1100; 1150; 1200; 1250 OC).
5. Для каждого значения температуры по таблицам [2] вычисляются удельные молярные энтальпии (H) компонентов газообразного продукта сгорания топлива, кДж/кмоль:
HH2O, HRO2, HN2, HВ = f(tТ). (3.26)
6. По правилу смешения вычисляется энтальпия газообразного продукта сгорания топлива (H Г), кДж/кмоль:
H Г = rH2O • HH2O + rRO2 • HRO2 + rN2 • HN2 + rB • HВ. (3.27)
7. По результатам вычислений путем интерполяции, учитывая, что в молярных величинах HГ = HТ, находим приращение температуры горения, OC (рис. 10)
tТ = f(HГ = HТ). (3.28)
8. Вычисляем теоретическую температуру горения, OC
tТ = t2 В + tТ. (3.29)
9. Общий КПД камеры сгорания можно выразить в виде произведения:
зКС = зТКС • зАКС. (3.30)
Здесь: зТКС - тепловой КПД КС; зАКС - аэродинамический КПД КС.
Тепловой КПД камеры сгорания, учитывающий все тепловые потери [4], определяется по формуле:
зТКС = 1 - (QНС + QОХЛ) / (QРН)), (3.31)
где: QНС - потери теплоты от неполного сгорания топлива (химический и физический недожог), кДж/м3, (кДж/кг). У современных КС эти потери не должны превышать 15 % общего расхода теплоты при работе во всем диапазоне рабочих нагрузок и 13 % при работе на расчетной нагрузке;
QОХЛ - потери за счет отдачи теплоты в окружающее пространство нагретой поверхностью КС и примыкающих к ней трубопроводов. Эти потери обычно бывают не более 0,5 % расхода теплоты, кДж/м3, (кДж/кг).
В существующих камерах сгорания тепловой КПД при работе на расчетном режиме:
зТКС = 0,97 … 0,98. (3.32)
Полные потери давления в камере сгорания складывается из следующих составляющих:
а) гидравлические (аэродинамические) потери, которые возникают без подвода теплоты в камере в результате потерь на трение при прохождении газового потока и наличия местных сопротивлений от воздухонаправляющих рёбер, завихрителей и т.д. Эти потери определяются при холодной продувке камеры:
б) дополнительные потери давления, вызванные нагревом газа при сгорании топлива в камере. Плотность газа в этом случае уменьшается, а скорость газового потока увеличивается (за счет увеличения объёмного расхода газа).
Потери давления в камере сгорания снижают КПД ГТУ. Это можно учесть с помощью аэродинамического КПД камеры сгорания зАКС, который обычно составляет:
зАКС = 0,98 … 0,99. (3.33)
Общий КПД современных камер сгорания согласно [5, 6]:
зКС = зТКС • зАКС = 0,95 … 0,98. (3.34)
Расчет действительной температуры горения базируется на алгоритме расчета теоретической температуры горения. Последовательность вычислений аналогична вычислению теоретической температуры горения. Расширенный алгоритм вычислений приводится в молярных величинах.
1. Вычисляется молярная энтальпия газообразного продукта сгорания топлива, которой соответствует действительная температура горения (действительная температура газов перед первой ступенью ГТ: t3 = tД):
HД(tД) = H3(t3) = QKС • зТП = [QРН • VН + HB(t2 В) + HПГ(tПГ)] • зТП, (3.35)
где: Q - кДж/м3; H - кДж/кмоль; VН = 22,414 м3/кмоль - объем 1 кмоля любого газа при НФУ; зТП = 0,97 … 0,98 - КПД теплового процесса горения, учитывающий тепловые потери КС в окружающую среду, так как корпус ГТД не имеет тепловой изоляции и применяется тепловое укрытие с применением вентиляторов отсоса воздуха из него, вследствие чего имеют место увеличенные потери в окружающую среду (рис. 11).
2. Задаемся рядом возможных теоретических температур горения (например, t3 = 1000; 1050; 1100; 1150; 1200; 1250 OC).
3. Для каждого значения температуры по таблицам [2] вычисляются удельные молярные энтальпии (H) компонентов газообразного продукта сгорания топлива, кДж/(кмоль):
HH2O, HRO2, HN2, H3В = f(t3). (3.36)
4. По правилу смешения вычисляется энтальпия газообразного продукта сгорания топлива (H Г), кДж/(кмоль•К):
H3 = rH2O • HH2O + rRO2 • HRO2 + rN2 • HN2 + rB • H3В. (3.37)
5. По результатам вычислений путем интерполяции, учитывая, что H3 = HД, находим приращение температуры горения, OC
tД = f(H3 = HД). (3.38)
6. Вычисляется действительная температура горения, OC
t3 = t2 В + tД. (3.39)
Рис. 9. Процесс в h,s-диаграмме изменения давления воздуха в камере сгорания
Рис. 10. H, ?t-диаграмма продукта сгорания топлива для определения теоретической и действительной температур горения
Рис. 11. Схема теплового укрытия ГТЭ-110: К - компрессор; КС - камера сгорания; ГТ - газовая турбина; Диф. - диффузор соединения ГТ с котлом-утилизатором
4. РАСЧЕТ ПАРАМЕТРОВ ПРОЦЕССА В ГАЗОВОЙ ТУРБИНЕ
Расчет параметров процесса ГТ осуществляется с учетом влияния на характер процесса расширения газов в ГТ воздуха, подводимого из компрессора на охлаждение элементов проточной части (рис. 12).
Особенностью работы ПГУ утилизационного типа является необходимость поддержания практически постоянной температуры выхлопных газов за газовой турбиной, т.е. перед котлом-утилизатором (КУ) в диапазоне нагрузок от 100 до 70% от номинальной нагрузки [15].
Для турбин ГТД-110 производства НПО «САТУРН» значение температуры газов за ГТ определено заводом-изготовителем на уровне 517оС [3]. Контролируемая величина этой температуры определяется как среднее значение, рассчитываемое программно-техническим комплексом (ПТК) ГТУ по 20-и датчикам, размещенных на выходе отработанных газов из ГТ. Количество точек измерения температуры соответствует количеству жаровых труб, расположенных в камере сгорания.
Поддержание контролируемого значения температуры газов за ГТ согласно инструкции по эксплуатации ГТД-110, работающей на КУ [3], производится путем автоматического изменения положения входного направляющего аппарата (ВНА), установленного перед первой ступенью компрессора. ВНА служит для управления расходом воздуха, необходимого для поддержания заданной температуры за ГТ [3]. Нагрузка ГТУ в диапазоне 100 60% от номинальной мощности и температура газов на выходе из ГТ регулируется путем изменения расхода воздуха через компрессор с помощью ВНА и расхода топлива регулирующим топливным клапаном (РКТ) при примерно постоянной температуре выхлопных газов ГТ. Этот диапазон нагрузок является наиболее экономичным, так как в нем незначительно изменяется КПД ГТУ.