CC BY
115
С. Г. АЛЕКСЕЕВ, канд. хим. наук, доцент, начальник отдела Уральского института ГПС МЧС РФ, г. Екатеринбург, Россия А. С. АВДЕЕВ, начальник сектора ГУ Судебно-экспертное учреждение ФПС "Испытательная пожарная лаборатория по Пермскому краю", г. Пермь, Россия Н. М. БАРБИН, д-р техн. наук, канд. хим. наук, заведующий кафедрой Уральского института ГПС МЧС РФ, г. Екатеринбург, Россия
С. А. ТИМАШЕВ, д-р техн. наук, профессор, директор Научно-инженерного центра "Надежность и ресурс больших систем и машин" УрО РАН, г. Екатеринбург, Россия Е. С. ГУРЬЕВ, канд. техн. наук, доцент, заместитель директора Научно-инженерного центра "Надежность и ресурс больших систем и машин" УрО РАН, г. Екатеринбург, Россия
УДК 614.84:665.74
МЕТОДЫ ОЦЕНКИ ВЗРЫВОПОЖАРООПАСНОСТИ ТОПЛИВОВОЗДУШНЫХ СМЕСЕЙ НА ПРИМЕРЕ КЕРОСИНА МАРКИ РТ. III. СП 12.13130.2009*
Проведен расчет параметров взрыва топливовоздушной смеси по методике, предложенной в СП 12.13130.2009, на примере керосина РТ. Выявлены возможности и недостатки данной методики. Ключевые слова: взрыв; керосин; топливовоздушная смесь; избыточное давление.
В данной работе продолжено рассмотрение возможностей существующих методик по оценке взрыво-пожароопасности топливовоздушных смесей (далее — TBC) на примере авиационного керосина марки РТ. Объектом настоящего исследования явилась методика расчета, представленная в СП 12.13130.2009 [1] (далее — СП), который пришел на смену НПБ 105-2003 [2]. Отметим, что данный СП, разработанный сотрудниками ВНИИПО МЧС России, предназначен для подтверждения требований Технического регламента по пожарной безопасности [3]. Сфера применения СП [1] допускает возможность рассмотрения и нашей гипотетической аварийной ситуации, связанной с разливом керосина марки РТ**.
Интенсивность испарения керосина W (кг/(м2-с)) при 35 °С определена по формуле [1]:
W = 10~6л/М • Рн = 1,15 • 10~4, (1)
где М — молекулярная масса жидкости, кг/кмоль; М = 151,2 кг/кмоль для керосина РТ [4, 5]; Рн — давление насыщенных паров, кПа; для керосина Т-1*** Рн = 9,33 кПа при 35 °С [6]. Расчетное значение интенсивности испарения авиакеросина марки РТ (0,58 ■ 10-4 кг/(м2 с)), полу-
* Продолжение. Начало см. в журнале "Пожаровзрывобез-
опасность", № 5 за 2010 г. и № 1 за 2011 г. * * Условия гипотетической аварийной ситуации и свойства авиационного керосина марки РТ приведены в работах [4, 5].
*** физико-химические свойства авиакеросинов Т-1 и РТ близки между собой [7], поэтому при расчете использовались данные для керосина Т-1.
ченное по формуле (1) для температуры 20 °С, хорошо согласуется с экспериментальными данными (0,57 ■ 10-4 кг/(м2-с)) при скорости воздушного потока 0,35-0,40 м/с. Эксперимент был проведен в вытяжном шкафу с использованием весов ВЛК-500 (рис. 1). Таким образом, правомочность применения уравнения (1) нашла независимое подтверждение.
Массу керосина, перешедшего в паровоздушную смесь за максимальное расчетное время испарения жидкости (Т = 3600 с), тп (кг) находим по уравнению
тп = WFиТ = 1,15 • 10~4 • 100 • 3600 = 41,4. (2)
Отметим, что для нашей расчетной аварийной ситуации значение массы паров керосина, полученное по уравнению (2), на порядок отличается от ре-
■ll'Stlr ' * _ •
__
Рис. 1. Общий вид установки для определения интенсивности испарения керосина марки РТ
© Алексеев С. Г., Авдеев А. С., Барбин Н. М., Тимашев С. А., Гурьев Е. С., 2012
зультата, который дает метод, предложенный в РД 03-409-01 [5] (435,44 кг), а с результатом, полученным по методике РБ Г-05-039-96 [4] (961,25 кг), расходится более чем в 23 раза.
Вычисляем приведенную массу паров авиакеросина тпр (кг) по формуле
Таблица 1. Избыточное давление и импульс волны давления взрыва TBC в зависимости от расстояния
m
пр
= (ßcr/ Qo) m п ^ = 39,3,
(3)
где Qcr — удельная теплота сгорания газа или пара, Дж/кг; Qcr = 42,9 ■ 106 Дж/кг [4]; Q0 — константа, Дж/кг; Q0 = 4,52 ■ 106 Дж/кг [1]; Z — коэффициент участия горючих газов и паров в горении; Z =0,1 [1]. С помощью уравнений (4) и (5) определяем зависимости избыточного давления АР и импульса волны давления i, развиваемые при сгорании ТВС, от расстояния геометрического центра топливоздуш-ного облака r:
АР = P
0,8тпр Vr+3тпр 6/r 2 + 5m Jr3 i = ш^р66 /r.
(4)
(5)
Результаты расчетов представлены в табл. 1.
Анализ полученных результатов показывает, что в радиусе 3мот геометрического центра облака ТВС образуется избыточное давление взрыва (вспышки) более 1 МПа. Это ставит под сомнение достоверность результатов расчета, так как известно, что для углеводородно-воздушных смесей избыточное давление взрыва в помещении обычно не превышает 0,6-0,8 МПа [8]. По В. Маршаллу, максимальное избыточное давление взрыва парового облака составляет 0,1-0,2 МПа, а уровень избыточного давления при детонации — порядка 1,8 МПа* [9].
Принимая во внимание, что существующие различные методы расчета давления взрыва дают порой несопоставимые результаты [11], для оценки так называемой нерабочей области ("мертвой зоны") уравнений (4) и (5) следует пользоваться "родственными" уравнениями. Наиболее подходящей формулой для этого является уравнение (6) из СП [1]. Данная формула используется для нахождения категорий помещений по взрывоопасности исходя из избыточного давления АР, которое определяется как разность между давлением взрыва Рвзр и атмосферным давлением Р0 для рассматриваемой гипотетической аварийной ситуации. При этом СП [1] не накладывает ограничений на объем помещения. Следовательно, между большим или очень большим помещением можно поставить знак приближенного равенства с открытым пространством. Хотя в СП [1] не рассматривается геометрия формы облака паро-
* Данное замечание совпадает с данными параметров детонации стехиометрических смесей ТВС низших углеводородов в верхней точке Чепмена - Жуге [10].
r, м АР, кПа i, Пас r, м АР, кПа i, Пас r, м АР, кПа i, Пас
1 23535,6 1387,4 16 35,2 86,7 80 4,0 17,3
2 3470,9 693,7 17 31,8 81,6 85 3,7 16,3
3 1205,3 462,5 18 29,0 77,1 90 3,5 15,4
4 591,6 346,9 19 26,6 73,0 100 3,1 13,9
5 349,8 277,5 20 24,6 69,4 105 2,9 13,2
6 232,0 231,2 23,4 19,4 59,3 110 2,8 12,6
7 166,4 198,2 25 17,6 55,5 115 2,6 12,1
8 126,1 173,4 30 13,6 46,2 120 2,5 11,6
9 99,6 154,2 35 11,0 39,6 125 2,4 11,1
9,24 94,6 150,2 40 9,2 34,7 130 2,3 10,7
10 81,2 138,7 45 7,9 30,8 135 2,2 10,3
11 67,8 126,1 50 7,0 27,7 140 2,1 9,9
12 57,8 115,6 55 6,2 25,2 145 2,0 9,6
12,88 51,0 107,7 60 5,6 23,1 150 2,0 9,2
13 50,1 106,7 65 5,1 21,3 151 2,0 9,2
14 44,1 99,1 70 4,6 19,8 152 1,9 9,1
15 39,2 92,5 75 4,3 18,5 155 1,9 9,0
воздушной смеси, но, принимая во внимание очень близкую аналогию между СП [1] и ГОСТ Р 12.3.047-98 [12] и приведенный в этой статье анализ по этому вопросу, можно сделать вывод, что форма паровоздушной смеси не изменяется при переходе от помещения большого объема к открытому пространству. Таким образом, выявленные аналогии свидетельствуют о возможности применения уравнения (6) для нахождения давления взрыва на границе паровоздушного облака ТВС и оценки "мертвой зоны" уравнений (4) и (5):
тп бсгР02
АР = Рвзр - Р0 =
Усв р вС TK н
(6)
где Р0 — атмосферное давление, кПа; Р0 =101,3 кПа [4, 5]; г =0,3 [1];
Усв — свободный объем помещения, м3; рв — плотность воздуха при 35 °С, кг/м3; Рв= 1,147 кг/м3;
Рв = М, /у (1 + 0,00367(10 - 273,15)]; Мв — молекулярная масса воздуха, кг/кмоль; Мв = 29 кг/кмоль;
Ур — киломольный объем газа при Т0 и Р0, м3/кмоль; Ур = 25,29 м3/кмоль*; Т0 — начальная температура воздуха, К; Т0 = 308,15 К [4, 5];
Ср — теплоемкость воздуха, Дж/(кг-К); Ср = 1,01103 Дж/(кг К) [1];
1 Киломольный объем Vtp найден с помощью универсального газового закона.
Кн — коэффициент, учитывающий негерметичность помещения и неадиабатичность процесса горения; Кн = 3 [1].
Учтем, что отношение тп/¥св в сущности представляет собой массовую концентрацию горючих паров в облаке ТВС. Для взрыва (вспышки) облака ТВС авиакеросина необходимо, чтобы концентрация его горючего пара находилась в области воспламенения, а для нашей задачи можно ограничиться двумя граничными условиями, когда концентрация горючих паров керосина равна НКПВ (СНКПВ = 0,7 % [4, 5]) и стехиометрической концентрации (Сстх = = 1,3 % [4, 5]). Для перевода объемной концентрации горючего вещества в облаке ТВС в массовую С' (кг/м3) воспользуемся уравнением
C = СМ/ (100 Vp). (7)
Таким образом, при Сстх = 7,77 ■ 10-2 кг/м3 максимальное избыточное давление взрыва (вспышки) облака ТВС составит 94,6 кПа, а при СНКПВ = = 4,18 ■ 10-2 кг/м3 — 51,0 кПа. Следовательно, "мертвая зона" уравнений (4) и (5) находится в диапазоне 9,23-12,87 мот геометрического центра облака ТВС.
СП [1] предлагает уравнение для определения радиуса горизонтальной взрывопожарной зоны с концентрацией горючего газа или пара выше НКПВ:
R нкпв = 3,1501 VK [-р-
V C НКПВ
0,813/ \0,333
т„
Р п Рн
(8)
где ЯНКПВ — радиус горизонтальной взрывопожар-ной зоны с концентрацией горючего газа или пара выше НКПВ, м; ^НКПВ = 23,4 м; К = ¿/3600 [1];
? — время испарения авиакеросина, с; для нашей задачи ? = 3600 с;
рп — плотность паров авиакеросина марки РТ при 35 °С, кг/м3; рп = 5,978 кг/м3; Рп = М/[^(1 + 0,00367(То - 273,15)]. Следует отметить, что СП [1] не предполагает приведения НКПВ (СНКПВ) к заданной температуре, в то время как ГОСТ 12.1.044-89 [13] позволяет это сделать с помощью следующей формулы:
Тп - 298,15"
СНКПВ^ = СНКПВ
1-
1250
= 0,694 %. (9)
Как видно из уравнения (9), переход к температуре 308,15 К не сильно отражается на изменении НКПВ и, как следствие, на радиусе RHKnBt (23,6 м), который увеличивается незначительно. Поэтому в нашем случае приведением НКПВ к заданным температурным условиям можно пренебречь.
Если даже допустить, что в радиусе от геометрического центра облака ТВС (RHKnB) концентрация горючих паров керосина равна НКПВ, то по логике вещей получается, что избыточное давление
взрыва (вспышки) на границе этой зоны должно составлять 51,0 кПа. Однако при расчете по формуле (4) для этого расстояния получается всего лишь 19,4 кПа избыточного давления (см. табл. 1). Таким образом, налицо несогласованность уравнений (4) и (8).
В СП [1] в отличие от методики расчета пожарного риска [14] и ГОСТ Р 12.3.047-98 [12] не дано определение высоты зоны паров ZНКПВ, в которой их концентрация не ниже НКПВ. Сравнение уравнения (10) [12, 14] с уравнением (8) показывает, что первое является упрощенным вариантом второго, что позволяет использовать формулу (11а), которая является уточненным вариантом уравнения (11) [12, 14], для нахождения ZНКПВ:
R НКПВ = 3,2^К
С
НКПВ
-НКПВ
= 0,12л/К
н
С
^нкпв = 0,12л/К
С
НКПВ
НКПВ 0,813
0,8
0,8
Р п Рн
0,33
(10)
т
п
Р п Рн
0,333
0,33
; (11)
Р п Рн
= 0,9 м. (11а)
Отметим, что в нормативных документах [1, 14] вопрос о форме облака ТВС не рассматривается со-всем,авГОСТР 12.3.047-98 [12] облако газо-и паровоздушных смесей представляется в виде цилиндра. Зная массу паров тп и геометрические размеры* (RHKnB = 23,4 м, ZHKnB = 0,9 м), можно найти среднюю концентрацию горючих паров керосина в цилиндрическом облаке ТВС, которая составит 2,68 ■ 10-2кг/м3. Учитывая, что это значение в 1,56 раза меньше НКПВ (4,18 ■ 10-2 кг/м3), можно еще раз констатировать, что формулы (8) и (11а) непригодны для оценки взрывоопасности ТВС. Если принять во внимание данные табл. 1 (RHKnB = 12,88 м), то высота цилиндрического облака ТВС составит 1,9 м.
Заметим, что в РБ Г-05-039-96 [15] приземное облако ТВС рассматривается как полусфера. Для расчета геометрии облака ТВС руководство предлагает применять следующие формулы:
V
ТВС
2240 х тк Т Q MCстх • 273,15 ;
(12)
V
2240 т п Т 0
ТВС
МСнкпв • 273,15
= 988,5 м3; (12а) = 0,78 #ТВС = 7,77 м, (13)
где х — доля массы керосина, переходящего в облако топливовоздушной смеси;
тк
масса исходного керосина.
* В нашем случае увеличением ЯНкпВ за счет радиуса пролива авиакеросина можно пренебречь.
Если учесть, что произведение %шк — это масса паров керосина тп, и вместо стехиометрической концентрации подставить значение НКПВ, то получится уравнение (12а), которое позволяет определить объем и радиус полусферы облака ТВС со средней концентрацией паров керосина, равной НКПВ. Несмотря на то что полученный объем УТВС с учетом погрешности формулы (12а) практически совпадает со значением объема паровоздушной смеси, определенным из соотношения тп/СНКПВ= 990,4 м3, полученное значение радиуса полусферы ТВС (7,77 м) меньше значения радиуса (12,88 м), приведенного в табл. 1, поэтому полусферическое облако ТВС можно исключить.
В отличие от ранее рассмотренных методик расчета пожаровзрывоопасности ТВС [4, 5] СП предлагает алгоритм определения интенсивности теплового излучения д по формуле
где ¥д
д = Еу¥д т, (14)
угловой коэффициент облученности;
т — коэффициент пропускания атмосферы.
Значение среднеповерхностной плотности теплового излучения пламени Еу принимается на основе имеющихся экспериментальных данных или по табл. 2. В связи с тем что в ней отсутствуют интересующие нас данные и что керосин является промежуточной фракцией между бензином и дизельным топливом при ректификации нефти, можно допустить, что значения Еу и удельной массовой скорости выгорания топлива Мвыг лежат между показателями бензина и дизельного топлива. На основании этого допущения получается, что для керосина Мвыг = 0,05 кг/(м2-с), что совпадает с литера-
Таблица 2. Среднеповерхностная плотность теплового излучения пламени в зависимости от диаметра пролива (очага) Ь и удельная массовая скорость выгорания Мвыг для некоторых жидких углеводородных топлив [1]
Топливо Еу, кВт/м2, при й, м кг/(м2-с)
10 20 30 40 50
СПГ (метан) 220 180 150 130 120 0,08
СУГ (пропан-бутан) 80 63 50 43 40 0,1
Бензин 60 47 35 28 25 0,06
Дизельное топливо 40 32 25 21 18 0,04
Нефть 25 19 15 12 10 0,04
Примечания: 1. СПГ — сжиженный природный газ; СУГ — сжиженные углеводороды. 2. Для диаметров очага менее 10 или более 50 м следует принимать значения Еу такими же, как и для очагов диаметром соответственно 10 и 50 м.
турными данными для удельной массовой скорости выгорания керосина (0,0483 кг/(м2-с) [16]).
Таким образом, показатель Еу может быть определен как среднее значение между данными бензина и дизельного топлива. Для эффективного диаметра пролива й = 11,287 м, вычисленного по формуле (15), среднеповерхностная плотность теплового излучения пламени Еу составит 48,65 кВт/м2. С помощью уравнения (16) находим высоту пламени Н (м):
Н = 42й
й = = 11,287 м;
Мвыг/(Р в4Ёй)
(15)
0,61
= 5,012 м, (16)
где ¥ — площадь пролива, м2; ¥ = 100 м2 [4, 5]; Я — ускорение свободного падения, м/с2; Я = 9,81 м/с2.
Угловой коэффициент облученности ¥д определяем с помощью формул:
¥д =
4
¥у ="
— аг^
И
Н ;
(17)
А
л/А2^
1;
(
агС^
И | агс1ё
5 -1 5 + 1
(А + 1)(5 - 1)
(А - 1)(5 + 1)
51 = 2г/й; И = 2Н/й; А = (И2 + 52 + 1)/(25);
(18)
(19)
(20) (21)
¥н ="
(В - 15)
.71
аг^
Л2 - 1 (А - 15)
(В + 1)(5 - 1) (В - 1)(5 + 1)
7^
агС^
1
(А + 1)(5 - 1) (А - 1)(5 + 1)
В = (1 + 52)/(25).
(22)
(23)
Коэффициент пропускания атмосферы т определяем по формуле
т = ехр[-7,0 ■ 10-4(г- 0,5й)]. (24)
Результаты расчетов интенсивности теплового излучения, углового коэффициента облученности и коэффициента пропускания атмосферы в зависимости от расстояния г представлены в табл. 3.
В СП [1] вводится такой параметр, как радиус воздействия высокотемпературных продуктов сгорания газо- или паровоздушной смеси в открытом пространстве Я¥, который равняется 1,2^НКПВ. В нашем случае Я¥ формально равен 31,2 м, но, принимая во внимание, что фактически ^НКПВ составляет не 23,4 м, а 12,88 м, принимаем Я¥ = 15,46 м.
1
Таблица 3. Расчетные данные интенсивности теплового излучения, углового коэффициента облученности и коэффициента пропускания атмосферы
r, м q, кВт/м2 Fq т
6 22,02 0,4527 0,9998
7 13,12 0,2700 0,9991
8 9,09 0,1872 0,9984
9 6,62 0,1364 0,9977
10 4,95 0,1021 0,9970
10,78 4,00 0,0826 0,9964
15 1,50 0,0310 0,9935
20 0,62 0,0129 0,9900
25 0,31 0,0065 0,9865
30 0,18 0,0037 0,9831
35 0,11 0,0023 0,9797
40 0,07 0,0015 0,9762
45 0,05 0,0011 0,9728
50 0,04 0,0008 0,9694
55 0,03 0,0006 0,9660
60 0,02 0,0004 0,9627
65 0,02 0,0003 0,9593
70 0,01 0,0003 0,9559
75 0,01 0,0002 0,9526
100 0,00 9E-05 0,9361
В СП [1] приведена заимствованная из НПБ 105-2003 [2] пробит-функция Pr1:
Pr1 = 5-0,26 lnV1, (25)
где V1 = (1750/ДР)8'4 + (290//)9'3.
Ранее нами было показано [5], что связь вероятностей событий с пробит-функциями описывается следующей формулой:
Р, = 0,4039(Pr,)5 - 10,092(Pr,)4 + 95,266(Pr,)3 -
- 420,16(Pr,)2 + 879,75(Pr,) - 707,94 (26) (r2 = 0,9999).
На основе решения уравнений (25) и (26) определена зависимость вероятности поражения человека избыточным давлением при сгорании ТВС керосина Р1 от расстояния r (рис. 2). При нахождении этой зависимости была исключена минимальная "мертвая зона" уравнений (4) и (5). Результаты расчетов вероятности поражения людей под действием избыточного давления по СП, представленные на рис. 2, вступают в противоречие с литературными данными [4, 8, 9, 17-19], что ставит под сомнение правомочность применения формулы (25).
Следует также отметить, что уравнение (25), взятое из РД 03-409-01 [20], в отличие от НПБ 105-2003 [2] и СП [1] используется не для нахождения услов-
9 10 И 12 13 14 г,м
Рис. 2. Вероятность поражения человека избыточным давлением Р1 при взрыве (вспышке) керосиновоздушной смеси от расстояния r
ной вероятности поражения человека избыточным давлением при сгорании газо-, паро- и пылевоздуш-ных смесей на расстоянии r от эпицентра, а для определения вероятности повреждений стен промышленных зданий, при которых возможно восстановление их без сноса*.
В СП [1] предлагается пробит-функция Pr2 для определения условной вероятности поражения человека тепловым излучением при пожаре пролива горючей жидкости, пожаре твердого материала или "огненном шаре" (Р2):
Pr2 = -12,8 + 2,56 ln (tq1'33); (27)
t = t0 + x/u, (28)
где t — эффективное время экспозиции, с;
q — интенсивность теплового излучения (см. табл. 3);
to — характерное время обнаружения пожара, с; допускается принимать to = 5 с [1]; x — расстояние от места расположения человека до зоны, где интенсивность теплового излучения не превышает 4 кВт/м2, м; x = 10,78 м (см. табл. 3);
u — скорость движения человека, м/с; допускается принимать u = 5 м/с [1]. Для нашей задачи условная вероятность поражения человека тепловым излучением, найденная по уравнениям (26)-(28), равна нулю. Полученный результат вступает в противоречие с данными расчетов по нахождению радиуса воздействия высокотемпературных продуктов сгорания газо- или паровоздушной смеси в открытом пространстве (RF) (см. выше). Заметим, что в ГОСТ P 12.3.047-98 также используется уравнение (27), однако первое слагаемое этого уравнения равно не минус 12,8, а минус 14,9 [12].
В заключение отметим, что СП [1] может применяться при прогнозировании взрывопожароопас-ности аварийных ситуаций, связанных с образованием ТВС. В качестве недостатка данной методики следует отметить плохую согласованность ее различных разделов между собой.
* Пример определения Prt по РД 03-409-01 дан в работе [5].
СПИСОК ЛИТЕРАТУРЫ
1. СП 12.13130.2009. Определение категорий помещений, зданий и наружных установок по взры-вопожарной и пожарной опасности (в ред. изм. № 1) : утв. приказом МЧС России от 09.12.2010 № 643 : введ. 01.05.2009. [Электронный ресурс]. Доступ из справочно-правовой системы "Кон-сультантПлюс".
2. НПБ 105-2003. Определение категорий помещений, зданий и наружных установок по взрыво-пожарной и пожарной опасности : утв. приказом МЧС России от 18.06.2003 № 314. [Электронный ресурс]. Доступ из справочно-правовой системы "КонсультантПлюс"
3. Технический регламент о требованиях пожарной безопасности : Федер. закон от 22.07.2008 № 123-ФЗ : принятГД РФ 04.07.2008 ;одобр. Советом Федерации 11.07.2008. [Электронныйресурс]. Доступ из справочно-правовой системы "КонсультантПлюс".
4. Алексеев С. Г., Авдеев А. С., Барбин Н. М. и др. Методы оценки взрывопожароопасности топливо-воздушных смесей на примере керосина марки РТ. I. РБ Г-05-039-96 // Пожаровзрывобезопас-ность. — 2010. — Т. 19, № 5. — С. 37-47.
5. Алексеев С. Г., Авдеев А. С., Барбин Н. М. и др. Методы оценки взрывопожароопасности топли-вовоздушных смесей на примере керосина марки РТ. II. РД 03-409-01 // Пожаровзрывобезопас-ность. — 2011. —Т. 20, № 1. —С. 21-27.
6. Варгафтик Н. Б. Справочник по теплофизическим свойствам газов и жидкостей. — М. : Наука, 1972.— С. 691.
7. Авиакеросин. [Электронный ресурс]. URL : Ьйр://га.'шМрей1а.о^/'шк1/Авиакеросин (дата обращения 15.08.2011).
8. Таубкин С. И. Пожар и взрыв, особенности их экспертизы. — М. : ВНИИПО, 1999. — 600 с.
9. Маршалл В. Основные опасности химических производств. — М. : Мир, 1989. — 672 с.
10. Бейкер У., Кокс П., Уэстайн П. и др. Взрывные явления. Оценка и последствия. — М. : Мир, 1986.—Кн. 1. —С. 36.
11. Акинин Н. И. Методические основы прогнозирования и предотвращения взрывов легковоспламеняющихся и взрывчатых материалов на опасных производственных объектах металлургических и коксохимических предприятий : дис. ... д-ра техн. наук. — М., 2005. — С. 133-146.
12. ГОСТ Р 12.3.047-98. ССБТ. Пожарная безопасность технологических процессов. Общие требования. Методы контроля. — Введ. 01.01.2000. — М. : Изд-во стандартов, 1998.
13. ГОСТ 12.1.044-89*. ССБТ. Пожаровзрывоопасность веществ и материалов. Номенклатура показателей и методы их определения. — Введ. 01.01.91. — М. : Изд-во стандартов, 1989; ИПК "Изд-во стандартов", 1996; 2001.
14. Об утверждении методики определения расчетных величин пожарного риска на производственных объектах : Приказ МЧС России от 10.07.2009 г. № 404 : зарег. в Минюсте РФ 17.08.2009 г., рег. № 14541. URL: http://www.mchs.gov.ru/law/index.php?ID=23378 (дата обращения 20.08.2011).
15. РБ Г-5-05-039-96. Руководство по анализу опасности аварийных взрывов и определению параметров их механического воздействия. — М.: НТЦ ЯРБ Госатомнадзора России, 2000. — 40 с.
16. Абдурагимов И. М., Говоров В. Ю., Макаров В. Е. Физико-химические основы развития и тушения пожаров. — М. : ВИПТШ МВД СССР, 1980. — С. 246.
17. Бейкер У., Кокс П., Уэстайн П. и др. Взрывные явления. Оценка и последствия. — М. : Мир, 1986.—Кн. 2.— С. 182.
18. Физика взрыва / Под ред. Л. П. Орленко. — М. : Физматлит, 2004. — Т. 1. — С. 608-612.
19. Комаров А. А. Прогнозирование нагрузок от аварийных дефлаграционных взрывов и оценка последствий их воздействия на здания и сооружения: дис.... д-ратехн. наук. —М., 2001. — С. 300.
20. РД 03-409-01. Методика оценки последствий аварийных взрывов топливно-воздушных смесей: утв. постановлением Госгортехнадзора РФ от 26.06.2001 № 25. [Электронный ресурс]. Доступ из справочно-правовой системы "КонсультантПлюс".
Материал поступил в редакцию 12 октября 2011 г.
Электронные адреса авторов: А1ех$коте@тш1.ги;
NMBarbin@yandex.ru; Тта$Неу8@швк1ги; sec@wekt.ru.
