Тепловое излучение гетерогенных продуктов сгорания при сжигании торфа Текст научной статьи по специальности «Нанотехнологии»

»ЭНЕРГЕТИЧЕСКОЕ« МТАЛЛУРГ^ЕСЖИ:

%%*ХИМИЧЕСКОЕ*%*

МАШИНОСТРОЕНИЕ

УДК [535.233:662.613:662.641]:621.18

ТЕПЛОВОЕ ИЗЛУЧЕНИЕ ГЕТЕРОГЕННЫХ ПРОДУКТОВ СГОРАНИЯ

ПРИ СЖИГАНИИ ТОРФА

В.А. КУЗЬМИН, И.А. ЗАГРАЙ, Р.В. РУКАВИШНИКОВА

Вятский государственный университет, г. Киров

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

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

Введение

Россия обладает 47 % мировых запасов торфа и имеет будущее для решения проблем местной энергетики. По запасам торфа Кировская область занимает второе место в Приволжском федеральном округе, что составляет 25 % (около 950 млн.тонн) от общих запасов [1, 2]. Стратегией социально-экономического развития Кировской области на период до 2020 года торфяная отрасль отнесена к отраслям опережающего роста [3]. Сравнительный анализ по использованию газа, мазута, угля и торфа показывает, что при соответствующей научной и технологической проработке использование торфа в энергетике не только не уступает, но и имеет более низкую себестоимость производимой энергии и сопутствующих продуктов.

По сравнению с другими регионами Кировская область лидирует по количеству котельных на торфе в сегменте ЖКХ, число которых составляет 51 единицу. Потенциал потребления биотоплива коммунальными котельными, с учетом планируемого перевода на это топливо существующих угольных и мазутных котельных, находящихся в территориальной близости к местам производства, оценивается в 1,7 млн. тонн. В настоящее время большинство существующих энергетических установок имеют заниженный КПД из-за слабо проработанной технологии перевода котлов с угля и мазута на торф.

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

Расчет теплового излучения (радиационных характеристик и характеристик излучения) продуктов сгорания является одним из способов оценки неполноты © В.А. Кузьмин, И.А. Заграй, Р.В. Рукавишникова Проблемы энергетики, 2015, № 3-4

сгорания топлива (недожога). В данной работе проводится комплексное опытно-теоретическое исследование теплового излучения гетерогенных продуктов сгорания (ГПС) при сжигании торфа.

Методика определения химического состава, оптических свойств и дисперсности зольных частиц

Исследование характеристик торфяной золы, образующейся после сжигания торфа, проводилось экспериментально на базе лаборатории современных методов физико-химического анализа, входящей в состав научно-образовательного центра «Нанотехнологии» Вятского государственного университета (ВятГУ). Работа по данному направлению является продолжением НИР [4], выполняемой по программе стратегического развития ВятГУ на период 2012-2016 годы.

Объектами исследования являлись образцы торфяной золы, полученные после сжигания топлива (кусковой и фрезерный торф, торфяные брикеты и т.д.) [5]. В качестве примера в работе подробно рассматривается образец золы после сжигания торфяного брикета (низинный торф, месторождение «Дымное», Верхнекамский р-н, пос. Светлополянск, Кировская область) в лабораторных условиях при температуре 600 °С.

Для изучения химического состава торфяной золы использовался рентгенофлуоресцентный спектрофотометр ЕВХ-ЭООНБ фирмы «Shimadzu» (Япония). Действие прибора основано на сборе и последующем анализе спектра, полученного путём воздействия на исследуемый материал рентгеновским излучением.

Программное обеспечение аппарата позволяет не только выявить элементный состав, но также имеет функцию пересчёта элементного состава в оксидную форму. Результаты по содержанию предельных оксидов: БЮ2 - 46,571 %, СаО - 25,293 %, Ре2Оз - 9,576 %, ЛЬОз - 7,181 %, БОз - 6,334 %, Р2О5- 2,039 %, К2О - 1,083 %, ТЮ2 -0,750 %, М§О - 0,633 %, СГ2О3- 0,204 %, 2пО - 0,141 %, БгО - 0,108 %, МпО - 0,081 %.

Полученные нами данные по составу и массовым долям компонентов торфяной золы близки по значениям к результатам из других работ [2].

Определение оптических свойств частиц торфяной золы, образовавшей после сгорания топлива, является важным для расчета тепловых потоков и теплообмена. Комплексный показатель преломления т = щ — п2г определяет оптические свойства вещества частиц и является важнейшим исходным параметром при расчете и моделировании радиационных характеристик. Действительная часть п1 (показатель преломления) определяет скорость распространения излучения в среде. Мнимая часть п2 (показатель поглощения) характеризует затухание амплитуды электромагнитных колебаний в веществе, обусловленное поглощением.

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

На основании собранной информации определялись эквивалентные показатели преломления и поглощения частиц торфяной золы известного химического состава по принципу аддитивности в спектральном интервале=0,4 -10 мкм (рис. 1). В основном оптические свойства золы определяются оптическими свойствами оксида кремния, который составляет около половины состава вещества.

Показатель преломления щ монотонно убывает в спектральном интервале Х=0,5.. .8 мкм, резкий скачок наблюдается при Х=9,2 .. .10 мкм. Показатель поглощения медленно возрастает при длинах волн X от 0,4 до 9,3 мкм, а затем наблюдается спад до Х=10 мкм.

3.6

2.4

1.8

1.2

0.6

щ1 щ2

1 1 - щ 4 *

2 * г * 1 Л

V » г г ✓ ч

4 ^ а * »

* ^ р *

у' г 0 # / "

10

0.1

0.01

0.001

ю-4

4 5 6 X, мкм

10 10

-5

Рис. 1. Оптические свойства частиц торфяной золы

Экспериментальная регистрация размеров частиц торфяной золы первоначально проводилась с помощью сканирующего электронного микроскопа JEOL ^М-6510Ь¥ [7, 8]. При данном виде анализа поверхность объекта сканируется электронным пучком с последующим получением карты рельефа проанализированной зоны. На рис. 2 показана микрофотография, на которой видны формы и размеры исследуемых частиц.

Рис. 2. Микрофотография частиц торфяной золы, полученная с помощью микроскопа ХЕОЬ 1БМ-6510

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

количество n заносились в систему MathCad, имеющую ряд встроенных функций для расчетов числовых статистических характеристик рядов случайных данных. Гистограммы счетного распределения частиц по размерам строились с использованием встроенной в MathCad функции hist, которая автоматически подсчитывает число частиц Дп, попавших в каждый интервал Дx. На рис.3 показано счетное распределение частиц золы торфа по размерам. Несмотря на то, что данный микрофотографический метод обладает хорошей наглядностью в силу визуализации объекта исследований, следует учитывать неточности при обнаружении отдельных мелких частиц, образующих агломераты в результате слипания.

fix), 1/мкм Q

0.08-г

30 40

X, мкм

Рис. 3. Гистограмма и функция распределения частиц по размерам

Для уточнения результатов по размерам частиц применялся лазерный анализатор микрочастиц «Ласка-2К» фирмы «ЛЮМЭКС» (г. Санкт-Петербург, Россия). Определение функции распределения частиц по размерам fx) производилось методом малоуглового светорассеяния [7, 8]. Через кювету с суспензией частиц пропускается световой пучок когерентного излучателя (лазера). С помощью многоэлементного фотоприемника измеряется радиальное распределение интенсивности света /(0), рассеянного ансамблем микрочастиц. Из данных радиального распределения интенсивности вычисляется функция распределения частиц по размерам fx). По способу реализации метод относится к интегрально-оптическим методам.

Графические данные гранулометрического анализа могут быть выражены в форме интегрального и дифференциального гранулометрического состава (Q и f(x) соответственно). Диапазон показаний анализатора составляет 0,5-100 мкм. Полученные данные по размерам частиц торфяной золы представлены на рис. 3. Из рис. 3 видно, что функция распределения имеет двухмодальный вид, но средний размер частиц совпадает с измерениями, полученными с помощью электронного микроскопа. Это говорит о большей чувствительности метода лазерного анализа. Однако недостатком метода является ограниченный диапазон измерений размеров частиц.

Методика математического моделирования радиационных характеристик и характеристик излучения

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

Радиационных характеристики полидисперсных систем сферических частиц конденсированной фазы рассчитывались по теории Ми. Для определения радиационных характеристик индивидуальных частиц (РХИЧ) - сечений ослабления о0сл, поглощения о погл и рассеянияо рас - в качестве исходных данных необходимы дисперсия комплексного показателя преломления т, параметр дифракции р (р = лхА) и диаметр частиц х [9, 10]:

2 2 - пх к ( ) _пхГ _ _

Сосл _ ^ Косл(т,р), страс _ ^ Крас(т,р), стпогл _стосл страс,

где Косл и Крас - факторы эффективности (безразмерные), определяемые по формулам:

2 да 2 да I 2 2\

Косл _— Е(2и + 1) Ке (ап + Ъп ) , Крас (m, р)_— Е(2п + ЩЫ + \Ъп\ ) • Р п_1 Р п_1

Здесь ап и Ъп - амплитуды парциальных волн (коэффициенты Ми). Определение теплового излучения гетерогенных продуктов сгорания определяется эффектами поглощения, рассеяния и излучения. Они, в свою очередь, носят статистический характер и определяются числовой концентрацией частиц N.

Количество частиц конденсата в единице объема продуктов сгорания определяется выражением

N ^рг (Т )ЖТ

(V) (1 _ 7 )рч (Т )ЖЧ

Здесь г - массовая доля конденсата; рч и рг - плотности частиц и газовой фазы (зависят от температуры Т); и Жг - скорости частиц и газа; (V} - средний объем

частиц, который вычисляется по формуле (для полидисперсной системы)

да

..3

V _1 п| х3 / ( х) йх.

6

0

Если среда, на которую падает пучок излучения, содержит облако сферических полидисперсных частиц одинакового состава, то радиационные характеристики единичного объема (РХЕО) - коэффициенты ослабления к^, поглощения а^ и рассеяния р^ - могут быть вычислены по формулам [9, 10, 11]:

да да да

к\ _ N • |сосд (х) / (х) dx, ах _ N • |апогл (х) / (х) йх, _ N • |арас (х) / (х) йх. 0 0 0 Индикатриса рассеяния

да

_ |т0 (х)/(х)йх 0

представлялась в ряд по полиномам Лежандра [11].

Для ГПС при сжигании торфа, как светорассеивающей среды, рассматривалось интегро-дифференциальное уравнение (ИДУ) переноса энергии излучения [9, 10]:

( /\ ^

(ПУ)4(г,П) + кх1х(г,| 1х(г\

(4п)

Г, Г , ПО

й ® + а, 1ах( Г ).

Здесь Q - направление; Iх - спектральная интенсивность излучения, r - координата; ю - телесный угол; Iox - спектральная интенсивность излучения абсолютно черного тела (АЧТ); - рассеяние назад.

*

Спектральный коэффициент поглощения а^ представлялся суммой

*

коэффициентов поглощения конденсированной и газовой фазы: а^ = а-к + аг^ .

Коэффициент поглощения газовой фазы аг^ вычислялся с использованием

информационной системы SPECTRA (http://spectra.iao.ru), включающей банки данных Hitran, Hitemp и др.

*

Коэффициент ослабления единичного объема ГПС находился так: к^ = а^ + вх .

В работе для решения ИДУ использовался метод сферических гармоник в Р3-приближении, который является естественным математическим аппаратом для приближенного решения задач теории переноса. Простота решения задачи для одномерной геометрии позволяет отслеживать и устанавливать основные закономерности процесса теплового излучения в зависимости от исходных факторов и параметров.

Спектральные и интегральные плотности потоков (Fx и F) через единицу площади поверхности, перпендикулярной направлению нормали, определялись по формулам:

<х>

FX = \ 4 (rQ) QndО, F = J Fx dX .

Q 0

Спектральные и интегральные степени черноты (ех и е) находились:

<ю <ю

SX = Fx/FA4T X , s = J FXdX/J FA4T X dX .

0 0

Расчет газообразных продуктов сгорания при сжигании торфа

В данной работе расчет процесса горения проводился по нормативным данным [12] исходя из химического состава горючей и рабочей массы топлива с помощью обобщенных теплотехнических характеристик топлива. При расчете объемов воздуха и продуктов сгорания принимали, что все горючие вещества рабочей массы торфа полностью окисляются с образованием оксидов исключительно с наивысшей степенью окисления по стехиометрическим уравнениям: 1) С + О2 = СО2; 2) S + О2 = S02; 3) H2 + 0,5 О2 = Н2О.

Горючая масса торфа представляла собой смесь органических веществ с количественным составом: Сг - 58,1 %; Ог - 32,6 %; Нг - 5,9 %; Nr - 2,8 %; Sr - 0,6 %. Содержание золы и влаги в топливе, соответственно, Ас = 46,88 % и Wp = 12,53 %. Коэффициент избытка воздуха для данного вида топлива выбран = 1,2. Влагосодержание воздуха (среднее влагосодержание атмосферного воздуха на расчетный период) приняли d„ = 10 г/кг.

В ходе расчета определили рабочую массу топлива с учетом зольности и влажности, получили следующие соотношения: Ар - 6,656 %; Ср - 23,583 %; Ор -13,232 %; Нр - 2,395 %; Nр - 1,137 %; Sр - 0,244 %. Определили расчетную низшую теплоту сгорания топлива QK = 1880 ккал/кг. Теоретически необходимые для горения количество сухого воздуха У°в = 2,299 м3/кг и атмосферного воздуха У°ввл = 2,335 м3/кг.

При этом действительное количество атмосферного воздуха, необходимого для сжигания топлива, составило Va = 2,802 м3/кг.

На основе полученных данных с помощью расчетных формул определили объемы по газовым составляющим: VN2 = 2,223 м3/кг; Vh2o = 0,894 м3/кг; Vco2 = 0,437 м3/кг; V02 = 0,097 м3/кг; Vso2 = 1,705-10-3 м3/кг.

Расчет продуктов сгорания с целью оптимизации производился по формулам в системе MathCad.

Исходные данные для проведения вычислительного эксперимента

С помощью математического моделирования исследовались радиационные характеристики и характеристики излучения ГПС газоторфяной смеси.

Выбранные: давление p = 105 Па, температура продуктов сгорания T = 1273 K и толщина излучающего слоя L = 7,5 м соответствуют наиболее характерным параметрам при сжигании торфа в котле БКЗ-210-140Ф на ТЭЦ-4 г. Кирова.

Массовая доля конденсата z = 0,11054 и состав газовой фазы N2 = 0,53887, H20 = 0,21799, C02 = 0,10866, О 2 = 0,02332, S02 = 0,00041 определены на основании теплотехнических характеристик исходного топлива.

С помощью данных по химическому составу золы и исходя из найденной доли золы в продуктах сгорания определили молярную массу продуктов сгорания ц = 0,0329 кг/моль и среднюю плотность зольных частиц рч = 3,278 г/см3.

Выбранный спектральный диапазон = 0,4 - 10 мкм ЛХ = 0,005 мкм) соответствует максимальной доле излучения, в него попадающей.

Результаты расчетов теплового излучения

Расчет РХИЧ показывает, что эффективные сечения оосл, опогл и орас изменяются не монотонно (рис.4). Максимумы функций оосл и орас соответствуют длинам волн X = 7,7 мкм и 9,25 мкм. Для данной системы частиц ослабление проходящего излучения через среду определяется в большей степени эффектами рассеяния, чем поглощения.

500 400 300

200

100 о

01 23456789 10

X, мкм

Рис. 4. Радиационные характеристики индивидуальных частиц торфяной золы

Величина опогл превышает орас в спектральном интервале X = 8,0 - 8,7 мкм. Результаты расчетов РХЕО показывают, что коэффициенты , а^ и р^ зависят от длины волны и пропорциональны концентрации частиц N (рис.5). На правой шкале рис. 6 приведен суммарный коэффициент поглощения газовой фазы аг^ для

рассматриваемых условий. Наибольший вклад в суммарный агх вносят молекулы С02 (в диапазоне X = 4,2 - 4,6 мкм), Н20 (во всем диапазон^ = 0,4 - 10 мкм) и БО2 (в диапазоне X = 7,1 - 7,7 мкм). Роль компонентов N2 и О2 пренебрежимо мала (величина аг^ находится в районе 10-9 1/мм и меньше). Влияние газовой фазы (а именно С02)

обуславливает ярко выраженный максимум у коэффициентов и а^ в области

X = 4,2 - 4,6 мкм.

к%, а^ , Рх > 1/мм агх , 1/мм

X, мкм

Рис. 5. Радиационные характеристики единичного объема продуктов сгорания при сжигании газоторфяной смеси

Вт/(см2-мкм)

/ 1 у АЧТ

1 V

1' 1 \

/ \ и и — ——

01 2 3 4 5 6 7 В 9 10

X, мкм

Рис. 6. Спектральная плотность потока излучения продуктов сгорания при сжигании газоторфяной смеси. Интегральная плотность потока Е = 7,694 Вт/см2

1

0.8 0.6 0.4 0.2 0

01 23 456 789 10

X, мкм

Рис. 7. Спектральная степень черноты продуктов сгорания при сжигании газоторфяной смеси. Интегральная степень черноты s = 0,5427

На рис. 7, 8 показаны характеристики излучения (спектральная плотность потока и спектральная степень черноты) продуктов сгорания при сжигании газоторфяной смеси. Спектральные распределения величин Fx и s х имеют селективный характер. Спектральная плотность потока исследуемой среды близка к излучению АЧТ в полосах поглощения газовой фазы (рис. 6). Это видно по спектральной степени черноты (рис.8). При X > 5,4 мкм sx имеет значения, превышающие 0,6, а при > 6,5 мкм значения, превышающие 0,8.

Выводы

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

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

Summary

Chemical composition, dispersivity and optical properties of peat ash particles are investigated in this paper. Calculation of gas phase components is based on thermal characteristics of the original fuel. Methodology and results of the calculation of the thermal radiation of heterogeneous combustion products at peat burning are given.

Key words: thermal radiation, combustion products, peat, ash particles, chemical composition, dispersivity, distribution function ofparticle size, optical properties.

Литература

1. Уланов Н.А. Торфяные и выработанные почвы. Киров: ОАО «Дом печати - Вятка», 2005.

287 с.

sx

Л \ и 1г WV* VN

1 1 J 1

1 \ 1

{

А * \)

2. Разработка концепции программы развития торфяной отрасли Кировской области на перспективу до 2013 года (ред. 21.03.2007): Отчет по выполнению I этапа работ / ОАО «ТГК-5»; рук. Э.Ю. Смелов. Киров, 2007. 141 с.

3. Стратегия социально-экономического развития Кировской области на период до 2020 года: принята Постановлением Правительства Кировской области от 06.12.2009 № 33/432. Киров: 2009. 202 с.

4. Разработка технических рекомендаций по повышению эффективности работы теплоэнергетических установок промышленных и жилищно-коммунальных предприятий: отчет о НИР (заключ.) / Вятский гос. ун-т; рук. В. А. Кузьмин; исполн.: Е.И. Маратканова, И.А. Заграй, Н.А. Кутергина, Н.Н. Миклин, Р.В. Некрасова. Киров, 2012. 205 с. № ГР 01201270369.

5. Кузьмин В.А., Заграй И.А., Рукавишникова Р.В. Исследование химического состава частиц торфяного шлама в продуктах сгорания энергетических установок // Материалы Всерос. науч.-практ. конф. «Общество, наука, инновации» / Вят. гос. ун-т; отв. ред. С.Г. Литвинец. Киров, 2014. 1 электрон. опт. диск (CD-ROM). - (ЭТФ. Секция «Физика и теплотехника», ст. 2). С. 1868-1871.

6. Кузьмин В.А., Заграй И.А., Рукавишникова Р.В. Исследование оптических свойств частиц торфяного шлама продуктов сгорания энергетических установок // Материалы Всерос. науч.-практ. конф. «Общество, наука, инновации» / Вят. гос. ун-т; отв. ред. С.Г. Литвинец. Киров, 2014. 1 электрон. опт. диск (CD-ROM). (ЭТФ. Секция «Физика и теплотехника», ст. 3). С.1872-1876.

7. Кузьмин В.А., Заграй И.А., Рукавишникова Р.В., Маратканова Е.И. Исследование оптических свойств и дисперсности частиц продуктов сгорания энергетических установок, работающих на торфе // Материалы Шестой Российской национальной конференции по теплообмену (РНКТ-6) / МЭИ. М., 2014. 1 электрон. опт. диск (CD-ROM). (Секция 9. «Радиационный и сложный теплообмен», ст. 7).

8. Кузьмин В.А., Заграй И.А., Рукавишникова Р.В., Миклин Н.Н. Определение средних размеров и функции распределения частиц по размерам продуктов сгорания энергетических установок, работающих на торфе // Всерос. ежегод. науч.-практ. конф. «Общество, наука, инновации»: сб.материалов / Вят. гос. ун-т; отв. ред. С.Г. Литвинец. Киров, 2014. 1 электрон. опт. диск (CD-ROM). (ЭТФ. Секция «Физика и теплотехника», ст. 4). С.1877-1882.

9. Кузьмин В.А. Тепловое излучение в двигателях и энергетических установках. Киров: ООО «Фирма «Полекс», 2004. 231 с.

10. Кузьмин В.А., Заграй И.А. Радиационные характеристики дисперсных систем энергетических установок // Монография. Германия: LAP Lambert Academic Publishing GmbH & Co. KG - ISBN: 978-3-659-51152-3. 2013. 140 с.

11. Кузьмин В.А., Пяткова (Заграй) И.А. Радиационные характеристики частиц конденсированной фазы гетерогенных продуктов сгорания в энергетических установках // Известия вузов. Проблемы энергетики. 2010. №1-2. С.14-20.

12. Белоусов В.Н., Смородин С.Н., Смирнова О.С. Топливо и теория горения 4.II. Теория горения // учебное пособие Спб: СПбГТУРП, 2011. 142 с.

Поступила в редакцию 16 декабря 2014 г

Кузьмин Владимир Алексеевич - д-р техн. наук, профессор кафедры «Физика» Вятского государственного университета (ВятГУ), г. Киров. Тел: 8(912)3734010. E-mail: [email protected].

Заграй Ираида Александровна - канд. техн. наук, докторант, доцент кафедры «Физика» Вятского государственного университета (ВятГУ), г. Киров. Тел: 8(905)8703213. E-mail: [email protected].

Рукавишникова Рената Валерьевна - аспирантка кафедры «Физика» Вятского государственного университета (ВятГУ), г. Киров. Тел: 8(909)1306055. E-mail: [email protected].