Примечания:
1. Здесь и далее под полным разрушением подразумевается утечка с диаметром истечения, соответствующим максимальному диаметру подводящего или отводящего трубопровода, или разрушения резервуара, емкости, сосуда или аппарата.
2. При определении частоты разгерметизации фильтров и кожухотрубных теплообменников указанное оборудование можно рассматривать как аппараты под давлением.
3. Аппараты воздушного охлаждения можно рассматривать как участки технологических трубопроводов, длина которых соответствует суммарной длине труб в пучках теплообменника.
Таблица 1.2
Частоты утечек из технологических трубопроводов
Таблица 1.3
Частота утечек из магистрального трубопровода
Таблица 1.4
Частоты возникновения пожара
для некоторых зданий и сооружений
Логическое дерево событий предназначено для графического отображения общего характера развития возможных аварийных ситуаций и аварий с отражением причинно-следственной взаимосвязи событий в зависимости от специфики опасности объекта оценки риска с учетом влияния на них имеющихся защитных мероприятий и является основой для оценки риска. Сценарий возникновения и развития аварийной ситуации и аварии на логическом дереве отражается в виде последовательности событий от исходного до конечного события (ветвь дерева событий).
Данный метод позволяет проследить развитие возможных аварийных ситуаций и аварий, возникающих вследствие реализации событий, инициирующих аварийную ситуацию. Анализ дерева событий представляет собой "осмысливаемый вперед" процесс, то есть процесс, при котором исследование развития аварийной ситуации начинается с исходного события с рассмотрением цепочки последующих событий, приводящих к возникновению аварии.
При построении логических деревьев событий учитываются следующие положения:
- выбирается аварийная ситуация, которая может повлечь за собой возникновение аварии с пожаром (взрывом) с дальнейшим его развитием;
- развитие аварийной ситуации и аварии должно рассматриваться постадийно с учетом места ее возникновения на объекте оценки риска, уровня потенциальной опасности каждой стадии и возможности ее локализации и ликвидации. На логическом дереве событий стадии развития аварийной ситуации и аварии могут отображаться в виде прямоугольников или других геометрических фигур с находящимися в них краткими названиями этих стадий;
- переход с рассматриваемой стадии на новую определяется возможностью либо локализации аварии на рассматриваемой стадии, либо развития аварии, связанного с вовлечением расположенных рядом технологического оборудования, помещений, сооружений и т.п. в результате влияния на них опасных факторов пожара (взрыва), возникших на рассматриваемой стадии. Вероятности переходов аварии со стадии на стадию одной ветви или с ветви на ветвь определяется, исходя из величин зон поражения опасными факторами пожара (взрыва), объемно-планировочных решений и конструктивных особенностей оборудования и сооружений промышленного предприятия. При этом каждой стадии иногда присваивается идентификационный номер, отражающий последовательность переходов со стадии на стадию;
- переход со стадии на стадию отображается в виде соединяющих линий со стрелками, указывающими направления развития пожароопасной ситуации и последующего пожара. При этом соединения стадий должны отражать вероятностный характер события с выполнением условий "И", "ИЛИ" и "И/ИЛИ" (условие "И" - взаимосвязь событий, условие "ИЛИ" - независимость событий, а условие "И/ИЛИ" - возможность частичной реализации либо условия "И", либо "ИЛИ" ввиду многообразия путей дальнейшего развития). Отображение взаимосвязи событий по условию "И" выполняется сплошной линией без ответвлений, по условию "ИЛИ" - сплошной линией с ответвлениями к взаимоисключающим друг друга событиям (стадиям), по условию "И/ИЛИ" - пунктирной линией с ответвлениями к потенциально возможным стадиям;
- для каждой стадии должен устанавливаться уровень ее опасности, характеризующийся возможностью перехода аварийной ситуации или аварии на соседние с пожароопасным участки объекта;
- при повторении одним из путей части другого пути развития аварии для упрощения построения логического дерева событий иногда вводят обозначение, представляющее собой соответствующую линию со стрелкой и надпись "на стадию (код последующей стадии)". При этом в случае перехода со стадии с более высоким уровнем опасности на стадию с меньшим уровнем следует дополнять код стадии с меньшим уровнем опасности кодом, соответствующим стадии с более высоким уровнем опасности, указываемым в скобках.
При анализе логических деревьев событий руководствуются следующими положениями:
- возможность предотвращения дальнейшего развития аварийной ситуации и аварии зависит от количества стадий/времени их протекания (т.е. от длины пути развития аварийной ситуации и аварии). Это объясняется большей вероятностью успешной ликвидации аварийной ситуации и аварии, связанной с увеличением времени на локализацию аварийной ситуации и аварии и количеством стадий, на которых эта локализация возможна;
- наличие у стадии нескольких разветвлений по принципу "И" или "И/ИЛИ" свидетельствует о приобретении аварийной ситуацией и аварией "цепного" характера развития, т.е. одновременного развития аварии по нескольким путям с дальнейшим увеличением их количества (так называемый "эффект домино"), что в значительной мере затрудняет успешную локализацию и ликвидацию пожара;
- наличие у стадии разветвлений по принципу "ИЛИ", одно из которых приходит на стадию локализации аварийной ситуации или аварии (например, тушение очага пожара, своевременное обнаружение утечки и ликвидация пролива, перекрытие запорной арматуры и т.п.), свидетельствует о возможности предотвращения дальнейшего развития аварийной ситуации и аварии по этому пути. Отсутствие стадии локализации аварийной ситуации и аварии на разветвлениях по принципу "ИЛИ" свидетельствует о невозможности приостановления дальнейшего развития аварийной ситуации и аварии на разветвляемой стадии;
- чем больше возможных путей развития аварийной ситуации и аварии приходит на одну стадию, тем больше вероятность возникновения этой стадии.
Для оценки вероятности перехода аварии со стадии на стадию, прежде всего, необходимо определение условной вероятности реализации различных ветвей дерева событий. Далее необходимо определение вероятностей срабатывания соответствующих средств предотвращения или локализации аварии. Кроме того, необходима оценка вероятности поражения расположенного в зоне возникновения аварии технологического оборудования и сооружений объекта в результате воздействия на них опасных факторов пожара.
В табл. 2.1 приведены условные вероятности мгновенного воспламенения и воспламенения с задержкой по времени в зависимости от массовой скорости истечения газа, двухфазной среды или жидкости при разгерметизации типового технологического оборудования на промышленных предприятиях.
Таблица 2.1
Условная вероятность мгновенного воспламенения
и воспламенения с задержкой
Условные вероятности воспламенения при различных диаметрах отверстия истечения при разгерметизации технологического оборудования на промышленных предприятиях приведены в табл. 2.2 и 2.3. При этом под задержкой воспламенения следует понимать время более одной минуты.
Таблица 2.2
Условная вероятность воспламенения при различных диаметрах
отверстия истечения
Таблица 2.3
и воспламенения с задержкой при различных диаметрах
отверстия истечения
Для особо опасных легковоспламеняющихся жидкостей (ЛВЖ) (температура вспышки
) следует использовать условные вероятности воспламенения как для двухфазной среды.В табл. 2.4 приведены условные вероятности различных сценариев развития аварий с выбросом из технологического оборудования сжиженных углеводородных газов при полном разрушении оборудования.
Таблица 2.4
Условная вероятность различных сценариев развития аварий
С ПОЖАРАМИ И ВЗРЫВАМИ
В настоящем приложении представлены методики оценки опасных факторов, реализующихся при различных сценариях аварий с пожарами на промышленных предприятиях с наличием горючих газов, легковоспламеняющихся и горючих жидкостей, основанные на современных отечественных и зарубежных теоретических и экспериментальных данных, нормативных документах и руководствах.
В случае, если при анализе сценариев аварий возникает необходимость учесть сценарии, методы оценки опасных факторов для которых в настоящем документе не изложены, или возникает необходимость в отдельных случаях применять методы, отличающиеся от описанных ниже, возможно использование иных методик при соответствующем обосновании необходимости их применения.
Рассматривается резервуар, находящийся в обваловании (рис. 3.1).
![]() Рис. 3.1. Схема для расчета истечения жидкости
из отверстия в резервуаре
Вводятся следующие допущения:
- истечение через отверстие однофазное;
- резервуар имеет постоянную площадь сечения по высоте;
- диаметр резервуара много больше размеров отверстия;
- размеры отверстия много больше толщины стенки;
- поверхность жидкости внутри резервуара горизонтальна;
- температура жидкости остается постоянной в течение времени истечения.
Массовый расход жидкости G
через отверстие во времени t (с) описывается соотношением , (3.1)где
, (3.2)где
g - ускорение свободного падения
;Зависимость высоты столба жидкости в резервуаре h (м) от времени t описывается формулой
. (3.3)Условия перелива струи жидкости (при
) через обвалование может быть приближенно записано следующим образом:где H - высота обвалования, м;
L - расстояние от стенки резервуара до обвалования, м.
Количество жидкости m (кг), перелившейся через обвалование за полное время истечения, описывается выражением
, (3.5)где
Величина
, (3.6)где a, b, c - параметры, описываемые формулами:
, , , м. (3.9)В случае, если жидкость в резервуаре находится под избыточным давлением
должна быть описана выражением . (3.10)Для определения количества жидкости, перелившейся через обвалование, и времени перелива следует проинтегрировать соответствующую систему уравнений, где величина
Массовая скорость истечения сжатого газа из резервуара описывается следующими соотношениями:
докритическое истечение
; (3.12)сверхкритическое истечение
(3.13)где G - массовый расход,
Массовая скорость истечения паровой фазы
описывается формулой , (3.15)где
M - молярная масса,
;R - универсальная газовая постоянная, равная
; - безразмерное давление сжиженного газа в резервуаре;Массовую скорость истечения паровой фазы можно также рассчитывать по формулам (3.11) - (3.14).
Массовая скорость истечения жидкой фазы G
описывается формулой , (3.16)где
- безразмерная температура сжиженного газа;T - температура сжиженного газа в резервуаре, К.
разрушении резервуара
Под квазимгновенным разрушением резервуара следует понимать внезапный (в течение секунд или долей секунд) распад резервуара на приблизительно равные по размеру части.
При таких аварийных ситуациях часть хранимого в резервуаре продукта может перелиться через обвалование.
Ниже представлена математическая модель, позволяющая оценить долю жидкости, перелившейся через обвалование при квазимгновенном разрушении резервуара. Приняты следующие допущения:
- рассматривается плоская одномерная задача;
- время разрушения резервуара много меньше характерного времени движения гидродинамической волны до обвалования;
- жидкость является невязкой;
- трение жидкости о поверхность земли отсутствует;
- поверхность земли является плоской, горизонтальной.
Система уравнений, описывающих движение жидкости, имеет вид
где h - высота столба жидкости над фиксированным уровнем, м;
u - средняя по высоте скорость движения столба жидкости,
x - координата вдоль направления движения жидкости, м;
t - время, с;
g - ускорение свободного падения
.Граничные условия с учетом геометрии задачи (рис. 3.2) имеют вид
; (3.19) ; (3.20)где a - высота обвалования.
![]() при квазимгновенном разрушении резервуара:
- уровень начального столба жидкости; - уровень жидкости в промежуточный моментвремени (результаты расчета)
Массовая доля жидкости Q (%), перелившейся через обвалование к моменту времени T, описывается выражением
, (3.22)где
R - ширина резервуара, м.
Один из возможных способов решения системы (3.17) с граничными условиями (3.18) - (3.21) - численное решение с использованием схемы Мак-Кормака, представляющей собой частный случай явной схемы типа "предиктор-корректор".
График указанной функции вместе с литературными экспериментальными данными представлен на рис. 3.3. Для практических целей представляется достаточным использование этого графика.
![]() Рис. 3.3. Зависимость доли перелившейся через
обвалование жидкости Q от параметров
1 - расчет; 2 - эксперимент
поступающих в окружающее пространство в результате
возникновения аварийных ситуаций
Количество поступивших в замкнутое или свободное пространство веществ, которые могут образовать взрывоопасные паровоздушные смеси или проливы горючих жидкостей, определяется, исходя из следующих предпосылок:
а) происходит расчетная авария одного из резервуаров или трубопроводов;
б) все содержимое резервуара (трубопроводов) или часть продукта (при соответствующем обосновании) поступает в замкнутое или свободное пространство. При этом в случае наличия на объекте нескольких резервуаров (с различными объемами) расчет следует проводить для вариантов с разгерметизацией каждого единичного резервуара;
в) при разгерметизации резервуара происходит одновременно утечка веществ из трубопроводов, питающих резервуар по прямому и обратному потоку в течение времени, необходимого для отключения трубопроводов. Расчетное время отключения трубопроводов определяется в каждом конкретном случае, исходя из реальной обстановки, и должно быть минимальным с учетом паспортных данных на запорные устройства и их надежности, характера технологического процесса и вида расчетной аварии. Под "временем отключения" следует понимать промежуток времени от начала возможного поступления горючего вещества из трубопровода (перфорация, разрыв, изменение номинального давления и т.п.) до полного прекращения поступления жидкости в окружающее пространство;
г) при расчетах принимается нулевая подвижность окружающего воздуха (безветрие);
д) в качестве расчетной температуры при аварийной ситуации с наземным расположением оборудования принимается максимально возможная температура воздуха в соответствующей климатической зоне, а при аварийной ситуации с подземным расположением оборудования - температура грунта, условно равная максимальной среднемесячной температуре окружающего воздуха в наиболее теплое время года;
е) длительность испарения жидкости с поверхности пролива принимается равной времени ее полного испарения, но не более 3600 с. Для относительно небольших проливов топлива (до 20 кг) время испарения допускается принимать равным 900 с, поскольку столь небольшие проливы могут быть достаточно эффективно удалены обслуживающим персоналом. Кроме того, в запас надежности идет неучет подвижности воздуха и уменьшение скорости испарения жидкости со временем вследствие ее охлаждения.
При проведении расчетов допускается использование справочных данных, опубликованных головными научно-исследовательскими организациями в области пожарной безопасности или выданных Государственной службой стандартных справочных данных. Допускается использование показателей пожаровзрывоопасности для смесей веществ и материалов по наиболее опасному компоненту.
Ниже приводятся основные расчетные формулы для определения масс горючих веществ, поступающих в открытое или замкнутое пространство в результате аварийных ситуаций.
Масса жидкости, поступившей в окружающее пространство при разгерметизации резервуара, определяется по формуле
, (3.23)где
Масса жидкости, поступившей самотеком при полном разрушении трубопровода, выходящего из резервуара, определяется по формуле
; (3.24) ; (3.25) , (3.26)где
n - число участков трубопроводов, связанных с местом разгерметизации;
g - ускорение свободного падения
.через дыхательную арматуру
В случае наполнения резервуара массу паров вычисляют по формуле
; (3.27) , (3.28)где
M - молярная масса топлива,
; ;поверхности в резервуаре
Массу паров ЛВЖ при испарении со свободной поверхности в резервуаре вычисляют по формуле
, (3.29)где
; (3.30)W - интенсивность испарения ЛВЖ,
.Радиус
и высота зоны, ограничивающие область концентраций, превышающих нижний концентрационный предел распространения пламени (НКПР), при неподвижной воздушной среде рассчитывают по формулам:для горючих газов (ГГ)
; (3.31) , (3.32)где
для паров ЛВЖ
; (3.33) , (3.34)где
K = T/3600;
T - продолжительность поступления паров в открытое пространство, с;
За начало отсчета горизонтального размера зоны принимают внешние габаритные размеры пролива.
при взрыве облака топливно-воздушной смеси
Методика количественной оценки параметров воздушных ударных волн при взрыве топливно-воздушных смесей (ТВС) распространяется на случаи выброса топлива в атмосферу из наружных установок.
Основными структурными элементами алгоритма расчетов являются:
- определение ожидаемого режима взрывного превращения облака ТВС;
- расчет максимального избыточного давления и импульса фазы сжатия воздушных ударных волн для различных режимов;
- определение дополнительных характеристик взрывной нагрузки;
- оценка поражающего воздействия.
Исходными данными для расчета параметров ударных волн при взрыве облака ТВС являются:
- вид топлива, содержащегося в облаке;
- агрегатное состояние смеси (газовая или гетерогенная);
- концентрация горючего в смеси
- стехиометрическая концентрация горючего с воздухом
- масса топлива, содержащегося в облаке
- удельная теплота сгорания топлива
- скорость звука в воздухе
);- информация об окружающем пространстве;
- эффективный энергозапас горючей смеси (E), вычисляемый по соотношению
(3.35)При расчете параметров взрыва облака, расположенного на поверхности земли, величина эффективного энергозапаса удваивается.
Ожидаемый режим взрывного превращения зависит от типа топлива и геометрических характеристик окружающего пространства.
Вещества, способные к образованию горючих смесей с воздухом, по своей чувствительности к возбуждению взрывных процессов разделены на четыре класса.
Класс 1. Особо чувствительные вещества (размер детонационной ячейки менее 2 см). К этому классу относятся такие вещества как
, Класс 2. Чувствительные вещества (размер детонационной ячейки лежит в пределах от 2 до 10 см). Типичными представителями этого класса являются
Класс 3. Среднечувствительные вещества (размер детонационной ячейки лежит в пределах от 10 до 40 см). Типичными представителями этого класса являются гексан
, октан, изооктан, пары и распылы бензина, LNG.Класс 4. Слабочувствительные вещества (размер детонационной ячейки больше 40 см). Типичными представителями этого класса являются метан, декан, распылы дизтоплива, керосина, бензола.
Классификация наиболее распространенных в промышленном производстве горючих веществ приведена в табл. 3.1. В случае, если вещество не внесено в классификацию, его следует классифицировать по аналогии с имеющимися в списке веществами, а при отсутствии информации о свойствах данного вещества его следует отнести к классу 1, т.е. рассматривать наиболее опасный случай.
Таблица 3.1
(ударно-волнового) поражения
При оценке масштабов фугасного (ударно-волнового) поражения необходимо учитывать различие химических соединений по теплоте горения, используемой для расчета полного запаса энерговыделения. В расчетах энерговыделения не учитывается теплота конденсации водяных паров. Для типичных углеводородов берется в расчет значение удельной теплоты сгорания
. Для иных горючих в расчетах используется удельное энерговыделение . Здесь Таблица 3.2
В связи с тем, что характер окружающего пространства в значительной степени определяет скорость взрывного превращения облака топливовоздушной смеси и, следовательно, параметры ударной волны, геометрические характеристики окружающего пространства также разделены на четыре класса в соответствии со степенью их опасности.
Класс I. Наличие длинных труб, полостей, каверн, заполненных горючей смесью, при сгорании которой возможно ожидать формирование турбулентных струй продуктов сгорания размером не менее трех размеров детонационной ячейки данной смеси. Если размер детонационной ячейки для данной смеси не известен, то минимальный характерный размер струй принимается равным 5 см для веществ класса 1; 20 см - для веществ класса 2; 50 см - для веществ класса 3 и 150 см - для веществ класса 4.
Класс II. Сильнозагроможденное пространство: наличие полузамкнутых объемов, высокая плотность размещения технологического оборудования, лес, большое количество повторяющихся препятствий.
Класс III. Среднезагроможденное пространство: отдельно стоящие технологические установки, резервуарный парк.
Класс IV. Слабозагроможденное и свободное пространство.
Для оценки действия взрыва возможные взрывные режимы превращения топливовоздушной смеси разбиты на шесть классов по диапазонам скоростей их распространения.
Класс 1. Детонация или горение со скоростью фронта пламени
.Класс 2. Дефлаграция, скорость фронта пламени
.Класс 3. Дефлаграция, скорость фронта пламени
.Класс 4. Дефлаграция, скорость фронта пламени
.Класс 5. Дефлаграция, скорость фронта пламени определяется соотношением
где
M - масса топлива, содержащегося в облаке горючей смеси, кг.
Класс 6. Дефлаграция, скорость фронта пламени определяется соотношением
, (3.37)где
M - масса топлива, содержащегося в облаке горючей смеси, кг.
Ожидаемый режим взрывного превращения определяется с помощью экспертной табл. 3.3, в зависимости от класса топлива и класса окружающего пространства.
Таблица 3.3
При определении максимальной скорости фронта пламени для взрывных процессов 2 - 4 классов дополнительно рассчитывается видимая скорость фронта пламени по соотношению (3.36). В том случае, если полученная величина больше максимальной скорости, соответствующей данному классу, она принимается за верхнюю границу диапазона ожидаемых скоростей взрывного процесса в топливовоздушной смеси.
Для дальнейших расчетов необходимо оценить агрегатное состояние топлива смеси. Предполагается, что смесь гетерогенная, если более 50% топлива содержится в облаке в виде капель. Провести такие оценки можно исходя из величины давления насыщенных паров топлива при данной температуре и времени формирования облака. Для летучих веществ, таких как пропан, при температуре +20 °C смесь можно считать газовой, а для веществ с низким давлением насыщенного пара (распыл дизтоплива при +20 °C) расчеты проводятся в предположении гетерогенной топливовоздушной смеси.
и импульса фазы сжатия воздушных ударных волн
После того, как определен вероятный режим взрывного превращения, рассчитываются параметры воздушных ударных волн (избыточное давление
Детонация газовых и гетерогенных ТВС
Для вычисления параметров воздушной ударной волны на заданном расстоянии от центра облака, при детонации облака ТВС, предварительно рассчитывается соответствующий безразмерный радиус по соотношению
где R - расстояние от центра облака;
E - эффективный энергозапас смеси.
Если соотношения записаны в функции аргумента
, величина E задается в МДж, то . После этого рассчитываются величины безразмерного давления В случае детонации газовой смеси расчет проводится по формулам:
Зависимости (3.39, 3.40) справедливы для значений
, в случае, если , то .В случае детонации облака гетерогенной топливовоздушной смеси
или
; (3.41') . (3.42')Зависимости (3.41, 3.42) справедливы для значений
, в случае, если , то .После определения безразмерных величин давления и импульса фазы сжатия вычисляются соответствующие им размерные величины:
; (3.43) . (3.44)Дефлаграция газовых и гетерогенных ТВС
В случае дефлаграционного взрыва облака ТВС к параметрам, влияющим на величины избыточного давления и импульса положительной фазы, добавляются величины скорости видимого фронта пламени (u) к степени расширения продуктов сгорания
.Как и для случая детонации, предварительно рассчитывается соответствующее безразмерное расстояние
Затем рассчитываются величины безразмерного давления
где
.Выражения (3.45, 3.46) справедливы для значений
, в случае, если , в выражения (3.45, 3.46) вместо (3.47)с перегретой жидкостью или сжиженным газом при воздействии
на него очага пожара
Избыточное давление
; (3.48) , (3.49)где
- приведенная масса, кг; (3.50)r - расстояние от центра резервуара, м;
; , (3.51)где k - доля энергии ударной волны (допускается принимать равной 0,5);
;m - масса ЛВЖ, ГЖ или СУГ, содержащаяся в резервуаре, кг;
T - температура жидкой фазы, К;
При наличии в резервуаре предохранительного устройства (клапана или мембраны) величина T определяется по формуле
, (3.52)где
A, B,
В настоящем разделе приводятся методы расчета интенсивности теплового излучения от пожара пролива и огненного шара, а также радиуса воздействия продуктов сгорания паровоздушного облака в случае пожара-вспышки.
Интенсивность теплового излучения q,
, для пожара пролива ЛВЖ, ГЖ или СУГ вычисляется по формулегде
;Значение
.Таблица 3.4
Среднеповерхностная плотность теплового излучения пламени
в зависимости от диаметра очага и удельная массовая
скорость выгорания для некоторых жидких
углеводородных топлив
Примечание. Для диаметров очага менее 10 или более 50 м следует принимать
Угловой коэффициент облученности
, (3.54)где
(3.55) (3.56) ; (3.57) ; (3.58)S = 2r/d; (3.59)
h = 2H/d, (3.60)
где r - расстояние от геометрического центра пролива до облучаемого объекта, м;
d - эффективный диаметр пролива, м;
H - высота пламени, м.
Эффективный диаметр пролива d, м, рассчитывается по формуле
, (3.61)где F - площадь пролива, м2.
Высота пламени H, м, вычисляется по формуле
, (3.62)где m - удельная массовая скорость выгорания топлива,
;g - ускорение свободного падения, принимаемое равным
.Коэффициент пропускания атмосферы
. (3.63)Интенсивность теплового излучения q
для огненного шара вычисляется по формуле (3.53).Величина
.Значение
, (3.64)где H - высота центра огненного шара, м;
r - расстояние от облучаемого объекта до точки на поверхности земли непосредственно под центром огненного шара, м.
Эффективный диаметр огненного шара
, (3.65)где m - масса продукта, поступившего в окружающее пространство, кг.
Величину H допускается принимать равной
Время существования огненного шара
. (3.66)Коэффициент пропускания атмосферы
. (3.67)паровоздушного облака в случае пожара-вспышки
В случае образования паровоздушной смеси в незагроможденном технологическим оборудованием пространстве и его зажигании относительно слабым источником (например, искрой) сгорание этой смеси происходит, как правило, с небольшими видимыми скоростями пламени. При этом амплитуды волны давления малы и могут не приниматься во внимание при оценке поражающего воздействия. В этом случае реализуется так называемый пожар-вспышка, при котором зона поражения высокотемпературными продуктами сгорания паровоздушной смеси практически совпадает с максимальным размером облака продуктов сгорания (т.е. поражаются в основном объекты, попадающие в это облако). Радиус воздействия высокотемпературных продуктов сгорания паровоздушного облака при пожаре-вспышке
, (3.68)где
;Интенсивность испарения W
для ненагретых жидкостей с удовлетворительной точностью может быть описана полуэмпирическим выражением , (3.69)где
M - молярная масса жидкости,
;Таблица 3.5
Значения коэффициента
При выбросе СУГ из поврежденного оборудования, в котором жидкость находится под давлением, часть продукта за счет внутренней энергии мгновенно испаряется, образуя с капельками жидкости облако аэрозоля. Массовую долю мгновенно испарившейся жидкости
, (3.70)где
; .Принимается, что при
вся масса жидкости, находящаяся в оборудовании, за счет взрывного характера испарения переходит в парокапельное облако.При
оставшаяся часть жидкости испаряется с поверхности пролива за счет потока тепла от подстилающей поверхности и воздуха.Интенсивность испарения жидкости со свободной поверхности W
описывается выражением , (3.71)где
; ;t - текущее время с момента начала испарения, с (но не менее 10 с);
u - скорость воздушного потока над поверхностью испарения,
d - характерный диаметр пролива, м;
Расчет геометрических параметров факелов при истечении паровой фазы СУГ проводится на основе схемы, показанной на рис. 3.4.
![]() Рис. 3.4. Схема для расчета геометрических параметров
факелов при истечении паровой фазы СУГ
Приведенный метод применим, когда скорость ветра больше
.При меньших значениях скорости ветра этот угол следует считать равным нулю.
Геометрические параметры газовых факелов в случае истечения паровой фазы СУГ вычисляются по формулам (для скорости ветра не менее
): ; (3.72) ; (3.73) ; (3.74) ; (3.75) ; (3.76) ; (3.77) ; (3.78) , (3.79)где
; ;Длина и диаметр факела, образующегося при истечении жидкой фазы СУГ, описываются формулами:
; (3.80) ; (3.81) ; (3.82) , (3.83)где
F - площадь отверстие истечения жидкой фазы, м2;
g - ускорение свободного падения,
B - диаметр факела, м.
Длина
; (3.84) , (3.85)где
G - расход газа,
k = 12,3 - коэффициент пропорциональности.
Размеры факела принимаются независимыми от направления истечения газа.
При истечении газа вдоль ограничивающей поверхности (земля, оборудование и т.д.) расчетную длину факела следует увеличить на 25%.
Оценку интенсивности теплового излучения пламени горящего резервуара (горение жидкости по всей площади поперечного сечения резервуара) можно получить по методу, изложенному в п. 3.6, принимая в качестве диаметра пролива диаметр резервуара.
ОПАСНЫМИ ФАКТОРАМИ АВАРИЙ С ПОЖАРАМИ И ВЗРЫВАМИ
Для аварий с пожарами и взрывами на промышленных объектах наиболее характерными поражающими факторами являются ударная волна, тепловое излучение и осколочное воздействие при взрывах технологических аппаратов.
Ниже рассмотрены различные подходы к определению критериев поражения людей, зданий и оборудования перечисленными поражающими факторами.
Известные критерии поражения объектов ударной волной (УВ) можно условно разделить на детерминированные и вероятностные. Детерминированные критерии приближенно показывают значения параметров УВ, при которых наблюдается тот или иной уровень поражения (разрушения). Вероятностные критерии показывают, какова условная вероятность того или иного уровня поражения (разрушения) при заданном значении поражающего фактора (например, избыточного давления и импульса положительной фазы УВ). Дальше рассматриваются как те, так и другие критерии.
поражающего действия УВ
Для описания поражения различных объектов УВ принято использовать так называемые P/i диаграммы. Для кратковременных импульсных воздействий зачастую приемлемым является использование только критических значений импульса i УВ, для относительно длительных воздействий, типичных для дефлаграционных взрывов паровоздушных облаков - критических значений избыточного давления
Воздействие УВ на конструкции во многом определяется величиной
воздействие определяется величиной - импульсом УВ. Собственные частоты некоторых сооружений приведены в табл. 4.1.Таблица 4.1
Значения критического давления
Таблица 4.2
Значения критического давления
Таблица 4.4
Значения критического давления
Таблица 4.5
поражающего действия УВ
В качестве вероятностного критерия поражения людей и/или зданий и сооружений используется понятие пробит-функции (probit-function). В общем случае пробит-функция
, (4.1)где a, b - константы, зависящие от степени поражения и вида объекта;
S - интенсивность воздействующего фактора.
Соотношения между величиной
Таблица 4.6
Значения пробит-функции
Для воздействия УВ на человека выражения для пробит-функции имеют вид
; (4.2) ; (4.3) , (4.5)где m - масса тела человека, кг;
i - импульс УВ, Па·с;
Пробит-функции для разрушения зданий имеют вид:
для тяжелых разрушений
; (4.7)для полного разрушения
; (4.8)Для практических целей целесообразно использовать пробит-функцию для тяжелых разрушений зданий (выражение 4.6), поскольку наиболее вероятный сценарий гибели людей от взрыва - это поражение их частями разрушающегося здания, расположенного вблизи взрывоопасного объекта.
Для теплового излучения могут быть использованы как детерминированные, так и вероятностные критерии поражения.
тепловым излучением
При анализе воздействия теплового излучения следует различать случаи импульсного и длительного воздействия. В первом случае следует говорить о дозе излучения D (случай огненного шара), во втором - о критической интенсивности теплового излучения
Величины
Таблица 4.7
Детерминированные критерии поражения технологического оборудования и строительных конструкций определяются, как правило, из условия потери ими несущей способности.
Основным критерием потери несущей способности является достижение на необогреваемой поверхности (например, внутренней поверхности резервуара, вблизи которого находится очаг пожара) критической температуры. Некоторые данные о значениях критических температур строительных материалов приведены ниже в табл. 4.9.
Оценка поражающего действия теплового излучения может быть проведена путем решения уравнения теплопроводности для облучаемого объекта (например, резервуара) с определением момента достижения температурой конструкционных материалов критических значений, приведенных в табл. 4.9.
Для приближенной оценки поражающего действия теплового излучения можно использовать метод, основанный на оценке критических величин плотностей тепловых потоков, вызывающих нагрев конструкционных материалов до критических температур. При этом величины критических тепловых потоков могут определяться как экспериментальным, так и расчетным путем.
Имеются экспериментальные данные по минимальной плотности теплового потока
В литературе описан метод расчетно-аналитического определения противопожарных разрывов между резервуарами с СУГ и ЛВЖ. Предполагается, что горение резервуара с ЛВЖ происходит во всем обваловании резервуара. В качестве критической температуры несмоченной стенки резервуара с СУГ принимается величина 427 °C. Расчеты показали, что плотность падающего на стенки резервуара теплового потока не должна превышать
в отсутствии ветра и при скорости ветра .Исходя из вышеизложенного, могут быть сделаны следующие практические рекомендации для учета поражающего действия теплового излучения по отношению к технологическому оборудованию и резервуарам. В случае полного охвата пламенем стенок резервуара представляется достаточным решение тепловой задачи в обыкновенных дифференциальных уравнениях с использованием критических температур из табл. 4.10. При несимметричном нагреве возможно использование понятия критической плотности теплового потока
.поражения тепловым излучением
В качестве вероятностного критерия оценки поражения тепловым излучением целесообразно использовать описанное выше понятие пробит-функции.
Для смертельного поражения человека величина пробит-функции описывается следующим выражением
; (4.11) . (4.12)Величина эффективного времени экспозиции t может быть вычислена по формулам:
для огненного шара
; (4.13)для пожара пролива
, (4.14)где m - масса горючего вещества, участвующего в образовании огненного шара, кг;
x - расстояние от места расположения человека до безопасной зоны (зона, где интенсивность теплового излучения меньше
);u - средняя скорость движения человека к безопасной зоне (может быть принята
).Величина дозы излучения D в случае пожара-вспышки для смертельного поражения человека может быть вычислена следующим образом. Вначале определяется эффективный радиус продуктов сгорания R по формуле
, (4.15)где
Величина D
определяется по следующим формулам в зависимости от расстояния r от геометрического центра паровоздушного облака:
Величины D для промежуточных значений r следует определять с помощью линейной интерполяции.
Оценку индивидуального и социального риска для зданий промышленных предприятий следует проводить в соответствии с ГОСТ Р 12.3.047-98 (приложение Ш).
Оценку потенциального, индивидуального, социального и коллективного рисков для территории предприятия следует проводить по изложенным ниже методикам.
Величина потенциального риска P(a)
в определенной точке местности a, где расположено предприятие, определяется с помощью соотношения , (5.1)где I - число сценариев развития аварии (ветвей логического дерева событий возникновения и развития аварии);
Условные вероятности поражения человека
Величина P(a) определяется посредством наложения зон поражения опасными факторами с учетом частоты реализации каждого сценария развития аварии на карту местности с привязкой их к соответствующему событию, инициирующему аварию (элементу оборудования, технологической установке) и ориентированию зоны поражения в соответствии с метеорологическими условиями (для струйного горения, пожара-вспышки, образования и взрывного превращения газопаровоздушного облака). При расчете риска рассматриваются различные метеорологические условия с типичными направлениями ветров и ожидаемой частотой их возникновения.
Процедура расчета риска предусматривает рассмотрение различных аварийных ситуаций и определение зон поражения опасными факторами пожара и взрыва, и частот их реализации. Для удобства расчетов территорию местности делят на зоны, внутри которых величины P(a) полагаются одинаковыми.
В необходимых случаях оценка условной вероятности поражения человека проводится с учетом совместного воздействия более чем одного опасного фактора (для ветвей со стадиями с условием перехода "И"). Так, например, для расчета условной вероятности поражения человека при реализации сценария развития аварии, связанного со взрывом резервуара с ЛВЖ под давлением, находящегося в очаге пожара, необходимо учитывать, кроме теплового излучения огненного шара, воздействие ударной волны и осколков.
Условная вероятность поражения человека
, (5.2)где h - число рассматриваемых опасных факторов аварии;
- условная вероятность поражения k-м опасным фактором.Результаты расчетов потенциального риска отображаются на карте (ситуационном плане) предприятия и прилегающих районов в виде замкнутых линий равных значений (изолинии функции P(a)).
Изолинии функции P(a) называются контурами риска. Их физический смысл состоит в том, что они разделяют территорию предприятия (так же, как и местность вокруг предприятия) на области, в которых ожидаемая частота возникновения опасных факторов аварии, приводящих к гибели людей, заключена в определенных, указанных на рисунке, пределах.
Контуры риска не зависят от количества работающих на предприятии или их должностных обязанностей, а определяются исключительно используемой технологией и надежностью применяемого оборудования. Потенциальный риск используется как мера (критерий допустимости/недопустимости) уровня пожарной безопасности объекта.
Для любого работника предприятия существует возможность гибели при возникновении аварии. Потеря жизни в течение определенного периода времени (года) является случайным событием, зависящим от рода его профессиональной деятельности, в том числе от продолжительности нахождения работника в областях, отвечающих различным контурам риска при его перемещениях по промышленной площадке предприятия в течение рабочей смены.
Для целей управления безопасностью персонала используется количественная мера возникновения этого случайного события - частота поражения опасными факторами пожара (взрыва) определенного человека, называемую индивидуальным риском.
Таким образом, индивидуальный риск определяется как ожидаемая частота поражения определенного работника предприятия опасными факторами аварий в течение года.
Области, на которые разбита территория предприятия, обозначается
J, j = 1,..., J.
Для удобства описания расчетов работники предприятия нумеруются
m = 1,..., M.
Текущий номер работника m однозначно определяет наименование должности работника, его категорию и другие особенности его профессиональной деятельности, необходимой для оценки безопасности.
Величина индивидуального риска
для работника предприятия m определяется с помощью соотношения , (5.3)где P(j) - величина потенциального риска в j-й области территории предприятия,
Доля времени, в течение которого работник находится в определенной области территории предприятия, рассчитывается на основе решений по организации эксплуатации и технического обслуживания оборудования.
Для анализа воздействия промышленных аварий на людей, а также для установления приемлемости того или иного уровня пожарной или промышленной безопасности используется понятие социального риска.
Социальный риск задается с помощью функции, значениями которой являются величины, определяющие, что в аварии с пожаром погибло не менее определенного количества человек.
Социальный риск S
определяется по формуле , (5.4)где L - число сценариев развития аварии, для которых выполняется условие
; .Ожидаемое число погибших в результате реализации i-го сценария развития аварии можно оценить по следующей формуле
, (5.5)где J - количество областей, на которые разбита территория предприятия и территория, прилегающая к предприятию (j - номер области);
Социальный риск S является интегральной величиной. В то же время социальный риск можно рассматривать и как векторную величину, компоненты которой имеют размерность
. (5.6)Для персонала предприятия в целом имеется ненулевая вероятность гибели части работников при возникновении аварии.
Количество погибших в течение определенного периода времени (года) является случайной величиной, зависящей от опасности производства, количества работающих и ряда других факторов.
Для целей управления безопасностью персонала используется математическое ожидание этой случайной величины. Указанную характеристику называют коллективным риском персонала от аварий.
Величина коллективного риска персонала C
определяется с помощью соотношения . (5.7)Связь индивидуального риска работника от аварий
. (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 Inlel 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, N 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. - The Hague, 1997. - 836 p.
10. Sallet D.W. Critical two-phase mass flow rates of liquefied gases//Journal of Loss Prevention in the Process Industries. - 1990. - Vol. 3, N 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, N 1. - p. 179 - 192
13. Шебеко Ю.Н., Шевчук А.П., Смолин И.М. Расчет влияния обвалования на растекание горючей жидкости при разрушении резервуара//Химическая промышленность. - 1994. - N 4. - С. 230 - 233
14. Андерсон Д., Таинненхил Дж., Плетчер Р. Вычислительная гидромеханика и теплообмен. - М: Мир, 1990. - Т. 1. - 179 с.
15. НПБ 105-03. Определение категорий помещений, зданий и наружных установок по взрывопожарной и пожарной опасности
16. Пожаровзрывоопасность объектов хранения сжиженного природного газа. Процессы испарения и формирования пожаровзрывоопасных облаков при проливе жидкого метана. Методики оценки параметров/И.А. Болодьян, В.П. Молчанов, О.И. Дешевых и др.//Пожарная безопасность. - 2000. - N 4. - С. 108 - 121
17. Пожарная опасность шаровых резервуаров для хранения под давлением сжиженных углеводородных газов/А.П. Шевчук, Ю.Н. Шебеко, И.М. Смолин и др.//Химическая промышленность. - 1992. - N 6. - С. 328 - 332
18. Борисов А.А., Гельфанд Б.Е., Цыганов С.А. О моделировании ударных волн давления, образующихся при детонации и горении газовых смесей//Физика горения и взрыва. - 1985. - Т. 21, N 2. - С. 90 - 97
19. Методика оценки последствий аварийных взрывов топливно-воздушных смесей. - М.: НТЦ "Промышленная безопасность". - 45 с.
20. Шебеко Ю.Н., Шевчук А.П., Смолин И.М. Расчет параметров ударных волн, образующихся при взрыве резервуара со сжиженным углеводородным газом в очаге пожара//Химическая промышленность. - 1993. - N 9. - С. 451 - 453
21. Mudan K.S. Thermal radiation hazards from hydrocarbon pool fires//Progress in Energy and Combustion Science. - 1984. - Vol. 10, N 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, N 3. - P. 329 - 365
23. Орлов В.Я., Никитин А.Г. Определение скорости испарения легковоспламеняющихся жидкостей при категорировании производств: Вопросы охраны труда при проектировании и строительстве зданий. - Казань: КХТИ, 1978. - С. 55 - 57
24. Пчелинцев В.А., Никитин А.Г., Хузиахметов Р.А. Оценка взрывопожароопасности производств, связанных с применением легковоспламеняющихся жидкостей в нагретом состоянии//Журн. ВХО им. Д.И. Менделеева. - 1985. - Т. 30, N 1. - с. 68 - 74
25. Маршалл В. Основные опасности химических производств. - М.: Мир, 1989. - 671 с.
26. Bais A.F., Zerefos C.S., Ziomas I.C. 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.C., Zerefos G.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. - N 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, N 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, N 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, N 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. Broeckmann B., Schecker H.G. Heat transfer mechanisms and boilover in burning oil-water systems//Journal of Loss Prevention in the Process Industries. - 1995. - Vol. 8, N 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, N 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, N 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. - N 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. - 326 p.
46. Pietersen C.M. Consequences of accidental releases of hazardous material//Journal of Loss Prevention in the Process industries. - 1990. - Vol. 3, N 1. - 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, N 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, N 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, N 1 - 3. - P. 357 - 374.
Вернуться в "Каталог нормативных документов"
Источник информации: https://internet-law.ru/documents/prod/metodika/6/rukov.html
На правах рекламы:
|
|||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||