Научная статья на тему 'Оптимальная степень сжатия воздуха в газотурбинной установке с промежуточным воздухоохладителем'

Оптимальная степень сжатия воздуха в газотурбинной установке с промежуточным воздухоохладителем Текст научной статьи по специальности «Химические технологии»

CC BY
1430
77
i Надоели баннеры? Вы всегда можете отключить рекламу.
Ключевые слова
ГАЗОТУРБИННАЯ УСТАНОВКА / GAS TURBINE PLANT / РАБОЧЕЕ ТЕЛО / WORKING MEDIUM / РЕАЛЬНЫЙ ГАЗ / REAL GAS / СТЕПЕНЬ СЖАТИЯ / COMPRESSION RATIO / ОПТИМИЗАЦИЯ / OPTIMIZATION / ВОЗДУХООХЛАДИТЕЛЬ / AIR COOLER / ТЕРМИЧЕСКИЙ КПД / COEFFICIENT OF THERMAL EFFICIENCY / ПОЛЕЗНАЯ МОЩНОСТЬ / USEFUL POWER / ТЕПЛОЁМКОСТЬ / HEAT CAPACITY

Аннотация научной статьи по химическим технологиям, автор научной работы — Хасанов Нариман Гаязович, Шигапов Айрат Багаутдинович

Разработана модель расчёта необратимых процессов расширения и сжатия реального газа. Приведён алгоритм для научно обоснованного выбора изоэнтропных КПД каскадов двухкаскадного компрессора при тепловом расчёте ГТУ. Доказана необходимость применения модели реального газа при поиске оптимальной степени сжатия воздуха в компрессоре по критериям максимальных термического КПД и полезной мощности.

i Надоели баннеры? Вы всегда можете отключить рекламу.

Похожие темы научных работ по химическим технологиям , автор научной работы — Хасанов Нариман Гаязович, Шигапов Айрат Багаутдинович

iНе можете найти то, что вам нужно? Попробуйте сервис подбора литературы.
i Надоели баннеры? Вы всегда можете отключить рекламу.

THE OPTIMUM OVERALL PRESSURE RATIO OF THE INTERCOOLED GAS TURBINE UNIT

In the article developed a mathematical model expansion and contraction of a real gas. Article shows the relationship between the isentropic efficiency of the compressor and the isentropic efficiency of the compressor stages. It was proved the need for a real gas model to find the optimal degree of compressed air in the compressor.

Текст научной работы на тему «Оптимальная степень сжатия воздуха в газотурбинной установке с промежуточным воздухоохладителем»

© Н.Г. Хасанов, А.Б. Шигапов УДК 621.438.082.2

ОПТИМАЛЬНАЯ СТЕПЕНЬ СЖАТИЯ ВОЗДУХА В ГАЗОТУРБИННОЙ УСТАНОВКЕ С ПРОМЕЖУТОЧНЫМ ВОЗДУХООХЛАДИТЕЛЕМ

Н.Г. Хасанов, А.Б. Шигапов

Казанский государственный энергетический университет, Россия

hasanovnariman @mail. ru

Резюме: разработана модель расчёта необратимых процессов расширения и сжатия реального газа. Приведён алгоритм для научно обоснованного выбора изоэнтропных КПД каскадов двухкаскадного компрессора при тепловом расчёте ГТУ. Доказана необходимость применения модели реального газа при поиске оптимальной степени сжатия воздуха в компрессоре по критериям максимальных термического КПД и полезной мощности.

Ключевые слова: газотурбинная установка, рабочее тело, реальный газ, степень сжатия, оптимизация, воздухоохладитель, термический КПД, полезная мощность, теплоёмкость.

THE OPTIMUM OVERALL PRESSURE RATIO OF THE INTERCOOLED GAS

TURBINE UNIT

N. G. Khasanov, A.B. Shigapov

Kazan State Power Engineering University, Russia

hasanovnariman@mail.ru

Abstract: In the article developed a mathematical model expansion and contraction of a real gas. Article shows the relationship between the isentropic efficiency of the compressor and the isentropic efficiency of the compressor stages. It was proved the need for a real gas model to find the optimal degree of compressed air in the compressor.

Keywords: gas turbine plant, working medium, real gas, compression ratio, optimization, air cooler, coefficient of thermal efficiency, useful power, heat capacity

Введение

Газотурбинные установки (ГТУ) с промежуточным охладителем (ПО) воздуха, рис. 1, находят применение в промышленности и энергетике [1]. Воздух в компрессоре сжимается при относительно низких температурах и высоких давлениях (например, степень сжатия пк = p2 /p в ГТУ GE ZMS100-42:1 [2]). Более низкие, чем в установке простого цикла, конечные температуры сжатия обуславливают высокое содержание горючего в продуктах сгорания в турбине. Таким образом, рабочее тело подобной установки не может быть достоверно рассчитано по законам идеального газа [3]. Основой теплового расчёта ГТУ, проводимого на ранних этапах проектных изысканий, является поиск оптимальной степени

сжатия воздуха в компрессоре опт по критериям максимальных полезной мощности Ьп :

Ьп = Ьт - Ьк,

где Ьт , - мощность газовой турбины и потребная мощность на привод компрессоров, МВт, и термического КПД цикла

т = Ьп 1 бкс ,

где Qкс - количество теплоты в камере сгорания МВт.

по

Рис. 1. Тепловая схема ГТУ с ПО: КНД, КВД - компрессоры низкого и высокого давления;

КС - камера сгорания; ГТ - газовая турбина; ПО - промежуточный охладитель

Во всех актуальных статьях поиск якопт осуществлён в модели идеального газа,

пример [4-9]. Цель настоящей работы - провести поиск як опт ГТУ с ПО в зависимости от

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

1. Модель процессов сжатия и расширения реального газа

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

- влияние давления на энтальпию рабочего тела в уравнении изоэнтропы (адиабаты):

врёТ = Т(дуIдТ)рёр, (для идеального газа - СрёТ = ТЯ(йрIр)); (1)

- влияние давления на теплоёмкость вр, показатель изоэнтропы к, коэффициент сжимаемости рабочего тела г : ср, к, 2 = /(Т; р);

- переменность и нелинейность зависимости вр, к, 2 рабочего тела от температуры;

- отличие уравнения состояния рабочего тела от уравнения идеального газа ру = ЯТ.

Анализ литературы показывает отсутствие аналитических зависимостей, комплексно

учитывающих проявления неидеальности газа, так как требование простого интегрирования реального уравнения адиабаты идёт вразрез с необходимостью учёта переменности термодинамических свойств рабочего тела от температуры и давления [10].

Принято решение использовать в статье численный способ расчёта. Далее будет приводиться пример для процесса сжатия газа. Рабочий процесс разбивается на п малых участков по давлению таким образом, чтобы пренебречь внутри участка изменением теплоёмкости и коэффициента сжимаемости (рис. 2), которые определяются при начальных параметрах участка Тш- и Р^, г = 1...п. Уравнение процесса, определяющее конечную температуру Ткг-, выведено из эмпирического уравнения состояния

pv = zRT

и уравнения адиабаты реального газа. Уравнения для отдельного участка имеют вид:

(Т, / Ttf) = (рк, / Рн," 'dT'PT)R''» ,

уравнение изоэнтропной работы реального газа при Пр, z = const из формул (1) и (2):

^из, =с РТн [(Ла /prn) Р Р •

рн Изменение давления Р/Рн рк

Рис. 2. Схема разбиения процесса сжатия на малые участки по давлению

Граничные условия всех участков совмещены по температуре и давлению -ТкI = Тн(1+1), ркI = рн(1+1), степень сжатия распределена равномерно: рк г- = рнг- ■ПуЧ,

степень сжатия в участке ЛуЧ = • Исходя из совмещения граничных условий и

формулы (2), Тнг и Рнг участков меняются, поэтому сумма определённых, при постоянных

теплоёмкостях, работ участков даст уже общую работу при переменной теплоёмкости

Еп

¡-Х^тг • Число п выбирается избыточно: п «1000. Таким образом, проявления

неидеальности воздуха учитываются комплексно.

Представленные в литературе аппроксимационные уравнения зависимостей термодинамических свойств рабочих тел от температуры и давления выполнены в узком диапазоне температур, соответствующих рабочим режимам установки простого цикла. Разработаны уравнения Ср,к,2 = /(Т;р) для воздуха в интервалах ДТ =220...1190 °К,

Ар =0,1.20 МПа. Построение происходило по табулированным данным национального института стандартов США [11]. Максимальная относительная погрешность 8 <2,2%. Построены зависимости для природного газа осреднённого состава [12]: Ср = /(Т; а),

к = /(Т;а), 2 = /(Т;а), где а - коэффициент избытка воздуха. ДТ =350.1850 °К, Да =1.10. Относительная погрешность <0,9%. В качестве модели камеры сгорания разработано аппроксимационное уравнение зависимости температуры на входе в газовую турбину Т3 от температуры нагретого в компрессоре воздуха Т2 (ДТ2 =400.1200 °К) и а (Да =1.9). Исходные данные получены в программном комплексе TDSOFT [13]. Для аппроксимации функции двух переменных применялся численный способ, а также метод множественной нелинейной регрессии в программном комплексе STATISTICA. 2. Учёт необратимости процессов сжатия и расширения

При тепловом расчёте потери энергии вследствие несовершенства проточной части агрегата учитываются изоэнтропным (внутренним) КПД:

"Лиз = ^из 1 ^ ,

- отношение изоэнтропной работы сжатия к фактической. Известно, что в компрессоре сумма изоэнтропных работ ступеней больше общей изоэнтропной работы компрессора [14]

(в турбине - наоборот). Элементарный участок можно считать ступенью компрессора, поэтому в рамках разработанной модели расчёта не выполняется условие

ш

Ь = [Хы ¿из/ )1Л из . (3)

Учёт необратимости необходимо вести внутри элементарного участка процесса:

Ьг = ¿из/1 Луч = сргТнг[(рк/ /рн/) Р Р -1] 1 Лиз.уч ,

где Лиз.уч - изоэнтропный КПД участка процесса, его можно распределить равномерно между участками. Тогда величину КПД участка можно вынести из под знака суммы Ь = (¿из/1Лиз.уч), уравнение (3) запишется правильно:

Ь = ^^/=1 ¿из/ ^1Л из.уч.

Это значит, стоит задача связать общий изоэнтропный КПД процесса и изоэнтропный КПД его участков. Луч и Луч принимаются постоянными. Часто

встречается ошибочное мнение, что произведение изоэнтропных КПД ступеней определяет общий изоэнтропный КПД агрегата, Лиз.уч = ПЛиз . В таком случае, с ростом числа участков п Лиз.уч должен сходиться к единице, но на самом деле Лиз.уч сходится к величине так называемого политропного КПД Лпол при п ( лш <Лпол < 1) [15]. По выкладкам [14], а так же формулам газовой динамики выведена связь Лиз.уч и Лиз :

_ пек-1

Лизуч = ^ [(%-1)1 Лиз ] +1 -1, (4)

, к

Ср/к

ек , ср

где ек =пкр , ср - средняя эффективная теплоёмкость. Так как перед началом расчёта

перепад действительных температур не известен, Ср допустимо рассчитать при средней температуре изоэнтропного процесса. Для процессов расширения

_1 - [1 -Лиз (1 -11 ет XI11 п

Лизуч 1 -1/п/Г ■

Проверка формул (4), (5) осуществляется следующим образом: при расчёте компрессора Лизуч подбирается вручную так, чтобы отношение общей изоэнтропной

работы к сумме действительных работ в участках равнялось Лиз . Подобранное вручную значение сравнивается с расчётом по (4), должно обеспечиваться равенство.

3. Математическая модель теплового расчёта ГТУ с ПО

Необходимо выбрать усреднённые, наиболее характерные и приближённые к действительности значения величины недоохлаждения воздуха тк = IТ (см. рис. 1) и относительного аэродинамического сопротивления в ПО ст = (1-р^/р21)'100, которые войдут в математическую модель постоянными. На основании анализа рекомендаций по тепловому расчёту ГТУ обширного круга авторов, показателей действующих установок в диапазоне мощностей 5-100 МВт, математического моделирования воздухоохладителя ГТ-100-750 в диапазоне степеней сжатия предвключённого компрессора 2...24 МПа, принято ст уравнять 3% от входного давления. Принято допущение об отсутствии недоохлаждения относительно температуры окружающего воздуха (тк =1).

Перед расчётом естественно задать в исходных данных общий изоэнтропный КПД компрессора Лиз . При этом в литературе отсутствуют указания и алгоритмы, как научно

обосновать выбор изоэнтропных КПД каскадов компрессора Ли1 и Лиз . То есть, как связаны изоэнтропные КПД компрессора в целом, и части процесса с определённой степенью сжатия и начальными условиями. Если выбрать Ли1 и Лиз произвольно или неправильно, нарушится энергетический баланс и изменится изначально заданный лиз.

Предложен алгоритм:

1) задаётся количество элементарных участков сжатия в компрессоре п , и по формуле (4) вычисляется изоэнтропный КПД элементарного участка процесса Лиз.уч при

общей степени сжатия пк ;

2) определяется степень сжатия в элементарной ступени Яуч - постоянная величина;

3)определяется число элементарных ступеней, укладывающихся в первый каскад компрессора: п1 = пк1-*;

4) решается обратная п.1 задача: по уже известному КПД элементарного участка процесса, степени сжатия в первом каскаде лк1) и числу ступеней п1 первого каскада определяется изоэнтропный КПД первого каскада по формуле

с<1)/Я П1

ли) =-К1)Р /[(1+

сW/R -К1) Р -1 n1

)n1 -1];

"Лиз.уч

5) определяется число элементарных участков второго каскада n2 = n - n1;

6) обратная п. 1 задача решается повторно по як2), n2, для определения Лиз . Проверка достоверности алгоритма: двухкаскадный компрессор рассчитывается при

отсутствии гидравлического сопротивления в теплообменнике и без охлаждения. При

правильно распределённых Лиз и Лиз потребная мощность на привод такого компрессора не должна отличаться от мощности простого однокаскадного компрессора с той же як, Лиз, и начальными параметрами воздуха.

4. Поиск оптимальной степени сжатия пк опт в ГТУ с ПО

Решение данной задачи для установки простой схемы ГТУ широко представлено в литературе как аналитическими, так и численными методами. Для установки с ПО задача усложняется: на величину "t и Lj влияет не только як, но и выбор степени сжатия в

первом каскаде компрессора . В силу большого количества рекуррентных зависимостей, даже аналитическое решение данной проблемы требует вычисления 3-4 циклов итерации. Применяются упрощённые зависимости идеального газа. Усложнение расчётных формул невозможно. Численным способом в модели полуидеального газа

решена только задача выбора я^опт при як = const [15].

Разработанная авторами данной статьи математическая модель на языке FORTRAN с применением библиотеки численных методов IMSL решает задачу методом перебора -анализирует с определённым шагом все совокупности як и tc^-* , Дяк =5.200, и выбирает те, которые обеспечивают максимум "t и Ln. Все дискретные зависимости аппроксимируются. Целевая функция представляет собой зависимость "t, рис. 3, и L от

iНе можете найти то, что вам нужно? Попробуйте сервис подбора литературы.

як при условии, что в рамках каждого постоянного як, як1* определён оптимально. Проведён поиск як опт в двух предположениях:

а) по уравнениям идеального газа, учитывая различие теплоемкостей воздуха и продуктов горения, определяя работу газа и конечную температуру процесса по теплоёмкости, найденной при средней температуре Т процесса. Учитывается методом последовательных приближений рекуррентная зависимость Ср = /(Т). Такой подход

является классическим при определении технико-экономических показателей ГТУ;

б) при учете неидеальных свойств газа, описанных в разделе 1 данной статьи.

Рис. 3. Зависимость термического КПД ГТУ и оптимальной степени сжатия в первом каскаде от общей степени сжатия воздуха. Начальная температура продуктов сгорания 1200 К

Исходные данные расчёта: расход воздуха Св=300 кг/с; Т1 =290 К ; р1 =0,1 МПа; Т12 =290 К ; ст =3%; Т3 =1100 * 1600 К ; Лиз.к =0,87; изоэнтопный КПД турбины Лиз.т =0 89. Часть переменных варьируется для изучения их влияния на положение лк опт и лкопт.

В результате расчёта наблюдаются высокие лк опт, пример для модели реального газа - в табл. 1. Это обусловлено низкой потребной работой компрессора Ьк в ГТУ с ПО, так как момент, когда скорость изменения роста Ьк с увеличением ^к превысит скорость роста работы турбины Ьт, запаздывает. Высокие лк опт по обеспечены эффективностью сжатия с высокой конечной температурой по затратам топлива в камере сгорания.

Таблица 1

Зависимость оптимальной степени сжатия от начальной температуры продуктов сгорания

Тз, К 1100 1200 1300 1400 1500 1600

_ /тшахч дкопт(¿п ) 16 20 25 31 38 47

„„ /„шах-, я копт(Л t ) 31 38 50 68 87 120

Зависимости Л1 = /(лк) , пример на рис. 3, и Ьп = /(%к) имеют плато в районе оптимума ^к (область, где значительное изменение ^к незначительно меняет Ьп и л? ). Несмотря на небольшое уточнение (в пределах 5% относительных) величин Ьи и л , за счёт учёта реальных свойств газа, ^копт меняется существенно. Данная особенность проявляется при всех начальных температурах продуктов сгорания Т3 . Идеальногазовая модель даёт завышенные значения ^копт (табл. 2). Причиной разницы, в первую очередь, служит учёт влияния давления на с„ воздуха - с ростом Лк, в модели реального газа, ср растёт быстрее не

только под влиянием растущей температуры, но и растущего давления, соответственно, быстрее растёт работа сжатия ср (Т2 - Т) • Разница в рассчитанных лкопт ощутимо меняет

эксплуатационные свойства компрессора (металлоёмкость, потребное количество охлаждающей воды, надёжность и др.). Очевидно, поиск як опт ГТУ с ПО необходимо проводить в модели

реального газа. Идеальногазовая модель кроме прочего даёт завышенные значения и

^тах . Наличие плато на графике ставит вопрос о возможности снижения як от оптимального значения без потери ^ и ¿п, как следствие, целесообразности применения регенерации совместно с ПО, включения ГТУ с ПО в парогазовую установку. На положение як опт, и

соответствующее л^, оказывают влияние разные факторы, качественная сторона тенденций схематично приведена в табл. 3.

Таблица 2

Абсолютная ошибка определения як / относительная погрешность по ¿Л^ и л™аХ при

применении модели идеального газа (в отличие от реального) в зависимости от характерных

параметров ГТУ

Оптимизация по мощности Оптимизация по термическому КПД

т} 1100 1200 1300 1400 1500 1600 1100 1200 1300 1400 1500 1600

Г] 223 6/3,2 10/3,4 17/3,7 28/3,9 34/4,1 45/4,4 19/4,2 30/4,3 40/4,7 52/2,7 64/3,0 64/3,3

273 1/2,5 3/2,5 5/2,5 7/2,5 11/2,6 14/2,8 8/3,6 8/3,3 14/3,3 21/3,5 33/0,2 47/3,3

323 1/1,3 1/1,2 1/1,3 2/1,3 3/1,4 5/1,5 2/2,1 3/1,9 3/2,0 6/2,1 7/2,3 16/2,6

тк I 1/1,67 2/1,2 3/1,7 4/1,8 7/1,9 10/2,1 4/2,6 7/2,4 10/2,5 13/2,7 24/3,0 39/3,3

1,1 1/1,4 1/1,4 2/1,4 3/1,5 5/1,6 8/1,7 5/2,4 5/2,2 7/2,2 7/2,3 15/2,4 23/2,7

1,2 0/1,1 1/1,2 2/1,2 2/1,2 3/1,3 4/1,3 5/1,7 10/1,8 6/1,9 6/2,0 8/2,1 18/2,2

<7 6 2/1,7 2/1,7 3/1,7 4/1,8 7/1,9 11/2,1 3/2,6 6/2,4 8/2,5 12/2,7 19/2,9 32/3,1

9 1/1,7 2/1,7 3/1,7 4/1,8 7/1,9 10/2,1 4/2,6 4/2,4 7/2,5 12/2,7 19/2,9 29/3,1

Пт.к 0,83 1/1,7 2/1,7 2/1,7 4/1,8 6/1,9 8/2,0 5/2,4 5/2,3 7/2,5 11/2,6 14/2,8 31/3,2

0,в9 2/1,7 2/1,7 3/1,8 6/1,9 8/2,0 12/2,1 4/2,7 7/2,5 12/2,6 16/2,8 29/-4,1 38/-7.2

Таблица 3

Характер влияния переменных на оптимальную степень сжатия - общую як опт и в первом

каскаде якопт

\ Измене Рост\Ние фактора п . ттах 11о критерию ьп ; о „тах По критерию ц, :

к к опт л(1) А' опт тт ах 71 11,к опт л(1) "■к опт Т1 гх

Гт3 Тл к опт к опт ^ ттах 1 1- л к опт 1 л(1) к опт Т п,тах

и ^ опт X л(1) » 11 к опт 1 .тах 1 л 1 к опт - X п,тах

к ^ опт Тл(1) |1к опт 1 гтах X /-„ ^ Кк опт Т л(|) 1 к опт X пГх

Та X Л * к опт - X /тах УГ * ' У! X л * к опт Т л(1) 1 Ч" опт X г|,тах

1Ч,3 * к опт 1 л(1) 1 11 к опт ф г тах 1 к опт X Л(1) ж к опт Г л,тах

Положение якош. определяет баланс: по критерию ¿Л^ - потребной работы на сжатие и работы турбины; по критерию Лшах - ¿п

и расхода топлива От опосредованно через конечную температуру сжатия. Разница в як опт, определённых по моделям

реального и идеального газа (табл. 2), также является переменной величиной, так что увеличивается с падением температур и ростом давлений. В связи с влиянием конечной температуры сжатия на Лг актуальна задача поиска оптимальной степени охлаждения воздуха в ПО. Анализ дал отрицательный результат, Лг , оценённый при оптимальных параметрах сжатия, неуклонно снижается при росте тк (табл. 4).

Таблица 4

Влияние степени охлаждения в ПО т^ на термический КПД лг

ь 1000 1200 1300 1400 1500 1600

л, т*=1 0,402 0,44 0,47 0,5 0,52 0.54

тк=1,1 0,38 0.42 0.46 0.49 0,51 0,53

тк=1.2 0,37 0.41 0,45 0,48 0.504 0,52

В литературных источниках встречается мнение, что при одинаковых степенях сжатия термический КПД ГТУ с ПО ниже, чем у установки простого цикла. Проведено сравнение величин и Лг для этих двух схем при оптимально рассчитанных степенях сжатия и прочих равных условиях. Благодаря применению описанной в разделе 3 модели распределения изоэнтропных КПД каскадов компрессора, Лизк остаётся постоянно

равным для обеих схем ГТУ, независимо от соотношения я^ / я^ в схеме с ПО. Предварительный расчёт як опт установки простого цикла показал необходимость применения модели реального газа при оптимизации по Лг. В интервале температур 1200 * 1600 °К мощность ГТУ с ПО выше мощности ГТУ простой схемы в среднем на 22,9%. Термический КПД ГТУ с ПО выше в среднем на 8% (относительных). При этом як опг ГТУ с ПО в обоих случаях значительно выше, разница растёт с ростом Т3.

Проведён анализ положения я^от-, соответствующей якопт. Изучено влияние

различных факторов на положение я^о^, результаты приведены в табл. 3. Для

оптимизации по Лг характерно высокотемпературное сжатие с низким значением я^ощ и изменением его в узком диапазоне (для снижения затрат на топливо в камере сгорания). При оптимизации по , я^ определяется так, чтобы обеспечить равенство потребных работ

каскадов. я^ и як2) примерно равны при Т2 = Т, при Т2 > Т степени сжатия я^ и я^ перераспределяются так, чтобы минимизировать добавку работы в высокотемпературной части зависимости ср = /(Т; р) (вследствие повышения средней температуры сжатия во

втором каскаде компрессора). Рассмотрен характер зависимостей ¿п = /(я^) и лг = / (як1))

при неоптимальных значениях як . я^ всегда, не зависимо от як, имеет оптимум по обоим

критериям расчёта между двумя предельными неэффективными случаями: я^ =1 (т.е. ГТУ

простой схемы) и я^ = як . Это опровергает тезис о том, что при як =сош1 и прочих равных условиях, эффективность цикла с ПО ниже, чем у установки простого цикла. На положение

лк опт влияет баланс двух основных факторов: затраты работы на привод компрессора и конечной температуры сжатия.

Выводы

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

2) Усовершенствованная численная модель расчёта по малым интервалам давления учитывает совокупно все проявления неидеальности рабочего тела, а также диссипацию энергии посредством изоэнтропного КПД агрегатов или показателя политропы.

3) Методика распределения изоэнтропных КПД каскадов компрессора пропорционально их степеням сжатия, термодинамическим параметрам, общей степени сжатия и заранее заданному общему изоэнтропному КПД теоретизирует задачу выбора изоэнтропных КПД каскадов компрессора ГТУ с ПО при тепловом расчёте.

4) Литературный анализ аналитических и численных моделей оптимизации степени

сжатия в ГТУ с ПО по критериям , Lr^ показал их некорректность. Разработанная численная модель позволяет строить зависимости Lп, % = f (як) при условии одновременно

оптимально распределённых и л42) .

5) Результаты расчётов показали высокие значения оптимальных степеней сжатия по Li1 ax и особенно по в ГТУ с ПО. Поиск якопт такой схемы должен проводиться

только в модели реального газа. Также для простой схемы ГТУ оптимизация по ^max должна проводиться с учётом неидеальных свойств. При оптимальных як и одинаковых T3 ГТУ с ПО имеют существенно более высокую экономичность и полезную мощность по сравнению с установкой простого цикла.

Литература

1. Saidi A. Intercoolers in gas turbine systems and combi-processes for production of electricity //ASME Turbo Expo 2000: Power for Land, Sea, and Air. - American Society of Mechanical Engineers, 2000. - 10 p.

2. Carl Knopf. Modeling, Analysis and Optimization of Process and Energy Systems. John Willey & Sons, New Jersey, 2012. 488p.

3. Манушин Э.А., Михальцев В.Е., Чернобровкин А.П. Теория и проектирование газотурбинных и комбинированных установок. М.: Машиностроение, 1977. 447 с.

4. Yang B. et al. Exergy performance analyses of an irreversible two-stage intercooled regenerative reheated closed Brayton CHP plant //International Journal of Exergy. — 2014. — I. 14. — №. 4. — p. 459-483.

5. Al-Sood M. M. A., Matrawy K. K., Abdel-Rahim Y. M. Optimum operating parameters of an irreversible gas turbine cycle. // Journal of Engineering Sciences, Assiut University, Vol. 40 -2014- № 6 p.1695-1714.

6. Sarath R. Numerical analysis for the prediction on the effect of heat transfer characteristics of combined cycle gas turbine using inter cooler // Interntional Journal of Engineering Research and Reviews. -2016. - Vol. 4, I. 2, p. 20-32.

7. Иванов В. А. Повышение эффективности стационарных и судовых газотурбинных установок // Вестник АГТУ. Серия: Морская техника и тех-нология. 2012. №2. C. 76-80.

8. Волощук В.А., Очков В.Ф., Орлов К.А. Термодинамическая оптимизация циклов некоторых ГТУ и ПГУ при помощи современных информационных технологий (Часть 1) //Новое в российской электроэнергетике. 2011. №. 7. С. 23-42.

9. Milind S. Patil at al. Thermal Performance of Reheat, Regenerative, Inter Cooled Gas Turbine Cycle // International Journal of Research in Mechanical Engineering & Technology. - 2015 - Vol. 5, Issue 2 p. 28-33.

10. Истомин В.А. Обобщенные показатели изоэнтропы реального газа. //Теплофизика высоких температур. 1998. Том 36, выпуск 5. С. 732-739.

11. P.J. Linstrom. NIST Chemistry WebBook, NIST Standard Reference Database Number 69, National Institute of Standards and Technology, Gaithersburg MD, 20899, http://webbook.nist.gov, (retrieved October 12, 2015).

12. Шигапов А.Б. Стационарные газотурбинные установки тепловых электрических станций. Казань: Казан. гос. энерг. ун-т, 2009. 416 с.

13. Назырова Р.Р. Исследование операций в оценке термодинамических характеристик. Казань: Изд. АБАК, 1999. 197 с.

14. Стечкин Б. С. и др. Теория реактивных двигателей. М.: Оборонгиз, 1956. 534 с.

15. Арсеньев Л.В., Тырышкин В.Г. Стационарные газотурбинные установки: справочник. Ленинград: Машиностроение, 1989. 543 с.

16. Гулина С.А., Орлова Г.М., Орлов М.Ю. Особенности конвертирования авиационного двигателя в газотурбинный привод центробежного нагнетателя для магистрального газопровода // Вестник. Самар. гос. техн. ун-та. Сер. технические науки. 2014. № 1 (41). С. 152-158.

Авторы публикации

Хасанов Нариман Гаязович - аспирант кафедры «Энергетическое машиностроение» (ЭМ) Казанского государственного энергетического университета (КГЭУ). E-mail: hasanovnariman@mail.ru .

Шигапов Айрат Багаутдинович - д.т.н., профессор кафедры «Энергетическое машиностроение» (ЭМ) Казанского государственного энергетического университета (КГЭУ). E-mail: shigapov38@mail.ru.

References

1. Saidi A. Intercoolers in gas turbine systems and combi-processes for production of electricity //ASME Turbo Expo 2000: Power for Land, Sea, and Air. - American Society of Mechanical Engineers, 2000. - 10 p.

2. Carl Knopf. Modeling, Analysis and Optimization of Process and Energy Systems. John Willey & Sons, New Jersey, 2012. 488p

3. Manushin Je.A., Mihal'cev V.E., Chernobrovkin A.P. Teorija i proektirovanie gazoturbinnyh i kombinirovannyh ustanovok. M.: Mashinostroenie, 1977. 447 p.

4. Yang B. et al. Exergy performance analyses of an irreversible two-stage in-tercooled regenerative reheated closed Brayton CHP plant //International Journal of Exergy. — 2014. — I. 14. — №. 4. — p. 459-483.

5. Al-Sood M. M. A., Matrawy K. K., Abdel-Rahim Y. M. Optimum operating parameters of an irreversible gas turbine cycle. // Journal of Engineering Sciences, Assiut University, Vol. 40 -2014- № 6 p.1695-1714.

6. Sarath R. Numerical analysis for the prediction on the effect of heat transfer characteristics of combined cycle gas turbine using inter cooler // Interntional Journal of Engineering Research and Reviews. -2016. - Vol. 4, I. 2, p. 20-32.

iНе можете найти то, что вам нужно? Попробуйте сервис подбора литературы.

7. Ivanov V. A. Povyshenie jeffektivnosti stacionarnyh i sudovyh ga-zoturbinnyh ustanovok // Vestnik AGTU. Serija: Morskaja tehnika i teh-nologija. 2012. №2. p. 76-80.

8. Voloshhuk V. A., Ochkov V. F., Orlov K. A. Termodinamicheskaja optimi-zacija ciklov nekotoryh GTU i PGU pri pomoshhi sovremennyh informacionnyh tehnologij (Chast' 1) //Novoe v rossijskoj jelektrojenergetike. — 2011. — №. 7. — p. 23-42.

9. Milind S. Patil at al. Thermal Performance of Reheat, Regenerative, Inter Cooled Gas Turbine Cycle // International Journal of Research in Mechanical Engineering & Technology. - 2015 - Vol. 5, Issue 2 p. 28-33.

10. Istomin V. A. Obobshhennye pokazateli izojentropy real'nogo gaza. Teplo-fizika vysokih temperatur, 1998, tom 36, vypusk 5, p. 732-739.

11. P.J. Linstrom. NIST Chemistry WebBook, NIST Standard Reference Database Number 69, National Institute of Standards and Technology, Gaithersburg MD, 20899, http://webbook.nist.gov, (retrieved October 12, 2015).

12. Shigapov A.B. Stacionarnye gazoturbinnye ustanovki teplovyh jelektricheskih stancij Kazan': Kazan. gos. jenerg. un-t, 2009., 416 p.

13. Nazyrova R.R. Issledovanie operacij v ocenke termodinamicheskih harakteristik. Kazan. ABAK. 1999. 197 p.

14. Stechkin B. S., Kazandzhan P. K., Alekseev L. P. Teorija reaktivnyh dvigatelej. M.: Oborongiz. 1956. 534 p.

15. Arsen'ev L.V., Tyryshkin V.G. Stacionarnye gazoturbinnye ustanovki: spravochnik. Leningrad: Mashinostroenie, 1989. 543 p.

16. Gulina S.A., Orlova G.M., Orlov M.Ju. Osobennosti konvertirovanija aviacionnogo dvigatelja v gazoturbinnyj privod centrobezhnogo nagnetatelja dlja magistral'nogo gazoprovoda. Vestn. samar. gos. tehn. unta. Ser. tehnicheskie nauki. 2014. № 1 (41), p. 152-158.

Authors of the publication

Khasanov Nariman Gayazovich - graduate student of department «Power Engineering» of Kazan State Power Engineering University.

Shigapov Ayrat Bagautdinovich- Doc. Sci. (Techn.), Assoc. Prof. of department «Power Engineering» of Kazan State Power Engineering University.

Поступила в редакцию 29 декабря 2016 г.

i Надоели баннеры? Вы всегда можете отключить рекламу.