Impact of steam flow into a combustion chamber of a contact gas-steam installation on its energy characteristics
- Authors: Galashov N.N.1, Tubolev A.A.1, Boldushevsky E.S.1, Minor A.A.2
-
Affiliations:
- National Research Tomsk Polytechnic University
- Tomsk Generation JSC
- Issue: Vol 335, No 2 (2024)
- Pages: 48-59
- Section: Articles
- URL: https://journal-vniispk.ru/2500-1019/article/view/267047
- DOI: https://doi.org/10.18799/24131830/2024/2/4436
- ID: 267047
Cite item
Full Text
Abstract
Relevance. Reduction of natural gas consumption and emissions of harmful substances into the environment based on introduction of water vapor into a combustion chamber of a contact gas-steam installation.
Aim. To carry out numerical studies on the influence of relative steam flow into the combustion chamber of the contact gas-steam installation on its energy characteristics.
Objects. Contact gas-steam installations based on gas turbines with steam injection into the combustion chamber.
Methods. Numerical methods based on material and energy balances of systems and elements of gas-steam installations.
Results. Based on the calculation of the thermal circuit of the contact gas-steam installation, the authors have studied the influence of the relative steam flow into the combustion chamber on its energy characteristics. It was determined that the absolute electrical efficiency of the contact gas-steam installation increases linearly with growth of relative steam flow into the combustion chamber. The range of changes in the relative steam flow into the combustion chamber strongly depends on the temperature of the gases behind the combustion chamber and the compression ratio in the air compressor; the smaller these parameters, the greater the range of changes. The maximum efficiency of 56% for all options is achieved at the maximum relative steam flow into the combustion chamber. It was established that the excess air coefficient, depending on the relative steam flow rate, decreases linearly, and the higher the temperature of the gases behind the combustion chamber and the compression ratio in the air compressor, the greater the rate of decline and the smaller the range of changes in the relative steam flow rate. It was revealed that the efficiency coefficient strongly depends on the relative steam flow into the combustion chamber, the temperature of the gases behind it and the degree of compression in the air compressor; with increasing these parameters, it increases linearly. It was determined that the temperature of the gases at the outlet of the gas turbine also strongly depends on the relative flow of steam into the combustion chamber, the temperature of the gases at its outlet and the compression ratio in the compressor. With an increase in the relative flow of steam into the combustion chamber, this temperature increases linearly from 600 to 700°C, while the higher the temperature of the gases at the outlet of the combustion chamber and the compression ratio in the compressor, the higher the temperature of the gases at the outlet of the gas turbine. The authors revealed the dependence of useful work on a gas turbine shaft on the relative steam flow into the combustion chamber. With an increase in the relative steam flow, the useful work on the gas turbine shaft increases along the branch of the parabola. The higher the temperature of the gases behind the combustion chamber and the compression ratio in the compressor, the steeper the branch of the parabola, but the smaller the range of changes in the relative steam flow. It was established that with an increase in the relative steam flow, the gas flow to the gas turbine decreases according to a hyperbola. Moreover, the lower the temperature of the gases behind the combustion chamber and the compression ratio in the compressor, the more the gas flow to the gas turbine drops.
Full Text
Введение
Быстрый рост стоимости природных энергоресурсов, таких как уголь, природный газ и нефть, требует разработки высокоэффективных энергосберегающих технологий снижения их расхода при сжигании. Как показано в [1]: «техническое перевооружение отечественных ТЭС, работающих на природном газе, с использованием высокоэкономичных парогазовых установок позволит более чем в 1,5 раза снизить удельные расходы топлива на производство электроэнергии».
В настоящее время для выработки электроэнергии с максимальным КПД используется комбинация газотурбинного цикла Брайтона и паротурбинного цикла Ренкина.
При этом возможны два способа комбинации этих циклов: раздельный – при котором сначала в газотурбинной установке (ГТУ) в цикле Брайтона при сжигании жидкого или газообразного топлива в камере сгорания (КС) теплота сгорания топлива преобразуется в газовой турбине (ГТ) в механическую мощность, которая в воздушном компрессоре (ВК) служит для подачи в камеру сгорания воздуха с необходимым давлением, а в газовой турбине – для привода ВК и генератора электрического тока. Дальше продукты сгорания с температурой 550…700 °С из ГТ входят в котел-утилизатор, где позволяют получить перегретый водяной пар с температурой 500…650 °С, который в цикле Ренкина вырабатывает мощность в паровой турбине для привода генератора электрического тока. Для подачи топлива в КС служит топливный компрессор (ТК). Этот способ хорошо освоен на современных мощных парогазовых установках (ПГУ) и позволил получить электрический КПД до 64 % [2, 3]. Недостатком способа является сложность схемы, зависимость режимов работы цикла Ренкина от цикла Брайтона. Сложность пуско-остановочных операций.
Второй способ заключается в том, что оба термодинамических цикла работают в одной контактной газопаровой установке (КГПУ). При этом в камеру сгорания вводится под давлением вода или водяной пар, где за счет тепловой энергии сжигаемого топлива они переходят в состояние перегретого пара. В турбине работу совершает смесь продуктов сгорания и перегретого водяного пара. Поскольку удельная теплоемкость и плотность водяного пара выше, чем у продуктов сгорания, то 1 кг смеси создает мощность больше, чем 1 кг продуктов сгорания. При этом сокращается массовый расход воздуха на кг сжигаемого топлива, что снижает затраты мощности на привод воздушного компрессора. Ввод пара или воды сильно влияет на параметры и тепловое состояние камеры сгорания. Одно из первых исследований камеры сгорания с вводом пара было выполнено В.А. Зысиным и Б.В. Турчаниновым на экспериментальной камере сгорания ГТ-700 НЗЛ [4]. Испытания показали, что при относительном массовом расходе пара 0,0875 к расходу воздуха температура стенки жаровой трубы снизилась в среднем на 200–250 °С, а химический недожег топлива практически отсутствовал. Основы работы установок по второму способу разработаны В.А. Зысиным в [5], где автор назвал такие установки «контактными газопаровыми установками». В это же время КГПУ начали разрабатываться за рубежом, где получили название STIG (Stеam Injected Gas Turbine).
Дальнейший теоретический анализ работы КГПУ по второму способу приведен в книге [6], где авторы обосновали влияние ввода пароводяного рабочего тела в камеру сгорания ГТУ. Были определены оптимальные и предельные параметры продуктов сгорания, снижение коэффициента избытка воздуха в камере сгорания и затраты мощности на привод ВК, уменьшение выбросов NOx в результате понижения температуры горения топлива, повышение примерно в два раза полезной удельной мощности ГТ. Также показано, что ввод воды в камеру сгорания может не только повысить КПД ГТУ, но и понизить из-за сжигания дополнительного расхода топлива на парообразование и перегрев пара, поэтому более эффективен ввод пара.
В ряде работ [7–19] показано, что ввод пара в КС позволяет снизить температуру металла КС и увеличить надежность ее работы, а также существенно уменьшить выбросы NOx и СО. В работах [14–19] определено, что оптимальный расход вводимого пара по влиянию на образование NOx и СО находится в пределах 1,5…2,5 кг пара/кг топлива. Для увеличения КПД требуется расход пара в КС больше 2,5 кг пара/кг топлива, для этого расход пара делят на две части: экологический вводят в зону горения топлива для поддержания допустимых выбросов оксидов азота NOx, а энергетический вводят за зоной горения для поддержания необходимой температуры на выходе КС.
Достоинством второго способа является более простая и дешевая схема КГПУ при меньших капитальных затратах. К недостаткам второго способа следует отнести повышенные потери теплоты с уходящими газами, а также вводимого в КС пара. Чтобы уменьшить эти потери за ГТ устанавливают котел-утилизатор, а на его выходе – теплоутилизатор для конденсации пара из парогазовой смеси продуктов сгорания, который используется для подготовки добавочной воды [20–23].
Расчетная схема контактной газопаровой установки
Целью данной работы является проведение исследований по влиянию относительного расхода пара в камеру сгорания контактной газопаровой установки на ее энергетические характеристики. Для расчета взята простая одновальная схема КГПУ с разомкнутыми циклами Брайтона и Ренкина, которая приведена на рис. 1.
Рис. 1. Схема газопаровой установки: КВОУ – комплексная воздухоочистная установка; ВК – воздушный компрессор; КС – камера сгорания; ТК – топливный компрессор; ГТ – газовая турбина; ГЭ – генератор электрического тока
Fig. 1. Cycle diagram of combined-cycle plant: CAPP – complete air purification plant; AС – air compressor; CC – combustion chamber; FC – fuel compressor; FT – fuel turbine; EG – electric generator
Числа на схеме определяют точки процессов на входе и выходе элементов.
Схема работает следующим образом. В КВОУ происходит очистка воздуха от загрязнений, при этом протекает изотермический процесс 0–1 со снижением давления от Р0 до Р1. ВК в политропном процессе 1–2 сжимает воздух от давления Р1 до Р2 с увеличением температуры от Т1 до Т2. ТК в политропном процессе 3–4 сжимает газообразное топливо от давления Р3 до Р4 с увеличением температуры от Т3 до Т4. Водяной пар в состоянии 5 подается в КС при давлении Р5 и температуре Т5. В КС в процессе окисления углеводородного топлива кислородом воздуха выделяется теплота и образуются продукты сгорания: СО2 и Н2О, которые с N2, избыточным воздухом и введенным паром с температурой Т6 идут на ГТ, где, расширяясь в политропном процессе 6–7, выходят с температурой Т7 и давлением Р7.
Методика расчета тепловой схемы КГПУ
В данной работе в качестве основы использована методика расчета схемы одновальной ГТУ без ввода дополнительного пара в КС, описанная в [24]. Эта методика дополнена определением количества вводимого дополнительного пара в КС, а также расчетом параметров потоков не как идеальных газов, а как реальных веществ, что отличает ее от методики, рассмотренной в [25], в которой также не учитывается изменение средней теплоемкости в процессах сжатия и расширения при изменении температуры. Удельные расходы определяются не по отношению к расходу воздуха, а по отношению к расходу топлива.
В методике приняты следующие обозначения: В, b – абсолютный и удельный расход топлива; D, d – абсолютный и удельный расход вводимого пара; G, g – абсолютные и удельные расходы воздуха и продуктов сгорания; N – мощность; h – КПД; Н – удельная работа.
Приняты следующие обозначения параметров: температура: (Т), К; (t), °С; давление (Р), МПа; удельная энтальпия (h), кДж/кг; удельная энтропия (s), кДж/(кг·К).
Параметры воды и водяного пара, воздуха, топлива и продуктов сгорания рассчитываются по функциям, определяемым в лицензированной библиотеке «CooProp» [26], которая встроена в надстройки пакета Excel и позволяет с помощью функции PropsSi определить параметры для 110 веществ и их смесей. Эта функция позволяет по двум независимым параметрам определить для жидкости и перегретого пара, а также на пограничных кривых: удельные энтальпию, энтропию, объем, внутреннюю энергию, теплоемкость при постоянном давлении и объеме, функцию Гиббса, давление, температуру, плотность, вязкость, теплопроводность, поверхностное натяжение, критические давление и температуру. При расчете параметров на пограничных кривых используется второй параметр Q=0 для жидкости и Q=1 для насыщенного пара. Достоинством функции является то, что большинство параметров в ней имеет общепринятое обозначение, а названия веществ записываются в латинском написании.
В следующем примере определяется энтальпия воды в кДж/кг. Термодинамическое состояние задается двумя переменными: Р и T, которые устанавливаются как 1 бар (100000 Па) и 15 °C (288,15 К) соответственно:
h=PropsSI("H";"T";15+273,25;"P";100000;"Water")/1000.
При расчете параметров с помощью функции PropsSI за начало расчета принята температура 0 К. Так как состав топлива приводится при нормальных условиях, т. е. при температуре 0 °С и давлении 101325 Па, то при этих параметрах расчетным путем определяется теплота сгорания топлива, поэтому для согласования с теплотой сгорания необходимо приведение параметров к нормальным условиям, для этого при расчете параметров необходимо дополнительно вычитать функцию PropsSI при температуре 273,15 К и давлении 101325 Па. Для приведенного выше примера это будет записано как
h=PropsSI("H";"T";15+273,25;"P";100000;"Water")/1000–PropsSI("H";"T";273,25;"P";101325;"Water")/1000.
Подстрочные индексы у параметров обозначают: числа – вход и выход для соответствующих элементов (рис. 1); «вх» – вход; «вых» – выход; «вк» и «тк» – воздушный и топливный компрессор; «кс» – камера сгорания; «гт» – газовая турбина; «в» – воздух; «гпс» – смесь продуктов сгорания и добавленного пара; «нв» – наружный воздух; «изб» – избыточный; «п» – пар; «s» – изоэнтропийный процесс; «н» – низшая теплота сгорания топлива; «м» – механический; «эг» – электрогенератор; «эп» – электропривод; «э» – электрическая; «е» – эффективная; «ут» – утечки. Комплексные подстрочные индексы записываются через точку или дефис.
Блок-схема расчета КГПУ приведена на рис. 2.
Рис. 2. Блок-схема расчета КГПУ
Fig. 2. Block diagram of contact gas-steam installation calculation
Исходные данные:
- наружный воздух: Рнв, Tнв, jнв;
- природный газ: СН4, С2H6, C3H8 и т. д.;
- природный газ в газопроводе: Pт, Tт;
- потери давления: dРквоу, dРвк.вых, dРтк.вых, dРкс, dРкс-гт, dРгт.вых;
- КПД: hвк, hтк, hм.вк, hм.тк, hгт, hгэ, hэп, hкс;
- потери от утечек в ВК αут;
- степень сжатия воздушного компрессора e;
- температура на выходе КС Tгпс;
- относительный расход вводимого пара d;
- электрическая мощность генератора Nэ.
Расчет параметров воздуха, топлива и продуктов сгорания
Для расчета параметров воздуха, топлива и продуктов сгорания используются уравнения [27]. По ним для заданного состава топлива рассчитываются: теоретическая масса воздуха для сжигания 1 кг топлива L0, кг/кг; удельная низшая теплота сгорания топлива Qн, кДж/кг; удельные массы продуктов сгорания кг/кг топлива: gN2; gCO2; gH2O. Расчет ведется для нормальных условий: Р=101325 Па; T=273,15 K. При этих условиях известны плотность r, молекулярная масса m и газовая постоянная R всех рассчитываемых веществ (табл. 1).
Таблица 1. Плотность, молекулярная масса и газовая постоянная
Table 1. Density, molecular weight and gas constant
Параметры Parameters | воздух air | пар steam | СО2 carbon dioxide | N2 nitrogen | CH4 methane | C2H6 ethane | C3H8 propane | C4H10 butane | топливо fuel |
r, кг/м3 density, kg/m3 | 1,293 | 0,804 | 1,977 | 1,25 | 0,717 | 1,356 | 2,004 | 2,703 | 0,725 |
m, г/моль molecular weight, g/mol | 28,96 | 18,016 | 44,01 | 28,013 | 16,04 | 30,07 | 44,097 | 58,12 | 16,22 |
R, Дж/(кгК) gas constant, J/(kgK) | 281,10 | 461,50 | 188,22 | 296,81 | 518,36 | 276,50 | 188,55 | 143,06 | 512,52 |
Расчет параметров в процессе сжатия в компрессоре
Алгоритм расчета компрессора:
- Принимаем mк¢=0,265.
- По заданным Рвх, Рвых, Па и Твх, К определяем:
- Энтальпию на входе
hвх=PropsSI("H";"T";Твх;"P";Рвх;"вещество")/1000.
- Изобарную теплоемкость Ср=R/mк¢, где R – газовая постоянная вещества.
- Температуру на выходе компрессора Твых=Твх(1+(emк¢–1)/hк, где hк – КПД воздушного или топливного компрессора.
- Энтальпию на выходе
hвых=PropsSI("H";"T";Твых;"P";Рвых;"вещество")/1000.
- Среднюю изобарную теплоемкость
Ср.ср=(hвых–hвх)/(Твых–Твых).
- Уточняем mк=R/Ср.ср.
- Если (abs(mк – mк¢)/mк)100>0,1, то mк¢=mк, идем на 4.
- Работа, затраченная на сжатие кг воздуха или топлива, Hк=hвых–hвх.
Расчет камеры сгорания
На рис. 3 приведена схема расходов и параметров КС.
Рис. 3. Схема расходов и параметров КС
Fig. 3. Scheme of costs and parameters of СС
Давления на входе и выходе КС:
Ркс.вх=Р2·(1–dРвк.вых); Ркс.вых= Ркс.вх·(1–dРкс).
Материальный баланс КС для абсолютных расходов
Gв+В+D=Gгпс=GN2+GCO2+GH2O+Gв.изб+D,
и для относительных
αL0+1+d=gN2+gCO2+gH2O+(α–1)L0+d.
На основе уравнения материального баланса можно получить коэффициент связи расхода воздуха через компрессор с расходом газа на выходе КС с учетом утечек в компрессоре αут:
b=Gв·(1+αут)/Gгпс=α·L0·(1+αут)/(α·L0+(1+d)).
Тепловой баланс КС для абсолютных расходов
В·Qн·hкс+Gв·h2+В·h4+D·h5=Gгпс·hгпс==GN2·h6N2+GCO2·h6CO2+GH2O·h6п+Gв.из·h6в+D·h6п,
и для относительных
Qн·hкс+α·L0·h2+h4+d·h5=gN2·h6N2++gCO2·h6CO2+gH2O·h6п+(α–1)·L0·h6в+d·h6п,
где h определяются по Т и P для соответствующих точек и веществ с помощью функции PropsSI.
Из последнего уравнения находим коэффициент избытка воздуха
α=(Qн·hкс+h4+d·(h6п–h5)+L0·h6в––gN2·h6N2–gCO2·h6CO2–gH2O·h6п)/(L0·(h6в–h2)).
Энтальпия газопаровой смеси на выходе КС
hгпс=d·h6п+(α·L0–1)·h6в+gN2·h6N2+gCO2·h6CO2+gH2O·h6п.
Газовая постоянная газопаровой смеси на выходе КС
Rгпс=d·Rп+(α·L0–1)·Rв+gN2·RN2+gCO2·RCO2+gH2O·Rп.
Расчет процесса расширения в газовой турбине
Алгоритм расчета турбины:
Давление на входе и выходе ГТ:
Ргт.вх=Ркс.вых·(1–dРгт.вх); Ргт.вых= Рнв·(1+dРгт.вых).
- Принимаем mгт¢=0,22.
- Определяем изобарную теплоемкость
Ср=Rгпс/mгт¢.
- Температура на выходе ГТ Т7=Т6(1–(1–e–mгт¢))hгт.
4.Для всех веществ газопаровой смеси по Т7 и Ргт.вых определяется энтальпия с помощью функции PropsSI.
- Энтальпия газопаровой смеси на выходе ГТ
h7гпс=d·h7п+(α·L0–1)·h7в+gN2·h7N2+gCO2·h7CO2+gH2O·h7п.
- Средняя изобарная теплоемкость
Ср.ср=(h6г–h7г)/(Т6–Т7).
- Уточняем mгт=Rгпс/Ср.ср.
- Если (abs(mгт – mгт¢)/mгт)100>0,1, то mгт¢=mгт, идем на 2.
- Работа, полученная при расширении кг газов,
Hгт=h6г–h7г.
10. Полезная работа на валу ГТ, переданная на ГЭ,
Не=(Нгт–b·Hвк)·hм.
Расчет показателей установки на заданную электрическую мощность Nэ
Gгпс=Nэ·hгэ/Не;
Gвк¢=b·Gгпс;
Gвк=Gвк¢ (1– αут);
В=(Gгпс–Gвк)/(1+d);
D=d·В;
Nгт=Gгпс·Hгт;
Nвк=Gвк¢·Hвк;
j=(Nгт–Nвк)/Nгт;
hэ=Gгпс·Не·hгэ/(Qн·В);
t7=T7–273,15;
Nэ.тк=В·Hтк/hэп.
На основе представленного алгоритма написана программа расчета схемы КГПУ в пакете электронных таблиц EXCEL с расчетом параметров веществ по уравнениям базы CoolProp.
Результаты расчета и их анализ
Расчеты проведены для исходных данных:
- наружный воздух: Рнв=0,1013 МПа, tнв=15 °С, jнв=0,6;
- природный газ: CH4=98,72 %, С2Н4=0,12 %, С3Н8=0,01 %, C4H10=0,009 %, СО2=0,14 %, N2=1 %;
- природный газ на входе ТК: Pтк.вх=0,2 МПа, tтк.вх=5 °С;
- потери давления: dРквоу=0,01, dРк.вых=dРкс==dРтк-кс=dРтпк-ст=0,03;
- КПД: hвк=0,86, hтк=0,9, hм.вк=0,995, hкс=0,995, hм.тк=0,995, hгт=0,93, hэг=0,982, hэп=0,99;
- степень сжатия компрессора e: 20, 30, 40 и 50;
- температура газопаровой смеси на выходе КС tгпс: 1300, 1400, 1500 и 1600 °С;
- относительный расход пара в КС d, кг/кг топлива: 0…13;
- электрическая мощность генератора Nэ=100 МВт.
Температура вводимого пара в КС принята на 40 °С ниже температуры газопаровой смеси на выходе из ГТ. Расчеты проведены при изменении относительного расхода ввода пара в КС для вариантов: 1) tгпс=1300 °C, e=20; 2) tгпс=1400 °C, e=30; 3) tгпс=1500 °C, e=40; 4) tгпс=1600 °C, e=50. Для каждого варианта максимальная величина вводимого расхода пара d определялась по минимальному коэффициенту избытка воздуха в КС 1,05.
Результаты расчетов для наиболее значимых показателей КГПУ приведены на рис. 4–9.
Рис. 4. Электрический коэффициент полезного действия в зависимости от относительного расхода пара, температуры газопаровой смеси за камерой сгорания и степени сжатия в компрессоре
Fig. 4. Electrical efficiency depending on the relative steam flow, gas temperature behind the combustion chamber and the compression ratio in the compressor
Рис. 5. Коэффициент избытка воздуха в зависимости от относительного расхода пара, температуры газопаровой смеси за камерой сгорания и степени сжатия в компрессоре
Fig. 5. Excess air coefficient depending on the relative steam flow, gas temperature behind the combustion chamber and the compression ratio in the compressor
Рис. 6. Коэффициент полезной работы в зависимости от относительного расхода пара, температуры газопаровой смеси за камерой сгорания и степени сжатия в компрессоре
Fig. 6. Efficiency factor depending on the relative steam flow, gas temperature behind the combustion chamber and the compression ratio in the compressor
Рис. 7. Температура на выходе газовой турбины в зависимости от относительного расхода пара, температуры газопаровой смеси за камерой сгорания и степени сжатия в компрессоре
Fig. 7. Temperature at the gas turbine outlet depending on the relative steam flow, gas temperature behind the combustion chamber and the compression ratio in the compressor
Рис. 8. Полезная работа на валу газовой турбины в зависимости от относительного расхода пара, температуры газопаровой смеси за камерой сгорания и степени сжатия в компрессоре
Fig. 8. Useful work on the gas turbine shaft depending on the relative steam flow, gas temperature behind the combustion chamber and the compression ratio in the compressor
Рис. 9. Расход газов на газовую турбину в зависимости от относительного расхода пара, температуры газопаровой смеси за камерой сгорания и степени сжатия в компрессоре
Fig. 9. Gas consumption at the gas turbine depending on the relative steam consumption, gas temperature behind the combustion chamber and the compression ratio in the compressor
Зависимость электрического КПД КГПУ от d для представленных выше вариантов приведена на рис. 4. Видим, что КПД КГПУ линейно возрастает с увеличением d: для первого варианта диапазон изменения d=0…12,5; для второго – 0…10; для третьего – 0…8 и для четвертого – 0…6,5. Максимальный КПД на уровне 56 % для всех вариантов достигается при максимальном d.
На рис. 5 видим, что коэффициент избытка воздуха в зависимости от относительного расхода пара линейно падает, чем выше температура газопаровой смеси за КС и степень сжатия в ВК, тем больше темп падения и меньше диапазон изменения d.
На рис. 6 показано изменение коэффициента полезной работы в зависимости от относительного расхода пара в КС. Для рассмотренных выше вариантов видим, что, чем больше температура газопаровой смеси за КС и степень сжатия в ВК, тем более круто увеличивается коэффициент полезной работы, но при этом его максимальное значение меньше, чем у вариантов с меньшей температурой газов за КС и степенью сжатия в ВК. Так, у варианта 1 максимальный коэффициент полезной работы равен 0,77, а у варианта 4 – 0,69.
На рис. 7 показана зависимость температуры газопаровой смеси на выходе ГТ от d. С увеличением d температура газопаровой смеси на выходе ГТ линейно возрастает. Для варианта 1 от 600 до 665 °С, для варианта 2 от 595 до 660 °С, для варианта 3 от 620 до 680 °С и для варианта 4 от 630 до 700 °С.
Зависимость полезной работы на валу ГТ от d приведена на рис. 8. С увеличением d полезная работа на валу ГТ возрастает по ветви параболы. Чем больше температура газопаровой смеси за КС и степень сжатия в ВК, тем круче ветвь параболы, но меньше диапазон изменения d. Для варианта 1 полезная работа на валу ГТ возрастает в допустимом диапазоне изменения d в 2 раза – от 450 до 900 кДж/кг; для варианта 2 в 1,9 раз – от 510 до 975 кДж/кг; для варианта 3 в 1,8 раз – от 580 до 1045 кДж/кг; для варианта 4 в 1,7 раз – от 655 до 1130 кДж/кг.
На рис. 9 приведена зависимость расхода газопаровой смеси из КС на ГТ от d. Видим, что с увеличением d расход газопаровой смеси на ГТ падает по гиперболе. При этом для варианта 1 с изменением d от 0 до 12 расход газопаровой смеси на ГТ падает в 1,86 раз – от 220 до 118 кг/с; для варианта 2 в 1,75 раз – от 193 до 110 кг/с; для варианта 3 в 1,66 раз – от 170 до 102 кг/с; для варианта 4 в 1,59 раз – от 153 до 96 кг/с.
Направление дальнейших исследований
Исследования показали, что в КГПУ с вводом пара в КС можно получить электрический КПД до 55 %. При этом не рассматривался вопрос способа получения вводимого в КС пара. Из рис. 6 следует, что температура уходящей из ГТ газопаровой смеси составляет 595…700 °С, что позволяет использовать ее в котле-утилизаторе для получения пара. Этот высокотемпературный пар можно непосредственно вводить в КС или, для повышения эффекта охлаждения КС и ГТ, вводить через паровую турбину с выходом из нее на ввод пара в КС, это позволит выработать на ней дополнительную электрическую мощность. Необходимо определить, что более эффективно. Также необходимо рассмотреть вопрос снижения потерь теплоты и пара в окружающую среду путем установки конденсационного теплоутилизатора на выходе котла-утилизатора, который позволит понизить температуру газопаровой смеси до 40…50 °С и сконденсировать из нее большую часть влаги. При этом низкотемпературная тепловая энергия уходящих газов может быть использована для отпуска теплоты потребителям [28] или в органическом цикле Ренкина для дополнительной выработки электроэнергии [29] с применением воздушных конденсаторов [30].
Выводы
1. Исследования показали, что:
- электрический коэффициент полезного действия контактной газопаровой установки линейно возрастает с увеличением d. Диапазон изменения d сильно зависит от температуры газопаровой смеси за камерой сгорания и степени сжатия в воздушном компрессоре, чем они меньше, тем больше диапазон изменения d. Максимальный коэффициент полезного действия на уровне 56 % для всех вариантов достигается при максимальном d;
- коэффициент избытка воздуха в зависимости от d линейно падает, при этом чем выше температура газопаровой смеси за камерой сгорания и степень сжатия в воздушном компрессоре, тем больше темп падения и меньше диапазон изменения d;
- коэффициент полезной работы сильно зависит от относительного расхода пара в камеру сгорания, температуры газопаровой смеси за ней и степени сжатия в воздушном компрессоре, с ростом этих параметров он линейно возрастает;
- температура уходящих газов из газовой турбины также сильно зависит от d. С увеличением d эта температура линейно возрастает от 600 до 700 °С, при этом чем выше температура газопаровой смеси на выходе камеры сгорания и степень сжатия в воздушном компрессоре, тем больше температура газов на выходе газовой турбины;
- полезная работа на валу газовой турбины с увеличением d возрастает по ветви параболы. Чем больше температура газопаровой смеси за камерой сгорания и степень сжатия в компрессоре, тем круче ветвь параболы, но меньше диапазон изменения d;
- расход газов из камеры сгорания на газовую турбину сильно зависит от d, с увеличением d расход газов падает по гиперболе, чем ниже температура газопаровой смеси на выходе камеры сгорания и степень сжатия в воздушном компрессоре, тем больше снижение расхода газа на газовую турбину.
2. Повысить эффективность контактной газопаровой установки с вводом пара в камеру сгорания можно установкой на выходе газовой турбины котла-утилизатора для подготовки вводимого пара и конденсационного теплоутилизатором на его выходе для утилизации теплоты и влаги из уходящих газов. При этом низкотемпературная тепловая энергия уходящих газов может быть использована для отпуска теплоты потребителям или в органическом цикле Ренкина для дополнительной выработки электроэнергии.
About the authors
Nikolay N. Galashov
National Research Tomsk Polytechnic University
Author for correspondence.
Email: gal@tpu.ru
ORCID iD: 0009-0005-5351-3584
Cand. Sc., Associate Professor
Russian Federation, TomskAlexander A. Tubolev
National Research Tomsk Polytechnic University
Email: Tubolev@tpu.ru
Cand. Sc., Associate Professor
Russian Federation, TomskEvgeny S. Boldushevsky
National Research Tomsk Polytechnic University
Email: franky575@rambler.ru
Postgraduate Student
Russian Federation, TomskAlexander A. Minor
Tomsk Generation JSC
Email: aam29061994@gmail.com
Chief Expert
Russian Federation, TomskReferences
- Olkhovskiy G.G. Сombined cycle power plants in domestic thermal power engineering. Elektricheskie stantsii, 2020, vol. 1, pp. 21–28. (In Russ.)
- Manushin E.A. Modern high-capacity gas-turbine power-generating and combined-cycle plants of the world companies. Gazoturbinnye tekhnologii, 2020, vol. 3, pp. 2–8. (In Russ.)
- Olkhovskii G.G. The most powerful power-generating GTUS (review). Thermal Engineering, 2021, vol. 68, pp. 490–495. (In Russ.)
- Zysin V.A, Turchaninov B.V. Operation of gas-turbine combined-cycle plants in steam-and-gas cycle with a waste heat boiler. Energomashinostroenie, 1960, vol. 9, pp. 18–21. (In Russ.)
- Zysin V.A. Combined and combined-cycle gas installations and cycles. Moscow, Leningrad, GEI Publ., 1962. 186 p. (In Russ.)
- Arsenyev L.V., Tyryshkin V.G. Combined installations with gas turbines. St Petersburg, Mashinostroenie Publ., 1982. 247 p. (In Russ.)
- Kolp D.A., Meller D.Zh. Commissioning of the world's first full-cycle gas turbine STIG based on the LM 5000 gas generator (Simpson Paper Company). Sovremennoe mashinostroyenie. Seriya: A, 1989, no. 11, pp. 1–14. (In Russ.)
- Galashov N.N., Tupolev A.A., Minor A.A., Boldushevskiy E.S. Influence of the temperature of steam injection into the combustion chamber of a gas-steam plant on its energy characteristics. Bulletin of the Tomsk Polytechnic University. Geo Assets Engineering, 2023, vol. 334, no. 5, pp. 27–36. (In Russ.)
- Yuefei Xiong, Jiang Qin, Kunlin Cheng, Youyin Wang. Influence of water injection on performance of scramjet engine. Energy, 2020, vol. 201, pp. 117477–117490.
- Abubaker A.M., Darwish Ahmad A., Magableh M.N.A., Najjar Y.S.H. Efficiency boosting and steam saving for a steam-injected gas turbine engine: optimization study of the running conditions. Journal of Energy Engineering – ASCE, 2021, vol. 147 (1), pp. 732–748.
- Van Der Spek M., Bonalumi D., Manzolini G., Ramirez A., Faaij A.P.C. Techno-economic comparison of combined cycle gas turbines with advanced membrane configuration and MEA solvent at part load conditions. Energy and Fuels, 2018, vol. 32 (1), pp. 625–645.
- Stathopoulos P., Rähse T., Vinkeloe J., Djordjevic N. Steam injected Humphrey cycle for gas turbines with pressure gain combustion. Energy, 2019, vol. 188, p. 116020.
- Ziółkowski P., Kowalczyk T., Lemański M., Badur J. On energy, exergy, and environmental aspects of a combined gas-steam cycle for heat and power generation undergoing a process of retrofitting by steam injection. Energy ConverSI with Con and Management, 2019, vol. 192, pp. 374–384.
- Chmielewski M., Niszczota P., Gieras M. Combustion efficiency of fuel-water emulsion in a small gas turbine. Energy, 2020, vol. 211, pp. 118961–118985.
- Ivanov A.A., Ermakov A.N., Shlyakhov R.A. On the deep suppression of NOx and CO emissions in gas turbines with water or steam injection. Izvestiya RAN. Energetika, 2010, vol. 3, pp. 119–128. (In Russ.)
- Gordin K.A., Maslennikov V.M., Filimonova E.A. Assessment of the emission level of nitrogen oxides when steam with natural gas is supplied to the combustion chamber of a gas turbine installation. Thermophysics of high temperatures, 2013, vol. 51, no. 6, pp. 937–944. (In Russ.)
- Datsenko V.V., Zeygarnik Yu.A., Kosoy A.S. Experience in the use of water and steam to ensure environmental standards in conversion gas turbine engines. Teploenergetika, 2014, vol. 4, pp. 49–56. (In Russ.)
- Zaiguo Fu, Huanhuan Gao, Zhuoxiong Zeng, Jiang Liu, Qunzhi Zhu. Generation characteristics of thermal NOx in a double-swirler annular combustor under various inlet conditions. Energy, 2020, vol. 200, pp. 117487–117501.
- Farokhipour A., Hamidpour E., Amani E. A numerical study of NOx reduction by water spray injection in gas turbine combustion chambers. Fuel, 2018, vol. 212, pp. 173–186.
- Galashov N., Tsibulskiy S., Melnikov D., Kiselev A., Gabdullina A. Efficiency of utilization of heat of moisture from exhaust gases of heat HRSG of CCGT. MATEC Web of Conferences. Tomsk, 2017. рр. 01027–01031.
- Mokhtari H., Ahmadisedigh H., Ameri M. The optimal design and 4E analysis of double pressure HRSG utilizing steam injection for Damavand power. Energy, 2017, vol. 118, pp. 399–413.
- Gabdullina A.I., Galashov N.N., Tsibulskiy S.A., Melnikov D.V., Asanov I.A., Kiselev A.S. Promising direction of perfection of the utilization combine cycle gas turbine units. MATEC Web of Conferences. Tomsk, 2016. pp. 01004–01008.
- Machácˇková A., Kocich R., Bojko M., Kuncˇická L., Polko K. Numerical and experimental investigation of flue gases heat recovery via condensing heat exchanger. International Journal of Heat and Mass Transfer, 2018, vol. 124. pp. 1321–1333.
- Kostyuk A.G., Frolov V.V., Bulkin A.E., Truhnij A.D. Steam and gas turbines for power plants. Moscow, MEI Publ. house, 2008. 561 p. (In Russ.)
- Galashov N.N., Tupolev A.A., Minor A.A., Bannova A.I. Parametric analysis of the gas-steam installation scheme using a mathematical model. Bulletin of the Tomsk Polytechnic University. Geo Assets Engineering, 2021, vol. 332, no. 12, pp. 124–135. (In Russ.)
- Bell Ian H., Wronski J., Quoilin S., Lemort V. Pure and pseudo-pure fluid thermophysical property evaluation and the open-source thermophysical property library CoolProp. Industrial & Engineering Chemistry Research, 2014, vol. 53, no. 6, pp. 2498–2508.
- Thermal calculation of boilers (standard method). 3rd ed. St Petersburg, NPO CKTI Publ., 1998. 257 p.
- Galashov N.N., Tupolev A.A., Bespalov V.V., Minor A.A., Boldushevskiy E.S. Calculation of the parameters of the scheme of a gas-steam plant with deep utilization and heat release. Bulletin of the Tomsk Polytechnic University. Geo Assets Engineering, 2022, vol. 333, no. 12, pp. 124–135. (In Russ.)
- Galashov N.N., Tsibulsky S.A. Analysis of the efficiency of steam-gas plants of the trinary type. Bulletin of the Tomsk Polytechnic University, 2014, vol. 325, no. 4, pp. 33–38. (In Russ.)
- Galashov N.N., Tsibul’ski S.A. Numerical analysis of the condensation characteristics of different heat-transfer media in an air-cooled condenser. Power Technology and Engineering, 2016, vol. 49, no. 5, pp. 365–370.
Supplementary files
