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

Аннотация

Разработана математическая модель топки НТКС, которая позволяет получить переходные характеристики топки при переменных параметрах топлива и режимах работы, а также в дальнейшем осуществить синтез системы автоматического управления топочными процессами.

Основной текст статьи

В настоящее время весьма актуальна проблема экономии энергетических ресурсов, и одним из вариантов ее решения является использование нетрадиционных технологий, например таких, как сжигание высокозольного твердого топлива в низкотемпературном кипящем слое (НТКС). В частности, топка НТКС может быть использована как источник теплоносителя для автономного газовоздушного воздухоподогревателя, предназначенного для обогрева ствола в зимний период. Однако одним из сдерживающих факторов внедрения подобных установок является несовершенство системы автоматизированного управления данным технологическим объектом, которая представляет собой ряд обособленных контуров управления, выполненных на основе устаревших регуляторов типа Р-25. Одним из необходимых условий для синтеза качественно новой системы автоматизации является установление математических зависимостей топки кипящего слоя

К настоящему времени имеется значительное количество работ, посвященных данной тематике. Обзор и анализ исследований, выполненных до 1989 года [1], а также более поздних работ [2, 3, 4] позволили выявить следующие недостатки существующих математических моделей:

• предполагается известной величина массового расхода топлива, тогда как в настоящее время технически реализуемо только объемное дозирование;
• параметры топлива: зольность, влажность, распределение частиц топлива по размерам предполагают заранее известными и для получения динамических характеристик системы необходимо производить каждый раз новый цикл вычислительных операций;
• не учитывается наличие термического разрушения и механического истирания частиц топлива в кипящем слое;
• сложность корректировки параметров модели при получении новых экспериментальных данных или ее использования для различных марок топлива ограничивает область применения полученных зависимостей.

Таким образом, целью работы является синтез на основе существующих аналитических зависимостей математической модели топки НТКС с решением следующих задач:

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

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

1. c_g , c_сл — удельная теплоемкость воздуха и материала слоя (Дж/(кг*К));
2. ρ_g , ρ_сл — плотность воздуха и насыпная плотность материала слоя (кг/м³);
3. S_сл — площадь зеркала горения (м²);
4. H — высота слоя (м);
5. T_о , T_сл — температура воздуха и материала слоя (К);
6. U_о — скорость газа через слой;
7. q₃ , q₄ — доля химического и механического недожога;
8. I_з , I_л — потери теплоты с отводимой золой и излучением.

Однако такой подход имеет целый ряд недостатков:

• величины с_g и с_сл при изменениях T_сл в пределах существования слоя не являются константами, а следовательно имеет смысл вместо теплоемкостей пользоваться удельными энтальпиями;
• расход теплоты при прогреве твердого топлива не учитывается;
• при определении теплоты, выносимой из слоя смесью газов, не учитывается химический состав поступающей в слой и исходящей из слоя смеси;
• массовое дозирование топлива возможно только в лабораторных условиях, а в реальных условиях применимо лишь объемное дозирование;
• при изменении расхода топлива не учитывается , что мгновенное значение теплоты, получаемое от сгорания, существенно зависит от его гранулометрического состава.

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

• масса инертного материала МСЛ, находящегося в слое, является постоянной величиной, с учетом непрерывно работающей системы золоудаления;
• топливо, поступающее в слой, имеет полидисперсный состав, который подчиняется нормальному закону распределения, переменной составляющей кторого является математическое ожидание эквивалентного диаметра частиц d_mo [5];
• исходящая из слоя смесь газов имеет температуру, равную T_сл [6];
• при работе котлоагрегата НТКС возможен переход с одной марки топлива на другую с изменением основных параметров (зольности, влажности и т.д.).

Исходя из изложенного, уравнение теплового баланса имеет вид:

1. I_сл — удельная энтельпия инертного материала слоя, ккал/кг.

В свою очередь, зависимости I_сл = f(Т_сл) и Т_сл = f(I_сл) при моделировании задаются массивом данных.

Рассмотрим последовательно все составляющие правой части выражения (2).

Количество теплоты в единицу времени — Q_gt, получаемое от сгорания жидкого топлива при розжиге:

1. V_gt — объемный расход жидкого топлива;
2. ρ_gt — плотность жидкого топлива;
3. Q_gt — низшая теплота сгорания жидкого топлива, ккал/кг.

Количество теплоты в единицу времени - Qv , вносимое в топку дутьевым воздухом:

1. T_v — температура дутьевого воздуха;
2. Q_dv — производительность дутьевого вентилятора;
3. I_v — удельная энтальпия воздуха.

Количество теплоты в единицу времени QT1 ,QT2 и QTn1 ,QTn2, вносимое в топку при сгорании сухой горючей массы топлив 1 и 2, и забираемое при нагреве этих масс до температуры слоя соответственно.

Подробнее рассмотрим определение величин QT1 и QT2. Для величин QTn1, и QTn2 алгоритм определения будет аналогичным. С целью упрощения выражений индексы топлив 1 и 2 опущены для всех переменных, кроме dt1 и dt2 – объемных долей соответствующих топлив в производительности забрасывателя. Первоначально определяется массовый расход сухой горючей составляющей топлива:

1. ρ_u — плотность сухой горючей массы;
2. ρ_p — плотность породы;
3. V_z — объемная производительность забрасывателя;
4. A_p — зольность топлива.

При моделировании предполагается, что сухая горючая масса распределяется по ряду ячеек, каждая из которых соответствует определенному диапазону диаметров частиц dmax,…dj,…dmin. Нумерация ячеек производится начиная с минимального диаметра. Текущий массовый расход для каждой ячейки определяется, как:

1. F(d,d_mo) — функция распределения частиц по диаметрам в зависимости от текущего математического ожидания.

Полагается, что все частицы, попавшие в одну ячейку, имеют одинаковый расчетный диаметр:

Во время нахождения в j-ой ячейке частицы последовательно проходят следующие стадии: прогрев до температуры слоя, выход и горение летучих, выгорание коксового остатка до , после чего происходит переход коксового остатка в ячейку j-1 [1]. Отбор теплоты на прогрев горючей массы из одной ячейки и получение теплоты от горения летучих и коксового остатка определяются из выражений:

1. t_rj,t_lj,t_gj — функция распределения частиц по диаметрам в зависимости от текущего математического ожидания;
2. Q_n — функция распределения частиц по диаметрам в зависимости от текущего математического ожидания;
3. G_trj(t) — функция распределения частиц по диаметрам в зависимости от текущего математического ожидания;
4. Q_d1, Q_d2,Q_dgt, Q_vo — оличество тепла в единицу времени, отбираемое из топки продуктами горения твердых и жидкого топлив, а также не прореагировавшим остатком дутьевого воздуха.

С учетом того, что массы выгорающих в единицу времени топлив известны, определение значений не составляет сложности. Таким образом, имеется возможность определения из выражения (2) текущей энтальпии слоя, а следовательно, и его температуры, т.е. основного технологического параметра. Предложенная математическая модель реализована в прикладном пакете MATLAB

На рис. 1а и 1в приведены полученные в разработанной математической модели реакции температуры слоя на ступенчатое изменение производительности забрасывателя для d_mo = 3 мм и d_mo = 10 мм соответственно. Моделирование производилось для следующих исходных данных: высота кипящего слоя Н = 0,4 м, площадь S = 2 м²,зольность топлива А = 50 %, влажность W = 8 %, производительность вентилятора Q_vd = 2,5 м³/с, марка угля — жирный.

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

Библиографический список

1. Махорин К. Е., Хинкис П. А. Сжигание топлива в псевдоожиженном слое. — К.: Наукова думка, 1989. — 204 с.
2. Корчевой Ю. П., Пацков В. П., Редькин В. Б., Майстренко А. Ю. Расчет выгорания частиц твердого топлива в кипящем слое с учетом внутрипористого реагирования// Тепломассообмен ММФ-92. Тепломассообмен в дисперсных системах: Т.5.-Минск. АНК"ИТМО им. А. В. Пылова", АНБ. — 1992.— С. 168–170.
3. Бубенчиков A. M., Старченко А. В., Стропус В. В. Математическое моделирование аэродинамики и тепломассопереноса в устройствах с циркулирующим кипящим слоем // Теплоэнергетика. — 1995. — № 9. — с. 37–41.
4. Рохман Б. Б., Шрайбер А. А., Чернявский Н. В. Инженерная методика расчета cгорания твердых топлив в реакторе с циркулирующим кипящим слоем применительно к пилотной установке по технологии фирмы "Лурги" // Пром. теплотехника. — 2004. — т. 26. — № 4. — с. 40–47.
5. Бородуля В. А., Гупало Ю. П. Математические модели химических реаторов с кипящим слоем. — Мн.: Наука и техника, 1976, 208 с.
6. Забродский С. С. Гидродинамика и теплообмен в псевдоожиженном (кипящем) слое. М.–Л., Госэнергоиздат, 1963. 488 с.
7. Бородуля В. А., Пальченок Г. И., Васильев Г. Г., Дрябин В. А., Галерштейн Д. М. Тепломассообмен и кинетика горения твердого топлива в кипящем слое // Проблемы тепло- и массообмена в современной технологии сжигания и газификации твердого топлива. Материалы международной школы-семинара. Минск, 1988, ч. 2. — с. 3–23