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

Ларионов В.И.

(кандидат технических наук, доцент, Центр исследований экстремальных ситуаций, г. Москва),

Акатьев В.А.

(кандидат технических наук, профессор, Центр исследований экстремальных ситуаций, г. Москва),

Александров А.А.

(инженер МГТУ им. Н.Э. Баумана, г. Москва)

ЗОНИРОВАНИЕ ТЕРРИТОРИИ В РАЙОНЕ АВТОЗАПРАВОЧНОЙ СТАНЦИИ ПО РИСКУ ВЗРЫВОВ БЕНЗИНО-ВОЗДУШНОЙ СМЕСИ ПРИ ЗАПРАВКЕ ЕМКОСТЕЙ

«Большие дыхания» происходят при заполнении резервуара, когда вся бензино-воздушная смесь (БВС) из резервуара вытесняется через дыхательный клапан в атмосферу. Учитывая большую интенсивность выброса БВС и сравнительно малое время заправки, «большие дыхания» можно рассматривать как залповые выбросы БВС. Для оценки взрывоопасности АЗС необходимо знать объем выброса и концентрацию паров бензина в БВС.

Для оценки концентрации паров бензина в газовой полости резервуара достаточно оценить минимальные (в герметичном объеме) и максимальные (с учетом «малых дыханий») ее значения. Общее давление в замкнутом газовом объеме резервуара будет суммироваться из парциального давления воздуха (атмосферное давление 101,3 кПа) и парциального давления паров бензина (давление насыщения). Зная отношение парциального давления паров бензина к общему давлению, можно определить концентрацию паров бензина в БВС в газовой полости резервуара /1/. Данные о давлении насыщения приняты по Н.Б. Варгафтику /2/.

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

сп = 100 Рн 7 (Рн + РоХ % об- (1)

где рн - давление состояния насыщения при температуре хранения, кПа; ро - атмосферное давление воздуха, кПа. Давление насыщения бензина и концентрация паров бензина в герметичной емкости, полученные для разных температур его хранения, представлены в табл. 1.

В действительности, максимальная концентрация паров бензина в газовой полости

резервуара при хранении значительно выше значений, представленных в табл. 1. Это связано с испарением бензина в резервуаре и многочисленными «малыми дыханиями». Объем выброса БВС при «малом дыхании» составляет до 2% объема газовой полости. По этой причине концентрация (табл. 2) может возрасти в 2 раза.

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

Плотность паров бензина рп может быть определена по формуле

М и Т

и о

рп = V т , кг/м3’ (2)

и

где Ми = 97 - молярная масса паров бензина, кг/кмоль;

Уи = 22,4 - молярный объем паров бензина, м3/кмоль;

То = 273 - стандартная температура, К;

Т = (273+1) - температура хранения бензина, К;

X - температура хранения, 0 С.

При 300 С плотность паров бензина равна

Ми То 97 273 рп = Уи Т = 22,4 303 _ 3,9 кг/м3-

Масса паров бензина в газовом пространстве резервуара Мб будет пропорциональна

Таблица 1

Температура хранения бензина, 0С +30 (лето) +5 25 (зима)

Давление насыщения бензина, кПа 16,8 5,3 2

Концентрация паров бензина в БВС, % об. 15 5 2

Таблица 2

Температура наружного воздуха, 0С 30 5 -25

Концентрация паров бензина в БВС, % об. 20-35 5-15 2-12

объемной концентрации сп, объему газового пространства Уп и плотности рп , т. е.

М = р У С /100, кг, (3)

п ~п п п ’’ у-/

где Уг - 0,8У;

У - объем резервуара, м3.

Результаты расчетов потерь бензина от одного «большого дыхания» резервуара приведены в табл. 3.

Таблица 3

Емкость резервуара, м3 Лето Весна-осень Зима

20 12-21 3.3-10 1.5-8.9

40 24-42 6.6-20 3.0-18

Потери бензина Мп от «больших дыханий» можно представить в виде линейной зависимости вида

М = кУ,

п 1 ’

где У - объем резервуара, м3;

к1 - удельная масса БВС для 1-го сезона, кг /м3 (табл. 4).

Таблица 4

Наименование Лето Весна-осень Зима

Удельная масса БВС, кь кг/м3 1.05 0.50 0.45

Радиус зоныг загазованности. В открытом воздушном пространстве при неподвижной воздушной среде радиус зоны загазованности с концентрацией, соответствующей нижнему концентрационному пределу распространения пламени (НКПР), в зависимости от массы выброса паров бензина может быть определен по формуле /3/

Я = 3,2 К1/2 (р /с )08 (М /(р р ))033, м, (4)

нкпр ’ чгн нкпру 4 пч'пгн// ’ ’V/

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

Мп - масса поступивших при дыхании паров бензина, кг;

рп - плотность паров бензина, кг/м3;

рн - давление насыщенных паров бензина

при расчетной температуре, кПа;

снкпр - концентрация, соответствующая

НКПР, % об. (для паров бензина снкп = 0,75

% об.);

К - коэффициент (К = Т/3600);

Т - продолжительность поступления паров ЛВЖ в открытое пространство, с.

Для емкости 50 м3, заправка которой проходит 1800 с, получаем следующую величину радиуса распространения паров бензина при безветрии

Янкпр = 3,2 (1800/3600)12 (16,8/0,75)08 (50/4 * 16,7)033 = 25 м.

Дрейф облака БВС. Диаметр облака БВС равен Б = 2Я = 50 м.

нкпр нкпр

При ветре возможен дрейф облака на расстояние до 150 м с вероятностью 0,7.

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

Частота взрыгва БВС. Известны аварии на АЗС со взрывом БВС при заправках. Например, 24 января 2003 г. при взрыве БВС на АЗС в г. Бугульме полностью разрушены автомойка и магазин, а телесные повреждения получили 3 человека.

По данным статистики за последние 10 лет в г. Москве на АЗС произошла одна авария с возгоранием БВС при «большом дыхании» емкости. Принимая, что к взрыву приводит каждый второй случай возгорания облака БВС, определим частоту аварии с взрывом Н на АЗС:

Н = п-Р/(№Т) = 1-0,5/(500-10) = 10-4, год1, (5) где п - число аварий, связанных с «большими дыханиями», на АЗС в Москве за 10 лет (п = 1);

Р - вероятность аварии со взрывом при условии, что произошло возгорание облака (Р = 0,5);

N - число АЗС в Москве (К = 500);

Т - время наблюдения (Т = 10 лет).

При дефлаграционном режиме избыточное давление взрывной волны ЛР р Па, определяется по формуле /5/

ДР^с /Р. = 1,788/ - 0,65/4, (6)

Р. = Ре(Уг • С0)2 • (о-1)/ о, (1Л >1 а.кР = 0,73251),

где = К/3/Ё7Р0;

Уг - скорость видимого фронта пламени (в расчетах принята равной 200 м/с);

С0 - скорость звука в воздухе (С0 = 340 м/с); а - степень расширения продуктов сгорания (а = 7).

Если Ка<Ка кр , то параметры взрывной волны не зависят от координаты, и в расчетах используются те же формулы с подстановкой значения К, = К,

а а.кр

Импульс фазы сжатия I Па-с, вычисляется из выражения

1^ /I. = 0,1293/ + 0,0464/ 12Л -0,025/ ,

і. = тЧодРоГ^-^ ---1

Со -

1 - 0,4

V, --1

Со

-

. (7)

Время фазы сжатия находится из уравнения

^эф = 21<1ес / ДРс!еС , с. (8)

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

Функция распределения дрейфа облака (рис. 1, а) определяет вероятность события, что облако будет на расстоянии меньше заданной величины Ь

_ Б(Ь) = Р(Ь < Ь), (9)

где Ь - случайные значения расстояния по направлению ветра.

Учитывая, что Р’(Ь) = А(Ь), построим функцию плотности распределения вероятностей (рис. 1, б) дрейфа облака. Тогда вероятность появления облака БВС в пределах отрезка ЛЬ на расстоянии Ь от емкости АЗС при заданном направлении ветра определяется из выражения

Р(Ь) = А(Ь)-ЛЬ. (10)

Функция плотности распределения вероятностей повторяемости ветров для г. Москвы приведена на рис. 2. Тогда вероятность появления облака БВС в пределах сектора аУ определяется из выражения

Р(У) = ф(У)ёУ. (11)

Для проведения расчетов вокруг емкости АЗС формируется матрица. Каждая ячейка матрицы представляет площадку с размерами 5х5 м2. При оценке вероятностей и интенсивностей

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

Расчетная схема к определению вероятности появления поражающего фактора различной интенсивности в пределах элементарной площадки с координатами (х, у) приведена на рис. 3.

Персонал на АЗС и клиенты могут находиться как в зданиях, так и на открытой местности. Принимается, что в здании находится 20% людей. На АЗС чаще применяются здания из легкого каркаса со стенами из панелей типа «Сэндвич».

ЗаконыI поражения людей. Под законами поражения людей понимают зависимость вероятности поражения людей от интенсивности поражающего фактора. Законы поражения людей, размещенных в зданиях, получены на основе теоремы полной вероятности. В расчетах учитывается, что событие С (общие, безвозвратные, санитарные потери) может произойти при получении зданием одной из степеней повреждения, образующих полную группу несовместных событий. Расчеты проводятся по формуле /7/

Р(Ф) = ЁРк(Ф) • Р(С^/В;), (12)

1=1

где Р(Ф) - вероятность поражения людей от воздействия поражающего фактора Ф; РВ1(Ф) - вероятность наступления 1-й степени повреждения здания при заданном значении поражающего фактора (закон разрушения здания);

Р^^ ) - вероятность получения людьми _)-й степени поражения при условии, что наступила 1-я степень повреждения здания;

б)

Г(Ц

Рисунок 1. Графики функции распределения ¥(Ь) (а) и функция плотности распределения вероятностей /(Ь) (б) дрейфа облака

п - рассматриваемое число возможных состояний повреждения здания.

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

импульс (при т+ £ 3Т/8 ) или

избыточное давление (при т+> 3Т/8), (13) где т+ - время действия фазы сжатия;

Т - период собственных колебаний конструкции (для одноэтажных панельных зданий можно принять Т = 0,07 с).

Степени разрушений зданий с ограждающими конструкциями из панелей типа «Сэндвич» от воздействия поражающего фактора даны в табл. 5.

V, град

Рисунок 2. Плотность распределения вероятностей ] (V) повторяемости ветров для г. Москвы

Таблица 5

Степени разрушения слабая | средняя | сильная | полная

При воздействии удельного импульса, Па-с 100-700 | 700-1000 | 1000-1300 | 1300-1700

При воздействии избыточного давления, кПа 10-30 | 30-50 | 50-65 | 65-105

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

а=к-Ф , (14)

расч 5 4 7

где Фрасч - величина расчетной нагрузки, вызывающая разрушения, при которых эксплуатация сооружения невозможна; к - коэффициент, зависящий от степени разрушения сооружения, принимаемый равным: для слабой степени - 0,1; средней - 0,2; сильной - 0,3; полной - 0,4.

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

Рисунок 3. Расчетная схема к определению вероятности появления поражающего фактора различной интенсивности в пределах элементарной площадки с координатами (х, у): 1 - зона действия взрыва;

2 - направление смещения облака; 0 - расположение емкости АЗС. Ь - рассматриваемое значение перемещения облака; W - точка взрыва; Я - расстояние от точки взрыва до центра элементарной площадки; Дх и Ду - размеры площадки; У - угол между вертикальной осью и направлением ветра;

Я8 - радиус зоны воздействия взрыва

Таблица 6

Степени разрушения и параметры нормального закона, Па-с

слабая средняя сильная полная

м а М а М а М а

400 85 850 170 1150 255 1800 340

Вероятность получения зданием не менее определенной степени разрушения выражается формулой

°(Ф) = 71Р ^ ехр(-^ •и2 )ди, (15)

где Z = (Ф-М)/а;

Ф - рассматриваемое значение поражающего фактора - импульса или давления; и - переменная интегрирования.

Значения интеграла (15) могут быть определены по формуле разложения

°(Ф) =1 Ф(х) + 2-; х = -^2;

<р(х)=7П

х--

1

. +—. 2!

I х_

3!" 7

1

+ — -4!

. (16)

Эти формулы следует ограничить диапазоном Ф <Ф<Ф . Соответствующие граничные

Ш1п тах •' г

значения составляют для аргумента Z = |3| и вероятностей 0(-3) = 0,0014 и 0(3) = 0,9986.

После воздействия здание может быть в одном из Ш несовместных событий: оказаться целым (событие В0), получить слабую В1, среднюю В2, сильную В3 или полную В4 степени разрушения (рис. 4). Вероятности наступления определенных степеней разрушения зданий могут быть представлены в виде зависимости:

Рвп(ф)=Рдп(ф);

Рві(ф)=Рді(ф) - Рлі+і(ф);

Рвд(Ф) = Рлд(ф) - Рл2(ф); Рво(Ф) = Рло(Ф) -Рлд(Ф),

(17)

где РА1(Ф), РА1+1(Ф) - вероятности наступления не менее 1-й, (1+1)-й степеней разрушения зданий.

При этом учитывается теорема о полной группе событий

£ Рві(ф)=д,

(18)

где Ш - число рассматриваемых событий;

РВ1 - вероятность наступления /-й степени

разрушения зданий.

Г рафики законов разрушения приведены на рис. 4.

В соответствии с изложенным методическим подходом могут быть получены законы разрушения зданий на воздействие давления. Вероятности наступления _)-й степени поражения людей в зданиях, при условии получения зданиями 1-й степени поражения, приняты по данным работы /8/ (табл. 7).

Таблица 7

Структура потерь Вероятность поражения людей при степени разрушения зданий

слабая средняя сильная полная

Общие 0,05 0,2 0,6 1

Безвозвратные 0 0,05 0,4 0,95

Санитарные 0,05 0,15 0,2 0,05

Законы поражения людей на открытой местности при воздействии взрывной волны приведены на рис. 6 /7/.

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

Тогда математическое ожидание потерь ЛМ(К на площадке с координатами х, у, при известной интенсивности поражающего фактора, можно определить по формуле

ДМ(х,у) = Р[Ф(х,у)] - ¥(х,у) • Дх • Ду , (19)

где у(х, у) - плотность людей в пределах площадки;

Лх и Лу - размеры, площадки;

Р[Ф(х, у)] - вероятность поражения людей в точке с координатами х, у при воздействии поражающего фактора Ф. Ожидаемое число пострадавших в целом в районе АЗС можно определить, суммируя потери по площадкам

М(К) = ДР[Ф(х,у)] • ¥(х,у) • ёх • ёу , (20)

где 8 - область интегрирования.

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

2р Ьтах

М(К) = Д | |Р[Ф(х,у)] ¥(х,у) •

8 0 0

• С(Ь) ф(У) • ёЬ • ёУ • ёх • ёу, (21)

где п = 3,14;

А(Ь) и ф(У) - функции плотности распределения вероятностей соответственно дрейфа облака и повторяемости направлений ветра за год.

В зависимости от решаемой задачи определяются общие, безвозвратные или санитарные потери. В этом случае ожидаемое число пострадавших с 1-ми условиями размещения при _)-м режиме взрывного горения М^(К) определяется по формуле

Законы поражения людей в поврежденных зданиях АЗС, полученные по формуле (15), приведены на рис. 5.

М(К) = £^£а;Му(К), (22)

]=1 1=1

где а. - доля людей в 1-х условиях размещения

а)

б)

I, Па-с

Рисунок 4. Законы разрушения зданий АЗС на воздействие импульса: а) вероятность возникновения не менее определенных степеней разрушения зданий; б) вероятность возникновения определенных степеней разрушения зданий Степени разрушения: 1 - слабая; 2 - средняя; 3 - сильная; 4 - полная

I, Па-с

/, Па с ЛРФ, кПа

Рисунок 5. Законы поражения людей в зданиях АЗС: а) на воздействие импульса; б) на воздействие давления;

1 - общие потери; 2 - безвозвратные; 3 - санитарные

ДРф, кПа

Рисунок 6. Законы поражения людей на открытой местности при воздействии взрывной волны:

1 - общие потери; 2 - безвозвратные; 3 - санитарные

(в расчетах принято а1= 0,8 - на открытой местности; а2= 0,2 - в зданиях);

^ - вероятность _)-го режима взрывного горения (в расчетах принимается ^= ё2= 0,5). С учетом изложенного, вероятность поражения людей в районе АЗС можно вычислить по формуле

P = N Я I IW^)]^).

• f(L) j(V) • dL • dV • dx • dy , (23)

а индивидуальный риск гибели человека из уравнения

Я = N Я Я |Р[Ф(х,у)] *(х,у) •

^ Б 0 0

| • С(Ь) ф(У) • ёЬ • ёУ • ёх • ёу , год-1, (24)

где N - общее число людей на рассматриваемой площади (области интегрирования по площади) АЗС;

Н - частота взрывного горения, год 1 .

2р Lm

Рисунок 7. Схема «типовой» АЗС:

1 - заглубленные резервуары; 2 - операторная; 3 - топливораздаточные колонки; 4 - автомобили клиентов; 5 - площадка парковки;

6 -трансформаторная подстанция

Моделирование аварий по сценарию «взрыв при большом дыгхании».

В качестве типовой принята АЗС, в состав которой входят: 4 емкости по 30 м3 каждая; здание операторной; 4-8 колонок; 4-8 автомобилей клиентов (рис. 7).

Размеры АЗС в среднем составляют 50х50 м2. Общая численность работающих на АЗС 12-15 чел. (8 операторов, работающих по 2 чел. через трое суток на четвертые; 1 уборщик территории; 1-2 слесаря; 1-3 заправщика; руководитель АЗС). Итак, на территории АЗС может находиться 10-17 чел.

В целях прогнозирования рисков разработана компьютерная программа «RISK».

Ожидаемое число погибших на АЗС в случае аварии со взрывом облака БВС может составить 1-2 чел., а общее число пострадавших -3-4 чел.

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

При решении задачи учитывается, что взрывоопасные объекты находятся на расстояниях 4-12 м до источника воздействия. Поэтому необходимо принимать во внимание взаимодействие ударно-волновых полей. Выход из строя соседних машин со взрывом бензобаков примем равным 100 кПа. Взаимодействие ударно-волновых полей ориентировочно учтем по формулам отраженных ударных волн.

На рис. 8 показаны ударноволновые поля при моделировании эффекта «домино» на АЗС. Расчеты показывают, что в случае взрыва облака БВС и при наличии машин клиентов на заправке детонация от первичного взрыва может распространиться на взрывоопасные мишени, находящиеся на расстоянии 12 м от емкости и 6 м друг от друга. При вторичных взрывах могут быть повреждены 4 легковых

Рисунок 8. Ударно-волновое поле при аварии на АЗС с учетом эффекта «домино»: 1 - емкость: 2 - автомобиль клиента

Рисунок 9. Поле вероятностей общего поражения людей при аварии на АЗС с учетом эффекта «домино»:

1 - емкость: 2 - автомобиль клиента

автомобиля с 20 литрами бензина. В этом случае площадь поражающего потенциала ударноволнового поля увеличи - лась в три раза.

На рис. 9 показано поле вероятностей общего поражения людей при вторичном взрыве трех легковых автомашин, в баках которых было 15 кг топлива. Произошел эффект «домино», в результате все объекты на АЗС взорвались. Площадь, в пределах которой возможно поражение людей, значительно увеличилась.

Зонирование территории вблизи АЗС по величине риска гибели человека было проведено с использованием ГИС-технологии.

Расчеты для «типовой» АЗС показали, что риск гибели человека при аварии со взрывом облака БВС может составить Я = 7,4-6 1/год, что превышает допустимое значение, равное Я = 10-6 1/год.

При зонировании территории вблизи АЗС по риску вначале проводилась оценка индивидуального риска гибели человека для каждой площадки. Затем точки с одинаковыми значениями риска соединялись изолиниями. При картировании область рисков разделена на четыре диапазона (рис. 10):

1) пренебрежительно малый риск 10-8 - 10-7;

2) малый риск 10-7 - 10-6 ;

3) повышенный риск 10-6 - 10-5;

4) высокий риск более 10-5.

Для зонирования территории г. Москвы по риску был разработан ГИС-проект. Структура ГИС включает базы данных по составу и размещению АЗС, математическим моделям и отчетной документации.

База данных включает картографическую М 1:10000 и семантическую информацию в це-

Условные обозначения М 1: 200 000

— 10-8 - 10-7 - пренебрежительно малый риск

— 10-7 - 10-6 - малый

— 10-7 - 10-5 - повышенный

— более 10-5 - высокий

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

смесей на АЗС в районе проспекта Мира

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

Блок математических моделей разделен на четыре группы: законы воздействия; законы поражения; модели прогнозирования ущерба; модели количественной оценки риска.

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

фактора проводилась с учетом влияния возможного дрейфа облака БВС и направления ветра.

Отображение средствами ГИС зон с различной величиной риска дает интегральное представление об опасности и об устойчивости элементов риска к поражающим факторам.

Выводы

В соответствии с ГОСТ Р 12.3.о47-98 к зонам с повышенным риском относят территорию, в пределах которой вероятность гибели человека составляет более Ю-6 год1. Следовательно, зона вокруг АЗС является зоной с повышенным риском. Учитывая, что в городах эксплуатируется большое количество АЗС, такое зонирование является актуальным.

Список использованной литературы:

1. Глинка Н.Л. Общая химия / Под ред. А.И. Ермакова. - М.: Интеграл-Пресс, 2002. - 728 с.

2. Варгафтик Н.Б. Справочник по теплофизическим свойствам газов и жидкостей. - М.: Наука, 1972. - 720 с.

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

4. Методика оценки последствий аварий на пожаро-, взрывоопасных объектах / Бодриков О.В., Елохин А. Н., Рязанцев Б.В. и др. - М.: МЧС России, 1994.

5. Методики оценки последствий аварий на опасных производственных объектах: Сборник документов. Серия 27. Выпуск 2. Гельфанд Б.Е., Дорофеев С.Б., Сидоров В.И. и др. - ГУП «НТЦ «Промышленная безопасность», 2001. - 224 с.

6. Теоретические основы реагирования на чрезвычайные ситуации. Механика разрушения / Под ред. В.И. Ларионова. - М.: 1999. - 276 с.

7. Акатьев В.А., Ларионов В.И., Ткачев В.Д. и др. Оперативное прогнозирование инженерной обстановки в чрезвычайных ситуациях / Под ред. С.К. Шойгу. - М.: Папирус, 1998. - 176 с.

8. Котляревский В.А., Ларионов В.И., Сущев С.П., Акатьев В.А. и др. Аварии и катастрофы. Предупреждение и ликвидация последствий. Кн. 6 / Под ред. В.А. Котляревского. - М.: Изд. АСВ, 2003. - 403 с.