34Вестник Командно-инженерного института МЧС Республики Беларусь, №1 (9), 2009 УДК. 614.8
ОЦЕНКА ИНДИВИДУАЛЬНОГО РИСКА ОТ ОБЪЕКТОВ С АММИАКОМ И ХЛОРОМ
Гаврилюк Д.АД Ильюшонок A.B., к.ф.-м.н., доцент, Катков В Л.*, д.т.н., профессор, Лешенюк Н.С., д.ф.-м.н., профессор, Отчик B.C., к.ф.-м.н., Попов В.М.**, к.ф.-м.н.
* Объединенный институт проблем информатики HAH Беларуси, **Институт физики им. Б.И. Степанова HAH Беларуси
Рассмотрена методология оценка индивидуального риска от химико-технологических объектов. Описан программно-информационный комплекс РИСК-1, созданный в рамках ГППНИ «Снижение рисков чрезвычайных ситуаций». Приводятся примеры работы комплекса - карты уровней концентраций, токсодоз и индивидуальных рисков, полученные в тестовых расчетах для некоторых опасных химических объектов г. Минска.
Республика Беларусь обладает разветвленной сетью предприятий, которые производят, транспортируют, хранят или используют в технологических процессах опасные химические вещества (ОХВ). Наиболее распространенными группами опасных химических объектов являются аммиачные установки и склады с хлором. На территории республики примерно на 500 объектах используется аммиак, примерно на 100 - хлор.
Общепринятой количественной мерой уровня опасности, создаваемой технологическим объектом, служит оценка риска [1,2]. Такая оценка позволяет дифференцировать техногенные объекты по уровню потенциальной опасности, которую они представляют для человека и окружающей природной среды.
Наиболее общим показателем риска является произведение вероятности аварии на ущерб от нее [1-4]. Если в течение года на объекте может произойти к опасных событий, то выражение для оценки риска имеет вид:
где 0\ - вероятность наступления в течение года /-го опасного собьгшя;
- возможный ущерб от /-го события.
В случае оценки индивидуального риска от химико-технологического объекта под величиной ущерба М7 чаще всего понимается вероятность гибели человека в результате токсического воздействия попавших в атмосферу ОХВ.
Для нахождения вероятности токсического поражения человека используется пробит-функция вида:
где а, Ъ и п - константы, свои для каждого конкретного ОХВ; С{х, у) - концентрация вещества в точке (х, у) зоны заражения; I - время действия ОХВ [5].
В случае п ф 1 для величины Сп1 часто используется термин «токсическая нагрузка». Если концентрация вещества изменяется в течение времени воздействия, то токсическая нагрузка будет определяться как:
(Поступила в редакцию 26 мая 2008 г.)
к
(1)
Pr = a + b-]n[Cn(x,y)-t],
(2)
ВТ = )сп{х,у,1)(11, (3)
О
где /о - время экспозиции.
Итак, для оценки риска, создаваемого в точке (х, у) опасным химическим объектом, необходимо знать вероятность возникновения на объекте аварии с выбросом (или проливом) ОХВ и токсическую нагрузку.
Нормативная документация [6] рекомендует определять вероятность возникновения аварии на опасном объекте на основе статистических данных, для чего требуется репрезентативная статистическая выборка о количестве аварий на том или ином объекте (группе объектов) за предшествующий период времени. Считая аварии независимыми событиями, распределенными по биномиальному закону, можно показать, что относительная ошибка в при определении вероятности возникновения аварии по ее частоте определяется следующим выражением [7]:
8 = 1/фйж, (4)
где X - вероятность аварии на объекте в течение года;
т - количество объектов, за которыми осуществляется наблюдение; п - число лет, в течение которых это наблюдение происходит.
Если положить, что вероятность аварии на объекте, использующем аммиак (т « 500), составляет 10"4, то для оценки этой вероятности даже со 100%-й ошибкой требуются статистические данные за 20 лет. Следует обратить внимание, что рассчитанное на основании таких данных значение X не будет нести значимой информации о вероятности аварии на конкретном объекте, потому что наблюдаемые объекты не являются однотипными. При выборе информации об однотипных объектах, их количество будет существенно меньше 500, и, следовательно, необходимы статистические данные за временной интервал, превышающий 20 лет в несколько раз. Такие данные отсутствуют уже по той причине, что срок эксплуатации исследуемых объектов составляет примерно 30-40
Отсутствие репрезентативной статистики по наиболее распространенным в Республике Беларусь группам опасных химических объектов требует другого подхода для определения вероятности аварии на них. Для таких объектов вероятность рассчитывается, исходя из структуры объекта и логики развития физико-химических процессов в ходе аварии (теоретико-вероятностный метод оценки риска) [2]. Для аммиачных установок и складов с хлором наибольшую опасность представляют аварии, связанные с полной разгерметизацией резервуаров и трубопроводов, расположенных на открытых площадках [7]. Обобщенная информация по частоте отказов элементов на объектах с аммиаком и хлором приведена в таблице.
К числу наиболее распространенных методов расчета концентрации опасных веществ, поступающих в атмосферу в результате промышленных аварий, относятся гауссовские модели и модели рассеяния, базирующиеся на интегральных законах сохранения массы и энергии. В основе гауссовских моделей лежит предположение о нормальном (гауссовском) распределении вредных веществ в атмосфере. В них учитываются в основном два процесса, происходящие с выбросом в атмосфере: перемещение примеси под действием ветра и рассеяние ее за счет атмосферной турбулентности, которая определяется с помощью эмпирических коэффициентов.
Таблица - Обобщенная информация по частоте отказов элементов оборудования на объектах с аммиаком и хлором
Событие Частота, год'1 Выбрасываемая (выливаемая) масса и данные для ее расчета
Разгерметизация контейнера с хлором 1,0-Ю-4 Все содержимое контейнера. Емкость контейнера
Разгерметизация изотермы с аммиаком 1,0-10'6 Все содержимое изотермы. Емкость изотермы
Разгерметизация резервуара с аммиаком под давлением 1,0-10'4 Все содержимое резервуара. Емкость резервуара
Разгерметизация трубопровода с аммиаком (до 90 % случаев - выброс содержимого через отверстие в стенке трубопровода; 10 % случаев - полный разрыв трубопровода) 5-Ю'6 •м"1 Истечение до отключения трубопровода. Длина, диаметр трубопровода, рабочее давление, производительность насоса, время срабатывания системы автоматики отключения
На распространение выбросов с плотностью, большей плотности воздуха, влияет их растекание под действием сил плавучести. Для аммиака это не имеет существенного значения, игнорирование же процесса гравитационного растекания в гауссовой модели для хлора, плотность паров которого превышает плотность воздуха, приводит к завышению его концентрации вблизи источника.
На изложенных выше подходах базируется работа созданного в рамках ГППНИ «Снижение рисков чрезвычайных ситуаций» программно-информационного комплекса РИСК-1, предназначенного для оценки индивидуальных рисков от объектов с аммиаком и хлором. В случае мгновенного выброса примесей из точечного источника, расположенного над подстилающей поверхностью в точке с координатами (0, 0, /?), поле концентрации ОХВ описывается формулой [8]:
(
Мех р
(х-Ш)2 у2 2сг 2а2
ехр
2а2
+ ехр
' (2 + /02Л 2а:
(5)
где М- масса ОХВ в выбросе;
Ко - радиус облака выброса в начальный момент времени;
и - скорость ветра;
аЛ, Оу, <з2 - стандартные отклонения облака примеси по направлениям соответствующих осей координат.
Направление оси х выбирается вдоль ветра. Входящие в формулу (5) коэффициенты дисперсии ог, Оу, ог вычисляются с помощью эмпирических зависимостей Смита-Хоскера, отражающих возрастание дисперсии с увеличением расстояния X = и 4 от источника выброса до центра облака [1]. Конкретный вид этих выражений зависит от класса вертикальной устойчивости атмосферы.
Для расчета поля концентрации по формуле (5), а следовательно, и токсической нагрузки, требуется знать ряд метеорологических параметров: скорость и направление ветра, температуру атмосферного воздуха, класс его вертикальной стратификации. При оценке риска невозможно спрогнозировать значения этих параметров на момент аварии, поэтому
необходимо учитывать весь спектр возможных состояний атмосферы в пределах характерного периода изменений [4]. Республиканский гидрометеорологический центр на регулярной основе предоставляет необходимую информацию по месяцам и годам.
При построении полей рисков важным является вопрос о выборе величины временного интервала, за который определяются средние значения метеопараметров. При осреднении за год они оказываются сильно «сглаженными». Например, среднегодовые температуры в Минске получаются положительными, в результате чего месяцы с отрицательными температурами как бы «выпадают» из расчета. Уменьшение временного интервала, за который определяются средние значения параметров среды, приводит к повышению точности оценки, но увеличивает объем вычислений. Поскольку расчеты проводятся на компьютере, это увеличение не носит принципиального характера. С учетом сказанного в программно-информационном комплексе РИСК-1 в качестве временного интервала осреднения принят месяц.
Исходными данными для определения пространственного распределения полей рисков с учетом возможных состояний окружающей среды являются: средняя скорость и повторяемость ветра в различных направлениях (по 8-румбовой схеме) по месяцам; среднемесячная температура воздуха; среднее за месяц число ясных и пасмурных дней; среднее за месяц количество дней со снежным покровом.
Риск от /-й аварии в точке с координатами (х, у) вычисляется по формуле:
Я (*, УУ= %= £ N. X 1\г] ± Рвт ж (*, ^ 5, , Тт ), (6)
где Ыт - количество дней в т-м месяце;
Рут>] - вероятность направления ветра в т-м месяце в/-м румбе;
Рвт,и - относительная частота наблюдения той или иной вертикальной устойчивости атмосферы (я = 1 - инверсия, я = 2 - изотермия, £ = 3 - конвекция) в течение т-го месяца;
Ку - средние значения скорости ветра в /-м румбе;
Тт - средняя за т-й месяц температура воздуха.
Результирующее поле риска /-й аварии получается наложением 12x9x3 = 324 пространственных «сеток».
Программно-информационный комплекс работает в трех режимах. Оценка - оценка величины индивидуального риска аварии на опасном химическом объекте по метеонаблюдениям за определенный период. Анализ- анализ конкретной чрезвычайной ситуации, случившейся на опасном объекте, с расчетом полей концентраций и токсодоз. Регламент - проведение регламентных работ по поддержанию комплекса РИСК-1 в актуальном состоянии, в основном пополнение баз данных комплекса новой метеорологической информацией.
По окончании сеанса работы комплекса РИСК-1 молено автоматически получить типовой отчет по оценке риска, исходящего от опасного объекта, или отчет с анализом конкретной чрезвычайной ситуации. В отчет включаются: электронные карты с полем концентрации ОХВ (режим Анализ, рисунок 1), полем токсодоз (режим Анализ, рисунок 2), полем индивидуальных рисков (режим Оценка, рисунок 3); графики концентрации и токсодоз по оси «факела» выброса и др. Каждая карта сопровождается краткой легендой, в которой указаны адрес опасного объекта, дата и тип аварии, метеоусловия, характеристики выброса, прогнозное число людей, пострадавших при аварии, и т. д. На приведенных рисунках ради экономии места легенды не указываются, а в подписях под рисунками дается часть информации из легенды и текстовой части итогового отчета комплекса после сеанса его работы.
Разрушение емкости с 1 т жидкого хлора. Начальный размер облака выброса 2,5 м. Площадь пролива жидкого хлора - 14.0 м2. Ветер юго-западный, 2,1 м/с. Масштаб карты 1 : 15 ООО
Рисунок 1 - Поле приземной концентрации хлора через 1000 с после выброса при гипотетической аварии на Минской водоочистной водопроводной станции
Пороговая токсодоза - 0,6 мг-мин/л, летальная токсодоза - 6 мг-мин/л. Возможное число пострадавших от пороговой токсодозы составит 2 чел., от летальной токсодозы - 5 чел.; всего пострадает 7 чел. из общего числа 140 чел. работающей смены
Рисунок 2 - Поле токсодоз хлора через 1000 с после начала выброса на очистной водопроводной станции г. Минска
Рисунок 3 - Поле индивидуальных рисков для очистной водопроводной станции Минска в случае средних метеоусловий за январь - декабрь 2005 г. при 8-румбовой схеме задания ветра
В программно-информационном комплексе РИСК-1 можно использовать также аэрокосмические снимки с двоякой целью: для уточнения расположения опасных объектов на местности, а также для выдачи результатов расчета прямо на снимок. Пример такой выдачи для объекта БЕЛМЕДПРЕПАРАТЫ (6 т аммиака), расположенного на пересечении улиц Аранской и Маяковского в Минске, приведен на рисунке 4. Видно, что при аварии на
этом объекте в случае северо-восточного ветра может сформироваться поле токсодозы аммиака, накрывающее Червенский рынок, на котором в определенные часы может находиться много людей.
Созданный программный комплекс можно использовать в Министерстве по чрезвычайным ситуациям:
• при районировании территории республики по величине риска аварий на опасных объектах;
• для периодической (плановой) оценки рисков аварий на опасных химических объектах;
• при расследовании причин и последствий чрезвычайных ситуаций;
• при подготовке планов мероприятий по уменьшению рисков и снижению ущерба от аварий, для мониторинга и прогнозирования чрезвычайных ситуаций природного и техногенного характера.
Также этот комплекс можно использовать в Министерстве промышленности при подготовке паспортов безопасности предприятий, имеющих химически опасные объекты и производства; в Министерстве природных ресурсов и охраны окружающей среды: при осуществлении государственного контроля по охране окружающей среды, при организации работ по прогнозированию изменения экологической обстановки вследствие антропогенной нагрузки в районах размещения опасных производств и интенсивной хозяйственной деятельности; в высших учебных заведениях: при изучении дисциплин «Охрана окружающей среды», «Моделирование ветрового переноса загрязнений», «Опасные факторы чрезвычайных ситуаций» и др.
В дальнейшем развитии комплекса РИСК предполагается реализовать в нем методику «Токси», ред. 3.1 [9].
Минск, опасный объект БЕЛМЕДПРЕПАРАТЫ (рядом с Червенским рынком)
Рисунок 4 - Наложение рассчитанного поля токсодозы аммиака на космический снимок
ЛИТЕРАТУРА
1. Количественная оценка риска химических аварий / Колодкин В.М. [и др.]; под ред. В.М. Колодкина. - Ижевск: изд. дом «Удмурдский университет», 2001. - 228 с.
2. Акимов, В.А. Природные и техногенные чрезвычайные ситуации: опасности, угрозы, риски / В.А. Акимов, В.Д. Новиков, H.H. Радаев. - М.: «Деловой экспресс», 2001. - 344 с.
3. Владимиров, В.А. Оценка риска и управление техногенной безопасностью / В.А. Владимиров, В.И. Измалков, A.B. Измалков. - М: «Деловой экспресс», 2002. - 184 с.
4. Шахраманьян, М.А. Новые информационные технологии в задачах обеспечения национальной безопасности России (природно-техногенные факторы) / М.А. Шахраманьян. - М.: Информ.-изд. центр ФЦ ВНИИ ГОЧС, 2003. - 397 с.
5. Лошадкин, H.A. Пробит-метод в оценке физиологически активных веществ при низких уровнях воздействия / H.A. Лошадкин [и др.] // Рос. хим. журн. (ЖРХО). - 2002. - Т. 16, № 6. - С. 63-67.
6. Нормы пожарной безопасности Республики Беларусь. Категорирование помещений, зданий и наружных установок по взрывопожарной и пожарной опасности: НПБ 5-2005. - Введ. 01.07.06. - Минск: ОДО «Друк-С», 2006. - 42 с.
7. Ильюшонок, A.B. Оценка вероятности возникновения аварии на объектах химико-технологического комплекса/ A.B. Ильюшонок, С.А. Лодята // Чрезвычайные ситуации: предупреждение и ликвидация. - 2007. - № 1 (21). - С. 36-44.
8. Методика оценки последствий химических аварий (Методика «Токси». Редакция 2.2) // Методики оценки последствий аварий на опасных производственных объектах: сб. док. / ГУП «НТЦ Промышленная безопасность». - М., 2001. - С. 121-204.
9. Методика оценки последствий аварийных выбросов опасных веществ (Методика «Токси». Редакция 3.1): проект. - М.: ФГУПНТЦ «Промышленная безопасность», 2005. - 67 с.
