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

I. N. Kartsan, E. S. Zhukova, S. V Litoshik

CHARACTERISTIC FEATURES OF TARGETING VIA AURBORNE LASER SPACE RADIO NAVIGATION DEVICES OF GLONASS SYSTEM

Characteristic features of usage of airborne laser systems for periodic adjustment of radio channels for range measurements of control stations of GLONASS satellite radio navigation system are considered in the article.

Keywords: airborne laser system, onboard control complex, ground control complex, radio frequency channel, the quantum-optical station.

© Карцан И. Н., Жукова Е. С., Литошик С. В., 2010

УДК630.43:662.519.6

В. С. Коморовский, Г. А. Доррер

МЕТОДИКА РАСЧЕТА ПАРАМЕТРОВ ЛЕСНЫХ ПОЖАРОВ КАК ДИНАМИЧЕСКИХ ПРОЦЕССОВ НА ПОВЕРХНОСТИ ЗЕМЛИ С ИСПОЛЬЗОВАНИЕМ ДАННЫХ КОСМИЧЕСКОГО МОНИТОРИНГА

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

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

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

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

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

Задача охраны лесов от пожаров является весьма актуальной во всем мире в настоящее время. Так, в Российской Федерации по данным информационной системы дистанционного мониторинга федерального агентства лесного хозяйства (ИСДМ-Рослесхоз) за пожароопасные сезоны 2008 и 2009 гг. произошло 35 337 пожаров, огнем

пройдено 15 565 510 га, в том числе 9 772 782 га покрытых лесом площадей [1]. Поскольку лес является одним из важнейших стратегических ресурсов, задача охраны лесов от пожаров, и в том числе моделирование и прогнозирование распространения пожаров, является чрезвычайно важной задачей. К настоящему времени разработано достаточное количество математических моделей распространения лесных пожаров и базирующихся на их основе методов прогнозирования. Эти модели и методы содержат разные подходы к рассмотрению процесса горения в природном пожаре, однако их общей особенностью является высокая требовательность к информационному обеспечению моделирования, т. е. условием построения адекватной модели и, соответственно, прогноза, является точное знание множества характеристик местности, на которой происходит пожар, состояния растительных горючих материалов, погодных условий. На текущий момент создать информационную базу для обеспечения построения прогнозов на основе таких моделей в России не представляется возможным.

Наиболее полная информация о лесопожарной обстановке на данный момент содержится в ИСДМ-Рослес-хоз. Эта информационная система существует и развивается с 1995 г. Источниками информации, представленной в ИСДМ-Рослесхоз, являются данные космического, авиационного и наземного мониторинга лесных пожаров. Причем, ввиду обширности зон наблюдения, методы космического мониторинга играют ведущую роль. В сборе данных участвуют спутники серий NOAA, TERRA, AQUA, SPOT, LANDSAT ETM+, МСУ-Э [1].

Наличие достоверного прогноза распространения и развития лесного пожара позволяет оценить угрозу природной среде, объектам экономики и населенным пунктам, принять необходимые меры по предотвращению ущерба, спланировать работу противопожарных сил. Одним из важнейших моментов при прогнозировании лесного пожара является построение его контуров. В настоящей работе представлена методика расчета контуров крупных (площадью более 200 га) лесных пожаров на основе ограниченной информации, доступной в ИСДМ-Рослесхоз [2; 3].

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

Данная методика может быть использована для моделирования других динамических процессов, упомянутых выше.

Исходные допущения. Были приняты следующие допущения:

1. Динамика изменения площади пожара определяется выражением

Б(Г) = к0(г - , (1)

где t - текущее время, сут; t0 - время возникновения пожара, сут; а - показатель динамики пожара; к0 - постоянный коэффициент имеющий размерность га/сута. Как ясно из геометрических соображений и будет показано ниже, с показателем а связано также изменение скорости фронта пожара: при а = 2 скорость фронта пожара постоянна, при а < 2 эта скорость со временем уменьшается, а при а > 2 - скорость возрастает.

2. Скорость фронта пожара в соответствии с принципом Гюйгенса представляется как у(ф, т) = У0(/)Х(ф), где у0(() -подлежащая определению максимальная скорость распространения фронта пожара (например, в направлении ветра); Х(ф) (1X1 < 1) - индикатриса полной скорости фронта, определяющая конфигурацию пожара; ф - направление распространения (0 < ф < 2 л). Таким образом, в общем случае

п0 = Пс(V, t), Х = Х( ^, Ф),

где V - величина скорости ветра, м/с. Там, где это не требуется, зависимость от V опускаем.

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

Формулы для расчета скорости фронта пожара. Для упрощения формул в данном разделе положим t0 = 0. Рассмотрим элементарное приращение площади й®. На рисунке показана кромка пожара в два близких момента времени t и t + Л. Из рисунка следует, что йБ = ппйШ, где пп - нормальная скорость фронта; й1 - приращение длины контура, которое вычисляется по формуле [4]

Й1 = 1йф = й(ф/х2(Ф) + Х2(Ф)|у0(т)йт . (2)

В свою очередь, нормальная скорость связана с полной скоростью соотношением

у(ф, t) = у0У )Х 2(ф)

Vп (ф, ) =

1, (уфt))2 2(ф)+Х'2(ф) '

V ^(ф, tу

Отсюда

I

йБ = йфdtv0(t)X2(ф)|у0(т)йт .

Приращение площади по всем направлениям распространения

2л л t

Й52л = | йБйф = 2dtv0(t)|йф|2(ф)|у0(т)йт.

0 0 0

Скорость прироста всей площади пожара

7^г Л t

— = 2^)(0|х2(ф)йф|Vo(т)йт , (3)

й 0 0

а динамика общей площади примет вид

Л t л

5(I) = 2|Х2(ф)йф|^(т)|^©йХйт .

0 0 0

Рассмотрим теперь закон изменения скорости фронта пожара. Будем искать функцию у0(^ в виде п0®(/р^, где величины р и п0® подлежат определению. Тогда интеграл по времени в уравнении (3) будет равен

|у0( т)й т = р У0/+\ 0 р

а скорость прироста площади

1 г 2|М |х2(ф)й ф.

С другой стороны, из первоначального предположения Б(() = к^ следует

— = к0аа-1

Л 0 ,

откуда получаем уравнение

2 12Р+' | Х2(ф)й ф= к0 аtа-1. (4)

Приравняв в уравнении (4) показатели степени при t, получим 2р + 1 = а - 1, откуда

Р = “-1.

2

Далее, приравняв множители при t, определим величину V

( ак0(Р+1) \1/2 _ а\к0 _ аД/к0

V0S _ ( р ) _ р _ „ г,,.. ч , (5)

где

2|Х2(Ф)^ ф 2 |х2(Ф)й ф

2Р( щ)

Формула (5) есть искомая формула для расчета скорости пожара.

На основе формулы (2) можно вычислить также периметр горящей кромки пожара. В частном случае, когда пожар развивается из очага малого радиуса (точечного очага)

р і Р _________________

Ь(і) _ 2|1йф _ 2|.0(т)йт|-у/х2 (щ Ф) + X'2(V,

0 0 0

Формула для скорости прироста периметра пожара

примет вид

йЩ)

йі

_ 2.0 (і)0(щ),

(6)

где

Є (V) _^Х2( V, Ф)+Х'2(щ,

Индикатрисы скорости фронта. Рассмотрим конкретные выражения для индикатрис.

Экспоненциальная индикатриса равна [4]

Х(ф) = ехр(а(^)(со8(/)-1)), где а(^) - коэффициент, зависящий от скорости ветра: а(^) = 0,785^ - 0,06^2.

Формула справедлива при условиях 0 < V < 3 м/с. Тогда интегралы, стоящие в выражениях (5) и (6), примут вид

Л

Р(м?) = |х2(^, ф)й ф =

0

л

= ехр(-2а( V)) | ехр(2а( V) оо8(ф))й ф,

Є( щ) _ ]уіХ 2( щ, Ф)+Х'2(щ

_ | XV ФЫ1 + 8Іп2 (2Ф)й Ф.

0 V 4

Вычисленные величины этих интегралов для некоторых значений скорости ветра приведены в табл. 1.

Эллиптическая индикатриса равна

Х(Ф) _ 1 - е<»> ,

1 - е(щ)СОЄ(Ф)

где е(щ) - эксцентриситет эллипса, зависящий от скорости ветра.

Для оценки зависимости эксцентриситета эллиптической индикатрисы от скорости ветра были аппроксимированы эллипсами графики индикатрис, приведенные в работе Ф. Альбини [5], для большого диапазона скоростей ветра и вычислены соответствующие эксцентриситеты. Была получена следующая аппроксимация зависимости эксцентриситета от скорости ветра:

е(щ) _ 1 - ехр(-0,4щ).

Затем были вычислены значения интеграла Р(щ) и Є(щ). Результаты расчетов приведены в табл. 2.

Сравнивая данные табл. 1 и 2, можно увидеть, что значения функций Р(щ) и Є(щ) для обеих индикатрис заметно различаются при одних и тех же скоростях ветра. При этом следует учитывать, что в первом случае скорость ветра вычислялась на высоте двух метров от земли, а во втором

- на высоте десяти метров, и этим также объясняется более вытянутая форма эллиптической индикатрисы.

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

Пусть погрешность определения площади пожара составляет 5? га. Рассмотрим вызванные этим погрешности параметров модели (1). Погрешность в оценке ко -эффициента к

д?

5? _—5к _ і “5к, ді

откуда

5к _^. іа

(7)

Погрешность в оценке времени распространения пожара:

ж

55 = —8t = kat a-l5t, дt

откуда

0

0

Таблица 1

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

V, м/с 0 0,1 0,2 0,5 1,0 2,0 3,0

а(щ) 0 0,078 0,155 0,378 0,725 1,33 1,815

Р(щ) 3,142 2,705 2,361 1,1695 1,178 0,818 0,686

6И 3,142 2,911 2,708 2,234 1,758 1,455 1,535

Таблица 2

Значения интегралов квадрата эллиптической индикатрисы в зависимости от скорости ветра

V, м/с 0 1,25 2,5 5 10 15

е(щ) 0 0,393 0,632 0,865 0,982 0,998

Р(у>) 3,142 1,469 0,782 0,588 0,539 0,193

6И 3,142 2,151 1,60 1,424 1,379 1,194

St = -

SS

k ata-1

Погрешность в оценке показателя a:

дS

SS = — Sa = kta in t Sa, дa

(8)

откуда

Sa = -

SS

Sa = -

= 0,07.

V0S = '

= 426 м/сут.

Р(ю) \ 1,178 С учетом возможной погрешности 5к величина скорости фронта будет находиться в следующем диапазоне: 275 < п05 < 577 м/сут.

По формуле (6) можно оценить скорость прироста периметра пожара: йЦі)

dt

- = 2v0 (t )Q(w) = 2 • 426 -1,758 = 1498 + 531 м/сут.

ке 1п Г (9)

Из формул (7)-(9) видно, что погрешности всех переменных, от которых зависит модель (1), зависят от времени и уменьшаются с ростом t.

Численный пример расчета скорости распространения пожара. Зарегистрированный в системе ИСДМ-Рос-лесхоз пожар К-1002 действовал в Идринском лесничестве Красноярского края с 10.05.2008 г. по 16.05.2008 г. Пожар был зарегистрирован на площади 393 га, площадь ликвидации 1163 га. Температура воздуха была 22,7 град., скорость ветра - 1 м/с, направление ветра - 180 град.

Примем, что погрешность определения площади составляет 55 = 100 га, Дt = 16 - 10 = 6 сут. Принимаем показатель скорости роста площади пожара а = 2, т. е. предполагаем, что скорость фронта пожара постоянна.

Вычислим коэффициент к и его погрешность в модели: К = ДУД^ = 770/36 = 21,4 га/сут2,

5к = 100/36 = 2,7 га/сут2.

Оценка погрешности коэффициента а по формуле (9) дает величину

100

21,4• 62 • 1п6

Поскольку скорость ветра невелика, воспользуемся экспоненциальной индикатрисой. Для скорости ветра

V = 1 м/с по табл. 1 определяем множители Р(V = 1,178 и б(^) = 1,758.

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

4к /21,4-104

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

Представленная в статье простая методика оценки параметров процесса распространения, основанная на использовании геометрической модели, позволяет осуществлять прогнозирование распространения лесных пожаров, исходя из информации, хранящейся в ИСДМ-Рослесхоз. В методике используются индикатрисы двух видов, в зависимости от скорости ветра в районе пожара: экспоненциальная и эллиптическая. Полученные формулы применены для численного расчета параметров реального пожара, оценены их погрешности. Описываемый метод расчета применяется для построения контуров пожаров на основании нейро-сетевых прогнозов в программном комплексе «Тайга-2» и в разработке информационной системы прогнозирования лесных пожаров для ИСДМ-Рослесхоз.

Библиографические ссылки

1. Информационная система дистанционного мониторинга федерального агентства лесного хозяйства [ Электронный ресурс]. Пушкино : ФГУ Авиалесоохрана, 2010 // URL: http://www.pushkino.aviales.ru/rus/main.sht.

2. Доррер Г. А. Проблемы прогнозирования крупных лесных пожаров // Материалы XI Всерос. науч.-практ. конф. : сб. ст. - Красноярск, 2009. - С. 5-8.

3. Коморовский В. С. Оценка возможности прогнозирования распространения лесных пожаров по данным «ИСДМ-Рослесхоз» [Электронный ресурс] // Современные проблемы информатизации в системах моделирования, программирования и телекоммуникациях : сб. ст. I Между-нар. науч. конф. // URL: http://econf.rae.ru/article/4679.

4. Доррер Г. А. Динамика лесных пожаров. Новосибирск : Изд-во СО РАН, 2008.

5. Albini F. A. Estimating wildfire behavior end effects // USDA Forest Service, Gen, Tec. Rep. INT-30, Ogden, 1976. (Intermountain Forest and Range Exp. Stn.).

V S. Komorovski, G. A. Dorrer

METHOD OF CALCULATION OF PARAMETERS OF FOREST FIRES, AS DYNAMIC PROCESSES IN THE EARTH SURFACE BASED ON SPACE MONITORING DATA

A simple method for parameters of forest fire calculation is suggested in this article. The method is based on the representation of a fire in terms of a mobile set on the surface of the Earth. The information base for the proposed methodology is the data of space monitoring offorest fires.

Keywords: forest fires, forecasting, space monitoring, dynamic processes.

© Коморовский В. С., Доррер Г. А., 2010