CC BY
26
Наука та прогрес транспорту. Вюник Дншропетровського нащонального унiверситету залiзничного транспорту, 2018, № 1 (73)
ЕКОЛОГ1Я НА ТРАНСПОРТ!
УДК 502.3:504.5
Н. Н. БЕЛЯЕВ1*, И. В. КАЛАШНИКОВ2*, В. А. КОЗАЧИНА3
2*
3*
'*Каф. «Гидравлика и водоснабжение», Днепропетровский национальный университет железнодорожного транспорта имени академика В. Лазаряна, ул. Лазаряна, 2, Днипро, Украина, 49010, тел. +38 (056) 273 15 09, эл. почта [email protected], ОЯСГО 0000-0002-1531-7882
2 Гос. предприятие «Проектно-изыскательный институт железнодорожного транспорта «Укрзализнычпроект», ул. Красноармейская, 7, Харьков, Украина, 61052, тел. +38 (057) 724 41 25, эл. почта [email protected],
ОЯСГО 0000-0002-2814-380Х
3 Каф. «Гидравлика и водоснабжение», Днепропетровский национальный университет железнодорожного транспорта имени академика В. Лазаряна, ул. Лазаряна, 2, Днипро, Украина, 49010, тел. +38 (056) 273 15 09,
эл. почта [email protected], ОЯСГО 0000-0002-6894-5532
РАСЧЕТ ТЕРРИТОРИАЛЬНОГО РИСКА ПРИ ТЕРАКТЕ: ЭКСПРЕСС МОДЕЛЬ
Цель. Работа предполагает разработку метода оценки территориального риска в случае теракта с использованием химического агента. Методика. Для описания процесса рассеивания в атмосфере химического агента, выброшенного в случае теракта, используется уравнение массопереноса примеси в атмосферном воздухе. Уравнение учитывает скорость ветрового потока, атмосферную диффузию, интенсивность эмиссии химического агента, наличие зданий возле места выброса химически опасного вещества. Для численного интегрирования моделирующего уравнения используется конечно-разностный метод. Особенностью разработанной численной модели является возможность оценки территориального риска в случае теракта при различных метеоусловиях и наличии зданий. Результаты. Разработана специализированная численная модель и пакет программ, которые могут быть использованы для оценки территориального риска как в случае терактов с применением химических агентов, так и в случае экстремальных ситуаций на химически опасных объектах и транспорте. Метод может быть реализован на компьютерах малой и средней мощности, что позволяет широко использовать его для решения задач рассматриваемого класса. Представлены результаты вычислительного эксперимента, позволяющие оценить возможности предложенного метода оценки территориального риска в случае теракта с использованием химического агента. Научная новизна. Предложен эффективный метод оценки территориального риска в случае теракта с применением химически опасного вещества. Метод может быть использован для оценки территориального риска в условиях городской застройки, что позволяет получать адекватные данные о возможных зонах поражения. Метод основан на численном интегрировании фундаментального уравнения массопереноса, выражающего закон сохранения массы в жидкой среде. Практическая значимость. Предложенный метод оценки территориального риска в случае теракта с применением химического агента может быть использован для расчета зон поражения возле административных зданий, центров и других социально значимых объектов.
Ключевые слова: территориальный риск; теракт; химическое загрязнение; численное моделирование; загрязнение воздушной среды
В настоящее время за рубежом повышенное внимание уделяется прогнозу последствий террористических актов с возможным применением химических или биологических агентов [1,
Введение
2, 3, 8, 9, 11]. Как отмечают зарубежные исследователи, сейчас высока доступность для террористов тех химически опасных веществ, которые используются в различных отраслях промышленности (например, хлор, аммиак и т.п.). В этой связи, за рубежом при изучении проблемы оценки последствий терактов часто
Наука та прогрес транспорту. Вюник Дншропетровського нащонального унiверситету залiзничного транспорту, 2018, № 1 (73)
рассматривается эмиссия именно таких химически опасных веществ и, реже, боевых отравляющих веществ [6, 7].
Известно, что теракты чаще происходят в условиях застройки (промышленной застройки или в селитебной зоне). В этой связи, для оценки последствий возможного теракта в создаваемых математических моделях обязательно следует учесть влияние зданий на рассеивание химического (биологического) агента. Однако, следует отметить, что такой учет крайне осложняет решение задачи, поскольку для прогноза возможных зон поражения следует знать поле скорости ветрового потока между зданиями, которое надо рассчитать для всех возможных метеоситуаций. Для решения данной задачи необходимо решить аэродинамическую задачу и определить это поле скорости.
Цель
Целью данной работы является разработка метода оценки территориального риска при террористической атаке, сопровождающейся эмиссией химического агента в условиях застройки. Метод расчета основывается на создании численной модели рассеивания химического агента при теракте.
Анализ публикаций. В настоящее время в Украине не разрабатываются методы по оценке последствий терактов (с применением химических или биологических агентов) в условиях застройки. Существуют публикации, где используются нормативные методики типа ОНД-86 или методика прогноза последствий аварий на химически опасных объектах и транспорте для оценки уровня загрязнения воздушной среды при выбросе химически опасных веществ. Для оценки такого уровня загрязнения при авариях на производствах используются методики типа «Токси» [4], «Аммиак», «SLAB» [6], реализующие модель Гаусса, или аналитическое решение уравнения массо-переноса, которые также не могут быть применены, поскольку не учитывают влияние зданий на рассеивание химического (биологического) агента. Одним из подходов к оценке территориального риска при экстремальных ситуациях на транспорте является применение численных кинематических моделей и др. [3, 4, 10, 12]. Следует отметить, что применение данных мо-
делеи ограничено только случаем рассеивания химически опасных веществ над ровноИ поверхностью. Следовательно, применять данные модели для оценки последствии терактов в условиях застроИки нельзя.
Методика
Методика - создание 2D численной модели для расчета аэродинамики обтекания зданий и оценки территориального риска при теракте на урбанизтрованной территории.
Моделирующие уравнения. Для создания адекватной модели прогноза последствий теракта с применением химического (биологического) агента необходимо учесть неравномерное поле скорости ветра при наличии застройки. В этой связи для расчета поля скорости ветрового потока при наличии зданий будем использовать модель безвихревого течения идеальной жидкости:
д2Р д2Р п
—Г+—Т = 0 ■
дх2 ду
(1)
где Р - потенциал скорости. Компоненты вектора скорости ветрового потока определяются соотношениями:
дР дР u =—, V = — дх ду
(2)
Для математического моделирования процесса распространения химического агента (биологического агента) в атмосферном воздухе применяется уравнение массопереноса, осред-ненное по высоте переноса примеси [2, 3, 5, 8]:
дС диС дчС , , 7„ч
--1---1---VоС = div\ugradL) +
дх ду
N
+Ев. ((Жх - X)5(у - у), (3)
1=1
где С - концентрация химического (биологического) агента в атмосферном воздухе; о - коэффициент, учитывающий распад агента; и, V - компоненты вектора скорости воздушного потока; ¡и = {¡,х ,иу ) - коэффициенты атмосферной турбулентной диффузии; в - интенсивность выброса химического (биологического) агента; д(х - х )(у - у ) - дельта-функция
Наука та прогрес транспорту. Вюник Дншропетровського нацюнального унiверситету залiзничного транспорту, 2018, № 1 (73)
Дирака; х,, у - координаты источника эмиссии агента; ^ - время.
Таким образом, для расчета вероятности попадания в зону поражения (то есть для оценки величины территориального риска) необходимо выполнить расчеты для конкретного рецептора по формуле (3). Необходимо подчеркнуть, что для решения поставленной задачи надо знать концентрацию опасного вещества в точке расположения конкретного рецептора и насколько эта величина превосходит некоторый пороговый уровень (например, ПДК или какую-либо другую величину). Следовательно, мы приходим к первоочередной задаче - определению поля концентрации опасного вещества для конкретного региона (территории) и конкретной метеоситуации.
Поле скорости ветрового потока в условиях застройки по моделирующему уравнению (3) определяется на основе решения аэродинамической задачи, т. е. путем решения уравнения (1). Постановка краевых условий для моделирующих уравнений (1), (2) рассмотрена в [2, 3, 5].
Для численного интегрирования уравнения (1) используется метод Либмана, а для численного интегрирования уравнения переноса примеси применяется неявная разностная схема расщепления [2, 5].
Расчет территориального риска. Для оценки территориального риска при теракте с применением химического агента будем учитывать, что каждому погодному состоянию )
отвечает конкретная зона загрязнения, характеризующаяся размерами и концентрацией опасного вещества. Вероятность реализации конкретной метеоситуации определяется по зависимости [1]:
Р(Щ )= кп/т, (4)
где Nп - число дней (часов), соответствующих определенной метеоситуации; Т - период наблюдений (прогноз метеоситуации).
Здесь под определенной метеоситуацией будем понимать конкретное значение скорости и направления ветра. Для оценки потенциального территориального риска необходимо оценить вероятность для человека, находящегося в каждой точке области, возле атакуемого объекта, оказаться под действием шлейфа (облака)
химически опасного вещества. При прогнозе нас интересуют те точки области в районе теракта, где концентрация химически опасного вещества превышает некоторое пороговое значение, при котором происходит та или иная степень поражения людей. Схематически на рис.1 показана ситуация попадания рецепторов А, В, С под действие источника эмиссии (баллона, в котором находится токсичный газ), причем отметим, что выброс происходит в условиях застройки.
скоростъ ПС 1)1» V
1>(W )
Рис. 1. Вероятность появления зон поражения для различных метеоусловий: 1 - здание
Fig. 1. Probability of occurrence of damage areas for various weather conditions: 1 - building
Как видно из рис. 1, человек, находящийся в точке А, попадает под воздействие шлейфа опасного вещества при любой метеоситуации. Человек, который расположен в точке В, попадает под воздействие шлейфа опасного вещества при метеоситуации PW2, PW3. Человек, находящийся в точке С, попадает под воздействие шлейфа опасного вещества только при метеоситуации PW3. Таким образом, вероятность попадания человека под действие шлейфа химически опасного вещества определится следующим образом:
n
p(wh =1 pW ), (5)
7=0
Наука та прогрес транспорту. Вюник Дншропетровського нащонального ушверситету з^зничного транспорту, 2018, N° 1 (73)
где P(W)
Z
суммарная вероятность всех рас-
сматриваемых метеоситуаций, при которых человек попадает в зону воздействия источника эмиссии и получает токсичное поражение.
Таким образом, при прогнозировании химической атаки террориста необходимо для конкретной точки конкретной области (поля), выполнить расчеты по формуле (5). Для их выполнения предварительно необходимо рассчитать значение концентрации химического агента в точке расположения конкретного человека, для конкретной метеоситуации, и определить, насколько эта величина превосходит заданный пороговый уровень.
Алгоритм решения. Оценку территориального риска (прогноз последствий) при теракте и при вероятной метеорологической ситуации РЩ будем осуществлять в такой последовательности [3]:
1) на первом этапе решения задачи формируется блок данных относительно инициирующего события (возможное место эмиссии химически опасного вещества, интенсивность эмиссии, режим эмиссии, вид химического агента);
2) на втором этапе формируется блок данных относительно вероятных метеоситуаций РЩ, характерных для области, где находится атакуемый объект;
3) на третьем этапе рассчитывается уровень химического заражения для вероятных метеоситуаций (на этом этапе проводится численное интегрирование уравнений (1) и (3) для конкретной метеоситуации);
4) на четвертом этапе определяются зоны, где концентрация превышает пороговое значение (например, смертельную концентрацию) при конкретной метеоситуации;
5) на пятом этапе осуществляется построении поля риска для рассматриваемого объекта.
Данный алгоритм оценки риска реализован в разработанном коде «СИЕМ-ШБК».
Результаты
Разработанный код был использован для решения задачи методом моделирования. На рис. 2, а представлена расчетная область. В этой области располагаются несколько зда-
ний. Предполагалось, что происходит теракт с эмиссией хлора. Выброс происходит между двумя зданиями. Объект атаки - фронтально расположенная зона возле последнего здания. Длительность выброса - 6 мин. Для рассматриваемого региона вероятные значения скорости ветра составляют 2 м/с (вероятность 54 %) и 5 м/с (вероятность 30 %). Стрелка на рисунке показывает направление ветра. Рассматриваются два сценария. Первый сценарий - отсутствие защиты возле объекта атаки, второй - присутствие защиты типа «трейлер», которая размещается перед объектом атаки.
На рис. 2, а, б показаны изолинии концентрации хлора для момента времени 190 с, при скорости ветра 5 м/с, без защиты и при защите. Из данных рисунков видно, что формируется обширная зона химического заражения, которая охватывает объект атаки.
а - а
В.ЭВ6Е»01 о B.9ZZE>81 г В.В5ЭЕ-01 d В.795Е.В1 i В.731Е+В1 n 9. GGSE+81 а Я.6В4М1 t B.S41E»B1 в В.47ГС.81 В.413Е>01 8.35@Е<01 у 0. 2S6E+01 В.223Е>81 В.159Е.81 B.S54E.B8 0.318Е+00
И i.
ÜJ]
б - b
Рис. 2. Зона химического заражения при химической атаке террориста возле объекта атаки:
а - без защиты, б - при защите
Fig. 2. Zone of chemical contamination during a chemical attack of a terrorist:
a - no protection near the target object, b - with protection
Наука та прогрес транспорту. Вюник Дншропетровського нащонального ушверситету залiзничного транспорту, 2018, N° 1 (73)
На рис. 3, а, б представлены матрицы потенциального территориального риска в районе химической атаки для момента времени 200 с, при возможных метеоситуациях и рассматриваемых сценариях (с защитой и без защиты перед объектом атаки). Вероятность поражения людей возле здания показана в процентах на каждом рисунке. Принимается, что поражение наступает, если концентрация хлора в расчетной точке превышает 3 мг/м3.
б — b
Как видно из представленных рисунков, значительному риску (до 99 %) подвергается пространство возле объекта атаки, если отсутствует защита. При наличии защиты - риск поражения людей возле объекта атаки снижается. Время расчета составляет 5 сек.
Научная новизна и практическая значимость
Разработана численная модель, позволяющая оценить величину потенциального территориального риска в случае теракта с применением химического агента.
Особенностью построенной модели является использование стандартной исходной информации, быстрота в получении прогнозных данных и удобство для анализа получаемых результатов.
Выводы
В работе представлен метод оценки территориального риска при атаке террориста с применением химически опасного вещества и рассеиванием химического агента в условиях застройки. В основу метода положен численный расчет загрязнения воздушной среды в условиях застройки с последующей оценкой размеров зон поражения. Дальнейшее совершенствование выбранного научного направления следует проводить в области создания 3D-модели для расчета территориального риска в случае химической атаки террориста.
а — а
Рис. 3. Матрица территориального риска при химической атаке террориста возле объекта атаки:
а - без защиты, б - при защите
Fig. 3. Matrix of territorial risk in the chemical attack of a terrorist:
a - no protection near the target object, b - with protection
СПИСОК ИСПОЛЬЗОВАННЫХ ИСТОЧНИКОВ
1. Алымов, В. Т. Техногенный риск: Анализ и оценка : учеб. пособие для вузов / В. Т. Алымов, Н. П. Тарасова. - Москва : Академкнига, 2004. - 118 с.
2. Беляев, Н. Н. Защита зданий от проникновения в них опасных веществ : монография / Н. Н. Беляев, Е. Ю. Гунько, Н. В. Росточило. - Днепропетровск : Акцент ПП, 2014. - 136 с.
3. Оценка техногенного риска при эмиссии опасных веществ на железнодорожном транспорте / Н. Н. Беляев, Е. Ю. Гунько, П. С. Кириченко, Л. Я. Мунтян. - Кривой Рог : Р. А. Козлов, 2017. - 127 с.
4. Стоецкий, В. Ф. Оценка риска при авариях техногенного характера / В. Ф. Стоецкий, В. И. Голинько, Л. В. Дранишников // Наук. вюн. НГУ. - 2014. - № 3. - С. 117-124.
Наука та прогрес транспорту. Вюник Дншропетровського нащонального ушверситету залiзничного транспорту, 2018, № 1 (73)
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 / О. 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. Produced water treatment by ceramic membrane: A numerical investigation by computational fluid dynamics / Horte'ncia Luma Fernandes Magalha~es, Antonio Gilson Barbosa de Lima, Severino Rodrigues de Farias Neto, Helton Gomes Alves and Josedite Saraiva de Souza // Advances in Mechanical Engineering. - 2017. -Vol. 9. - Iss. 3. - Р. 1-20. doi: 10.1177/1687814016688642.
11. Regucki, P. Numerical modelling of sulphate ion concentration in wastewater from a closed cooling system / Pawel Regucki, Barbara Janowska // E3S Web of Conferences. - 2017. - Vol. 17. - Р. 1-8. doi: 10.1051/e3sconf/20171700078.
12. Tashvigh, A. A. Soft computing method for modeling and optimization of air and water gap membrane distillation - a genetic programming approach / Akbar Asadi Tashvigh, Bahram Nasernejad // Desalination and Water Treatment. - 2017. - Vol. 76. - Р. 30-39. doi: 10.5004/dwt.2017.20696.
М. М. БШЯеВ1*, I. В. КАЛАШН1КОВ2*, В. А. КОЗАЧИНА3*
'*Каф. «Пдравлжа та водопостачання», Дтпропетровський нацюнальний унгверситет з^зничного транспорту iKcm акадетка В. Лазаряна, вул. Лазаряна, 2, Дшпро, Украша, 49010, тел. +38 (056) 273 15 09, ел. пошта [email protected], ORCID 0000-0002-1531-7882
2 Держ. тдприемство «Проектно-вишукувальний шститут залiзничного транспорту Украши «Укрзалiзничпроект», вул. Червоноармшська, 7, Харкв, Украша, 61052, тел. +38 (057) 724 41 25, ел. пошта [email protected],
ORCID 0000-0002-2814-380X
3 Каф. «Пдравлжа та водопостачання», Дтпропетровський нацюнальний ушверситет залiзничного транспорту iменi акадетка В. Лазаряна, вул. Лазаряна 2, Дтпро, Украша, 49010, тел. +38 (056) 273 15 09,
ел. пошта [email protected], ORCID 0000-0002-6894-5532
РОЗРАХУНОК ТЕРИТОР1АЛЬНОГО РИЗИКУ П1Д ЧАС ТЕРАКТУ: ЕКСПРЕС МОДЕЛЬ
Мета. Робота передбачае розробку методу оцшки терш^ального ризику в разi терористичного нападу з використанням хiмiчного агента. Методика. Для опису процесу розсшвання в атмосферi хiмiчного агента, викинутого у випадку теракту, використовуеться рiвняння масопереносу домшки в атмосферному повггрг Рiвняння враховуе швидшсть вирового потоку, атмосферну дифузш, iнтенсивнiсть емiсii х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чних агентiв, так i в разi екстрема-льних ситуацiй на хiмiчно небезпечних об'ектах i транспортi. Метод може бути реалiзований на комп'ютерах мало1' та середньо1' потужностi, що дозволяе широко використовувати його для виршення задач даного класу. Представленi результати обчислювального експерименту, що дозволяють оцiнити можливосп запропонова-ного методу оцiнки територiального ризику в разi терористичного нападу з використанням хiмiчного агента. Наукова новизна. Запропоновано ефективний метод оцшки територiального ризику в разi теракту зi застосуванням хiмiчно небезпечноi речовини. Метод може бути використаний для оцшки територiального ризику в умовах мiськоi забудови, що дозволяе отримувати адекватнi данi про можливi зони ураження. Метод засно-
Наука та прогрес транспорту. Вюник Дншропетровського нащонального ушверситету залiзничного транспорту, 2018, № 1 (73)
ваний на чисельному iнтегруваннi фундаментального рiвняння масопереносу, що виражае закон збереження маси в редкому середовищi. Практична значимкть. Запропонований метод оцiнки територiального ризику в разi терористичного нападу зi застосуванням хiмiчного агента може бути використаний для розрахунку зон ураження бшя адмiнiстративних будiвель, центрiв та шших соцiально значущих об'ектiв.
Ключовi слова: територiальний ризик; теракт; хiмiчне забруднення; чисельне моделювання; забруднення повiтряного середовища
M. M. BILIAIEV1*, I. V. KALASHNIKOV2*, V. A. KOZACHYNA3*
'*Dep. «Hydraulics and Water Supply», Dnipropetrovsk National University of Railway Transport named after Academician V. Lazaryan, Lazaryan St., 2, Dnipropetrovsk, Ukraine, 49010, tel. +38 (056) 273 15 09, e-mail [email protected], ORCID 0000-0002-1531-7882
2 State Enterprise «Design and Exploration Institute of Railway Transport of Ukraine «Ukrzaliznichproekt», Str. Red Army, 7, Kharkov, 61052, тел. +38 (057) 724-41-25, e-mail: [email protected], ORCID ID: 0000-0002-2814-380X
3 Dep. «Hydraulics and Water Supply», Dnipropetrovsk National University of Railway Transport named after Academician V. Lazaryan, Lazaryan St., 2, Dnipropetrovsk, Ukraine, 49010, tel. +38 (056) 273 15 09,
e-mail [email protected], ORCID 0000-0002-6894-5532
TERRITORIAL RISK ASSESMENT AFTER TERRORIST ACT: EXPRESS MODEL
Purpose. The paper involves the development of a method to assess the territorial risk in the event of a terrorist attack using a chemical agent. Methodology. To describe the process of chemical agent scattering in the atmosphere, ejected in the event of a terrorist attack, the equation of mass transfer of an impurity in atmospheric air is used. The equation takes into account the velocity of the wind flow, atmospheric diffusion, the intensity of chemical agent emission, the presence of buildings near the site of the emission of a chemically hazardous substance. For numerical integration of the modeling equation, a finite difference method is used. A feature of the developed numerical model is the possibility of assessing the territorial risk in the event of a terrorist attack under different weather conditions and the presence of buildings. Findings. A specialized numerical model and software package has been developed that can be used to assess the territorial risk, both in the case of terrorist attacks, with the use of chemical agents, and in case of extreme situations at chemically hazardous facilities and transport. The method can be implemented on small and medium-sized computers, which allows it to be widely used for solving the problems of the class under consideration. The results of a computational experiment are presented that allow estimating possibilities of the proposed method for assessing the territorial risk in the event of a terrorist attack using a chemical agent. Originality. An effective method of assessing the territorial risk in the event of a terrorist attack using a chemically hazardous substance is proposed. The method can be used to assess the territorial risk in an urban environment, which allows you to obtain adequate data on possible damage areas. The method is based on the 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 assessing the territorial risk in the event of a terrorist attack using 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
REFERENCES
1. Alymov, V. T., & Tarasova, N. P. (2004). Tekhnogennyy risk: Analiz i otsenka: Uchebebnoe posobie dlya vuzov. Moscow: Akademkniga. (in Russian).
2. Biliaiev, N. N., Gunko, E. Yu., & Rostochilo, N. V. (2014). Zashchita zdaniy ot proniknoveniya v nikh opasnykh veshchestv: Monografiya. Dnepropetrovsk: Aktsent PP. (in Russian).
3. Biliaiev, N. N., Gunko, E. Yu., Kirichenko, P. S., & Muntian, L. Y. (2017). Otsenka tekhnogennogo riska pri emissii opasnykh veshchestv na zheleznodorozhnom transporte. Krivoy 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., & Biliaiev, N. N. (1997). Chislennoe modelirovanie rasprostraneniya zagryazneniya v okruzhayushchey srede. Kyiv: Naukova dumka. (in Russian).
Наука та прогрес транспорту. Вюник Дншропетровського нащонального ушверситету залiзничного транспорту, 2018, № 1 (73)
6. Barret, A. M. (2009). Mathematical Modeling and Decision Analysis for Terrorism Defense: Assessing Chlorine Truck Attack Consequence and Countermeasure Cost Effectivness. (Dissertation) 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. (undated). Retrived from: http://www.era.europa.eu/Document-Register/Documents/Cefic%20guidance%20on%20risk% 20assessment.pdf. (In English)
10. Horte'ncia Luma Fernandes Magalha~es, Antonio Gilson Barbosa de Lima1, Severino Rodrigues de Farias Neto, Helton Gomes Alves, & Josedite Saraiva de Souza. (2017). Produced water treatment by ceramic membrane: A numerical investigation by computational fluid dynamics. Advances in Mechanical Engineering, 9(3), 1-20. doi: 10.1177/1687814016688642. (in English)
11. Regucki, P., & Janowska, B. (2017) Numerical modelling of sulphate ion concentration in wastewater from a closed cooling system. E3S Web of Conferences, 17, 1-8. doi: 10.1051/e3sconf/20171700078. (in English)
12. Tashvigh, A. A., & Nasernejad, B. (2017) Soft computing method for modeling and optimization of air and water gap membrane distillation - a genetic programming approach. Desalination and Water Treatment, 76, 30-39. doi:10.5004/dwt.2017. (In English)
Статья рекомендована к публикации д.физ.-мат.н., проф. C. А. Пичуговым (Украина);
д.т.н., проф. С. З. Полищуком (Украина)
Поступила в редколегию: 22.09.2017
Принята к печати: 10.01.2018
