CC BY
149
ВЕСТНИК ТОМСКОГО ГОСУДАРСТВЕННОГО УНИВЕРСИТЕТА
2008 Математика и механика № 3(4)
УДК 521.534:536.245.022
А.М. Гришин О МАТЕМАТИЧЕСКОМ МОДЕЛИРОВАНИИ ТОРФЯНЫХ ПОЖАРОВ
Дается краткий обзор работ и общая система уравнений для математического моделирования торфяных пожаров. Предлагается программа дальнейших научных исследований этих пожаров для анализа условий их возникновения и распространения, а также создания новых способов борьбы с ними.
Ключевые слова: торф, торфяной пожар, физическая модель, пиролиз и сушка торфа, горение, процессы переноса, атмосфера, математическое моделирование.
Обзор физического и математического моделирования торфяных пожаров
К отдельному типу лесных пожаров относятся почвенные и торфяные пожары, которые возникают в результате действия очагов горения при низовых лесных пожаров. Под влиянием излучения и конвективного теплового потока от фронта низового лесного пожара, горящих пней и валежника происходит нагревание, высушивание, пиролиз, а затем зажигание подстилки - слоя полуразложившихся листвы, хвои, травы и тонких веточек. Плотность подстилки в абсолютно сухом состоянии изменяется в пределах от 30 до 300 кг/м3 [1 - 5], что выше, чем плотность опада. Влагосодержание подстилки тоже выше, чем у напочвенного покрова, и по данным [1-5] изменяется, как правило, в пределах от 470 до 560 %. Поэтому подстилка высыхает значительно медленнее, чем напочвенный покров, и лишь в очень сухую погоду достигает минимального значения - 14 % [1 - 5]. Если под подстилкой расположен слой торфа, то процесс заглубления горения может продолжаться, и может возникнуть торфяной пожар. Согласно [1 - 5], торфом называется продукт неполного разложения отмершей растительной массы, образующийся в условиях избыточной влажности и недостаточной аэрации. Как правило, залежи торфа (торфяники) имеются на болотах. Торфяники образуются в результате так называемого процесса заболачивания озер на равнинах, покрытых лесом. Например, значительная часть территории Томской области покрыта Васюганскими болотами. Болота имеются и в европейской части России (Московская, Владимирская, Тверская, Рязанская области).
Вопрос о пиролизе торфа, по-видимому, впервые обсуждался в работе [1], где и дана величина энергии активации этой реакции Е = 35,1 кДж/моль, а вопрос об энергетическом балансе беспламенного горения торфяных слоев лесной почвы -в работах [2 - 5].
На основе исследования горения верхних торфянистых слоев лесной почвы в [3, 4] было установлено, что процесс распространения торфяного пожара в основном определяется горением коксика - конденсированного продукта пиролиза торфа. Согласно [3], углерод, содержащийся в торфе, почти полностью сгорает до CO2. Утверждается [3], что одна из причин беспламенного горения торфа заключается в его низкой теплопроводности. Высказана гипотеза о том, что горение торфа при его высоком влагосодержании объясняется перемещением влаги в на-
правлении теплового потока от фронта горения [3]. В работе [4] с использованием термопар и газоанализатора ВТИ-2 экспериментально исследован газообмен при почвенных пожарах и установлено, что беспламенное горение торфа происходит вблизи границы раздела слоев кокса и золы, а скорость горения определяется скоростью диффузии кислорода через слой золы в очаг горения. Через этот же слой происходит выброс газообразных и дисперсных продуктов сгорания (СО2, СО, частицы зол и др.) в атмосферу. Утверждается, что процесс горения торфа происходит с химическим недожогом [4] и при сгорании кокса образуются СО2 и СО, а изменение объема газовой фазы не происходит. К недостатком этой работы следует отнести игнорирование процесса пиролиза торфа и фильтрации газов при горении торфа, а также зависимости параметров состояния реагирующей среды от времени и вертикальной пространственной координаты.
Механизм сгорания торфа при почвенных (подземных пожарах) обсуждается в монографии [5]. В ней утверждается, что слой рыхлой реагирующей среды - торф с малой степенью разложения - горит в беспламенном режиме. Отдельные частицы слоя (стебельки мха и других растений) соприкасаются друг с другом лишь в отдельных точках. При этом общая структура сохраняется и слой реагирующей среды сгорает как целый кусок, покрываясь сверху золой. Скорость горения рыхлого слоя торфа согласно [5] лимитируется скоростью распространения «раскаленных точек» по отдельным частицам торфа. Утверждается, что необходимо учитывать перенос теплоты конвекцией и излучением. Э.В. Коневым на основе закона сохранения энергии предложена полуэмпирическая модель горения рыхлого торфяного слоя, основанная на использовании только одного закона природы -закона сохранения энергии. Наряду с моделью горения «рыхлого слоя» в [5] предложена модель беспламенного горения плотного слоя. В рамках этой модели учтена диффузия кислорода через слой золы и использован закон сохранения массы сгорающего материала в алгебраической форме [5]
ри (1 - т - п) = рЁ га,
где р - плотность сгоревшего материала, рш и у - плотность газообразных продуктов сгоревшего материала и скорость выброса этих продуктов в атмосферу, я -относительная площадь пор в слое золы, т - доля кокса, п - зольность материала, и - скорость продвижения фронта сгорания (тления). В итоге в [5] разработаны две полуэмпирические модели горения торфа. Недостатком моделей, предложенных Э.В. Коневым, является игнорирование закона сохранения количества движения, многокомпонентности среды и многокомпонентной диффузии, отсутствие конкретных формул для скоростей химических реакций и уравнения состояния смеси газов, образующейся в торфе при его сгорании.
Математическая модель пористого реагирующего тела с учетом фильтрации и процессов тепло- и массообмена была предложена в [6 - 8] в связи с исследованием термохимического разрушения тел при их входе в атмосферу с гиперзвуковой скоростью. В дальнейшем в работе [9] на основе обработки результатов наблюдения за реальными торфяными пожарами в Томской области, с использованием результатов [1 - 8] и методов механики сплошных реагирующих сред была предложена новая общая математическая модель торфяных пожаров. При создании этой модели использованы законы сохранения массы, количества движения (с учетом фильтрации) и энергии, а также уравнение состояния для смеси реагирующих газов. Среда считалась двухтемпературной и были учтены пиролиз торфа, много-
компонентная диффузия и фильтрация газообразных продуктов пиролиза и воздуха в порах. В справочнике [10] дана информация о свойствах торфа, которая представляет интерес для математического моделирования торфяных пожаров. В нем изложены итоги экспериментальных исследований различных свойств торфа и предложена одномерная система уравнений, описывающая одновременный перенос теплоты и воды при сушке торфа. Недостатком этой модели является игнорирование переноса влаги вследствие фильтрации. Кроме того, в этой системе уравнений не учитывается многофазность среды и нет выражения для скорости испарения воды. Интересные экспериментальные и теоретические результаты были получены в статье [11]. К недостаткам этой работы следует отнести игнорирование многофазности реагирующей среды при горении торфа и теории фильтрации. В дальнейшем был выполнен цикл работ по математическому моделированию торфяных пожаров [12 - 15], которые подтвердили физические основы математической модели [8]. Кроме того, в работе [13] была сделана попытка создания базы данных для математической модели. Затем в [16] была предложена уточненная математическая модель торфяных пожаров второго поколения, в рамках которой учитываются двухтемпературность пористой среды, частицы пепла, сажи, дыма, свободной воды и влияние многокомпонентности газовой фазы, а в работе
[17] были получены некоторые аналитические решения задач о зажигании слоя торфа.
В статье [18] рассмотрены вопросы глобального мониторинга и оценки экологических последствий природных пожаров.
Физическая модель почвенных и торфяных пожаров
Структура болота представлена на рис. 1, а на рис. 2 предложена схема основных физико-химических процессов, характеризующих механизм распространения торфяного пожара. Она получена в результате анализа экспериментальных и натурных данных о торфяных пожарах [1 - 5].
В дальнейшем эта схема использовалась при создании первого варианта общей математической модели торфяных пожаров [9] и при конкретном математическом моделировании возникновения и распространения торфяных пожаров в работах [12 - 17].
Рис. 1. Структура болота: 1 - подстилка; 2 - слой торфа; 3 - слой воды; 4 - дно водоема; 5 - границы зоны горения торфа; 6 - конвективная колонка над очагом торфяного пожара
Рис. 2. Схема основных физико-химических процессов во фронте торфяного пожара
В работе [7] обсуждались экологические последствия торфяных пожаров. В июле 2002 г. в окрестностях Москвы загрязнение окружающей среды привело к снижению видимости до расстояния менее 100 метров. В воздухе находились дискретные частицы диаметром меньше микрона, формальдегид, полицикличе-ские ароматические углеводороды и окись углерода.
Схема основных физико-химических процессов и экологических последствий торфяных пожаров с учетом физико-химических процессов в зоне торфяного пожара представлена на рис. 3.
Рис. 3. Схема физико-химических процессов и экологических последствий, возникающих при торфяных пожарах
Стрелки показывают последовательность событий и символизируют процессы переноса теплоты и массы компонентов в слое торфа и контрольном объеме приземного слоя атмосферы над торфяником.
В итоге при торфяных пожарах образуется рыхлый слой серого пепла и облако (конвективная колонка) газообразных и дисперсных продуктов неполного сгорания торфа, в том числе оксида углерода и частиц сажи, которые оказывают негативное влияние на людей и животных в зоне действия торфяного пожара. В сентябре 2002 г. были случаи госпитализации некоторых москвичей в результате действия смога, образовавшегося в некоторых районах Московской области [18], так как задымление этой области в результате действия торфяных пожаров, длившееся с конца июля по начало сентября 2002 г., оказало вредное влияние на здоровье москвичей, поскольку дым торфяных пожаров содержит вредные для здоровья людей вещества [18] (частицы сажи и пепла) с диаметром менее микрона, формальдегид, по-лициклические ароматические углеводороды и окись углерода).
Схему процессов, представленную на рис. 3, можно рассматривать как физическую модель торфяного пожара второго поколения. В ее рамках более полно учитываются реальные физико-химические процессы в очаге торфяного пожара и окружающей среде.
Пожары на торфяниках затухают при прекращении доступа кислорода и при сильном увлажнении торфа, в результате чего увеличиваются тепловые потери от фронта торфяного пожара. Таким образом, предельные условия распространения торфяных пожаров имеют вид
се < с*, Ж > Ж*, (1)
где се - равновесная концентрация кислорода в торфяной залежи, а индекс * приписывается критическим значениям концентрации кислорода и влагосодержания, при которых горение торфа прекращается.
Система уравнений для математического моделирования процессов, возникающих в слое торфа при торфяных пожарах
Известно, что наряду с лесными пожарами огромный ущерб окружающей среде и людям приносят торфяные пожары. Для их тушения отсутствуют эффективные методы, и поэтому привлечение больших материальных и людских сил для тушения торфяников в Московской области летом 1972 г. и летом 2002 г. не привело к успеху.
В связи с этим представляет интерес математическое моделирование торфяных пожаров для исследования закономерностей их распространения и определения предельных (критических) условий, при создании которых горение торфа прекращается.
В составе торфа 56 % углерода, 36 % кислорода, 6 % водорода и 2 % азота [3]. Энергетический баланс при горении торфяников обсуждался в работах [12, 13]. Очевидно, что представляют интерес уравнения состояния среды, баланса массы различных компонентов, количества движения и энергии. В отличие от верховых лесных пожаров при горении торфяников в порах торфа, очевидно, не образуются достаточно крупные горящие частицы, а дисперсная среда состоит из частичек дыма, размер которых значительно меньше радиуса пор. Вместе с тем, имея ввиду особенности процессов переноса при пожарах на торфяниках, можно сделать следующие упрощающие предположения:
1. Торфяник является двухтемпературной горючей недеформируемой пористодисперсной средой, состоящей из следующих фаз: ф! - сухое горючее органиче-
ское вещество из лигнина и целлюлозы; ф2 - вода, связанная с этим веществом; ф3 - коксик (из твердого продукта пиролиза торфа); ф4 - пористый слой золы; ф5 -газовая фаза; ф6 - частицы сажи; ф7 - частицы дыма.
2. Скорость и температура дисперсных частиц совпадают со скоростью и температурой несущей газовой фазы в соответствующей точке пространства.
3. Газодисперсная фаза состоит из тех же компонентов, что и в общей физикоматематической модели лесных пожаров [16].
4. Силы инерции в уравнении движения газовой фазы малы по сравнению с силой давления, и уравнения сохранения количества движения можно использовать в форме квадратичного закона Дарси [10].
5. Влиянием обмена энергией между газовой фазы в порах дисперсными частицами дыма можно пренебречь по сравнению с теплообменом со стенками мак-ропор, а последний учитывается выбором коэффициента объемной теплоотдачи.
6. Температура газовой фазы совпадает с температурой дисперсных частиц.
С учетом сделанных допущений и результатов, полученных в [1 - 6, 11, 16], имеем следующую систему уравнений в частных производных для описания физико-химических процессов, протекающих в слое торфа.
Эр дру .
-£+ я = І = !’2’3’ ді дх ,•
Ф ёгайр = -\ к + Р*^=
V + р^, ф = 2 Фі ,
і=5
аТ 7
і=5
(
Цфф і
дТ
дх
і
дТ_
дх ,•
N
І
а =1
дС
а
дх
(2)
(3)
і J
+аV (і; - Т) + срі (Г5 - Т)(1 -ас)+
(4)
+ Ср7 (Ті Т)К2І + Ч3іК3 + Ч52К52 + Ч53^53 + Ч54К54 + Ч2і (К8- К 8+ ) ;
$ а С а
Р = * о
аі
4 дТ д
^ ріф/ Срі ^ « і =1 д дХі
дх,
Р^эфф о
дС
иса
дх,
- Са Є, а = 1,2,...,N ;
+ Чі$ Чгі ^2і + ?3і ^3і +аг(Т Т
р 5фі = Л й —-- л
15 , р2
І = -л
ді
Рз-
дфз
З/
Мс
М,
5ф4
5/
- а Я - ■
- а4R3s Я7 ;
7 7 N г
Е Са = 1 Е Фі = 1 Р = Рг Е
а=1 і =1 а =1
(5)
(6)
(7)
(8) (9)
б - (1 - ас )Л1* + Л2* + "ГГ" Л3 + ЛІ
м,
-Л
(*)
(10)
Здесь г - время, г - радиус-вектор любой точки в рассматриваемом слое торфа; х1
- декартовы координаты этой точки; 2 - массовая скорость образования газодис-
8
Л
персной фазы, к, - коэффициент фильтрации; срі, ср7 , ср - теплоемкости при постоянном давлении отдельных фаз, водяного пара и газообразных продуктов пиролиза; у,, IV | - компоненты и модуль осредненной скорости газа и дисперсных частиц; р - истинная плотность г-й фазы; Я£ = ^8+ - ^ ; ^8—, ^8+'> - массовые скорости конденсации паров и испарения свободной воды в газодисперсной фазе;
Цъ, 9з, - теплоты пиролиза, испарения связанной воды и горения коксика; д52, д53, д54 - теплоты горения С02, Н2, СО; р - давление газа в порах; са - массовая концентрация а-компонента в газодисперсной среде; N - количество компонентов в газодисперсной среде; N. = N - 3 - количество компонентов газовой фазы; Яа - массовая скорость образования а-компонента газовой фазы в результате пиролиза торфа, испарения воды, гетерогенных и гомогенных химических реакций; Я, Я, Я- массовые скорости образования частиц сажи, дыма и
7
капелек воды при конденсации ее паров; рг - X Раг - плотность газовой фазы;
а=1
раг - парциальные плотности компонентов газовой фазы (индекс 1 соответствует кислороду, 2 - СО; 3 - Н2; 4 - СН4; 5 - СО2; 6 - N2 7 - парам воды Н20); саг - раг / рг, а-1,2...; 7 - массовые концентрации для смеси газов;
са = рО, /р, а = 1,2...; 10 - массовые концентрации компонентов газодисперсной смеси, в которую входят все 7 компонентов газовой фазы и частиц пепла; р,0 - парциальные плотности отдельных компонентов газодисперсной фазы
(а = 1,2,_10; первые семь значений р,0 соответствуют компонентам газовой фазы,
а последние 8 - частицам сажи; 9 - частицам дыма; 10 - капелькам воды); р - р5 ф5 + р6 ф6 + р7 ф7 + р8 ф8 - плотность газодисперсной смеси; р5 - р8 - истинные плотности газовой фазы, частиц сажи, частиц дыма и капелек воды в мак-ропорах торфа; X - коэффициент молекулярной теплопроводности газовой фазы; Оа эфф - эффективный коэффициент диффузии а-компонента, в том числе частиц сажи, золы и капелек воды; аг - коэффициент внутреннего теплообмена; Я5а -массовые скорости изменения компонентов газовой фазы (индекс 1 соответствует скорости изменения - О2, 2 соответствует горению СО, 3 - Н2, 4 - СН4); Я& - массовые скорости образования (исчезновения ) вещества конденсированных фаз (индекс 1 соответствует скорости пиролиза торфа, 2 - испарению воды, связанной с торфом, 3 - горению коксика); Т и Т - температуры газовой и конденсированной фаз; X, - коэффициент теплопроводности пористой конденсированной фазы;
- тепловые эффекты химических реакций; к =1, 2, 3, где 1 соответствует реакции пиролиза, 2 - массовой скорости испарения связанной воды, 3 - скорости горения коксика; аС - коксовое число торфа; Я - универсальная газовая постоянная; Мс и М\ - атомарная и молекулярная массы углерода и кислорода; Ма- молекулярная масса а-компонента газовой фазы; как правило, Р = § - ускорение силы тяготения; Хэфф и рэфф - эффективные коэффициенты теплопроводности и вязкости газа; ау = яа, - коэффициент объемного теплообмена в слое торфа; я - удельная поверхность торфа в данной торфяной залежи; а, - коэффициент теплообмена типичного элемента торфа с окружающей средой; в (10) нижний индекс «уя» соответствует реакциям с участием конденсированных веществ; верхний индекс «я» -характеристикам дисперсной фазы.
О математическом моделировании процессов переноса в приземном слое атмосферы в окрестности очага торфяного пожара
Над очагом торфяного пожара возникает конвективная колонка (облако дисперсных и газообразных продуктов горения торфа). Математическое описание параметров состояния этой колонки (скорости, давления, температуры, концентраций компонентов) представляет сложную задачу математической физики. Эта система уравнений может быть получена как частный случай из общей математической модели лесных пожаров [17]. Для упрощения этой системы были сделаны следующие допущения:
1. В результате торфяных пожаров над поверхностью торфяника возникает конвективная колонка из нагретого очагом пожара воздуха, а также продуктов пиролиза и горения торфа.
2. Для упрощения постановки задачи конвективная колонка считается осесимметричной.
3. Соответствующая система уравнений для моделирования процессов переноса в атмосфере над очагом торфяного пожара как частный случай может быть получена из системы уравнений, предложенной в [17] для описания состояния газодисперсной среды, возникающей при лесных пожарах.
4. На границе раздела атмосферы и очага торфяного пожара должны выполняться соотношения, выражающие законы сохранения массы, количества движения и энергии, примеры которых даны в [17, 18].
Кроме того, для математического моделирования такой нестационарной задачи, как зажигание и горение торфа, необходимы начальные условия, аналогичные условиям, представленным в [13].
База данных общей математической модели лесных и торфяных пожаров и методики численного решения задач
В основную систему уравнений (2) - (10) входят теплофизические коэффициенты, скорости сушки и реакций пиролиза и горения и множество эмпирических постоянных, совокупность которых представляет базу данных общей математической модели. Компоненты базы данных определялись, с одной стороны, из литературных источников, а с другой - на основе экспериментальных исследований с помощью решения обратных задач механики реагирующих сред. Очень важно, чтобы все компоненты базы данных были согласованы друг с другом. Это осуществлялось методом их варьирования в пределах их точности и выбора из условия наилучшего согласования соответствующих расчетных и экспериментальных данных. Очень важно, при этом, чтобы решение прямых и обратных задач математической теории лесных пожаров при численном решении определялось с контролируемой точностью, что представляет собой сложную математическую задачу, так как основная система уравнений (2) - (10) нелинейна. Для численного решения прямых задач теории лесных пожаров использовались разработанные в ТГУ итерационно-интерполяционный метод и специальные методики численного счета, основанные на применении метода Патанкара - Сполдинга.
В систему уравнений (2) - (10), представленную выше, входит ряд параметров, которые характерны только для торфяных пожаров (к5, - коэффициенты теп-
лопроводности и диффузии, теплофизические и термокинетические коэффициенты). К сожалению, из-за ограниченности полноценных экспериментальных дан-
ных по возникновению и распространению торфяных лесных пожаров значительной базы данных для математической модели (2) - (10) до сих пор нет. Поэтому расчеты проводились с использованием правдоподобных допущений о порядках величин упомянутых выше величин и коэффициентах тепло- и массообмена, использованных при записи граничных условий. В качестве основных методов решения задач теории торфяных пожаров использовались итерационно-интерполяционный метод и некоторые приближенные аналитические методы.
Некоторые результаты исследования торфяных пожаров
Основные результаты научных исследований, полученных в Томском государственном университете по проблеме моделирования торфяных пожаров, состоят в следующем:
1. Разработаны общие математические модели торфяных лесных пожаров первого и второго поколения, в рамках которых учитываются известные экспериментальные данные по исследованию торфяных пожаров [6, 17].
2. Созданы упрощенные двух- и трехфронтовые модели для одномерного нестационарного распространения процесса горения на торфяниках [17, 18].
3. В результате математического моделирования исследованы основные закономерности физико-химических процессов в пористых реагирующих средах [116] Показано, что при инициировании горения торфяника источником повышенной температуры для одних значений коэффициентов тепло- и массообмена, описывающих взаимодействие слоя горящего торфа с атмосферой и с инертной почвой, на которой расположен этот слой, может возникнуть поверхностное горение торфа вблизи границы раздела «атмосфера - горящий слой торфа», а для других -волна тления уходит в глубь торфяника [18].
4. С использованием упрощенной математической модели торфяных пожаров
[18] была решена задача о стационарном распространении фронта торфяного пожара. Было установлено, что структура фронта торфяного пожара в качественном отношении совпадает со структурой фронта лесного пожара и горение носит диффузионный характер, т.е. лимитируется притоком окислителя. Расчеты показали, что скорость распространения ю сильно зависит от начальной пористости торфа и коэффициента массообмена с окружающей средой и с ростом этих величин ю растет вследствие наличия запаса кислорода в порах и притока 02 во фронт торфяного пожара из приземного слоя атмосферы.
5. Для создания баз данных, характеризующих свойства реагирующих пористых сред, разработаны методы решения обратных задач математической физики.
Выводы и предложения по развитию теории торфяных пожаров
Проблема математического и физического моделирования торфяных пожаров далеко не исчерпана. Представляет интерес исследование следующих принципиальных вопросов:
1. Создание экспериментальных установок и методик для определения теплофизических характеристик для определения коэффициентов теплопроводности и диффузии в слое торфа в обычных условиях и при распространении торфяных пожаров.
2. Создание экспериментальных установок и методик для определения термокинетических постоянных, характеризующих процессы сушки и пиролиза торфа.
3. Создание детерминированно-вероятностной модели для прогноза возникновения торфяных пожаров в конкретных метеорологических и техногенных условиях с использованием геоинформационных технологий.
4. Создание компьютерных программ для прогноза развития процесса горения торфяных залежей, а также экологических последствий горения торфяников с использованием геоинформационных технологий.
5. Исследование предельных условий распространения торфяных пожаров и разработка новых способов и устройств для борьбы с ними.
Решение этих вопросов позволит в первом приближении завершить создание математической теории торфяных пожаров и использовать ее для создания как способов и средств для борьбы с торфяными пожарами, так и прогнозов экологических последствий торфяных пожаров. На основании имеющегося опыта исследования задач теории лесных пожаров [17 - 19] можно утверждать, что для их решения необходимо сочетать экспериментальные и теоретические методы исследования. В качестве методической основы теории торфяных пожаров целесообразно использовать механику многофазных реагирующих сплошных сред и теорию горения и метод математического моделирования с применением технологии распараллеливания вычислительных операций на суперкомпьютерах, что позволит существенно уменьшить затраты машинного времени и получать поля температур и другие характеристики торфяного пожара в режиме, опережающем реальное время развития торфяного пожара. Для сокращения числа эмпирических входных данных, например эмпирических коэффициентов тепло-и массообмена, целесообразно использовать теорию сопряженных задач механики реагирующих сред и методы решения некорректных задач математической физики [20].
Поскольку в рамках общей системы уравнений (2) - (10) учитывались основные физико-химические процессы, определяющие распространение пожара, то с помощью нее можно давать прогноз экологических последствий торфяных пожаров.
Кроме того, в рамках этой модели можно определить поля плотности излучения над очагом торфяного пожара в различные моменты времени, что в принципе позволяет с использованием метода акад. А.Н. Тихонова [21] и математической теории катастроф создать новую методику обнаружения и диагностики торфяных пожаров по данным аэрокосмического мониторинга. Очевидно, что при этом целесообразно использовать современные геоинформационные технологии и спутниковый мониторинг торфяных пожаров в реальном времени для принятия своевременных решений по борьбе с торфяными пожарами.
ЛИТЕРАТУРА
1. Фарберов И.Л. О процессах пиролиза // Теоретический и окислительный пиролиз топлив и высокополимерных материалов: Сб. статей. М.: Наука, 1966. С. 3 - 8.
2. Курбатский Н.П. Исследование количества и свойств лесных горючих материалов // Вопросы лесной пирологии. Красноярск: ИЛиД СО АН СССР, 1970. С. 3 - 5.
3. Гундар С.В. Об энергетическом балансе беспламенного горения органической части почвы при лесных пожарах // Вопросы лесной пирологии. Красноярск, 1974. С. 74 - 82.
4. Гундар С.В. О газообмене при почвенных пожарах // Проблемы лесной пирологии: Сб. статей. Красноярск: Институт леса и древесины им. В.Н. Сукачева СО АН СССР, 1974. С. 137 - 146.
5. Конев Э.В. Физические основы горения растительных материалов. Новосибирск: Наука, 1977. 239 с.
