CC BY
62
УДК 621.311.22:662.756
И.А. Ростунцева
МОДЕЛИРОВАНИЕ ПРИРОДООХРАННЫХ ТЕХНОЛОГИЙ В ТЕПЛОЭНЕРГЕТИКЕ С УЧЕТОМ МАКРОКИНЕТИЧЕСКИХ ПРОЦЕССОВ
Представлена аналитическая модель оценки экологического воздействия теплоэнергетического оборудования на окружающую среду. В основу разрабатываемой модели положен макрокинетический закон протекания брутто-реакции образования оксидов азота в зоне высоких температур, отражающий закономерности образования оксидов азота по термическому механизму. Получены концентрации кислорода и азота в зоне реакции в зависимости от режимов работы горелок и топки котла, а также от возможного их изменения за счет ввода добавочной влаги. Определены аналитические зависимости скорости реакции, концентрации оксидов азота и коэффициента выхода оксидов азота при сжигании топлива с учетом ввода влаги.
Теплоэнергетическое оборудование, концентрация оксидов азота, аналитическая модель, скорость реакции образования оксидов азота, константы реакции
I.A. Rostuntseva
MODELING OF ENVIRONMENTAL TECHNOLOGIES THERMAL ENGINEERING BASED MACROKINETIC PROCESSES
Presented an analytical model for evaluating the environmental impact of power equipment on the environment. The basis of the model of the law laid macrokinetic flow gross reaction of nitrogen oxides in the zone of high temperatures, reflecting the regularities of formation of nitrogen oxides on the thermal mechanism. Received oxygen and nitrogen in the reaction zone, depending on the modes of operation of burners and furnaces of the boiler, as well as their possible changes due
to entry of additional moisture. Defined analytical reaction rate, concentration of nitrogen oxides and exit rates of nitrogen oxides from fuel combustion, taking into account the input of moisture.
Power equipment, the concentration of nitrogen oxides, an analytical model, the rate of reaction of nitrogen oxides, the constants of the reaction
При проектировании нового и модернизации действующего теплоэнергетического оборудования ТЭС необходима теоретическая оценка его экологического воздействия на окружающую среду, позволяющая быстро и с достаточной точностью определять допустимое загрязнение атмосферного воздуха и степень снижения вредных выбросов при внедрении природоохранных технологий.
Первым этапом в разработке аналитической модели образования вредных выбросов в факеле энергетической установки является исследование физико-химических механизмов образования вредных выбросов, в частности одного из наиболее токсичных компонентов в составе дымовых газов - оксидов азота (NOx) при сжигании различных видов топлива.
В рамках данной разработки проведено теоретическое исследование возможности образования термических, «быстрых» и топливных оксидов азота на различных уровнях прохождения дымовых газов по тракту энергетической установки, а также факторов, влияющих на количественную характеристику выхода оксид азота при сжигании топлива.
Из проведенного анализа установлено, что преобладающее влияние на формирование концентраций NOx в зоне активного горения оказывают термические оксиды азота, для снижения, образования которых целесообразно регулировать температурный уровень и концентрацию свободного кислорода, что, возможно, осуществлять путем ввода влаги в топочные процессы.
Исходя, из последнего разработана аналитическая модель исследования образования вредных выбросов с учетом ввода влаги в зону горения.
В основу разрабатываемой модели положен макрокинетический закон протекания брутто-реакции образования оксидов азота в зоне высоких температур. Методика отражает закономерности образования NOx в основном по термическому механизму Я.Б.Зельдовича, исходя из избирательного воздействия дополнительно вводимой влаги преимущественно на термические оксиды азота. В методике определено влияние, оказываемое на образование оксидов азота температурного уровня, концентрации реагирующих веществ и время нахождения их в зоне высоких температур.
При разработке модели образования оксидов азота с учетом ввода влаги использовались известные положения, изложенные в работах Я.Б.Зельдовича, И.Я.Сигала, Т.Б. Эфенди-ева и ряда других авторов [1-5].
Для определения концентраций оксидов азота, образующихся в процессе горения, использовался метод разложения экспоненты, согласно которому концентрация NO находится из выражения
dC no
Cno= (—N°)-T NO, (1)
dx
dCNO
где (---NO) - скорость реакции образования оксидов азота при максимальной температуре
dx
факела Тт , % об./с; Т NO - условное время реакции, с.
С учетом макрокинетического закона выход оксидов азота будет описываться уравнением
dCNO
(_N°) = 2 . КП ■ Со2 ■ CN2 , (2)
dx
где Кц - константа скорости реакции, (% об./с)-1; С02 и Ск2 - объемные концентрации соответственно атомарного кислорода и молекулярного азота в зоне максимальной температуры Тт , % об.
Концентрация атомарного кислорода в зависит от максимальной температуры в зоне горения. При температурах выше (1800-2000) 0К ее можно найти из условия равновесия при горении с избытком окислителя по формуле
Со = [К] ■ д/СО2 , (3)
где [К] - константа равновесия , (% об./с)-0,5 , С02- концентрация кислорода в продуктах сгорания, % об.
Константы Кц и [К] могут быть определены по выражениям:
КП = 6,1-107 ■ exp ■ [-R5400 ], (4)
R ■ Tm
[К]= 0,75 ■ 104 ■ exp ■ [ 59300 ], (5)
K ■ Tm
где R - универсальная газовая постоянная (R=8,326 Дж/моль ■ К ), Tm - максимальная температура факела, 0К.
Подставив (3)-(5) в (2), получим следующее соотношение для определения скорости образования оксидов азота:
;d^NO)= 2 ■ Kii ■ VC07-Cn2= 9,15 ■ 1011 ■ 7CO2-Cn2 ■ exp ■ (-^^; dT v 02 v O2 r . Tm
Концентрации кислорода и азота в зоне реакции определяются согласно принятому в методике допущению при окислении азота воздуха. При этом указанные концентрации зависят от параметров, характеризующих режимы работы горелок и топки котла, а также от возможного их изменения за счет ввода добавочной влаги.
21 ■ у0- (а + 0,5 ■ Да -1)
Со2=---------------------------------------------------------------гор-т- , (7)
[У0 + (а + 0,5 ■ Да -1) ■ у0+Ф ■ g ]
г гор т п впр
79 ■ у0 ■ (а + 0,5 ■ Да -1)
С _ У гор т (8)
Ст_~0-----------------------------------------------------------0- , (8)
[У0 + (а + 0,5 ■ Да -1) ■ у0 + £ ■ g ] г гор т у п впр
где агор - коэффициент избытка воздуха в горелках; Дат - присосы воздуха в топке;
У° и У0- теоретические объемы соответственно продуктов сгорания и воздуха, м3/м3
3 3 3
(м /кг); ■& - удельный объем пара м /кг; g - водотопливное отношение кгводы/м топлива п впр
(кгводы/кгтоплива).
Подставляя (7), (8) в (6) и сделав ряд преобразований, получим зависимость изменения скорости образования оксидов азота от количества вводимой влаги:
(^^N0^ = 9,15 ■ 1011 ■ , ^ = ■ ехр[-^ ], (9)
ёт ^7+^ $ Тт
1 2 п впр
где L =121 ■ V0■ (а + 0,5 ■ Да -1) -[79 ■ 0- (а + 0,5 ■ Да -1)] , (10)
1 V V гор т V гор т
т = V0 + (а + 0,5 ■ Да -1) ■ V0 , (11)
2 г гор т
134700
L3= -g—. U2)
При расчете уровня максимальных температур Тт с учетом ввода влаги в зону горения использовалась методика, в основу которой положено определение адиабатической температуры в топке и критерия Больцмана.
Адиабатическая температура горения рассчитывается по известной зависимости, когда все тепловыделение в топке воспринимается продуктами сгорания с учетом поправки на дополнительно вводимую влагу.
Получены зависимости адиабатической температуры горения Та , критерия Больцмана и максимальной температуры дымовых газов в ядре факела Тт от условий теплообмена и с учетом влияния ввода влаги.
В итоге получены аналитические зависимости для определения концентраций оксидов азота CgNo2 (г/м3), учитывающая ввод дополнительной влаги в топку котла и коэффициента выхода оксидов азота при сжигании топлива также с учетом ввода влаги.
g 11 L, ■ Ат L, T4 - (T - ДТ) 4 (13)
С No2 = 286-99 ■ 10 ■ , 1 Т 3 ■ exp[ - Г3 ] ■ -т------------------т-^. ( '
(L, + »п ■ g„„„> Тт Т„4,- (Т m - ДТ)
'2 П 6 впр m m
Расчет теоретической зависимости образования оксидов азота от количества вводимой влаги в топку котла
Таблица 1
Расчетные параметры Водотопливное отношение g^, %
G 1G 2G 3G 4G 5G
Максимальная температура факела Тт, 0К 2G1G 1997 1983 1969 1956 1942
Объемная концентрация атомарного кислорода СО2, % об. 1,3170 1,3030 1,2890 1,2750 1,2639 1,2470
Объемная концентрация молекулярного азота См2, % об. 71,602 70,575 69,547 68,520 68,150 66,465
Скорость образования оксидов азота (ёСш2Мт), % об. 0,1577 0,1272 0,0966 0,0661 0,0585 0,0049
Концентрация оксидов азота в топке Січ02 , г/м3 0,2500 0,2040 0,1580 0,1392 0,1204 0,1114
Используя зависимость (13), можно оценить коэффициент выхода оксидов азота при сжигании топлива с вводом влаги:
,,вых 286,99 ■ 1011 L. ■ Ат L Т4 -(Т - ДТ)4
KNO2 =-------g=0-----, 1 Т ■ exp[--Л.] . „ ' m-(I4)
СNO2 L I Тт Т4 ■ (Т - ДТ)4
NO2 ^(l2 +tfn ■ gВпр) т k т '
Аналитическая проверка зависимостей (13) и (14) выполнена для энергетического котла БКЗ-320-140 ГМ, для которого проводились опытные замеры концентраций оксидов азота с впрыском воды в тракт горячего воздуха.
Экспериментальная зависимость изменения концентраций оксидов азота от вводимой влаги в топку котла приведена в табл. 2 (по опытным данным [6]).
Таблица 2
Экспериментальная зависимость изменения концентраций оксидов азота в зависимости от вводимой влаги в топку
Расчетные параметры Водотопливное отношение Йвпр, %
0 10 20 30 40 50
Концентрация оксидов азота в топке Смо2 , г/м3 0,2530 0,2065 0,1600 0,1410 0,1220 0,1130
Как видно из анализа табл.1, 2, получена достаточная сходимость экспериментальных результатов и теоретических данных, рассчитанных по предлагаемой методике. Максимальная величина отклонений расчетной величины концентрации NO2 от экспериментальных значений не превышает 1,5 %.
Для учета реальных условий эксплуатации котлов аналитическая модель дополнена математической моделью расчета процесса горения и расчета топочной камеры.
В математической модели представлен сложный комплекс процессов, происходящих при факельном горении, включая аэродинамику топочной камеры и теплообмен с поверхностями нагрева. В данной модели производится позонный расчет максимальной температуры горения, скорости изменения температуры горения с учетом всех возможных режимов работы топочной камеры, что позволяет проследить влияние темпов охлаждения, начального и конечного состава продуктов сгорания, регулируемых путем ввода активных добавок (воды) с целью снижения концентраций вредных выбросов.
С использованием прикладного языка программирования разработан программный комплекс макрокинетической модели снижения образования выбросов при проектировании природоохранных мероприятий, учитывающий совместно аналитическую (кинетическую) модель расчета образования вредных выбросов в топке котла, а также математическую модель расчета процесса горения и расчета топочной камеры, с помощью которого просто и оперативно, а также с достаточной точностью можно произвести определение выбросов оксидов азота котельными установками.
ЛИТЕРАТУРА
1. Зельдович Я.Б., Садовников П.Я., Франк-Каменецкий Д. А. Окисление азота при горении. М., 1946. 245 с.
2. Сигал И.Я. Защита воздушного бассейна при сжигании топлива. Л.: Недра. 1988. 312 с.
3. Эфендиев Т.Б. Образование оксидов азота в парогенераторах // Теплоэнергетика. 1975. №9. С.20-23.
4. Росляков Л.П. Механизм влияния добавок воды и влагосодержания топлива на образование термических и топливных оксидов азота // Изв. вузов. Энергетика.1988. №7.С.59-63.
5. Рейсиг В.А., Сигал И.Я. Кинетические особенности образования оксидов азота в теплоэнергетических установках // Теплоэнергетика.1993. №1.С.28-31.
6. Шупарский А.И., Голубь Н.В., Ерофеева В.И., Ростунцева И. А. Снижение выбросов оксидов азота при вводе воды в воздуховоды котлов. // Теплоэнергетика .1991. № 8.С.104-107.
BIBLIOGRAPHY
1. Zeldovich J.B., Sadovnikov P.Ya., Frank-Karnenetsky D.A. oxidation of nitrogen at burning. M., 1946.245 with.
2. Sigal I.JA. Protection of air pool at fuel burning. Л: Bowels. 1988. 312 with.
3. Efendiev ^B. Forrnation oxidation of nitrogen in stearn and gas generators // Teploenergeti-ka.1975. №9. С.20-23.
4. Roslyakov L.P. Mechanisrn of influence of additives of water and a fuel rnoisture content on obra-zova-nie therrnal and fuel оксидов nitrogen // Izv. High schools Energetika.1988. №7. С.59-63.
5. Rejsig V. A, Sigal I.JA. Kinetic of feature of forrnation oxidation of nitrogen in heat-power installations // Teploenergetika.1993. №1. С.28-31.
6. Decrease of ernissions oxidation of nitrogen at water input in air lines of coppers / Shuparskij A.I., Golub N.V., Yerofeev V. I, Rostuntseva I.A. // Power systern.1991. № 8. С.104-107.
Ростунцева Ирина Алексеевна -
ассистент кафедры «Тепловые электрические станции» Саратовского государственного технического университета
Rostuntseva Irina Alexeyevna -
assistant of the Department «Thermoelectric power stations» Saratov State Technical University
