The application of thermodynamic calculation for the monitoring of the characteristics of gas-distributing-unit gas-turbine drives

Full Text

В настоящее время газовое топливо широко используется в наземных газотурбинных установках. Для получения достоверных термодинамических и эксплуатационных характеристик ГПА необходим учет влияния состава сжигаемого топливного газа. Приближенные расчеты их характеристик обычно производятся с использованием математических моделей, предусматривающих постоянство компонентного состава природного газа. Для точного определения параметров цикла невозможно использовать модели процессов на основе молекулярно-кинетической теории идеальных газов. Как показал анализ используемых методов расчета характеристик ГТД [1, 5], их достоверность в расчетных моделях может быть обеспечена только при учете температуры атмосферного воздуха, фактического состава топливного газа и температуры при его сжигании. Необходимость учета переменности теплофизических свойств воздуха и продуктов сгорания в процессе термодинамических расчетов циклов достаточно очевидна, поскольку неучет в расчетах изменения свойств рабочего тела, его зависимости от температуры и величины коэффициента избытка воздуха α в камере сгорания дает ошибку в вычислениях более 5 % [5]. Расчеты параметров ГТД значительно усложняются при их работе на природном газе произвольного состава, зависящем от места его добычи. Как правило, используемый в них топливный газ не имеет постоянного химического состава: так, например, содержание в нем метана СН4 для различных месторождений изменяется от 80 до 99 % при переменном составе остальных компонентов - предельных углеводородов (этан, пропан, бутан, пентан), примесей азота, серы, углекислого газа, сероводорода и водорода. Соответственно этому изменяется состав продуктов сгорания, коэффициент избытка воздуха α, мощность Ne и экономичность двигателя. Особенность состава рабочего тела ГТД определяет его термодинамические свойства, зависящие от переменности теплоемкости воздуха и продуктов сгорания, в диапазоне рабочих значений давлений и температур. Для различных составов природного газа (в зависимости от месторождения) определены химико-физические данные компонентов газовой смеси и рассчитано теоретическое количество воздуха, необходимое для сгорания 1 м3 газа (стехиометрический коэффициент): . Вышеуказанные данные сведены в таблицу. Теплофизические показатели газового топлива Газовое топливо (месторождение) CH4 Др. углевод. CO CO2 N2 H2О H2 O2 , мДж/кг L0, м3/м3 Водород 100 % 111,667 2,380 Метан 100 % 46,602 9,520 Ямбург 98,2 1,12 0,67 0,01 46,183 9,523 Вуктылское 81,8 12,84 0,3 5,06 43,432 10,325 Оренбургское 84,6 5,05 1,3 9,05 41,638 9,080 Покровское 65,3 8,4 0,3 26,0 34,101 8,002 Сухой отбенз. газ 92,7 4,7 0,8 1,7 0,1 45,245 9,868 Биогаз 71,7 - 24,1 4,1 33,428 6,833 ПНГ 40 50 5,0 1 4 32,072 11,353 Доменный газ 0,1 - 30,1 10 55,1 2,3 2,2 5,335 0,778 В зависимости от месторождения природного газа за счет различного содержания примесей и тяжелых углеводородов теплота его продуктов сгорания изменяется в среднем на 27 %, а стехиометрический коэффициент L0 - на 22 %. Эти показатели определяют значение коэффициента избытка воздуха α для термодинамического цикла ГТД. Поэтому для проведения инженерных термодинамических расчетов ГТД оказалось необходимым и целесообразным создать программу определения теплофизических свойств рабочего тела, учитывающую характеристику и компонентный состав газа и продуктов его сгорания. В этой программе использованы широко применяемые на практике электронные таблицы Excel с включенными в них массивами данных для любого состава топливного газа. Они включают характеристики компонентов газового топлива (CmHn; Н2S; CO2; O2; CO; H2; H2O; N2), продуктов сгорания газообразного топлива (N2, CO2, O2, H2O) и воздуха. Теплофизические параметры продуктов сгорания и воздуха (- теплоемкость, - энтропия, - энтальпия) определялись по зависимостям, рекомендованными авторами [2], в виде аппроксимирующего полинома шестого порядка: ; ; , где αj - коэффициенты соответствующих полиномов, приводятся в [2]; 1000 - масштабный коэффициент. Для воздуха расчеты проводятся в диапазоне температур от -25 до 1650 °С, для продуктов сгорания - от 650 до 1650 °С. Теплофизические параметры компонентов смеси продуктов сгорания (энтальпия и логарифм относительного давления ) при их стехиометрическом составе (коэффициенте избытка воздуха α = 1) рассчитывается по формуле , ; , где r - объемная доля компонента, %. Дополнительно рассчитываются показатель изоэнтропы k(Т) и газовая постоянная R (Т). Для остальных значений коэффициента избытка воздуха α > 1 определяется объемная доля воздуха в продуктах сгорания rв с учетом влажности воздуха на входе в двигатель: , где V0 - суммарный объем продуктов сгорания; - влагосодержание. Зависимость объемной доли воздуха rв от коэффициента избытка воздуха при α от 1 до ∞ для природного газа месторождения Ямбург приведена на рис. 1. Рис. 1. Зависимость объемной доли воздуха rв от коэффициента избытка воздуха α (для природного газа месторождения Ямбург) В соответствии с рекомендациями [2] молекулярная масса продуктов сгорания для любого содержания воздуха в продуктах сгорания табулируется с помощью линейной интерполяции, значения энтальпии и относительного давления данной смеси определяется линейной интерполяцией между значениями при a = 1 и a = ∞, что соответствует rв = 0% и rв = 100% : ; . Все рассчитанные параметры для продуктов сгорания и воздуха табулируются в функции температуры с шагом в 1°С. Для рассчитанных термодинамических параметров смеси продуктов сгорания газового топлива применена равномерная шкала температур Т с интервалом в 1 °С и неравномерные шкалы энтальпий h(Т) и относительных давлений π(Т). Эти шкалы равноправны с точки зрения рассмотрения их в качестве аргумента или функции. В среде пакета Microsoft Excel с помощью функции «Просмотр» искомое значение (рассчитанный или заданный в ходе термодинамического расчета цикла параметр) определяется в массиве аргумента «Просматриваемый вектор». Форма массива позволяет просматривать искомое значение из списка рассчитанного параметра в строке или столбце указанного аргумента, находит указанное значение и возвращает значение из аналогичной позиции столбца или строки искомого параметра. По значению имеющегося аргумента, допустим рассчитанному значению энтальпии h, находится искомый параметр - температура в цикле Т или относительное давление , или наоборот. Погрешность определения искомого значения составляет 0,01 %. Используемая модель теплового процесса и разработанная программа позволяют рассчитывать теплофизические параметры рабочего тела для любого состава топливного газа и автоматически использовать их при проведении теплового расчета двигателя. Модель теплового процесса ГТД составлена на основе объединенного закона термодинамики: , где - изменение теплоты в цикле; - изменение энтропии в цикле при сохранении условия зависимости от температуры. Авторы [4] предлагают связь между температурой и давлением в изоэнтропическом процессе сжатия или расширения определять в виде отношения , где функция , где , - значение энтропии при соответствующих давлениях. Термодинамическая функция π(Т), являющаяся относительным давлением, зависит только от одного термодинамического параметра - температуры Т. Параметр относительного давления определяется для условий изоэнтропического процесса. Для изоэнтропического процесса и , где е - основание натурального логарифма; = 8,3142 - мольная газовая постоянная. Энтропия определяется: . (1) Для изоэнтропического процесса имеем: , (2) где - значение при температуре газа, соответствующей концу изоэнтропического процесса; - значение при температуре газа, соответствующей началу изоэнтропического процесса. Для реального процесса, протекающего при тех же давлениях P1 и P2 и той же начальной температуре, , (3) где - значение при температуре газа, соответствующей концу реального процесса. Вычитая из уравнения (2) уравнение (3), получаем , т. е. изменение энтропии в реальном процессе равно разности значений , соответствующих температурам газа в конце реального и изоэнтропического процессов (при одинаковом отношении давлений). Последнее равенство после подстановки значений из уравнения (1) примет вид , (4) т. е. изменение энтропии определяется через значения логарифма относительного давления. Изменение энтальпии с учетом зависимости ср от температуры Т определится: , (5) относительное давление: . Использование уточненных значений теплофизических параметров рабочего тела (энтальпий h (Т), относительного давления πis(Т)) позволяет выполнять термодинамический расчет цикла двигателей, работающих на топливном газе любого состава. Разработанная модель по сравнению с существующими моделями позволяет более детально учесть переменность теплофизических свойств рабочего тела в зависимости от процентного содержания основных компонентов и примесей. На базе разработанной математической модели для определения термодинамических показателей ГТД составлен алгоритм термодинамического расчета двигателя для газообразного топлива произвольного состава (рис. 2). Расчет цикла можно выполнить для любого диапазона температуры окружающей среды. Использование разработанного алгоритма термодинамического расчета не требует дополнительных итерационных расчетов, достоверные значения термодинамических параметров ГТД получаются за счет введения точных составов топливного газа. Алгоритм и программа расчета параметров ГТД разработаны с использованием самого доступного вычислительного расчетного инструмента Excel, стандартного Microsoft Office, что облегчает инженерные расчеты и позволяет корректировать исходные данные. Использование данной программы позволяет выполнять термодинамический расчет цикла для вновь создаваемых двигателей с целью получения оптимальных параметров цикла для конкретных условий эксплуатации. Оценка адекватности созданной программы термодинамического расчета ГТД проведена на основе экспериментального исследования газотурбинного двигателя MS 5002В, работающего на природном газе в составе газоперекачивающего агрегата ГТК-25И. Выполнено сравнение результатов расчета мощности с использованием разработанной модели и программы (тепловой способ) с результатами измерения мощности ГТД с помощью торсионного вала и с результатами определения мощности на муфте «ГТУ - нагнетатель» по производительности приводимого нагнетателя в соответствии с ГОСТ 20440 (эксплуатационный метод). Теплотехнические измерения проводились на трех режимах, характеризующихся различными оборотами свободной турбины. Состав природного газа, его теплофизические свойства представлены специализированной лабораторией. На основе результатов построена графическая зависимость полного давления за турбиной высокого давления (ТВД) от оборотов силового вала (рис. 3). Для трех режимов нанесены точки, соответствующие измеренному давлению и давлению, рассчитанному по вышеизложенному алгоритму и программе. Рис. 2. Алгоритм термодинамического расчета, основанный на уточненном определении теплофизических параметров рабочего тела Рис. 3. Зависимость давления за ТВД на 3 режимах Расхождение между величинами давления составляет не более 1 %. Расхождение между результатами мощности, определенной вышеуказанными способами, составляет от 0,5 до 1,8 % для всех режимов работы ГТД. С помощью разработанной программы выполнены многочисленные вариантные расчеты известных ГТД (НК-12СТ, НК-14СТ, НК-16СТ, НК-38СТ, АЛ-31СТ и др.), эксплуатируемых на компрессорных станциях магистральных газопроводов. Результаты термодинамического расчета близки к показателям их паспортных данных. При этом было установлено, что изменение удельного расхода топлива при эксплуатации ГТД в зависимости от содержания метана в топливном газе составило в среднем 12 %. Разработанная методика и программа термодинамического расчета параметров ГТД может быть применена для вновь создаваемых двигателей с возможностью определения оптимальных параметров цикла в конкретных условиях их эксплуатации. Также методика и программа применимы для мониторинга характеристик приводных ГТД, находящихся в эксплуатации на компрессорных станциях газопроводов, с определением их мощности во всем диапазоне рабочих режимов, с использованием штатных средств измерения. Это позволяет реализовать инженерные мероприятия и планировать работу этих двигателей в режимах, позволяющих повысить их эффективность. При этом получены следующие результаты: 1. Создана математическая модель учета теплофизических параметров продуктов сгорания топливного газа, которая обеспечивает получение более точных данных показателей эффективности ГТД. 2. Разработанные алгоритм и программа уточненного определения теплофизических параметров рабочего тела ГТД в широком диапазоне составов топливного газа могут быть использованы в случаях адаптации ГТД к различным видам газового топлива - природного, попутного, доменного газа, биогаза. 3. Рассчитанные значения показателей эффективности ГТД хорошо согласуются с экспериментальными данными за счет учета точного состава топливного газа, что позволяет выполнять термодинамический расчет цикла как для вновь создаваемых двигателей, так и для двигателей, находящихся в эксплуатации. Разработанная программа термодинамического расчета ГТД реализована в самом распространенном расчетном инструменте Excel и может практически применяться инженерами - разработчиками ГТД, инженерами компрессорных станций магистральных газопроводов и студента вузов.

About the authors

Svetlana A Gulina

Samara State Technical University

Email: kr_oeg@mail.ru
(Ph.D. (Techn.)), Associate Professor 244, Molodogvardeyskaya st., Samara, 443100, Russia

Anastasiya S Gulina

Samara State Technical University

Student 244, Molodogvardeyskaya st., Samara, 443100, Russia

References

Supplementary files