Способ управления нагревом заготовок в многозонной пламенной печи Советский патент 1993 года по МПК F27B9/40 

венное расположение расчетных участков, б) схема массопереноса (стрелки) и положение факела; на фиг. 2 температурные поля вдоль свода печи: кривые 1, 2, 3 - для факелов 5, 7 и 9 м соответственно при тепловой нагрузке 800 м /ч и производительности печи 120 т/ч; 1, 2, 3 - то же, при тепловой нагрузке 1050 м3/ч; 1, 2, 3 - то же, при тепловой нагрузке 550 м/ч; на фиг. 3 - профили удельных поглощенных поверхностью нагреваемых заготовок потоков излучения от рабочего пространства по длине зоны печи: 5, 4, б - для факелов длиной 5, 7, 9, соответственно при тепловой нагрузке 800 м3/ч и производительности печи 120т/ч; на фиг. 4 - функции изменения параметра настройки модели (К/L, ф/ для длины факелов 5 м (кривые 7, 9) и 9 м (кривые 8, 10) и различного места установки термопар: кривые 7, 8 для термопары ТПП. кривые 9, 10 для термопары ТПП.

Согласно фиг. 1 математическая модель теплотехнической зоны печи содержит 65 расчетных участков, а именно: 17 объемных участков газа, 30 поверхностных участков кладки и 9 поверхностных и 9 объемных участков металла.

Модель позволяет учитывать следующие особенности теплообмена: неоднородность радиационных характеристик среды по объему, развитию и положение факела относительно тепловоспринимающих поверхностей, специфику аэродинамики газов и выгорания топлива по длине факела и объему рабочего пространства, распределение потерь тепла между отдельными элементами печи, конвективную радиационные составляющие теплообмена. Теплофизические свойства кладки и металла соответствуют материалам,применяемым в кладки печи и нагреваемому металлу. Степень черноты металла и кладки принималась 0,8. Коэффициенты конвективного обмена на поверхности принимались 25-35 вт/м . Теплота сгорания топлива 38,4 М Дж/м .

Длина факела определялась по степени топлива. Полагалось, что полная длина факела равна участку, на котором выгорает 98% поданного в факел топлива. Для описания процесса выгорания топлива и подсоса в факел воздуха использованы известные соотношения (см. стр. 122 ск. Лисиенко В. Г., Волков В. В., Гончаров А. Л., Математическое моделирование теплообмена в печах и агрегатах, Киев, Наукова думка, 1984. с. 224),

k 1 -e

(x In),

0)

к 0,,85, (), (2) 1ф In

(1 -e-mnx),(3) где к- степень выгорания топлива, х - расстояние от выходного сечения горелки до рассматриваемого сечения, м;

От - безразмерный подсос воздуха в зону горения, представляющий собой соотношение количества подсосного воздуха на пути до рассматриваемого сечения факела к теоретически необходимому при а 1;

Ого - коэффициент расхода воздуха в продуктах горения на выходе из рабочего пространства зоны. Коэффициент m определяется из условия, что в пределах длины пути подсоса выгорает 85% топлива, а коэффициент получается из условия, что топливо в факеле выгорает на 98% (ic 0,98), т.е. в конце полной длины факела недожог топлива равен 2%.

Коэффициенты поглощения объемных участков модели определяли из концентрации излучающих компонентов приходящихся на соответствующие участки. Парциальные давления углекислого газа и водяного пара( и ) по длине факела определяются по соотношению:

РХ Рп.г.

/cV0

FVo + (1 - А:) + (От - к) и,

(4)

где Р п.г. - парциальное давление излучающего компонента в продуктах полного сгорания при а 1 атм;

VQ - теоретический выход продуктов сгорания при при нормальных условиях, м /м ; - теоретически необходимое количество воздуха при нормальных условиях сжигания 1 м топлива.

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

п + m|

Б AijTi4+ 2 gijTi+ , i 1

где Ti - температуры расчетных зон, К;

Aij - коэффициенты радиационного теплообмена между зонами i и j, Вт/К ;

gij - коэффициент конвективного теплообмена или переноса тепла газами между зонами i и j, Вт/К :

QJ- внутреннее тепловыделение в зоне j или тепло,проводимое из другой зоны, Вт;

n,m - количество объемных и поверхностных зон модели.

Количество тепла, поглощенное каждой j-ой из К объемных зон металла, равно:

WMI TMI-WMJ TMJ

+ (Jnj TMj) FnJ

/kmЛм

где: WMI; WMJ - теплоемкость потока,входящего в зону из зоны i и выходящего из зоны

JJ

TMI, TMJ - среднемассовые температуры зон i и j металла, К;

Tnj, Tnj -температура и площадь поверхности зоны металла ограничивающей объемную зону металла j, К, м ;

50kj, Aokj - толщины и коэффициент теплопроводности слоя окалины, м, Вт/(м К);

дм. AMJ - то же для металла;

i/j- коэффициент формы массивного тела для слябов равен 1/3.

Для того, чтобы найти температурные поля в рабочем пространстве печи, необходимо решить систему (лп+п) уравнений вида (5) и К уравнений вида (6).

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

Анализ результатов вычислительных экспериментов позволит сделать следующие выводы:

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

2.4, На своде печи имеется область с центром в точке А (см, фиг. 2), температура в которой практически не зависит от длины факела, но хорошо отражает тепловую нагрузку зоны (см. серии кривых на фиг. 21,2, 3; 1, 2, 3 ; 1.2.3).

3. Температурное поле поверхности свода по длине зоны печи существенно неравномерно. Область наибольших температур находится над ядром факела. И именно эта неравномерность температурного поля кладки по длине зоны печи определяет существенную неравномерность распределения падающих тепловых потоков на поверхность нагреваемого металла - (см. фиг. 3).

54. Рассмотрение распределения температуры греющей среды по длине печи показывает, что для различных длин факелов это распределение также различно (см. фиг. 4). Поэтому однозначно определить распреде0 ление температуры греющей среды по длине зоны одной термопарой невозможно.

5. Кроме того, расчеты подтвердили сделанный ранее вывод, что при постоянной тепловой нагрузке эффективней работа на

5 коротком факеле,

Полученные результаты расчетов и сделанные выводы позволяют принять следующие технические решения; а) температуру теплотехнической воды печи следует конт0 ролировать в области,не зависящей от длины факела. Поскольку регулятор топлива поддерживает необходимую температуру в теплотехнической зоне, и, тем самым, обеспечивает конечную температуру нагрева за5 готовок, до датчик температуры зоны, обеспечивающий обратную связь регулятора топлива, в заявленном техническом решении мы называем основным;

б) высокотемпературная область кладки 0 над ядром факела, согласно заявленному техническому решению, контролируется вторым датчиком температуры - дополнительным. Этот датчик используется в обратной связи регулятора длины факела, 5 который подбирает оптимальную с точки зрения теплоусвоения длину факела;

в) сравнивая температуры,измеренные основной и дополнительной термопары, возможно однозначно определить длину.

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

г) поскольку основная термопара измеряет температурный урове нь зоны (см. п.а.) и определен профиль распределения грею0 щей среды по длине зоны (см. п. в), то несложно воспроизвести распределение температуры греющей среды по длине зоны.

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

горелки не обладают стабильностью длины факела во всем диапазоне работы.

Эффективность использования топлива достигается выбором минимальной длины факела. Однако, уменьшение длины факела ведет к повышению температуры кладки печи (см. фиг. 2). Это увеличение температуры может достигнуть и предельнодопустимых значений из условий стойкости кладки. В заявляемом способе предлагается регулировать длину факела, путем ее сокращения потока температура кладки не достигнет предельного значения. При достижении предельной температуры кладки длину факела увеличивают, тем самым растягивается процесс выгорания топлива и тепловыделение в зоне печи распределяется более равномерно. Кроме того, эффективность возрастает вследствие исключения ложных срабатываний регулятора расхода топлива при изменении длины факела, что обеспечивается специальным выбором места установки основного датчика температуры. Исключение ложных срабатываний регулятора продлевает срок службы исполнительных механизмов.

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

qpy(L) Ob(L)K(L, 1ф)Ттп4-Тм4/д, (7) где L- координата положения сляба в печи, м;

JB(L) - приведенный коэффициент внешнего теплообмена по длине печи, Вт/(м2 К4);

Ттп-температура в зоне печи, измеренная основной термопарой, К;

TM(L) - температура поверхности заготовки, находящейся в координате L, К;

ф - текущая длина факела в зоне печи, м;

K(L, ф) - функция,определяющая связь профиля распределения температуры греющей среды по координате с температурой измеренной термопарой в локальной области теплотехнической зоны.

Как видно из уравнения (7) настроечными параметрами являются овО-) и К/L, ф).

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

Для этого используя основные выражения для расчета коэффициентов радиационного теплообмена(Лисиенко В. Г., Волков В.

В., Гончаров А. Математическое моделирование теплообмена в печах и агрегатах - Киев: Наукова думка, 1984, 232 с) согласно зональному методу уравнение (5) возможности свести к виду (7) для каждой поверхностной зоны металла при условии:

оВ(1ф) Ан(1ф)(8)

25

К;(1ф) V ТАи(1ф)Т 4/А Юф)ТтгЛ (9)

где ой;(ф), К(ф), Ау(ф) - коэффициенты уравнения (5), (7) рассчитанные для фиксированной длины факела ф,

Ттп - температура расчетной зоны модели, которая совпадает с местом установки термопары контролирующей температурный уровень зоны нагрева.

Изменение параметра Kj(L, ф) подлине

печи зависит от выбора места установки датчика температурного уровня зоны, т.к. формула (9) определяет связь температуры измеряемой этим датчиком и падающим потоком излучения на поверхность металла в

конкретной координате. На фиг. 4 приведено изменение параметра K(L, 1ф) для различных мест установки термопар.

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

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

Блок-схема устройства для реализации заявленного способа представлена на фиг. 5. Система содержит зону 11 многозонной печи с торцовым отоплением, горелку 12 с регулируемой длиной факела, регулирующий орган 13 расхода топлива, регулятор 14 расхода топлива, основную термопару 15, математическую модель 16, блок 17 слежения за заготовками, регулирующий орган 18 длины факела, регулятор 19 длины факела, дополнительную термопару 20, задатчик 21 максимальной температуры кладки, блок 22 идентификации распределения температуры греющей среды.

Согласно заявляемому способу управ- ление нагревом металла осуществляется следующим образом.

В зоне 11 печи с торцовым отоплением устанавливают горелки 12 с регулируемой длиной факела и регулирующим органом 18. Предварительно определяют места установки термопар основной 12 и дополнительной 20. При этом основную термопару 15 устанавливают в область,не зависящую от длины факела, а дополнительную термопару 20 над ядром факела, т.е. в область мак- симальных температур кладки. Регулирование длины факела осуществляется регулятором 19 независимо от расхода топлива в зону 11 печи по рассогласованию между заданной максимально допустимой температурой кладки, заданной задатчиком 20 и текущей максимальной температурой кладки,измеренной дополнительной термопарой 20.

.В качестве горелки 12 с изменяемой длиной факела возможно использование известной горелки типа ФСГ-Р, в которой изменение длины факела достигается за счет изменения условий перемешивания топлива и воздуха путем перераспределения топлива между центральными и периферийными подводами.

Температурный уровень в зоне 11 печи поддерживается регулятором 12 расхода топлива, через исполнительный механизм 13, в зависимости от разности температуры измеренной основной термопарой 15 установленной в области А кладки зоны 11 печи и заданного значения температурного уровня зоны рассчитанного по математической модели 16. Математическая модель 16 рассчитывает температуры заготовок по их координатам в печи и определяет требуемый температурный уровень зоны согласно критерию управления. Информация о положении заготовок в зоне и их параметрах поступает в математическую модель из блока 17 слежения за заготовками.

Математическая модель 16 рассчитывает температуру заготовок в соответствии с распределением температуры греющей среды, профиль которой определяется блоком 12 идентификации распределения температуры греющей среды в зависимости от рассчитанной по разности температур основной 15 и дополнительной 20 термопар длины факела зоны (фиг. 4). Практическая

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

Экономический эффект при использовании заявляемого изобретения в составе АСУТП нагрева металла составит по оценкам сделанным с помощью автоматической модели 3% экономии топлива за счет оптимизации длины факела.

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

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

поддерживая их минимально возможными.