CC BY
17
Наука та прогрес транспорту. Вюник Дншропетровського нащонального унiверситету залiзничного транспорту, 2018, № 3 (75)
ЕКОЛОГ1Я НА ТРАНСПОРТ!
УДК 502.3:504.5
Н. Н. БЕЛЯЕВ1*, И. В. КАЛАШНИКОВ2*, И. В. КЛИМЕНКО3*, В. А. КОЗАЧИНА4*
1 Каф. «Гидравлика и водоснабжение», Днепропетровский национальный университет железнодорожного транспорта имени академика В. Лазаряна, ул. Лазаряна, 2, Днипро, Украина, 49010, тел. +38 (056) 273 15 09,
эл. почта water.supply.treatment@gmail.com, ОКСГО 0000-0002-1531-7882
2 ГП «Проектно-изыскательный институт железнодорожного транспорта «Укрзализнычпроект», ул. Конарева, 7, Харьков, Украина, 61052, тел. +38 (057) 724 41 25, эл. почта uzp38@ukr.net, ОКСГО 0000-0002-2814-380Х
3 Каф. «Высшая математика», Днепропетровский национальный университет железнодорожного транспорта имени академика В. Лазаряна, ул. Лазаряна, 2, Днипро, Украина, 49010, тел. +38 (056) 273 15 09,
эл. почта water.supply.treatment@gmail.com, ОКСГО 0000-0002-1746-9118
4*Каф. «Гидравлика и водоснабжение», Днепропетровский национальный университет железнодорожного транспорта имени академика В. Лазаряна, ул. Лазаряна, 2, Днипро, Украина, 49010, тел. +38 (056) 273 15 09, эл. почта v.kozachyna@gmail.com, ОКСГО 0000-0002-6894-5532
3Б ЧИСЛЕННАЯ МОДЕЛЬ ДЛЯ ОЦЕНКИ ТЕРРИТОРИАЛЬНОГО РИСКА ПРИ ТЕРАКТЕ
Цель. Работа предполагает разработку 3Б численной модели для оценки территориального риска в случае теракта с использованием химического агента. Методика. Для описания процесса рассеивания в атмосфере химического агента, выброшенного в случае теракта, используется трехмерное уравнение массопере-носа примеси в атмосферном воздухе. Уравнение учитывает поле скорости ветрового потока, атмосферную диффузию, интенсивность эмиссии химического агента, наличие зданий возле места выброса химически опасного вещества. Для численного интегрирования моделирующего уравнения используется конечно -разностный метод. Особенностью разработанной численной модели является возможность оценки территориального риска в случае теракта при различных метеоусловиях и наличии зданий. Для расчета поля скорости ветрового потока в условиях застройки применяется трехмерное уравнение для потенциала скорости. Результаты. Разработанная численная модель и пакет программ могут быть использованы для оценки территориального риска как в случае терактов с применением химических агентов, так и в случае экстремальных ситуаций на химически опасных объектах и транспорте. Построенная численная модель может быть реализована на компьютерах малой и средней мощности, что позволяет широко использовать ее для решения задач рассматриваемого класса, при разработке плана ликвидации аварийной ситуации. Представлены результаты вычислительного эксперимента, позволяющие оценить возможности предложенного метода оценки территориального риска в случае теракта с использованием химического агента. Научная новизна. Предложена эффективная численная модель для оценки территориального риска в случае теракта с применением химически опасного вещества. Метод может быть использован для оценки территориального риска в условиях городской застройки, что позволяет получать адекватные данные о возможных зонах поражения. Метод основан на численном интегрировании фундаментального уравнения мас-сопереноса, выражающего закон сохранения массы в жидкой среде. Практическая значимость. Предложенный метод оценки территориального риска в случае теракта с применением химического агента может быть использован для расчета зон поражения возле административных зданий, центров и других социально значимых объектов.
Ключевые слова: территориальный риск; теракт; химическое загрязнение; численное моделирование; загрязнение воздушной среды
Наука та прогрес транспорту. Вюник Дншропетровського нащонального унiверситету залiзничного транспорту, 2018, № 3 (75)
Введение
Теракты с применением химических (биологических) агентов представляют собой уже реальную угрозу для населения разных стран [1-4, 6, 7-9, 11]. В качестве агентов могут быть использованы как химически опасные вещества, применяемые в различных отраслях промышленности (хлор, аммиак и т. под.), так и специальные, разработанные для боевого применения. Если рассматривать теракт в селитебной зоне, то особенностью его является ограниченный по массе и по времени выброс агента, т. е. эмиссия - полунепрерывная. Это связано с тем, что незаметно поставить в селитебную зону большие емкости с химическим (биологическим) агентом достаточно сложно. Например, емкостью для опасного вещества может являться баллон для топлива, используемый на автозаправочных станциях, но заполненный поражающим агентом.
Если теракт имеет место на промышленном объекте или в транспорте, то выброс опасного вещества является более мощным и происходит значительно дольше по времени. В случае теракта, сопровождающегося эмиссией химического (биологического) агента, в атмосфере формируется зона загрязнения, интенсивность которой в различных точках области различна. В определенных местах и в определенные моменты времени значение концентрации агента будет выше некоторого порогового значения (например, смертельной концентрации). Эта поражающая концентрация будет зависеть в первую очередь от мощности выброса и от метеоусловий.
Актуальной проблемой является прогноз риска поражения людей в селитебной зоне при возможных терактах с применением химических (биологических) агентов. Такой прогноз может быть выполнен только методом математического моделирования.
Цель
Основной целью данной работы является разработка 3Б численной модели для оценки территориального риска при террористической атаке, сопровождающейся эмиссией химического агента в условиях застройки.
В настоящее время в Украине для оценки размеров и интенсивности зон химического за-
ражения используют методику ОНД-86, методику прогноза последствий аварий на химически опасных объектах и транспорте, модель Гаусса. За рубежом оценку зон поражения проводят чаще всего по модели Гаусса [4, 10], методом СББ-моделирования [2, 3, 12]. Следует подчеркнуть, что применение СББ-моделирования дает возможность получать прогнозные данные с учетом такого важного фактора, как наличие зданий в зоне рассеивания опасного вещества. За рубежом СББ-модели основаны на применении уравнений Навье-Стокса совместно с определенной моделью турбулентности, что приводит к значительным затратам компьютерного времени, даже при реализации моделей на мощных зарубежных компьютерах. Расчет одного варианта задачи может длиться несколько суток, что неудобно при проведении серийных расчетов.
Методика
Создание 3Б СББ-модели необходимо для расчета аэродинамики обтекания зданий и оценки территориального риска при теракте на урбанизированной территории.
При оценке территориального риска будем учитывать, что зона заражения различна для различных метеоситуаций Р (^ ) . Под определенной метеоситуацией будем понимать конкретное значение скорости и направления ветра.
Вероятность реализации конкретной метеоситуации в регионе известна и определяется по зависимости [1]:
Р(Щ) = Нп/Г , (1)
где НП - число дней (часов), соответствующих определенной метеоситуации; Т - период наблюдений.
Для оценки потенциального территориального риска необходимо оценить вероятность для человека, находящегося в каждой точке области возле атакуемого объекта, оказаться под действием шлейфа химически опасного вещества.
Вероятность попадания человека под действие шлейфа химически опасного вещества определяется следующим образом:
Наука та прогрес транспорту. Вюник Дншропетровського нащонального унiверситету залiзничного транспорту, 2018, № 3 (75)
P P (W
(2)
где P(W),
суммарная вероятность всех рас-
сматриваемых метеоситуации, при которых человек попадает в зону воздействия источника эмиссии и получает токсичное заражение.
Расчет зоны заражения для каждого (вероятного) метеоусловия выполняется на базе трехмерного уравнения массопереноса [2, 3, 5, 8]:
д(w - wS)C
+1 а (()5( X - хг (Г) )5( .у - у, (Г) )5( 7 - ), (3)
где С - концентрация химического (биологического) агента в атмосферном воздухе; с -коэффициент, учитывающий распад агента; и, V, - компоненты вектора скорости воздушного потока; - скорость гравитационного оседания примеси; ц, = (ц,х, цу, ) - коэффициенты атмосферной турбулентной диффузии; а - интенсивность выброса агента;
5( х - х,)(у - у,)( г - г,) - дельта-функция Дирака; х,, у,, г, - координаты источника эмиссии агента; ^ - время.
Поскольку рассматривается теракт в селитебной зоне или в промышленной зоне вблизи зданий, то возникает необходимость расчета поля скорости ветрового потока (параметры и, V, ^). Для этого расчета при наличии зданий будем использовать модель потенциального течения:
д2 P д2 P д2 P
дк2 ду2
дг2
= 0,
(4)
где Р - потенциал скорости.
После определения поля потенциала скорости осуществляется расчет компонент вектора скорости воздушного потока:
u = -
ÖP_
дк '
V = -
ÖP_ ду
w = -
ÖP_
Постановка краевых условий для моделирующих уравнений (3, 4) рассмотрена в [2, 3, 5].
Для численного интегрирования уравнения (4) используется метод Ричардсона, а для численного интегрирования уравнения переноса примеси применяется неявная разностная схема расщепления [2, 5].
Расчет территориального риска. Таким образом, при прогнозировании химической атаки террориста необходимо для конкретной точки конкретной области (поля) выполнить расчеты по формуле (5). Для их выполнения предварительно необходимо рассчитать значение концентрации химического агента в точке расположения конкретного человека, для конкретной метеоситуации, и определить, насколько эта величина превосходит заданный пороговый уровень.
Алгоритм решения. Оценку территориального риска (прогноз последствий) при теракте и при вероятной метеорологической ситуации РЖ будем осуществлять в такой последовательности [3]:
1) на первом этапе решения задачи формируется блок данных относительно инициирующего события (возможное место эмиссии химически опасного вещества, интенсивность, режим эмиссии, вид химического агента);
2) на втором этапе формируется блок данных относительно вероятных метеоситуаций РЖ,, характерных для области, где находится атакуемый объект;
3) на третьем этапе рассчитывается уровень химического заражения для вероятных метеоситуаций (на этом этапе проводится численное интегрирование уравнений (1) и (3) для конкретной метеоситуации);
4) на четвертом этапе определяются зоны, где концентрация превышает пороговое значение (например, смертельную концентрацию) при конкретной метеоситуации;
5) на пятом этапе осуществляется построении поля риска для рассматриваемого объекта.
Данный алгоритм оценки риска реализован в разработанном коде «Ы8К-3Б Q».
Результаты
Разработанный код был использован для решения следующей модельной задачи. Рассматривается сценарий атаки, когда возле двух
зданий происходит эмиссия хлора в двух местах. Полагается, что на момент атаки возможны такие значения скорости ветра: 3 м/с и 6 м/с.
На рис. 1 показаны изолинии концентрации хлора для момента времени 110 с, при скорости ветра 6 м/с (уровень ъ = 2м). Из данного рисунка видно, что при эмиссии химически опасного вещества формируется зона заражения, которая охватывает оба здания и создает угрозу для людей, находящихся на любой стороне зданий.
0.340Е+01 coordiïiate х 0.167Е+03
Рис. 1. Зона химического заражения при атаке террориста
Fig. 1. Zone of chemical contamination after the chemical attack of a terrorist
Это особенно хорошо видно из рис. 2, где представлена матрица территориального риска в районе химической атаки (момент времени 110 с) при вероятных метеоситуациях. Вероятность заражения людей возле здания показана в процентах. Принимается, что заражение наступает, если концентрация хлора в расчетной точке превышает 3 мг/м3.
Как видно из рис. 2, зона заражения (риск до 99 %) формируется не только возле зданий, но и за ними. То есть для рассматриваемого сценария здания не выполняют защитную функцию барьера.
Отметим, что время расчета составляет порядка 10 сек.
Наука та прогрес транспорту. Вюник Дншропетровського нащонального ушверситету з&тзничного транспорту, 2018, N° 3 (75)
Рис. 2. Матрица территориального риска (уровень 3 м) при химической атаке террориста
Fig. 2. Matrix of territorial risk (level z=2m) after the chemical attack of a terrorist
Научная новизна и практическая значимость
Разработана 3D численная модель, позволяющая оценить величину территориального риска в случае теракта с применением химического агента в условиях застройки (селитебная зона или промышленная площадка).
Особенностью построенной модели является использование стандартной исходной информации, быстрота в получении прогнозных данных и удобство для анализа получаемых результатов.
Выводы
В работе представлена 3D-модель для оценки территориального риска при атаке террориста с применением химически опасного вещества и рассеиванием химического агента в условиях застройки. В основу расчета риска положено численное моделирование рассеивания опасного вещества в условиях застройки с последующей оценкой размеров зон поражения для каждой вероятной метеоситуации.
Дальнейшее совершенствование выбранного научного направления следует проводить в области создания SD-модели для расчета территориального риска на базе уравнений, описывающих вязкое движение воздушного потока.
СПИСОК ИСПОЛЬЗОВАННЫХ ИСТОЧНИКОВ
1. Алымов, В. Т. Техногенный риск. Анализ и оценка : учеб. пособие для вузов / В. Т. Алымов, Н. П. Тарасова. - Москва : Академкнига, 2004. - 118 с.
2. Беляев, Н. Н. Защита зданий от проникновения в них опасных веществ : монография / Н. Н. Беляев, Е. Ю. Гунько, Н. В. Росточило. - Днепропетровск : Акцент ПП, 2014. - 136 с.
Наука та прогрес транспорту. Вюник Дншропетровського нащонального унiверситету залiзничного транспорту, 2018, N° 3 (75)
3. Оценка техногенного риска при эмиссии опасных веществ на железнодорожном транспорте / Н. Н. Беляев, Е. Ю. Гунько, П. С. Кириченко, Л. Я. Мунтян. - Кривой Рог : Р. А. Козлов, 2017. - 127 с.
4. Стоецкий, В. Ф. Оценка риска при авариях техногенного характера / В. Ф. Стоецкий, В. И. Голинько, Л. В. Дранишников // Наук. вюн. НГУ. - 2014. - № 3. - С. 117-124.
5. Численное моделирование распространения загрязнения в окружающей среде / М. З. Згуровский, В. В. Скопецкий, В. К. Хрущ, Н. Н. Беляев. - Киев : Наук. думка, 1997. - 368 с.
6. Barret, A. M. Mathematical Modeling and Decision Analysis for Terrorism Defense: Assessing Chlorine Truck Attack Consequence and Countermeasure Cost Effectiveness : Degree of Doctor of Philosophy / Anthony Michael Barret ; Carnegie Mellon University. - Pittsburg, Pennsylvania, 2009. - 123 p.
7. Berlov, O. V. Atmosphere protection in case of emergency during transportation of dangerous cargo / O. V. Berlov // Наука та прогрес транспорту. - 2016. - № 1 (61). - С. 48-54. doi: 10.15802/stp2016/60953
8. Biliaiev, M. M. Numerical Simulation of Indoor Air Pollution and Atmosphere Pollution for Regions Having Complex Topography / M. M. Biliaiev, M. M. Kharytonov // NATO Science for Peace and Security. Series C: Environmental Security. - Dordrecht, 2012. - P. 87-91. doi: 10.1007/978-94-007-1359-8_15
9. Cefic Guidance on safety Risk Assessment for Chemical Transport Operations J Verlinden [Электронный ресурс]. - Режим доступа: http://www.era.europa.eu/Document-Register/Documents/Cefic%20guidance%20on%20risk%20assessment.pdf - Загл. с экрана. - Проверено : 29.01.2018.
10. Effect of barriers on the status of atmospheric pollution by mathematical modeling / Zahra Naserzadeh, Farideh Atabi, Faramarz Moattar, Naser Moharram Nejad // Bioscience Biotechnology Research Communication. - 2017. - Vol. 10 (1). - P. 192-204.
11. Protective Action Criteria. A Review of Their Derivation, Use, Advantages and Limitations [Электронный ресурс] / Environmental Public Health Science Unit, Health Protection Branch, Public Health and Compliance Division, Alberta Health. - Edmonton, Alberta, 2017. - Режим доступа: http://open.alberta.ca/publications/9781460131213 - Загл. с экрана. - Проверено : 14.06.2018.
12. The analysis of the use of mathematical modeling for emergency planning purposes / Ondrej Zavila, Pavel Dobes, Jakub Dlabka, Jan Bitta // Bezpecnostni vyzkum. The Science for Population Protection. - 2015. -№ 2. - P. 1-9.
М. М. ЫЛЯСВ1*, I. В. КАЛАШНЖОВ2*, I. В. КЛИМЕНКО3*, В. А. КОЗАЧИНА4*
1 Каф. «Пдравлжа та водопостачання», Дшпропетровський нацюнальний ушверситет залiзничного транспорту iменi академжа В. Лазаряна, вул. Лазаряна, 2, Дшпро, Украша, 49010, тел. +38 (056) 273 15 09,
ел. пошта water.supply.treatment@gmail.com, ORCID 0000-0002-1531-7882
2 Держ. шдприемство «Проектно-вишукувальний шститут залiзничного транспорту Украши «Укрзалiзничпроект», вул. Конарева, 7, Харкгв, Украша, 61052, тел. +38 (057) 724 41 25, ел. пошта uzp38@ukr.net,
ORCID 0000-0002-2814-380X
3 Каф. «Вища математика», Дшпропетровський нацюнальний утверситет залiзничного транспорту iменi академжа В. Лазаряна, вул. Лазаряна, 2, Дншро, Украша, 49010, тел. +38 (056) 273 15 09,
ел. пошта water.supply.treatment@gmail.com, ORCID 0000-0002-1746-9118
4 Каф. «Пдравлжа та водопостачання», Дшпропетровський нацюнальний ушверситет залiзничного транспорту iменi академжа В. Лазаряна, вул. Лазаряна 2, Дншро, Украша, 49010, тел. +38 (056) 273 15 09,
ел. пошта v.kozachyna@gmail.com, ORCID 0000-0002-6894-5532
3D ЧИСЕЛЬНА МОДЕЛЬ ДЛЯ ОЦ1НКИ ТЕРИТОР1АЛЬНОГО РИЗИКУ ПРИ ТЕРАКТ
Мета. Робота передбачае розробку 3D чисельно! моделi для оцшки територiального ризику в разi теро-ристичного нападу з використанням хiмiчного агента. Методика. Для опису процесу розсшвання в атмос-ферi хiмiчного агента, викинутого у випадку теракту, використовуеться рiвняння масопереносу домшки в атмосферному повггрг Рiвняння враховуе швидшсть вггрового потоку, атмосферну дифузш, штенсившсть емюп хiмiчного агента, наявшсть будiвель б™ мюця викиду хiмiчно небезпечно! речовини. Для чисельного штегрування моделюючого рiвняння використовуеться кшцево^зницевий метод. Особливютю розроблено! чисельно! моделi е можливють оцшки територiального ризику в разi теракту при рiзних метеоумовах та ная-вносп будiвель. Для розрахунку поля швидкосл впрового потоку в умовах забудови використовуеться три-вимiрне рiвняння для потенщалу швидкостi. Результати. Розроблена спецiалiзована чисельна модель та
Наука та прогрес транспорту. Вюник Дншропетровського нащонального ушверситету залiзничного транспорту, 2018, № 3 (75)
пакет програм можуть бути використанi для оцiнки TeprnopiaibHoro ризику як у випадку теракпв i3 застосу-ванням хiмiчних агенпв, так i в paзi екстремальних ситуaцiй на xiMi4HO небезпечних об'ектах i тpaнспopтi. Метод може бути peaлiзoвaний на комп'ютерах мало! та середньо! пoтужнoстi, що дозволяе широко викори-стовувати його для виршення задач даного класу. Пpeдстaвлeнi результати обчислювального експерименту, що дозволяють оцшити мoжливoстi запропонованого методу oцiнки тepитopiaльнoгo ризику в paзi терорис-тичного нападу з використанням xiмiчнoгo агента. Наукова новизна. Запропоновано ефективний метод оцшки тepитopiaльнoгo ризику в pasi теракту зi застосуванням xiмiчнo небезпечно! речовини. Метод може бути використаний для оцшки тepитopiaльнoгo ризику в умовах мюько! забудови, що дозволяе отримувати aдeквaтнi дaнi про мoжливi зони ураження. Метод заснований на чисельному iнтeгpувaннi фундаментального piвняння масопереносу, що виражае закон збереження маси в рщкому сepeдoвищi. Практична значимкть. Запропонований метод oцiнки тepитopiaльнoгo ризику в paзi терористичного нападу зi застосуванням xiмiчнoгo агента може бути використаний для розрахунку зон ураження б™ адмшют-ративних будiвeль, цeнтpiв та iншиx сoцiaльнo значущих oб'ектiв.
Ключовi слова: тepитopiaльний ризик; теракт; xiмiчнe забруднення; чисельне моделювання; забруднення пoвiтpянoгo середовища
M. M. BILIAIEV1*, I. V. KALASHNIKOV2*, I. V. KLYMENKO3*, V. A. KOZACHYNA4*
1 Dep. «Hydraulics and Water Supply», Dnipropetrovsk National University of Railway Transport named after Academician V. Lazaryan, Lazaryan St., 2, Dnipro, Ukraine, 49010, tel. +38 (056) 273 15 09, e-mail Water.supply.treatment@gmail.com, ORCID 0000-0002-1531-7882
2 State Enterprise «Design and Exploration Institute of Railway Transport of Ukraine «Ukrzaliznichproekt», Konarev St., 7, Kharkiv, 61052, tel. +38 (057) 724 41 25, e-mail uzp38@ukr.net, ORCID 0000-0002-2814-380X
3 Dep. «Higher Mathematics», Dnipropetrovsk National University of Railway Transport named after Academician V. Lazaryan, Lazaryan St., 2, Dnipro, Ukraine, 49010, tel. +38 (056) 273 15 09, e-mail Water.supply.treatment@gmail.com, ORCID 0000-0002-1746-9118
4 Dep. «Hydraulics and Water Supply», Dnipropetrovsk National University of Railway Transport named after Academician V. Lazaryan, Lazaryan St., 2, Dnipro, Ukraine, 49010, tel. +38 (056) 273 15 09,
e-mail v.kozachyna@gmail.com, ORCID 0000-0002-6894-5532
3D NUMERICAL MODEL FOR TERRITORIAL RISK ASSESMENT AFTER TERRORIST ACT
Purpose. The paper deals with the development of 3D numerical model for a territorial risk assessment in a case of a terrorist attack with usage of a chemical agent. Methodology. To describe the process of chemical agent dispersion in the atmosphere, emitted in a case of a terrorist attack, authors used three-dimensional equation of mass transfer of an admixture in the atmospheric air. The equation takes into account a velocity of the wind flow, atmospheric diffusion, an intensity of chemical agent emission, a buildings presence near a place of a chemically hazardous substance release. For numerical integration of the modeling equation, a finite difference method is used. A feature of the developed numerical model is possibility of a territorial risk assessment in a case of a terrorist attack under different weather conditions and buildings presence. For calculation of wind flow velocity field in a case of building presence 3D equation for potential of velocity is used. Findings. The specialized numerical model and software package can be used for a territorial risk assessment, both in a case of terrorist attacks with usage of chemical agents and in a case of extreme situations at chemically hazardous facilities and transport. The method can be implemented on small and medium-powered computers, which allows it to be widely used for solving this class of problems, when developing an emergency response plan. The results of a computational experiment are presented that allow estimating possibilities of the proposed method for assessing a territorial risk in a case of a terrorist attack with usage of chemical agent. Originality. The effective method of a territorial risk assessment in a case of a terrorist attack with usage of a chemically hazardous substance is proposed. The method can be used to assess a territorial risk in an urban environment, which allows obtaining adequate data about possible affected areas. The method is based on numerical integration of the fundamental mass transfer equation, which expresses the law of conservation of mass in a liquid medium. Practical value. The proposed method for a territorial risk assessment in a case of a terrorist attack with usage of a chemical agent can be used to calculate the affected areas near administrative buildings, centers and other socially significant facilities.
Keywords: territorial risk; terrorist act; chemical pollution; numerical modelling; air pollution
Наука та прогрес транспорту. Вюник Дншропетровського нащонального ушверситету залiзничного транспорту, 2018, № 3 (75)
REFERENCES
1. Alymov, V. T., & Tarasova, N. P. (2004). Tekhnogennyy risk. Analiz i otsenka: Uchebebnoe posobie dlya vuzov. Moscow: Akademkniga. (in Russian).
2. Belyaev, N. N., Gunko, Y. Y., & Rostochilo, N. V. (2014). Zashchita zdaniy ot proniknoveniya v nikh opasnykh veshchestv: Monografiya. Dnepropetrovsk: Aktsent PP. (in Russian).
3. Belyaev, N. N., Gunko, Y. Y., Kirichenko, P. S., & Muntyan, L. Y. (2017). Otsenka tekhnogennogo riska pri emissii opasnykh veshchestv na zheleznodorozhnom transporte. Krivoi Rog: Kozlov R. A. (in Russian).
4. Stoetsky, V. F., Golinko, V. I., & Dranishnikov, L. V. (2014). Risk assessment in man-caused accidents. Scientific Bulletin of National Mining University, 3, 117-124. (in Russian).
5. Zgurovskiy, M. Z., Skopetskiy, V. V., Khrushch, V. K., Belyaev, N. N. (1997). Chislennoe modelirovanie rasprostraneniya zagryazneniya v okruzhayushchey srede. Kiev: Naukova dumka. (in Russian).
6. Barret, A. M. (2009). Mathematical Modeling and Decision Analysis for Terrorism Defense: Assessing Chlorine Truck Attack Consequence and Countermeasure Cost Effectivness. (Dissertation of Doctor of Philosophy). Carnegie Mellon University, Pittsburg, Pennsylvania, USA. (in English)
7. Berlov, O. V. (2016). Atmosphere protection in case of emergency during transportation of dangerous cargo. Sciance and Transport Progress, 1(61), 48-54. doi: 10.15802/stp2016/60953 (in English)
8. Biliaiev, M. M., & Kharytonov, M. M. (2012). Numerical Simulation of Indoor Air Pollution and Atmosphere Pollution for Regions Having Complex Topography. NATO Science for Peace and Security. Series C: Environmental Security. doi: 10.1007/978-94-007-1359-8_15 (in English)
9. Cefic Guidance on safety Risk Assessment for Chemical Transport Operations J Verlinden. (n.d.). Retrived from http://www.era.europa.eu/DocumentRegister/Documents/Cefic%20guidance%20on%20risk%20asses sment.pdf (in English)
10. Zahra Naserzadeh, Farideh Atabi, Faramarz Moattar, & Naser Moharram Nejad. (2017). Effect of barriers on the status of atmospheric pollution by mathematical modeling. Bioscience Biotechnology Research Communication, 10(1), 192-204. (in English)
11. Government of Alberta. Protective Action Criteria: A Review of Their Derivation, Use, Advantages and Limitations. Environmental Public Health Science Unit, Health Protection Branch, Public Health and Compliance Division, Alberta Health. Edmonton, Alberta. Retrived from http://open.alberta.ca/publications/9781460131213 (in English)
12. Ondrej Zavila, Pavel Dobes, Jakub Dlabka, & Jan Bitta. (2015). The analysis of the use of mathematical modeling for emergency planning purposes. Bezpecnostni vyzkum. The Science for Population Protection, 2, 1-9. (in English)
Статья рекомендована к публикации д.физ.-мат.н., проф. C. А. Пичуговым (Украина)
Поступила в редколегию: 12.02.2018
Принята к печати: 30.05.18
