Известия РАН. Энергетика, 2020, № 2, стр. 133-140

Исследование эффективности дожигания топлива в дополнительной камере сгорания ГТУ, имеющей контактный теплообменник для подогрева подпиточной сетевой воды

Е. Л. Степанова 1*, П. В. Жарков 1

1 Институт систем энергетики им. Л.А. Мелентьева СО РАН
Иркутск, Россия

* E-mail: step@isem.irk.ru

Поступила в редакцию 17.01.2020
После доработки 25.02.2020
Принята к публикации 27.02.2020

Полный текст (PDF)

Аннотация

В работе выполнено исследование энергетической и экономической эффективности двух теплофикационных газотурбинных установок, технологические схемы которых включают одну и две камеры сгорания топлива, соответственно. Обе установки оснащены контактными теплообменниками для подогрева подпиточной сетевой воды. Разработаны математические модели конструкторско-поверочного расчета газотурбинных установок, которые позволяют проводить один конструкторский и несколько поверочных расчетов в режимах с различными температурами наружного воздуха. Модели были созданы с помощью разработанного в Институте систем энергетики им. Л.А. Мелентьева СО РАН программно-вычислительного комплекса Система машинного построения программ. Проведены оптимизационные расчеты режимов работы газотурбинных установок по критерию минимума цены электроэнергии. Исследования были выполнены при одинаковых климатических условиях работы, температурных графиков тепловой сети, диапазонах электрических мощностей. Полученные результаты позволили сделать вывод о том, что при работе данных теплофикационных газотурбинных установок при одинаковых климатических условиях, электрических нагрузках и необходимости обеспечения потребителей тепловой нагрузкой определенного уровня, применение дожигания топлива во второй камере сгорания оправданно с точки зрения экономической и энергетической эффективности.

Ключевые слова: газотурбинные установки, математическое моделирование энергоустановок, оптимизационные расчеты режимов работы энергоустановок, утилизация тепла дымовых газов, энергетическая и экономическая эффективность

ВВЕДЕНИЕ

В настоящее время большое количество строящихся и реконструируемых энергоблоков в мире, в том числе и России, содержат в своем составе газотурбинные установки (ГТУ). Относительно традиционных паросиловых установок, ГТУ имеют основные преимущества в виде высокой маневренности и, как правило, меньшей удельной стоимости. Следует отметить, что наибольшую энергетическую и экономическую эффективности возможно достичь при работе ГТУ в режиме когенерации, а также с помощью использования технологий утилизации тепла уходящих газов [14]. ГТУ мощностью до 60 МВт не теряют на сегодняшний день своей актуальности для энергетической отрасли промышленности. Наиболее простую компоновку имеют одновальные ГТУ с одной камерой сгорания [5]. В определенных случаях в целях стабилизации параметров рабочего тела (например, сетевой воды), а также увеличения тепловой мощности ГТУ возможно применение дожигание топлива в среде выходных газов перед входом в котел-утилизатор. В этом случае важно комплексно оценить целесообразность применения той или иной технологической схемы установки, с учетом ее энергетической и экономической эффективности. Провести оценку вариантов технологических схем ГТУ, их оптимальных конструктивных характеристик, тепловых и электрических нагрузок в различных режимах, эффект от внедрения технологий утилизации теплоты уходящих газов можно с помощью математического моделирования и оптимизации. В Институте систем энергетики им. Л.А. Мелентьева (ИСЭМ) СО РАН разработан программно-вычислительной комплекс Система машинного построения программ (СМПП), который позволяет создавать математические модели и проводить конструкторские и поверочные расчеты энергоустановок [6]. Для решения поставленной задачи были разработаны схемы двух теплофикационных ГТУ. Принципиальные технологические схемы (рис. 1) включают в себя воздушный компрессор, газовую турбину, котел-утилизатор, контактный теплообменник, водоводяной подогреватель и насосные установки. Отличие схем заключается в следующем: первая – имеет одну камеру сгорания топлива, вторая – две камеры сгорания. Схемы не имеют пикового источника тепла, регулирование тепловой нагрузки ГТУ осуществляется за счет байпасирования части расхода уходящих газов помимо котла-утилизатора и контактного теплообменника.

Рис. 1.

Принципиальные технологические схемы теплофикационных ГТУ. Схемы (а), (б): 1 – воздушный компрессор; 2 – камера сгорания топлива; 3 – газовая турбина; 4 – котел-утилизатор; 5 – контактный теплообменник; 6 – водоводяной подогреватель; 7-9 – насосы. Схема (а): 10 – дополнительная камера сгорания топлива.

Созданы математические модели конструкторско-поверочного расчета обеих теплофикационных ГТУ. Разработка велась на основе моделей конструкторского расчета, позволяющих определить геометрические размеры теплообменников, номинальные расходы, давления газа на входе и выходе из газовой турбины и моделей поверочных расчетов, которые при полученных конструктивных характеристиках определяют параметры теплоносителей (газа, воздуха и воды). Модели конструкторско-поверочного расчета ГТУ позволяют проводить один конструкторский и несколько поверочных расчетов в режимах с различными температурами наружного воздуха. Конструкторский расчет осуществляется в номинальном режиме с определением конструктивных характеристик оборудования ГТУ. Поверочные расчеты выполняются для характерных режимов с различными тепловыми нагрузками и продолжительностями в течение года. На рисунке 2 приведена блок-схема расчета теплофикационных ГТУ.

Рис. 2.

Блок-схема расчета теплофикационной ГТУ.

В качестве критерия оптимизации при проведении оптимизационных исследований принималась цена электроэнергии при заданных цене отпускаемого тепла и внутренней норме возврата капиталовложений, определяющей уровень экономической эффективности инвестиционного проекта [7]. В данной работе использован новый эффективный метод оптимизации, который обеспечивает лучшую сходимость, более точное приближение к оптимуму и имеет большее быстродействие по сравнению с ранее использованными авторами методами [8, 9]. Основная идея метода состоит в том, что решение систем уравнений, описывающих ТЭУ, передается с уровня математической модели установки на уровень решения оптимизационной задачи. В результате точное решение системы уравнений достигается лишь в оптимальной точке [10].

Математическая постановка решаемой задачи имеет следующий вид:

(1)
$\mathop {\min }\limits_{{{Q}^{p}},{{x}_{k}},{{B}_{k}},{{x}_{i}},{{B}_{i}}} {{C}^{{{\text{элек}}}}}$
при условиях
(2)
${{S}_{k}} = f\left( {{{x}_{k}},{{B}_{k}},{{\gamma }_{k}},{{Q}^{p}}} \right),$
(3)
${{G}_{k}}\left( {{{x}_{k}},{{B}_{k}},{{\gamma }_{k}},{{Q}^{p}}} \right) \geqslant 0,$
(4)
${x_{k}^{{{\text{min}}}} \leqslant {{x}_{k}} \leqslant x_{k}^{{{\text{max}}}}},$
(5)
${{N}_{k}} = f\left( {{{Q}^{p}},{{B}_{k}},{{x}_{k}},{{\gamma }_{k}}} \right),$
(6)
${{{N}^{{{\text{min}}}}} \leqslant {{N}_{k}} \leqslant {{N}^{{{\text{max}}}}}},$
(7)
${{G}_{i}}\left( {{{x}_{i}},{{B}_{i}},{{\gamma }_{i}},{{S}_{k}},{{Q}_{i}}} \right) \geqslant 0,$
(8)
${x_{i}^{{{\text{min}}}} \leqslant {{x}_{i}} \leqslant x_{i}^{{{\text{max}}}}},\,\,\,\,{i = 1,...,n},$
(9)
${{K}_{{{\text{об}}}}} = \theta \left( {{{S}_{k}},{{d}_{{{\text{уд}}}}}} \right),$
(10)
${{K}_{{{\text{сум}}}}} = {{K}_{{{\text{об}}}}} + {{K}_{{{\text{ст}}}}},$
(11)
${{Q}_{i}} = {{q}_{i}}\left( {{{Q}^{p}},{{\gamma }_{i}}} \right),$
(12)
${{N}_{i}} = f\left( {{{Q}_{i}},{{B}_{i}},{{x}_{k}},{{x}_{i}},{{\gamma }_{i}}} \right),$
(13)
${{Q}_{{{\text{год}}}}} = \sum\limits_{i = 1}^n {{{Q}_{i}}{{T}_{i}}} ,\,\,\,\,{{Э}_{{{\text{год}}}}} = \sum\limits_{i = 1}^n {\left( {{{N}_{i}} - N_{i}^{{{\text{сн}}}}} \right){{T}_{i}}} ,\,\,\,\,{{B}_{{{\text{год}}}}} = \sum\limits_{i = 1}^n {{{B}_{i}}{{T}_{i}}} .$
(14)
$IRR\left( {{{B}_{{{\text{год}}}}},{{Q}_{{{\text{год}}}}},{{Э}_{{{\text{год}}}}},{{C}^{{{\text{теп}}}}},{{C}^{{{\text{элек}}}}},{{K}_{{{\text{сум}}}}},\varphi } \right) = IR{{R}_{Z}},$
где Cэлек – цена электрической энергии; xk – вектор независимых оптимизируемых параметров, определяющих конструктивные характеристики установки (параметры цикла, конструктивные параметры элементов и режимные параметры в номинальном режиме); xi – вектор оптимизируемых режимных параметров в i-м режиме (индексом i обозначаются параметры, относящиеся к характерным режимам, в которых проводятся поверочные расчеты); Bk – часовой расход топлива газовым котлом в номинальном режиме; Bi – расход топлива газовым котлом в i-м режиме; C теп – цена тепловой энергии; Стоп – стоимость топлива; Qp – расчетная тепловая нагрузка; R – постоянная составляющая в цене тепла; Sк – вектор конструктивных характеристик установки; γк – вектор исходных данных, определяющих внешние условия работы установки в номинальном режиме; Gklk-мерная векторная функция ограничений-неравенств в номинальном режиме; Nk – полная электрическая мощность в номинальном режиме; N min, Nmax – минимальное и максимальное значение электрической мощности ПТУ в номинальном режиме; Gilp-мерная векторная функция ограничений-неравенств в i-м режиме; γi – вектор исходных данных, определяющих внешние условия работы в i-м режиме; Kоб – капиталовложения в оборудование; dуд – вектор удельных стоимостей элементов оборудования; Kсум – суммарные капиталовложения; Kст – капиталовложения, учитывающие непредвиденные затраты и затраты на строительство; Bгод – годовой расход топлива; Qгод – годовая теплопроизводительность установки; Эгод – годовой полезный отпуск электроэнергии; Ti – продолжительность i-го режима; Qi – отпуск тепловой энергии в i-м режиме; Ni – полная электрическая мощность в i-м режиме; $N_{i}^{{{\text{сн}}}}$ – расход электроэнергии на собственные нужды в i-м характерном режиме; IRR – внутренняя норма возврата капиталовложений; IRRZ – заданное значение внутренней нормы возврата капиталовложений; φ – вектор параметров, определяющих условия кредитования и налогообложения; $x_{k}^{{\min }},$ $x_{k}^{{\max }},$ $x_{i}^{{\min }},$ $x_{i}^{{\max }}$ – векторы минимальных и максимальных значений xk и xi соответственно; n – число режимов, в которых проводятся поверочные расчеты (общее число характерных режимов n +1).

Для проведения оптимизационных исследований были назначены оптимизируемые параметры цикла в номинальном режиме, конструктивные и режимные оптимизируемые параметры, общее количество которых для каждой ГТУ составляло 45 и 48 параметров соответственно. Список оптимизируемых параметров конструкторского расчета включал: расчетную тепловую нагрузку потребителя; расход топлива; внутренние, наружные диаметры и шаги труб поверхностей нагрева котла-утилизатора и теплообменников; массовые скорости воды; температура газа на входе в газовую турбину; долю от общего расхода уходящих газов после газовой турбины на котел-утилизатор. Список режимных оптимизируемых параметров в характерных режимах (поверочные расчеты) включал: расход топлива; расход греющей воды подогревателя подпитки сетевой воды; температуру газа на входе в газовую турбину; долю от общего расхода уходящих газов после газовой турбины на котел-утилизатор. При оптимизационных исследованиях были наложены ограничения на: максимально допустимую температуру газа на входе в газовую турбину, неотрицательность температурных напоров в котле-утилизаторе, механические напряжения металла труб теплообменников (в количестве 244 ограничений для ГТУ, имеющей одну камеру сгорания и 249 – для ГТУ с двумя камерами сгорания).

РЕЗУЛЬТАТЫ ОПТИМИЗАЦИОННЫХ ИССЛЕДОВАНИЙ

Оптимизационные расчеты режимов работы обеих теплофикационных ГТУ проводились для следующих условий: расчетная температура наружного воздуха принята равной –25°С; отопительный период разбит на три периода времени, каждому периоду соответствовала своя температура наружного воздуха – –20°С, –10°С, +1.5°С и свой характерный режим с определенными тепловыми нагрузками; продолжительности отопительных периодов составляли 418 ч, 1316 ч, 3176 ч, неотопительного периода – 3850 ч. Диапазон изменения электрической мощности ГТУ в номинальном режиме был принят равным 50–60 МВт. Температурный график тепловой сети: 120/70°С. Предельно допустимое значение температуры газа перед газовой турбиной принято равным 1300°С. При расчете капиталовложений была принята следующая исходная информация [11]. Удельные стоимости составляли для: газовой турбины – 70 долл./кВт (на единицу максимальной полной мощности газовой турбины), воздушного компрессора – 50 долл./кВт, электрогенератора турбины, насосов – 60 долл./кВт, теплообменников, сделанных из низкоуглеродистой стали – 7 тыс. долл./т, из стали 20–21 тыс. долл./т, электротехнического оборудования 0.192 тыс. долл./кВт., систем, зависящих от расхода топлива 240 тыс. долл./(т/ч). Внутренняя норма возврата капиталовложений – 0.15. Исследования проводились при цене топлива равной 100 долл./т у.т. и цене тепла – 4.78 долл./ГДж.

В таблицах 1 и 2 представлены основные показатели, полученные в результате оптимизационных расчетов номинального и четырех характерных режимов работы теплофикационной ГТУ с одной камерой сгорания и ГТУ с дожиганием топлива в среде выходных газов во второй камере сгорания соответственно.

Таблица 1.  

Основные показатели работы теплофикационной ГТУ с одной камерой сгорания в номинальном и характерных режимах

Основные показатели Номинальный режим Номер характерного режима
1 2 3 4
Полезная электрическая мощность ГТУ, МВт 60.0 58.0 57.0 56.1 53.0
Тепловая нагрузка потребителей, ГДж/ч 278.0 258.3 203.5 150.0 88.0
Расход топлива ГТУ, т у.т./ч 21.5 20.8 20.1 19.6 17.0
Температура воздуха на выходе из воздушного компрессора, °С 275.9 313.7 336.3 361.5 397.3
Расход воздуха на выходе из воздушного компрессора, т/ч 792.7 667.8 655.5 643.3 623.5
Давление газа на входе в газовую турбину, МПа 1.42 1.42 1.41 1.40 1.33
Давление газа на выходе из газовой турбины, МПа 0.1 0.1 0.1 0.1 0.1
Температура газа на входе в газовую турбину, °С 849.3 1055.1 1090.9 1129.5 1190.2
Температура газа на выходе из газовой турбины, °С 465.9 511.4 534.6 560.3 587.2
Количество тепла, полученного за счет подогрева подпиточной сетевой воды, ГДж/ч 102.8 58.5 18.6 13.3 6.27
Площадь теплообменной поверхности котла-утилизатора, м2 5014.0
Площадь теплообменной поверхности водоводяного подогревателя подпитки сетевой воды, м2 1632.0
Таблица 2.  

Основные показатели работы теплофикационной ГТУ с двумя камерами сгорания в номинальном и характерных режимах

Основные показатели Номинальный режим Номер характерного режима
1 2 3 4
Полезная электрическая мощность ГТУ, МВт 60.0 58.0 57.1 56.5 56.1
Тепловая нагрузка потребителей, ГДж/ч 420.0 396.2 310.1 217.7 124.9
Расход топлива ГТУ, т у.т./ч 30.5 24.7 22.1 20.9 19.0
Температура воздуха на выходе из воздушного компрессора, °С 296.4 317.6 340.7 365.8 400.5
Расход воздуха на выходе из воздушного компрессора, т/ч 614.1 573.4 561.6 552.2 541.4
Давление газа на входе в газовую турбину, МПа 1.45 1.45 1.45 1.44 1.43
Давление газа на выходе из газовой турбины, МПа 0.1 0.1 0.1 0.1 0.1
Температура газа на входе в газовую турбину, °С 881.2 1104.0 1160.7 1195.6 1227.3
Температура газа на выходе из газовой турбины, °С 478.5 539.8 575.2 596.7 619.2
Количество тепла, полученного за счет подогрева подпиточной сетевой воды, ГДж/ч 237.4 76.9 68.5 65.2 50.9
Площадь теплообменной поверхности котла-утилизатора, м2 4298.0
Площадь теплообменной поверхности водоводяного подогревателя подпитки сетевой воды, м2 2594.0

В таблице 3 приведены основные годовые показатели расчетов режимов работы обеих теплофикационных ГТУ, на основании которых можно провести сравнительный анализ энергетической и экономической эффективности установок при работе в одинаковых климатических условиях и при электрических нагрузках одного уровня.

Таблица 3.  

Таблица основных годовых показателей расчета режимов работы теплофикационных ГТУ

Основные показатели ГТУ с одной камерой сгорания ГТУ с двумя камерами сгорания
Цена тепла, долл./ГДж 4.78
Цена электроэнергии, цент./кВт 4.27 3.57
Удельные капиталовложения, долл./кВт 617.5 583.8
Годовой расход топлива, тыс. т у.т.год 170.7 175.0
Годовой отпуск электроэнергии, тыс. МВт ч 491.5 500.8
Годовой отпуск тепла, тыс. ГДж 1148.6 1741.8
Удельный расход условного топлива на отпуск тепла, кг у.т./ГДж 36.74
Удельный расход условного топлива на отпуск электроэнергии, г у.т./кВт ч 262.0 220.0

Как видно, из приведенных результатов оптимизационных исследований (табл. 1, 2) вся тепловая нагрузка ГТУ покрывается нагревом сетевой воды в котле-утилизаторе без привлечения пикового источника. Максимальная электрическая полезная мощность обеих ГТУ достигает в номинальном режиме 60 МВт, полученная тепловая нагрузка потребителей ГТУ с одной камерой сгорания равна 278.0 ГДж/ч, для ГТУ с двумя камерами сгорания – 420.0 ГДж/ч. Количество дополнительного тепла, полученного за счет подогрева подпиточной сетевой воды с помощью контактного теплообменника находится в диапазоне 6.27–102.8 ГДж/ч для ГТУ с одной камерой сгорания и в диапазоне 50.9–237.4 ГДж/ч для ГТУ с двумя камерами сгорания. Результаты исследований показали, что ГТУ с двумя камерами сгорания при цене топлива равной 100 долл./т у.т. и цене тепла – 4.78 долл./ГДж имеет меньшую цену электроэнергии (3.57 цент./кВт), чем ГТУ с одной камерой сгорания (4.27 цент./кВт), больший отпуск электрической и тепловой энергии, и более низкий расход условного топлива на отпуск электроэнергии (220.0 г у.т./кВт ч), чем расход условного топлива на отпуск электроэнергии ГТУ с одной камерой сгорания (262.0 г у.т./кВт ч).

ЗАКЛЮЧЕНИЕ

Выполнено исследование энергетической и экономической эффективности двух теплофикационных ГТУ, технологические схемы которых включают одну и две камеры сгорания топлива соответственно. Проведены оптимизационные расчеты их режимов работы по критерию минимума цены электроэнергии. Полученные результаты позволяют сделать вывод о том, что при работе данных теплофикационных установок в одинаковых климатических условиях, с электрическими нагрузками одного уровня, а так же при необходимости обеспечения потребителей большей тепловой нагрузкой, чем максимально-возможная тепловая нагрузка ГТУ с одной камерой сгорания (без пикового источника), применение дожигания топлива в среде выходных газов во второй камере сгорания ГТУ является достаточно эффективным в энергетическом и экономическом плане мероприятием.

Работа выполнена в рамках научного проекта III.17.1.1. программы фундаментальных исследований СО РАН, рег. № АААА-А17-117030310433-6.

Список литературы

  1. Фаворский О.Н. Выбор тепловой схемы и профиля отечественной мощной энергетической ГТУ нового поколения и ПГУ на ее основе // Теплоэнергетика, 2010. № 2. С. 2–7.

  2. Ольховский Г.Г. Теплоэнергетические технологии в период до 2030 г. // Изв. РАН. Энергетика. 2008. № 6. С. 79–94.

  3. Najjar Y.S.H. Gas turbine cogeneration systems: a review of some novel cycles // Applied Thermal Engineering. 2000. P. 179–197.

  4. Bade M., Bandyopadhyay S. Analysis of gas turbine integrated cogeneration plant: process integration approach // Applied Thermal Engineering. 2015. P. 118–128.

  5. Цанев С.В., Буров В.Д., Земцов А.С., Осыка А.С. Газотурбинные энергетические установки. М.: Издательский дом МЭИ, 2011. 428 с.

  6. Клер А.М., Тюрина Э.А. и др. Оптимизационные исследования энергетических установок и комплексов. – Новосибирск: Академическое издательство “Гео”, 2016. 298 с.

  7. Клер А.М., Степанова Е.Л., Максимов А.С. Исследование эффективности теплофикационного энергоблока с противодавленческой паровой турбиной и утилизацией тепла уходящих газов // Теплофизика и аэромеханика. 2018. Т. 25. № 6. С. 963–973.

  8. Клер А.М., Деканова Н.П., Скрипкин С.К. Математическое моделирование и оптимизация в задачах оперативного управления тепловыми электростанциями. Новосибирск: Наука. Сиб. предприятие РАН. 1997. 120 с.

  9. Клер А.М., Деканова Н.П., Тюрина Э.А. Теплосиловые системы: оптимизационные исследования. Новосибирск: Наука, 2005. 236 с.

  10. Kler A.M., Zharkov P.V., Epishkin N.O. Parametric optimization of supercritical power plants using gradient methods // Energy. 2019. V. 189.

  11. Клер А.М., Захаров Ю.Б., Потанина Ю.М. Оценка эффективности впрыска воды между частями воздушного компрессора для ГТУ и ПГУ // Теплофизика и аэромеханика. 2017. Т. 24. № 3. С. 497–505.

Дополнительные материалы отсутствуют.