автореферат диссертации по безопасности жизнедеятельности человека, 05.26.03, диссертация на тему:Модели оценки последствий распространения опасных веществ в воздушной среде

На правах рукописи

Сумской Сергей Иванович

МОДЕЛИ ОЦЕНКИ ПОСЛЕДСТВИЙ РАСПРОСТРАНЕНИЯ ОПАСНЫХ ВЕЩЕСТВ В ВОЗДУШНОЙ СРЕДЕ

Специальность 05 26 03 - Пожарная и промышленная безопасность (по проблемам нефти и газа)

АВТОРЕФЕРАТ диссертации на соискание ученой степени кандидата технических наук

Москва - 2007

003176052

Работа выполнена в Автономной некоммерческой организации «Агентство исследований промышленных рисков»

Научный руководитель: доктор технических наук Печеркин Андрей Станиславович

Официальные оппоненты:

доктор технических наук Буйновский Станислав Николаевич (ОАО «НТЦ «Промышленная безопасность») кандидат технических наук Ковалев Сергей Андреевич (ООО «ВНИИГАЗ»)

Ведущая организация

Российский государственный университет нефти и газа имени И М Губкина

Защита состоится 09 ноября 2007 в 14 часов 00 минут

на заседании диссертационного совета Д 417 001 01 при НТЦ «Промышленная безопасность» по адресу 105082, Москва, Переведеновский переулок, дом 13, строение 14, а/я 38, тел 620-47-50

С диссертацией можно ознакомиться в библиотеке НТЦ «Промышленная безопасность»

Автореферат разослан «_» октября 2007 года

диссертационного совета, д т н

Ученый секретарь

М В Лисанов

ОБЩАЯ ХАРАКТЕРИСТИКА РАБОТЫ Актуальность темы. Выброс опасных веществ (ОВ) в атмосферу является одним из самых серьезных вариантов аварии на опасном производственном объекте (ОПО) При выбросе ОВ в атмосферу может возникать целый спектр поражающих факторов токсическое поражение, образование в воздухе топливно-воздушных смесей (ТВС), снижение концентрации кислорода, генерация ударных волн (УВ), термическое поражение, химический ожег Размеры зон поражения могут достигать нескольких сотен и даже тысяч метров

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

Среди аварий и катастроф, происшедших за историю человечества и имеющих в своей основе химическую опасность, максимальное количество погибших наблюдалось именно при распространении ОВ через воздушную среду Авария, вызванная утечкой паров токсичного метилизоцианата в Бхопале (Индия) 3 декабря 1984 года, стала крупнейшей за всю историю мировой промышленности В ней погибло около 3 тысяч человек и пострадало около 200 тысяч человек В 1989 году на железной дороге под Уфой при аварии на продуктопроводе с широкой фракцией легких углеводородов произошел взрыв облака ТВС В результате аварии погибли и серьезно пострадали 1224 человека из 1284, оказавшихся в зоне поражения Это крупнейшая из всех катастроф, связанная с распространением в атмосфере пожаровзрывоопасного облака

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

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

2 Разработка моделей распространения и рассеяния ОВ в атмосфере

3 Исследование с помощью численного моделирования развития и последствий аварии при взрыве ТВС в тоннеле, где проходят трубопроводы с нефтью и нефтепродуктами

4 Исследование с помощью численного моделирования развития и последствий аварии на типовом объекте (в тч на типовом объекте нефтеперерабатывающих производств) - складе хлора

Научная новизна

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

2 Для описания рассеяния ОВ в атмосфере разработаны две модели

Первая модель описывает рассеяние в рамках интегрального метода и предназначена для расчета рассеяния выброса с плотностью больше плотности воздуха (т н «тяжелых газов») Модель отличается тем, что в ней используется переменная по пространству концентрация как в первичном, так и во вторичном облаках Так же отличительной особенностью является рассмотрение не единичных облаков, а цуга облаков, соответствующих различным стадиям аварии

Вторая универсальная модель рассматривает распространение и рассеяние в атмосфере неравновесной двухфазной среды («жидкие капли-газ») При этом учитывается как наведенная в ходе выброса, так и изначально присущая атмосфере турбулентность Отличие от других подходов состоит в том, что используется существенно измененная к-Е модель турбулентности, что позволяет более точно описать турбулентность в течениях с большими градиентами

3 Предложен методологический подход для анализа взрыва ТВС в тоннеле, предназначенном для транспортировки нефти и нефтепродуктов (по проложенным в тоннеле трубопроводам) Новизна подхода состоит в частности в том, что при оценке параметров генерируемых УВ, образующихся даже при детонации ТВС, предлагается учитывать влияние турбулентности Проведено численное исследование такого типа аварии Получены аналитические зависимости для различных нагрузок, действующих

на стенки тоннеля при детонации в нем TBC Исследованы параметры поля давления при выбросе продуктов детонации (ПД) из тоннеля Показано, что размер зоны поражения волнами давления при взрыве TBC в тоннеле существенно отличается от размера зон поражения при взрывах на открытых площадках Это отличие наблюдается как в форме зоны поражения, так и в характере спада давления с расстоянием Выяснено, что при взрывах в тоннелях большую опасность, нежели воздействие УВ, может представлять распространение в атмосфере горячих ПД

4 Предложен методологический подход для анализа аварии «полное разрушение контейнера со сжиженным хлором в закрытом помещении» Особенностью предложенного подхода является то, что учитываются такие факторы, как генерация волн давления в процессе вскипания хлора с учетом временной длительности этого процесса, забор хлора из помещения с помощью аварийной вентиляции, эмиссия хлора наружу из закрытого помещения, срабатывание водяных завес Проведено численное исследование такой аварии на примере контейнера с 1 т хлора По результатам численного исследования определено, что перепад давления во фронте хлорной волны при ее подходе к стенкам и потолочному перекрытию не превосходит 5 кПа что является условием сохранения целостности помещением склада после разрушения баллона В условиях сохранении целостности помещения склада и своевременного включения аварийной вентиляции путем численного моделирования показано, что размеры зон поражения при авариях рассматриваемого типа по сравнению с выбросом того же количества хлора на открытом пространстве значительно меньше (зона смертельного поражения в 5,8 раз, а зона порогового поражения в 3,5 раза)

Практическая ценность

Результаты работы могут быть использованы для прогнозирования последствий аварий, связанных с выбросом OB, например при разработке деклараций промышленной безопасности ОПО, мероприятий по защите персонала и населения от возможных аварий, локализации и ликвидации аварийных ситуаций, при оценке воздействия аварийных выбросов на окружающую природную среду, при проектировании объектов, на которых получаются, используются, перерабатываются, образуются, хранятся, транспортируются, уничтожаются опасные вещества, при обосновании условий страхования ответственности опасных объектов Основные положения работы реализованы в нормативных документах - Методика оценки последствий химических аварий (Методика "Токси" Вторая редакция) («ТОКСИ-2»),

- Методика оценки последствий химических аварий (Методика "Токси" Третья редакция) («ТОКСИ-3»)

- Правила безопасности для наземных складов жидкого аммиака

Методики реализованы в программном комплексе «ТОКСИ+»

Результаты работы нашли конкретное применение при оценки размеров опасных зон на таких объектах как «Сахалин-1», «Сахалин-2», перевалочные нефтебазы «Шесхарис» и «Грушовая», ряд РНУ «Транснефти», предприятие «Горводоканал» г Пензы и т д

Положения, выносящиеся на защиту В соответствии с перечисленными целями и задачами исследования на защиту выносятся следующие положения

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

2 Модель, описывающая рассеяние залповых и продолжительных выбросов OB в атмосфере в приближении «тяжелого газа» с учетом пространственного изменения концентрации в облаке

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

4 Методологический подход к анализу аварии в протяженных тоннелях, по которым осуществляется транспортировка нефти и нефтепродуктов, подход апробирован путем численного моделирования детонации облака TBC в тоннеле и последующего выброса из тоннеля продуктов детонации

5 Методологический подход к анализу аварии на типовом складе хлора, подход апробирован путем численного моделирования параметров хлорной волны и последствий воздействия на людей

Апробация работы. Материалы диссертации докладывались на Международной конференции студентов и молодых ученых "Современная физика и экология" (Москва 27 апреля - 1 мая, 1992, МИФИ), Научной сессии МИФИ-99 (Москва, 1999), 16th International Colloquium on the Dynamics of explosions and Reactive Systems (ICDERS) (Cracow, Poland, August 3-8, 1997), Второй международной конференции

"Безопасность трубопроводов" (28-31 августа, 1997), Третьей международной конференции "Безопасность трубопроводов" (6-10 сентября, 1999), 17th International Colloquium on the Dynamics of Explosions and Reactive Systems (ICDERS) (Heidelberg, 1999), научной сессии МИФИ (2000), VIII Всероссийском съезде по теоретической и прикладной механике (Пермь, 2001), Второй всероссийской конференции «Безопасность в нефтегазовой промышленности» (Москва, 2001), тематическом семинаре «Об опыте декларирования промышленной безопасности и развития методов оценки риска опасных производственных объектов» (Москава, 2002), Н-ой научно-практической конференции «Проблемы прогнозирования чрезвычайных ситуаций» (Москва, 23 октября 2002 г), Международной научно-практической конференции «Промышленная безопасность» (Москва, 2003 г), Международной научной школе «Моделирование и Анализ Безопасности и риска в сложных системах» (Санкт-Петербург, 2003, 2006), пятом тематическом семинаре «Об опыте декларирования промышленной безопасности и страхования ответственности» (Москва, 26-27 октября 2004), 6-ой научно-технической конференции, посвященной 75-летию Российского государственного университета нефти и газа им И М Губкина "Актуальные проблемы состояния и развития нефтегазового комплекса России" (Москва, 2005), тематическом семинаре «Об опыте декларирования промышленной безопасности Оценка риска аварий на опасных производственных объектах» (Москва, 2005), международной конференции «Управление рисками и устойчивое развитие единой системы газоснабжения России (Москва, 2006)

Публикации Основные результаты диссертации опубликованы в 27 печатных работах, список которых приведен в конце автореферата

Структура и объем диссертации Диссертационная работа состоит из введения, пяти глав и заключения Общий объем работы составляет 165 страниц, включая 6 таблиц, 40 рисунков, одно приложение и библиографический список литературы из 95 наименований

СОДЕРЖАНИЕ РАБОТЫ

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

Модель развития аварии строится для двух вариантов хранения OB в оборудовании

- в газообразном виде,

- в жидкой фазе (в т ч с наличием над поверхностью жидкости объема, заполненного газовой фазой)

При этом ОВ в жидкой фазе может находиться и в перегретом, и в захоложенном состоянии

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

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

Два для ОВ, находящегося в технологическом оборудовании в газообразном состоянии

сценарий 1 полное разрушение оборудования, содержащего ОВ в газообразном состоянии,

сценарий 2 нарушение герметичности (частичное разрушение) оборудования, содержащего ОВ в газообразном состоянии Для ОВ, находящегося в технологическом оборудовании в жидком состоянии сценарий 3 полное разрушение оборудования, содержащего ОВ в жидком состоянии,

сценарий 4 нарушение герметичности (частичное разрушение) оборудования, содержащего ОВ в жидком состоянии

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

Действительно, с точки зрения предельной общности можно выделить три типа оборудования

- емкостное оборудование (резервуары, сосуды под давлением, цистерны, реакторы, теплообменное оборудование и т д ),

- трубопроводы и установленный на них фитинг,

- насосно-компрессорное оборудование

Комбинация вышеприведенных типов оборудования позволяет смоделировать большую часть гипотетических аварий

При этом в общем случае можно выделить восемь возможных стадий развития аварии

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

- стадия 2 истечение жидкой фазы до отсечения аварийного участка,

- стадия 3 истечение жидкой фазы из аварийного участка после его отсечения,

- стадия 4 истечение газа при наличии пролива жидкой фазы и испарение с пролива,

- стадия 5 истечение газа из разрушенного оборудования при отсутствии пролива жидкой фазы,

- стадия 6 испарение ОВ с пролива при отсутствии истечения его жидкой или газообразной фаз из разрушенного оборудования,

- стадия 7 испарение ОВ из емкости при отсутствии пролива,

- стадия 8 ликвидация аварии (ликвидация отверстия разгерметизации и пролива)

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

выброса ОВ

На каждой стадии аварии формируются свои облака ОВ в атмосфере либо первичное (стадия 1), либо вторичные (стадии 2-7)

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

- размеры облака,

- содержание ОВ (масса для первичных облаков, расход для вторичных) в облаке с разбиением по фракциям,

- плотность среды в облаке,

- длительность поступления ОВ в окружающую среду и, в частности, в атмосферу

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

Модель рассеяния «тяжелого газа» в рамках интегрального подхода учитывает такие факторы, как

- изменение концентрации ОВ в облаке,

- движение облака в переменном по высоте ветре,

- гравитационное растекание,

- рассеяние облака в вертикальном направлении за счет атмосферной турбулентности с учетом наличия сред разной плотности,

- рассеяние облака в горизонтальном направлении (за счет атмосферной турбулентности и за счет гравитационного растекания),

- теплообмен за счет подмешивания воздуха,

- наличие аэрозолей и фазовые переходы в облаке,

- теплообмен облака с подстилающей поверхностью,

- частичное изменение скорости (по направлению)

Распределение концентрации опасного вещества в первичном облаке сй(х,у,г,1) описывается зависимостями (случай наличия в облаке «ядра» с постоянной концентрацией)

Г _ -,я "1 Г 222

(*-*„) +У -г

с°{х,у,г,г) = с^ ехр

2 р*

— - ехр

j

К2]

при (*-*„) +у'>Г

(1)

Основные уравнения для описания движения первичного облака в этом случае выглядят следующим образом

Сохранение массы выброшенного ОВ <2

2

б =С ц 71 ^ф Яэфф

Изменение массы облака <2сум1

^фф Рвозд 71 ^ф Рвозд Уподм ^"[^фф]

Гравитационное растекание облака

1-

Рэфф

Л1- -1 ""Мх1 Боковое рассеяние выброса за счет атмосферной диффузии

(2)

(3)

(4)

(5)

Сохранение энергии в облаке ^Эфф

¿А*]-« ^>ФФ еВ03И +

2 71 ^>фф Н^ Рво-щ Уподм ^ [**ф ] ево т+п ^Ф Е.

Положение центра облака л;

В уравнения (1-7) входят следующие величин!,! Д:,фф - эффективный радиус облака, Д>ФФ - эффективная высота облака, кэфф- эффективная скорость движения облака, Сц -концентрация О В в центре облака; X. — вертикальная дисперсия; —

определяющий скорость подмешивания через боковую поверхность, величина этого коэффициента выбрана равной 0,63 исходя из наилучшего совпадения расчетных и экспериментальных данных; — коэффициент, равный 1,151; — скорость

подмешивания воздуха в облако за счет диффузии в верти [сальном направлении; £„„„ — удельный тепловой поток в облако от поверхности земли; с — удельная внутренняя энергия подмешиваемого воздуха; ртт — плотпость воздуха в окружающей среде; — дисперсия вдоль оси у (в поперечном направлении).

Данная система при замыкании ее эмпирическими соотношениями для £П01, г;"'";,, а также при задании связи между плотностью р3(|,ф, температурой 7;фф и давлением в окружающей среде моле«' быть решена с помощью соответствующих численных методов.

Аналогичная модель построена и для продолжительного выброса. Обе модели были верифицированы но опытам Торней Айленд2. Между результатами расчета и экспериментальными данными имеет место удовлетворительное совпадение, что можно видеть на рис. 1, где для одного из опытов показана зависимость концентрации от расстояния.

горизонтальная дисперсия,

Уподм — коэффициент пропорциональности,

Гис. I, Зависимость ■ияк'ги.мдльио досгигасмой в Точке концентрации (ут расстояния при залповом ныПросс смет «ают-фрсон-П» (опыт 01)7 серии Торней Анлснд)

К.Р. Ьее5п [.«не РгеуепГюп Т.^' Ргосс^5 '' I: : I йачгоп, '^Л 1

МсО^а .1. оГНагаг[1оч^ Ммепа1. 16, р. 1—8, 1М7.

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

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

Чтобы при описании турбулентного движения учесть как изначально имеющуюся в атмосфере турбулентность, так и возникающую в процессе развития наведенных течений в турбулентной вязкости V, выделяются две составляющие и у„

(8)

У„ - турбулентная составляющая вязкости, обусловленная наведенным течением, У„ - турбулентная составляющая вязкости, обусловленная наличием турбулентности в атмосфере

Согласно ¿-^-модели турбулентная вязкость определяется соотношением

, (9)

е

где к и е - энергия турбулентности и скорость диссипации энергии турбулентности

Согласно стандартному подходу = 0 093 Однако такое значение коэффициента справедливо для условий равновесного потока, когда скорости диссипации и генерации турбулентности равны Использование Си = 0 09 в расчетах высокоскоростных потоков при наличии градиентов и УВ неоправданно завышает генерацию турбулентности

Исходя из вышесказанного предлагается использовать следующее выражение для

Сц

3 Селезнев В Е, Алешин В В, Клишин Г С, Фотин С В Численный анализ пожарной опасности магистральных газопроводов М Едиториал УРСС, 2004,328 с

В этой записи учтено взаимное влияние на коэффициент Си скоростей генерации (Р) и диссипации (г) турбулентной энергии, что значительно улучшает качество прогноза

Коэффициенты 0 7 и 6 8 в формуле (10) выбраны таким образом, чтобы, с одной стороны, совпадать со стандартной к-s моделью в условиях PIs = 1, когда Сц= 0 09, а с другой стороны соответствовать экспериментальным данным4, согласно которым в дальних следах струй, где генерация турбулентности невелика Р/е-0, Сц= 0 7

Для описания процессов кипения и испарения использовалось предположение о лимитирующей роли давления насыщения для мгновенного вскипания и температуры кипения для испарения

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

Численное моделирование аварий, связанных с взрывными превращениями в тоннеле. Одной из важнейших практических задач является оценка последствий аварии в тоннелях, по которым происходит транспортировка углеводородов, в т ч нефти и газа В данной работе рассматривается следующая задача В тоннеле длиной 3,4 км и диаметром 5 м проложены нефтепроводы В результате утечки нефти и последующего ее испарения во внутреннем пространстве тоннеля образуется облако TBC, которое инициируется на расстоянии 1,7 км от входа в тоннель При решении задачи рассматривалась ситуация, когда весь тоннель заполнен стехиометрической смесью паров нефти с воздухом и предполагался самый опасный режим энерговыделения - детонация После окончания детонации в тоннеле происходит струйный выброс ПД наружу, сопровождающийся распространением УВ

Задача о развитии аварии в тоннеле решалась в два этапа На первом этапе рассчитывалось течения до окончания детонации в тоннеле На втором - истечение ПД из тоннеля и формирование воздушной УВ

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

4 W Rodi, "The prediction of free turbulent boundary layers by use of a two-equation model of turbulence", Ph D Thesis University of London 1972

На основе этого решения записаны аналитические выражений для пространственно-временных зависимостей давления, скорости и температуры

Здесь л" - пространственная переменная, расстояние вдоль оси тоннеля, нулевая координата (0) соответствует середине тоннеля, месту инициирования детонации, г -пространственная координата, расстояние по радиусу, х. - положение фронта детонации, на момент выхода волны из тоннеля х.=Ы2, и - время, за которое детонация проходит расстояние х., £ - длина тоннеля, 1=3400 м, с1 — диаметр тоннеля, сЬ= 5 м, х, - точка разделяющая область покоя, присутствующую в тоннеле, от области движения ПД, их, и, - составляющие скорости ПД соответственно вдоль оси х и г, с - местная скорость звука в ПД, р,р,Т - плотность, давление и температура ПД, ц — молекулярный вес ПД, Л — универсальная газовая постоянная, у - показатель адиабаты ПД, £> - скорость детонации Индекс "О" в формулах (11) означает, что параметры относятся к фронту детонации, т е это параметры Чепмена-Жуге На рис 2 для двух моментов времени 1= 0,944 с (момент выхода детонации из тоннеля) и 0,472 с (время прохода детонационной волной Ы4) приведены распределения давления по длине тоннеля

Из полученных результатов следует, что поражающие факторы в тоннеле имеют чрезвычайно высокую интенсивность амплитуда давления на фронте волны составляет около 19 атм, температура повышается и держится продолжительное время (не менее 1-2 сек) на уровне 2500-2800 К, скорость движущихся по тоннелю ПД составляет несколько сот метров

(И)

На втором этапе для моделирования выброса ПД из тоннеля использовалась система уравнений Рейнольд са, что позволяет учесть влияние турбулентности па параметры генерируемой УВ.

-S

S

i 100 ■

Гис. 2 l':ni[p<'iv']fm'. я i;inictini[ пи длине ......с. i ti 1'и..-. 3 И sói.1 i очны'' даиления, достигаемые при

í.ivi i ну \ рамичк(.1\ момента иреиени 0,472 с. н ACI аннинн сферического облака его 0,944 с n-ф. i .i i р:Н]ни (2) н при детонации и тоннеле (3),

ч зависимое*и от расстояния до границы облаков.

Данные по достигаемым давлениям на оси тоннеля представлены на рис. 3.. На рисунке представлена величина избыточного давления в УВ в зависимости от расстояния. Для сравнения приведены давления в УВ при детонации и дефлаграции того же количества топлива. Вблизи от места выхода детонации из облака (на расстоянии в нескольких десятков мэтров) УВ при выходе из тоннеля спадает медленнее, однако затем этот спад становится быстрее, чем в УВ от сферического заряда.

Еще большее расхождение наблюдается при сравнении форм зон поражения давлением при детонации облака в тоннеле и на открытом пространстве. Гак при взрыве сферического облака зона с избыточным давлением 14 кПа (соответствует уровню разрушения оконных проемов и легко сбрасываемых конструкций) представляет собой круг радиусом около 200 метров, а при детонации облака в тоннеле - сектор с радиусом также 200 м, но полушириной 40 м.

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

Рассматривается авария в закрытом помещении (герметизированное здание без око л с плотными железными воротами) — разрушение контейнера с I т сжиженного

1.5

хлора При подобной аварии на открытой площадке зоны поражения при наихудших условиях составят 290 м - смертельное поражение и 830 м - пороговое воздействие При аварии в помещении рассматриваемого типа имеет место ряд особенностей, способных повлиять на размеры зоны поражения

В первую очередь была рассмотрена роль помещения С помощью численных расчетов показано, что помещение сохраняет целостность при повышении давления в результате аварии и не разрушается Этот вывод базируется на анализе характера ударно-волновых нагрузок на стены помещения Из рис 4, где приведена типовая зависимость давления от времени, видно, что повышение давления на фронте УВ не превосходит 5 кПа

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

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

Рис. 4 Типичная зависимость давления от Рис. 5 Распределение ингаляционной токсодозы

времени в волнах давления, воздействующих на (на оси выброса, на поверхности земли) в стены и потолок хлораторной зависимости от расстояния

В этом случае размер зоны смертельного поражения (на открытом пространстве) уменьшается до 50 м, а размер зоны пороговых воздействий до 240 м, что видно из рис 5 , где представлена зависимость токсодозы от расстояния на оси выброса При этом для людей, находящихся в помещении, зона пороговых воздействий составит 50 м

В заключении приводятся основные результаты и выводы по диссертационной работе

В приложении приводятся данные об использованных в работе коэффициентах турбулентной диффузии в атмосфере

Основные результаты диссертационной работы

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

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

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

4 Предложен методологический подход для анализа опасности типового объекта -склада хлора В подходе учитываются особенности распространения выброса хлора в помещении, наличие аварийной вентиляции, фильтрация наружу, наличие водяных завес Проведено численное моделирование распространение хлорного облака в замкнутом помещении и на открытом пространстве при эмиссии из помещения Показано, что своевременное включение аварийной вентиляции (с нейтрализацией забираемого хлора) приводит к сокращению в 3,5-5,8 раз размеров зон поражения по сравнению с ситуациями, когда выброс происходит на открытой площадке или при эмиссии из помещения, но с отключенной аварийной вентиляцией

5 Основные результаты работы реализованы в нормативных документах «Методика оценки последствий химических аварий (Методика "Токси" 2-я редакция», «Методика оценки последствий химических аварий (Методика "Токси" 3-я редакция», «Правила безопасности для наземных складов жидкого аммиака» Разработанные методики реализованы в программном комплексе «ТОКСИ+»

Публикации по теме диссертации

1 Башлачев В К , Григоров С И , Емельянов Е Н , Паталаха А Е , Сидоров В И , Сумской С И , Сидоров В И Печеркин А С , Шаргатов В А Расчет наземных концентраций сероводорода при возможных аварийных ситуациях на астраханском газовом комплексе Химическая промышленность 1994 №3 С 15-20

2 A A Borisov, SI Sumskoy. Modeling of accelerating gaseous flames in ducts Conference Proceedings 16-th International Colloquium on the Dynamics of explosions and Reactive Systems (ICDERS) Cracow University of Mining and Metalurgy 1997 P 2223

3 Башлачев В К , Григоров С И , Емельянов Е Н , Паталаха А Е , Сидоров В И , Сумской С И , Сидоров В И Печеркин А С, Шаргатов В А Оценка опасности производственных объектов, содержащих окисляющие вещества Химическая промышленность 1999 № 6 С 49-52

4 Губин С А , Маклашова И В , Лыков С М , Сидоров В И, Печеркин А С , Сумской С И Верификация методик для расчета рассеяния токсических выбросов Химическая промышленность 1999 № 10 С 58(662)-66(670)

5 Сумской С И , Губин С А , Маклашова И В Расчет рассеяния тяжелого газа в атмосфере при авариях на трубопроводах и емкостях // Научная сессия МИФИ-2000 Сб научных трудов Т5 М МИФИ 2000 С 46

6 Гражданкин А И, Дегтярев Д В , Лисанов М В, Печеркин А С, Сидоров В И , Сумской В И Анализ риска аварий на нефтепроводных системах КТК-Р и БТС Безопасность жизнедеятельности №6 2002 С 17-22

7 Дадонов Ю А , Гражданкин А И , Лисанов М В , Печеркин А С , Васильева Е В , Сумской С И Теория и практика оценки риска магистральных нефтепроводов В кн «Промышленная безопасность», сб избранных статей под общей редакцией В И Сидорова М ГУП «Научно-технический центр по безопасности в промышленности Госгортехнадзора России» 2002 С 276-283

8 Лыков С М, Гражданкин А И , Лисанов М В , Печеркин А С , Сумской С И , Ханин Е В Анализ риска газонаполнительной станции. В кн «Промышленная безопасность», сб избранных статей под общей редакцией В И Сидорова М Государственное унитарное предприятие «Научно-технический центр по безопасности в промышленности Госгортехнадзора России» 2002 С 300-317

9 Лисанов М В , Пчельников А В , Карманович А И, Сумской С И Оценка риска аварий на объектах хранения нефти и нефтепродуктов В кн «Об опыте декларирования промышленной безопасности и развития методов оценки риска опасных производственных объектов материалы тематического семинара» Под общ ред В И Сидорова М «Государственное унитарное предприятие «Научно технический центр по безопасности в промышленности Госгортехнадзора России» 2002 С 51-57

10 Лисанов MB, Пчельников А В., Сумской С И Моделирование и оценка последствий аварийных выбросов опасных веществ В кн «Об опыте декларирования промышленной безопасности и развития методов оценки риска опасных производственных объектов материалы тематического семинара» Под общ ред В И Сидорова М «Государственное унитарное предприятие «Научно

технический центр по безопасности в промышленности Госгортехнадзора России», 2002 С 51-57

11 Гражданкин А И , Дегтярев Д В , Лисанов М В , Печеркин А С , Сидоров В И , Сумской С И Анализ риска аварий и прогнозирование чрезвычайных ситуаций на нефтепроводных системах КТК-Р и БТС В кн Проблемы прогнозирования чрезвычайных ситуаций II научно-практическая конференция 23 октября 2002 г Доклады и выступления М Едиториал УРСС 2003 С 276-290

12 Дегтярев Д В , Лисанов М В , Печеркин А С , Пчельников А В , Сумской С И Совершенствование методического обеспечения оценки последствий аварий на опасных производственных объектах В кн Проблемы прогнозирования чрезвычайных ситуаций II научно-практическая конференция 23 октября 2002 г Доклады и выступления -М Едиториал УРСС 2003 С 317-323

13 Дадонов Ю А , Дегтярев Д В , Кручинина И А, Лисанов М В , Сумской С И Декларирование промышленной безопасности и анализ риска магистрального аммиакопровода «Тольятти-Одесса» ОАО «Трансаммиак» В кн «Международная научно-практическая конференция «Промышленная безопасность» тезисы докладов» Под общ ред В М Кульечева М «Государственное унитарное предприятие «Научно технический центр по безопасности в промышленности Госгортехнадзора России» 2003 С 109-112

14 Агапов А А , Бородовский С Я , Печеркин А С , Уткин Ю И , Сумской С И Программные средства для оценки последствий аварий с выбросом токсичных веществ В кн «Международная научно-практическая конференция «Промышленная безопасность» тезисы докладов» Под общ ред В М Кульечева М «Государственное унитарное предприятие «Научно технический центр по безопасности в промышленности Госгортехнадзора России» 2003 С 204-208

15 Лисанов М В , Пчельников А В , Сумской С И Оценка риска аварий на объектах хранения нефтепродуктов В кн Моделирование и Анализ Безопасности и риска в сложных системах Труды международной Научной Школы СПб Изд-во СПбГУАП 2003 С 403-409

16 Пчельников А В , Гражданкин А И , Кручинина И А , Дадонов Ю А , Лисанов М В , Сумской С И Оценка риска аварий на объектах хранения и перевалки нефти и нефтепродуктов Безопасность труда в промышленности №6 2004 С 33-37

17 Пчельников А В , Лисанов М В , Симакин В В , Кручинина И А , Сумской С И Количественная оценка риска аварий на объектах хранения нефтепродуктов Технологии ТЭК №4(17) 2004 С 72-75

18 Пчельников А В , Сумской СИ Об основных положениях методики оценки последствий выброса опасных веществ на основе модели «тяжелого газа» Об опыте декларирования промышленной безопасности и страхования ответственности Развитие методов оценки риска аварий на опасных производственных объектах Пятый тематический семинар Москва, 26-27 октября 2004 года Тезисы докладов М Федеральное государственное унитарное предприятие «Научно технический центр по безопасности в промышленности Госгортехнадзора России» 2004 С 106-110

19 Шаталов А А , Лисанов М В , Печеркин А С, Пчельников А В , Сумской С И Методика расчета распространения аварийных выбросов, основанная на модели «тяжелого» газа Безопасность труда в промышленности № 9 2004 С 46-52

20 Пчельников А В , Сумской С И, Лисанов М В Количественная оценка риска аварий на объектах перевалки нефти и нефтепродуктов В сб "Актуальные проблемы состояния и развития нефтегазового комплекса России" Тезисы докладов Том 2 Москва-2005 С 505

21 Пчельников А В , Лисанов М В , Печеркин А С , Шаргатов В А Сумской С И Верификация методик оценки последствий аварийных выбросов газа от источников продолжительного действия Безопасность труда в промышленности 2005 № 8 С 28-35

22 Пчельников А В , Лисанов М В, Сумской С И Моделирование рассеяния выбросов опасных веществ в атмосфере Российский химический журнал (Журнал Российского химического общества им Д И Менделеева) 2005 Т XLIX №4 стр 1828

23 Агапов А А , Буйко Е В , Гельфанд Б Е , Дадонов Ю А , Ильин А М , Карабанов Ю Ф , Лисанов М В , Сумской С И , Сидоров В И Печеркин А С , Шаталов А А Методики оценки последствий аварий на опасных производственных объектах (сборник документов) Серия 27 Декларирование промышленной безопасности и оценка риска Выпуск 2 М ГУП «Научно-технический центр по безопасности в промышленности Госгортехнадзора России», 2001 204 С

24 Лисанов М В , Гражданкин А И , Пчельников А В , Савина А В Сумской С И Анализ риска аварий на нефтепроводных системах БТС и МН «Дружба» Безопасность труда в промышленности 2006 №1 С 34-40

25 Агапов А А, Лисанов М В , Печеркин А С , Пчельников А В , Сумской С И , Шушунов В В Моделирование аварийных ситуаций на опасных производственных объектах Программный комплекс ТОКСИ+ (Версия 3 0) Серия 27 Декларирование промышленной безопасности и оценка риска Выпуск 5 М ООО «Научно-технический центр по безопасности в промышленности», 2006 252 С

26 Сумской С И, Пчельников А В , Шамшин И.О, Викторов С Б Моделирование взрывного превращения топливновоздушных смесей в протяженных тоннелях Безопасность труда в промышленности 2007 №3 С 46-50

27 Сумской С И , Пчельников А В , Шамшин И О, Викторов С Б Моделирование последствий взрывного превращения топливновоздушных смесей в протяженных тоннелях Безопасность труда в промышленности 2007 №4 С 42-48

Тираж 100 экз Заказ 234

Отпечатано с готовых оригинал-макетов на ротапринье

ОАО «НТЦ «Прочыш пенная безопасность»

105082 Москва, Переведеновский переулок, дом 13, строение 14,

Формат 60x90 1/16 Объем 1,0 п л

ВВЕДЕНИЕ.

ГЛАВА 1. МОДЕЛЬ РАЗВИТИЯ АВАРИЙНОЙ СИТУАЦИИ ПРИ

ВЫБРОСЕ ОПАСНОГО ВЕЩЕСТВА.

1.1. Основные события, инициирующие аварийную ситуацию.

1.2. Основные стадии развития аварийной ситуации.

1.3. Моделирование основных стадий развития аварий.

ГЛАВА 2. ИНТЕГРАЛЬНАЯ МОДЕЛЬ РАССЕЯНИЯ ВЫБРОСА ОПАСНОГО ВЕЩЕСТВА ДЛЯ СЛУЧАЯ ОТРИЦАТЕЛЬНОЙ ПЛАВУЧЕСТИ.

2.1. Модель рассеяния первичного облака (1-ая стадия аварийной ситуации).

2.2. Модель рассеяние вторичного облака (2-7-ые стадии аварийной ситуации).

2.3. Верификация модели.

ГЛАВА 3. МОДЕЛЬ РАСПРОСТРАНЕНИЯ ВЫБРОСА В РЕАЛЬНОЙ АТМОСФЕРЕ С УЧЕТОМ ВОЗНИКАЮЩИХ В ХОДЕ АВАРИИ ГРАДИЕНТОВ ДАВЛЕНИЯ, СКОРОСТИ И ПЛОТНОСТИ, А ТАКЖЕ С УЧЕТОМ НАЛИЧИЯ АЭРОЗОЛЕЙ И ИХ ФАЗОВЫХ ПЕРЕХОДОВ.

3.1. Уравнения движения.

3.2. Модель турбулентности.

3.3. Модель испарения и кипения капель.

3.4. Верификация модели.

ГЛАВА 4. МОДЕЛИРОВАНИЕ АВАРИЙ, СВЯЗАННЫХ СО

ВЗРЫВНЫМИ ПРЕВРАЩЕНИЯМИ В ТОННЕЛЯХ.

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

ЗАКРЫТЫХ СКЛАДАХ ХЛОРА.

5.1. Постановка задачи.

5.2. Моделирование распространение облака хлора на открытом ровном пространстве.

5.3. Особенности моделирование распространение хлора на рассматриваемых типовых объектах (складах хлора).

5.4. Моделирование распространение облака хлора в помещении (без учета работы вентиляции).

5.5. Моделирование распространение облака хлора при эмиссии из помещения (без учета работы вытяжной вентиляции).

5.6. Моделирование распространение облака хлора при эмиссии из помещения (с учетом работы вытяжной вентиляции).

5.7. Моделирование влияния водяных завес и защиты людей в помещении при распространение облака хлора.

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

При выбросе ОВ в атмосферу может возникать целый спектр поражающих факторов:

- токсическое поражение;

- образование в воздухе топливно-воздушных смесей (ТВ С), обладающих способностью к энерговыделению в различных режимах;

- снижение концентрации кислорода в атмосфере до уровня опасного для жизни;

- генерация ударных волн при расширении выбрасываемых в атмосферу веществ;

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

- химический ожег.

Приведенный перечень включает лишь химические опасности и может быть дополнен при введение в рассмотрение иных типов опасности, например, биологических или радиационных.

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

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

Бхопале (Индия) 3 декабря 1984 года, стала крупнейшей за всю историю мировой промышленности. В ней погибло около 3000 человек и пострадало около 200000 человек [1].

Катастрофы подобного масштаба могут возникать и в результате природных катаклизмов. Например, при выбросе углекислого газа на озере Ниос (Камерун) в августе 1986 года погибло около 1800 человек и пострадало около 1000 человек. Облако газа распространилось на расстояние до 25 км. Причиной гибели стало снижение концентрации кислорода до опасного уровня.

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

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

Выше были приведены два примера самых тяжелых последствий при токсическом поражении и при асфексии. Сопоставимую угрозу представляют и выбросы пожаровзрывоопасных веществ. Облака пожаровзрывоопасных веществ могут дрейфовать на сотни метров, сохраняя при этом способность к воспламенению. Воспламенение дрейфующих облаков вблизи мест сосредоточения людей может привести к многочисленным жертвам. Размеры зон, охватываемых пламенем (или детонационной волной), распространяющимся в ТВС, могут составлять сотни метров. Памятен такой случай. В 1989 году на железной дороге под Уфой при аварии на продуктопроводе с широкой фракцией легких углеводородов (ШФЛУ) произошел взрыв облака ТВС. Место утечки находилось на расстоянии около 900 м от железной дороги. В результате аварии погибли и серьезно пострадали 1224 человека из 1284, оказавшихся в зоне поражения [2]. Это крупнейшая из всех катастроф, связанная с распространением в атмосфере пожаровзрывоопасного облака.

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

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

К настоящему времени сложились три основных подхода для количественного описания процесса рассеяния выброса в атмосфере:

- гауссовские модели рассеяния, называемые также дисперсионными моделями;

- модели рассеяния, базирующиеся на интегральных законах сохранения либо в облаке в целом (залповый выброс), либо в поперечном сечении облака (продолжительный выброс); иногда такие модели называют моделями с сосредоточенными параметрами, одна из подгрупп этих моделей называется моделью рассеяния «тяжелого газа»;

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

Гауссовские модели, построенные на распределении Гаусса, описывают пространственно-временное изменение концентрации газа. Для практической реализации гауссовских моделей в них вводятся эмпирические коэффициенты, описывающие атмосферную турбулентность. Гауссовские модели учитывают только два процесса, происходящие с выбросом в атмосфере, — перемещение в поле ветра и рассеяние за счет атмосферной турбулентности. Естественно, что приемлемую точность оценок эти модели дают для тех зон, где доминируют указанные процессы: либо на большом удалении от места выброса, либо во всем диапазоне расстояний, но при слабом выбросе. В России для решения практических задач в области промышленной безопасности, долгое время широко использовались именно простые гауссовские модели. Примеры таких моделей многократно описаны в литературе [3,4, 5].

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

Большинство газовых выбросов при промышленных авариях в начальный момент имеют большую плотность, чем воздух. Причины этому могут быть различные:

- высокая молярная масса, более 29,5 г/моль (молярная масса воздуха), при нормальных условиях к «тяжелым газам» относятся пропан, бутан, пары нефти, бензина, хлор, сероводород и т.д.;

- низкая температура (например, метан при нормальных условиях имеет плотность ниже плотности воздуха, а при температуре, равной температуре кипения, его плотность становится выше плотности воздуха, и тогда метан уже является «тяжелым газом», поэтому процесс рассеяния метана от пролива следует рассчитывать с учетом этого факта);

- наличие аэрозолей (например, аммиак и при нормальных условиях, и при температуре кипения — «легкий газ», соответственно при смешении с воздухом образующиеся смеси имеют плотность меньше, чем воздух, однако в случае капельных включений в смеси она становится «тяжелым газом»);

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

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

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

Необходимость учета указанных эффектов с целью существенного улучшения прогностических возможностей методического аппарата (на сотни и даже тысячи процентов) и стала главной причиной массового перехода от моделей гауссовского типа к моделям рассеяния «тяжелого газа»: попытки модифицировать гауссовские модели к условиям рассеяния «тяжелых газов» в целом не дали приемлемой точности [6-8].

Для описания поведения в атмосфере «тяжелого газа» были разработаны специальные модели — «модели рассеяния тяжелого газа». Эти модели созданы в рамках т.н. интегрального подхода. Первоначально этот подход был разработан для турбулентных струй [9,10,11]. Суть интегрального метода состоит в том, что для описания течения используются интегральные характеристики. Для струй в качестве таких характеристик выбираются потоки массы, импульса, энергии и т.д. через поперечное сечение. В рамках интегрального метода определяются именно эти характеристики, а более детальная структура течения определяется по ним с помощью априорно заданных профилей.

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

Подробный обзор моделей «тяжелого газа» дан в [1,26]. Первая модель рассеяния «тяжелого газа» была предложена в 1974 г. ван Ульденом (van Ulden) [12]. В ней рассматривалась только одна интегральная характеристика - объем облака. Методика была предназначена для расчета только залпового выброса с однородным распределением и только в спокойной атмосфере.

Позднее эта модель неоднократно пересматривалась [13-18], однако ее ограниченность полностью преодолена не была.

Позднее был предложен еще ряд моделей: «Бритиш Газ» [19-22] (с постоянным профилем концентрации в облаке и постоянной скоростью перемещения облака), модели «Денз/Кранч» [23] (рассматриваются две интегральные характеристики объем и температура), «СЛЭБ» [24], «Хегадас/Хегабох» [25]. Последняя пара моделей нашла наибольшее распространение и фактически стала своеобразным «каноном» для моделей рассеяния «тяжелого газа». В ней рассматриваются облака-шлейфы от проливов с переменной концентрацией по объему и залповые выбросы с постоянной концентрацией по объему.

На сегодня существует и ряд компьютерных реализаций моделей рассеяния «тяжелого газа»: методика Всемирного банка [27], свод методик HGSYSTEM [28], методики, созданные при поддержке ведущих зарубежных исследовательских организаций TNO — The Netherlands Organization of Applied Scientific Research (Нидерланды) [29], Det Norske Veritas (DNV Technica) (Норвегия) [30].

В нашей стране значительный опыт в разработке моделей рассеяния, в т.ч. в приложении к «тяжелым газам», имеется в "НТЦ "Промышленная безопасность" [31-40]. Результаты этих работ нашли широкое приложение при анализе риска на целом ряде объектов [41-49], а также в ряде нормативных документов - «ТОКСИ 2» [50], «ТОКСИ 3» [53], «Правилах безопасности для наземных складов аммиака» [51] - и программном комплексе «ТОКСИ+»[52,53].

Помимо этого в России имелся единичный опыт разработок подобного рода - ГОСТ Р 12.3.047-98 [54]. В ГОСТ Р 12.3.047-98 дается методика расчета залповых выбросов на основе модели «тяжелого газа» в приложении к сжиженным углеводородным газам. К сожалению, конспективность изложения и ряд грубых опечаток делают фактически невозможным ее использование на практике. Для продолжительного выброса в [54] модель «тяжелого газа» вообще не используется.

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

Следует отметить, что количество разработок подобного рода (моделей «тяжелого газа») в мире сравнительно невелико и если оценивать их в целом, то можно отметить определенные недостатки: недостаточная верификация используемых в моделях коэффициентов; и, в некоторых случаях, недостаточная доведенность моделей до логической законченности. Например, для моделирования залповых выбросов практически все модели используют подход, при котором облако рассматривается однородным по всему объему с постоянной концентрацией.

Кроме того, в уже упомянутых работах [12-30] отсутствуют унифицированные модели развития аварийных ситуаций, привязанные к конкретному типу оборудования, а вместо этого рассматриваются отдельные физические процессы.

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

Как и любые модели, модель рассеяния «тяжелого газа» обладает рядом ограничений. Прежде всего, это связано с необходимостью учета пространственных факторов, в частности, особенностей рельефа местности и наличия застройки. В наибольшей степени учесть эти факторы позволяют модели, основанные на численном решении полной системы уравнений газодинамики. Это сложный, но наиболее перспективный подход. В настоящее время он реализован в нескольких программных продуктах, например, в вычислительных комплексах PHOENIX [55], AutoReaGas [56], в компьютерных программах, разработанных во ВНИИГАЗ [57,69], в "НТЦ "Промышленная безопасность" [70,71], в кодах «Сигмет» [58] и «ФЕМЗ» [59]. Однако не всегда, такие модели учитывают как реальные свойства атмосферы, так и сложные особенности наведенных газодинамических потоков. В частности, [57], [69] не позволяют рассматривать многофазные потоки, [55], [57] и [69] не учитывают особенности развития турбулентности в областях с сильными градиентами, в [56] не учитывается реальные характеристики атмосферы.

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

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

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

1) развитие и последствия аварийных ситуаций при взрывных превращениях топливно-воздушных смесей в протяженных тоннеля;

2) изучение влияния геометрии и формы оболочки гетерогенного заряда на характеристики взрывных волн, определение наиболее опасного сценария взрыва гетерогенного заряда;

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

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

2. Интегральная модель рассеяния залповых и продолжительных выбросов «тяжелого газа» с учетом переменной по пространству концентрацией ОВ и наличием аэрозольных включений.

3. Универсальная модель распространения и рассеяния выброса ОВ в атмосфере с учетом атмосферной и наведенной турбулентности, многофазности потока, фазовых переходов вещества. Универсальная модель включает в качестве составляющих частей модифицированную к-£ модель турбулентности (модификация касается учета влияния на турбулентность соотношения скорости ее генерации и диссипации) и модель фазовых переходов лимитируемых давлением и температурой в газовой фазе.

4. Методологический подход к анализу развития и последствий аварийных ситуаций при взрыве TBC в тоннеле.

5. Результаты исследования последствий аварийных ситуаций в тоннелях при взрыве в них TBC.

6. Методологический подход к анализу развития и последствий аварийных ситуаций при разрушении баллона с хлором на закрытом складе.

7. Результаты исследования последствий аварийных ситуаций на закрытых складах хлора.

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

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

Диссертация состоит из введения, пяти глав, заключения и приложения.

ЗАКЛЮЧЕНИЕ

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

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

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

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

4. Предложен методологический подход к анализу последствий взрывных превращений в тоннеле. Он состоит в рассмотрении двух стадий аварии: распространение взрывного превращения в тоннеле и выброс продуктов из тоннеля. Для первой стадии процесса на основе автомодельного решения задачи о распространении детонационной волны получены аналитические зависимости для задания силовых, скоростных и тепловых нагрузок на стенки тоннеля. Эти зависимости можно использовать для определения картины возможных разрушений самого тоннеля и оценки масштабов возможных обрушений и завалов в нем, а также для решения задачи об истечении продуктов детонации из тоннеля. Для оценки последствий взрывных процессов в тоннеле с точки зрения генерируемых в этом случае ударных волн предложен подход на основе решения уравнений Рейнольдса, а не Эйлера, обычно используемая при моделировании ударных волн. Показано, что зоны поражения ударной волной при авариях в тоннелях определяются не только и не столько количеством топлива, участвовавшим в том или ином взрывном превращении. Эти зоны заметно меньше, чем зоны поражения при сгорании (взрыве) того же количества топлива на открытом пространстве и имеют форму сектора, а не круга. Также показано, что в случае, когда быстро выделившаяся при взрыве энергия передается в окружающую среду длительное время, существенную роль может играть генерируемая в атмосфере турбулентность. Показано, что при оценке последствий взрывов в тоннелях может возникнуть необходимость оценки зон поражения не только от ударной волны, но и от других поражающих факторов, в первую очередь от распространения горячих продуктов взрыва. В решенной задаче о взрыве в тоннеле зона поражения человека ударной волной оказалась меньше, чем зона распространения горячих продуктов.

5. Предложен методологический подход для анализа опасности типового объекта склада хлора. В подходе учитываются особенности распространения выброса хлора в помещении. Также учитывается наличие аварийной вентиляции, фильтрации наружу через неплотности помещение, наличие водяных завес и нахождение людей в помещении. Показано, что волны давления, образующиеся при разрыве одного однотонного контейнера со сжиженным хлором в помещении объемом 2000 м3 не способны вызвать разрушение здания. Повышение давления в такой волне не превосходит 5 кПа. Показано, что при сохранении целостности здания склада и своевременном включение аварийной вентиляции с нейтрализацией хлора зона смертельного поражения сокращается до 50 м, а порогового воздействия до 240 м; при выбросе на открытом пространстве зона смертельного поражения составляет 290 м, а порогового 830 м; аналогичные зоны поражения наблюдаются и при выбросе в помещении при отключенной вентиляции. В случае если люди укрываются в помещении зона порогового воздействия поражения составит 50 м.

Результаты выполненных расчетов и их качественное и количественное согласие с имеющимися экспериментальными данными позволяют сделать общий вывод о том, что предложенные модели способны адекватно описывать распространение выброса в атмосфере, а разработанные на их основе методики «Токси-2» [5], «Токси 3» [60], «Правила для наземных складов аммиака»[51] и программный комплекс «ТОКСИ +» могут быть рекомендованы для использования в расчетах различных задач промышленной безопасности и экологии.

1. Маршалл В. Основные опасности химических производств. М.: Мир, 1989, 672 с.

2. Makhviladze G.M., Yakush S.E. Proc. of the 29th Int. Symp. on Comb., 2002, p. 195— 210.

3. Горский В.Г., Моткин Г.А., Петрунин B.A., Терещенко Г.Ф., Шаталова А.А., Швецова-Шиловская Т.Н. Научно-методические аспекты анализа аварийного риска М.: Экономика и информатика, 2002, - 260 с.

4. F. Pasquill, Atmosferic Diffusion, Ellis Horwood, Chichester, 2nd edn., 1974.

5. Burgess D., Zabetakis M.G. Fire and explosion hazards of LNG. In: "Fire Protection Manual for Hydrocarbon Processing Plants" (ed. Vervalin C.V.), Huston: Gulf, 1964, p. 49.

6. Clancey V.J. Liquid and vapour emission and dispersion. Course on Loss Prevention in the Process Industries. Dept. of Chem. Engng. Loughborourgh Univ. of Tecnol., 1976.

7. Germeles A.E, Drake E.M. Gravity spreading and atmospheric dispersion of LNG vapour clouds. Transport of Hazardous Cargoes, 4,1975.

8. Абрамович Г. H., Гиршович Т. А., Крашенников С. Ю., Секундов А. Н., Смирнова И. П. Теория турбулентных струй. М.: Наука, 1984, 714 с.

9. Бруяцкий Е.Н. Турбулентные стратифицированные течения. Киев: Наук. Думка, 1986,294 с.

10. Cox R.A., Roe D.R. A model of the dispersion of dense vapour clouds. Loss Prevention and Safety Promotion, 2, p. 359,1977.

11. Cox R.A., Carpenter R.J. Further development of a dense vapour cloud dispersion model for hazard analysis. In: Heavy Gas Dispersion and Risk Assessment. Dordrecht: Reidel, 1980, p. 55-73.

12. Bradley C.I., Carpenter R.J., Waite P.J., Ramsay C.G., English M.A. Recent development of a simple box-type model for dense vapour cloud dispersion. In: Heavy Gas Dispersion and Risk Assesment-2, Dordrecht: Reidel, 1980, p. 77-98.

13. Carpenter R.J., Cleaver R.P., Waite P.J. English M.A. The calibration of simple model for dense gas dispersion using the Thorney Island Phase 1 trials data. J. Haz. Materials, 16, 293-307,1987.

14. Kaiser G.D., Walker B.C. Releases of anhydrous ammonia from pressurized containers -the importance of denser-than-air mixtures. Atmos. Environment, 12, p. 2289-2304, 1978.

15. Ermak D.L., Chan ST. A study of heavy gas effect on atmospheric dispersion of dense gases. In: Air Pollution Modelling and its Application (ed. de Wispelaere C.), 5, p. 723, NY: Plenum Press, 1986

16. Puttock J.S. The development and use of HEGABOX/HEGADAS dispersion models for hazard analysys. In: Vapour Cloud Modeling, p. 317, 1987.

17. F.P. Lees, Loss Prevention in the Process Industries, 2nd edn., Butterworth-Heinemann, Boston, MA, 1996.

18. Руководство по оценке индустриальных опасностей (Techniques for Assessing Industrial Hazards: a Manual) // World Bank Tech. Paper 55, 1988.

19. The HGSYSTEM version 3.0 technical reference manual/ Shell Internationale Research Maatschappij BV. — Hague, 1994.

20. Methods for the calculation of physical effects CPR 14E. Part 1. TNO Yellow book. 3-rd edition. TNO, 1997.

21. Witlox H.W.M. Unified Dispersion Model (UDM Version 6.0) Consequence Modelling Documentation: Theory Manual — Det Norske Veritas, 2000.

22. Губин С.А., Маклашова И.В., Лыков С.М., Сидоров В.И., Печеркин A.C., Сумской С.И. Верификация методик для расчета рассеяния токсических выбросов. Химическая промышленность. 1999. № 10. С. 58(662)-66(670).

23. Сумской С.И., Губин С.А., Маклашова И.В. Расчет рассеяния тяжелого газа в атмосфере при авариях на трубопроводах и емкостях. // Научная сессия МИФИ-2000. Сб. научных трудов. Т.5. М.: МИФИ. 2000. С. 46.

24. Шаталов A.A., Лисанов М.В., Печеркин A.C., Пчельников A.B., Сумской С.И. Методика расчета распространения аварийных выбросов, основанная на модели «тяжелого» газа. Безопасность труда в промышленности. № 9.2004. С. 46-52.

25. Пчельников A.B., Лисанов М.В., Печеркин A.C., Шаргатов В.А. Сумской С.И. Верификация методик оценки последствий аварийных выбросов газа от источников продолжительного действия. Безопасность труда в промышленности. 2005. № 8. С. 28-35.

26. Пчельников A.B., Лисанов М.В., Сумской С.И. Моделирование рассеяния выбросов опасных веществ в атмосфере. Российский химический журнал (Журнал Российского химического общества им. Д.И.Менделеева). 2005. Т. XLIX. №4. стр. 18-28.

27. Гражданкин А.И., Дегтярев Д.В., Лисанов М.В., Печеркин A.C., Сидоров В.И., Сумской В.И. Анализ риска аварий на нефтепроводных системах КТК-Р и БТС. Безопасность жизнедеятельности. №6. 2002. С. 17-22.

28. Лисанов М.В., Пчельников A.B., Сумской С.И. Оценка риска аварий на объектах хранения нефтепродуктов. В кн.: Моделирование и Анализ Безопасности и риска в сложных системах: Труды международной Научной Школы. СПб.: Изд-во СПбГУ АП. 2003. С. 403-409.

29. Пчельников A.B., Гражданкин А.И., Кручинина И.А., Дадонов Ю.А., Лисанов М.В., Сумской С.И. Оценка риска аварий на объектах хранения и перевалки нефти и нефтепродуктов. Безопасность труда в промышленности. № 6. 2004. С. 33-37.

30. Пчельников A.B., Лисанов М.В., Симакин В.В., Кручинина И.А., Сумской С.И. Количественная оценка риска аварий на объектах хранения нефтепродуктов. Технологии ТЭК. № 4(17). 2004. С. 72-75.

31. Пчельников A.B., Сумской С.И., Лисанов М.В. Количественная оценка риска аварий на объектах перевалки нефти и нефтепродуктов. В сб.: "Актуальные проблемы состояния и развития нефтегазового комплекса России". Тезисы докладов. Том 2. Москва 2005. С. 505.

32. Методика оценки последствий химических аварий (Методика "Токси". Вторая редакция) («Методика ТОКСИ-2»)

33. ПБ 09-579-03 Правила безопасности для наземных складов жидкого аммиака.

34. М.: «Государственное унитарное предприятие «Научно технический центр по безопасности в промышленности Госгортехнадзора России». 2003. С. 204-208.

35. ГОСТ 12.3.047-98 ССБТ Пожарная безопасность технологических процессов. Общие требования. Методы контроля.

36. Иванов А.В., Мастрюков Б.С. Изв. вузов. Черная металлургия, 1999,1 11, с. 64—68.

37. Century Dynamics/TNO (1999). AutoReaGas User Documentation. Version 3.0.

38. Едигаров A.C., Сулейманов B.A. Математическое моделирование, 1995,1 4, с. 37— 52.

39. Ermak D.L., Chan S.T., Morgan D.J., Morris L.K., A comparison of dense gas dispersion model simulation with Burro series LNG spill tests results. J. Haz. Materials, 6 (1\2), p/ 129-143,1982.

40. Методика оценки последствий аварийных выбросов опасных веществ (Методика «Токси», редакция 3.1). Проект. М.:ФГУП НТЦ «Промышленная безопасность», 2005.-67 с.

41. В.Г. Мацак, Л.К. Хоцянов, "Гигиеническое значение скорости испарения и давления пара токсических веществ, применяемых в производстве", М.:Медгиз, 1959,231 с.

42. Общих правил взрывобезопасности для взрывопожароопасных химических, нефтехимических и нефтеперерабатывающих производств (ПБ 09-540-03)

43. Н.К. Fauske, М. Epstein "Source term consideration in connection with chemical accidents and vapor cloud modeling", J. Loss Prev. Process Ind., 1988, Vol. 1, April, pp. 75-83.

44. ESDU Data Sheet 74031, Characteristics of atmospheric turbulence near the ground -Part II. Single point data for strong winds (neutral atmosphere), Engineering Sciences Data Unit, London, 1974

45. McQuaid, J., "Some experiments on the structure of stably-stratified shear flows", Tech. Paper p. 21, Safety in Mines Research Establishment, Sheffield, UK, 1976

46. Kranenburg, C., "Wind-induced entrainment in a stably stratified fluid", J. Fluid Mech. 145, pp.253-273, 1984

47. McQuaid J. J. of Hazardous Material, 16, p. 1—8,1987.

48. Mercer A., Nussey C. J. of Hazardous Material, 16, p. 9—20, 1987.

49. Селезнев В.Е., Алешин В.В., Клишин Г.С., Фотин С.В. Численный анализ пожарной опасности магистральных газопроводов. М.: Едиториал УРСС, 2004, 328 с.

50. Сумской С.И., Пчельников А.В., Шамшин И.О., Викторов С.Б. Моделирование взрывного превращения топливновоздушных смесей в протяженных тоннелях Безопасность труда в промышленности. 2007. №3. С. 46-50

51. Сумской С.И., Пчельников А.В., Шамшин И.О., Викторов С.Б. Моделирование последствий взрывного превращения топливновоздушных смесей в протяженных тоннелях. Безопасность труда в промышленности. 2007. №4. С. 42-48.

52. Нигматулин Р.И. Динамика многофазных сред. 4.1. М.: Наука, 1987.

53. W. Rodi, "The prediction of free turbulent boundary layers by use of a two-equation model of turbulence", Ph. D. Thesis, University of London, 1972.

54. Шамшин И. О. Моделирование течений при взрывах многофазных сред. Диссертация на соискание ученой степени кандидата физико-математических наук. М. 2003.

55. R.B. Miles, J. Connors, Е. Markovitz, P. Howard, and G. Roth, "Instantaneous supersonic velocity profiles in an underexpanded sonic air jet by oxygen flow tagging", Phys. Fluids A, vol. 1,№2, p.389-393,1989

56. D. Pawel, P. J. Van Tiggelen, H. Vasatko, H. Gg. Wagner. Initiation of detonation in Various Gas Mixtures. Combustion & Flame. 15. p. 173-177, 1970.

57. Гостинцев Ю.А., Губин С. А., Шаргатов В. А. Численное моделирование процессов при горении открытых объемов перемешанных газовых смесей. Препринт. Черноголовка: ИХФ АН СССР, 1986.

58. G.W.Koroll, S.R.Mulpuru. Dilution with Steam on the Burning Velocity and Structure of Premixed Hydrogen Flames. Twenty-first Symposium (Int.) on Combustion , 1986, pp. 1811 -1819.

59. Fujimoto M. Chemical reaction in shock wave I: The ignition delay of a hydrogen-oxygen mixture in a shock tube. Bull Chem Soc Japan , 1963, vol. 36, p. 1233 1236.

60. Волков O.M., Проскуряков Г.А. Пожарная безопасность на предприятиях транспорта и хранения нефти и нефтепродуктов. М., Недра, 1981,256 с.

61. ПБ 09-540-03 Общие правила взрывобезопасности для взрывопожароопасных химических, нефтехимических и нефтеперерабатывающих производств

62. РД 03-409-01 Методика оценки последствий аварийных взрывов топливновоздушных смесей (Методика ТВС) (с изменениями и дополнениями).

63. Щелкин К.И, Трошин Я. К. Газодинамика горения. М.: Изд-во АН СССР, 1963. 255 с.

64. Зельдович Я. Б., Баренблатт Г. И., Либрович В. Б., Махвиладзе Г. М. Математическая теория горения и взрыва. М.: Наука, 1980.

65. Зверев И.Н., Смирнов Н.Н. Газодинамика горения. М.: Изд-во Моск. ун-та, 1987.

66. Викторов С.Б., Губин С.А., Маклашова И.В., Пепекин В.И. Прогнозирование детонационных характеристик безводородных взрывчатых веществ // Химическая физика. 2005. Т. 24. N. 12. С. 22-45.

67. Годунов С.К., Забродин А.В., Иванов М.Я., Крайко А.Н., Прокопов Г.П. Численное решение многомерных задач газовой динамики. М.: Наука, 1976.

68. Колган В.П. Применение принципа минимальных значений производной к построению конечно-разностной схемы для расчета разрывных решений газовой динамики // Ученые записки ЦАГИ. 1972. Т. 3. № 6. С. 68.

69. Техника безопасности при хранении, транспортировании и применении хлора.// А.Ф. Тимофеев, Б.Ю. Ягуд//М.,1996. 519 с.

70. Methods for the calculation of physical effects CPR 14E (Part 1, 2) "TNO Yellow book", 3rd edition, The Netherlands Organisation of Applied Scientific Research, The Netherlands, 1997

71. Methods for the determination of possible damage to people and objects resulting from releases of hazardous materials CPR 16E "TNO Green book", 1st edition, The Netherlands Organisation of Applied Scientific Research, The Netherlands, 1992

72. McRae G.J., Goodin W.R., Seinfield J.H. Atmospheric Environment, 16, 4, p.679-696, 1982.

73. Businger J.A., Wyngaard J.C., Isumi Y., Bradley E.F. J. Atmos. Sci., 28, 2, p. 181-189, 1971.

74. Doran J.C., Horst T.W., Nicola P.W. Validation in measured values of lateral diffusion parameter. J. Appl. Met., 17, p.825-831, 1978.