CC BY
54
УДК 621.311.22
Оценка величины потери теплоты с химическим недожогом в камере сжигания дополнительного топлива
Е.С. Малков
ФГБОУВПО «Ивановский государственный энергетический университет имени В.И. Ленина»,
г. Иваново, Российская Федерация E-mail: admin@tes.ispu.ru
Авторское резюме
Состояние вопроса: В настоящее время теоретически обоснована возможность применения камеры сжигания дополнительного топлива и газового подогревателя сетевой воды при их последовательной установке за газовым подогревателем конденсата в газоходе котла-утилизатора и предложен способ получения тепловой энергии для нужд теплофикации на ПГУ-КЭС. Условия для протекания процесса горения в предложенной компоновке отличаются от условий в схемах, применяемых в настоящее время в энергетике. Поэтому для комплексной оценки эффективности такого способа получения тепловой энергии необходимо оценить величину потери теплоты с химическим недожогом.
Материалы и методы: Решение поставленной задачи осуществлено двумя способами - с применением нормативного метода расчета котельных агрегатов и с использованием программного комплекса Flow Vision (ООО «ТЕСИС») для проведения многовариантных расчетов характеристик газового потока в котле-утилизаторе. Результаты: На основе многовариантных расчетов в программном комплексе Flow Vision получена зависимость величины потери теплоты с химическим недожогом от характеристик газового потока на входе в камеру сжигания дополнительного топлива в диапазоне температур 90-150оС и объемной концентрации кислорода 12-17 %. Выводы: Разные способы оценки величины потери теплоты с химическим недожогом позволили получить практически одинаковые результаты. Результаты работы необходимы для комплексного анализа эффективности применения камеры сжигания дополнительного топлива для нужд теплофикации в условиях расположения ее за газовым подогревателем конденсата котла-утилизатора и экономичности работы парогазового энергоблока.
Ключевые слова: котел-утилизатор, камера сжигания дополнительного топлива, потеря теплоты с химическим недожогом, газовый подогреватель сетевой воды.
Measuring of heat losses under incomplete combustion in supplementary combustion chambers
E.S. Malkov
Ivanovo State Power Engineering University, Ivanovo, Russian Federation E-mail: admin@tes.ispu.ru
Abstract
Background: It is currently proved theoretically that supplementary combustion chambers and network water gas heaters can be used in recovery boilers if placed in tandem after the condensate gas heater in the boiler gas duct. The paper presents a method of obtaining thermal power for district heat supply provided by combined cycle gas turbine units. The combustion conditions in the proposed arrangement differ from those in the schemes currently used in power engineering. Therefore, to give a comprehensive assessment of the efficiency of the new method of obtaining thermal power, it is necessary to calculate the heat losses under incomplete combustion.
Materials and methods: The problem was solved in two ways - using the standard method of calculating boilers and applying the Flow Vision software (OOO «TESIS») to the multivariate calculation of the recovery boiler gas flow characteristics.
Results: The multivariate calculations made in the Flow Vision software allowed the authors to obtain the dependence of incomplete combustion heat losses on the characteristics of the gas flow entering the supplementary combustion chamber within the temperature range of 90-150 °C and at the oxygen volume concentration of 12-17 %.
Conclusions: The two different methods of incomplete combustion heat loss calculation had almost identical results. The research results can be employed in comprehensive analysis of the efficiency of the supplementary combustion chamber used for district heat supply provision if the chamber is set after the condensate gas heater, and the efficiency of combined-cycle power plants.
Key words: recovery boiler, supplementary combustion chamber, incomplete combustion heat loss, network water gas heater.
Объемная концентрация кислорода О2 в утилизатором (КУ) дух составляет 100-140 °С.
выхлопных газах энергетических газовых тур- По данным Всероссийского теплотехнического
бин 13-16 %, а их температура за котлом- института (ВТИ), в тепловых схемах парогазо-
вых установок (ПГУ) эксплуатация горелочных устройств КУ возможна при О2 > 12 %, коэффициенте избытка воздуха а > 2 и температуре газов &газ > 100 оС [1]. Поэтому выхлопные газы газотурбинной установки (ГТУ) можно использовать в качестве малоактивного окислителя процесса горения. Данный потенциал потока газов предложено использовать для увеличения тепловой мощности ПГУ-КЭС за счет установки камеры сжигания дополнительного топлива (КСДТ) и газового подогревателя сетевой воды (ГПСВ) за газовым подогревателем конденсата (ГПК) КУ.
При горении топлива в потоке газов, забалластированных продуктами сгорания и обедненных кислородом, снижаются полнота сгорания и скорость распространения пламени. Поэтому для анализа изменения значения КПД КУ пкУ при сжигании дополнительного топлива необходимо оценить величину потери теплоты от химической неполноты сгорания q3.
При отсутствии опытных данных и аналитической информации данная задача была решена двумя способами - с использованием нормативного метода теплового расчета котельных агрегатов1 (далее - Нормативный метод) и программного продукта Flow Vision.
Состав топливно-газовой смеси в пред-пламенной зоне КСДТ близок к составу топ-ливно-воздушной смеси при сжигании доменного газа с характеристиками, представленными в табл. 1.
Таблица 1. Характеристики балласта доменного газа и отдельных составляющих продуктов сгорания
Наименование характеристик Значения
Концентрация в топливе:
азота М2, % 55
диоксида углерода СО2, % 12,5
Теоретические объемы:
воздуха vе, м3/м3 0,78
газов Vг, м3/м3 1,63
водяных паров ^20, м3/м3 0,07
Коэффициент избытка воздуха в зоне 1,1
горения а
Теплота сгорания онс, МДж/м3 3,78
Потеря теплоты от химической непол- 1,5
ноты сгорания д3, %
Тепловое напряжение топочного объе- 0,23
ма ду, МДж/м3
Значения теоретических объемов отдельных составляющих продуктов сгорания природного газа «Саратов - Н.Новгород» при коэффициенте избытка воздуха аГТУ = 1,0 представлены в табл. 2. Для данного топлива при высокой суммарной концентрации горючих компонентов (95,3 %) характерны высокие значения теплоты сгорания и объемов газов.
Согласно данным табл. 2, при сжигании в камере сгорания ГТУ природного газа при а = 2,8 удельный объем уходящих из ГТУ газов составляет \/ухгаз = 36,8 м3/м3.
Таблица 2. Значения теоретических объемов отдельных составляющих продуктов сгорания природного газа при а = 1,0
Наименование характеристик
Теоретические объемы: воздуха vе, м3/м3 диоксида углерода ур*о2, м3/м3 азота ум2, м3/м3 водяных паров ун20, м3/м3 Суммарный объем газов V--, м3/м3 Теплота сгорания онс, МДж/м3
Значения
9,57
1,03
7,59
2,13
10,75
36,1
Удельный расход уходящих из ГТУ газов, используемых в качестве окислителя в КСДТ,
при объеме кислорода в них ^02 = 4,41 м3/м3 и
а = 1,1 составляет
0,21аГУ0 ,
ксдт
V ух vo2
V-Т, = 21,1 м3/ м3.
Суммарная концентрация азота и диоксида углерода в топливно-окислительной смеси (в предпламенной зоне) при сжигании топлива в количестве иТОПЛ = 1 м составляет:
• для доменного газа
к, 0,01(М2+С02) + 0,79оУ ° л„, ,„„, Ы2+СО2 = -—^-^-:-г—100% = 73%;
2 2 ^0пл +агу0
• для природного газа, сжигаемого в КСДТ, при а = 1,1
0,01(М2+СР2) + 0,79УКСДТ 1(
n2+co2 =
1 + V
100% = 75%.
ксдт
1 Тепловой расчет котлоагрегатов (норм. метод) / под ред.
Н.В. Кузнецова. - М.: Энергия, 1998.
Полученные результаты свидетельствуют о том, что обе топливно-окислительные смеси перед сжиганием инертными газами забалластированы в равной мере.
Таким образом, для оценки величины потери теплоты от химической неполноты сгорания д3 при сжигании в КСДТ в первом приближении можно ориентироваться на данные при сжигании доменного газа.
Если химический недожог представлен в виде оксида углерода, трудносжигаемого компонента, то в случае сжигания доменного газа с располагаемой теплотой топлива Одг = 3,86 МДж/м3 при удельном объеме сухих
газов Усдхг = 1,71 м3/м3 его концентрация в продуктах сгорания, %, составляет
дг 0,01 дд г Ордг С0д = -—43 дгр = 0,268,
осо Ч-ух
где ОСО = 12,64 МДж/м - теплота сгорания оксида углерода.
При одинаковых значениях концентрации инертных газов в топливно-окислительной смеси полнота сгорания топлива зависит от
вероятности контакта кислорода с горючими элементами в предпламенной зоне. Поэтому концентрация оксида углерода в продуктах сгорания за КСДТ может быть определена величиной, обратно пропорциональной содержанию кислорода в предпламенной зоне.
Относительное содержание кислорода в топливно-воздушной смеси для 1 м3 топлива составляет:
• перед сжиганием доменного газа
= 02ЧГ = 0,097 м3/м3;
1 + агуо
• перед сжиганием природного газа в КСДТ
0ксдт = 15,5агц<сдт = о 9 м3/м3, 2 100(1 + агц<сдт)
где 15,5 - концентрация кислорода в газах, поступающих из ГТУ в КСДТ, %.
Тогда концентрация оксида углерода в продуктах сгорания за КСДТ, %, определяется как
СОкедт = срд.г одг
г о2сдт
= 0,175.
При удельном объеме сухих уходящих газов Ус = 17,3 м3/м3, температурах газов за КУ Эку = 98-115°С и располагаемой теплоте топлива С^'1" = 36,3 МДж/м3 потери теплоты от
химической неполноты сгорания природного газа в КСДТ, %, составляют
12,64 СОКСДТУс[у|х г100%
ц3 =—-г-сух-= 1,06%.
Српр.г -
На основании полученных результатов можно сделать вывод, что использование нормативного метода расчета котельных агрегатов в сравнении с условиями сжигания доменного газа позволяет получить величину потери теплоты от химической неполноты сгорания при сжигании в КСДТ природного газа примерно равной д3 = 1,0 %.
Ввиду активного развития компьютерной техники и высокой стоимости проведения натурных экспериментов, основным подходом к решению большинства технических задач стало математическое моделирование рабочих процессов, позволяющее глубоко и полно исследовать влияние конструктивных и режимных факторов на основные показатели работы оборудования.
Для решения задач газовой динамики в настоящее время широко применяются численные методы, которые сводятся к тому, что область решения задачи, называемая моделью, разбивается на систему подобластей (или узлов), в каждой из которых набор газодинамических параметров заменяется постоянными функциями. При этом газодинамическая и физико-химическая задачи не разделимы. В
подавляющем большинстве случаев газодинамические параметры, характеризующие динамику горения, зависят от внутренних свойств пламени.
В КСДТ, в основном, применяются диф-фузионно-стабилизаторные горелки, сжигание топлива в которых осуществляется в турбулентном следе. В таких горелках используются уголковые стабилизаторы, расположенные вертикально в одной плоскости. Природный газ вводится в зону рециркуляции системой струй через специальный трубчатый коллектор или через отверстия в тыльной стороне стабилизатора. Выхлопные газы ГТУ поступают в зону горения из обтекающего стабилизатор потока (рис. 1). Такой способ сжигания топлива отличается высокой интенсивностью процесса смесеобразования при малой длине факела по потоку. Режим горения рассматриваемого диффузионного факела близок к режиму турбулентного горения гомогенных смесей [1].
Наиболее сложная из известных моделей КСДТ [2] представляет собой несколько рядов фронтовых устройств диффу-зионно-стабилиза-торного типа, включающих уголковые стабилизаторы и находящиеся внутри них газовые коллекторы со струйными форсунками подачи топлива. Фронтовые устройства собираются из модулей, в которых на кромках уголковых стабилизаторов установлены поперечные уголки с короткими топливными коллекторами.
Решение поставленной задачи осуществлялось с помощью программного продукта Flow Vision, который в настоящее время широко используется для научно-исследовательских работ и решения прикладных задач газовой динамики.
Реализованная в программном комплексе Flow Vision модель массопереноса «Горение» предполагает одну брутто-реакцию с одним или двумя конечными продуктами, скорость которой определяется конкретной моделью горения. При моделировании решаются полные (без упрощений) уравнения Навье-Стокса, уравнение энергии, записанное через полную энтальпию, и уравнения турбулентных переменных [3].
Расчет произведен на основе математической модели Аррениуса-Магнуссена. Для
всех компонентов, кроме окислителя, решаются конвективно-диффузионные уравнения. Окислитель является нерасчетным компонентом. Его массовая доля определяется из соотношения [4]
y = I Yo - Af • /1, если Yo > Af • ib
'o = 1 „
10 иначе,
(1)
М = У* - Уг,
где У0, У0, Уf, Ус - восстановленные и истинные массовые доли окислителя и горючего соответственно; ¡1 - стехеометрический коэффициент при окислителе. Восстановленные массовые доли - это параметры, полученные в результате смешения компонентов (без процесса горения), а истинные - с учетом процесса горения.
Для горючего решаются два уравнения: однородное для восстановленной массовой доли (2) и неоднородное для истинной массовой доли (3):
д(у;)
дт
)
дт
7(pYfV )-7 (р YfV )-
VJ;
t,ett
= 0;
VJ
t,ett
■■-Wf
(2)
(3)
где р - плотность, кг-м- ; V - векторный дифференциальный оператор (оператор Гамильтона); т - время, с; V- скорость, м-с-1; -диффузионный поток горючего (эффективное значение), кг-м-2-с-1; Щ - скорость реакции горения, кг-м-3-с-1.
Скорость реакции горения определяется уровнем пульсаций восстановленной массовой доли горючего [5]:
1 1 -у у Щ ~ +
Wan Wub Wkin = ATanbse-BTabs (р Yf )n_ (р Yo)n_o;
\ 0,25
P ь I Yf
Wub = C
Ц8
(
Y = min
Pk 2
1 Jf
' Yf'
s ■ (v Yo P—min I Yf,— k { f 1
(4)
(5)
(6) (7)
где Щкп - кинетическая скорость процесса горения; Щигь - скорость турбулентного смещения горючего и окислителя; у - параметр, определяющий вес турбулентной реакции горения в модели «Аррениус-Магнуссен»; А, В, п, п_о, п_Я - параметры констант кинетической реакции; ТэЬв - абсолютная температура, К; Уо - истинная массовая доля окислителя; С -параметр константы турбулентной реакции; ц - молекулярная динамическая вязкость, Па-с; е - скорость диссипации турбулентной энергии, м2-с-3; к - турбулентная энергия, м2-с-2; дя - дисперсия восстановленной массовой доли горючего.
Для дисперсии восстановленной массовой доли горючего gf в данной модели решается конвективно-диффузионное уравнение [3, 4]:
df + V(pVgf) = IV(TVgf) + 2,8ц (gf )2 -2p —,
(8)
где ц - турбулентная динамическая вязкость, Па-с.
Трехмерная модель участка газохода, в котором установлена КСДТ, представляет собой параллелепипед размерами 11630x9590 мм и длиной 15 м. В данной модели принято последовательное расположение ГПК, КСДТ и ГПСВ. Температура входящего потока газов в КСДТ в различных режимах работы ГТУ изменяется в диапазоне от 100 до 140 °С, а коэффициент избытка воздуха а - от 2,8 до 5,7. Характеристики топлива (природного газа) выбраны из существующей базы веществ Flow Vision. Входящий поток выхлопных газов ГТУ создан в пользовательской базе веществ. Для решения поставленной задачи исследовательской работы проведены многовариантные расчеты с различным составом газов потока ГТУ (содержанием окислителя) и его температурой.
Первый этап расчетов - определение термодинамических характеристик (скорость потока, распределение окислителя и топлива по потоку) газовой смеси без активации процесса горения. Расчеты на данном этапе завершаются при условии стабилизации массовых долей окислителя и топлива в выходном сечении модели.
Второй этап - в потоке газовой смеси с помощью модификатора «поджег» в расчетной модели инициируется процесс горения. Исследовательская модель оснащена постоянно работающей дежурной зоной горения, обеспечивающей устойчивость пламени. Расчет завершается при условии стабилизации искомых параметров в выходном сечении.
На основании результатов расчетов разработана номограмма изменения величины q3 в зависимости от количества сжигаемого топлива и параметров потока газов ГТУ: его температуры и объемной концентрации кислорода на входе в КСДТ (рис. 2). В правой части рис. 2 выделена область параметров работы ГТЭ-110 в диапазоне относительных нагрузок пГТУ = 0,3-1,0 и температур наружного воздуха £нв = -30-30 оС, в левой части - диапазон расходов топлива в КСДТ для выработки тепловой мощности 34,2 МВ.
Анализ полученных зависимостей (рис. 2) показывает, что при стремлении объемной концентрации О2 на входе в КСДТ к 12 % величина q3 увеличивается в большей степени. Причем чем меньше О2 во входящем потоке газов, тем более заметно становится влияние температуры газов &газ.
Например, для параметров Эгаз = 130оС и Вксдт = 0,9 кг/с уменьшение объемной концентрации О2 с 17 до 16 % вызовет увеличение д3 на 0,11 %, а уменьшение объемной концентрации О2 с 14 до 13 % приведет к увеличению д3 на 0,23 %. Это обусловлено уменьшением количества кислорода, приходящегося на каждую долю топлива. Увеличение температуры газов на входе в КСДТ Эгаз снижает величину д3, причем данное влияние наиболее выражено при меньших значениях объемной концен-
трации О2. Это объясняется увеличением скорости процесса горения. Так, при постоянном значении объемной концентрации О2 = 16 % и Вксдт = 0,9 кг/с возрастание Эгаз на 60оС (с 90 до 150оС) уменьшит д3 на 0,13 %, а при объемной концентрации О2 = 14 % и тех же самых условиях уменьшение д3 составит 0,36 %.
Выделенные области рабочих параметров позволяют сделать вывод, что величина д3 изменяется от 0,5 до 1,04 %.
Рис. 2. Зависимости величины потери тепла с химической неполнотой сгорания д3 от количества сжигаемого топлива, температуры газов и объемной концентрации кислорода на входе в КСДТ
По полученной номограмме (рис. 2) при тех же параметрах входящего потока, которые использовались при расчетах с применением нормативного метода, величина q3 равна 0,8 %. Отличие результатов расчетов по нормативному методу и с использованием компьютерной модели Flow Vision составляет 0,26 %. Данная разница является несущественной, что говорит о возможности применения обоих методов. Таким образом, величина потери теплоты от химической неполноты сгорания в указанных ранее условиях работы для расчетов эффективности работы КУ с КСДТ может быть принята равной 0,5-1,0 %.
При установке КСДТ в переходном газоходе, соединяющем диффузор ГТУ с КУ, рекомендовано [1] располагать первую по ходу газов теплообменную поверхность на расстоянии 4-5 м от горелочных устройств для выравнивания температурного поля и стабилизации потока. Такое расстояние обусловлено высокими скоростями потока на входе в КУ (свыше 12 м/с). Скорость входящего в КСДТ потока газов при расположении её за ГПК во всем рабочем диапазоне ГТУ не превышает 7 м/с. При указанных ранее параметрах работы КСДТ длина пламени в разработанной модели газохода в программном продукте Flow Vision со-
ставляет не более 0,7 м. Выравнивание температурного поля и стабилизация потока осуществляется на расстоянии 1,3-1,7 м. На основании разработанной модели газохода при данной компоновке поверхностей нагрева и КСДТ за ГПК КУ первый ряд труб ГПСВ рекомендуется располагать на расстоянии 2 м от горелочных устройств.
Список литературы
1. Цанев С.В., Буров В.Д., Ремезов А.Н. Газотурбинные и парогазовые установки тепловых электростанций: учеб. пособие для вузов / под ред. С.В. Цанева. - М.: Изд-во МЭИ, 2002. - 584 с.
2. Пат. 2447364 Российская Федерация МПК F23G7/06. Камера дожигания / Д.Ю. Бантиков, В.И. Васильев, В.Н. Лавров, Ю.И. Цыбизов, В.П. Чикин, Б.Ю. Шариков; опубл. 10.12.2011.
3. Маркова Т.В, Жлуктов С.В. Развитие моделей горения в ПК FlowVision ООО «Тесис», г. Москва [Электронный ресурс]. - Режим доступа: http://www.tesis.com.ru/ infocenter/downloads/flowvision/fv_es12_tesis3.pdf
4. Работа с документацией Flow Vision Help [Электронный ресурс]. - Режим доступа: http://www.flowvision.ru/ webhelp/ fvru_30804/
5. Аксенов, А.А, Похилко, В.И., Тишин, А.П. Исследование двухступенчатого сжигания метана в вихревой горелке [Электронный ресурс]. - Режим доступа: http://www.tesis.com.ru/infocenter/downloads/flowvision/fv_me than.pdf
References
1. Tsanev, S.V., Burov, V.D., Remezov, A.N. Gazo-turbinnye i parogazovye ustanovki teplovykh elektrostantsiy [Gas turbine and combined-cycle power units of thermal power plants]. Moscow, Izdatel'stvo MEI, 2002. 584 p.
2. Bantikov, D.Yu., Vasil'ev, V.l., Lavrov, V.N., Tsybizov, Yu.l., Chikin, V.P., Sharikov, B.Yu. Kamera dozhi-ganiya [Afterburner combustion chambers]. Patent RF, no. 2447364, 2011.
3. Markova, T.V, Zhluktov, S.V. Razvitie modeley goreniya v PK FlowVision OOO «Tesis», g. Moskva [The development of combustion models in the PC Flow Vision]. Available at:
http://www.tesis.com.ru/infocenter/downloads/flowvision/fv_es12_t esis3.pdf
4. Rabota s dokumentatsiey Flow Vision Help [Working with Flow Vision Help documents]. Available at: http://www.flowvision.ru/ webhelp/ fvru_30804/
5. Aksenov, A.A, Pokhilko, V.I., Tishin, A.P. Issle-dovanie dvukhstupenchatogo szhiganiya metana v vik-hrevoy gorelke [Investigation of two-stage methane combustion in swirl burners]. Available at: http://www.tesis.com.ru/infocenter/downloads/flowvision/fv_ methan.pdf
Малков Евгений Сергеевич,
ФГБОУВПО «Ивановский государственный энергетический университет имени В.И. Ленина»,
аспирант, инженер кафедры тепловых электрических станций,
адрес: г. Иваново, ул. Рабфаковская, д. 34, кор. В, ауд. 408,
телефон (4932) 26-99-31,
e-mail: admin@tes.ispu.ru
