В свою очередь, окисленность шлака (рис. 3 б) в конечный период плавки растет быстрее, так как при снижении содержания углерода ниже [С] = 0,2 ± 0,05 % доля кислорода, идущего на реакцию с углеродом, уменьшается, а на окисление железа, соответственно, увеличивается, что может увеличить потери железа со шлаком.
Кривые скорости нагрева шлако-металлического расплава от обезуглероживания металла (рис. 3 в) имеют вид экспоненциальной зависимости, максимальны в начальный период и непрерывно снижаются в связи с постоянным увеличением массы шлако-металлического расплава.
Таким образом, применение математической модели обезуглероживания металла в ДСП и программы расчета составляющих уС при непрерывной подаче ЖМО по ходу плавки с интенсификацией процесса шлакообразования кислородом от ТКГ позволяет рассчитывать кривые содержания углерода в металле, окисленности шлака и скорости нагрева металла в условиях обезуглероживания в течение расчетного периода электроплавки стали с учетом непрерывной загрузки ЖМО при подаче кислорода через ТКГ, что позволяет использовать модель и ее программу для расчета управляемого окислительного рафинирования при непрерывном плавлении окатышей в ванне дуговой печи.
Проверку модели на адекватность осуществляли путем сравнения расчетных данных уС и (БеО) с данными опытных плавок из паспортов [2], [5] и результаты сравнения (разброс данных не более 10 %) свидетельствует о возможности практического использования модели (рис. 1) для управления режимом плавления ЖМО в ДСП.
Литература
1. Бакакин, В.А. Исследование кинетических закономерностей обезуглероживания металлического расплава /
B.А. Бакакин, С.И. Филиппов, С.Н. Падерин, М.М. Волков // Изв. ВУЗов «Черная металлургия». - 1972. - № 5. - С. 5 - 7.
2. Бартенева, О.И. Исследование процессов нагрева и обезуглероживания металла в 150-т дуговой печи с переменной массой ванны / О.И. Бартенева, Э.Э. Меркер // Изв. ВУЗов «Черная металлургия». - 2001. - № 9. - С. 65 - 66.
3. Григорян, В.А. Теоретические основы электросталеплавильных процессов / В.А. Григорян, Л.Н. Белянчи-ков, А.Я. Стомахин. - М., 1987.
4. Меркер, Э.Э. Анализ процесса самораскипания сталеплавильной ванны при продувке металла кислородом / Э.Э. Меркер // Изв. ВУЗов «Черная металлургия». - 2001. -№ 6. - С. 62 - 64.
5. Меркер, Э.Э. Энергосберегающие условия плавления металлизованных окатышей в ванне дуговой печи / Э.Э. Меркер // Бюллетень «Черная металлургия». - 2008. Вып. 1 (1297). - С. 35 - 39.
6. Падерин, С.Н. Термодинамика и расчеты процесса глубокого обезуглероживания стали / С.Н. Падерин, Е.П. Падерина // Изв. ВУЗов «Черная металлургия». - 2005. -№ 10. - С. 19 - 24.
7. Падерин, С.Н. Термодинамическое моделирование окислительных процессов при обезуглероживании стали /
C.Н. Падерин, П.С. Падерин, И.В. Кузьмин // Изв. ВУЗов «Черная металлургия». - 2003. - № 5. - С. 6 - 11.
8. Трахимович, В.И. Использование железа прямого восстановления при выплавке стали / В. И. Трахимович, А.Г. Шалимов. - М., 1982.
9. Филиппов, С.И. Теория металлургических процессов / С.И. Филиппов. - М., 1967.
10. Филиппов, С.И. Теория процесса обезуглероживания стали / С.И. Филиппов. - М., 1956.
11. Хмелева, С.Л. Математическое моделирование процесса обезуглероживания в электродуговой печи / С. Л. Хмелева, А.Г. Падалко // Сборник трудов Всероссийской научно-практической конференции «Моделирование, программное обеспечение и наукоемкие технологии в металлургии». - Новокузнецк, - 2011. - С. 347 - 355.
УДК 621.778.04.
Р.А. Юдин, Н.И. Шестаков, В.Р. Аншелес, И.Р. Юдин, Н.А. Тувалин ОСНОВНЫЕ ЗАКОНОМЕРНОСТИ НЕПОЛНОГО СГОРАНИЯ ПРИРОДНОГО ГАЗА
В статье приведены аналитические исследования, позволившие выявить основные закономерности неполного сжигания природного газа. Показано, что сумма объемов оксида и диоксида углерода в продуктах неполного сгорания равна стехио-метрическому объему диоксида углерода: СО2тах = CO + CO2, а сумма объемов водорода и водяных паров - стехиометриче-скому объему водяных паров: H2Omax = H2 + H2O. Получены аналитические формулы, позволяющие рассчитывать химический состав продуктов неполного сгорания при различных режимных параметрах.
Оксид и диоксид углерода, водород и водяные пары, коэффициент расхода первичного воздуха.
The paper presents the analytical studies of the basic laws of incomplete combustion of the natural gas. It is shown that the sum of the amounts of oxide and carbon dioxide in the products of incomplete combustion is the stoichiometric volume of carbon dioxide is СО^^ = CO + CO2, and the sum of the amounts of hydrogen and water vapor is equal to the stoichiometric amount of water vapor: H2Omax = H2 + H2O. The analytical formulas of calculating the chemical composition of the products of incomplete combustion at various regime parameters are received.
Oxide and carbon dioxide, hydrogen and water vapor, flow rate of primary air.
В металлургии и машиностроении для обеспечения мало- и безокислительной термохимической обработки стали используют печи открытого нагрева с двухстадийным сжиганием природного газа. На первой стадии его сжигают при существенных недостатках воздуха, отвечающим коэффициентам расхода первичного воздуха а! = 0,45 - 0,60. На второй стадии продукты неполного сгорания дожигают соответственно с коэффициентами расхода вторичного воздуха а2 = 0,60 - 0,45. Указанным значениям отвечает суммарный коэффициент расхода воздуха аЕ = 1,05, близкий к стехиометрическому значению.
Следует отметить, что в машиностроении эти печи именуют безокислительными, так как окисление поверхности металлопродукции удовлетворяет требованиям точной штамповки. Однако в металлургии такой подход неприемлем. Это связано с тем, что, например, при нагреве проволоки, несмотря на наличие на поверхности только цветов побежалости, измеряемых в ангстремах, в непрерывном технологическом процессе после термохимической обработки необходимо дополнительно ее подтравливать, поэтому в металлургии такой нагрев называют малоокислительным.
Преимуществом открытого малоокислительного нагрева является то, что при сжигании природного газа одновременно с реализацией заданных температурных технологических режимов достигается эффективная защита поверхности металлопродукции от окисления и обезуглероживания. Природный газ является наиболее приемлемым естественным топливом, позволяющим решить указанную задачу.
Существуют проходные и протяжные печи непрерывного действия с последовательными и параллельными зонами неполного сжигания и дожигания [6] - [7]. В печах с последовательными зонами продукты сгорания движутся противотоком к нагреваемой металлопродукции. Такая организация ее термохимической обработки позволяет обеспечить окончательный нагрев до заданных технологических температур в малоокислительной атмосфере продуктов неполного сгорания, содержащей помимо окислительных и нейтральных восстановительные компоненты: оксид углерода (СО) и водород (Н2). Из зоны окончательного нагрева продукты неполного сгорания противотоком к металлу поступают в зону предварительного нагрева, где их дожигают вторичным воздухом. При этом предварительный нагрев металла осуществляют до температуры не более 600°С, что позволяет минимизировать его окисление.
Печи с параллельными зонами неполного сжигания и дожигания разделены по высоте теплопроводным промежуточным подом, выполненным, например, из карбида кремния, что практически позволяет полностью исключить контакт металлопродукции с окислительными продуктами полного сгорания.
Однако для повышения эффективности работы печей открытого малоокислительного металлопродукции с двухстадийным сжиганием природного газа весьма важно не только грамотно решать задачи конструирования и автоматизации, но эффективно рас-
считывать химический состав продуктов неполного сгорания при различных коэффициентах а1 и других теплотехнических параметрах, например, при различных температурах нагрева первичного воздуха. В общепринятой инженерной практике расчеты ведут, допуская, что в продуктах неполного сгорания полностью отсутствуют кислород, несгоревшие углеводородные газы и сажистый углерод.
Общепринятая методика расчета химического состава продуктов неполного сгорания газообразного углеводородного топлива разработана В.Ф. Копыто-вым [2]. Наиболее детально эта методика и пример ее использования представлены в монографии Б.С. Мастрюкова [3, с. 15 - 16, с. 27 - 28]. Методика В.Ф. Копытова основана на том, что согласно термодинамическому принципу компоненты продуктов сгорания находятся в термодинамическом равновесии по следующей обратимой физико-химической реакции:
СО2 + Н2 « СО + Н2О. (1)
Константа термодинамического равновесия этой реакции является табулированной величиной, значения которой в зависимости от температуры определяют по формуле:
Хро * и н
к = —^ = у (/), (2)
усо * 2н2о
где х, у, и, 2 - объемы означенных в формуле компонентов в продуктах неполного сгорания в м3/м3.
Для определения химического состава продуктов неполного сгорания составляют материальный баланс химических элементов, участвующих в реакции
(1). Для этого первоначально составляют систему, состоящую из трех линейных алгебраических уравнений, соответственно для углерода, водорода и кислорода. В левой части каждого из этих уравнений в принятых обозначениях неизвестных переменных записывают количество данного элемента в продуктах неполного сгорания, а в правой - в исходном топливе и воздухе, подаваемом на горение при различных коэффициентах расхода первичного воздуха а1. Четвертым уравнением этой системы является формула константы термодинамического равновесия
(2), при подстановке в которую трех неизвестных образуется квадратное уравнение. Данное уравнение следует решать относительно одного из компонентов реакции (1).
Решение системы уравнений по методике В.Ф. Копытова трудоемко тогда, как в монографии [1, с. 101 - 106] А.Е. Лившицем показано, что искомые неизвестные являются функцией трех стехиометри-ческих параметров процессов горения: стехиометри-ческих объемов диоксида углерода, водяных паров и кислорода, необходимого для полного сгорания топлива, что позволяет существенно уменьшить объем расчетов.
Однако в монографии [4, с. 129 - 130] показано, что для расчета химического состава продуктов не-
полного сгорания необходимо использовать только два из указанных параметров, так как значение любого из трех указанных параметров несложно выразить по двум остальным. Предпочтительнее использовать стехиометрические объемы диоксида углерода и водяных паров, так как именно эти параметры являются суммой соответствующих продуктов неполного сгорания. Эта закономерность позволяет записать рассматриваемую систему линейных алгебраических уравнений в следующем виде:
хСО2 = СО2 + СО2шах " УсО
МИ2 + 2НгО = Н2Отах (3)
2ХС02 + 2НгО = 2СО2 " 2О2 +
+ а, (2СО2шах + Н2Ошах + 28О2 - СО) - уго.
Здесь СО2шах, Н2Ошах - максимальные содержания диоксида углерода и водяных паров в продуктах сгорания, соответствующие стехиометрическим условиям; СО2, О2 - содержание диоксида углерода и кислорода в исходном топливе; 8О2 - содержание сернистого газа в продуктах сгорания, соответствующее содержанию сероводорода в исходном топливе в м3/м3.
Эту систему несложно решить относительно одного из искомых параметров, например, уСО, используя или формулы Крамера, или метод Лагранжа. Здесь предпочтительнее использовать формулы Крамера, так как содержание хСОг определяется
только объемом уСО, а определитель системы (3) равен единице. Действительно:
Д =
1 0 0 0 1 1 2 0 1
= 1.
Поэтому хСО = Дх, иН = Ди , а гН О = Дг . Тогда:
1 0 СО2 + СО21пах - УСО
0 1 Н.Ошах
2 0
2СО2 -2О2 +а1(2СО2шах + ^О^) +
+а1(2Н28 - СО) - Усо
= а (2СО2шах +Н2Ошах )-2СО2шах + + а1 (2£Ог - СО) - 2О2 + уСО.
(4)
к--
Подставив хСО , иН , гН О в формулу константы равновесия (2) получим:
._СОгш„[ + СОг - уСО__
уСО [аДгСО^ + НгОшах) - 2СОгш!|1[ + а1(28Ог - СО) - 2Ог + уСО ]'
х [(1 - а )(2СО2шах + НгО„„) - а. (28Ог - СО) + 2Ог - уш.
(6)
Для топлива, содержащего только углеводородные горючие газы и негорючие компоненты, эта формула имеет следующий вид:
к =
(СОгшах + СОг - Усо) • [(1 -а1) • (гСОгшах + НгОшх) - Усо]
Усо[а1(2СОгтш( + НгОшх) - гСОгшах + Усо]
(6.1)
Расчеты неполного сгорания природного газа следует проводить, используя метод итераций. Это обусловлено тем, что константа термодинамического равновесия, как и значения теплоемкостей компонентов продуктов сгорания, является функцией температуры. Для уменьшения числа итераций еще в большей степени, чем при полном сгорании топлива, целесообразно использовать хотя бы приблизительные зависимости искомых температур от коэффициента расхода первичного воздуха а!. Получение таких зависимостей является весьма актуальной задачей.
Теплотехнические параметры метана весьма близки с характеристиками природного газа различ-
ных месторождений. Для метана СО2
= 1 м3/м3
Н2Ошах = 2 м /м, а низшая теплота сгорания Qн= 35820 кДж/м3. Для природного газа, поступавшего по магистральному газопроводу, например, 12.04.12 г. в г. Череповец 0НР= 35997 кДж/м3, СО2шах = 1,0074 м3/м3, Н2Ошах = 1,9988 м3/м3. В связи с этим получение и использование табулированных и графических зависимостей теплотехнических параметров процессов неполного сгорания метана от коэффициента расхода первичного воздуха а! при расчетах неполного сжигания природного газа вполне приемлемо.
При неполном сжигания метана (0,50 < а! < 1,0) формула (6) принимает следующий вид:
к =
(1 - Усо) • [4(1 -а1)- Усо] Усо(4а1 - 2 + Усо)
(6.2)
Согласно второму уравнению системы получим:
При неполном сжигании природного газа (0,50 < а1 < 1,0), поступающего в г. Череповец, формула (6) принимает следующий вид:
иНг = НгОшах - ^ = ^О^ - а1 (СОг^ + ЩО^ )--а1(2Нг8-СО) + 2Ог +2СОгшах - Усо = = (1 - а1) (гСОгшах + НгОшах ) + 2Ог - а1 (28Ог - СО) - уго.
(5)
к = (1,0078 - Усо) • [4,0136 • (1 -а1)) - Усо] (6 3) Усо [4,0136а1 - 2,0148 + Уго] ' '
Формулу (6.2) несложно привести к следующим
квадратным уравнениям, корнем которых является содержание оксида углерода в продуктах неполного сгорания - >>сО:
(1 - к) уСо - [3 - 2(2а1 -1)(1 - к)] усо + 4(1 - а,),
к,уСо - [3 - 2к1 (2а1 -1)]усо + 4(1 - а1). (7)
Очевидно, что константа термодинамического равновесия к1 = 1 - к, однако квадратное уравнение (7) имеет более компактный вид. Зависимость к от температуры представлена в таблице.
Таблица
Зависимость от температуры к1 = (1 - к) = Г
Получить компактные квадратные уравнения из формулы (6.3) для Череповецкого природного газа и тем более из формулы (6) для топлива произвольного химического состава без введения дополнительных алгебраических комплексов как, например, в монографии [1, с. 101 - 106] весьма затруднительно. В связи с этим необходимые уравнения следует получать, предварительно рассчитав СО2тах, Н2Отах, например, как в (6.3), и, задавая коэффициент расхода первичного воздуха а1. Квадратное уравнение (7) имеет следующее решение:
_ 3 - 2к1(2а1 -1) - ^[3 - 2к1 (2а1 -1)]2 - 16к1 (1 - а1) усо _ 2к
(8)
Объемы других компонентов продуктов неполного сгорания метана следует рассчитывать по сомножителям числителя и знаменателя формулы (6.3). Следует отметить, что компактность (7 - 8) во многом обеспечивается тем, что для метана, как показано выше, используемые стехиометрические параметры являются целыми числами, определение которых не представляет трудностей. Однако основная сложность при расчетах химического состава продуктов неполного сгорания заключается в том, что, как следует из таблицы, он зависит не только от коэффициента расхода первичного воздуха аь но и от температур горения топлива.
При полном сгорании, как показано в работе [5], химический состав продуктов сгорания зависит исключительно от коэффициента избытка воздуха as и не зависит от указанных температур. В связи с этим весьма актуальной является разработка современного алгоритма расчета химического состава продуктов неполного сгорания и соответствующих ему температур горения топлива при различных режимных параметрах, в котором используются полученные закономерности.
Литература
1. Ващенко, А.И. Окисление и обезуглероживание стали / А.И. Ващенко, А.Г. Зеньковский, А.Е. Лившиц, Л.А. Шульц. - М., 1972.
2. Копытов, В. Ф. Сжигание газа с недостатком воздуха / В.Ф. Копытов // Сталь. - 1954. - № 3. - С. 57 - 62.
3. Мастрюков, Б.С. Теплотехнические расчеты промышленных печей / Б.С. Мастрюков. - М., 1972.
4. Юдин, Р.А. Научно-технические основы процессов горения / Р.А. Юдин. - Череповец, 2013.
5. Юдин, Р.А. Особенности горения топлива произвольного химического состава / Р.А. Юдин, И.Р. Юдин // Вестник ЧГУ. - 2013. - № 1. - С. 15 - 18.
6. Юдин, Р.А. Повышение эффективности работы печей метизного производства / Р.А. Юдин // Производство проката. - 1998. - № 5. - С. 36 - 46.
7. Юдин, Р.А. Совершенствование конструкций и способов отопления печей-ванн горячего оцинкования / Р.А. Юдин, Н.И. Шестаков, И.Р. Юдин // Сталь. - 2004. - № 9. -С. 45 - 46.
Температура, °С 1200 1300 1400 1500 1600 1700 1800
Константа равновесия, к 0,607 0,655 0,703 0,74 0,758 0,776 0,789
Температура, °С 1900 2000 2100 2200 2300 2400
Константа равновесия, к 0,801 0,813 0,825 0,837 0,841 0,845