Список рекомендуемой литературы

1. Основы пожарного дела. Теребнев В.В., Артемьев Н.С., Шадрин К.В. – М.: Центр Пропаганды, 2006. – 328 с.

2. Справочник спасателя-пожарного. Теребнев В.В., Артемьев Н.С., Грачев В.А. – М .: Центр Пропаганды, 2006. – 528 с.

3. Система стандартов безопасности труда. ГОСТ Р 12.3.047 – 98. Пожарная безопасность технологических процессов. Общие требования. Методы контроля. – ВНИИПО, 1998.

4. Основные процессы и аппараты химической технологии. Касаткин А.Г. – М.:ООО ТИД «Альянс», 2005.– 753 с.

5. Химическая технология: Учеб. Пособие для студ. высш. учеб. Заведений: В 2т. – М.: Гуманит. изд. центр ВЛАДОС, 2003. – Т.1– 368 с.

6. Федеральный закон Российской Федерации от 22 июля 2008 г. № 123 – ФЗ «Технический регламент о требованиях пожарной безопасности»

.

Приложение 1

Метод расчета индивидуального и социального риска

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

Согласно расчетной методике уровень обеспечения безопасности людей при пожарах отвечает требуемому, если:

,                                 (1.1)

где  — нормируемый индивидуальный риск,  = 10^-6 год^-1; Q_в — расчетный индивидуальный риск.

Расчетный индивидуальный риск Q_в в каждом здании (помещении) рассчитывают по формуле

Q_в = Q_п P_пp (1 – Р_э) (1 – P_п.з).                                     (1.2)

где Q_п — вероятность пожара в здании в год; Р_пр — вероятность присутствия людей в здании, при работе: 0,33 — в одну смену; 0,67 — в две смены; 1,00 — в три смены; Р_э — вероятность эвакуации людей; Р_п.з — вероятность эффективной работы технических решений противопожарной защиты.

Вероятность эвакуации Р_э рассчитывают по формуле

Р_э = 1 - (1 - Р_э.п)(1 - Р_д.в).                   (1.3)

где Р_э.п — вероятность эвакуации по эвакуационным путям; Р_д.в — вероятность эвакуации по наружным эвакуационным лестницам, переходам в смежные секции здания

Вероятность Р_э.п рассчитывают по формуле

(1.4)

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

Расчетное время эвакуации людей из помещений и зданий устанавливают по расчету времени движения одного или нескольких людских потоков через эвакуационные выходы от наиболее удаленных мест размещения людей. При расчете весь путь движения людского потока подразделяют на участки (проход, коридор, дверной проем, лестничный марш, тамбур) длиной li и шириной di. Начальными участками являются проходы между рабочими местами, оборудованием, рядами кресел и т. п. При определении расчетного времени длину и ширину каждого участка пути эвакуации принимают по проекту. Длину пути по лестничным маршам, а также по пандусам измеряют по длине марша. Длину пути в дверном проеме принимают равной нулю. Проем, расположенный в стене толщиной более 0,7 м, а также тамбур следует считать самостоятельными участками горизонтального пути, имеющими конечную длину li. Расчетное время эвакуации людей tр следует определять как сумму времени движения людского потока по отдельным участкам пути ti по формуле

t_p = t₁ + t₂ + t₃ +,...- + t_i,                                      (1.5)

где t₁ — время движения людского потока на первом (начальном) участке, мин; t₁, t₂, t₃ ,..., t_i — время движения людского потока на каждом из следующих после первого участка пути, мин.

Время движения людского потока по первому участку пути ti, мин, рассчитывают по формуле

,                                 (1.6)

где l₁ — длина первого участка пути, м;  — скорость движения людского потока по горизонтальному пути на первом участке, м/мин (определяют по таблице 1.1 в зависимости от плотности D).

Плотность людского потока на первом участке пути D₁ рассчитывают по формуле

                                     (1.7)

где N₁ — число людей на первом участке, чел; f— средняя площадь горизонтальной проекции человека, м², принимаемая равной

0,100 — взрослого в домашней одежде;

0,125 — взрослого в зимней одежде;

0,070— подростка;  — ширина первого участка пути, м.

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

,                                  (1.8)

где ,  — ширина рассматриваемого i-го и предшествующего ему участка пути, м; ,  — интенсивности движения людского потока по рассматриваемому i-му и предшествующему участкам пути, м/мин [интенсивность движения людского потока на первом участке пути q = определяют по таблице 1.1 по значению D₁, установленному по формуле (1.7)].

Если значение  определяемое по формуле (1.8), меньше или равно q_max, то время движения по участку пути ti, мин, равно:

                                        (1.9)

при этом значения q_max, м/мин, следует принимать равными: 16,5 — для горизонтальных путей; 19,6 — для дверных проемов; 16,0 — для лестницы вниз; 11,0— для лестницы вверх.

Таблица 1.1 — Интенсивность и скорость движения людского потока при различной плотности на разных участках путей эвакуации

Плотность потока D, м2/м2	Горизонтальный путь		Дверной проем, интенсивность q, м/мин	Лестница вниз		Лестница вверх
	Скорость v, м/мин	Интенсивность q, м/мин		Скорость v, м/мин	Интенсивность q, м/мин	Скорость v, м/мин	Интенсивность q, м/мин
0,01	100	1,0	1,0	100	1,0	60	0,6
0,05	100	5,0	5,0	100	5,0	60	3,0
0,10	80	8,0	8,7	95	9,5	53	5,3
0,20	60	12,0	13,4	68	13,6	40	8,0
0,30	47	14,1	16,5	52	16,6	32	9,6
0,40	40	16,0	18,4	40	16,0	26	10,4
0,50	33	16,5	19,6	31	15,6	22	11,0
0,70	23	16,1	18,5	18	12,6	15	10,5
0,80	19	15,2	17,3	13	10,4	13	10,4
0,90 и более	15	13,5	8,5	8	7,2	11	9,9
Примечание — Интенсивность движения в дверном проеме при плотности потока 0,9 и более, равная 8,5 м/мин, установлена для дверного проема шириной 1,6 м и более, а при дверном проеме меньшей ширины d интенсивность движения следует определять по формуле q = 2,5 + 3,75 d

Если значение q_i, определенное по формуле (1.8), больше q_max то ширину S, данного участка пути следует увеличивать на такое значение, при котором соблюдается условие

q_i  > q_max                     (1.10)

При невозможности выполнения условия (1.10) интенсивность и скорость движения людского потока по участку i определяют по таблице 1.1 при значении D = 0,9 и более.

Приложение 2

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

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

а) происходит расчетная авария одного из аппаратов;

б) все содержимое аппарата поступает в помещение;

в) происходит одновременно утечка веществ из трубопроводов, питающих аппарат по прямому и обратному потоку в течение времени, необходимого для отключения трубопроводов.

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

- времени срабатывания системы автоматики отключения трубопроводов согласно паспортным данным установки, если вероятность отказа системы не превышает 0,000001 в год или обеспечено резервирование ее элементов (но не более 120 с);

- 120 с, если вероятность отказа системы автоматики превышает 0,000001 в год и не обеспечено резервирование ее элементов;

- 300 с при ручном отключении.

г) происходит испарение с поверхности разлившейся жидкости. Площадь испарения при разливе на пол определяют (при отсутствии справочных данных), исходя из расчета, что 1 л смесей и растворов, содержащих 70 % и менее (по массе) растворителей, разливается на площади 0,5 м², а остальных жидкостей — на 1 м² пола помещения;

д) происходит также испарение жидкостей из емкостей, эксплуатируемых с открытым зеркалом жидкости, и со свежеокрашенных поверхностей;

е) длительность испарения жидкости принимают, равной времени ее полного испарения, но не более 3600 с.

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

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

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

Избыточное давление D_р, кПа, для индивидуальных горючих веществ, состоящих из атомов С, Н, О, N, Cl, Вr, I, F, рассчитывают по формуле

,                               (2.1)

где p_max — максимальное давление, развиваемое при сгорании стехиометрической газовоздушной или паровоздушной смеси в замкнутом объеме. При отсутствии данных допускается принимать p_max  равным 900 кПа; р₀ — начальное давление, кПа (допускается принимать равный 101 кПа); т — масса горючего газа (ГГ) или паров легковоспламеняющихся (ЛВЖ) и горючих жидкостей (ГЖ), вышедших в результате расчетной аварии в помещение, кг; Z— коэффициент участия горючего при сгорании газопаровоздушной смеси. Допускается принимать Z по таблице 2.1; V_св — свободный объем помещения, м³; ρ _{г, п} — плотность газа или пара при расчетной температуре t_р, кг/м³;

В качестве расчетной температуры следует принимать максимально возможную температуру воздуха в данном помещении в соответствующей климатической зоне или максимально возможную температуру воздуха по технологическому регламенту с учетом возможного повышения температуры в аварийной ситуации. Если такого значения расчетной температуры tр определить не удается, допускается принимать ее равной 61 °С;

С_ст — стехиометрическая концентрация ГГ или паров ЛВЖ и ГЖ, % (об.), вычисляемая по формуле

,                                   (2.2)

где  — стехиометрический коэффициент кислорода в реакции сгорания; n_с, n_н, n_о, n_х — число атомов С, Н, О и галоидов в молекуле горючего;

К_н — коэффициент, учитывающий негерметичность помещения и неадиабатичность процесса горения. Допускается принимать К_н равным трем.

Таблица 2.1

Вид горючего вещества	Значение Z
Водород и нагретые выше температуры вспышки высокотемпературные органические теплоносители	1,0
Горючие газы	0,5
Легковоспламеняющиеся и горючие жидкости, нагретые до температуры вспышки и выше	0,3
Легковоспламеняющиеся и горючие жидкости, нагретые ниже температуры вспышки, при наличии возможности образования аэрозоля	0,3
Легковоспламеняющиеся и горючие жидкости, нагретые ниже температуры вспышки, при отсутствии возможности образования аэрозоля	0,0

Расчет избыточного давления для индивидуальных веществ и смесей

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

                                      (2.3)

где H_т — теплота сгорания, Дж/кг; ρ_в — плотность воздуха при начальной температуре Т₀, кг/м3; С_р — теплоемкость воздуха, Дж/(кг·К) [допускается принимать равной 1,01·103 Дж/(кг·К)]; Т₀ — начальная температура воздуха, К.

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

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

Данные для расчета

В помещение со свободным объемом V_св = 160 м³ при аварийной разгерметизации аппарата поступает 117,9 кг паров ацетона (определенных в соответствии с приложением И). Максимально возможная температура для данной климатической зоны t_р = 36 °С. Молярная масса ацетона М = 58,08 кг/кмоль. Химическая формула ацетона С₃Н₆O. Максимальное давление при сгорании стехиометрической паровоздушной смеси ацетона в замкнутом объеме Р_max = 572 кПа.

Расчет

Стехиометрический коэффициент кислорода в реакции сгорания ацетона равен

Стехиометрическая концентрация паров ацетона составит

%(об.).

Плотность паров ацетона ρ_п при расчетной температуре t_р равна

кг/м³

Тогда избыточное давление Dр при сгорании паровоздушной смеси ацетона для расчетной аварии составит

кПа

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

Данные для расчета

Через помещение, свободный объем которого V_св = 200 м³, проходит трубопровод с проходным сечением диаметром d_тр = 50 мм, по которому транспортируется водород Н₂ с максимальным расходом q = 5·10-3 м³/с при нормальных условиях и с максимальным давлением р_т = 150 кПа. Трубопровод оснащен системой автоматического отключения с временем срабатывания 2 с и с обеспечением резервирования ее элементов. Задвижки системы установлены перед стеной помещения в месте ввода трубопровода и за стеной данного помещения в месте вывода трубопровода. Длина отсекаемого участка трубопровода L_тр = 10м. Максимально возможная температура для данной климатической зоны t_p = 39 °С. Плотность водорода ρ_в при данной t_p равна 0,0787 кг/м³. Молярная масса водорода М = 2,016 кг/кмоль. Максимальное давление при сгорании стехиометрической газовоздушной смеси водорода в замкнутом объеме р_max = 730 кПа.

Расчет

Объем водорода, поступившего в помещение в результате аварийной разгерметизации трубопровода, будет равен

V_в = V_1т + V_2т = 0,01 + 0,02945 = 0,03945 м³ ,

V_1т = q · T = 5 · 10-3 · 2 = 0,01 м³,

м³.

Масса водорода, поступившего в помещение при расчетной аварии, составит m_в = V_в ρ_в = 0,03945 · 0,0787 = 3,105 · 10^-3 кг.

Стехиометрический коэффициент кислорода в реакции сгорания водорода равен

Стехиометрическая концентрация водорода составит

%(об.).

Избыточное давление D_р при сгорании водородовоздушной смеси, образующейся в результате расчетной аварии, равно

кПа

Расчет избыточного давления при сгорании пылевоздушной смеси в помещении

Избыточное давление при сгорании пылевоздушной смеси D_p, кПа, рассчитывают по формуле

,                                 (2.4)

где М— расчетная масса взвешенной в объеме помещения горючей пыли, образовавшейся в результате аварийной ситуации, кг; Н_т — теплота сгорания пыли, Дж/кг; p₀— начальное атмосферное давление, кПа (допускается принимать равным 101 кПа); Z — доля участия взвешенной горючей пыли при сгорании пылевоздушной смеси; V_св — свободный объем помещения, м3; ρ_В — плотность воздуха до сгорания пылевоздушной смеси при начальной температуре Т_о, кг/м³; С_р — теплоемкость воздуха, ДжДкг·Л) [допускается принимать равной 1010 Дж/(кг·А)]; Т_о — начальная температура воздуха, К, К_н — коэффициент, учитывающий негерметичность помещения и неадиабатичность процесса горения. Допускается принимать К_н равным трем.

К пылям, способным образовывать горючие пылевоздушные смеси, относят дисперсные материалы, характеризующиеся наличием показателей пожарной опасности: нижним концентрационным пределом распространения пламени, максимальным давлением, развиваемым при сгорании пылевоздушной смеси (более 50 кПа), и скоростью его нарастания, минимальным пожароопасным содержанием кислорода (менее 21 %).

Z рассчитывают по формуле

Z= 0,5 F,                                              (2.5)

где F— массовая доля частиц пыли размером менее критического, с превышением которого аэровзвесь становится взрывобезопасной, т. е. неспособной распространять пламя.

В отсутствие возможности получения сведений для расчета Z допускается принимать Z = 0,5.

Пример расчета избыточного давления при сгорании пыли

Пример 3. Определение избыточного давления, возникающего при сгорании полиэтиленовой пыли в помещении

Данные для расчета

В помещении горит полиэтиленовая пыль при следующих исходных условиях:

М_вз = 10 кг; M_ав = 90 кг; F = 0,3; H_т = 47·106 Дж/кг; V_св = 2000 м3; V_ав = 20 м³; Р_в = 1,2 кг/м³; T_о = 298 К; p_ст = 0,1 кг·м³.

Определяем Z

Z = 0,5F= 0,5 · 0,3 = 0,15.

Определяем М

отсюда следует, что М = 14 кг.

Принимая К_н = 3 и подставляя исходные данные в выражение для расчетного избыточного давления при сгорании пылевоздушной смеси, получим:

 кПа.

Приложение 3

МЕТОД РАСЧЕТА РАЗМЕРОВ ЗОН РАСПРОСТРАНЕНИЯ ОБЛАКА ГОРЮЧИХ ГАЗОВ И ПАРОВ ПРИ АВАРИИ

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

Мгновенный выброс СУГ может происходить при повреждении резервуара или иного аппарата, в котором СУГ находится под давлением. За счет внутренней энергии СУГ его массовая доля d мгновенно испаряется, образуя с капельками жидкости облако аэрозоля. За счет больших скоростей вихревых потоков происходит быстрое вовлечение в облако воздуха и быстрое испарение оставшейся части СУГ.

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

Ма₀ = ( l - d) M_gL_g / ( C_p.a ( T_a - T_g ) + X_wL_w),                       (3.1)

где М_g — масса выброшенного СУГ, кг; С_р.a — удельная теплоемкость воздуха, Дж/(кг·К); L_g — удельная теплота парообразования СУГ, Дж/кг; T_a — температура окружающего воздуха, К; Т_g — температура кипения СУГ при атмосферном давлении, К; Х_w — массовая доля водяных паров в воздухе; L_w — удельная теплота парообразования воды, Дж/кг.

d определяют из соотношения

d = 1 - ехр (- С_p.g (T_a - T_g) / L_g),                              (3.2)

где C_p.g — удельная теплоемкость СУГ, Дж/(кг · К).

Принимают, что образовавшееся облако дрейфует по ветру со скоростью v_d = 0,6n_в (n_в — скорость ветра) и имеет в начальный момент форму цилиндра, высота которого равна его радиусу. С течением времени высота облака уменьшается, а радиус растет. Изменение во времени радиуса, высоты облака и концентрации газа в нем в этой фазе (называемой фазой падения) определяется путем решения методом Рунге-Кутта системы обыкновенных дифференциальных уравнений:

dM_a / dt = ρ_a p r² a₂ a₃ n_в Ri-1 + 2 ρ_a а₁ (dr / dt) p r h,

dT / dt =((dMa / dt) C_p.a ( Ta - T ) + p r² ( Tg_r - T )1,333 ) / ( Ma C_p.a + Mg C_p.g ),

(3.3)

dr / dt = a₄ (gh (ρ_g.a - ρ_a) / r_g.a) 0,5,

где M_a — масса воздуха в облаке, кг; ρ_a — плотность воздуха, кг/м³; r — радиус облака, м;

а₁, a₂, a₃, a₄ — коэффициенты ( а₁ = 0,7, а₂ = 0,5, a₄ = 1,07, a₃ = 0,3 для классов устойчивости А—В (классы устойчивости даны по Паскуиллу, таблица 3.1); 0,24 — для С—В; 0,16 — для E—F);

Ri — число Ричардсона, определяемое из соотношения

Ri = (5,88 h₀,48 g / ( a32 nв₂ )) ( ρ_g.a - ρ_a ) / ρ_a ;

h — высота облака, м; Т— температура облака. К; Т_gr — температура земной поверхности. К; ρ_g.a  — плотность паровоздушного облака, кг/м³.

Таблица 3.1— Классы устойчивости атмосферы по Паскуиллу

Класс по Паскуиллу	Типичная скорость ветра, м/с	Описание погоды	Вертикальный градиент температуры, К/м
А	1	Безоблачно	>>>0,01
В	2	Солнечно и тепло	>>0,01
С	5	Переменная облачность в течение дня	>0,01
D	5	Облачный день или облачная ночь	»0,01
Е	3	Переменная облачность в течение ночи	<0,01
F	2	Ясная ночь	Инверсия (отрицательный градиент)

Решением системы вышеуказанных уравнений являются зависимости

M_a = M_a(t), Т= Т(t), ρ= ρ(t).

Для решения системы уравнений необходимы дополнительные соотношения

ρ_g.a = (M_a  + M_g ) / (M_a / ρ_a + Mg / ρ_g ) ( T_a / T ). (3.4)

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

(ρ_g.a - ρ_a ) / ρ_g.a  < 10-3 .                                   (3.5)

Зависимость h = h(t) находим из соотношения

h(t)=(M_a / ρ_a + M_g /ρ_g ) (T /T_a )( 1/(p r(t)²).                              (3.6)

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

              (3.7)

где s_y, s_z — среднеквадратичные отклонения, зависящие от величины xc - x0 ; х_c — координата центра облака в направлении ветра, м x₀ — координата точки окончания фазы падения, м; s_y (x_c - x₀); s_z (x_c - x₀) зависят от класса устойчивости по Паскуиллу. При x_c = x₀ принимается s_y0 = r / 2,14, s_z0 = h / 2,14; при x_c > x₀

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

(3.8)              

где Q = т· t _j,— масса СУГ в j-м элементарном объеме, кг; т — массовая скорость истечения СУГ, кг/с; x_j— координата центра j-го элементарного объема, м;  — среднеквадратичные отклонения распределения концентраций в j-м элементарном объеме, м.  - определяют аналогично  в 3.7

Пример расчета динамики паровоздушного облака в открытом пространстве

Пример 1. Для расчета динамики паровоздушного облака (движения в пространстве границы облака, определяемой НКПВ) принимается, что в некоторый момент времени t₀ начинается истечение пропана с массовой скоростью 1,3 кг/с, скорость ветра составляет 1 м/с, градиент температуры составляет 0,667 К/м.

t₀ — время начала истечения Рисунок 3.1 — Границы паровоздушного облака по НКПВ на различные моменты времени от начала истечения

Результаты расчета границы облака для двух значений времени t₀ + 10 с и t₀ + 300 с представлены на рисунке 3.1. Процедура расчета, реализованная на ПЭВМ, представлена на блок-схеме (рисунок 3.2).

Рисунок 3.2 — Алгоритм расчета параметров паровоздушного облака

Приложение 4

МЕТОД РАСЧЕТА ИНТЕНСИВНОСТИ ТЕПЛОВОГО ИЗЛУЧЕНИЯ

Расчет интенсивности теплового излучения «огненного шара» q, кВт/м2, проводят по формуле

q = E_f·F_q·t,                                       (4.1)

где E_f — среднеповерхностная плотность теплового излучения пламени, кВт/м2; F_q — угловой коэффициент облученности; t  - коэффициент пропускания атмосферы

E_f определяют на основе имеющихся экспериментальных данных. Допускается принимать E_f  равным 450 кВт/м².

F_q рассчитывают по формуле

,            (4.2)

где Н— высота центра «огненного шара», м; D_s — эффективный диаметр «огненного шара», м; r — расстояние от облучаемого объекта до точки на поверхности земли непосредственно под центром «огненного шара», м.

Эффективный диаметр «огненного шара» D_s рассчитывают по формуле

D_s =5,33 m 0,327,                              (4.3)

где т — масса горючего вещества, кг.

H определяют в ходе специальных исследований. Допускается принимать H равной D_s/2.

Время существования «огненного шара» t_s, с, рассчитывают по формуле

t_s = 0,92 m 0,303.                                              (4.4)

Коэффициент пропускания атмосферы т рассчитывают по формуле

t = ехр [-7,0 · 10-4 (- Ds / 2)].                          (4.5)

Пример определения время существования «огненного шара» и интенсивность теплового излучения от него на расстоянии 500 м при разрыве сферической емкости с пропаном объемом 600 м³ в очаге пожара.

Данные для расчета

Объем сферической емкости 600 м³. Плотность жидкой фазы 530 кг/м³. Степень заполнения резервуара жидкой фазы 80 %. Расстояние от облучаемого объекта до точки на поверхности земли непосредственно под центром «огненного шара» 500 м.

Расчет

Находим массу горючего т в «огненном шаре» по формуле

т = Vra = 600 · 530 · 0,8 = 2,54 · 105 кг,

где V— объем резервуара, м³ (V =  600 м³); r — плотность жидкой фазы, кг/м³ (r = 530 кг/м³); a — степень заполнения резервуара (a = 0,8).

По формуле (4.3) определяем эффективный диаметр «огненного шара» D_s

D_s = 5,33 (2,54 · 105)0,327 = 312 м.

По формуле (4.2), принимая H = Ds /2 = 156 м, находим угловой коэффициент облученности F_q

По формуле (4.5) находим коэффициент пропускания атмосферы t: t = ехр [-7,0 · 10-4 (  - 312/2)] = 0,77.

По формуле (4.1), принимая E_f= 450 кВт/м², находим интенсивность тепловогоизлучения q

q = 450 · 0,037 · 0,77 = 12,9 кВт/м².

По формуле (4.4) определяем время существования «огненного шара» t_s

t_s = 0,92 (2,54 · 105)0,303 = 40 с.

Приложение 5

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

Исходя из рассматриваемого сценария аварии, определяют массу т, кг, горючих газов и (или) паров, вышедших в атмосферу из технологического аппарата (приложение 2).

Избыточное давление D_p, кПа, развиваемое при сгорании газопаровоздушных смесей, рассчитывают по формуле

,             (5.1)

где р₀ — атмосферное давление, кПа (допускается принимать равным 101 кПа); r — расстояние от геометрического центра газопаровоздушного облака, м; m_пp — приведенная масса газа или пара, кг, рассчитанная по формуле

m_пр = (Q_сг / Q₀)m_г,п Z, (5.2)

где Qсг — удельная теплота сгорания газа или пара, Дж/кг; Z— коэффициент участия, который допускается принимать равным 0,1; Q₀— константа, равная 4,52 · 106 Дж/кг; m_г,п — масса горючих газов и (или) паров, поступивших в результате аварии в окружающее пространство, кг.

Импульс волны давления i, Па · с, рассчитывают по формуле

.                                       (5.3)

Пример 1. Рассчитать избыточное давление и импульс волны давления при выходе в атмосферу пропана, хранящегося в сферической емкости объемом 600 м³, на расстоянии 500 м от нее.

Данные для расчета

Объем емкости 600 м³. Температура 20 °С. Плотность сжиженного пропана 530 кг/м³. Степень заполнения емкости 80 % (по объему). Удельная теплота сгорания пропана 4,6 · 10⁷ Дж/кг. Принимается, что в течение времени, необходимого для выхода сжиженного газа из емкости, весь пропан испаряется.

Расчет

Находим приведенную массу m_пр  по формуле (5.2):

m_пр  = 4,6 · 107 / 4,52 · 106 ·(0,8· 530 · 600) · 0,1 = 2,59 · 105 кг.

Находим избыточное давление D_p по формуле (5.1)

D_p = 101 [0,8 (2,59 · 105)0,33 / 500 + 3 (2,59 · 105) 0,66 / 5002 + 5 (2,59 · 105) / 5003] = 16,2 кПа.

Находим импульс волны давления i по формуле (5.3):

i = 123 (2,59 · 105)0,66 / 500 = 1000 Па · с.

Приложение 6

СПРАВОЧНОЕ

Таблица 6.1— Пороговое количество веществ для технологических производств

Наименование веществ или групп веществ	Категория опасности по ГОСТ 19433	Пороговое количество веществ, т, не менее
Газы горючие сжатые, сжиженные и растворенные под давлением
Ацетилен С2Н2	231	50
Водород Н2	231	50
Сернистый водород H2S	241	50
Оксид этилена (СН2)2О	241	50
Аммиак NH3	241	500
Все остальные сжатые, сжиженные и растворенные под давлением горючие газы	231,232,241,911	200
Легковоспламеняющиеся и горючие жидкости
Оксид пропилена С3Н6О	311	50
Все остальные легковоспламеняющиеся и горючие жидкости	311, 312, 321, 322, 324, 325, 335, 314, 315, 323, 331	200
Твердые вещества
Саморазлагающиеся вещества	415,416,417,418, 521, 522, 523	10
Окисляющие вещества
Кислород жидкий 02	212	2000
Нитрат аммония NН4NО3	511	5000
Хлорат натрия NaClО3	511	250
Пероксид метилэтилкетона (концентрация более 60%)	523	250
Пероксид метилизобутилкетона (концентрация более 60 %)	523	50
Надуксусная кислота (концентрация более 60 %) СН3СОООН	522	50
Хлор Сl2	222	50
Бром Вr2	832	500
Оксиды азота N2О, N2O2	512	50
Примечания 1 Наименования классов (подклассов) опасных веществ приведены в соответствии с ГОСТ 19433. 2 Пороговые количества опасных веществ, относящихся по ГОСТ 19433 к категориям 434 и 437 (самовозгорающиеся твердые вещества), классу 8 (едкие и коррозионно-активные вещества), категориям 436 и 438 (твердые легковоспламеняющиеся вещества), категориям 425, 913 и 921 (вещества, горючие газы при взаимодействии с водой), должны согласовываться со специально уполномоченными государственными органами по безопасности в промышленности. 3 При совместном использовании несовместимых по ГОСТ 12.1.004 веществ их пороговое количество должно быть согласовано со специально уполномоченными государственными органами по безопасности в промышленности

Приложение 7