Агентный подход при моделировании лесопожарных ситуаций

УДК 630.43 Дата подачи статьи: 14.04.16

DOI: 10.15827/0236-235X.115.101-108

АГЕНТНЫЙ ПОДХОД ПРИ МОДЕЛИРОВАНИИ ЛЕСОПОЖАРНЫХ СИТУАЦИЙ

С.В. Яровой, аспирант, ach_bask@mail.ru

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

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

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

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

Природный пожар - неконтролируемый процесс горения, стихийно возникающий и распространяющийся в природной среде. Данный термин объединяет лесные, торфяные и степные пожары. Природные пожары наносят огромный экономический и экологический ущерб регионам, в которых они возникают. Для управления эффективной борьбой с природными пожарами разработано большое количество информационных систем, позволяющих с различной степенью точности моделировать их поведение (системы BehavePlus, FARSITE, WFDS [1-3] и др.). Данные системы основаны на математических моделях, имеющих различную точность, время расчета и требовательность к вычислительным ресурсам.

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

Агентное моделирование (agent-based modeling, ABM) представляет собой новый подход к имита-

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

Для моделирования лесопожарных ситуаций особый интерес также представляют геоинформационные системы (ГИС), обладающие возможностями хранения и визуализации пространственной информации. Совместное использование агентного моделирования и ГИС позволяет объединить преимущества обоих направлений, применить перспективные средства для моделирования процессов, протекающих в условиях реального ландшафта.

Идея моделирования кромки природного пожара с помощью агентного подхода не нова [5, 6]. Агентные модели такого рода создавались с помощью специализированных систем для программного описания деятельности агентов и состояния среды (агентных платформ), таких как NetLogo, StarLogo, Repast Simphony и др. Данный подход позволяет упростить программную реализацию

имитационной модели, но значительно затрудняет ее интеграцию с ГИС. Кроме того, в работах [5, 6] моделирование свободно распространяющегося пожара производится без моделирования противопожарных сил и средств.

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

Описание среды моделирования

Среда моделирования представляет собой участок карты реальной местности с нанесенным на нее слоем, описывающим растительные горючие материалы (РГМ). Таким образом, имеется возможность совмещать данные ГИС и агентной модели. В модели используется информация о населенных пунктах, реках, озерах, дорогах, оврагах, лесных насаждениях и др. Позиционирование на карте производится заданием пространственных координат, различных для каждого объекта. В качестве модели горючего используется модель тонкого слоя [7]. В данной модели лесной покров рассматривался как однородный слой на поверхности земли с усредненными характеристиками РГМ, определяемыми на основании пирологического описания выделов карты лесов.

Таким образом, для каждого местоположения на карте заранее задается тип поверхности, определяющий наличие основных проводников горения (ОПГ) и их тип. В зависимости от типа ОПГ определяются базовая скорость распространения кромки пожара УОПГ (м/мин.) и теплота воспламенения слоя Q (Дж/м2).

В качестве основных показателей среды также принимаются следующие параметры:

- величина скорости ветра по данным ближайшей метеостанции Уш;

- вектор направления ветра Ш;

- класс пожарной опасности.

Данные параметры в значительной степени влияют на скорость и направление распространения пожара, а также на интенсивность горения.

В модели выделены два типа агентов: участок фронта распространения пожара (агенты типа А) и противопожарная команда (агенты типа В).

Агенты типа А используются для моделирования поведения лесного пожара и в совокупности представляют контур пожара на карте. Модель допускает наличие одновременно нескольких таких контуров, распространяющихся независимо друг от друга в различных частях карты.

Агенты типа B имитируют действия противопожарных сил и оказывают воздействие на агентов типа A.

Свойства и характеристики агентов типа A

Агент, представляющий участок фронта распространения пожара, характеризуется следующими параметрами:

- текущие пространственные координаты (X, Y);

- вектор внешней нормали к контуру n = [nx, ny]T;

- вектор скорости движения и = [Ux, Uy]T;

- значение интенсивности горения I;

- значение уровня задымления M;

- ссылки L и R на соседних агентов для создания единого контура;

- текущее состояние горения onFire=[true, false]. Основная цель, которой агенты A руководствуются при исполнении, - перемещение по среде моделирования. Расчет пространственных координат данных агентов на каждом шаге моделирования осуществляется по заранее заданному закону. Также во время исполнения модели агенты A могут порождать новых агентов или погибать при определенных условиях. Область видимости агентов A ограничена расстоянием до соседних по контуру агентов за счет ссылок L и R.

Расчет пространственных координат. Расчет новых координат AX, AY на каждой итерации моделирования выполняется посредством основанного на методе подвижных сеток (МПС) алгоритма [8, 9]. Согласно данному методу, расчетная сеточная область не задается заранее, а продвигается и развивается вместе с контуром пожара, представленным в виде набора точек: Ci(t)={xi(t), yi(t), pi(t), L(i), R(i)}, ieN(t), где x,(f), y(t) - координаты точки на плоскости; p,(t) - вектор внешней нормали к контуру; L(i) и R(i) - номера точек контура, расположенных, соответственно, слева и справа от Ci; N(t) - множество номеров агентов, образующих контур в момент t.

Для расчета следующего состояния контура используются скорость распространения пожара Uo и индикатриса нормальной скорости ^(ф,ю).

Приращение координат в точке Ci(t) производится по формулам:

xi(t+At) = xi(t)+Uo^n^i,rayAtxos(ai), yi(t+At) = yi(t) +ио^п(ф>,ю) •At- sin(ai), где ф,- - угол между нормалью к контуру в точке Ci(t) и направлением ветра ю; ai - угол между нормалью к контуру в точке Ci(t) и горизонтальной осью.

Для использования данного алгоритма сопоставили каждую точку Ci(t) с агентами типа A в нашей модели. В качестве входных данных для алгоритма используются характеристики агентов (X,, Y, n, и, L, R), а также параметры среды, такие как скорость и направление ветра, класс пожарной опасности, тип

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

На рисунке 1 визуально представлен контур пожара при расчете, а также схематично указано положение агента Аг на разных итерациях процесса моделирования t и t+1. Стрелка п- это направление нормали к контуру для г агента, г е N(0 Вектор Ш - направление ветра.

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

Расчет интенсивности горения РГМ. Согласно [10], интенсивность I рассчитывается по формуле I = Q•Vi•pbe, где и - значение скорости распространения пламени, м/мин.; рЬе - эффективная плотность горючего, кг/м3, то есть количество горючего материала (из общего запаса, содержащегося в слое), которое участвует в воспламенении; Q - теплота воспламенения (количество тепла, необходимое для доведения единицы веса горючего материала до воспламенения), Дж/кг.

Поскольку в предлагаемой модели используется модель тонкого слоя горючего, была использована эмпирическая формула, в которой в качестве Q используется теплота сгорания слоя, МДж/м2 [11].

Таким образом, интенсивность горения агента Аг, на которого не оказывается воздействие агентом В, при моделировании рассчитывается следующим образом: I = фи + Тад + 1ад)/3, где I - интенсивность горения агентаАг; Q- поверхностная теплота сгорания слоя; и - значение скорости движения

агента Ai, 1щ) - интенсивность горения агента Аад; Il(í) - интенсивность горения агента Аод.

То есть при расчете интенсивности горения текущего агента Ai учитывается аналогичный показатель соседних агентов.

Также на каждом шаге моделирования показатель интенсивности I может быть изменен случайным образом в положительную или отрицательную сторону с вероятностью Pnf. Это необходимо для имитации неоднородности таксационных выделов, ошибок в определении типов леса и других факторов, которые могут служить причиной ошибок прогноза [11].

Расчет задымленности. Задымление является одним из отрицательных факторов, увеличивающих ущерб от лесного пожара, а также значительно затрудняющих борьбу с ним. В настоящее время существует большое количество работ, посвященных моделированию дымового шлейфа [12, 13], однако реализация этих моделей требует высоких вычислительных затрат, а сами модели избыточны для данного исследования. Таким образом, была выбрана концепция локального пламени [8]. Согласно данной концепции, несколько огрубляя физическую картину распространения дымового облака, принимается, что каждое локальное пламя (в нашем случае каждый агент типа А) порождает локальный дымовой поток, конфигурацию которого в приземном слое можно определить следующим образом: ls(x, y, t, ф<) = иДх, y, <a)(t - ф< - ф), где (x, y) - координаты точки в области активного горения; фш - направление ветра, рад; ф - текущий угол, рад, 0 < ф < 2п; t - текущее время, сек., t0 < t < t0+r; t0 - время воспламенения горючего в точке (x, y), сек.; r - продолжительность горения горючего в точке (x, y); %s(w, фш - ф) - индикатриса локального дымового потока, определяющая его форму; uo - горизонтальная скорость распространения дымового облака в приземном слое.

Согласно [13], us0 вычисляется с помощью выражения Uso = 0,28 + 0,607w - 0,0217w2, где w - скорость ветра над пологом леса, м/с.

В предлагаемой модели на каждом шаге моделирования агенты типа A перемещаются в новую точку пространства, поэтому r всегда будет равна времени, эквивалентному одному шагу моделирования (при реализации модели использовался шаг = 5 мин.).

Сам уровень задымления M непосредственно около агента A принят в виде безразмерной величины, зависящей от интенсивности горения I, и может принимать значение от 0 (при I = 0) до 1 (при I ^ max).

Добавление новых агентов в процессе моделирования. В процессе моделирования агенты А перемещаются, увеличивая размер контура пожара в десятки и сотни раз. Расстояние между агентами также растет, что приводит к очень грубому представлению контура пожара. Таким образом, суще-

/°о р / / Ai ¿ w b(,)(t) . .A(')(t+1^ n-

AnQ

О \ \ \ °.....о /AR(i)(t) V Л

Рис. 1. Представление контура пожара

Fig. 1. Fire edge presentation

ствует необходимость упорядочения расстояния между соседними агентами Ai(t) и AR(i)(t), i е N(t). Для этого рассчитывается расстояние l(i) между каждой парой агентов. Если l(i)>lmax, порождается новый агент Aj(t), который располагается между парой Ai(t) и AR(i)(t). В качестве характеристик нового агента Aj(t) записываются средние значения соответствующих характеристик агентов Ai(t) и AR(i)(t). Данный процесс схематично представлен на рисунке 2.

M)

Рис. 2. Добавление нового агента типа A Fig. 2. Adding a new A agent

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

Удаление агентов в процессе моделирования. Удаление (гибель) агентов типа А производится в двух случаях.

Во-первых, удаление агентов происходит, если расчетное расстояние 1(г) между агентами Аг(1) и Ая({() меньше заданного значения 1тт. В этом случае агент Аг(1) удаляется, а ссылки Я и Ь у агентов Ая({) и Ащ) корректируются следующим образом: Я(Ь(1))=Я(1), Ь(Я(1))=Ь(1).

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

Во время данного анализа вычисляются точки пересечения каждого из отрезков контура, образуемых агентами Аг и Ащ), со всеми остальными отрезками Аг и Ая(г).

Пересечение контуров, в том числе и самопересечение, имеет место, если для некоторой пары отрезков Аг, Ая(г) и А], Аяу) (г фу, г ф Я(), у Ф Я(г)) точка их пересечения £ лежит внутри каждого из них. В этом случае изменяются ссылки у агентов Аг и Аяуу. Я(г)=Я(!'), Ь(Я('))=1. Агенты, оказавшиеся

внутри области получившегося контура, удаляются.

Переход состояний агента. Как уже было указано, для агентов А возможны два состояния, характеризующие текущее горение. Обозначим эти состояния как опР (агент «горит») и по1Р (агент «не горит»).

Переход агента из состояния опР в по£Р возможен в нескольких случаях.

Во-первых, агент может самостоятельно перейти в по£Р, когда попадает в область, в которой отсутствуют ОПГ. Во-вторых, он может быть «потушен» агентами типа В, оказывающими на них воздействие, снижающее интенсивность горения I.

Переход агента из состояния по1Р в опР в данной модели не предусмотрен: если агент однажды перешел в состояние по£Р, он продолжает в нем находиться до конца процесса моделирования. Таким образом, можно сформулировать два простых правила перехода состояния для агентов типа А.

Правило агента А № 1. Если УОПГ = 0, переходи в состояние по£Р («не горит») и прекращай движение по среде, иначе - продолжай движение.

Правило агента А № 2. Если I = 0, переходи в состояние по1Р («не горит») и прекращай движение по среде, иначе - продолжай движение.

Исходя из данных правил, получим таблицу переходов состояний агентов типа А (см. таблицу).

Переходы состояний агентов типа А

State transitions of A agents

Состояние Копг I Переход состояния

onF 0 0 onF ^ notF

onF >0 0 onF ^ notF

onF 0 >0 onF ^ notF

onF >0 >0 onF ^ onF

notF любое любое notF ^ notF

Свойства и характеристики агентов типа B

Агент, моделирующий противопожарную команду, характеризуется следующими параметрами:

- текущие пространственные координаты (X, Y);

- значение скорости движения по местности и;

- направление обхода контура пожара а = [R, L];

- значение скорости тушения кромки пожара VT (мера воздействия агента B на интенсивность горения I агента A);

- номер fc ближайшего по расстоянию агента типа A ;

- значение психологического состояния P;

- текущее состояние борьбы с пожаром inFight [true, false].

При исполнении агенты типа B руководствуются только одной целью - «потушить» все очаги возгорания. Для этого агент B перемещается по среде к ближайшему агенту Ai и воздействует на его показатель интенсивности горения I. Как только данный показатель становится равным 0, агент Ai меняет свое состояние на «не горит», а агент B перемещается к следующему агенту AR(i).

Перед запуском модели выбирается количество агентов типа B, затем задаются характеристики каждого агента (и, а, Vt, P). Начальное расположение агентов на карте также задается заранее в зависимости от выбранной тактики тушения низового лесного пожара (схемы организации тушения). Схем тушения может быть четыре [14, 15]:

- по всему периметру (для слабых по интенсивности пожаров);

- охватом с фронта (для слабого и среднего по интенсивности лесного пожара);

- охватом с фланга (для пожаров средней интенсивности);

- сведением на клин (для сильного лесного пожара).

Область видимости агентов B ограничена только самой средой моделирования. Где бы ни находился данный агент на карте, он всегда будет передвигаться в сторону ближайшего очага пожара (агента A).

В процессе моделирования расчет новых координат BX, BY производится различными способами в зависимости от текущего состояния борьбы с пожаром. Если inFight = true, то по мере «тушения» в качестве координат X, Y агента B принимаются текущие координаты агентов A. Если inFight = false, производится расчет расстояния Sji от текущего агента Bj до каждого агента Ai. Затем выбирается наименьшее из этих расстояний и производится перезапись Bfc=i. Высчитывается нормаль направления движения N от Bj до A,-. Сами координаты BX, BY рассчитываются следующим образом:

Bx(t+1)= Bx(t)+Nx-v,

By(S+1)= Bx(t)+NY и

Добавление новых агентов в процессе моделирования. Автоматическое добавление новых агентов типа B в процессе моделирования не предусмотрено. Однако возможно добавление новых агентов вручную, если этого требуют цели и задачи исследования.

Удаление агентов в процессе моделирования. Удаление агентов типа B производится только в случае «гибели» противопожарной команды, которая может произойти в процессе «тушения» агентов A.

Взаимодействие агентов типа A и B. Рассмотрим более подробно взаимодействие агентов A и B, имитирующее процесс тушения участка кромки пожара противопожарной командой. Если агент B справится со своей задачей, то агент A перейдет в состояние notF («не горит»). В противном

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

Скорость тушения кромки пожара VT определяется перед моделированием с учетом профессиональной подготовки, физического состояния, оснащенности и количества пожарных в команде. При тушении пожара эффективность труда во многом зависит также от психологического состояния пожарных P, которое обусловлено большим числом факторов, таких как напряженность работ, задым-ленность, высокая температура и др. [15]. Данный показатель может иметь значение от 0 до 1.

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

Помимо рассмотренных факторов, на работоспособность влияет и большое количество факторов вероятностного характера. Обозначим совокупность этих факторов как F. Они могут как снизить работоспособность (характер растительности, состояние погоды, неуверенные действия руководителя), так и повысить ее (кратковременный отдых, смена режима работы, уверенные действия руководителя). Значение F рассчитывается случайно в заданных пределах и напрямую увеличивает или уменьшает показатель P.

Таким образом, интенсивность горения агента Ai, находящегося под воздействием агента Bj, рассчитывается по формуле Ii(t+1) = Ii(t)-(Vy  Pj), где Ii(t+\) - интенсивность горения агента Ai в период времени t+1; I,(t) - интенсивность горения агента Ai в период времени t; VTj - скорость тушения кромки пожара агентом Bj; Pj - психологическое состояние пожарных агента Bj.

Переход состояний агента. Для агентов B возможны два состояния, характеризующие текущее состояние борьбы с пожаром. Обозначим эти состояния как fight (агент «тушит кромку») и move (агент «перемещается по карте»).

Переход агента из состояния move в fight возникает, когда пространственное расстояние от агента B до агента Afc равно нулю, то есть противопожарная команда достигла кромки пожара.

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

Таким образом, сформулируем два простых правила перехода состояния для агентов типа B.

Правило агента B № 1: если твое состояние move и расстояние от агента B до агента Afc равно 0, переходи в состояние fight и начинай взаимодействие с агентом Afc, иначе - продолжай движение.

Правило агента B №2: если твое состояние fight и количество непотушенных агентов A в контуре равно 0, то переходи в состояние move и двигайся к ближайшему горящему агенту A, иначе -продолжай борьбу.

Реализация агентной модели

Представим два примера работы созданной системы.

Пример № 1.

Параметры среды, выставленные перед началом моделирования:

- тип ОПГ - лишайниковый;

- класс пожарной опасности - 2;

- скорость ветра под пологом леса - 1 м/с (Ю/З).

Согласно [14, 15], интенсивность представленного пожара определяется как слабая, поэтому выбрали схему тушения по всему периметру. Расположили 3 противопожарные команды по 6 человек вокруг пожара и задали направление обхода против часовой стрелки. Площадь очага пожара при обнаружении ~ 4 га. На рисунке 3 представлены исходная позиция пожарных команд и исходный контур пожара, а на рисунке 4 - положение пожарных команд и контур пожара спустя 4 часа (48 итераций моделирования).

Как видно из представленного примера, после 4 часов локализации пожар продолжает развиваться и потушено менее 50 % текущей кромки пожара. Одновременно с этим общая длина горящей кромки пожара на 4-м часе тушения приблизительно равна длине исходной кромки пожара, что говорит о том, что 3 отряда пожарных по 6 человек не способны справиться с данным пожаром.

Рис. 3. Пример работы имитационной системы № 1

Fig. 3. Example of № 1 simulation system operation

Рис. 4. Пример работы имитационной системы № 1

Fig. 4. Example of №1 simulation system operation

Пример № 2.

Параметры среды, выставленные перед началом моделирования:

- тип ОПГ - лишайниковый;

- класс пожарной опасности - 3;

- скорость ветра - 1,5 м/с (С/В).

Согласно [14, 15], интенсивность представленного пожара определяется как средняя, поэтому выбрали схему тушения охватом с фронта. Расположили 2 противопожарные команды по 12 человек около фронта в одной точке и задали направление обхода в противоположных направлениях. Площадь очага пожара при обнаружении ~ 3 га. На рисунке 5 представлены исходная позиция пожарных команд и исходный контур пожара, на рисунке 6 - положение пожарных команд и контур пожара спустя 3 часа (36 итераций моделирования). Как видим, спустя 3 часа локализации пожара фронт полностью потушен, а пожарные команды продолжают успешно продвигаться по флангам в сторону тыла пожара. Поскольку тушение фронта

I

Рис. 5. Пример работы имитационной системы № 2

Fig. 5. Example of № 2 simulation system operation

Рис. 6. Пример работы имитационной системы № 2

Fig. 6. Example of № 2 simulation system operation

пожара является наиболее опасным и трудоемким процессом относительно тушения других элементов пожара, то уже на этом этапе можно сказать, что для борьбы с данным пожаром была выбрана правильная схема тушения, а 2 отряда пожарных по 12 человек способны данную схему реализовать (см. http://www.swsys.ru/uploaded/image/2016_3/ 2016-3-аор/5оря).

В результате проведенного исследования автором была разработана мультиагентная модель для имитации прямого метода тушения низового лесного пожара немеханизированными противопожарными отрядами. На основе данной модели со-

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

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

Литература

1. Andrews P.L., Bevins C.D., Seli R.C. BehavePlus fire modeling system, version 2.0: User's Guide. USDA Forest Service Gen. Techn. Rep. RMRS-GTR-106WWW. Ogden, 2003, 45 p.

2. Finney M.A. FARSITE: Fire are simulator model, development and evaluation. USDA Forest Service, Res. Paper RMRS-RP-4. Ogden, 1998, 47 p.

3. Mell W., Jenkins M.A., Gould J., Cheney Ph. A physics-based approach to modeling grassland fires. Intern. Journ. of Wildland Fire. 2007, vol. 16, pp. 1-22.

4. Аксенов К.А., Гончарова Н.В. Динамическое моделирование мультиагентных процессов преобразования ресурсов. Екатеринбург: Изд-во УГТУ-УПИ, 2006. 311 с.

5. Кухта В.Б. Метод моделирования распространения низового пожара в лесных насаждениях с использованием агент-ного подхода // Вестн. Московского гос. ун-та леса - Лесной вестник. 2014. № 5. С. 92-97.

6. Dai D., Zhang Y. Simulating fire spread in a community using an agent-based model. Proc. 12th Intern. Conf. on GeoComputa-tion. Wuhan, China, LIESMARS Wuhan Univ., 2013, pp. 130-132.

7. Shatalov P.S., Dorrer G.A. Parallel computation of forest fire and its interaction with infrastructure objects. Proc. 4th Fire Behave and Fuels Conf., St. Petersburg, 2013, p. 3.

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

9. Годунов С.К., Прокопов Г.П. Об использовании подвижных сеток в газодинамических расчетах // Журн. вычислит. матем. и математ. физики. 1972. Т. 12. № 2. С. 429-439.

10. Rothermel R.C. A mathematical model for fire spread predictions in wildland fuels. USDA Forest Service Research Paper INT-115, Ogden, 1972, 40 p.

11. Волокитина А.В. Защита населенных пунктов от чрезвычайных ситуаций, связанных с природными пожарами (прак-тич. рекоменд.). Красноярск: Изд-во СО РАН, 2002. 63 с.

12. Brown D.F., Dunn W.E., Lazaro M.A., Policastro A.J. The smoke-fireplume model: tool for eventual application to prescribed burns and wildland fires. Joint Fire Sc. Program Conf. and Workshop, Boise, 1999, 12 p.

13. Мальбахов В.М., Шлычков В.А., Лежнин А.А., Дубровская О.А. Численная модель распространения дымового шлейфа при лесных пожарах с параметрическим учетом процессов горения // География и природные ресурсы. Спец. вып., 2004. С. 170-174.

14. Иванов В.А., Иванова Г.А., Москальченко С.А. Справочник по тушению природных пожаров. Красноярск: Изд-во ПРООН/МКИ, 2011. 130 с.

15. Щетинский Е.А. Тушение лесных пожаров: пособие для лесных пожарных. М.: Изд-во ВНИИЛМ, 2002. 104 с.

DOI: 10.15827/0236-235X.115.101-108 Received 14.04.16

FOREST FIRE SIMULATION USING AN AGENT-BASED MODEL 1Yarovoy S.V., Postgraduate Student, ach_bask@mail.ru

1 Siberian State Technological University, Mira Ave. 82, Krasnoyarsk, 660049, Russian Federation

Abstract. The article considers the task of modeling forest fire situation. To solve it the authors propose to use the agent-based approach with the geographic information systems (GIS) technology that allows combining the advantages

of both. As a result, an agent-based model has been developed. This model allows predicting the spread of a forest fire edge and simulating the direct method of fighting grass-roots forest fires.

Any portion of the real terrain maps can be used as an environment for agents. The proposed model uses two types of agents. Agents of the first type (A agents) are areas of fire front propagation and are used to model the behavior of forest fire. Taken together, they represent forest fire edge on a map. The model allows multiple edges, which independently spread in different parts of the map. Agents of the second type (B agents) are used for simulating firefighting actions and affect A agents.

Based on this model the authors developed a software package for a wide range of tasks. First, the system is useful for determining the operational parameters of grass-roots forest fires spread and development of measures to prevent and fight it. Second, this system can serve as a basis for forest firefighting training systems for forestry workers and the Ministry of Emergency Situations as a part of advanced professional training, as well as for teaching students of universities and colleges.

Keywords: wildfire, wildfire control, simulation modeling, geographic information system, agent-based model, software.

References

1. Andrews P.L., Bevins C.D., Seli R.C. BehavePlus fire modeling system, ver. 2.0: User's Guide. USDA Forest Service Gen. Techn. Rep. RMRS-GTR-106WWW. Ogden, 2003, 45 p.

2. Finney M.A., FARSITE: Fire are simulator model, development and evaluation. USDA Forest Service, Res. Paper RMRS-RP-4. Ogden, 1998, 47 p.

3. Mell W., Jenkins M.A., Gould J., Cheney Ph. A physics-based approach to modeling grassland fires. Int. Journ. of WildlandFire. 2007, vol. 16, pp. 1-22.

4. Aksenov K.A., Goncharova N.V. Dinamicheskoe modelirovanie multiagentnykhprotsessovpreobrazovaniya resursov [Dynamic Simulation of Multi-Agent Processes of Resource Conversion]. Ekaterinburg, UGTU-UPI Publ., 2006, 311 p. (in Russ.).

5. Kukhta V.B. Method of modeling ground fire spread in forest plantations using an agent-based approach. VestnikMos-kovskogo gos. universiteta lesa - Lesnoy vestnik [Journ. of the Moscow State Forest Univ. Forestry Journ.]. 2014, no. 5, pp. 92-97 (in Russ.).

6. Dai D., Zhang Y. Simulating fire spread in a community using an agent-based model. Proc. 12th Int. Conf. on Geo-Computation. LIESMARS Wuhan Univ. Publ., Wuhan, China, 2013, pp. 130-132.

7. Shatalov P.S., Dorrer G.A. Parallel computation of forest fire and its interaction with infrastructure objects. Proc. 4th Fire Behave and Fuels Conf. St. Petersburg, 2013, p. 3.

8. Dorrer G.A. Dinamika lesnykhpozharov [Forest Fire Dynamics]. Novosibirsk, SB RAS Publ., 2008, 404 p. (in Russ.).

9. Godunov S.K., Prokopov G.P. On the use of moving grids in the gas-dynamic calculations. Zhurn. vychislitelnoy ma-tematiki i matematicheskoy fiziki [Computational Mathematics and Mathematical Physics Journ.]. 1972, vol. 12, no. 2, pp. 429-439 (in Russ.).

10. Rothermel R.C. A mathematical model for fire spread predictions in wildland fuels. USDA Forest Service Research Paper INT-115. Ogden, 1972, 40 p.

11. Volokitina A.V. Zashchita naselennykh punktov ot chrezvychaynykh situatsy, svyazannykh s prirodnymi pozharami (prakticheskie rekomendatsii) [Protection of Inhabited Localities from Emergency Situations Connected with Wildfires (Practical Recommendations)]. Krasnoyarsk, SB RAS Publ., 2002, 63 p. (in Russ.).

12. Brown D.F., Dunn W.E., Lazaro M.A., Policastro A.J. The smoke-fireplume model: tool for eventual application to prescribed burns and wildland fires. Joint Fire Science Program Conf. and Workshop. Boise, 1999, 12 p.

13. Malbakhov V.M., Shlychkov V.A., Lezhnin A.A., Dubrovskaya O.A. Numerical model of smoke plume from forest fires with parametrical account of burning processes. Geografiya i prirodnye resursy [Geography and Natural Resources]. 2004, spec. iss., pp. 170-174 (in Russ.).

14. Ivanov V.A., Ivanova G.A., Moskalchenko S.A. Spravochnikpo tusheniyuprirodnykhpozharov [Forest Fire Guide]. Krasnoyarsk, PROON/MKI Publ., 2011, 130 p. (in Russ.).

15. Shchetinsky E.A. Tushenie lesnykh pozharov: posobie dlya lesnykh pozharnykh [Forest Fire Fighting: A Guide for Forest Firefighters]. Moscow, VNIILM publ., 2002, 104 p.

Примеры библиографического описания статьи

1. Яровой С.В. Агентный подход при моделировании лесопожарных ситуаций // Программные продукты и системы. 2016. Т. 29. № 3. С. 101-108.

2. Яровой С.В. Агентный подход при моделировании лесопожарных ситуаций // Программные продукты и системы. 2016; DOI: 10.15827/0236-235X.115.101-108.

3. Yarovoy S.V. Forest fire simulation using an agent-based model. Programmnye produkty i sistemy [Software & Systems]. 2016, vol. 29, no. 3, pp. 101-108 (in Russ.); DOI: 10.15827/0236-235X. 115.101-108.