CC BY
97
УДК 004.021
О. С. Колясев, А. М. Сверчков *, Т. В. Савицкая
Российский химико-технологический университет им. Д.И. Менделеева, Москва, Россия 125480, Москва, ул. Героев Панфиловцев, д. 20, корп. 1 *e-mail: andrey_sverchkov@mail.ru
МОДЕЛИРОВАНИЕ ПОСЛЕДСТВИЙ АВАРИЙ НА МАГИСТРАЛЬНЫХ ТРУБОПРОВОДАХ С ИСПОЛЬЗОВАНИЕМ ПРОГРАММНОГО КОМПЛЕКСА ТОХ1+М8К
Проанализирован ряд методик по расчету интенсивности истечения газа при разрушении газопровода, расчету геометрических размеров пламени котлованного и струевого типов и расчету теплового потока с поверхности пламени. Проведен анализ результатов, полученных при использовании модуля "Струевое горение газа" входящего в состав программного комплекса TOXI+Risk.
Ключевые слова: опасный производственный объект, магистральные газопроводы, пламя колонного типа, струевое пламя, пробит-функция.
Эксплуатация опасных производственных объектов (ОПО), в том числе объектов систем газораспределения и газопотребления, связана с высоким риском возникновения аварий, экономические потери от которых составляют десятки миллионов рублей. Так, например, за последние 5 лет, по данным Федеральной службы по экологическому, технологическому и атомному надзору на объектах газораспределения и газопотребления произошло 204 аварии, общий ущерб от которых составил порядка 2,6 млрд. рублей [1].
Наиболее распространенным фактором антропогенного воздействия, приводящего к отрицательным последствиям, является загрязнение окружающей среды и, в частности, перенос вредных веществ на большие расстояния.
В целях повышения уровня промышленной безопасности ОПО транспортировки газа ООО «Научно-исследовательский институт природных газов и газовых технологий - Газпром ВНИИГАЗ» разработало методические указания по проведению анализа риска для ОПО газотранспортных предприятий [2].
Данный документ содержит подробное описание последовательных этапов
количественного анализа риска применительно к технологическим объектам транспорта газа. Также он содержит ряд приложений, включающих расчетные методики, например:
- Методика расчета интенсивности истечения и количества выбрасываемого газа при разрушениях газопроводов;
- Методика расчета пространственно-временного распределения выбрасываемых продуктов с учетом их физико-химической трансформации;
- Методика расчета физических параметров и масштабов распространения поражающих факторов аварий и др.
На основе вышеуказанных приложений в закрытом акционерном обществе «Научно-технический центр исследований проблем промышленной безопасности» (ЗАО НТЦ ПБ) разработан ряд программных модулей, включенных в состав программного комплекса TOXI+Risk [3]. Он используется на кафедре компьютерно-интегрированных систем в химической технологии в учебном процессе и предоставлен на безвозмездной основе ЗАО НТЦ ПБ [4, 5].
В программном комплексе TOXI+Risk начиная с версии 4.3 реализован модуль «Струевое горение газа». Рассмотрим его работу и проанализируем результаты расчета. В данном модуле реализованы несколько методик [2]:
1. Методика расчета интенсивности истечения газа при разгерметизации однониточного газопровода;
2. Методики расчета геометрических размеров пламени котлованного и струевого типов;
3. Методики расчета теплового потока с поверхности пламени (пламя колонного типа в виде вертикального цилиндра и пламя струевого типа в виде лежащих полуцилиндров).
Самой трудоемкой, в плане реализации расчетов и получаемых форм изолиний интенсивности теплового излучения, методикой является расчет теплового потока с поверхности пламени в виде двух лежащих полуцилиндров.
Рассмотрим последовательность данных вычислений. Расчет угловых коэффициентов излучения проводится для точек поверхности земли расположенных на лучах Л1 - Л6 от боковой поверхности панели, на лучах Лс1 и Лс2 от торцевых поверхностей полуцилиндра (рис. 1). Лучи Л1, Л3, Л4 и Л6 проводится по нормали к оси цилиндра через боковую поверхность цилиндра, лучи Л2 и Л5 - по нормали к оси цилиндра на половине его длины. Лучи ЛС1 и ЛС2 проводятся по оси цилиндра в направлении от торцевых
поверхностей, начало этих лучей лежит на торцевых поверхностях цилиндра.
Для лучей Ль..Лб угловой коэффициент излучения от лежащего полуцилиндра для вертикально ориентированных площадок-приемников излучения фв и для горизонтального расположения площадок фг, расположенных на поверхности грунта, рассчитывается по формулам:
*=гЬщ ■ И* (т^г) - ■ (ЦЩ)+
, ¿1 <Х1+21) ^ Н--. ■ аг&д
2- п ■ Нг
2 ■ \ ■
агад ^ ^ - 1г ■ аггЛд
12-(Х2 + 72) / |ся3- 1)-А'2
Н--. ■ аг^д
и (Л2+1)Г
Д, + 1
+Н{-1
г_л / |(Я2 1)■ Х2
• агс1д
12 +Н21
(Я2+ 1) 2г I У
ЭКЕ ЭКЕ
Х2 = (Яг + 1)г + 1г = (Яг - 1)г + 1\.
где Яэкв - эффективный диаметр факела, м; х -расстояние по горизонтали от оси цилиндра до места расположения приемника излучения, м; Ьф1 - длина полуцилиндра пламени с одной стороны от луча, проведенного по нормали к оси полуцилиндра, м; Ьф2 - длина полуцилиндра пламени с другой стороны от луча, проведенного по нормали к оси полуцилиндра, м. Для лучей ЛС1 и ЛС2 расчет угловых коэффициентов излучения при х > Яэкв проводится по формуле:
Ч>
= п1'т*Я\—Гх^ж)
В результате данных расчетов определяют так называемый геометрический фактор. Умножая его на коэффициент поглощения излучения атмосферой и на интенсивность излучения с единицы поверхности пламени, получаем искомую величину - радиационное тепловое
воздействие пожаров на газопроводах на прилегающие объекты.
Первым действием при работе с модулем (рис. 2) является выбор типа пламени, который будем рассматривать (пламя колонного типа или струевое пламя). Далее необходимо выбрать вещество, после чего произойдет автоматическое заполнение всех необходимых для расчета величин. Стоит отметить, что вещества берутся из общей для всех модулей базы данных веществ. Данная база имеет несколько десятков самых распространенных веществ. Пользователь может вносить изменения в данную базу, как в существующие записи, так и добавлять новые вещества. В данном модуле поля ввода для свойств вещества являются активными, что позволяет, при необходимости, изменить те или иные свойства. Далее необходимо задать диаметр трубопровода и интенсивность выброса газа при аварии, причем, если выбрано пламя колонного типа, то значение интенсивности выброса газа будет одно. Если же выбрано струевое пламя, то значений интенсивности выброса газа будет уже два. Одно для нагнетающего участка трубопровода, второе для всасывающего участка трубопровода. Для задач определения расстояния от точки выброса по заданной интенсивности теплового излучения и определения интенсивности теплового излучения в заданной точке пространства никаких дополнительных исходных данных больше не нужно, можно приступать к расчету.
Для задач определения расстояния от точки выброса по заданной вероятности поражения тепловым излучением и определения вероятности поражения тепловым излучением в заданной точке пространства необходимо выбрать методику, пробит-функция из которой будет применяться в расчетах. Если пользователю не известно значение интенсивности выброса газа при аварии, то рядом с соответствующим полем ввода есть кнопка с изображением калькулятора, при нажатии на которую откроется окно определения интенсивности выброса со своим перечнем исходных данных.
После заполнения всех необходимых исходных данных в левой части программного модуля перейдем к правой. Существует шесть детерминированных критериев, рекомендуемых [б, 7], численные значения интенсивности излучения которых равны, например: 1400Вт/м2 -без негативных последствий; 4200Вт/м2 -безопасно для человека в брезентовой одежде; 13900Вт/м2 - воспламенение древесины и др. Пользователю предоставляется возможность расчета неограниченного количества
интересующих его значений в рамках одной задачи. Для этого необходимо нажать на знак "+" в соответствующем типе расчетов и заполнить критерий расчета.
Рис. 1. Модель излучателя в виде двух лежащих полуцилиндров
-Г Струе вого горение газа ВНЕ
Тип пламени при аварии (•" Пламя колонного типа С Струевое пламя Г" В - Без негативны* последствий 14ш Ву/м2 _ 15зэ.п4 [— _ |Безопаснодля человека в брезентовой одежде 4200 Вт/м2 — 386 67 м »
Вещество
Метан -w Г" ш - |Нв™Рв"ссимая 6оль чвевз 20 с«к 7000 Вт/м2 - 673.1 м
Молярная масса
0,016 кг/моль - „ Непереносимая боль через 3-5 сек , 0500 Вт/м2 , 531 85
Показатель адиабаты Г" Н - | Воспламенение древесины , 3300 Вг/„2 - 445.78
¡1.42 Низшая теплота сгорания Г" ш - роспламанани. реаины 14в00 ¡Вт/м2 - 427.83 м
[50100000 Дж/кг . ^ |. | Определение интенсивности теплового излучения в заданной точке пространства
Температура транспортировки ¡283 К |Пользовагельскии критерии 5(ю . м . 11605.78 ВтЛ,2 -
Плотность при НУ . + I — I Определение расстояния от точки выброса по заданной интенсивности теплового излучения
0,71 Б8 кг/мЗ Г" ■ - IПользовательский критерии 00 Вт,м2 ^ 5447,56 м
Параметры разгерметизации Интенсивность выброса газа при аварии + _ Определение расстояния от точки выброса по заданной вероятности поражения тепловым излучением
|— щ ^ | Смертельное поражение эо ^ 742,27 м
9111,78 кг/с - h Mil Внутренний диаметр M Г 1400 мм , ^ | | Определение вероятности поражения тепловым излучением в заданной точке пространства
\ Смертельное поражение 300 м -^100 Ч -
П робит-функция Г" ГОСТ 12.3.047-98 Методика определения расч. величин пожарного риска на ПО,2009 СТО Газпром 2-2.3-400-2009
Рис. 2. Главное окно программ
Рассмотрим результаты расчета на конкретном примере, исходными данными которого являются: внутренний диаметр трубопровода,
транспортирующего природный газ - 1400мм; температура транспортировки - 283К; рабочее давление - 7,5МПа; тип труб - стальные новые; время отключения компрессорной станции - 270с; время перекрытия задвижек - 360с; расстояние от места разрыва до задвижки по ходу движения газа и против - по 10000м; расстояние от места разрыва до компрессорной станции по ходу движения газа и против - по 60000м.
В ходе моделирования аварийной ситуации (гильотинного разрыва газопровода) проведены
го модуля "Струевое горение газа"
расчеты для различных видов пламени и для различных видов пробит-функций, на основе которых определялись вероятности термического поражения. Кроме того, были проведены исследования для конкретной удаленной от эпицентра аварии точки для разных сценариев. Результаты показали, что интенсивность излучения на расстоянии в 500 м для пламени колонного типа - 11605,8 Вт/м2, а для струевого пламени в направлении, перпендикулярном оси струи для нагнетающего отрезка трубопровода -7762,9 Вт/м2. Сценарий с пламенем колонного типа для конкретной удаленной точки является более опасным, чем сценарий со струевым
пламенем, т. к. в первом случае массовый расход газа суммируется для обоих концов разрушенного трубопровода.
Расчеты также показывают, что одна и та же вероятность поражения, например 90%, достигается на расстоянии: 447,8 м по методике [6]; 559,62 м по методике [7] и 712,27 м в соответствии с [2]. Таким образом на результаты сильное влияние оказывает вид пробит-функции.
Результаты расчета зон поражения тепловым излучением, полученные из данного модуля, участвуют в построении полей потенциального риска в программном комплексе TOXI+Risk, так как методология риск-анализа является эффективным инструментом управления показателями безопасности сложных
технологических объектов.
Колясев Олег Сергеевич, студент факультета информационных технологий и управления РХТУ им. Д. И. Менделеева, Россия, Москва
Сверчков Андрей Михайлович, аспирант кафедры КИС ХТРХТУ им. Д. И. Менделеева, Россия, Москва
Савицкая Татьяна Вадимовна, д.т.н., профессор кафедры КИС ХТ РХТУ им. Д. И. Менделеева, Россия, Москва
Литература
1. Пантюхова Ю.В. Методика оценки уровня промышленной безопасности опасных производственных объектов систем газораспределения и газопотребления: дис. канд. техн. наук. - М, 2011. - 127 с.
2. Методические указания по проведению анализа риска для опасных производственных объектов газотранспортных предприятий ОАО «Газпром»// Газпром экспо. - 2009. - 310 с.
3. Агапов А.А., Лазукина И.О., Марухленко А.Л. и др. Использование программного комплекса ТОКСИ+Risk для оценки пожарного риска// Безопасность труда в промышленности. - 2010. - № 1. - С. 46-52.
4. Егоров А.Ф., Савицкая Т.В., Агапов А.А. Опыт использования программного комплекса ТОКСИ+Risk для подготовки специалистов// Безопасность труда в промышленности. - 2012. -№ 5. - С. 22-25.
5. Савицкая Т.В., Егоров А.Ф., Запасная Л.А. и др. Сравнительный анализ результатов моделирования последствий химических аварий с использованием программного комплекса ТОКСИ+Risk// Безопасность труда в промышленности. - 2012. - № 8. - С. 78-83.
6. ГОСТ Р 12.3.047-98. Система стандартов безопасности труда. Пожарная безопасность технологических процессов. Общие требования. Методы контроля.// ИПК Издательство стандартов. - 1998. - 85 с.
7. Методика определения расчетных величин пожарного риска на производственных объектах (утв. приказом Министерства Российской Федерации по делам гражданской обороны, чрезвычайным ситуациям и ликвидации последствий стихийных бедствий 10.07.09 № 404). -2009. - 63 с.
Kolyasev OlegSergeevich, Sverchkov Andrey Mikhaylovich*, Savitskaya Tatyana Vadimovna
D.I. Mendeleev University of Chemical Technology of Russia, Moscow, Russia. *e-mail: andrey_sverchkov@mail.ru
MODELING THE CONSEQUENCES OF ACCIDENTS ON TRUNK GAS PIPELINES USING SOFTWARE, TOXI+RISK
Abstract
Analyzed a number of methods for the calculation of the intensity of the gas flow in the destruction of the pipeline, the calculation of geometrical sizes of flame of columned and of jet types and calculation of the heat flux from the surface of the flame. The analysis of the results obtained using the module "jet-burning gas" is part of the software package TOXI+Risk.
Key words: Hazardous production facility, the trunk gas pipelines, flame of columned type, flame of jet type, probit-func
