ГОСУДАРСТВЕННЫЙ КОМИТЕТ РОССИЙСКОЙ ФЕДЕРАЦИИ
ПО ОХРАНЕ ОКРУЖАЮЩЕЙ СРЕДЫ
ПРИКАЗ
от 5 марта 1997 г. N 90
ОБ УТВЕРЖДЕНИИ МЕТОДИК РАСЧЕТА ВЫБРОСОВ ЗАГРЯЗНЯЮЩИХ
ВЕЩЕСТВ В АТМОСФЕРУ
С целью обеспечения единого подхода к расчету выбросов
загрязняющих веществ в атмосферный воздух приказываю:
Утвердить Методику расчета выбросов от источников горения при
разливе нефти и нефтепродуктов (Приложение 1), Методику
определения и расчета выбросов загрязняющих веществ от лесных
пожаров (Приложение 2).
Председатель Государственного
комитета Российской Федерации
по охране окружающей среды
В.И.ДАНИЛОВ-ДАНИЛЬЯН
Приложение 1
к Приказу Государственного
комитета Российской Федерации
по охране окружающей среды
от 5 марта 1997 г. N 90
МЕТОДИКА
РАСЧЕТА ВЫБРОСОВ ОТ ИСТОЧНИКОВ ГОРЕНИЯ ПРИ РАЗЛИВЕ
НЕФТИ И НЕФТЕПРОДУКТОВ
Область применения
Настоящая Методика устанавливает общие требования к расчету
выбросов вредных газообразных и дисперсных веществ в атмосферу при
горении нефти и нефтепродуктов, разлитых на различных типах
подстилающих поверхностей (вода; инертная почва с буграми и
впадинами; почва, покрытая растительностью, в том числе и лесной;
болото).
Настоящая Методика может использоваться для определения
экологического ущерба в результате неконтролируемого горения нефти
и нефтепродуктов, разлитых на различных типах подстилающих
поверхностей.
1. Нормативные ссылки
В настоящем документе использованы ссылки на следующие
официальные издания:
ГОСТ 8.310-78 ГСИ Государственная служба стандартных
справочных данных. Основные положения.
ГОСТ 8.417-81 ГСИ Единицы физических величин.
ГОСТ 2.321-84 ЕСКД Обозначения буквенные.
2. Общие положения
2.1. Цель Методики - расчет итоговых и текущих значений масс
выбросов загрязняющих веществ и тепла в атмосферу при
неконтролируемом горении нефти и нефтепродуктов в открытом
пространстве на различных типах подстилающей поверхности.
2.2. При создании настоящего документа использовались
численные значения коэффициентов генерации поллютантов, полученных
различными авторами на основе лабораторных экспериментов и анализа
последствий реальных пожаров. В том случае, если экспериментальные
данные отсутствовали, использовались экспертные оценки.
2.3. Горение представляет собой быстро протекающее химическое
превращение, сопровождающееся выделением теплоты, света и вредных
веществ в атмосферу. Различают организованное (контролируемое)
горение в топках паровых котлов и различных двигателях и
неконтролируемое горение. При организованном горении соединяется
расчетное количество Н и НП и кислорода. Неконтролируемое горение
имеет место при пожарах в открытом пространстве, возникающих в
результате аварий на нефтебазах и нефтехимических производствах
или на трубопроводах. Оно представляет собой сложный физико -
химический процесс, на скорость которого влияет не только
химическая реакция, но и неконтролируемый приток окислителя из
окружающей среды.
В результате неконтролируемого горения разлитой нефти и
нефтепродуктов возникает конвективная колонка - струя нагретых
продуктов полного и неполного сгорания топлива, которые
выбрасываются благодаря этой колонке в приземный слой атмосферы.
Высота конвективной колонки тем больше, чем большее количество
тепла выделяется при горении, т.к. основная движущая сила
продуктов сгорания - сила Архимеда. Очаг пожара имеет сложную
структуру и включает в себя зону пиролиза углеводородного топлива,
зону догорания газообразных и конденсированных продуктов пиролиза.
Горение нефти и нефтепродуктов происходит при постоянном давлении
и имеет диффузионный характер, т.е. лимитируется поступлением
кислорода благодаря подсосу воздуха из окружающей среды. Любой
пожар имеет начало, стадию квазистационарного горения и стадию
потухания, когда горение прекращается из-за сгорания разлитой
нефти, в результате чего устанавливается новое термодинамическое
равновесие.
2.4. Исходными данными для Методики называют описание
предварительного обследования зоны аварии, которое включает:
1) карту - схему района аварии;
2) краткое природно-климатическое описание района, в котором
расположен очаг пожара;
3) метеорологические характеристики окружающей среды
(температура, скорость ветра, наличие осадков);
4) тип подстилающей поверхности;
5) запас техногенных и природных горючих материалов;
6) экологическое состояние окружающей среды (фоновые
концентрации загрязняющих веществ).
2.5. При расчете выбросов загрязняющих веществ в атмосферу
необходимо использовать следующий общий алгоритм:
а) обследование места аварии и получение исходных данных
(экологического паспорта зоны аварии), в том числе оценка массы
М(0) нефти и нефтепродуктов, разлитых на поверхности;
б) выбор типа подстилающей поверхности;
в) выбор математической модели для расчета выбросов вредных
веществ в атмосферу;
г) численное решение задачи об определении итоговых выбросов
загрязняющих веществ;
д) численное решение задачи об определении текущих выбросов
загрязняющих веществ, если это необходимо;
е) оценка экологического ущерба от выбросов загрязняющих
веществ.
3. Определения, обозначения и сокращения
В настоящей Методике используются следующие термины с
соответствующими определениями и обозначениями:
3.1. Методика расчета выбросов представляет собой совокупность
определений, физических и математических моделей и алгоритмов,
позволяющих найти массу выбросов вредных веществ в атмосферу в
результате горения нефти и нефтепродуктов для различных исходных
данных.
3.2. Физической моделью явления или среды называется описание
основных причинно-следственных связей, объясняющих исследуемое
явление или структуру среды. Она создается на основе наблюдений,
экспериментальных исследований и известных в естественных науках
законов природы. С использованием физической модели строится
математическая модель выбросов вредных продуктов горения в
атмосферу.
3.3. Математической моделью явления или среды называется
совокупность математических уравнений, адекватно отражающих
физическую модель явления или структуру среды, которые получены на
основе законов природы в результате удержания существенных
(родовых) и отбрасывания несущественных признаков исследуемого
явления или среды, что позволяет понимать сущность явления или
структуру среды и достаточно точно прогнозировать основные
характеристики исследуемого явления.
3.4. Масштаб экологической катастрофы при разливе нефти (Н) и
нефтепродуктов (НП) характеризуется начальной массой М(0) нефти
или нефтепродуктов, оказавшейся выброшенной в результате аварии в
окружающую среду, и площадью территории, покрытой ими - S(0). Эти
величины получают эксперты в результате обследования территории,
на которой произошла экологическая катастрофа.
3.5. Растительными горючими материалами (РГМ) называются
природные углеводородные топлива, к которым относятся тонкие
веточки, хвоинки или листья в кронах деревьев и опавшие на землю,
а также напочвенный покров (трава, кустарники, мох, лишайник).
3.6. Выбросом загрязняющих веществ в атмосферу называется
поступление за определенное время в воздух или образование в нем
физико-химических агентов и веществ, неблагоприятно воздействующих
на людей и окружающую среду. Выброс любого вредного вещества альфа
обозначается М(альфа) и измеряется в единицах массы (г, кг, моль).
3.7. Удельным выбросом (коэффициентом эмиссии) альфа-вещества
называется отношение:
К(альфа) = М(альфа) / М(г), (3.1)
где М(г) - масса сгоревшей нефти или нефтепродукта.
Если известна химическая формула Н и НП, то при организованном
горении можно определить коэффициенты эмиссии с помощью
стехиометрии - науки о количественных соотношениях, в которых
различные вещества вступают друг с другом в химическую реакцию.
Для неконтролируемого горения К(альфа) можно определить только
опытным путем.
3.8. Недожогом Н и НП называется масса несгоревшего в условиях
неконтролируемого горения топлива М(н). Величина М(н) измеряется в
кг. Известно, например, что если неконтролируемое горение имеет
место на водной подстилающей поверхности, то на поверхности воды
остается слой нефти толщиной 2 мм.
3.9. Количество сгоревшего углеводородного топлива М(г) можно
определить, используя закон сохранения массы:
М(г) = М(0) - М(н). (3.2)
3.10. При аварии трубопровода, фонтане нефти на
нефтепромыслах, аварии танкера или поезда с железнодорожными
цистернами Н и НП топливо разливается по поверхности, называемой
подстилающей. Существуют четыре типа подстилающей поверхности:
1) водная поверхность;
2) инертная почва с известной пористостью и проницаемостью;
3) почва, покрытая растительностью, которая, впитывая Н и НП,
сгорает вместе с углеводородным топливом;
4) болота, которые представляют совокупность живых болотных
растений, отмерших растений, а также свободной и связанной в
растениях воды.
3.11. Поллютант - вещество, загрязняющее среду обитания.
Русский синоним этого слова - загрязнитель. Масса выброса
поллютанта альфа-сорта, возникающего при горении Н и НП, на
основании (3.1) и (3.2) определяется по формуле:
М(альфа) = К(альфа) (М(0) - М(н)). (3.3)
Величина М(альфа) измеряется в кг.
3.12. Выбросом тепла в атмосферу называется количество
теплоты Q, выделяющееся при горении массы М(г) Н и НП. Выброс
тепла в атмосферу определяется по формуле:
Q = q (М(0) - М(н)). (3.4)
3.13. М(0) - масса нефти или нефтепродукта, разлитых на
поверхности в результате аварии, т.
3.14. Q - количество тепла, выброшенного в атмосферу в
результате горения нефти или нефтепродуктов, кДж.
3.15. q - тепловой эффект горения нефти или нефтепродукта,
кДж/кг.
3.16. S(0) - площадь территории, покрытая разлитыми нефтью или
нефтепродуктами, кв. м.
3.17. р(0) - плотность нефти или нефтепродукта, кг/куб. м.
3.18. h - толщина слоя нефти или нефтепродукта, разлитых на
поверхности, мм.
3.19. h(*) - величина критической толщины слоя нефти или
нефтепродукта, ниже которой горение прекращается, мм.
3.20. фи - пористость почвы (фи = V(n) / V, V(n) - объем пор в
единице объема пористой среды; V - объем пористой среды).
3.21. М(н) - масса не сгоревшего в результате пожара
нефтепродукта или нефти, кг.
3.22. М(г) - масса сгоревшего в результате пожара
нефтепродукта или нефти, кг.
3.23. К - коэффициент полноты сгорания нефти или
нефтепродукта, определяющий, какая часть исходной массы топлива
сгорела при пожаре:
К = М(г) / М(0). (3.5)
3.24. m - запас лесных горючих материалов, кг/кв. м.
3.25. W - влагосодержание лесных горючих материалов в
процентах.
3.26. W(*) - предельное значение влагосодержания, выше
которого ЛГМ не горит, в процентах.
3.27. К(в) - доля площади болота, занятая водой.
3.28. W' - влагосодержание грунта.
3.29. x, y - координаты контура нефтяного пожара - кривой на
подстилающей поверхности, которая охватывает зону горения и
отделяет ее от остальной части подстилающей поверхности.
3.30. a, b - большая и малая полуоси эллипса, соответствующего
контуру пожара.
3.31. x(0), y(0) - координаты центра эллипса.
3.32. V(е) - скорость ветра, м/с.
3.33. омега(n) - скорость распространения пламени по нормали к
контуру нефтяного пожара, м/с.
3.34. Тэта - угол между скоростью распространения пламени по
нормали к контуру и направлением ветра.
3.35. омега (А) - скорость распространения фронта пожара по
направлению ветра, м/с.
3.36. омега (В) - скорость распространения фронта пожара
против скорости ветра, м/с.
3.37. омега (С) - скорость распространения фронта пожара
перпендикулярно скорости ветра, м/с,
3.38. омега (Z) - линейная скорость послойного сгорания нефти,
м/с.
3.39. омега (Z(0)); омега (Z (~~)) <*> - линейная скорость
послойного сгорания нефти при нулевой и предельно большой скорости
ветра, м/с.
--------------------------------
<*> Здесь и далее бесконечность обозначается ~~.
3.40. t - текущее время горения, с.
3.41. F(*) - площадь нефтяного пятна, кв. м.
3.42. М(*) - масса несгоревших нефтепродуктов в момент времени
окончания горения t(*), кг.
3.43. t(*) - полное время горения нефти или нефтепродуктов, с.
3.44. М(Г(*)) - полная масса нефти, выгоревшая к моменту
времени t(*), кг.
3.45. М(альфа(*)) - полное значение выбросов поллютантов к
моменту времени t(*), кг.
3.46. Q(*) - полное значение выбросов тепла к моменту времени
t(*), кДж.
3.47. h(i) - глубина лунки, м.
3.48. r(i) - радиус основания лунки, м.
3.49. альфа(i) - полураствор угла при вершине лунки, град.
4. Описание методики расчета итоговых выбросов
вредных веществ и тепла в атмосферу при горении
топлива на водной поверхности
Известно, что Н и НП обладают меньшей плотностью, чем вода,
они не растворяются в воде и при аварии растекаются на водной
поверхности. Особенностью горения нефти на водной поверхности
является то, что на ней остается слой топлива h(*), который не
сгорает. Величина h(*) зависит от сорта Н и НП.
Для массы недожога М(н) в этом случае следует использовать
формулу:
М(н) = ро(0) S(0) h(*). (4.1)
Введем коэффициент недожога:
М(н)
К(н) = ----. (4.2)
М(0)
Тогда коэффициент полноты сгорания равен:
М(0) - М(н)
К = ----------- = 1 - К(н). (4.3)
М(0)
Очевидно, что по определению 0 < К(н) < 1 и 0 < К < 1. Для
водной поверхности:
ро(0) S(0) h(*)
К = 1 - ---------------. (4.4)
М(0)
Зная К, получим расчетные формулы для выброса поллютантов и
тепла при горении топлива на водной подстилающей поверхности:
М(альфа) = К К(альфа) М(0), (4.5)
альфа = 1...N; Q = q К М(0).
Конкретные значения К(альфа) приведены в таблице 4.1.
Таблица 4.1.
КОЭФФИЦИЕНТЫ ЭМИССИИ ПОЛЛЮТАНТОВ ПРИ ГОРЕНИИ
Н и НП И ЛГМ
----T--------------T------------------------------------T-----------¬
¦ N ¦ Поллютант ¦ К(альфа) для Н и НП (кг/кг) ¦ К(альфа) ¦
¦п/п¦ ¦ ¦ для ЛГМ ¦
¦ ¦ ¦ ¦ (кг/кг) ¦
¦ ¦ +-----------T------------T-----------+-----------+
¦ ¦ ¦ нефть ¦диз. топливо¦ бензин ¦ ¦
+---+--------------+-----------+------------+-----------+-----------+
¦1. ¦Оксид углерода¦ -2¦ -3 ¦ -1¦ -1¦
¦ ¦CO ¦8.40 x 10 ¦7.06 x 10 ¦3.11 x 10 ¦1.35 x 10 ¦
¦ ¦ ¦ ¦ ¦ ¦ ¦
¦2. ¦Диоксид ¦ ¦ ¦ ¦ -1¦
¦ ¦углерода CO(2)¦1.00 ¦1.00 ¦1.00 ¦1.35 x 10 ¦
¦ ¦ ¦ ¦ ¦ ¦ ¦
¦3. ¦Оксиды азота ¦ -3 ¦ -2 ¦ -2¦ -4¦
¦ ¦NO(x) ¦6.9 x 10 ¦2.61 x 10 ¦1.51 x 10 ¦4.05 x 10 ¦
¦ ¦ ¦ ¦ ¦ ¦ ¦
¦4. ¦Оксиды серы (в¦ ¦ ¦ ¦ ¦
¦ ¦пересчете на ¦ -2¦ -3 ¦ -3¦ -6¦
¦ ¦SO(2)) ¦2.78 x 10 ¦4.71 x 10 ¦1.20 x 10 ¦1.00 x 10 ¦
¦ ¦ ¦ ¦ ¦ ¦ ¦
¦5. ¦Сероводород ¦ -3¦ -3 ¦ -3¦ -6¦
¦ ¦(H(2)S) ¦1.00 x 10 ¦1.00 x 10 ¦1.00 x 10 ¦1.00 x 10 ¦
¦ ¦ ¦ ¦ ¦ ¦ ¦
¦ ¦ ¦ -1¦ -2 ¦ -3¦ -2¦
¦6. ¦Сажа (C) ¦1.70 x 10 ¦1.29 x 10 ¦1.47 x 10 ¦1.10 x 10 ¦
¦ ¦ ¦ ¦ ¦ ¦ ¦
¦7. ¦Синильная ¦ -3¦ -3 ¦ -3¦ -6¦
¦ ¦кислота (HCN) ¦1.00 x 10 ¦1.00 x 10 ¦1.00 x 10 ¦1.00 x 10 ¦
¦ ¦ ¦ ¦ ¦ ¦ ¦
¦8. ¦Дым (ультра- ¦ ¦ ¦ ¦ ¦
¦ ¦дисперсные ¦ ¦ ¦ ¦ ¦
¦ ¦частицы ¦ -6¦ -6 ¦ -6¦ -2¦
¦ ¦SiO(2)) ¦1.00 x 10 ¦1.00 x 10 ¦1.00 x 10 ¦5.50 x 10 ¦
¦ ¦ ¦ ¦ ¦ ¦ ¦
¦9. ¦Формальдегид ¦ -3¦ -3 ¦ -4¦ -6¦
¦ ¦(HCHO) ¦1.00 x 10 ¦1.18 x 10 ¦5.33 x 10 ¦1.00 x 10 ¦
¦ ¦ ¦ ¦ ¦ ¦ ¦
¦10.¦Органические ¦ ¦ ¦ ¦ ¦
¦ ¦кислоты (в ¦ ¦ ¦ ¦ ¦
¦ ¦пересчете на ¦ -2¦ -3 ¦ -4¦ -6¦
¦ ¦CH(3)COOH) ¦1.50 x 10 ¦3.65 x 10 ¦5.33 x 10 ¦1.00 x 10 ¦
L---+--------------+-----------+------------+-----------+------------
5. Описание методики расчета итоговых выбросов вредных
веществ при горении топлива на инертной почве
Инертная почва характеризуется количеством микровпадин и
впитыванием Н и НП в почву в результате фильтрации. Например,
скорость фильтрации для песчаной почвы W(ф) = 0.000065 м/с, а для
глинистой W(ф) = 0.00001 м/с.
Для инертной почвы для определения величины К(н) предлагается
следующая формула:
К(н) = фи W', 0 < фи < 1, 0 < W' < 1. (5.1)
Здесь фи - пористость грунта;
W' - влагосодержание грунта.
При разливе нефти на инертной почве горение сосредотачивается
в микровпадинах. Пусть на обследованной местности обнаружено
n микровпадин. Тогда для каждой отдельной впадины i имеем
следующие формулы для определения выбросов вредных веществ и
тепла:
М(альфа, i) = К(i) К(альфа, i) М(0, i),
Q(i) = q К(i) М(0, i),
альфа = 1...N, i = 1...n, (5.2)
где индекс i приписывается характеристикам i-ой впадины.
В целом для всей территории имеем следующие формулы:
n
М(альфа) = сумма К(i) К(альфа, i) М(0, i),
i=1
n
Q = q сумма К(i) М(0, i). (5.3)
i=1
Здесь и выше коэффициенты полноты сгорания К(i) определяются
эмпирически, в результате обследования местности после пожара.
Следует отметить, что для песчаной почвы значения К(i) при прочих
равных условиях меньше, чем для глинистой, т.к. песчаная почва
имеет большую пористость и проницаемость. Поэтому для нее более
значительная часть нефти впитывается в почву и не сгорает, а для
глинистой и каменистой почвы имеет место обратный эффект.
6. Описание методики расчета итоговых выбросов
вредных веществ при горении топлива на почве,
покрытой растительностью
Если почва покрыта растительностью (трава, лишайник Cladonia,
мох, кустарник) или отмершими частями растений (опад хвойных и
лиственных деревьев, подстилка торф), то наряду с нефтепродуктами
сгорают и растительные материалы.
Для сухой почвы, покрытой напочвенным покровом (трава, опад
хвои, опад листьев и др.), считается, что имеет место полное
сгорание и нефти и растительности. Тогда:
К(н) = 0. (6.1)
Лесные горючие материалы можно рассматривать как своеобразное
углеводородное топливо с известными коэффициентами эмиссии,
элементарным составом, химической формулой и тепловым эффектом.
Слой растительности характеризуется типом, запасом
растительных горючих материалов (РГМ) m (кг/кв. м) и
влагосодержанием W в %. В обычных условиях при влагосодержании
W > W(*), где W(*) - предельное значение влагосодержания, выше
которого ЛГМ не горит. Однако РГМ, смоченный нефтью, бензином или
керосином, сгорает даже при W > W(*). Поэтому для расчета выбросов
поллютантов и тепла при горении ЛГМ необходимо использовать
формулы, аналогичные выражениям 5.2:
(р) n(1) (р) (р)
М(альфа) = сумма К(i) К(альфа, i) S (0,i) m(i),
i=1
n(1) (р) (р)
Q(р) = сумма q(i) К(i) S(0, i) m(i). (6.2)
i=1
Здесь S(0,i) - площади, занятые растительностью, в кв. м;
m(i) - запас ЛГМ на i-ой площади в кг/кв. м;
n(1) - общее количество участков лесной (степной)
растительности на обследуемой территории;
(р)
q(i) - значения тепловых эффектов для разных типов ЛГМ в
кДж/кг.
Если слой растительности сплошным образом покрывает
обследуемую территорию, где имеет место разлив нефти, то
формулы (6.2) можно упростить:
(р) (р) (р)
М(альфа) = К К (альфа) S(0) m,
(р) (р) (р)
Q = q К S(0) m, (6.3)
(р) (р) (р)
где К , К(альфа) , q - средние значения упомянутых выше
коэффициентов;
S(0) - общая площадь леса в кв. м;
m - средний запас ЛГМ в кг/кв. м.
При совместном горении нефтепродуктов и ЛГМ выражения для
выбросов вредных веществ и тепла принимают вид:
(р)
М(альфа, Сигма)= М(альфа) + М(альфа) ,
(р)
альфа = 1...N, Q(Сигма) = Q + Q . (6.4)
Здесь М(альфа, Сигма) и Q(Сигма) - суммарные итоговые выбросы
вредных веществ и теплоты;
а N = N(0) + N(р) - полное количество поллютантов,
соответствующее одновременному горению нефтепродуктов и ЛГМ;
N(0), N(р) - количество поллютантов от горения нефти и
растительности соответственно;
М(альфа) и Q определяются по формуле (5.3).
7. Описание методики расчета выбросов поллютантов при
горении нефти и нефтепродуктов, разлитых на болоте
Для болота коэффициент недожога определяется по формуле:
К(н, в) k(b), M принадлежит B,
К(н) = {
1 - k(b), M не принадлежит B. (7.1)
Здесь k(b) - коэффициент водности болота, т.е. доля его
поверхности, занятая водой;
K(н, в) = р h(*) S(0) / М(0) - коэффициент недожога для водной
поверхности, определяемый по формуле (4.4);
В - область болота, покрытая водой.
Расчеты выбросов вредных веществ и тепла, согласно физико -
математической модели выбросов при неконтролируемом горении нефти
и нефтепродуктов, на болоте предполагают априорное знание величин
удельных выбросов (коэффициентов эмиссии) вредных веществ и
удельных теплот сгорания Н и НП (нефть, бензин, дизельное топливо)
и растительных горючих материалов (РГМ).
В продуктах горения Н и НП, помимо безопасных соединений,
таких как вода и углекислый газ, присутствуют значительные
концентрации вредных веществ. К ним относятся сернистый и серный
ангидриды, сероводород, окислы азота и углерода, сажа, синильная и
органические кислоты, предельные и циклические углеводороды,
формальдегид и 3,4 бенз(а)пирен.
Поскольку предлагаемая методика предназначена для расчета
выбросов при горении Н и НП, разлитых на болоте, то необходимо
знать величину удельных выбросов при горении РГМ. В результате
частичного или полного сгорания РГМ происходит дополнительный
вклад массы поллютантов в сумму выбросов при горении Н и НП.
Продукты горения растительных материалов (мох, торф, трава,
кустарники, а также деревья хвойных и лиственных пород) включают в
себя, помимо продуктов горения, присущих Н и НП, дым - дисперсные
частицы зольного остатка, состоящего, в основном, из окиси
кремния.
Запас РГМ может быть значительным, и, соответственно, выбросы
поллютантов при его горении вносят большой вклад в общие выбросы
загрязняющих веществ в результате горения Н и НП. Поэтому
игнорирование вклада продуктов горения РГМ в общую массу выбросов
при горении нефти и нефтепродуктов приведет к значительному
занижению величины экологического ущерба.
Значения удельных выбросов при горении Н и НП и ЛГМ приведены
в таблице 4.1.
Будем считать для определенности, что k(b) - доля территории
болота, покрытая водой, а (1 - k(b)) - кочками, на которых
произрастает мох сфагнум и другая болотная растительность. Тогда
для расчета выбросов вредных веществ и тепла М(альфа) и Q
необходимо использовать формулы:
М(альфа) = (М(0) - М(н)) К(альфа) k(b) S(0),
Q = q (М(0) - М(н)) k(b) S(0), (7.2)
(р) (p)
а для определения М(альфа) , Q - выражения:
(р) (р) (р)
М(альфа) = К К(альфа) S(0) m (1 - k(b)),
(р) (р) (р)
Q = q К S(0) (1 - k(b)). (7.3)
Здесь и выше S(0) - площадь, на которой разлита нефть на
болоте;
k(b) - доля этой площади, соответствующая водной поверхности;
(р)
К - коэффициент полноты сгорания растительности на болоте;
m - запас растительности в кочках на болоте, кг/кв. м;
(р)
К(альфа) - удельные выбросы вредных веществ при горении
растительности;
М(0) - масса нефти, разлитой на болоте;
М(н) = ро h(*) k(b) S(0) - масса недогоревшей нефти на болоте;
ро - плотность нефти;
h(*) - предельная высота слоя нефти, при которой нефть еще
горит на водной поверхности болота.
Приложение А
(обязательное)
МАТЕМАТИЧЕСКАЯ МОДЕЛЬ ДЛЯ РАСЧЕТА ТЕКУЩИХ ВЫБРОСОВ
ПОЛЛЮТАНТОВ ПРИ ГОРЕНИИ НЕФТИ НА ПОВЕРХНОСТИ ВОДЫ
Известно, что при крупных авариях (горение нефти при аварии
танкеров, горение нефти при ее разливе в результате разрыва
нефтепровода, горение нефти на нефтепромыслах) возникает
необходимость определения выбросов вредных веществ в текущий
момент времени. При организации тушения очага пожара это
необходимо и для прогноза масштабов экологического бедствия и
оценки времени горения. В данном разделе для тех же типов
подстилающих поверхностей даются формулы для изменения массы
выбросов поллютантов при изменении времени. Весь процесс горения и
выбросов поллютантов разбивается на стадию формирования очага
пожара и стадию догорания нефти и нефтепродуктов. В данном разделе
приводятся формулы для расчета выбросов поллютантов типичных
динамических режимов горения нефти и нефтепродуктов на поверхности
воды.
При разливе нефти на водной поверхности образуется нефтяное
пятно, по форме напоминающее эллипс, большая полуось которого
ориентирована в направлении ветра. Поэтому площадь горения
представляет собой площадь эллипса:
2
(x - x(0)) y
---------- + --- = 1, S = пи a b. (А.1)
2 2
a b
Здесь x и y - координаты контура (см. рис. 4);
направление оси x совпадает с направлением ветра v(e);
a = 1/2 (омега(А) + омега(В)) t, b = омега(С) t - большая и
малая полуоси эллипса;
x(0) = 1/2 (омега(А) - омега(В)) t, y = 0 - координаты центра
эллипса;
омега(А) - скорость распространения фронта горения по
направлению ветра;
омега(В) - скорость распространения фронта пожара против
скорости ветра;
омега(С) - скорость распространения горения перпендикулярно
скорости ветра.
Вид эллипса, отграничивающего площадь горения, представлен на
рис. А.1, где:
0 - точка, где произошло зажигание нефти или нефтепродуктов;
кси, эта - подвижная система координат, связанная с центром
эллипса (x, 0);
М - любая точка контура;
А, В, С, D - характерные точки контура нефтяного пожара.
В табл. А. 1 дана зависимость омега(n) от скорости ветра. Эти
данные для нефти можно аппроксимировать формулой:
омега(n) = омега(0) + (омега(~~) - омега(0)) [1 -
- exp (-бэта V(e) cos Тэта)], (А.2)
где омега(0) = 0.025 м/с, омега(~~) = 0.61 м/с,
бэта = 0.045 - эмпирические постоянные.
В общем случае величины омега(0), омега(~~), бэта зависят от
сорта нефти и вида нефтепродукта, от влагосодержания W горючего и
от угла Тэта между скоростью распространения пламени по нормали к
контуру и вектором скорости ветра V(е). Следуя аналогии между
распространением нефтяного пожара и лесного пожара, будем считать,
что при Тэта = пи/2 и 3 пи/2 величина омега(n) = омега(0).
Рис. А.1. Вид контура, отграничивающего площадь
горения нефти или нефтепродуктов, разлитых
на водной поверхности. <*>
--------------------------------
<*> Рисунок не приводится.
Наряду с распространением пламени по площади пятна разлитой
нефти имеет место ее послойное вертикальное сгорание. Линейную
скорость послойного сгорания можно рассчитать по формуле:
омега(z) = омега(z 0) + (омега(z ~~) - омега(z 0)) [1 -
- exp (-бэта(z) V(e))], (A.3)
где омега(z 0) и омега(z ~~) - линейные скорости послойного
сгорания нефти и нефтепродуктов при нулевой и предельно большой
скорости ветра, а бэта - эмпирическая постоянная.
В табл. А.1 даны значения упомянутых выше эмпирических
постоянных.
Таблица А.1
ВЛИЯНИЕ СКОРОСТИ ВЕТРА НА ЛИНЕЙНУЮ СКОРОСТЬ ПОСЛОЙНОГО
ВЫГОРАНИЯ НЕФТИ И НЕФТЕПРОДУКТОВ
---------------------T--------------T--------------T-------------¬
¦ Вид горючего ¦ бэта(z), ¦ омега(z 0), ¦ омега(z ~~),¦
¦ ¦ сек./м ¦ м/с ¦ м/с ¦
+--------------------+--------------+--------------+-------------+
¦Дизельное топливо ¦ 0.1 ¦ 0.000047 ¦ 0.0001762 ¦
¦ ¦ ¦5 ¦ ¦
¦ ¦ ¦ ¦ ¦
¦Тракторный керосин ¦ 0.1 ¦ 0.000033 ¦ 0.0001299 ¦
¦ ¦ ¦3 ¦ ¦
¦ ¦ ¦ ¦ ¦
¦Автомобильный бензин¦ 0.1 ¦ 0.000065 ¦ 0.0002925 ¦
¦ ¦ ¦0 ¦ ¦
L--------------------+--------------+--------------+--------------
Для мазута линейная скорость выгорания не меняется при
изменении скорости ветра от 0 до 2.8 м/с и составляет
0.0000166 м/с.
Скорость послойного горения сильно зависит от влажности нефти
или нефтепродукта. В частности, согласно, мазут сгорает только
тогда, когда его влажность не превышает 0.7%, а нефть - при
влажности, не превышающей 10%.
Следует отметить, что скорость распространения нефтяного
пожара омега(n) значительно превышает скорость послойного сгорания
омега(z). Например, при скорости ветра 2 м/с для нефти омега /
омега(z) = 9. Скорость распространения омега известна только для
нефти. Однако известно, что линейная скорость горения
нефтепродуктов растет с убыванием массы нефтепродукта. Поэтому
можно для определения омега(0) и омега(~~) в формуле (А.2)
использовать простые соотношения:
омега(0 Н) ро(Н)
омега(0 i) = ----------------,
ро(i) (А.4)
омега(~~ Н) ро(Н)
омега(~~ i) = -----------------,
ро(i)
где индексы Н и i соответствуют нефти и рассматриваемому
нефтепродукту.
Для определения омега(0 z) и омега(~~ z) имеем аналогичные
формулы:
омега( 0 z Н) ро(Н)
омега(z 0 i) = -------------------,
ро(i) (А.5)
омега(~~ z Н) ро(Н)
омега(~~ z i) = -------------------.
ро(i)
Зная площадь горения и скорость послойного сгорания, легко
определяем массовую скорость сгорания нефти или нефтепродукта:
dМ ¦
-- = -пи a b омега(z) ро(s), М¦ = М(0), (А.6)
dt ¦t=0
где ро(s) - плотность нефти или нефтепродукта, а М(0) -
величина выброса нефти или нефтепродукта на подстилающую водную
поверхность (количество разлитой нефти или нефтепродукта).
Уравнение (А.6) представляет собой закон сохранения массы
разлитой нефти или нефтепродукта при их горении на водной
подстилающей поверхности. Решая это уравнение, найдем массу М(г)
сгоревшей к моменту времени t нефти или нефтепродукта:
t
М(г) = М(0) - М = пи ро(s) интеграл a b омега(z) dt. (А.7)
0
Зная текущие значения М(г) и коэффициенты эмиссии К(альфа)
различных поллютантов при горении нефти и нефтепродуктов, легко
находим текущий выброс каждого поллютанта в отдельности и выброс
тепла Q в атмосферу:
М(альфа) = М(г) К(альфа), Q = q М(г), (А.8)
где q - тепловой эффект горения нефти или нефтепродукта.
Надо сказать, что площадь горения не может превысить площади
нефтяного пятна. Совокупность формул (А.1) - (А.8) является
упрощенной математической моделью для расчета выбросов поллютантов
и тепловой энергии в атмосферу.
Поскольку ро(s), омега(А), омега(В), омега(С), омега(z) можно
с удовлетворительной степенью точности считать постоянными
величинами, то величина М(г) легко определяется в результате
интегрирования правой части (А.6):
М(г) = М(0) - М = 1/6 пи ро(s) омега(z) (омега(А) +
3
+ омега(В)) омега(С) t . (А.9)
С помощью формулы (А.9) и таблиц значений К(альфа) находят
выбросы для любых поллютантов и теплоты в любой момент времени:
К(альфа)
М(альфа) = -------- пи ро(s) омега(z) (омега(А) +
6
3
+ омега(В)) омега(С) t ,
q 3
Q = --- пи ро(s) омега(z) (омега(А) + омега(В)) омега(С) t .
6
(А.10)
Легко видеть, что М(г), М(альфа) и Q растут с ростом омега(А),
омега(В), омега(С), омега(z) по линейному закону, а с ростом
времени - по степенному (кубическому) закону.
Следует отметить, что кубический закон горения реализуется
только для первых моментов времени (t < t(1)), когда имеет место
распространение пламени по поверхности разлитой нефти. Значение
времени t(1), при котором пламя охватывает всю площадь разлитой
нефти, легко получается, если известна площадь нефтяного пятна
F(*):
------------------------------------
¦ 2 F(*)
t(1) = v (---------------------------------). (А.11)
пи (омега(А) + омега(В)) омега(С)
При выводе этой формулы, так же, как и при выводе формул
(А.9) - (А.10), считалось, что омега(А), омега(В), омега(С) не
зависят от времени.
Подставляя (А.11) в (А.9), (А.10), получим значения:
М(1г) = М(г) (t(1)) = 1/3 ро(s) омега(z) t(1) F(*), (А.12)
М(1 альфа) = М(альфа) (t(1)) =
= 1/3 К(альфа) ро(s) омега(z) t(1) F(*), (А.13)
Q(1) = Q (t(1)) = 1/3 q ро(s) омега(z) t(1) F(*). (А.14)
При (t > t(1)) распространение фронта нефтяного пожара
прекращается и имеет место вторая стадия процесса - режим
догорания нефтяной пленки с послойной скоростью сгорания омега(z).
Известно, что нефть и нефтепродукты не сгорают полностью и на
водной поверхности остается пленка толщиной h = h(*) = 2 мм. Пусть
площадь пятна в момент окончания горения известна и равна F(*),
тогда масса несгоревшего горючего М(*) равна:
М(*) = ро(s) h(*) F(*). (А.15)
Полная масса топлива, сгоревшая к моменту t = t(*) прекращения
горения, равна:
М(г) (t(*)) = М(0) - М(*). (А.16)
Интегрируя (А.6) с учетом того, что при t > t(1) величина
пи a b = F(*) = const, получаем, используя (А.12) в качестве
начального условия:
М(г) = М(1 г) + ро(s) омега(z) F(*) (t - t(1)),
t >= t(1). (А.17)
Выражения для выбросов поллютантов и тепла при t > t(1) имеют
вид:
М(1 альфа) = [М(1 г) + ро(s) омега(z) F(*) (t -
- t(1))] К(альфа), (А.18)
Q(альфа) = [М(1 г) + ро(s) омега(z) F(*) (t -
- t(1))] q. (А.19)
Величина времени, в течение которого длится второй период
горения, легко определяется по формуле:
М(0) - М(1 г) - М(*)
t(2) = --------------------. (А.20)
ро(s) омега(z) F(*)
Полное время горения определяется по формуле:
t(*) = t(1) + t(2). (А.21)
Подставляя (А.11) и (А.20) в формулу (А.21), получаем полное
время горения:
------------------------------------
¦ 2 F(*)
t(*) = v (--------------------------------- +
пи (омега(А) + омега(В)) омега(С)
М(0) - М(1 г) - М(*)
+ --------------------. (А.22)
ро(s) омега(z) F(*)
По определению М(г *) имеем следующее выражение для полной
массы нефти, выгоревшей к моменту t(*):
М(г *) = М(0) - М(*) = М(0) - ро(s) h(*) F(*). (А.23)
Зная М(г *), легко находим полное значение выбросов
поллютантов и теплоты к моменту прекращения горения:
М(альфа *) = (М(0) - М(*)) К(альфа) = (М(0) -
- ро(s) h(*) F(*)) К(альфа), (А.24)
Q(*) = (М(0) - ро(s) h(*) F(*)) q. (А.25)
При подсчете Q(*) необходимо знать величину q. Для ее
определения следует использовать формулу Д.И. Менделеева:
р р р р р
q = 4.19 [81 C - 246 H - 26 (O - S ) - 6 W ], (А.26)
где q - низшая теплота сгорания в кДж/кг;
р р р р р
величины C , H , O , S , W характеризуют элементарное
содержание в % углерода, водорода, кислорода, серы и влаги;
верхний индекс р означает, что проценты берутся к рабочей
массе горючего.
Таким образом, аналитические формулы (А.10), (А.18) - (А.19)
позволяют получить конкретные текущие значения выбросов
поллютантов и теплоты в приземный слой атмосферы, если известны
М(0), h(*), F(*), ро(s), К(альфа).
Приложение Б
(обязательное)
МАТЕМАТИЧЕСКАЯ МОДЕЛЬ И АЛГОРИТМ РАСЧЕТА
ТЕКУЩИХ ВЫБРОСОВ ПРИ ГОРЕНИИ НЕФТИ И НЕФТЕПРОДУКТОВ
НА ИНЕРТНОЙ ПОЧВЕ
При горении нефти и нефтепродуктов, разлитых на инертной
почве, необходимо учитывать впадины и возвышенности. Очевидно, что
нефть и нефтепродукты сосредотачиваются во впадинах. Пусть
количество впадин равно N, а количество горючего в любой из них -
М(г i). Тогда для любого отдельного очага горения имеем следующее
выражение для выбросов поллютантов и тепла:
М(альфа i) = М(г i) К(альфа), Q(i) = q М(г i). (Б.1)
Общее количество поллютантов и теплоты, выброшенных в
атмосферу:
(г) N
М(альфа) = К(альфа) сумма М(г i),
i=1
(г) N
Q = q сумма М(г i). (Б.2)
i=1
Разлитая на почве нефть скапливается в отдельных лужах
(впадинах) и при зажигании горит в них так же, как и на водной
подстилающей поверхности (п. 4.2). Лунка имеет форму конуса с
высотой h(i), радиусом r(i) и полураствором угла альфа(i) при
вершине. Кроме того, поскольку почва смачивается нефтепродуктами,
то горение слоя нефти в лунке происходит до самого конца.
Рассмотрим горение нефтепродуктов в i-ой лунке.
Скорость их сгорания определяется решением дифференциального
уравнения сохранения массы нефтепродукта:
dМ(i) 2
----- = -ро(i) пи r(i) омега(z). (Б.3)
dt
С начальным условием:
М(i)(0) = М(i 0). (Б.4)
Величину М(i) можно определить по формуле:
3
пи r(i)
М(i) = ро(i) -------- ctg альфа(i). (Б.5)
3
Из (Б.5) находим:
3 М(i) 1/3
r(i) = [---------------------]. (Б.6)
пи ро(i) ctg альфа(i)
Подставляя (Б.6) в (Б.3) и интегрируя по t, получим величину
текущей массы нефтепродукта в i-ой лунке:
1/3
М(i)= [М(i0) -
ро(i) пи омега(zi) 3 2/3 3
- ------------------ (---------------------) t]. (Б.7)
3 пи ро(i) ctg альфа(i)
Очевидно, что время горения в i-ой лунке равно:
1/3 1/3 2/3
3 М(i0) (ctg альфа (i))
t(i г) = ------------------------------. (Б.8)
1/3 1/3
ро(i) пи омега(z i)
Величина массы топлива, сгоревшего к моменту времени t, равна:
1/3
М(г i) = (М(i 0) - М(i)) = М(i 0) - [М(i 0) -
ро(i) пи омега(z i) 3 2/3 3
- ------------------- (---------------------) t] . (Б.9)
3 пи ро(i) ctg альфа(i)
Зная коэффициенты эмиссии К(альфа), получаем выражение для
выброса альфа-поллютанта из i-ой лунки в любой момент времени:
1/3
М(альфа i) = (М(io) - М(i)) К(альфа) = {М(i 0) - [М(i 0) -
ро(i) пи омега(z i) 3 2/3 3
- ------------------ (---------------------) t] } К(альфа).
3 пи ро(i) ctg альфа(i)
(Б.10)
Суммируя выбросы альфа-компонента по всем лункам, получаем:
N N
М(альфа) = Сумма М(альфа i) = Сумма К(альфа) {М(i 0) -
i=1 i=1
1/3 ро(i) пи омега(z i) 3 2/3 3
- [М(i0) - ------------------- (---------------------) t] }.
3 пи ро(i) ctg альфа(i)
(Б.11)
Подставляя в (Б.11) выражение (Б.8), находим суммарный выброс
альфа-компонента:
N
М(альфа *) = Сумма К(альфа) М(i 0). (Б.12)
i=1
Выброс тепла от горения нефтепродуктов во всех лунках
определяется по формуле:
N
Q = Сумма q М(г i). (Б.13)
i=1
К моменту прекращения горения нефтепродуктов в атмосферу
выбрасывается следующее количество тепла:
N
Q(*) = Сумма q М(i 0). (Б.14)
i=1
Таким образом, с помощью формул (Б.11) - (Б.14) удается
получить как выбросы поллютантов, так и выбросы тепла.
Приложение В
(обязательное)
МАТЕМАТИЧЕСКАЯ МОДЕЛЬ И МЕТОДИКА ДЛЯ РАСЧЕТА
ТЕКУЩИХ ВЫБРОСОВ ПОЛЛЮТАНТОВ ПРИ ГОРЕНИИ НЕФТЕПРОДУКТОВ,
РАЗЛИТЫХ НА РАСТИТЕЛЬНОМ ПОКРОВЕ
Если почва покрыта растительностью (трава, лишайник Cladonia,
мох, кустарник) или отмершими частями растений (опад хвойных и
лиственных деревьев), то выбросы поллютантов будут обусловлены не
только горением нефти, но и сгоранием растительного напочвенного
покрова. Поэтому для определения выбросов поллютантов и тепла
имеют место формулы:
(г) (р)
М(альфа) = М(альфа) + М(альфа) ,
(г) (р)
Q = Q + Q . (В.1)
(р) (р)
Здесь М(альфа) , Q(альфа) - выбросы альфа-поллютанта и
тепловой энергии, обусловленные горением растительных горючих
материалов.
Следует считать, что распространение сложного пожара
осуществляется с линейной скоростью, присущей распространению
пламени по разлитому нефтепродукту. Тогда форма контура горения
определяется с использованием формул (А.1), (А.2). Допускается,
что скорости распространения фронта пожара и выгорания
напочвенного покрова, смоченного нефтью или нефтепродуктами,
определяются линейной скоростью распространения и линейной
скоростью выгорания нефтепродукта. Тогда законы сохранения массы
нефтепродукта и лесного горючего материала (ЛГМ) принимают вид:
dМ(2)
----- = -пи a b омега(z) ро(2), М(2)(0) = М(0 2), (В.1)
dt
dМ(1)
----- = -пи a b омега(z) ро(1), М(1)(0) = М(0 1). (В.2)
dt
Здесь М(1) - текущая масса нефти на подстилающей поверхности;
М(2) - текущая масса растительности.
Интегрируя (В.1), получают выражение для масс нефтепродуктов и
ЛГМ, сгоревших к моменту времени t на стадии распространения
пожара:
пи 3
М(г 1) = -- ро(1) омега(z) (омега(А) + омега(В)) омега(С) t ,
6
0 < t < t(1), (В.3)
пи 3
М(г 2) = -- ро(2) омега(z) (омега(А) + омега(В)) омега(С) t ,
6
0 < t < t(1). (В.4)
Выбросы поллютантов для 0 < t < t(1) при горении
нефтепродуктов и ЛГМ определяются по формулам:
пи К(альфа)
М(альфа 1) = ----------- ро(1) омега(z) (омега(А) +
6
3
+ омега(В)) омега(С) t , 0 < t < t(1), (В.5)
пи К(альфа)
М(альфа 2) = ----------- ро(2) омега(z) (омега(А) +
6
3
+ омега(В)) омега(С) t , 0 < t < t(1). (В.6)
Общая масса альфа-поллютанта при 0 < t < t(1):
М(альфа) = М(альфа 1) + М(альфа 2). (В.7)
Общее количество тепла, выбрасываемого в атмосферу при
0 < t < t(1), равно:
пи
Q = -- (q(1) ро(1) + q(2) ро(2)) омега(z) (омега(А) +
6
3
+ омега(В)) омега(С) t . (В.8)
Следует отметить, что так же, как и при горении нефти,
разлитой по водной поверхности, при t > t(1) имеет место режим
догорания. Однако, в отличие от упомянутого выше случая, будем
предполагать, что имеет место полное сгорание ЛГМ, смоченных
нефтепродуктами.
Время t(1) легко определить по формуле (А.11), а для
определения масс нефтепродукта и ЛГМ, сгоревших к моменту времени
t, где t(1) < t < t(*), имеем следующие формулы:
М(1 г) = М(11) + ро(1) омега(z) F(*) (t - t(1)), (В.9)
М(2 г) = М(21) + ро(2) омега(z) F(*) (t - t(1)). (В.10)
Здесь М(11) = М(1 г)(t(1)), М(21) = М(2 г)(t(1)) - массы
нефтепродуктов и ЛГМ, сгоревших к моменту времени t(1).
Величина времени, в течение которого длится второй период
(период догорания), определяется выражением:
М(1 0) + М(2 0) - М(11) - М(21)
t(2) = -------------------------------. (В.11)
(ро(1) + ро(2)) омега(z) F(*)
Выражение для выбросов поллютантов и теплоты для момента
времени во втором периоде горения имеет вид:
М(1 альфа) = М(1 г) К(1 альфа), Q(1) = М(1 г) q(1), (В.12)
М(2 альфа) = М(2 г) К(2 альфа), Q(2) = М(2 г) q(2). (В.13)
Полное время горения определяется по формуле:
t(*) = t(1) + t(2), (В.14)
где t(1), t(2) определяются по формулам (А.11) и (В.11).
С учетом допущения о полном сгорании масс нефтепродуктов и ЛГМ
за время t(*) имеем:
М(1*) = М(10), М(2*) = М(20). (В.15)
Таким образом, с помощью формул (В.6) - (В.15) полностью
решается вопрос о расчете выбросов в любой момент времени на
временном интервале 0 < t < t(*).
Приложение 2
к Приказу Государственного
комитета Российской Федерации
по охране окружающей среды
от 5 марта 1997 г. N 90
МЕТОДИКА
ОПРЕДЕЛЕНИЯ И РАСЧЕТА ВЫБРОСОВ ЗАГРЯЗНЯЮЩИХ
ВЕЩЕСТВ ОТ ЛЕСНЫХ ПОЖАРОВ
Область применения
Настоящая Методика устанавливает общие требования к расчету
выбросов загрязняющих газообразных и дисперсных веществ в
атмосферу при неконтролируемом горении лесных горючих материалов
на лесных пожарах различных типов (низовых, верховых и торфяных).
Настоящая Методика может использоваться для определения
экологического ущерба в результате неконтролируемого горения
лесных горючих материалов в открытом пространстве на различных
типах подстилающей поверхности.
1. Нормативные ссылки
В настоящем документе использованы ссылки на следующие
официальные издания:
ГОСТ 8.310-78 ГСИ Государственная служба стандартных
справочных данных. Основные положения.
ГОСТ 8.417-81 ГСИ Единицы физических величин.
ГОСТ 17.6.1.01-83 Охрана природы. Охрана и защита лесов.
Термины и определения. М.: Изд-во стандартов, 1984.
ГОСТ 2.321-84 ЕСКД Обозначения буквенные.
2. Общие положения
2.1. Цель Методики - расчет итоговых и текущих значений масс
выбросов загрязняющих веществ и тепла в атмосферу при
неконтролируемом горении лесных горючих материалов в открытом
пространстве.
2.2. При ее создании использовались численные значения
коэффициентов генерации поллютантов, полученных различными
авторами на основе лабораторных экспериментов и анализа
последствий реальных пожаров. В тех случаях, если
экспериментальные данные отсутствовали, использовались экспертные
оценки.
2.3. Горение представляет собой быстро протекающее химическое
превращение, сопровождающееся выделением теплоты, света и вредных
веществ в атмосферу. Различают организованное (контролируемое)
горение в топках паровых котлов и различных двигателях и
неконтролируемое горение. При организованном горении соединяется
расчетное количество горючего и кислорода. Неконтролируемое
горение имеет место при пожарах в открытом пространстве,
возникающих в лесах в результате воздействия сухих гроз и
неосторожного обращения человека с огнем (брошенные окурки,
непогашенные костры и т.п.). Оно представляет собой сложный
физико-химический процесс, на скорость которого влияет не только
химическая реакция, но и неконтролируемый приток окислителя из
окружающей среды. В результате неконтролируемого горения лесных
горючих материалов возникает конвективная колонка - струя нагретых
продуктов полного и неполного сгорания топлива, которые
выбрасываются благодаря этой колонке в приземный слой атмосферы.
Высота конвективной колонки тем больше, чем большее количество
тепла выделяется при горении, т.к. основная движущая сила
продуктов сгорания - сила Архимеда.
2.4. Исходными данными для Методики называют описание
предварительного обследования зоны лесного пожара, которое
включает:
1) карту - схему района пожара;
2) краткое природно-климатическое описание района, в котором
расположен очаг пожара;
3) метеорологические характеристики окружающей среды
(температура, скорость ветра, наличие осадков);
4) запас лесных горючих материалов (ЛГМ);
5) экологическое состояние окружающей среды (фоновые
концентрации загрязняющих веществ).
2.5. При расчете выбросов загрязняющих веществ в атмосферу
необходимо использовать следующий общий алгоритм:
а) обследование зоны лесного пожара и получение исходных
данных (экологического паспорта лесной территории в зоне пожара),
в том числе оценка первоначальной массы М(0) лесных горючих
материалов соответствующей площади лесной территории, пройденной
огнем;
б) заполнение карточки учета лесного пожара;
в) выбор математической модели для расчета выбросов
загрязняющих веществ в атмосферу;
г) численное решение задачи об определении итоговых выбросов
загрязняющих веществ;
д) численное решение задачи об определении текущих выбросов
загрязняющих веществ, если это необходимо;
е) оценка экологического ущерба от выбросов загрязняющих
веществ.
3. Определения, обозначения и сокращения
В настоящей Методике используются следующие термины с
соответствующими определениями и обозначениями:
3.1. Лесными горючими материалами (ЛГМ) называются природные
углеводородные топлива, к которым относятся тонкие веточки,
хвоинки или листья в кронах деревьев и опавшие на землю, а также
напочвенный покров (трава, кустарники, мох, лишайник), болотные
растения и торф.
3.2. Лесным пожаром называется явление неуправляемого
многостадийного горения в открытом пространстве на покрытой лесом
площади, в рамках которого имеют место взаимосвязанные процессы
конвективного и радиационного переноса энергии, нагревания, сушки
и пиролиза лесных горючих материалов (ЛГМ), а также горение
газообразных и догорание конденсированных продуктов пиролиза ЛГМ.
3.3. В любой момент времени можно выделить на территории,
покрытой лесом, достаточно большой контрольный объем среды - зону
пожара, внутри которой параметры состояния среды в результате
физико-химических превращений, обусловленных лесным пожаром,
отличаются от невозмущенных значений, определяемых погодными
условиями и типом растительности. Схема зоны лесного пожара
представлена на рис. 3.1, 3.2 <*>.
--------------------------------
<*> Рисунки не приводятся.
Схема физико-химических процессов в зоне лесного пожара
представлена на рис. 3.3. Химический состав газообразных и
дисперсных продуктов горения ЛГМ дается в таблице 4.1.
Рис. 3.3. Схема физико-химических процессов в зоне
лесного пожара.
СХЕМА ФИЗИКО-ХИМИЧЕСКИХ ПРОЦЕССОВ В ПРИЗЕМНОМ СЛОЕ
АТМОСФЕРЫ НАД ЗОНОЙ ЛЕСНОГО ПОЖАРА
--------------¬ -------------¬
¦ Конденсация ¦ ¦ Коагуляция ¦
¦ H(2)O ¦ ¦ частиц ¦
L-T------------ L----------T--
¦ ^ Всплытие продуктов горения ^ ¦
¦ ¦ -------------------------¬ ¦ ¦
¦ L-------+ Нагретые газообразные +--------- ¦
¦ ¦ и дисперсные продукты ¦ ¦
¦ L------------------------- ¦
V V
Выпадение Седментация
осадков частиц
СХЕМА ФИЗИКО-ХИМИЧЕСКИХ ПРОЦЕССОВ ВО ФРОНТЕ
ЛЕСНОГО ПОЖАРА
--------------¬ -------------¬
¦ Теплота ¦ ¦ Прогрев ¦
¦сгорания ЛГМ +------------------------------------>¦ ЛГМ ¦
L-------------- L------T------
^ ¦
¦ V
-------+------------¬ ------------¬ --------------¬
¦ Сгорание ¦ ¦ Пиролиз ¦ ¦ Сушка ¦
¦продуктов пиролиза ¦<------+ ЛГМ ¦<---------+ ЛГМ ¦
L-------------------- L------------ L--------------
3.4. Наиболее сильное изменение параметров состояния среды
происходит в некоторой части зоны лесного пожара, называемой
фронтом пожара, который распространяется с некоторой скоростью по
территории, покрытой лесом. Эта скорость определяется процессами
переноса массы и энергии, а также физико-химическими процессами -
сушкой, пиролизом ЛГМ, горением газообразных и конденсированных
продуктов пиролиза ЛГМ. Визуально фронт лесного пожара наблюдается
в виде светящейся зоны лесного пожара.
3.5. Поверхность 1, отделяющая фронт от несгоревших ЛГМ,
называется внешней кромкой фронта лесного пожара. Ее проекцию на
подстилающую поверхность будем называть контуром лесного пожара.
Очевидно, что последний отделяет ЛГМ, пораженные огнем, от
несгоревших материалов. Внешняя кромка, распространяющаяся по
ветру, называется передней, а против ветра - задней кромкой. На
схеме зоны лесного пожара (рис. 3.1 и 3.2) Г - граница зоны
лесного пожара, 1 - контур лесного пожара (рис. 3.1) и внешняя
кромка фронта лесного пожара (рис. 3.2), 2 - внутренняя кромка
фронта пожара (рис. 3.2) и проекция этой кромки на горизонтальную
подстилающую поверхность (рис. 3.1); горизонтальные стрелки на
рис. 3.2 обозначают скорость (по величине и направлению)
невозмущенного ветра для различных высот, а вертикальные w -
скорость вдува газообразных и дисперсных горючих продуктов горения
из фронта пожара в приземный слой атмосферы.
3.6. Скорость распространения лесного пожара направлена по
нормали к контуру лесного пожара, поэтому называется нормальной
скоростью распространения и определяется формулой:
r
дельта(n)
омега(n) = lim ---------, (3.1)
дельта(t)
r
где дельта(n) - приращение контура пожара по нормали за
промежуток времени дельта(t).
3.7. Поверхность 2 (рис. 3.1, 3.2), отделяющая фронт пожара от
сгоревших ЛГМ, называется внутренней кромкой лесного пожара. Она
может продвигаться как по ветру, так и против ветра и,
следовательно, с меньшей скоростью, чем передняя внешняя кромка,
которая распространяется в направлении ветра. В результате толщина
(ширина) фронта для различных пожаров может меняться в зависимости
от времени, скорости ветра и запаса ЛГМ, который мало отличается
от невозмущенных значений.
3.8. В зоне лесного пожара имеет место пористо-дисперсная
среда, которая называется средой пожара и состоит из сухого
органического вещества (объемная доля фи(1)), воды в связанном с
этим веществом состоянии (фи(2)), конденсированного продукта
пиролиза - коксика (фи(3)), конденсированного продукта горения -
золы (фи(4)), газовой фазы (фи(5)), включающей в себя компоненты
воздуха и газообразных продуктов пиролиза (фи(5)), и дисперсных
продуктов горения, которые состоят из частиц сажи (объемная
доля фи(6)), и частиц золы (фи(7)).
3.9. Границей зоны лесного пожара Г называется поверхность, на
которой параметры состояния среды (температура, скорость ветра,
концентрации компонентов) не отличаются от параметров состояния
невозмущенной среды. Очевидно, параметры состояния среды в зоне
пожара асимптотически стремятся к невозмущенным значениям. Поэтому
для конкретного определения границы зоны пожара необходимо априори
определять степень допустимых отличий Е невозмущенных и
возмущенных значений параметров состояния на границе зоны
пожара Г.
3.10. Выбросом загрязняющих веществ в атмосферу называется
поступление за определенное время в воздух или образование в нем
физико-химических агентов и веществ, неблагоприятно воздействующих
на людей и окружающую среду. Выброс любого вредного вещества альфа
обозначается М(альфа) и измеряется в единицах массы (г, кг, моль).
Различают итоговый и текущий выбросы загрязняющих веществ.
Итоговым выбросом называется выброс загрязняющих веществ за все
время горения. Текущим называется выброс, соответствующий
рассматриваемому моменту времени, которое меньше, чем полное время
горения. Поэтому текущий выброс всегда меньше итогового.
3.11. Удельным выбросом (коэффициентом эмиссии) альфа-вещества
при пожаре называется отношение:
m(альфа)
К(альфа) = --------, [К(альфа)] = кг/кг,
m(г)
альфа = 1, 2,..., N, (3.2)
где m(г) - масса ЛГМ на единице площади лесной территории,
сгоревшая при лесном пожаре;
m(альфа) - масса альфа-компонента, образованного при горении
ЛГМ на той же единице площади лесной территории;
индекс альфа изменяется от 1 до N, где N - общее число
поллютантов (вредных веществ), возникающих при лесном пожаре.
Если известна химическая формула ЛГМ, то при организованном
горении можно определить коэффициенты эмиссии с помощью
стехиометрии - науки о количественных соотношениях, в которых
различные вещества вступают друг с другом в химическую реакцию.
Для неконтролируемого горения К(альфа) вместо (3.2) удобнее
определять и опытным путем.
3.12. Недожогом ЛГМ называется масса не сгоревшего в условиях
неконтролируемого горения топлива m(н) на единицу площади. Наряду
с m(н) используется коэффициент недожога, определяемый по формуле:
m(н)
К(н) = ----. (3.3)
m(0)
3.13. m(o) - запас ЛГМ в абсолютно сухом состоянии, кг/кв. м.
3.14. Количество сгоревшего ЛГМ m(г) можно определить,
используя закон сохранения массы:
m(г) = m(0) - m(н), (3.4)
где m(г) измеряется в кг/кв. м.
3.15. Поллютант - вещество, загрязняющее среду обитания.
Русский синоним этого слова - загрязнитель. Масса выброса
поллютанта альфа-сорта, возникающего при горении единицы площади
растительного покрова ЛГМ, на основании (3.3) и (3.4) определяется
по формуле:
m(альфа) = К(альфа) (m(0) - m(н)). (3.5)
Величина m(альфа) измеряется в кг/кв. м.
3.16. Выбросом тепла в атмосферу называется количество
теплоты Q(n), выделяющееся при горении массы m(г) ЛГМ. Выброс
тепла в атмосферу определяется по формуле:
Q(n) = q (m(0) - m(н)). (3.6)
Величина Q(n) измеряется в Дж/кв. м.
3.17. q - тепловой эффект горения ЛГМ, Дж/кг.
3.18. S(i) - площадь лесной территории, пройденная огнем,
кв. м.
3.19. М(н) - масса не сгоревших в результате пожара ЛГМ на
площади S(i), кг.
3.20. М(г) - масса сгоревших в результате пожара ЛГМ на
площади S(i), кг.
3.21. К - коэффициент полноты сгорания, определяющий, какая
часть исходной массы топлива сгорела при пожаре:
m(г)
К = ----, К = 1 - К(н) (3.7)
m(0)
3.22. W - влагосодержание ЛГМ в процентах.
3.23. W(*) - предельное значение влагосодержания, выше
которого ЛГМ не горят, в процентах.
3.24. x, y - координаты контура пожара.
3.25. a, b - большая и малая полуоси эллипса, соответствующего
контуру пожара.
3.26. x(0), y(0) - координаты центра эллипса.
3.27. V(е) - скорость ветра, м/с.
3.28. омега(А) - скорость распространения фронта пожара по
направлению ветра, м/с.
3.29. омега(В) - скорость распространения фронта пожара против
скорости ветра, м/с.
3.30. омега(С) - скорость распространения фронта пожара
перпендикулярно скорости ветра, м/с.
3.31. t - текущее время горения, с.
3.32. t(*) - полное время горения ЛГМ, с.
3.33. М(альфа*) - полное значение выбросов поллютантов к
моменту времени t(*), кг.
3.34. Q(*) - полное значение выбросов тепла к моменту
времени t(*), Дж.
3.35. S(г) - открытая площадь горения торфяной залежи, через
которую происходит выброс поллютантов в атмосферу, кв. м.
4. Описание методики расчета итоговых выбросов вредных
веществ и тепла при распространении лесных пожаров
Работниками региональных авиационных баз охраны и защиты лесов
и региональных управлений лесами Федеральной лесной службы России
после каждого лесного пожара заполняется карточка учета лесного
пожара (Приложение Г). В ней указывается тип лесного пожара
(низовой, верховой или подземный), выгоревшая площадь и ряд других
данных, которые носят служебный характер.
Величина К, как следует из второй формулы (3.7), выражается
через коэффициент недожога К(н), который может быть определен в
результате осмотра последствий лесного пожара.
Имея карточку лесного пожара и зная К(альфа) и К, итоговый
выброс массы альфа-компонента и тепла Q(n) для любого типа лесного
пожара, следует определять по формуле:
М(альфа i) = S(i) К(i) К(альфа i) m(30 i),
Q(ni) = q(i) k(i) m(30 i),
i = 1, 2, 3. (4.1)
S(i) - площадь лесной территории, пройденная огнем;
Q(ni) - теплота, выделившаяся при пожаре в Дж;
q(i) - тепловой эффект горения ЛГМ в Дж/кг;
индекс i, равный
1 - соответствует параметрам низового лесного пожара;
2 - верховому лесному пожару;
3 - пожару на торфяниках.
Остальные буквенные обозначения имеют тот же смысл, что и в
предыдущих формулах.
Количество поллютантов определяется не только скоростью
горения и размерами очага горения, но и коэффициентами эмиссии
К(альфа) поллютантов, которые зависят от типа растительности и
условий горения. В таблице 4.1 приведены К(альфа).
Таблица 4.1
ОСРЕДНЕННЫЕ ЗНАЧЕНИЯ КОЭФФИЦИЕНТОВ ЭМИССИИ К(АЛЬФА)
-----------------------------------------------T-----------------¬
¦ Наименование поллютанта ¦Значение К(альфа)¦
¦ ¦ кг/кг ¦
+----------------------------------------------+-----------------+
¦Оксид углерода (окись углерода) ¦ 0.135 ¦
¦Углекислый газ ¦ 0.094 ¦
¦Оксид азота ¦ 0.000405 ¦
¦Сажа (элементарный углерод) при горении ¦ 0.0014 ¦
¦Дым (режим горения) ¦ 0.014 ¦
¦Дым (режим тления) ¦ 0.055 ¦
¦Сажа при тлении ¦ 0.011 ¦
¦Метан ¦ 0.075 ¦
¦Другие углеводороды ¦ 0.011 ¦
¦Озон ¦ 0.001 ¦
L----------------------------------------------+------------------
Приложение А
(обязательное)
МАТЕМАТИЧЕСКАЯ МОДЕЛЬ И АЛГОРИТМ РАСЧЕТА ТЕКУЩИХ
ВЫБРОСОВ ПОЛЛЮТАНТОВ В АТМОСФЕРУ ДЛЯ НИЗОВЫХ
ЛЕСНЫХ ПОЖАРОВ
Довольно часто встречается ситуация, когда пожары продолжаются
несколько недель и даже несколько месяцев. В результате возникает
задымленность лесных территорий, по причине которой прекращаются
авиарейсы на местных авиалиниях и плавание судов по рекам. Поэтому
представляет интерес прогноз выбросов поллютантов и тепла для
любых моментов времени.
На основании законов сохранения массы и энергии загрязняющих
компонентов для определения массы выбросов, теплоты и контура
лесного пожара необходимо решать уравнения:
dQ
-- = интеграл q K омега(n) m(3) ds, (А.1)
dt L
dM(альфа)
-------- = интеграл К(альфа) K омега(n) m(3) ds, (А.2)
dt L
dфи
--- + омега(n) ¦grad фи¦ = 0, (А.3)
dt
с соответствующими начальными условиями
М(альфа)(0) = М(альфа 0), Q (0) = Q(0),
фи¦ = фи(0) (x,y). (А.4)
¦t=0
Здесь (А.1) - (А.3) интегро-дифференциальные уравнения для
определения Q(t), М(альфа) (t) и контура лесного пожара
фи = фи(0)(x,y,t). Решение данной системы интегро-дифференциальных
уравнений представляет значительные математические трудности.
В данном Приложении дается упрощенная полуэмпирическая
математическая модель и методика расчета выбросов от низовых
лесных пожаров.
Считается, что очаг лесного пожара представляет собой плоский
источник поллютантов, который увеличивается с ростом времени.
Примем, что контур лесного пожара в любой момент времени в
неподвижной системе координат описывается эллипсом (рис. А.1) <*>:
--------------------------------
<*> Не приводится.
2 2
(x - x(0)) y (омега (А) + омега(В)) t
----------- + -- = 1, а = ------------------------,
2 2 2
а b
(омега(А) - омега(В)) t
b = омега(С) t, X(0) = -----------------------. (А.5)
2
Периметр и площадь этого эллипса определяются по формулам:
----
L = пи [1.5 (a + b) - va b], (А.6)
пи 2
S = пи a b = -- (омега (А) + омега(В)) омега(С) t . (А.7)
2
Известно, что сухие ЛГМ сгорают почти полностью, т.е. К = 1, а
К(н) = 0, в то время как при определенном (предельном)
влагосодержании W = W(*) прогресс горения прекращается, т.е.
К = 0, а К(н) = 1. Исходя из этих физических соображений, будем
считать, что:
W(*) - W
К = --------, (А.8)
W(*)
где W - влагосодержание ЛГМ.
Для низовых лесных пожаров величина W(*) = 0.13, а K примерно
0.5.
Если предположить, что запас m(3) не меняется по периметру
контура лесного пожара, а К(альфа) и скорость распространения не
зависят от времени, то, разбивая контур на N равных частей и
подставляя в (3.7) формулы (А.5) и (А.7), получаем после
интегрирования следующие выражения для массы альфа-загрязняющего
компонента:
(омега(А) + омега(В))
М(альфа) = К К(альфа) m(3) пи {1.5 [--------------------- +
2
-------------------------------
¦ (омега(А) + омега(В)) 2 N
+ омега(С)] - v --------------------- омега(С)} t Сумма омега(ni).
2 i=1
(А.9)
Здесь омега(ni) - значения скорости распространения лесного
пожара, соответствующее i-ой части периметра эллипса (рис. А.1).
Схема разбиения периметра лесного пожара на
элементарные части и направление нормальной
скорости распространения
Рис. А.1 <*>.
--------------------------------
<*> Рисунок не приводится.
При выводе (А.9) было учтено, что процесс распространения
лесного пожара симметричен относительно оси x.
В общем случае, когда m(3) и омега(n) несимметричны
относительно оси x, необходимо использовать формулу:
К(альфа) пи (омега(А) + омега(В))
М(альфа) = К ----------- {1.5 [--------------------- + омега(С)] -
2 2
------------------------------
¦ (омега(А) + омега(В)) 2 N
- v -------------------- омега(С)} t Сумма m(3) омега(ni),
2 i=1 (А.10)
где N - число точек по всему периметру контура лесного пожара.
Легко видеть, что выброс загрязняющего компонента в атмосферу
растет с ростом времени и скорости распространения по
квадратичному закону, а с ростом коэффициентов эмиссии - по
линейному закону. Для определения омега(n) в различных точках
контура низового лесного пожара удобно использовать формулу,
которая с учетом принятых выше обозначений имеет вид:
n
f(1) (ню(е), альфа) (m(3) / m(0))
омега(n) = омега(no) --------------------------------------------
[1 + G(г) (ро - ро(0)) / ро(0)] [1 + С (W -
-----------------------------------. (А.11)
- W(0) + D (T - T(0))] / f(2)(бэта)
-
¦ 1 + В(1) V(е) cos альфа + В(4) V(е),
¦ при 0 є V(е) є V(1);
¦
¦ пи
¦ В(1) (альфа є --) = В(1) (0),
¦ 2
¦
¦ пи
¦ В(1) (альфа > --) = В(1) (пи);
¦ 2
¦
¦ 1 + [В(1) (0) V(е) + В(2) (V(е) -
¦
¦ - V(1))] cos альфа + В(3) (V(е) - V(1)) +
¦
¦ + В(4) V(е);
¦
f(1) (V(е), альфа) ={ при V(1) є V(е) є V(2),
¦
¦ V(1) - V(е)
¦ альфа є альфа(кр) = arccos -----------;
¦ ню'
¦ (А.12)
¦
¦ 1 + В(1) (пи) V(е) cos альфа + В(3) (V(е) -
¦
¦ V(1)) + В(4) V(е);
¦
¦ при V(1) Є V(е) Є V(2), альфа Є альфа(кр);
¦
¦ b
¦ альфа (V(е) cos альфа) , при V(е) > V(2),
¦
¦ V(2)
¦ альфа є arccos ----;
L V(е)
-
¦ 1 + В(1) V(е) cos альфа + В(4) V(е),
¦ при 0 є V(е) є V(1);
¦
¦ пи
¦ В(1) (альфа є --) = В(1) (0),
¦ 2
¦
¦ пи
¦ В(1) (альфа > --) = В(1) (пи);
¦ 2
¦
¦ 1 + [В(1) (0) V(е) + В(2) (V(е) -
¦
¦ - V(1))] cos альфа + В(3) (V(е) - V(1)) +
¦
¦ + В(4) V(е);
f(1) (V(е), альфа) ={
¦ при V(1) є V(е) є V(2),
¦
¦ V(1) - V(е)
¦ альфа є альфа(кр) = arccos ----------- ;
¦ ню'
¦ (А.13)
¦
¦ 1 + В(1) (пи) V(е) cos альфа + В(3) (V(е) -
¦
¦ V(1)) + В(4) V(е);
¦
¦ при V(1) Є V(е) Є V(2), альфа Є альфа(кр);
¦
¦ b
¦ альфа (V(е) cos альфа) , при V(е) > V(2),
¦
¦ V(2)
¦ альфа є arccos ----;
L V(е)
-
¦ (1 - sin 8°) [1 - sin 2 (бэта - 8°)],
¦
¦ при бэта > 8°;
f(2) (бэта) = {
¦ (1 - sin 8°), при -18° є бэта є 8°; (А.14)
¦
¦ примерно 0,7 при бэта є -18°;
L
где омега(no) - скорость распространения низового пожара на
горизонтальной подстилающей поверхности;
V(е) - скорость ветра на стандартной высоте;
m(з) - запас лесных горючих материалов на подстилающей
поверхности, соответствующий скорости распространения омега =
= омега(no) (таблица А.2);
G(г) = 0.3;
ро(0) - плотность слоя на горизонтальной поверхности,
соответствующая скорости распространения омега = омега(no)
(таблица А.2);
n = (0.2 - 0.35) - эмпирический коэффициент;
ню' - пульсация скорости ветра (ню приблизительно равна
1 м/с);
бэта - угол наклона к горизонту рельефа местности;
альфа - угол между направлением скорости распространения
фронта пожара омега и скоростью ветра;
С, D, омега(no), a, b - эмпирические постоянные;
Т(0), W(0) - стандартные начальные температуры горизонтального
слоя ЛГМ и влагосодержание соответственно.
Скорость ветра в лесу V(е) определяется по формуле:
h(1)
V(е) = Ламбда ---- V(em), (А.15)
_
h(3)
где h(1) - высота слоя проводников горения;
_
h(3) - высота верхней границы полога леса;
Ламбда - эмпирический множитель, который при сомкнутости
полога 0; 0.2; 0.4; 0.6; 0.8; 1 принимает значения 0.77; 0.62;
0.47; 0.32; 0.165; 0.09 соответственно.
_
При отсутствии полога леса вместо h(3) и скорости ветра V(em)
по измерениям на метеостанции для высоты 10 м надо брать
_
h(3) = 2 м и скорость ветра по измерениям на высоте 2 м.
Числовые данные для конкретных типов подстилающей поверхности
приведены в таблице А.2.
Таблица А.2
ДАННЫЕ ДЛЯ ОПРЕДЕЛЕНИЯ СКОРОСТИ РАСПРОСТРАНЕНИЯ
ФРОНТА ЛЕСНОГО ПОЖАРА
----------T-------T-------------------T-------------------------T----------------¬
¦Проводни-¦ оме- ¦ 0 < V(e) < V(1) ¦ V(1) < V(e) < V(2) ¦ V(е) > V(2) ¦
¦ ки ¦га(no),+---------T---------+-------------------------+----------------+
¦горения ¦ м ¦ пи ¦ пи ¦ V(1) - V(e)¦ V(2)¦
¦ ¦ ----- ¦альфає-- ¦альфаЄ-- ¦альфа є arccos-----------¦альфаєarccos----¦
¦ ¦ мин. ¦ 2 ¦ 2 ¦ ню ¦ V(e)¦
¦ ¦ +---------+---------+----T----T----T----T-----+--------T-------+
¦ ¦ ¦В(1) (0) ¦В(1) (пи)¦V(1)¦В(2)¦В(3)¦В(4)¦ V(2)¦ a ¦ b ¦
¦ ¦ ¦ с/м ¦ с/м ¦м/с ¦с/м ¦с/м ¦с/м ¦ м/с ¦ ¦ ¦
+---------+-------+---------+---------+----+----+----+----+-----+--------+-------+
¦ низовые пожары ¦
+---------T-------T---------T---------T----T----T----T----T-----T--------T-------+
¦Отмершая ¦ 0.7 ¦ 0.52 ¦ 0.25 ¦ 1.6¦0.63¦0.25¦0.25¦ 4 ¦ 0.15 ¦ 2.5 ¦
¦трава ¦ ¦ ¦ ¦ ¦ ¦ ¦ ¦ ¦ ¦ ¦
¦(весна) ¦ ¦ ¦ ¦ ¦ ¦ ¦ ¦ ¦ ¦ ¦
+---------+-------+---------+---------+----+----+----+----+-----+--------+-------+
¦Вейник: ¦ ¦ ¦ ¦ ¦ ¦ ¦ ¦ ¦ ¦ ¦
¦ лето ¦ 0.5 ¦ 0.26 ¦ 0.13 ¦ 3.0¦0.37¦0.12¦0.12¦ 8.5 ¦ 0.15 ¦ 2.5 ¦
¦ осень ¦ 0.4 ¦ 0.90 ¦ 0.40 ¦ 1.6¦1.10¦0.40¦0.40¦ ¦ ¦ ¦
+---------+-------+---------+---------+----+----+----+----+-----+--------+-------+
¦Лишайник ¦ 0.22 ¦ 2.40 ¦ 0.80 ¦ 1.8¦1.20¦0.80¦0.80¦ 4 ¦ 0.15 ¦ 2.5 ¦
¦Cladonia ¦ ¦ ¦ ¦ ¦ ¦ ¦ ¦ ¦ ¦ ¦
+---------+-------+---------+---------+----+----+----+----+-----+--------+-------+
¦Мхи: ¦ ¦ ¦ ¦ ¦ ¦ ¦ ¦ ¦ ¦ ¦
¦ мох ¦ 0.24 ¦ 2.00 ¦ 0.75 ¦ 1.5¦1.55¦0.75¦0.75¦ 4 ¦ - ¦ - ¦
¦Шребера ¦ ¦ ¦ ¦ ¦ ¦ ¦ ¦ ¦ ¦ ¦
¦ сфагнум ¦ 0.24 ¦ 2.00 ¦ 0.75 ¦ 1.5¦1.55¦0.75¦0.75¦ ¦ ¦ ¦
+---------+-------+---------+---------+----+----+----+----+-----+--------+-------+
¦Опад ¦ ¦ ¦ ¦ ¦ ¦ ¦ ¦ ¦ ¦ ¦
¦листвы: ¦ ¦ ¦ ¦ ¦ ¦ ¦ ¦ ¦ ¦ ¦
¦ береза ¦ 0.18 ¦ 2.30 ¦ 1.00 ¦ 1.7¦2.00¦1.00¦1.00¦ 4 ¦ - ¦ - ¦
¦ осина ¦ 0.24 ¦ 1.80 ¦ 0.75 ¦ 1.7¦1.45¦0.75¦0.75¦ ¦ ¦ ¦
+---------+-------+---------+---------+----+----+----+----+-----+--------+-------+
¦Опад хвои¦ ¦ ¦ ¦ ¦ ¦ ¦ ¦ ¦ ¦ ¦
¦(сосна ¦ 0.11 ¦ 3.00 ¦ 1.70 ¦ 1.7¦3.50¦1.70¦1.70¦ 4 ¦ - ¦ - ¦
¦обыкно- ¦ ¦ ¦ ¦ ¦ ¦ ¦ ¦ ¦ ¦ ¦
¦венная) ¦ ¦ ¦ ¦ ¦ ¦ ¦ ¦ ¦ ¦ ¦
+---------+-------+---------+---------+----+----+----+----+-----+--------+-------+
¦ верховые пожары ¦
+---------T-------T---------T---------T----T----T----T----T-----T--------T-------+
¦Хвоя ¦ ¦ ¦ ¦ ¦ ¦ ¦ ¦ ¦ ¦ ¦
¦(листья) ¦ ¦ ¦ ¦ ¦ ¦ ¦ ¦ ¦ ¦ ¦
¦и тонкие ¦ ¦ ¦ ¦ ¦ ¦ ¦ ¦ ¦ ¦ ¦
¦веточки в¦ ¦ ¦ ¦ ¦ ¦ ¦ ¦ ¦ ¦ ¦
¦кронах ¦ ¦ ¦ ¦ ¦ ¦ ¦ ¦ ¦ ¦ ¦
¦деревьев ¦ ¦ ¦ ¦ ¦ ¦ ¦ ¦ ¦ ¦ ¦
¦(кустар- ¦ ¦ ¦ ¦ ¦ ¦ ¦ ¦ ¦ ¦ ¦
¦ника) ¦ 1.35 ¦ - ¦ - ¦ - ¦ - ¦ - ¦ - ¦ 4 ¦ 1.00 ¦ 2.5 ¦
L---------+-------+---------+---------+----+----+----+----+-----+--------+--------
Таблица А.2
(Продолжение)
ДАННЫЕ ДЛЯ ОПРЕДЕЛЕНИЯ СКОРОСТИ РАСПРОСТРАНЕНИЯ
ФРОНТА ЛЕСНОГО ПОЖАРА
----------T-----------------T-----T----T----------------------¬
¦Проводни-¦m(3)' є m(3) < ~~¦ро(0)¦ G ¦ 0 є W є W(*) ¦
¦ ки +-----T-----T-----+ ¦ +----T----T----T-------+
¦горения ¦m(3) ¦ n ¦m(3)'¦ кг /¦ ¦W(0)¦ C ¦W(*)¦ D ¦
¦ ¦кг / ¦ ¦кг / ¦ куб.¦ ¦ ¦ ¦ ¦ 3 -1¦
¦ ¦кв. м¦ ¦кв. м¦ м ¦ ¦ ¦ ¦ ¦x10 К ¦
+---------+-----+-----+-----+-----+----+----+----+----+-------+
¦ низовые пожары ¦
+---------T-----T-----T-----T-----T----T----T----T----T-------+
¦Отмершая ¦ ¦ ¦ ¦ ¦ ¦ ¦ ¦ ¦ ¦
¦трава ¦ ¦ ¦ ¦ ¦ ¦ ¦ ¦ ¦ ¦
¦(весна) ¦0.30 ¦ 0.34¦0.15 ¦1 - 5¦0.3 ¦0.13¦13 ¦0.18¦ 0.3 ¦
+---------+-----+-----+-----+-----+----+----+----+----+-------+
¦Вейник: ¦ ¦ ¦ ¦ ¦ ¦ ¦ ¦ ¦ ¦
¦ лето ¦0.30 ¦ 0.20¦ - ¦1 - 5¦0.3 ¦0.18¦ - ¦ - ¦ - ¦
¦ осень ¦0.50 ¦ - ¦0.10 ¦ 1 ¦ ¦0.13¦12.9¦0.25¦ 6.3 ¦
+---------+-----+-----+-----+-----+----+----+----+----+-------+
¦Лишайник ¦0.10 ¦ - ¦ - ¦ 4 ¦ ¦0.10¦ 7.7¦0.2 ¦ 3.7 ¦
¦Cladonia ¦ ¦ ¦ ¦ ¦ ¦ ¦ ¦ ¦ ¦
+---------+-----+-----+-----+-----+----+----+----+----+-------+
¦Мхи: ¦ ¦ ¦ ¦ ¦ ¦ ¦ ¦ ¦ ¦
¦ мох ¦ ¦ ¦ ¦ ¦ ¦ ¦ ¦ ¦ ¦
¦Шребера ¦0.25 ¦ - ¦ - ¦ 12 ¦ ¦0.06¦29 ¦0.08¦ 14 ¦
¦ сфагнум ¦0.20 ¦ - ¦ - ¦ - ¦ ¦ - ¦ - ¦ - ¦ ¦
+---------+-----+-----+-----+-----+----+----+----+----+-------+
¦Опад ¦ ¦ ¦ ¦ ¦ ¦ ¦ ¦ ¦ ¦
¦листвы: ¦ ¦ ¦ ¦ ¦ ¦ ¦ ¦ ¦ ¦
¦ береза ¦0.30 ¦ 0.32¦0.20 ¦ 11 ¦ - ¦0.11¦ 8.8¦0.18¦ 4.3 ¦
¦ осина ¦0.30 ¦ - ¦0.20 ¦ 9 ¦ - ¦0.11¦10 ¦0.17¦ 4.8 ¦
+---------+-----+-----+-----+-----+----+----+----+----+-------+
¦Опад хвои¦ ¦ ¦ ¦ ¦ ¦ ¦ ¦ ¦ ¦
¦(сосна ¦ ¦ ¦ ¦ ¦ ¦ ¦ ¦ ¦ ¦
¦обыкно- ¦ ¦ ¦ ¦ ¦ ¦ ¦ ¦ ¦ ¦
¦венная) ¦0.30 ¦ 0.30¦0.10 ¦ 30 ¦0.33¦0.09¦ 7.2¦0.13¦ 3.5 ¦
+---------+-----+-----+-----+-----+----+----+----+----+-------+
¦ верховые пожары ¦
+---------T-----T-----T-----T-----T----T----T----T----T-------+
¦Хвоя ¦ ¦ ¦ ¦ ¦ ¦ ¦ ¦ ¦ ¦
¦(листья) ¦ ¦ ¦ ¦ ¦ ¦ ¦ ¦ ¦ ¦
¦и тонкие ¦ ¦ ¦ ¦ ¦ ¦ ¦ ¦ ¦ ¦
¦веточки в¦ ¦ ¦ ¦ ¦ ¦ ¦ ¦ ¦ ¦
¦кронах ¦ ¦ ¦ ¦ ¦ ¦ ¦ ¦ ¦ ¦
¦деревьев ¦ ¦ ¦ ¦ ¦ ¦ ¦ ¦ ¦ ¦
¦(кустар- ¦ ¦ ¦ ¦ ¦ ¦ ¦ ¦ ¦ ¦
¦ника) ¦ - ¦ - ¦ - ¦ - ¦ - ¦ - ¦ - ¦ - ¦ - ¦
L---------+-----+-----+-----+-----+----+----+----+----+--------
Примечание. W(*) - критическое влагосодержание (при W > W(*)
*
горение прекращается); В(з) = В(1) (пи) примерно = В(з); ню' -
среднее значение пульсаций скорости, равное 1 м/с (по данным
киносъемки: Т(0) = 293 - 298 К).
Надо сказать, что, если известна выгоревшая площадь А, то,
используя (А.7), легко определить время горения:
------------------------------------
¦ 2 А
t(г) = v (---------------------------------). (А.16)
пи омега(С) (омега(А) + омега(В))
Приложение Б
(обязательное)
МАТЕМАТИЧЕСКАЯ МОДЕЛЬ И АЛГОРИТМ РАСЧЕТА
ТЕКУЩИХ ВЫБРОСОВ ПОЛЛЮТАНТОВ В АТМОСФЕРУ
ПРИ ВЕРХОВЫХ ЛЕСНЫХ ПОЖАРАХ
Для повальных верховых лесных пожарах справедлива та же общая
математическая модель (А.1) - (А1.3) для определения выбросов
загрязняющих веществ и теплоты, что и в Приложении А, если под
контуром S понимается контур верхового лесного пожара.
Повальный верховой пожар возникает, если выполняется следующее
условие:
-----------------
h(2) < h(n) = k v (m(з) омега (n),
u(~~) > 2.0 м/с, (Б.1)
где h(2) - высота нижней границы полога леса;
k = 1.6 - эмпирический коэффициент;
u(~~) - скорость ветра в кронах деревьев.
Скорость распространения верхового лесного пожара в том
случае, если направление распространения совпадает с направлением
ветра, определяется формулой:
----------------------------------------
_ ¦ _2 _
u(~~) [q + v (q + 8q эпсилон (пи + 1) (Тэта(0) - 1)
омега = ---------------------------------------------------,
_
2q (пи + 1)
4 ро(i) фи(i н) С(p i)
пи = Сумма --------------------. (Б.2)
i=1 ро(5~~) С(p5)
_ -1 q(з) М(1) альфа(с)
q = (ро(5~~) С(р 5) Т(~~) {[------------------ +
М(с)
+ q(5) (1 - альфа(с)) ню(т)] ро(1) фи(1н) -
- q(2) ро(2) фи(2н)},
альфа ламбда(т)
эпсилон = --------------------------------, (Б.3)
2
дельта h (ро(5~~) С(р5) u(*~~))
_
Т(0) k q(0) _
Тэта(0) = ----- = 1 + ------, q(0) = q¦
Т(~~) 1 + пи ¦фи(2н) = 0,
где альфа(с) - коксовое число;
М(1), М(с) - молекулярная и атомная массы кислорода и
углерода;
q(5) - теплота сгорания газообразных продуктов пиролиза;
ню(т) - массовая доля газообразных горючих продуктов пиролиза;
q(2) - теплота испарения связанной воды;
k = 1.6 - эмпирическая постоянная;
q(3) - теплота сгорания коксика;
ро(1), ро(2), ро(3) и ро(5) - плотности сухого органического
вещества, связанной с ним воды, коксика и воздуха;
альфа - коэффициент теплообмена полога леса с приземным слоем
атмосферы;
ламбда(т) - коэффициент турбулентной теплопроводности;
Т(0) - максимальная температура во фронте верхового лесного
пожара;
индекс (~~) приписывается величинам для нормальной среды в
дали от фронта;
н - начальным значениям параметров состояния в зоне лесного
пожара.
При применении формулы (Б.2) возможно, что подкоренное
выражение будет отрицательным. Данная ситуация отвечает большим
значениям влагосодержания ЛГМ и верховой пожар в этом случае не
возникает.
Формулы (Б.2), (Б.3) получены в результате аналитических и
численных исследований. При 2 м/с < u(*) < 4 м/с, ню(г) = 0.7,
альфа(с) = 0.06, 0% < W < 85% (W - влагосодержание ЛГМ),
0.15 кг/куб. м < ро фи < 0.3 кг/куб. м выражение (Б.3)
аппроксимирует результаты численных расчетов с погрешностью, не
превышающей 10%.
Для определения контура повального верхового лесного пожара
можно использовать ту же формулу (А.5), что и для низовых лесных
пожаров с полуосями a(в), b(в) и координатой x(0).
(В) (Н)
(омега(А) + омега(В) ) t (Н)
а(в) = --------------------------, b(В) = омега(С) t, (Б.4)
2
(В) (Н)
омега(А) - омега(В) ) t
x(0) = --------------------------.
2
Здесь ось x направлена в сторону ветра;
(В)
омега(А) - скорость верхового лесного пожара, определяемая по
формуле (Б.2);
(Н)
омега(С) - скорость низового лесного пожара,
распространяющегося перпендикулярно скорости ветра, и определяемая
по формуле (А.12);
(Н)
омега(В) - скорость низового лесного пожара, который
распространяется против скорости ветра ((1.13)).
Для повального верхового пожара при определении М(альфа) можно
использовать ту же методику, что и для низового лесного пожара.
Очевидно, что в этом случае сгорает не только напочвенный покров,
но и масса ЛГМ в кронах деревьев. Поэтому для определения выбросов
при верховых лесных пожарах необходимо использовать формулу:
(В) (Н)
К (альфа В) пи m(З в) (омега(А) + омега(В) )
М(альфа В) = --------------------- {1.5 [----------------------- +
2 N(*) 2
--------------------------------------
¦ (В) (Н)
(Н) ¦ (омега(А) + омега(В) ) (Н)
+ омега(С) ] - v (----------------------- омега (С) }
2
2 N
t Сумма омега(n) (альфа(j)). (Б.5)
j=1
В отличии от низового лесного пожара верховой не
распространяется против ветра. Поэтому в формуле (Б.5)
суммирование по j ограничивается N(*) < N, где N - полное число
разбиений для всего эллипса, моделирующего контур низового лесного
пожара. Известно, что верховой лесной пожар возникает и
распространяется лишь в том случае, когда равновесная скорость
ветра в пологе леса в одномерном случае удовлетворяет условию:
u(альфа) > u(альфа *) = 2 м/с. (Б.6)
В двумерном случае это условие принимает вид:
u(n альфа) (альфа(j)) > u(альфа *). (Б.7)
Из условия равенства:
u(n альфа) (альфа(j)) = u(альфа *) (Б.8)
легко находим предельный угол альфа(j *), при котором еще
имеет место распространение верхового лесного пожара. Этому углу и
соответствует значение N(*).
Общее количество поллютантов, выделяющихся при распространении
повального верхового лесного пожара равно:
М(альфа) = М(альфа В) + М(альфа Н), (Б.9)
где М(альфа Н) определяется по формуле (А.9), а М(альфа В) -
по (Б.5).
Легко видеть, что при повальном лесном пожаре М(альфа)
примерно = t, как и в случае низового лесного пожара. Площадь
повального верхового лесного пожара определяется по формуле:
пи (В) (Н) (Н) 2
А(В) = -- (омега(А) + омега(В) ) омега(С) t . (Б.10)
2
Если известна площадь верхового лесного пожара S(в), то время
горения легко определить по формуле:
-----------------------------------------
¦ 2 S(В)
t(г) = v (--------------------------------------). (Б.11)
(В) (Н) (Н)
пи (омега(А) + омега(В) ) омега(С)
Приложение В
(обязательное)
ФИЗИКО-МАТЕМАТИЧЕСКАЯ МОДЕЛЬ И МЕТОДИКА РАСЧЕТА
ТЕКУЩИХ ВЫБРОСОВ ВРЕДНЫХ ВЕЩЕСТВ В АТМОСФЕРУ
ПРИ ГОРЕНИИ ТОРФЯНИКОВ
К отдельному типу лесных пожаров относятся почвенные и
торфяные пожары, которые возникают под влиянием заглубления
горения ЛГМ при низовых лесных пожарах. Под влиянием нагрева от
фронта низового лесного пожара и притока окислителя из окружающей
среды высушивается, гидролизуется, а затем загорается подстилка -
слой полуразложившихся листвы, хвои, травы и тонких веточек.
Плотность подстилки в абсолютно сухом состоянии изменяется в
пределах от 30 до 300 кг/куб. м, что выше, чем плотность опада.
Влагосодержание подстилки тоже выше, чем у напочвенного покрова и
по данным /15/ изменяется, как правило, в пределах от 470% до
560%. Поэтому подстилка высыхает значительно медленнее, чем
напочвенный покров, и лишь в очень сухую погоду достигает
минимального значения - 14%. Если под подстилкой расположен слой
продуктов неполного разложения растительных материалов, называемых
торфом, то процесс заглубления горения может продолжаться и может
возникнуть торфяной пожар.
По условиям образования различают низинный, переходный и
верховой торфы. Низинный торф залегает на глубине 1,6 - 3,5 м и
имеет степень разложения 26 - 31%, зольность 6 - 12% и влажность в
естественном состоянии 86 - 92%, а верховой - глубину залегания
1,6 - 4,4 м, степень разложения 18 - 46%, зольность 2,7 - 4% и
влажность 89 - 93%. Под степенью разложения понимают отношение
массы гумуса (разложившейся части органики) ко всей массе торфа.
Низинный торф образуется на болотах в результате грунтового
питания болот водой, а верховой - при поступлении воды из
атмосферы. Влагосодержание торфов существенно больше
влагосодержания подстилки и составляет от 2500 до 3070%. С
увеличением степени разложения его влагосодержание уменьшается.
Как и у лесных горючих материалов, основными горючими элементами у
торфов являются углерод (52 - 56% от общей массы) и водород (5 -
6% от общей массы), кроме того, в составе торфа имеется от 30% до
40% атомов кислорода, связанного в молекулах химических веществ,
из которых состоит торф. Среднее значение величины теплотворной
способности торфа равно 5500 кКал/кг.
При нагревании торф высушивается, затем происходит его пиролиз
с образованием газообразных горючих компонентов и кокса, и если в
зоне нагрева имеется окислитель, то происходит воспламенение и
последующее горение торфа. Из наблюдений следует, что горение
торфа в естественных условиях после его заглубления под слой почвы
в условиях избыточного влагосодержания и недостатка кислорода
происходит в режиме тления. Линейная скорость распространения
фронта торфяного пожара составляет в среднем 7 мм/час. Горение
носит диффузионный характер, то есть лимитируется поступлением
окислителя.
Максимальная температура горения изменяется в пределах
623К є Т(г) є 673К. Фронт торфяного пожара неоднороден, то есть
горение носит очаговый характер. В результате горения образуются
продукты полного (CO(2), H(2) O и пепел), неполного окисления (CO)
и пиролиза торфа - метан (CH(4)), водород (H(2)), сажа, дым. Пепел
имеет белый цвет, рыхлую структуру и отваливается на вертикальных
участках. В процессе горения торфа первоначально плоская
поверхность фронта горения может принимать форму параболоида
вращения. В связи с тем, что горение при торфяных пожарах носит
подпочвенный характер, его очень трудно обнаружить. Над
поверхностью почвы при пожаре на торфянике стелется синеватый
дымок (угарный газ - газообразный продукт неполного горения торфа
в толстых слоях имеет синеватый цвет). В отличие от низового
лесного пожара, процессы переноса массы, энергии и количества
движения имеют существенно меньшую скорость, чем при низовых
пожарах. В частности, течение газа на торфяных пожарах принимает
характер просачивания газа через поры, называемое фильтрацией,
что, в конечном счете, существенно замедляет скорость поступления
окислителя к фронту торфяного пожара и предопределяет малую
скорость его распространения.
Таким образом, общая схема физико-химических процессов при
лесных пожарах, приведенная на рис. 3.3, сохраняет свою силу и для
пожаров на торфяниках, но, в отличие от низового и верхового
лесных пожаров, роль излучения при передаче тепловой энергии из
зоны горения торфяного пожара мала по сравнению с передачей
энергии путем теплопроводности и фильтрации продуктов горения.
Другим отличием является то, что фронт горения представляет собой
поверхность, положение которой в пространстве определяется
расположением торфяной залежи и приходом кислорода из приземного
слоя атмосферы. Поскольку горение при торфяных пожарах имеет
диффузионный характер, то фронт горения можно отождествлять с
поверхностью, расположенной внутри торфяной залежи, на которой
температура Т = Т(г), где Т(г) - температура горения, а плотность
диффузионного потока кислорода q(1) = q(1:), где q(1:) - значение
q, соответствующее Т = Т(г).
Пусть поверхность горения S(г), плотность торфа в торфяной
залежи ро(т) и скорость нормального распространения горения
омега(n) по торфянику известны.
Тогда для массовой скорости выброса альфа-поллютанта и тепла в
атмосферу имеем уравнение и начальные условия:
dМ(альфа)
--------- = интеграл интеграл К(альфа) К ро(т) омега(n) ds,
dt S(г)
альфа = 1,...,N, М(альфа)¦ = 0. (В.1)
¦t=0
dQ
-- = интеграл интеграл q К(альфа) К(омега(n)) ро(т) ds,
dt S(г)
Q¦ = Q(0), (В.2)
¦t=0
где К(альфа) - удельный выброс альфа-поллютанта в атмосферу;
К - коэффициент полноты сгорания;
N - общее количество поллютантов, возникающих при торфяном
пожаре.
Если считать, что поверхность горения - плоскость,
перпендикулярная вектору силы тяжести, а К(альфа), ро(т) и
омега(n) не зависят от координаты точки внутри торфяной залежи, то
уравнения (В.1), (В.2) упрощаются и принимают вид:
dМ(альфа)
-------- = К К(альфа) ро(т) омега(n) S(г), М(альфа)¦ = 0,
dt ¦t=0
(В.3)
dQ
-- = q K(альфа) К омега(n) ро(т) S(г), Q¦ = Q(0), (В.4)
dt ¦t=0
Значение полноты сгорания К в (В.3) и (В.4) определяется по
формуле (А.10), где величина предельного влагосодержания W(*) =
= 30.
Значения удельного выброса (коэффициента генерации) К(альфа)
даются в таблице Г.1 Приложения Г.
Если считать, что все величины в правых частях (В.3), (В.4) не
зависят от времени, то интегрируя (В.3), (В.4) по времени с учетом
начального условия, получаем значение выброса для любого
альфа-поллютанта и тепла:
М(альфа) = К К(альфа) ро(т) омега(n) S(г) t, (В.5)
Q = q К(альфа) К омега(n) ро(т) S(г) t. (В.6)
Выражение для скорости горения омега(n) имеет вид:
омега(n) = омега(0) (1 - W / W(*)). (В.7)
Здесь омега(0) = 0.002 мм/с, W и W(*) - влагосодержание и
предельное влагосодержание торфа.
Если известна масса торфяной залежи М(0), то, пользуясь
определением М(альфа), из (В.5), (В.6) получаем:
М = М(0) - ро(т) S(г) омега(n) t, (В.8)
где М - текущая масса торфа в торфяной залежи.
Очевидно, что в момент окончания горения М = 0 и из этого
условия находим время горения:
t(г) = М(0) / ро(г) S(г) омега(n). (В.9)
Подставляя (В.9) в (В.8), легко находим итоговый выброс
альфа-поллютанта при сгорании всей торфяной залежи:
М(альфа) (t(г)) = К К(альфа) М(0). (В.10)
Формулы (В.8) и (В.10) однозначно определяют текущий и
итоговый выброс альфа-поллютанта.
Приложение Г
(обязательное)
ПОНЯТИЕ О МАКСИМАЛЬНОЙ МОДЕЛИ ЛЕСА.
БАЗА ДАННЫХ ДЛЯ МОДЕЛЕЙ ВЫБРОСОВ ЗАГРЯЗНЯЮЩИХ
ВЕЩЕСТВ ПРИ ЛЕСНЫХ ПОЖАРАХ
Очевидно, что величина выбросов поллютантов в атмосферу
зависит от типа лесного пожара, запаса и типа лесного горючего
материала и его влагосодержания, типа лесных фитоценозов,
метеоусловий, рельефа местности и ряда других условий. Для оценки
сверху величины выбросов поллютантов в атмосферу при возгорании
лесных фитоценозов целесообразно иметь так называемую максимальную
модель леса.
Под максимальной моделью леса будем понимать гипотетический
лесной массив, в котором запас лесных горючих материалов для
каждого из ярусов леса максимален, а их влагосодержание
минимально. Очевидно, что при горении такого лесного массива
температура горения и интенсивность пожара (количество тепловой
энергии, выделившейся на единицу длины кромки пожара) будут
максимальны, чем и объясняется название - максимальная модель
леса. Как правило, при лесных пожарах сгорают тонкие (до 7 мм в
диаметре) веточки, хвоинки (в хвойных лесах) и листья (в
лиственных). На основании данных о запасе лесных горючих
материалов для различных ярусов леса и теплотворной способности, а
также вышеизложенного в качестве максимальной модели хвойного леса
предлагается следующая совокупность ярусов гипотетического леса:
1. Ярус мхов, лишайников с включениями из опавших хвоинок и
тонких веточек (нулевой слой). Высота его составляет 15 см,
плотность ро = 20 кг/куб. м, запас ЛГМ m(0)= 3,0 кг/кв. м,
теплотворная способность q = 19446 кДж/кг, влагосодержание W = 8%.
2. Первый ярус леса - травы и кустарники. Высота слоя - 2 м,
плотность ЛГМ в этом слое ро(с) = 0,8 кг/куб. м, запас ЛГМ в
первом ярусе m(0) = 1,6 кг/кв. м, теплотворная способность q =
= 17808 кДж/кг, влагосодержание W = 80%.
3. Второй ярус леса - подрост - представляет собой
совокупность деревьев высотой до 6 м. Высота слоя h(2) = 6 м,
толщина дельта(2) = 5 м, то есть нижняя граница полога подроста
(полог подроста - совокупность крон молодых деревьев) находится в
первом ярусе леса (первый и второй ярусы частично перекрываются).
Плотность слоя ро(с) = 0,4 кг/куб. м, запас ЛГМ m(o) =
= 2,0 кг/кв. м, теплотворная способность q = 21949 кДж/кг,
влагосодержание W = = 80%.
4. Третий ярус леса - совокупность крон деревьев. Высота
_
верхней границы полога леса h(3) = 22 м, высота нижней -
h(3) = 5 м, толщина слоя дельта = h(3) - h(3) = 17 м, то есть
второй и третий ярусы леса частично перекрываются. Плотность слоя
ЛГМ в пологе леса ро(с) = 0,3 кг/куб. м, запас ЛГМ m(з) =
= 5,1 кг/кв. м, теплотворная способность q = 21949 кДж/кг,
содержание воды W = 90%.
Используя лесотаксационные описания лесных массивов,
информацию о вертикально-фракционной структуре леса, данные о
запасе ЛГМ, можно построить максимальные модели леса для каждого
региона страны.
Кроме данных о запасе лесных горючих материалов, для расчета
М(альфа) необходимо знать величины К и К(альфа) для каждого типа
лесного пожара.
Для определения М(альфа) необходимо знать величины W(*) и К.
Эти данные приведены в таблице Г.1.
Таким образом, располагая данными о типе лесного фитоценоза,
типе лесного пожара, данными о запасе ЛГМ и его влагосодержании,
приведенными выше, можно оценить выбросы поллютантов в атмосферу в
результате горения ЛГМ при различных лесных пожарах.
Следует отметить, что в областных управлениях лесами и на
областных базах авиационной охраны лесов ведутся карточки учета
лесных пожаров, образец которой представлен ниже.
Располагая информацией из этой карточки учета, а в ней даны
тип пожара и размер выгоревшей площади, можно легко рассчитать
итоговые выбросы загрязняющих веществ и дать оценку величины
ущерба в результате попадания выбросов поллютантов в атмосферу.
Таблица Г.1
ЗНАЧЕНИЯ КОЭФФИЦИЕНТОВ ЭМИССИИ ДЛЯ РАЗЛИЧНЫХ
ТИПОВ ЛЕСНЫХ ПОЖАРОВ
--------------------------T--------------------------------------¬
¦Название поллютанта и его¦ К(альфа) для различных лесных ¦
¦ формула ¦ пожаров ¦
¦ +------------T------------T------------+
¦ ¦ низовой ¦ пожар на ¦ повальный ¦
¦ ¦ ¦ торфянике ¦ верховой ¦
+-------------------------+------------+------------+------------+
¦Оксид углерода CO ¦ 0.135 ¦ 0.135 ¦ 0.135 ¦
¦Диоксид углерода CO(2) ¦ 0.094 ¦ 0.094 ¦ 0.094 ¦
¦Оксиды азота NO(x) ¦ 0.000405 ¦ 0.000405 ¦ 0.000405 ¦
¦Сажа C ¦ 0.0062 ¦ 0.011 ¦ 0.0014 ¦
¦Дым (ультрадисперсные ¦ ¦ ¦ ¦
¦частицы SiO(2)) ¦ 0.0345 ¦ 0.055 ¦ 0.014 ¦
¦Метан CH(4) ¦ 0.075 ¦ 0.075 ¦ 0.075 ¦
¦Непредельные углеводороды¦ 0.011 ¦ 0.011 ¦ 0.011 ¦
L-------------------------+------------+------------+-------------
Рис. Г.1. Карточка учета лесного пожара
Приложение N 5
____________ авиабаза к приказу по Центральной авиабазе
____________ авиаотделение от "___" ________ 199__ г. N 59
КАРТОЧКА УЧЕТА ЛЕСНОГО ПОЖАРА N
Принадл. лесн. фонда _____________ Лесхоз ______________________ Леснич-во ____________________
Квартал (урочище)_________________ Тип леса _________________ га Район прим. сил ______________
Дата обнаружения _________________ Площ. обнар. ________________ Способ обнаруж. ______________
Причина пожара ___________________ Азимут и раст. ______________ Ширина _______________________
и долгота ________________________ Дата ликвидац. ______________ Площ.ликвидац. _______________
лесная _______________________ га, нелесная ________________ га, в т.ч. верх. ______________ га
подземный ____________________ га
-------T--------------------T--------T------T---------T-----------------------T------T--------¬
¦Число,¦ Площадь в га ¦Характе-¦ Вид ¦Состояние¦ Работает ¦ Треб.¦ Налет ¦
¦месяц +----T------T--------+ристика ¦пожара¦ или ¦ людей тех. средств¦допол-¦ ¦
¦ ¦лес-¦нелес-¦ в том ¦ ¦ ¦ причина +-----------T-----------+ ни- +---T----+
¦ ¦ная ¦ ная ¦ числе ¦ ¦ ¦неосмотра¦тип группы,¦ тип, ¦тельно¦тип¦час.¦
¦ ¦ ¦ ¦верховой¦ ¦ ¦ ¦количество ¦количество ¦ ¦ВС ¦мин.¦
L------+----+------+--------+--------+------+---------+-----------+-----------+------+---+-----
|
