МИНИСТЕРСТВО
РОССИЙСКОЙ ФЕДЕРАЦИИ ПО ДЕЛАМ ФЕДЕРАЛЬНОЕ ГОСУДАРСТВЕННОЕ УЧРЕЖДЕНИЕ РУКОВОДСТВО
|
Наименование оборудования |
Событие, инициирующее аварию |
Диаметр отверстия истечения, мм |
Частота разгерметизации, год-1 |
Резервуары, емкости, сосуды и аппараты под давлением |
Разгерметизация с последующим истечением жидкости, газа или двухфазной среды |
5 |
4,0 • 10-5 |
12,5 |
1,0 • 10-5 |
||
25 |
5,0 • 10-6 |
||
50 |
5,0 • 10-6 |
||
Полное разрушение |
2,0 • 10-6 |
||
Внешнее воздействие, приводящее к реализации огненного шара |
- |
2,5 • 10-5 |
|
Насосы (центробежные) |
Разгерметизация с последующим истечением жидкости или двухфазной среды |
5 |
4,3 • 10-3 |
12,5 |
6,1 • 10-4 |
||
25 |
5,1 • 10-4 |
||
50 |
2,0 • 10-4 |
||
Полное разрушение |
1,0 • 10-4 |
||
Компрессоры (центробежные) |
Разгерметизация с последующим истечением газа |
5 |
1,1 • 10-2 |
12,5 |
1,3 • 10-3 |
||
25 |
3,9 • 10-4 |
||
50 |
1,3 • 10-4 |
||
Полное разрушение |
1,0 • 10-4 |
||
Резервуары для хранения ЛВЖ и ГЖ при давлении, близком к атмосферному |
Разгерметизация с последующим истечением жидкости в обвалование |
Все типы |
1,0 • 10-4 |
Квазимгновенное разрушение |
- |
5,0 • 10-6 |
|
Резервуары с плавающей крышей |
Пожар в кольцевом зазоре по периметру резервуара |
- |
4,6 • 10 3 |
Пожар по всей поверхности резервуара |
- |
9,3 • 10-4 |
|
Резервуары со стационарной крышей |
Пожар на дыхательной арматуре |
- |
9,0 • 10-5 |
Пожар по всей поверхности резервуара |
- |
9,0 • 10-5 |
Примечания:
1. Здесь и далее под полным разрушением подразумевается утечка с диаметром истечения, соответствующим максимальному диаметру подводящего или отводящего трубопровода, или разрушения резервуара, емкости, сосуда или аппарата.
2. При определении частоты разгерметизации фильтров и кожухотрубных теплообменников указанное оборудование можно рассматривать как аппараты под давлением.
3. Аппараты воздушного охлаждения можно рассматривать как участки технологических трубопроводов, длина которых соответствует суммарной длине труб в пучках теплообменника.
Диаметр трубопровода, мм |
Частота утечек, (м-1 • год-1) |
|||
Малая (диаметр отверстия 12,5 мм) |
Средняя (диаметр отверстия 25 мм) |
Большая (диаметр отверстия 50 мм) |
Разрыв |
|
50 |
5,7 • 10-6 |
2,4 • 10-6 |
9,4 • 10-7 |
4,7 • 10-7 |
100 |
2,8 • 10-6 |
1,2 • 10-6 |
4,7 • 10-7 |
2,4 • 10-7 |
150 |
1,9 • 10-6 |
7,9 • 10-7 |
3,1 • 10-7 |
1,6 • 10-7 |
250 |
1,1 • 10-6 |
4,7 • 10-7 |
1,9 • 10-7 |
9,4 • 10-7 |
600 |
4,7 • 10-7 |
2,0 • 10-7 |
7,9 • 10-8 |
3,9 • 10-8 |
900 |
3,1 • 10-7 |
1,3 • 10-7 |
5,2 • 10-8 |
2,6 • 10-8 |
1200 |
2,4 • 10-7 |
9,8 • 10-8 |
3,9 • 10-8 |
2,0 • 10-8 |
Таблица 1.3
Диаметр трубопровода, мм |
Частота аварий (м-1 • год-1) при размерах дефектного отверстия в виде ромба со сторонами (мм) |
||
306 × 38 |
765 × 96 |
1530 × 190 |
|
1020 |
8,8 • 10-8 |
5,6 • 10-8 |
1,6 • 10-8 |
Таблица 1.4
Наименование объекта |
Частота возникновения пожара, (м-1 • год-1) |
Электростанции |
2,2 • 10-5 |
Склады химической продукции |
1,2 • 10-5 |
Склады многоменклатурной продукции |
9,0 • 10-5 |
Инструментально-механические цеха |
0,6 • 10-5 |
Цеха по обработке синтетического каучука и искусственных волокон |
2,65 • 10-5 |
Литейные и плавильные цеха |
1,89 • 10-5 |
Цеха по переработке мясных и рыбных продуктов |
1,53 • 10"5 |
Цеха горячей прокатки металлов |
1,89 • 10-5 |
Текстильные производства |
1,53 • 10-5 |
Логическое дерево событий предназначено для графического отображения общего характера развития возможных аварийных ситуаций и аварий с отражением причинно-следственной взаимосвязи событий в зависимости от специфики опасности объекта оценки риска с учетом влияния на них имеющихся защитных мероприятий и является основой для оценки риска. Сценарий возникновения и развития аварийной ситуации и аварии на логическом дереве отражается в виде последовательности событий от исходного до конечного события (ветвь дерева событий).
Данный метод позволяет проследить развитие возможных аварийных ситуаций и аварий, возникающих вследствие реализации событий, инициирующих аварийную ситуацию. Анализ дерева событий представляет собой «осмысливаемый вперед» процесс, то есть процесс, при котором исследование развития аварийной ситуации начинается с исходного события с рассмотрением цепочки последующих событий, приводящих к возникновению аварии.
При построении логических деревьев событий учитываются следующие положения:
- выбирается аварийная ситуация, которая может повлечь за собой возникновение аварии с пожаром (взрывом) с дальнейшим его развитием;
- развитие аварийной ситуации и аварии должно рассматриваться постадийно с учетом места ее возникновения на объекте оценки риска, уровня потенциальной опасности каждой стадии и возможности ее локализации и ликвидации. На логическом дереве событий стадии развития аварийной ситуации и аварии могут отображаться в виде прямоугольников или других геометрических фигур с находящимися в них краткими названиями этих стадий;
- переход с рассматриваемой стадии на новую определяется возможностью либо локализации аварии на рассматриваемой стадии, либо развития аварии, связанного с вовлечением расположенных рядом технологического оборудования, помещений, сооружений и т.п. в результате влияния на них опасных факторов пожара (взрыва), возникших на рассматриваемой стадии. Вероятности переходов аварии со стадии на стадию одной ветви или с ветви на ветвь определяется, исходя из величин зон поражения опасными факторами пожара (взрыва), объемно-планировочных решений и конструктивных особенностей оборудования и сооружений промышленного предприятия. При этом каждой стадии иногда присваивается идентификационный номер, отражающий последовательность переходов со стадии на стадию;
- переход со стадии на стадию отображается в виде соединяющих линий со стрелками, указывающими направления развития пожароопасной ситуации и последующего пожара. При этом соединения стадий должны отражать вероятностный характер события с выполнением условий «И», «ИЛИ» и «И/ИЛИ» (условие «И» - взаимосвязь событий, условие «ИЛИ» - независимость событий, а условие «И/ИЛИ» - возможность частичной реализации либо условия «И», либо «ИЛИ» ввиду многообразия путей дальнейшего развития). Отображение взаимосвязи событий по условию «И» выполняется сплошной линией без ответвлений, по условию «ИЛИ» - сплошной линией с ответвлениями к взаимоисключающим друг друга событиям (стадиям), по условию «И/ИЛИ» - пунктирной линией с ответвлениями к потенциально возможным стадиям;
- для каждой стадии должен устанавливаться уровень ее опасности, характеризующийся возможностью перехода аварийной ситуации или аварии на соседние с пожароопасным участки объекта;
- при повторении одним из путей части другого пути развития аварии для упрощения построения логического дерева событий иногда вводят обозначение, представляющее собой соответствующую линию со стрелкой и надпись «на стадию (код последующей стадии)». При этом в случае перехода со стадии с более высоким уровнем опасности на стадию с меньшим уровнем следует дополнять код стадии с меньшим уровнем опасности кодом, соответствующим стадии с более высоким уровнем опасности, указываемым в скобках.
При анализе логических деревьев событий руководствуются следующими положениями:
- возможность предотвращения дальнейшего развития аварийной ситуации и аварии зависит от количества стадий / времени их протекания (т.е. от длины пути развития аварийной ситуации и аварии). Это объясняется большей вероятностью успешной ликвидации аварийной ситуации и аварии, связанной с увеличением времени на локализацию аварийной ситуации и аварии и количеством стадий, на которых эта локализация возможна;
- наличие у стадии нескольких разветвлений по принципу «И» или «И/ИЛИ» свидетельствует о приобретении аварийной ситуацией и аварией «цепного» характера развития, т.е. одновременного развития аварии по нескольким путям с дальнейшим увеличением их количества (так называемый «эффект домино»), что в значительной мере затрудняет успешную локализацию и ликвидацию пожара;
- наличие у стадии разветвлений по принципу «ИЛИ», одно из которых приходит на стадию локализации аварийной ситуации или аварии (например, тушение очага пожара, своевременное обнаружение утечки и ликвидация пролива, перекрытие запорной арматуры и т.п.), свидетельствует о возможности предотвращения дальнейшего развития аварийной ситуации и аварии по этому пути. Отсутствие стадии локализации аварийной ситуации и аварии на разветвлениях по принципу «ИЛИ» свидетельствует о невозможности приостановления дальнейшего развития аварийной ситуации и аварии на разветвляемой стадии;
- чем больше возможных путей развития аварийной ситуации и аварии приходит на одну стадию, тем больше вероятность возникновения этой стадии.
Для оценки вероятности перехода аварии со стадии на стадию, прежде всего, необходимо определение условной вероятности реализации различных ветвей дерева событий. Далее необходимо определение вероятностей срабатывания соответствующих средств предотвращения или локализации аварии. Кроме того, необходима оценка вероятности поражения расположенного в зоне возникновения аварии технологического оборудования и сооружений объекта в результате воздействия на них опасных факторов пожара.
В табл. 2.1 приведены условные вероятности мгновенного воспламенения и воспламенения с задержкой по времени в зависимости от массовой скорости истечения газа, двухфазной среды или жидкости при разгерметизации типового технологического оборудования на промышленных предприятиях.
Таблица 2.1
Массовый расход истечения, кг • с-1 |
Вероятность мгновенного воспламенения |
Вероятность последующего воспламенения при отсутствии мгновенного воспламенения |
Вероятность сгорания с образованием избыточного давления при последующем воспламенении |
||||||||
Диапазон |
Номинальное среднее значение |
газ |
двухфазная смесь |
жидкость |
газ |
двухфазная смесь |
жидкость |
газ |
двухфазная смесь |
жидкость |
|
Малый (<1) |
0,5 |
0,005 |
0,005 |
0,005 |
0,005 |
0,005 |
0,005 |
0,080 |
0,080 |
0,050 |
|
Средний (1-50) |
10 |
0,035 |
0,035 |
0,015 |
0,036 |
0,036 |
0,015 |
0,240 |
0,240 |
0,050 |
|
Большой (>50) |
100 |
0,150 |
0,150 |
0,040 |
0,176 |
0,176 |
0,042 |
0,600 |
0,600 |
0,050 |
|
Полный разрыв |
Не определено |
0,200 |
0,200 |
0,050 |
0,240 |
0,240 |
0,061 |
0,540 |
0,540 |
0,100 |
|
Условные вероятности воспламенения при различных диаметрах отверстия истечения при разгерметизации технологического оборудования на промышленных предприятиях приведены в табл. 2.2 и 2.3. При этом под задержкой воспламенения следует понимать время более одной минуты.
Диаметр истечения, мм |
Условная вероятность воспламенения |
|
Газ/двухфазная среда |
Жидкость |
|
5 |
0,1 |
0,04 |
25 |
0,2 |
0,08 |
50 |
0,3 |
0,12 |
Разрыв |
0,5 |
0,2 |
Таблица 2.3
Условная вероятность мгновенного воспламенения и воспламенения с задержкой при различных диаметрах отверстия истечения
Диаметр истечения, мм |
Мгновенное воспламенение |
Воспламенение с задержкой |
5 |
0,1 |
0,9 |
25 |
0,2 |
0,8 |
50 |
0,3 |
0,7 |
Разрыв |
0,4 |
0,6 |
Для особо опасных легковоспламеняющихся жидкостей (ЛВЖ) (температура вспышки tвсп ≤ 28 °С) следует использовать условные вероятности воспламенения как для двухфазной среды.
В табл. 2.4 приведены условные вероятности различных сценариев развития аварий с выбросом из технологического оборудования сжиженных углеводородных газов при полном разрушении оборудования.
Таблица 2.4
Сценарий аварии |
Вероятность |
Сценарий аварии |
Вероятность |
Факел |
0,06 |
Сгорание с развитием |
|
Огненный шар |
0,70 |
избыточного давления |
0,01 |
Горение пролива |
0,03 |
Без горения |
0,03 |
Сгорание облака |
0,17 |
Итого |
1 |
В настоящем приложении представлены методики оценки опасных факторов, реализующихся при различных сценариях аварий с пожарами на промышленных предприятиях с наличием горючих газов, легковоспламеняющихся и горючих жидкостей, основанные на современных отечественных и зарубежных теоретических и экспериментальных данных, нормативных документах и руководствах.
В случае, если при анализе сценариев аварий возникает необходимость учесть сценарии, методы оценки опасных факторов для которых в настоящем документе не изложены, или возникает необходимость в отдельных случаях применять методы, отличающиеся от описанных ниже, возможно использование иных методик при соответствующем обосновании необходимости их применения.
3.1. Истечение жидкостей и газов
3.1.1. Истечение жидкости
Рассматривается резервуар, находящийся в обваловании (рис. 3.1).
Рис. 3.1. Схема для расчета истечения жидкости из отверстия в резервуаре
Вводятся следующие допущения:
- истечение через отверстие однофазное;
- резервуар имеет постоянную площадь сечения по высоте;
- диаметр резервуара много больше размеров отверстия,
- размеры отверстия много больше толщины стенки;
- поверхность жидкости внутри резервуара горизонтальна;
- температура жидкости остается постоянной в течение времени истечения.
Массовый расход жидкости G (кг • с-1) через отверстие во времени t (с) описывается соотношением
(3.1)
где G0 - массовый расход в начальный момент времени, кг • с-1, описываемый выражением
(3.2)
где ρ - плотность жидкости, кг • м-3;
g - ускорение свободного падения (9,81 м • с-2);
μ - коэффициент истечения;
Ahol -площадь отверстия, м2;
hhol - высота расположения отверстия, м;
AR - площадь сечения резервуара, м2;
h0 - начальная высота столба жидкости в резервуаре, м.
Зависимость высоты столба жидкости в резервуаре h (м) от времени t описывается формулой
(3.3)
Условия перелива струи жидкости (при h0 > hhol) через обвалование может быть приближенно записано следующим образом:
(3.4)
где Н - высота обвалования, м;
L - расстояние от стенки резервуара до обвалования, м.
Количество жидкости m (кг), перелившейся через обвалование за полное время истечения, описывается выражением
(3.5)
где tpour - время, в течение которого жидкость переливается через обвалование, с (т.е. время, в течение которого выполняется условие (3.4)).
Величина tpour описывается выражением
(3.6)
где а, b, с - параметры, описываемые формулами:
(3.7)
(3.8)
(3.9)
В случае, если жидкость в резервуаре находится под избыточным давлением ΔР (Па), величина мгновенного
массового расхода G0 (кг • с-1) должна быть описана выражением
(3.10)
Для определения количества жидкости, перелившейся через обвалование, и времени перелива следует проинтегрировать соответствующую систему уравнений, где величина ΔР может, вообще говоря, быть переменной.
3.1.2. Истечение сжатого газа
Массовая скорость истечения сжатого газа из резервуара описывается следующими соотношениями:
докритическое истечение
(3.11)
(3.12)
сверхкритическое истечение
(3.13)
(3.14)
где G - массовый расход, кг • с-1;
Ра - атмосферное давление, Па;
РV - давление газа в резервуаре, Па;
γ - показатель адиабаты газа;
Ahol - площадь отверстия, м2;
μ - коэффициент истечения (0,6-0,8);
ρV - плотность газа в резервуаре при давлении РV, кг • м-3.
3.1.3. Истечение сжиженного газа из отверстия в резервуаре
Массовая скорость истечения паровой фазы GV (кг • с-1) описывается формулой
(3.15)
где μ - коэффициент истечения;
Ahol - площадь отверстия, м2;
Рс - критическое давление сжиженного газа, Па;
М - молярная масса, кг • моль-1;
R - универсальная газовая постоянная, равная 8,31 Дж • К-1 • моль-1;
Тс - критическая температура сжиженного газа, К;
PR = РV/Рс - безразмерное давление сжиженного газа в резервуаре;
РV - давление сжиженного газа в резервуаре, Па.
Массовую скорость истечения паровой фазы можно также рассчитывать по формулам (3.11)-(3.14).
Массовая скорость истечения жидкой фазы G (кг • с-1) описывается формулой
(3.16)
где ρL - плотность жидкой фазы, кг • м-3;
ρV - плотность паровой фазы, кг • м-3;
TR = Т/Тс - безразмерная температура сжиженного газа;
Т - температура сжиженного газа в резервуаре, К.
3.1.4. Растекание жидкости при квазимгновенном разрушении резервуара
Под квазимгновенным разрушением резервуара следует понимать внезапный (в течение секунд или долей секунд) распад резервуара на приблизительно равные по размеру части.
При таких аварийных ситуациях часть хранимого в резервуаре продукта может перелиться через обвалование.
Ниже представлена математическая модель, позволяющая оценить долю жидкости, перелившейся через обвалование при квазимгновенном разрушении резервуара. Приняты следующие допущения:
- рассматривается плоская одномерная задача;
- время разрушения резервуара много меньше характерного времени движения гидродинамической волны до обвалования;
- жидкость является невязкой;
- трение жидкости о поверхность земли отсутствует;
- поверхность земли является плоской, горизонтальной. Система уравнений, описывающих движение жидкости, имеет вид
где h - высота столба жидкости над фиксированным уровнем, м;
hG - высота подстилающей поверхности над фиксированным уровнем, м;
u - средняя по высоте скорость движения столба жидкости, м • с-1;
х - координата вдоль направления движения жидкости, м;
t - время, с;
g - ускорение свободного падения (9,81 м • с-2).
Граничные условия с учетом геометрии задачи (рис. 3.2) имеют вид
(3.19)
(3.20)
(3.21)
где а - высота обвалования.
Массовая доля жидкости Q (%), перелившейся через обвалование к моменту времени Т, описывается выражением
(3.22)
где uN - средняя по высоте скорость движения столба жидкости при х = b, м • с-1;
hN - высота столба жидкости при х = b, м;
h0 - начальная высота столба жидкости в резервуаре, м;
R - ширина резервуара, м.
Рис. 3.2. Типичная картина движения жидкости в обваловании при квазимгновенном разрушении резервуара:
- - - - - - уровень начального столба жидкости;
______ - уровень жидкости в промежуточный
момент времени (результаты расчета)
Один из возможных способов решения системы (3.17) с граничными условиями (3.18)-(3.21) - численное решение с использованием схемы Мак-Кормака, представляющей собой частный случай явной схемы типа «предиктор-корректор».
График указанной функции вместе с литературными экспериментальными данными представлен на рис. 3.3. Для практических целей представляется достаточным использование этого графика.
Рис. 3.3. Зависимость доли перелившейся через обвалование жидкости Q от параметров a/h0:
1 - расчет; 2 - эксперимент
3.2. Количественная оценка массы горючих веществ, поступающих в окружающее пространство в результате возникновения аварийных ситуаций
3.2.1. Общие положения
Количество поступивших в замкнутое или свободное пространство веществ, которые могут образовать взрывоопасные паровоздушные смеси или проливы горючих жидкостей, определяется, исходя из следующих предпосылок:
а) происходит расчетная авария одного из резервуаров или трубопроводов;
б) все содержимое резервуара (трубопроводов) или часть продукта (при соответствующем обосновании) поступает в замкнутое или свободное пространство. При этом в случае наличия на объекте нескольких резервуаров (с различными объемами) расчет следует проводить для вариантов с разгерметизацией каждого единичного резервуара;
в) при разгерметизации резервуара происходит одновременно утечка веществ из трубопроводов, питающих резервуар по прямому и обратному потоку в течение времени, необходимого для отключения трубопроводов. Расчетное время отключения трубопроводов определяется в каждом конкретном случае, исходя из реальной обстановки, и должно быть минимальным с учетом паспортных данных на запорные устройства и их надежности, характера технологического процесса и вида расчетной аварии. Под «временем отключения» следует понимать промежуток времени от начала возможного поступления горючего вещества из трубопровода (перфорация, разрыв, изменение номинального давления и т.п.) до полного прекращения поступления жидкости в окружающее пространство;
г) при расчетах принимается нулевая подвижность окружающего воздуха (безветрие);
д) в качестве расчетной температуры при аварийной ситуации с наземным расположением оборудования принимается максимально возможная температура воздуха в соответствующей климатической зоне, а при аварийной ситуации с подземным расположением оборудования - температура грунта, условно равная максимальной среднемесячной температуре окружающего воздуха в наиболее теплое время года;
е) длительность испарения жидкости с поверхности пролива принимается равной времени ее полного испарения, но не более 3 600 с. Для относительно небольших проливов топлива (до 20 кг) время испарения допускается принимать равным 900 с, поскольку столь небольшие проливы могут быть достаточно эффективно удалены обслуживающим персоналом. Кроме того, в запас надежности идет неучет подвижности воздуха и уменьшение скорости испарения жидкости со временем вследствие ее охлаждения.
При проведении расчетов допускается использование справочных данных, опубликованных головными научно-исследовательскими организациями в области пожарной безопасности или выданных Государственной службой стандартных справочных данных. Допускается использование показателей пожаровзрывоопасности для смесей веществ и материалов по наиболее опасному компоненту.
Ниже приводятся основные расчетные формулы для определения масс горючих веществ, поступающих в открытое или замкнутое пространство в результате аварийных ситуаций.
3.2.2. Разгерметизаций надземного резервуара
Масса жидкости, поступившей в окружающее пространство при разгерметизации резервуара определяется по формуле
(3.23)
где mа - масса жидкости, кг;
ρL - плотность жидкости, кг • м-3;
VR - объем жидкости в резервуаре, м3.
Масса жидкости, поступившей самотеком при полном разрушении трубопровода, выходящего из резервуара, определяется по формуле
(3.24)
(3.25)
(3.26)
где GL - начальный расход жидкости, истекающей из резервуара через разгерметизированный трубопровод, кг • с-1;
μ - коэффициент истечения;
τ - расчетное время отключения трубопроводов, связанных с местом разгерметизации, с;
dP -диаметр трубопроводов, м (в случае различных диаметров трубопроводов, связанных с местом разгерметизации, объем выходящего топлива рассчитывается для каждого трубопровода в отдельности);
Li - длина i-го участка трубопровода от запорного устройства до места разгерметизации, м;
n - число участков трубопроводов, связанных с местом разгерметизации;
ΔPR - напор столба жидкости в резервуаре, Па;
hL - высота столба жидкости (от верхнего уровня жидкости в резервуаре до уровня места разгерметизации), м;
g - ускорение свободного падения (g = 9,81 м • с-2).
3.2.3. Масса паров ЛВЖ, выходящих через дыхательную арматуру
В случае наполнения резервуара массу паров вычисляют по формуле
(3.27)
(3.28)
где mV - масса выходящих паров ЛВЖ, кг;
ρV - плотность паров ЛВЖ, кг • м-3;
PS - давление насыщенных паров топлива при расчетной температуре, кПа, определяемое по справочным данным;
P0 - атмосферное давление (допускается принимать равным 101 кПа);
VR - геометрический объем резервуара, м3;
М - молярная масса топлива, кг • кмоль-1;
V0 - мольный объем, равный 22,413 м3 • кмоль-1;
t0 - расчетная температура, °С.
3.2.4. Масса паров ЛВЖ при испарении со свободной поверхности в резервуаре
Массу паров ЛВЖ при испарении со свободной поверхности в резервуаре вычисляют по формуле
(3.29)
где GV - расход паров ЛВЖ, кг • с-1, определяемый соотношением
GV = FRW; (3.30)
τE - время поступления паров из резервуара, с;
FR - максимальная площадь поверхности испарения ЛВЖ в резервуаре, м2;
W - интенсивность испарения ЛВЖ, кг • м-2 • с-1).
3.3. Максимальные размеры взрывоопасных зон
Радиус (Rнкпр, м) и высота (Zнкпp, м) зоны, ограничивающие область концентраций, превышающих нижний концентрационный предел распространения пламени (НКПР), при неподвижной воздушной среде рассчитывают по формулам:
для горючих газов (ГГ)
(3.31)
(3.32)
где mг - масса ГГ, поступившего в открытое пространство при аварийной ситуации, кг;
ρг - плотность ГГ при расчетной температуре и атмосферном давлении, кг • м-3;
Снкпр - нижний концентрационный предел распространения пламени ГГ, % (об.);
для паров ЛВЖ
, (3.33)
, (3.34)
где Мп - масса паров ЛВЖ, поступивших в открытое пространство за время полного испарения, но не более 3600 с, кг;
ρп - плотность паров при расчетной температуре, кПа;
Рн - давление насыщенных паров при расчетной температуре, кПа;
К = Т/3600;
Т - продолжительность поступления паров в открытое пространство, с;
Снкпр - нижний концентрационный предел распространения пламени паров, % (об.).
За начало отсчета горизонтального размера зоны принимают внешние габаритные размеры пролива.
3.4. Определение параметров волны давления при взрыве облака топливно-воздушной смеси
Методика количественной оценки параметров воздушных ударных волн при взрыве топливно-воздушных смесей (ТВС) распространяется на случаи выброса топлива в атмосферу из наружных установок.
Основными структурными элементами алгоритма расчетов являются:
- определение ожидаемого режима взрывного превращения облака ТВС;
- расчет максимального избыточного давления и импульса фазы сжатия воздушных ударных волн для различных режимов;
- определение дополнительных характеристик взрывной нагрузки;
- оценка поражающего воздействия.
Исходными данными для расчета параметров ударных волн при взрыве облака ТВС являются:
- вид топлива, содержащегося в облаке;
- агрегатное состояние смеси (газовая или гетерогенная);
- концентрация горючего в смеси (Сг);
-стехиометрическая концентрация горючего с воздухом (Сст);
- масса топлива, содержащегося в облаке (Мт);
- удельная теплота сгорания топлива (Еуд);
- скорость звука в воздухе С0 (обычно принимается равной 340 м • с-1);
-информация об окружающем пространстве;
- эффективный энергозапас горючей смеси (Е), вычисляемый по соотношению
(3.35)
При расчете параметров взрыва облака, расположенного на поверхности земли, величина эффективного энергозапаса удваивается.
3.4.1. Определение ожидаемого режима взрывного превращения
Ожидаемый режим взрывного превращения зависит от типа топлива и геометрических характеристик окружающего пространства.
Классификация горючих смесей по степени чувствительности
Вещества, способные к образованию горючих смесей с воздухом, по своей чувствительности к возбуждению взрывных процессов разделены на четыре класса.
Класс 1. Особо чувствительные вещества (размер детонационной ячейки менее 2 см). К этому классу относятся такие вещества как Н2, C2H2, C2H4O, С3Н6О, R-NO2.
Класс 2. Чувствительные вещества (размер детонационной ячейки лежит в пределах от 2 до 10 см). Типичными представителями этого класса являются С3Н8, С4Н10, С2Н6, С3Н6, С4Н8.
Класс 3. Среднечувствительные вещества (размер детонационной ячейки лежит в пределах от 10 до 40 см). Типичными представителями этого класса являются гексан (С6Н14), октан, изооктан, пары и распылы бензина, LNG.
Класс 4. Слабочувствительные вещества (размер детонационной ячейки больше 40 см). Типичными представителями этого класса являются метан, декан, распылы дизтоплива, керосина, бензола.
Классификация наиболее распространенных в промышленном производстве горючих веществ приведена в табл. 3.1. В случае, если вещество не внесено в классификацию, его следует классифицировать по аналогии с имеющимися в списке веществами, а при отсутствии информации о свойствах данного вещества, его следует отнести к классу 1, т.е. рассматривать наиболее опасный случай.
Класс 1 |
Класс 2 |
Класс 3 |
Класс 4 |
Ацетилен |
Акрилонитрил |
Ацетальдегид |
Бензол |
Винилацетилен |
Акролеин |
Ацетон |
Декан |
Водород |
Бутан |
Бензин |
Дизтопливо |
Гидразин |
Бутилен |
Винилацетат |
о-Дихлорбензол |
Изопропилнитрат |
Бутадиен |
Винилхлорид |
Додекан |
Метилацетилен |
1,3-Пентадиен |
Гексан |
Керосин |
Нитрометан |
Пропан |
Генераторный газ |
Метан |
Окись пропилена |
Пропилен |
Изооктан |
Метилбензол |
Окись этилена |
Сероуглерод |
Метиламин |
Метилмеркаптан |
Этилнитрат |
Этан |
Метилацетат |
Метилхлорид |
|
Этилен |
Метилбутилкетон |
Нафталин |
|
Эфиры: |
Метилпропилкетон |
Окись углерода |
|
диметиловый |
Метилэтилкетон |
Фенолхлорбензол |
|
дивиниловый |
Октан |
Эгиленбензол |
|
метилбутиловый |
Пиридин |
|
|
ШФЛУ |
Сероводород |
|
|
|
Спирты: |
|
|
|
метиловый |
|
|
|
этиловый |
|
|
|
пропиловый |
|
|
|
амиловый |
|
|
|
изобутиловый |
|
|
|
изопропиловый |
|
|
|
Циклогексан |
|
|
|
Этилформиат |
|
|
|
Этилхлорид |
|
|
|
LNG |
|
Классификация горючих смесей по масштабам фугасного (ударно-волнового) поражения
При оценке масштабов фугасного (ударно-волнового) поражения необходимо учитывать различие химических соединений по теплоте горения, используемой для расчета полного запаса энерговыделения. В расчетах энерговыделения не учитывается теплота конденсации водяных паров. Для типичных углеводородов берется в расчет значение удельной теплоты сгорания Q0 = 44 МДж • кг-1. Для иных горючих в расчетах используется удельное энерговыделение Q = βQ0. Здесь β - корректировочный параметр. Для условно выделенных классов опасных веществ величины параметра β представлены в табл. 3.2.
Классы опасных веществ |
β |
Классы опасных веществ |
β |
Класс 1 |
Класс 3 |
||
Ацетилен |
1,1 |
Этиламин |
0,80 |
Метилацетилен |
1,05 |
Метиламин |
0,70 |
Винилацетилен |
1,03 |
Спирты: |
|
Окись этилена |
0,62 |
метиловый |
0,45 |
Гидразин |
0,44 |
этиловый |
0,61 |
Изопропилнитрат |
0,41 |
пропиловый |
0,69 |
Этилнитрат |
0,30 |
амиловый |
0,79 |
Водород |
2,73 |
Циклогексан |
1 |
Нитрометан |
0,25 |
Ацетальальдегид |
0,56 |
Класс 2 |
Винилацетат |
0,51 |
|
Этилен |
1,07 |
Бензин |
1 |
Диэтилэфир |
0,77 |
Гексан |
1 |
Дивинилэфир |
0,77 |
Изооктан |
1 |
Окись пропилена |
0,7 |
СПГ |
1 |
Акролеин |
0,62 |
Пиродин |
0,77 |
Сероуглерод |
0,32 |
Циклопропан |
1 |
Бутан |
1 |
Класс 4 |
|
Бутилен |
1 |
Метан |
1,14 |
Бутадиен |
1 |
Трихлорэтан |
0,15 |
1,3-Пентадиен |
1 |
Метилхлорид |
0,12 |
Этан |
1 |
Бензол |
1 |
Диметилэфир |
0,66 |
Декан |
1 |
Диизопропилэфир |
0,82 |
Дизтопливо |
1 |
ШФЛУ |
1 |
Додекан |
1 |
Пропилен |
1 |
Керосин |
1 |
Пропан |
1 |
Метилбензол |
1 |
Класс 3 |
Метилмеркаптан |
0,53 |
|
Печной газ |
0,09 |
Нафталин |
0,91 |
Винилхлорид |
0,42 |
Окись углерода |
0,23 |
Сероводород |
0,34 |
Дихлорэтан |
0,24 |
Ацетон |
0,65 |
Дихлорбензол |
0,42 |
Кумол |
0,84 |
Трихлорэтан |
0,14 |
Классификация окружающей территории
В связи с тем, что характер окружающего пространства в значительной степени определяет скорость взрывного превращения облака топливовоздушной смеси и, следовательно, параметры ударной волны, геометрические характеристики окружающего пространства также разделены на четыре класса в соответствии со степенью их опасности.
Класс I. Наличие длинных труб, полостей, каверн, заполненных горючей смесью, при сгорании которой возможно ожидать формирование турбулентных струй продуктов сгорания размером не менее трех размеров детонационной ячейки данной смеси. Если размер детонационной ячейки для данной смеси не известен, то минимальный характерный размер струй принимается равным 5 см для веществ класса 1; 20 см - для веществ класса 2; 50 см - для веществ класса 3 и 150 см - для веществ класса 4.
Класс II. иСильнозагроможденное пространство: наличие полузамкнутых объемов, высокая плотность размещения технологического оборудования, лес, большое количество повторяющихся препятствий.
Класс III. Среднезагроможденное пространство: отдельно стоящие технологические установки, резервуарный парк.
Класс IV. Слабозагроможденное и свободное пространство.
Классификация взрывного режима
Для оценки действия взрыва возможные взрывные режимы превращения топливовоздушной смеси разбиты на шесть классов по диапазонам скоростей их распространения.
Класс 1. Детонация или горение со скоростью фронта пламени 500 м • с-1.
Класс 2. Дефлаграция, скорость фронта пламени 300-500 м • с-1.
Класс 3. Дефлаграция, скорость фронта пламени 200-300 м • с-1.
Класс 4. Дефлаграция, скорость фронта пламени 150-200 м • с-1.
Класс 5. Дефлаграция, скорость фронта пламени определяется соотношением
(3.36)
где k1 - константа, лежащая в диапазоне от 35 до 43;
М - масса топлива, содержащегося в облаке горючей смеси, кг.
Класс 6. Дефлаграция, скорость фронта пламени определяется соотношением
(3.37)
где k2 - константа, лежащая в диапазоне от 17 до 26;
М - масса топлива, содержащегося в облаке горючей смеси, кг.
Ожидаемый режим взрывного превращения определяется с помощью экспертной табл. 3.3, в зависимости от класса топлива и класса окружающего пространства.
Класс смеси |
Класс окружающего пространства |
|||
I |
II |
III |
IV |
|
1 |
1 |
1 |
2 |
3 |
2 |
1 |
2 |
3 |
4 |
3 |
2 |
3 |
4 |
5 |
4 |
3 |
4 |
5 |
6 |
При определении максимальной скорости фронта пламени для взрывных процессов 2-4 классов дополнительно рассчитывается видимая скорость фронта пламени по соотношению (3.36). В том случае, если полученная величина больше максимальной скорости, соответствующей данному классу, она принимается за верхнюю границу диапазона ожидаемых скоростей взрывного процесса в топливовоздушной смеси.
Для дальнейших расчетов необходимо оценить агрегатное состояние топлива смеси. Предполагается, что смесь гетерогенная, если более 50 % топлива содержится в облаке в виде капель. Провести такие оценки можно исходя из величины давления насыщенных паров топлива при данной температуре и времени формирования облака. Для летучих веществ, таких как пропан, при температуре +20 °С смесь можно считать газовой, а для веществ с низким давлением насыщенного пара (распыл дизтоплива при +20 °С) расчеты проводятся в предположении гетерогенной топливовоздушной смеси.
3.4.2. Расчет максимального избыточного давления и импульса фазы сжатия воздушных ударных волн
После того, как определен вероятный режим взрывного превращения, рассчитываются параметры воздушных ударных волн (избыточное давление ΔР и импульс фазы сжатия I+) в зависимости от расстояния от центра облака.
Детонация газовых и гетерогенных ТВС
Для вычисления параметров воздушной ударной волны на заданном расстоянии от центра облака, при детонации облака ТВС, предварительно рассчитывается соответствующий безразмерный радиус по соотношению
Rx = R/(E/P0)1/3 (3.38)
где R - расстояние от центра облака;
Р0 - атмосферное давление;
Е - эффективный энергозапас смеси.
Если соотношения записаны в функции аргумента λ = R0/E1/3, величина Е задается в МДж, то λ = 2,15 Rx. После этого рассчитываются величины безразмерного давления (Рх) и импульс фазы сжатия (Ix).
В случае детонации газовой смеси расчет проводится по формулам:
1n(Рх) = -1,124 - 1,66(1n(Rx)) + 0,260 (ln(Rx)2; (3.39)
1n(Iх) = -3,4217 - 0,898(1n(Rx)) - 0,0096 (ln(Rx))2. (3.40)
Зависимости (3.39, 3.40) справедливы для значений Rx больших величины Rk = 0,2, в случае, если Rx < Rk, то Рх полагается равным 18, а в выражение (3.40) вместо Rx подставляется величина Rxx = 0,14.
В случае детонации облака гетерогенной топливовоздушной смеси
(3.41)
Iх = 0,022/ Rx (3.42)
или
Рх = 0,271/λ + 0,62/λ2 + 0,236/λ3; (3.41')
Ix= 0,047/λ. (3.42')
Зависимости (3.41, 3.42) справедливы для значений Rx больших величины Rk = 0,2, в случае, если Rx < Rk, то Рх полагается равным 18, а в выражение (3.42) вместо Rx подставляется величина Rxx = 0,134.
После определения безразмерных величин давления и импульса фазы сжатия, вычисляются соответствующие им размерные величины:
(3.43)
(3.44)
Дефлаграция газовых и гетерогенных ТВС
В случае дефлаграционного взрыва облака ТВС, к параметрам, влияющим на величины избыточного давления и импульса положительной фазы, добавляются величины скорости видимого фронта пламени (u) к степени расширения продуктов сгорания (σ). Степень расширения продуктов сгорания дли газовых смесей принимается σ = 7, для гетерогенных - σ = 4. Для расчета параметров ударной волны при дефлаграционном взрыве гетерогенных облаков эффективная величина энергозапаса смеси домножается на коэффициент χ = (σ - 1)/ σ.
Как и для случая детонации, предварительно рассчитывается соответствующее безразмерное расстояние Rx от центра облака (3.38).
Затем рассчитываются величины безразмерного давления (Рх1) и импульса фазы сжатия (Ix1).
(3.45)
(3.46)
где
Выражения (3.45, 3.46) справедливы для значений Rx больших величины Rкp1 = 0,34, в случае, если Rx < Rкp1, выражения (3.45, 3.46) вместо Rx подставляется величина Rкp1. Далее вычисляются величины (Рх2) и (Ix2), которые для случая детонации газовой смеси рассчитываются по соотношениям (3.39, 3.40), а для детонации гетерогенной смеси - по соотношениям (3.41, 3.42). Окончательные значения Рх и Ix выбираются из условия
Px = min(Px1, Px2);
Ix = min (Ix1, Ix2). (3.47)
3.5. Параметры волны давления при взрыве резервуара с перегретой жидкостью или сжиженным газом при воздействии на него очага пожара
Избыточное давление ΔР и импульс i в ударной волне, образующиеся при взрыве резервуара с перегретой ЛВЖ, ГЖ или СУГ в очаге пожара, определяются по формулам:
(3.48)
(3.49)
где приведенная масса, кг; (3.50)
r - расстояние от центра резервуара, м;
H0 = 4,52·106 Дж·кг-1;
Eeff - эффективная энергия взрыва, рассчитываемая по формуле
Eeff = kCpm(T-Tb), (3.51)
где k - доля энергии ударной волны (допускается принимать равной 0,5);
Ср - удельная теплоемкость жидкости (допускается принимать равной 2 000 Дж • кг-1 • К-1;
m - масса ЛВЖ, ГЖ или СУГ, содержащаяся в резервуаре, кг;
Т - температура жидкой фазы, К;
Тb - нормальная температура кипения, К.
При наличии в резервуаре предохранительного устройства (клапана или мембраны) величина Т определяется по формуле
(3.52)
где Pval - давление срабатывания предохранительного устройства; А, В, СA - константы уравнения зависимости давления насыщенных паров топлива от температуры (константы Антуана), определяемые по справочной литературе. Единицы измерения Pval (кПа, мм рт. ст., атм) должны соответствовать используемым константам Антуана.
3.6. Интенсивность теплового излучения
В настоящем разделе приводятся методы расчета интенсивности теплового излучения от пожара пролива и огненного шара, а также радиуса воздействия продуктов сгорания паровоздушного облака в случае пожара-вспышки.
3.6.1. Пожар пролива
Интенсивность теплового излучения q, кВт • м-2, для пожара пролива ЛВЖ, ГЖ или СУГ вычисляется по формуле
где Ef - среднеповерхностная интенсивность теплового излучения пламени, кВт • м-2;
Fq - угловой коэффициент облученности;
τ - коэффициент пропускания атмосферы.
Значение Ef принимается на основе имеющихся экспериментальных данных или по табл. 3.4. При отсутствии данных для нефтепродуктов допускается принимать величину Ef равной 40 кВт • м-2.
Среднеповерхностная плотность теплового излучения пламени в зависимости от диаметра очага и удельная массовая скорость выгорания для некоторых жидких углеводородных топлив
Топливо |
Еf, кВт • м-2, при d, м |
m, кг-м-2 • с-1 |
||||
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 |
Примечание. Для диаметров очага менее 10 или более 50 м следует принимать Ef такой же, как и для очагов диаметром 10 и 50 м соответственно.
Угловой коэффициент облученности Fq определяется по формуле
(3.54)
где FV, FH - факторы облученности для вертикальной и горизонтальной площадок соответственно, определяемые с помощью выражений:
(3.55)
(3.56)
(3.57)
(3.58)
S = 2r/d; (3.59)
h = 2H/d, (3.60)
где r - расстояние от геометрического центра пролива до облучаемого объекта, м;
d - эффективный диаметр пролива, м;
Н - высота пламени, м.
Эффективный диаметр пролива d, м, рассчитывается по формуле
(3.61)
где F - площадь пролива, м2.
Высота пламени Н, м, вычисляется по формуле
(3.62)
где m - удельная массовая скорость выгорания топлива, кг • м-2 • с-1);
ρа - плотность окружающего воздуха, кг • м-3;
g - ускорение свободного падения, принимаемое равным 9,81 м • с-2.
Коэффициент пропускания атмосферы τ для пожара пролива определяется по формуле
τ = ехр[-7·10-4(r - 0,5d)]. (3.63)
3.6.2. Огненный шар
Интенсивность теплового излучения q(кВт • м-2) для огненного шара вычисляется по формуле (3.53).
Величина Еf определяется на основе имеющихся экспериментальных данных. Допускается принимать Еf равной 450 кВт • м-2.
Значение Fq определяется по формуле
(3.64)
где Н - высота центра огненного шара, м;
DS - эффективный диаметр огненного шара, м;
r - расстояние от облучаемого объекта до точки на поверхности земли непосредственно под центром огненного шара, м.
Эффективный диаметр огненного шара DS определяется по формуле
DS = 5,33m0,327, (3.65)
где m - масса продукта, поступившего в окружающее пространство, кг.
Величину H допускается принимать равной DS/2.
Время существования огненного шара tS, с, определяется по формуле
tS = 0,92m0,303. (3.66)
Коэффициент пропускания атмосферы τ для огненного шара рассчитывается по формуле
(3.67)
3.6.3. Определение радиуса воздействия продуктов сгорания паровоздушного облака в случае пожара-вспышки
В случае образования паровоздушной смеси в незагроможденном технологическим оборудованием пространстве и его зажигании относительно слабым источником (например, искрой) сгорание этой смеси происходит, как правило, с небольшими видимыми скоростями пламени. При этом амплитуды волны давления малы и могут не приниматься во внимание при оценке поражающего воздействия. В этом случае реализуется так называемый пожар-вспышка, при котором зона поражения высокотемпературными продуктами сгорания паровоздушной смеси практически совпадает с максимальным размером облака продуктов сгорания (т.е. поражаются в основном объекты, попадающие в это облако). Радиус воздействия высокотемпературных продуктов сгорания паровоздушного облака при пожаре-вспышке RF определяется приближенным соотношением
(3.68)
где Еi - объемный коэффициент расширения продуктов сгорания. Для нефтепродуктов величина Еi может быть принята равной 7, при этом = 1,82;
XLFL - горизонтальный размер взрывоопасной зоны, определяемый по п. 3.3.
3.7. Испарение жидкости и СУГ из пролива
3.7.1. Жидкость
Интенсивность испарения W (кг • м-2 • с-1) для ненагретых жидкостей с удовлетворительной точностью может быть описана полуэмпирическим выражением
(3.69)
где η - коэффициент, принимаемый по табл. 3.5 в зависимости от скорости и температуры воздушного потока над поверхностью испарения;
М - молярная масса жидкости, кг • кмоль-1;
РS - давление насыщенного пара при расчетной температуре жидкости, кПа.
Значения коэффициента η
Скорость воздушного потока, м • с-1 |
Значение коэффициента η при температуре t (°С) воздуха |
||||
10 |
15 |
20 |
30 |
35 |
|
0 |
1,0 |
1,0 |
1,0 |
1,0 |
1,0 |
0,1 |
3,0 |
2,6 |
2,4 |
1,8 |
1,6 |
0,2 |
4,6 |
3,8 |
3,5 |
2,4 |
2,3 |
0,5 |
6,6 |
5,7 |
5,4 |
3,6 |
3,2 |
1,0 |
10,0 |
8,7 |
7,7 |
5,6 |
4,6 |
3.7.2. Сжиженный углеводородный газ
При выбросе СУГ из поврежденного оборудования, в котором жидкость находится под давлением, часть продукта за счет внутренней энергии мгновенно испаряется, образуя с капельками жидкости облако аэрозоля. Массовую долю мгновенно испарившейся жидкости 8 определяют из соотношения
(3.70)
где Ср - удельная теплоемкость СУГ, Дж • кг-1 • К-1;
Та - температура окружающего воздуха, К;
Tg - температура кипения СУГ при атмосферном давлении, К;
Lg - удельная теплота парообразования СУГ, Дж • кг-1.
Принимается, что при δ ≥ 0,35 вся масса жидкости, находящаяся в оборудовании, за счет взрывного характера испарения переходит в парокапельное облако.
При δ < 0,35 оставшаяся часть жидкости испаряется с поверхности пролива за счет потока тепла от подстилающей поверхности и воздуха.
Интенсивность испарения жидкости со свободной поверхности W (кг • м-2 • с-1) описывается выражением
(3.71)
где λS - коэффициент теплопроводности материала, на поверхность которого разливается жидкость, Вт • м-1 • К-1;
CS - удельная теплоемкость материала, Дж • кг-1 • К-1);
ρS - плотность материала, кг • м-3;
T0 - начальная температура материала, К;
t - текущее время с момента начала испарения, с (но не менее 10 с);
λа - коэффициент теплопроводности воздуха при температуре T0;
u - скорость воздушного потока над поверхностью испарения, м • с-1;
d - характерный диаметр пролива, м;
va - кинематическая вязкость воздуха при T0, м2- с-1.
3.8. Размеры факела при струйном горении
3.8.1. Истечение паровой фазы СУГ
Расчет геометрических параметров факелов при истечении паровой фазы СУГ проводится на основе схемы, показанной на рис. 3.4.
Рис. 3.4. Схема для расчета геометрических параметров факелов при истечении паровой фазы СУГ
Приведенный метод применим, когда скорость ветра больше 5 м • с-1.
При меньших значениях скорости ветра этот угол следует считать равным нулю.
Геометрические параметры газовых факелов в случае истечения паровой фазы СУГ вычисляются по формулам (для скорости ветра не менее 5 м • с-1):
(3.72)
(3.73)
(3.74)
(3.75)
(3.76)
(3.77)
(3.78)
(3.79)
где Ua - скорость ветра, м • с-1;
Uj - скорость истечения паровой фазы СУГ, принимаемая равной скорости звука, м • с-1;
d0 - диаметр отверстия истечения, м;
ρj - плотность истекающей паровой фазы за срезом сопла (при атмосферном давлении), кг • м-3;
ρа - плотность воздуха, кг • м-3;
γ - показатель адиабаты истекающей паровой фазы (допускается принимать равным 1,2);
Р0 - атмосферное давление, Па;
Θ - угол наклона оси факела к вертикали, град;
Мj - молярная масса истекающей паровой фазы, кг • кмоль-1;
Ma - молярная масса воздуха, кг • кмоль-1;
LBV - высота центра верхнего основания факела над горизонтом, м;
W1, W2 - диаметры нижнего и верхнего оснований факела, м.
3.8.2. Истечение жидкой фазы СУГ
Длина и диаметр факела, образующегося при истечении жидкой фазы СУГ, описываются формулами:
(3.80)
(3.81)
(3.82)
(3.83)
где D1 - эффективный диаметр отверстия истечения, м;
ρL - плотность жидкой фазы, кг • м-3;
G1 - расход жидкой фазы через отверствие истечения, кг • с-1;
F - площадь отверстия истечения жидкой фазы, м2;
LF - высота факела, м;
g - ускорение свободного падения, м • с-2;
В - диаметр факела, м.
3.8.3. Истечение сжатого газа
Длина LF и диаметр dF - факела углеводородных газов вычисляют по формулам:
LF = kG0,4; (3.84)
(3.85)
где LF, dF - длина и максимальный диаметр факела, м;
G - расход газа, кг • с-1;
k = 12,3 - коэффициент пропорциональности.
Размеры факела принимаются независимыми от направления истечения газа.
При истечении газа вдоль ограничивающей поверхности (земля, оборудование и т.д.) расчетную длину факела следует увеличить на 25 %.
3.9. Тепловое излучение от горящего резервуара
Оценку интенсивности теплового излучения пламени горящего резервуара (горение жидкости по всей площади поперечного сечения резервуара) можно получить по методу, изложенному в п. 3.6, принимая в качестве диаметра пролива диаметр резервуара.
Для аварий с пожарами и взрывами на промышленных объектах наиболее характерными поражающими факторами являются ударная волна, тепловое излучение и осколочное воздействие при взрывах технологических аппаратов.
Ниже рассмотрены различные подходы к определению критериев поражения людей, зданий и оборудования перечисленными поражающими факторами.
4.1. Критерии поражения ударной волной
Известные критерии поражения объектов ударной волной (УВ) можно условно разделить на детерминированные и вероятностные. Детерминированные критерии приближенно показывают значения параметров УВ, при которых наблюдается тот или иной уровень поражения (разрушения). Вероятностные критерии показывают, какова условная вероятность того или иного уровня поражения (разрушения) при заданном значении поражающего фактора (например, избыточного давления и импульса положительной фазы УВ). Дальше рассматриваются как те, так и другие критерии.
4.1.1. Детерминированные критерии оценки поражающего действия УВ
Для описания поражения различных объектов УВ принято использовать так называемые Р/i диаграммы. Для кратковременных импульсных воздействий зачастую приемлемым является использование только критических значений импульса i УВ, для относительно длительных воздействий, типичных для дефлаграционных взрывов паровоздушных облаков - критических значений избыточного давления ΔР положительной фазы УВ.
Воздействие УВ на конструкции во многом определяется величиной τ/T, где τ - время воздействия положительной фазы УВ, T - период собственных колебаний конструкции. При τ/T > 2,5 воздействие определяется величиной ΔР, при τ/T < 0,1 - импульсом УВ. Собственные частоты некоторых сооружений приведены в табл. 4.1.
Наименование зданий и сооружений |
Т, с |
Двухэтажное кирпичное здание |
0,14 |
Башенный копер шахты высотой 70 м |
0,20 |
Крупнопанельное здание |
0,22÷0,32 |
Пятиэтажный каркасный производственный корпус с железобетонными диафрагмами |
0,25 |
Трехэтажное кирпичное здание заводоуправления |
0,34 |
Силосный корпус элеватора высотой 25 м |
0,36 |
Металлическая водонапорная башня высотой 20 м |
0,40 |
Каркасное трехэтажное здание фабрики высотой 17,5 м и размерами в плане 114 × 37 м |
0,85 |
Кирпичная дымовая труба |
2,86 |
Значения критического давления ΔР для разрушения УВ тех или иных элементов зданий приведены в табл. 4.2, для повреждений некоторых промышленных конструкций - в табл. 4.3.
Характер повреждений элементов зданий |
ΔР, кПа |
Разрушение остекления |
2÷7 |
Разрушение перегородок и кровли: |
|
деревянных каркасных зданий |
12 |
кирпичных зданий |
15 |
железобетонных каркасных зданий |
17 |
Разрушение перекрытий: |
|
деревянных каркасных зданий |
17 |
промышленных кирпичных зданий |
28 |
промышленных зданий со стальным и железобетонным каркасом |
30 |
зданий с массивными стенами |
42 |
Разрушение стен: |
|
шлакоблочных зданий |
22 |
деревянных каркасных зданий |
28 |
кирпичных зданий со стенами в 1,5 кирпича |
40 |
зданий с массивными стенами |
100 |
Разрушение фундаментов |
215÷400 |
Характер повреждений промышленных конструкций |
ΔР, кПа |
Незначительное повреждение стальных конструкций каркасов, ферм |
8÷10 |
Разрушение стальных каркасов, ферм и перемещение оснований |
20 |
Разрушение промышленных стальных несущих конструкций |
20 ÷ 30 |
Разрушение опорных структур резервуаров |
100 |
Перемещение цилиндрических резервуаров, повреждение трубопроводов |
50÷100 |
Повреждение ректификационных колонн |
35 ÷ 80 |
Незначительные деформации трубопроводных эстакад |
20 ÷ 30 |
Перемещение трубопроводных эстакад, повреждение трубопроводов |
35 ÷ 40 |
Разрушение трубопроводных эстакад |
40 ÷ 55 |
Значения критического давления ΔР, при которых реализуется поражение органов слуха человека, приведены в табл. 4.4.
Эффект |
Избыточное давление ΔР, кПа |
Граница временной потери слуха |
2,0 |
Нижний порог разрыва барабанной перепонки |
34,5 |
50 %-я вероятность разрыва барабанной перепонки |
103 |
100 %-я вероятность разрыва барабанной перепонки |
400 |
Значения критического давления ΔР для людей, находящихся в зданиях, приведены в табл. 4.5.
Вид воздействия |
Давление воздействия, кПа |
Люди, находящиеся в неукрепленных зданиях, погибнут в результате прямого поражения УВ, под развалинами зданий или вследствие удара о твердые предметы |
190 |
Наиболее вероятно, что все люди, находящиеся в неукрепленных зданиях, либо погибнут, либо получат серьезные повреждения в результате действия взрывной волны, либо при обрушении здания или перемещении тела взрывной волной |
69÷76 |
Люди, находящиеся в неукрепленных зданиях, либо погибнут или получат серьезные повреждения барабанных перепонок и легких под действием взрывной волны, либо будут поражены осколками и развалинами здания |
55 |
Обслуживающий персонал получит серьезные повреждения с возможным летальным исходом в результате поражения осколками, развалинами здания, горящими предметами и т.п. Имеется 10%-я вероятность разрыва барабанных перепонок |
24 |
Возможна временная потеря слуха или травмы в результате вторичных эффектов взрывной волны, таких, как обрушение зданий, и третичного эффекта переноса тела. Летальный исход или серьезные повреждения от прямого воздействия взрывной волны маловероятны |
16 |
С высокой надежностью гарантируется отсутствие летального исхода или серьезных повреждений. Возможны травмы, связанные с разрушением стекол и повреждением стен здания |
5,9÷8,3 |
4.1.2. Вероятностные критерии оценки поражающего действия УВ
В качестве вероятностного критерия поражения людей и/или зданий и сооружений используется понятие пробит-функции (probit-function). В общем случае пробит-функция Рr описывается выражением
Pr = a + blnS, (4.1)
где а, b - константы, зависящие от степени поражения и вида объекта;
S - интенсивность воздействующего фактора.
Соотношения между величиной Рr и условной вероятностью поражения человека приведено в табл. 4.6.
Таблица 4.6
Значения пробит-функции
Условная вероятность поражения, % |
Величина Рr |
|||||||||
0 |
1 |
2 |
3 |
4 |
5 |
6 |
7 |
8 |
9 |
|
0 |
- |
2,67 |
2,95 |
3,12 |
3,25 |
3,36 |
3,45 |
3,52 |
3,59 |
3,66 |
10 |
3,72 |
3,77 |
3,82 |
3,87 |
3,92 |
3,96 |
4,01 |
4,05 |
4,08 |
4,12 |
20 |
4,16 |
4,19 |
4,23 |
4,26 |
4,29 |
4,33 |
4,36 |
4,39 |
4,42 |
4,45 |
30 |
4,48 |
4,50 |
4,53 |
4,56 |
4,59 |
4,61 |
4,64 |
4,67 |
4,69 |
4,72 |
40 |
4,75 |
4,77 |
4,80 |
4,82 |
4,85 |
4,87 |
4,90 |
4,92 |
4,95 |
4,97 |
50 |
5,00 |
5,03 |
5,05 |
5,08 |
5,10 |
5,13 |
5,15 |
5,18 |
5,20 |
5,23 |
60 |
5,25 |
5,28 |
5,31 |
5,33 |
5,36 |
5,39 |
5,41 |
5,44 |
5,47 |
5,50 |
70 |
5,52 |
5,55 |
5,58 |
5,61 |
5,64 |
5,67 |
5,71 |
5,74 |
5,77 |
5,81 |
80 |
5,84 |
5,88 |
5,92 |
5,95 |
5,99 |
6,04 |
6,08 |
6,13 |
6,18 |
6,23 |
90 |
6,28 |
6,34 |
6,41 |
6,48 |
6,55 |
6,64 |
6,75 |
6,88 |
7,05 |
7,33 |
99 |
7,33 |
7,37 |
7,41 |
7,46 |
7,51 |
7,58 |
7,65 |
7,75 |
7,88 |
8,09 |
Для воздействия УВ на человека выражения для пробит-функции имеют вид
Rr = 5,0 - 5,74lnS; (4.2)
(4.3)
(4.4)
(4.5)
где m - масса тела человека, кг;
PS - избыточное давление в УВ, Па;
i - импульс УВ, Па • с;
Р0 - атмосферное давление, Па.
Пробит-функции для разрушения зданий имеют вид:
для тяжелых разрушений
Рr = 5,0 - 0,26lnV; (4.6)
(4.7)
для полного разрушения
Рr = 5,0 - 0,22lnV; (4.8)
Для практических целей целесообразно использовать пробит-функцию для тяжелых разрушений зданий (выражение 4.6), поскольку наиболее вероятный сценарий гибели людей от взрыва - это поражение их частями разрушающегося здания, расположенного вблизи взрывоопасного объекта.
4.2. Критерии поражения тепловым излучением
Для теплового излучения могут быть использованы как детерминированные, так и вероятностные критерии поражения.
4.2.1. Детерминированные критерии оценки поражения тепловым излучением
При анализе воздействия теплового излучения следует различать случаи импульсного и длительного воздействия. В первом случае следует говорить о дозе излучения D (случай огненного шара), во втором - о критической интенсивности теплового излучения qCR (например, пожар пролива).
Величины qCR для воспламенения некоторых горючих материалов приведены ниже в табл. 4.7, для различных степеней поражения человека - в табл. 4.8.
Таблица 4.7
Материалы |
qCR, кВт • м-2 |
Древесина (сосна влажностью 12 %) . |
13,9 |
Древесно-стружечные плиты (плотностью 417 кг • м-3) |
8,3 |
Торф брикетный |
13,2 |
Торф кусковой |
9,8 |
Хлопок-волокно |
7,5 |
Слоистый пластик |
15,4 |
Стеклопластик |
15,3 |
Пергамин |
17,4 |
Резина |
14,8 |
Уголь |
35,0 |
Рулонная кровля |
17,4 |
Картон серый |
10,8 |
Декоративный бумажно-слоистый пластик |
19,0-24,0 |
Металлопласт |
24,0-27,0 |
Плита древесно-волокнистая |
13,0 |
Плита древесно-стружечная |
12,0 |
Плита древесно-стружечная с отделкой «Полиплен» |
12,0 |
Плита древесно-волокнистая с лакокрасочным покрытием под ценные породы дерева |
12,0-16,0 |
Винилискожа обивочная пониженной горючести |
30,0-32,0 |
Кожа искусственная |
17,9-20,0 |
Стеклопластик на полиэфирной основе |
14,0 |
Лакокрасочные покрытия |
25,0 |
Обои моющиеся ПВХ на бумажной основе |
12,0 |
Линолеум ПВХ |
10,0-12,0 |
Линолеум алкидный |
10,0 |
Линолеум ПВХ на тканевой основе |
6,0-12,0 |
Покрытие ковровое |
4,0-6,0 |
Сено, солома (при минимальной влажности до 8 %) |
7,0 |
Легковоспламеняющиеся, горючие и трудногорючие жидкости при температуре самовоспламенения, °С: |
|
300 |
12,1 |
350 |
15,5 |
400 |
19,9 |
500 и выше |
28,0 и выше |
Таблица 4.8
Степень поражения |
Интенсивность излучения, кВт • м-2 |
Без негативных последствий в течение неограниченного времени |
1,4 |
Безопасно для человека в брезентовой одежде |
4,2 |
Непереносимая боль через 20-30 с |
7,0 |
Ожог 1 степени через 15-20 с |
10,5 |
Ожог 2 степени через 30-40 с |
|
Непереносимая боль через 3-5 с |
|
Ожог 1 степени через 6-8 с |
20,0 |
Ожог 2 степени через 12-16 с |
|
Мгновенные болевые ощущения через 4 с |
|
Летальный исход с вероятностью 50 % при длительности воздействия около 10 с |
44,5 |
Детерминированные критерии поражения технологического оборудования и строительных конструкций определяются, как правило, из условия потери ими несущей способности.
Основным критерием потери несущей способности является достижение на необогреваемой поверхности (например, внутренней поверхности резервуара, вблизи которого находится очаг пожара) критической температуры. Некоторые данные о значениях критических температур строительных материалов приведены ниже в табл. 4.9.
Таблица 4.9
Материал |
Критическая температура, °С |
Сталь марок Ст3, Ст5 в стальных конструкциях; арматурные стали классов: A-I, A-II, A-IV, Am-IV, A-V, Am-V |
500 |
Арматурная сталь классов: |
|
А-III |
600 |
Am-VI, Am-VII, B-I, Bp-I |
450 |
Высокопрочная арматурная проволока классов В-II, Вр-11 и арматурные канаты класса К-7 |
400 |
Сталь марок Ст3, Ст5 при жестком закреплении стальных конструкций и для металлических каркасов |
300 |
Конструкции из алюминиевых сплавов |
250 |
Гранит |
700 |
Известняк |
800 |
Тяжелый бетон с заполнителем |
600 |
Тяжелый бетон с известняковым заполнителем |
675 |
Кирпич глиняный обыкновенный |
900 |
Силикатный кирпич |
700 |
Оценка поражающего действия теплового излучения может быть проведена путем решения уравнения теплопроводности для облучаемого объекта (например, резервуара) с определением момента достижения температурой конструкционных материалов критических значений, приведенных в табл. 4.9.
Для приближенной оценки поражающего действия теплового излучения можно использовать метод, основанный на оценке критических величин плотностей тепловых потоков, вызывающих нагрев конструкционных материалов до критических температур. При этом величины критических тепловых потоков могут определяться как экспериментальным, так и расчетным путем.
Имеются экспериментальные данные по минимальной плотности теплового потока qmin, при которой возможен прогрев труб до температур самовоспламенения типичных нефтепродуктов (табл. 4.10).
Температура °C |
qmin (кВт • м-2) при времени прогрева (мин) |
|||
1 |
3 |
5 |
более 15 |
|
250 |
21,3 |
13,2 |
11,9 |
8,9 |
300 |
27,1 |
19,2 |
16,9 |
12,0 |
350 |
34,9 |
25,8 |
22,7 |
15,5 |
400 |
- |
34,9 |
30,2 |
19,9 |
500 и более |
- |
- |
34,9 |
27,9 |
В литературе описан метод расчетно-аналитического определения противопожарных разрывов между резервуарами с СУГ и ЛВЖ. Предполагается, что горение резервуара с ЛВЖ происходит во всем обваловании резервуара. В качестве критической температуры несмоченной стенки резервуара с СУГ принимается величина 427 °С. Расчеты показали, что плотность падающего на стенки резервуара теплового потока не должна превышать 13,5 кВт • м-2 в отсутствии ветра и 22 кВт • м-2 при скорости ветра 10 м • с-1.
Исходя из вышеизложенного, могут быть сделаны следующие практические рекомендации для учета поражающего действия теплового излучения по отношению к технологическому оборудованию и резервуарам. В случае полного охвата пламенем стенок резервуара представляется достаточным решение тепловой задачи в обыкновенных дифференциальных уравнениях с использованием критических температур из табл. 4.10. При несимметричном нагреве возможно использование понятия критической плотности теплового потока qmin. Если решается задача о возможности переноса пламени от горящего резервуара к негорящему, может быть использована табл. 4.10 и формула (4.9). Если решается задача о разрушении несмоченной стенки резервуара, то величина qmin может быть принята приближенно равной 20 кВт • м-2.
4.2.2. Вероятностные критерии оценки поражения тепловым излучением
В качестве вероятностного критерия оценки поражения тепловым излучением целесообразно использовать описанное выше понятие пробит-функции.
Для смертельного поражения человека величина пробит-функции описывается следующим выражением
Рr = -14,9 + 2,56ln(D); (4.11)
D = tq4/3. (4.12)
Величина эффективного времени экспозиции t может быть вычислена по формулам:
для огненного шара
t = 0,92m0,303; (4.13)
для пожара пролива
(4.14)
где m - масса горючего вещества, участвующего в образовании огненного шара, кг;
t0 - характерное время, за которое человек обнаруживает пожар и принимает решение о своих дальнейших действиях (может быть принято равным 5 с);
х - расстояние от места расположения человека до безопасной зоны (зона, где интенсивность теплового излучения меньше 4 кВт • м-2);
u - средняя скорость движения человека к безопасной зоне (может быть принята 5 м • с-1).
Величина дозы излучения D в случае пожара-вспышки для смертельного поражения человека может быть вычислена следующим образом. Вначале определяется эффективный радиус продуктов сгорания R по формуле
R = Rex(Ei-1)1/3, (4.15)
где Rex - радиус взрывоопасной зоны, м, вычисляемый в соответствии с изложенной в п. 3.3 методикой;
Ei - коэффициент расширения продуктов сгорания (может быть принят равным 7).
Величина D ((кВт • м-2)4/3 • с) определяется по следующим формулам в зависимости от расстояния r от геометрического центра паровоздушного облака:
Расстояние |
D, (кВт • м-2)4/3 • с |
r < 1,1R |
D = 3 000 |
1,1R ≤ r < 1,2R |
2 300 < D <3 000 |
1,2R ≤ r < 1,3R |
1 000 < D < 2 300 |
1,3R ≤ r < 1,4R |
200 < D < 1 000 |
r ≥ 1,4R |
D = 0 |
Величины D для промежуточных значений r следует определять с помощью линейной интерполяции.
Оценку индивидуального и социального риска для зданий промышленных предприятий следует проводить в соответствии с ГОСТ Р 12.3.047-98 (Приложение III).
Оценку потенциального, индивидуального, социального и коллективного рисков для территории предприятия следует проводить по изложенным ниже методикам.
5.1. Потенциальный риск
Величина потенциального риска Р(а) (год-1) в определенной точке местности а, где расположено предприятие, определяется с помощью соотношения
(5.1)
где I - число сценариев развития аварии (ветвей логического дерева событий возникновения и развития аварии);
Qdi(a) - условная вероятность поражения человека в определенной точке местности а в результате реализации i-го сценария развития аварии, отвечающего определенному событию, инициирующему аварию;
Q(Ai) - частота реализации в течение года i-го сценария развития аварии, год-1.
Условные вероятности поражения человека Qdi(a) определяются по значениям пробит-функций.
Величина Р(а) определяется посредством наложения зон поражения опасными факторами с учетом частоты реализации каждого сценария развития аварии на карту местности с привязкой их к соответствующему событию, инициирующему аварию (элементу оборудования, технологической установке) и ориентированию зоны поражения в соответствии с метеорологическими условиями (для струйного горения, пожара-вспышки, образования и взрывного превращения газопаровоздушного облака). При расчете риска рассматриваются различные метеорологические условия с типичными направлениями ветров и ожидаемой частотой их возникновения.
Процедура расчета риска предусматривает рассмотрение различных аварийных ситуаций и определение зон поражения опасными факторами пожара и взрыва, и частот их реализации. Для удобства расчетов территорию местности делят на зоны, внутри которых величины Р(а) полагаются одинаковыми.
В необходимых случаях оценка условной вероятности поражения человека проводится с учетом совместного воздействия более чем одного опасного фактора (для ветвей со стадиями с условием перехода «И»). Так, например, для расчета условной вероятности поражения человека при реализации сценария развития аварии, связанного со взрывом резервуара с ЛВЖ под давлением, находящегося в очаге пожара, необходимо учитывать, кроме теплового излучения огненного шара, воздействие ударной волны и осколков.
Условная вероятность поражения человека Qdi(a) от совместного независимого воздействия несколькими опасными факторами в результате реализации i-го сценария развития аварии определяется следующим образом:
(5.2)
где h - число рассматриваемых опасных факторов аварии;
Qk - вероятность реализации k-го опасного фактора;
- условная вероятность поражения k-м опасным фактором.
Результаты расчетов потенциального риска отображаются на карте (ситуационном плане) предприятия и прилегающих районов в виде замкнутых линий равных значений (изолинии функции Р(а)).
Изолинии функции Р(а) называются контурами риска. Их физический смысл состоит в том, что они разделяют территорию предприятия (так же, как и местность вокруг предприятия) на области, в которых ожидаемая частота возникновения опасных факторов аварии, приводящих к гибели людей, заключена в определенных, указанных на рисунке, пределах.
Контуры риска не зависят от количества работающих на предприятии или их должностных обязанностей, а определяются исключительно используемой технологией и надежностью применяемого оборудования. Потенциальный риск используется как мера (критерий допустимости / недопустимости) уровня пожарной безопасности объекта.
5.2. Индивидуальный риск
Для любого работника предприятия существует возможность гибели при возникновении аварии. Потеря жизни в течение определенного периода времени (года) является случайным событием, зависящим от рода его профессиональной деятельности, в том числе от продолжительности нахождения работника в областях, отвечающих различным контурам риска при его перемещениях по промышленной площадке предприятия в течение рабочей смены.
Для целей управления безопасностью персонала используется количественная мера возникновения этого случайного события - частота поражения опасными факторами пожара (взрыва) определенного человека, называемую индивидуальным риском.
Таким образом, индивидуальный риск определяется как ожидаемая частота поражения определенного работника предприятия опасными факторами аварий в течение года.
Области, на которые разбита территория предприятия, обозначается
J, j = 1, ... J.
Для удобства описания расчетов работники предприятия нумеруются
m = 1, ..., М
Текущий номер работника m однозначно определяет наименование должности работника, его категорию и другие особенности его профессиональной деятельности, необходимой для оценки безопасности.
Величина индивидуального риска Rm (год-1) для работника предприятия m определяется с помощью соотношения
(5.3)
где - величина потенциального риска в j-й области территории предприятия, год-1;
qjm - доля времени, в течение которого работник предприятия т находится в j-й области территории предприятия.
Доля времени, в течение которого работник находится в определенной области территории предприятия, рассчитывается на основе решений по организации эксплуатации и технического обслуживания оборудования.
5.3. Социальный риск
Для анализа воздействия промышленных аварий на людей, а также для установления приемлемости того или иного уровня пожарной или промышленной безопасности используется понятие социального риска.
Социальный риск задается с помощью функции, значениями которой являются величины, определяющие, что в аварии с пожаром погибло не менее определенного количества человек.
Социальный риск S (год-1) определяется по формуле
(5.4)
где L - число сценариев развития аварии, для которых выполняется условие Ni ≥ N0;
Ni - ожидаемое число погибших в результате реализации i-го сценария развития аварии;
N0 - число погибших, для которого оценивают величину социального риска. В настоящем документе принимается N0 = 10.
Ожидаемое число погибших в результате реализации i-го сценария развития аварии можно оценить по следующей формуле
(5.5)
где J - количество областей, на которые разбита территория предприятия и территория, прилегающая к предприятию (j - номер области);
- условная вероятность поражения человека, находящегося в j-й области, опасными факторами пожара (взрыва) при реализации i-го сценария развития аварии;
nj - среднее число людей, находящихся в j-й области.
Социальный риск S является интегральной величиной. В то же время социальный риск можно рассматривать и как векторную величину, компоненты которой имеют размерность год-1. При этом результаты расчетов социального риска могут быть представлены в виде так называемых F/N диаграмм, где по горизонтальной оси откладывают N - количество погибших в результате реализации i-гo сценария развития аварии, а по вертикальной оси - F - частоту реализации сценария, при котором погибло не менее N человек. Подобные зависимости могут быть аппроксимированы кривой - графиком непрерывной функции F(N). При этом указанная выше величина S описывается следующим выражением
(5.6)
5.4. Коллективный риск
Для персонала предприятия в целом имеется ненулевая вероятность гибели части работников при возникновении аварии.
Количество погибших в течение определенного периода времени (года) является случайной величиной, зависящей от опасности производства, количества работающих и ряда других факторов.
Для целей управления безопасностью персонала используется математическое ожидание этой случайной величины. Указанную характеристику называют коллективным риском персонала от аварий.
Величина коллективного риска персонала С (чел • год-1) определяется с помощью соотношения
(5.7)
Связь индивидуального риска работника от аварий Rm и коллективного риска персонала от аварий устанавливает соотношение
(5.8)
1. Corder I. The application of risk techniques to the design and operation of pipelines / Institution of Mechanical Engineers. Conference C502. Paper C502/016, 1995. - P. 56-69.
2. OREDA (Offshore Reliability Data handbook). - DnV, 1992. - P. 156-267.
3. Hydrocarbon Leak and Ignition Database: Report No. 11/4/180//E&P Forum, 1992. - P. 134-157.
4. Gas Pipeline Incident 1970-1997 EGIG: Third report. Document No. 98.R.0120 / EGIG (European Gas Pipeline Incident Data Group), 1998. - December. - 167 p.
5. Offtakes and Pressure-Regulating Installations for Inlet Pressures between 7 and 100 bar: Recommendations on Transmission and Distribution Practice: IGE/TD/9. Communication 1229 / Institution of Gas Engineers, 1986. - 234 p.
6. McConnel R.A. The Use of Slam Shut Valves on LCA Plants / Process Safety Progress. - 1997. - Summer. - Vol. 16, № 2. - P. 79-97.
7. Оценка риска аварий на линейной части магистральных нефтепроводов / М.В. Лисанов, А.С. Печеркин, В.И. Сидоров и др. - М: НТЦ «Промышленная безопасность», 2001. - 187 с.
8. ГОСТ Р 12.3.047-98. Пожарная безопасность технологических процессов. Общие требования. Методы контроля.
9. CPR 14E. Methods for the Calculation of Physical Effects (Yellow Book). Committee for the prevention of Disaster. - TheHague,1997.-836p.
10. Sallet D.W. Critical two-phase mass flow rates of liquefied gases // Journal of Loss Prevention in the Process Industries. - 1990. - Vol. 3, № 1. - P. 38-42.
11. Обеспечение пожарной безопасности объектов хранения и переработки СУГ: Рекомендации. - М: ВНИИПО, 1999. - 156 с.
12. Greenspan N.P., Young R.E. Flow over a containment dyke // Journal of Fluid Mechanics. - 1987. - Vol. 87, № 1. - P. 179-192.
13. Шебеко Ю.Н., Шевчук А.П., Смолин И.M. Расчет влияния обвалования на растекание горючей жидкости при разрушении резервуара // Химическая промышленность. - 1994. - № 4. - С. 230-233.
14. Андерсон Д., Таинненхил Дж., Плетчер Р. Вычислительная гидромеханика и теплообмен. - М.: Мир, 1990. - Т. 1.-179 с.
15. НПБ 105-03. Определение категорий помещений, зданий и наружных установок по взрывопожарной и пожарной опасности.
16. Пожаровзрывоопасность объектов хранения сжиженного природного газа. Процессы испарения и формирования пожаровзрывоопасных облаков при проливе жидкого метана. Методики оценки параметров / И.А. Болодъян, В.П. Молчанов, Ю.И. Дешевых и др. // Пожарная безопасность. - 2000. - № 4. - С. 108-121.
17. Пожарная опасность шаровых резервуаров для хранения под давлением сжиженных углеводородных газов / А.П. Шевчук, Ю.Н Шебеко, И.М. Смолин и др. // Химическая промышленность. - 1992. - № 6. - С. 328-332.
18. Борисов А.А., Гельфанд Б.Е., Цыганов С.А. О моделировании ударных волн давления, образующихся при детонации и горении газовых смесей // Физика горения и взрыва. - 1985. -Т. 21, № 2. - С. 90-97.
19. Методика оценки последствий аварийных взрывов топливно-воздушных смесей. - М.: НТЦ «Промышленная безопасность». - 45 с.
20. Шебеко Ю.Н., Шевчук А.П., Смолин И.М. Расчет параметров ударных волн, образующихся при взрыве резервуара со сжиженным углеводородным газом в очаге пожара // Химическая промышленность. - 1993. - № 9. - С. 451-453.
21. Mudan K.S. Thermal radiation hazards from hydrocarbon pool fires // Progress in Energy and Combustion Science. - 1984. - Vol. 10, № 1. - P. 59-80.
22. Leslie I.R.M., Birk A.M. State of the art review of pressure liquified gas container failure modes and associated projectile hazards // Journal of Hazardous Materials. - 1991. - Vol. 28, № 3. - P. 329-365.
23. Орлов В.Я., Никитин А.Г. Определение скорости испарения легковоспламеняющихся жидкостей при категорировании производств: Вопросы охраны труда при проектировании и строительстве зданий. - Казань: КХТИ, 1978. - С. 55-57.
24. Пчелинцев В.А., Никитин А.Г., Хузиахметов Р.А. Оценка взрывопожароопасности производств, связанных с применением легковоспламеняющихся жидкостей в нагретом состоянии // Журн. ВХО им. Д.И. Менделеева. - 1985. - Т. 30, № 1. - С. 68-74.
25. Маршалл В. Основные опасности химических производств. - М.: Мир, 1989. - 671 с.
26. Bais A.F., Zerefos C.S., Ziomas I.С. Design of a system for real-time modeling of the dispersion of hazardous gas releases in industrial plants. 1. Emissions from Industrial stacks // Journal of Loss Prevention in the Process Industries. - 1989. - July. - Vol. 2. - P. 155-160.
27. Ziomas I.С, Zerefos C.S., Bais A.F. Design of a system for real-time modeling of the dispersion of hazardous gas releases in industrial plants. 2. Accidental releases from storage installations // Journal of Loss Prevention in the Process Industries. - 1989. - October. - Vol. 2. - P. 194-199.
28. Моделирование распространения паровоздушного облака тяжелого газа при его мгновенном выбросе и непрерывном истечении / А.П. Шевчук, Ю.Н. Шебеко Л.В. Гуринович, И.М. Смолин // Химическая промышленность. - 1992. - № 10. - С. 622-625.
29. Some aspects of fire and explosion hazards of large LPG storage vessels / Yu.N. Shebeko, I.M. Smolin, A.Ya. Korolchenko et al. // Journal of Loss Prevention in the Process Industries. - 1995. - Vol. 8, № 3. - P. 163-168.
30. Van Ulden A.P. On the spreading of a heavy gas released near the ground // In: First International Symposium on Loss Prevention and Safety Promotion in the Process Industries. - 1974. - P. 221-226.
31. Попов П.С, Реутт В.Ч., Гришин В.В. Взрывоопасные зоны и горение компактных струй сжиженного газа при аварийном истечении // Экспресс-информация ВНИИПО. Серия: Пожарная профилактика в технологических процессах и строительстве. - М: ВНИИПО, 1974. - Вып. 36. - С. 35-47.
32. Мольков В.В., Некрасов В.П. Динамика сгорания газа в постоянном объеме при наличии истечения // Физика горения и взрыва. - 1981. - Т. 17, № 4. - С. 17-24.
33. Термогазодинамика пожаров в помещениях / В.М. Астапенко, Ю.А. Кошмаров, И.С. Молчадский, А.Н. Шевляков. - М: Стройиздат, 1988. - 448 с.
34. Методы расчета температурного режима пожара в помещениях зданий различного назначения: Рекомендации. - М.: ВНИИПО, 1988. - 53 с.
35. Fire and explosion risk assessment for LPG storages / Yu.N. Shebeko, A.Ya. Korolchenko, A.P. Shevchuk et al. // Fire Science and Technology. - 1995. - Vol. 15, № 1-2. - P. 37-45.
36 Chamberlain G.A. Developments in design methods for predicting thermal radiation from flares // Chem. Eng. Des. - 1987. - July. - Vol. 65. - P. 123-136.
37. Broechnann В., Schecker Н.G. Heat transfer mechanisms and boilover in burning oil-water systems // Journal of Loss Prevention in the Process Industries. - 1995. - Vol. 8, № 3. - P. 137-147.
38. Ramskill P.K. A description of the «ENGULF» computer codes-codes to model the thermal response of an LPG tank either fully or partially engulfed by fire // Journal of Hazardous Materials. - 1988. - Vol. 20, № 1-3. - P. 177-196.
39. Thermal response analysis of LPG tanks exposed to fire / N.U. Aydemir, V.K. Magapu, A.C.M. Sousa, J.E.S. Venart // Journal of Hazardous Materials. - 1988. - Vol. 20, № 1-3. - P.239-262.
40. Исаченко В.П., Осипова В.А., Сукомел А.С. Теплопередача. - М.: Энергия, 1975. - 486 с.
41. Взрывные явления: оценка и последствия / У. Бейкер, П. Кокс, П. Уэстайн и др. - Т. 1. - М.: Мир, 1986. - 319 с.
42. Взрывные явления: оценка и последствия / У. Бейкер, П. Кокс, П. Уэстайн и др. - Т. 2. - М.: Мир, 1986. - 384 с.
43. Clancey V.J. The effects of explosions // Institute of Chemical Engineering Symposium Series. - 1982. - № 71. - P. 87-108.
44.Кочегаров В.П. Нагрузки от ударных волн при детонации газовоздушных смесей: Охрана труда в строительстве. - М.: МИСИ, 1978. - С. 90-95.
45. CPR 16E. Methods for the determination of possible damage. Committee for the prevention of Disaster. - Voorburg,1989. - 326p.
46. Pietersen C.M. Consequences of accidental releases of hazardous material // Journal of Loss Prevention in the Process Industries. - 1990. - Vol. 3,№ l. - P. 136-141.
47. An analysis of the risks arising from the transport of liquefied gases in Great Britain / G. Puredy, H.S. Campbell, G.C Grint, L.M. Smith // Journal of Hazardous Materials. - 1988. - Vol. 20, № 1-3. - P. 335-355.
48. Волков О.М. Пожарная безопасность резервуаров с нефтепродуктами. - М.: Недра, 1984. - 234 с.
49. Martinsen W.E., Johnson D.W., Millsap S.B. Determining Spasing by Radiant Heat Limits // Plant Operations Progress. - 1989. - Vol. 8, № 1. - P. 25-28.
50. Risk assessment for Installations where liquefied petroleum gas (LPG) is stored in bulk vessels above ground / G.A. Clay, R.D. Fitzpatric, N.W. Hurst et al. // Journal of Hazardous Materials. - 1988. - Vol. 20, № 1-3. - P. 357-374.