Введение
Исследованиями прогноза возникновения и распространения лесных, степных и торфяных пожаров занимаются ученые во многих странах на протяжении уже более 50-ти лет. Наиболее актуальна проблема природных пожаров для Соединенных Штатов Америки, Канады, России, Франции, Испании, Греции, Австралии и других стран. Среди наиболее известных зарубежных систем оценки пожарной опасности можно отметить такие как канадская система CFFDRS (Canadian Forest Fire Danger Rating System) [1], национальная система США NFDRS (National Fire Danger Rating System) [2], австралийская Национальная система прогноза пожара NFBP (National Fire Behaviour Prediction system) [3]. В настоящее время в России прогнозирование пожарной опасности ограничивается оценкой пожарной опасности по условиям погоды в соответствии с показателем В.Г. Нестерова (GOST Р 22.1.09-99 [4]). В качестве общего недостатка разработанных на сегодняшний день и программно реализованных систем следует отметить их неуниверсальность, т.е. невозможность их использования в других климатических условиях без предварительной адаптации. Это связано в первую очередь с тем, что в эти системы заложены эмпирические данные, полученные для определенных лесных горючих материалов (ЛГМ) и климатических условий, которые резко отличаются в разных частях света. Основным недостатком является игнорирование реальной причины возникновения лесных пожаров, а именно, процесса сушки слоя ЛГМ на подстилающей поверхности, а также учет антропогенной нагрузки.
В связи с этим в работах [5-7] была предложена новая детерминированно-вероятностная методика прогноза лесной пожарной опасности, обоснованность которой была доказана в результате ретроспективной проверки. Для определения вероятностей появления природных факторов (молний) и антропогенной нагрузки использовались соответствующие частоты.
Целью моей работы является уточнение существующей модели с учетом формулы для оценки частоты появления источников загорания и индекса грозовой пожарной опасности.
Детерминированно-вероятностная модель прогноза катастроф
В настоящее время для прогноза катастроф широко применяется метод математического моделирования с использованием электронных вычислительных машин.
При математическом моделировании используются следующие типы математических моделей катастроф:
1) детерминированные;
2) вероятностные;
3) смешанные (детерминированно-вероятностные);
4) имитационные;
Согласно [5] наиболее эффективным инструментом познания катастроф являются детерминированные математические модели.
Детерминированной математической моделью физико-химического явления называется совокупность дифференциальных, интегральных, интегро-дифференциальных, трансцендентных и алгебраических уравнений, а также соответствующих граничных и начальных условий, которые адекватно описывают движение, деформацию и разрушение тел и поля физических величин (скорость, давление, плотность, температура, концентрация) для исследуемого катастрофического явления. Иными словами, каждая детерминированная математическая модель может рассматриваться как конкретная задача математической физики.
Наиболее плодотворно исследование различных проблем естествознания, в том числе и задач теории катастроф, с помощью смешанных математических моделей. Здесь необходимо отметить то, что если рассматриваются задачи о моделировании глобальных или региональных катастроф, например возникновения лесных пожаров, то вероятностный анализ этих проблем должен предшествовать использованию детерминированных математических моделей. Иными словами, результаты математического моделирования должны быть получены в режиме, опережающем время возникновения и развития катастроф.
Далее остановимся на одной из таких моделей, которая позволяет с определенной точностью осуществлять прогнозирование лесной пожарной опасности и обоснованность которой была доказана в результате ретроспективной проверки. Речь идет о детерминированно-вероятностной методике, предложенной на кафедре физической и вычислительной механики Томского государственного университета.
В работах [5-7] утверждается, что возникновение лесного пожара носит вероятностный характер и зависит не только от погодных условий и грозовой активности, но и от уровня антропогенной нагрузки, скорости ветра, влагосодержания растительных горючих материалов и реакционной способности этих материалов. В работе [5] предложена общая схема возникновения лесных пожаров
Рис.1. Схема возникновения лесного пожара.
Из схемы Рис.1. видно, что под действием осадков и при высушивании слоя растительных горючих материалов, достигается такое состояние РГМ, при котором влагосодержание слоя становится меньше критического значения. Именно в этом случае под воздействием сухих гроз, а также под воздействием антропогенной нагрузки происходит воспламенение растительных горючих материалов и возникает низовой лесной пожар.
Все эти факторы были учтены в предложенной детерминированно-вероятностной методике.
Далее возникает сложность в определении основного проводника горения. Считается, что для каждого индивидуального участка леса должен быть установлен свой основной проводник горения (ОПГ), но на практике имеет место огромное разнообразие образующихся смесей из растительных горючих материалов. Эта смесь может иметь признаки двух, трех и даже четырех ОПГ. Кроме того, при определении типов ОПГ необходимо учесть их сезонную динамику.
Здесь верхний индекс k как раз соответствует проводнику горения (ОПГ), в левой части неравенства находится текущее влагосодержание k-го ОПГ на i-м выделе в j-й временной интервал, в правой части – критическое влагосодержание k-го основного проводника горения на i-м выделе, при котором ОПГ не воспламеняется и лесной пожар не возникает.
Далее, после применения теории вероятностей и физических соображений, была получена следующая формула (динамическая модель) для оценки вероятности возникновения лесного пожара для j-го временного интервала лесопожарного сезона на i-й лесной территории [5-7]:
Здесь Pij – вероятность возникновения лесного пожара для j-го интервала (шага по времени) на i-й контролируемой лесной территории; Pi(А) – вероятность антропогенной нагрузки; Pi(ЛП/А) – вероятность возникновения пожара вследствие антропогенной нагрузки на площади Fi; Pi(М) – вероятность возникновения сухих гроз на площади Fi; Pi(ЛП/М) – вероятность возникновения лесного пожара от молнии при условии, что сухие грозы могут иметь место на площади Fi; Pij(C) – вероятность того, что влагосодержание слоя РГМ будет меньше критического (вероятность возникновения пожара по метеоусловиям).