Развитие нефтяной промышленности диктует жёсткие условия обеспечения комплексной безопасности объектов с целью защиты как коммерческих интересов, так и населения и территорий. Особую актуальность обеспечение безопасности приобретает сегодня из-за критического сближения опасных производственных объектов с территорией активной жизнедеятельности населения.
В этой связи анализ рисков возникновения ЧС необходим для обоснованного выбора методов предупреждения неблагоприятных последствий воздействия опасностей на объекты. На сегодняшний день не развита единая методологическая база по оценке риска возникновения ЧС, слабо развиты методические основы учета вариабельности последствий ЧС в резервуарных парках, в том числе эффекта домино и эффекта BLEVE. В этой связи необходимо развитие методических основ прогноза риска ЧС в резервуарных парках.
Первые разработки в оценке риска (на основе дерева отказов), появились в 60 годы прошлого столетия. В последующие 10 лет начали применять методы анализа причин и последствий. Широкого развития и использования достигли множество методов системного анализа для оценки безопасности объекта экономики, такие как дерево отказов, дерево событий, сети GERT, моделирование на основе орграфов и другие [3, 6, 7, 9].
При выборе метода анализа риска точность и адекватность оценки зависят от этапа функционирования рассматриваемого объекта, цели анализа, типа производственного объекта, характера опасности и т.п. Владение необходимой информацией является важным условием проведения оценки риска [3]. Оценка риска возникновения чрезвычайной ситуации в резервуарных парках включает анализ и обработку статистических данных об аварийности на объектах хранения нефтепродуктов на основе годовых отчетов о деятельности Федеральной службы по экологическому, технологическому и атомному надзору в период с 2005 по 2015 г. и прочих информационных ресурсов. Анализ литературных источников показал, что в среднем вероятность возникновения аварии на объектах нефтепродуктообеспечения составляет 3,6×10–3 год–1 [2].
На основе всей совокупности методов можно представить резервуарный парк как систему, состоящую из различных взаимосвязанных элементов, обладающих различными свойствами (рис. 1).
Анализ резервуарного парка как системы предполагает анализ взаимосвязей его отдельных элементов. Поэтому целесообразно использовать метод анализа на основе орграфов. Для резервуарного парка с помощью взвешенных орграфов можно определить влияние различных компонентов системы на резервуары, приводящие к ЧС (разрушающее действие). На рис. 2 представлен взвешенный орграф для резервуарного парка. Весовые значения влияния внешних и внутренних факторов, приводящих к разрушению резервуара с нефтепродуктом, получены методом экспертных оценок (20 экспертов, среднее значение). Анализ орграфа позволяет выявить как метод взвешенных орграфов также используется при анализе загрязнения окружающей среды. Так, в работе [8] представлен взвешенный орграф, описывающий экологическую ситуацию в городе, который позволяет определить степень влияния различных факторов на уровень загрязнения города.
При решении как экологических, так и техногенных задач установить основные причины и определить вероятность возникновения аварии можно, используя метод построения дерева отказов, который нашел наибольшую применимость при оценке риска.
Метод дерева отказов широко используется при анализе риска возникновения ЧС на опасных производственных объектах, при выявлении причин выхода из строя технического оборудования. Так, в работе [1] данный метод используется для выявления причин отказа технологического оборудования, применяемого на нефтеперерабатывающих и нефтехимических производствах.
При анализе риска возникновения ЧС помимо причин особое значение имеет оценка последствий. Анализ литературных источников показал, что наиболее эффективным методом является дерево событий, представляющее собой схему причинно-следственных закономерностей возникновения и развития аварий. В работе [4] приведено дерево событий для определения последовательности и сочетания различных событий, возникающих в результате разлива бензина на проезжую часть.
С помощью дерева событий можно определить вариацию развития пролива при прохождении ливневых осадков с низкой вероятностью воспламенения, (с преимущественной угрозой загрязнения литосферы и гидросферы) и без осадков (с высокой вероятностью образования пожара пролива). Важным этапом прогноза является учет возможности возникновения эффекта BLEVE – как элемента наиболее опасного развития сценария ЧС в резервуарном парке.
Рис. 1. Модель системы «Резервуарный парк»
Рис. 2. Взвешенный орграф для резервуарного парка
В то же время методическая основа оценки возможности эффекта BLEVE находится на начальном этапе. Как показал анализ исследований в данной области, при оценке возможности возникновения эффекта BLEVE применяют несколько методов.
Широкое использование нашла методика оценки BLEVE согласно ГОСТу Р 12.3.047-2012, в которой рассчитывают показатель δ, характеризующий возможность возникновения волн сжатия и зависящий от свойств хранящегося вещества и технологических особенностей предохранительных клапанов. Если δ < 0,35, BLEVE не происходит. При δ > 0,35 вероятность возникновения данного явления велика.
В то же время при расчете показателя δ не учитывается вероятность попадания резервуара в зону пожара пролива (принимается условно, что резервуар находится в зоне воздействия теплового излучения). В этом случае можно использовать вероятностный метод, который основан на представлении процессов и явлений как случайных событий [5]. Тогда для вероятностной оценки устойчивости резервуаров возможно использование подхода по оценке геометрической вероятности в условиях равновозможного проявления эффекта BLEVE на территории объекта. Здесь важное значение приобретает установление границ объектов, рассматриваемых при оценке.
При установлении границ площади территории, которая будет находиться в зоне действия теплового излучения пожара, необходимо учесть неравномерность пространственного распределения температуры и интенсивности Также необходимо определить границы основания резервуаров. При расположении резервуара или нескольких резервуаров в обваловании площадь границ пожара пролива может быть равна площади обвалования. Тогда вероятность возникновения эффекта BLEVE будет определяться по формуле
(1)
где n – количество резервуаров на рассматриваемой территории (территория, объединенная особенностями рельефа местности или единым обвалованием, шт. резервуаров; Sосн.рез – площадь основания резервуара, м2; Sрассм.тер – площадь рассматриваемой территории, м2; Sпож.прол – площадь пожара пролива, м2.
Рис. 3. Расчётная схема оценки вероятности образования эффекта BLEVE в резервуарном парке
Расчётная схема для определения площадей представлена на рис. 3.
Полученное значение вероятности при анализе риска может быть использовано как значение одного из звеньев «дерева отказов» либо как вероятность инициирующего события в «дереве событий», либо как вероятностная составляющая значения риска. В то же время такой подход не учитывает конструктивные особенности резервуаров и интенсивность теплового излучения пожара пролива.
Конструктивные особенности частично можно учесть исходя из оценки пределов огнестойкости материала резервуара и давления, возникающего внутри него. На основании прогноза обстановки рассчитывается предел огнестойкости для резервуаров, находящихся в зоне действия пожара, и определяется возможность разгерметизации резервуара в зависимости от величины интенсивности теплового излучения, температуры и времени их воздействия, с учетом материала, из которого сделан резервуар. Температура на поверхности резервуара определяется на основании закона излучения, по формуле
(2)
где Е – энергия передаваемая поверхности при пожаре пролива, Вт/м2; σ – постоянная Стефана – Больцмана, равна 5.67·10−8 Вт/м2К4; Ат – степень черноты, при горении нефтепродуктов допускается принимать 0,9.
Энергия, передаваемая поверхности при пожаре пролива, определяется как интенсивность теплового излучения. Если условно предположить, что давление в резервуаре подчиняется уравнению Менделеева – Клапейрона, то его величину можно определить по формуле
(3)
где P1 – атмосферное давление, кПа; Т1 – температура окружающей среды, °С; Т2 – температура после нагрева, составляет, °С:
С учетом условия прочности по гипотезе О. Мора максимальное давление в резервуаре равно
(4)
где d – диаметр резервуара, м; t – толщина стенки резервуара, равна м; [sp] – допускаемом напряжении на растяжение, МПа.
Если рассчитанное давление превышает допустимое значение, происходит разрушение резервуара и существует вероятность возникновения эффекта BLEVE.
Поскольку поведение паров нефтепродуктов сложно сравнить с поведением идеального газа, описываемым уравнением Менделеева – Клапейрона, то возможность возникновения эффекта BLEVE можно оценить, сопоставив время достижения нефтепродуктом внутри резервуара температуры кипения и время воздействия теплового потока на резервуар с учетом интенсивности теплового излучения.
Рис. 4. Алгоритм оценки возможности возникновения эффекта BLEVE
Для оценки продолжительности нагрева нефтепродукта в резервуаре необходимо рассчитать количество энергии, затрачиваемой на приращение температуры нефтепродукта до значения кипения. Оценка времени нагрева содержимого резервуара до температуры кипения проводится на основании расчётного значения интенсивности теплового излучения с учётом площади поглощающей поверхности резервуара:
(5)
где q – интенсивность теплового излучения, воздействующая на резервуар, Вт/м2; m1 – количество нефтепродукта в резервуаре, кг; Δt – изменение температуры, °С.
Если время достижения нефтепродуктом температуры кипения меньше времени воздействия теплового излучения пожара, то существует вероятность возникновения эффекта BLEVE.
Таким образом, обоснованный выбор подходов к оценке вероятности возникновения эффекта BLEVE можно представить в виде алгоритма (рис. 4).
Таким образом, существующие методы системного анализа и алгоритмы принятия решений позволяют учесть многие факторы влияния, визуально отобразить взаимосвязь причин и вариацию последствий различных ЧС, что повысит своевременность реагирования на изменение состояния безопасности. При анализе возможности возникновения эффекта BLEVE необходимо применение обоснованного выбора метода и/или нескольких методов оценки комплексного учета всех факторов влияния.