МЕТОД ОЦЕНКИ ЭФФЕКТИВНОСТИ ОГНЕЗАЩИТЫ СТАЛЬНЫХ КОНСТРУКЦИЙ НА ОБЪЕКТАХ НЕФТЕГАЗОВОГО КОМПЛЕКСА В УСЛОВИЯХ ОТКРЫТОГО ПОЖАРА
А.В. Федоров, доктор технических наук.
Санкт-Петербургский государственный университет информационных технологий, механики и оптики.
А.А. Кузьмин, кандидат педагогических наук, доцент; Н.Н. Романов, кандидат технических наук, доцент; Д.А. Минкин, кандидат технических наук, доцент. Санкт-Петербургский университет ГПС МЧС России
Описывается структура интегрированной диалоговой системы, включающей в себя электронную таблицу Microsoft Office Excel, семейство макросов на Visual Basic for Applications и программный комплекс ELCUT 6.3, которая позволяет на основе разветвленного алгоритма расчета теплового режима открытого пожара при разгерметизации резервуара с перегретой горючей жидкостью, а также газом под давлением или сжиженным газом на объекте нефтегазового комплекса, прогнозировать последствия воздействия интенсивного теплового излучения на стальные конструкции в зависимости от характеристик их огнезащиты.
Ключевые слова: нефтегазовый комплекс, стальные конструкции, огнезащита, тепловой режим, тепловое излучение, пожарная нагрузка, огненный шар, горючая жидкость, горючий газ, вспученный вермикулит
METHOD FOR ESTIMATION OF THE EFFICIENCY OF FIRE PROTECTION OF STEEL STRUCTURES ON OBJECTS OF OIL AND GAS COMPLEX IN CONDITIONS OF OPEN FIRE
A.V. Fedorov.
Saint-Petersburg state university of information technologies, mechanics and optics. A.A. Kuzmin; N.N. Romanov; D.A. Minkin.
Saint-Petersburg university of State fire service of EMERCOM of Russia
The article describes the structure of the integrated dialog system, including the Microsoft Office Excel spreadsheet, the Visual Basic for Applications macros family, and the ELCUT 6.3 software package, which allows, based on a branched algorithm for calculating the open fire's thermal mode for depressurizing a tank with an overheated combustible liquid, and As well as gas under pressure or liquefied gas at the oil and gas complex to predict the effects of intense heat radiation on steel structures, depending on From the characteristics of their fire protection.
Keywords: oil and gas complex, steel structures, fire protection, thermal regime, thermal radiation, fire load, fire ball, flammable liquid, combustible gas, expanded vermiculite
Развитие нефтегазового комплекса России в современных условиях сопровождается ростом числа аварийных ситуаций, сопровождающихся разрушением резервуаров и аппаратов, приводящих к возгоранию и даже взрыву горючих жидкостей (ГЖ), газов под давлением или сжиженных природных газов (СС111 ) [1]. Повышение
пожаровзрывобезопасности объектов нефтегазового комплекса становится важнейшей составной частью обеспечения защищенности государства от угроз техногенного и экологического характера.
Пожары на объектах нефтегазового комплекса, как правило, сопровождаются резким повышением температуры до 1 100 оС и более [2]. Для незащищенных строительных конструкций такое температурное воздействие может привести к утрате их несущей, ограждающей и теплоизолирующей способности.
Некоторые сооружения на объектах нефтегазового комплекса представляют собой легковозводимые сооружения, несущие конструкции которых изготавливаются, как правило, из углеродистой стали различных марок, которая отличается высокой теплопроводностью. Это приводит к тому, что в условиях пожара незащищенные стальные конструкции быстро прогреваются до температур, превышающих 400^500 оС. Под воздействием этих температур и нормативной нагрузки интенсивно развиваются температурные деформации, что ведет к быстрому разрушению всего сооружения (в пределах всего 0,12^0,25 ч). Повышение предела огнестойкости конструкций, то есть времени от начала воздействия пожара до наступления предельного состояния (для легковозводимых сооружений потери несущей способности стальных несущих конструкций) до требуемого уровня, достигается применением огнезащиты.
Применяемые для огнезащиты обычные (невспучивающиеся) материалы, как правило, занимают дополнительный объем и увеличивают за счет собственной массы нагрузку на защищаемую стальную конструкцию [3]. Кроме того, средства огнезащиты должны не только обеспечивать защиту стальной конструкции от внешнего огневого воздействия, но и обладать адгезией к подложке материала конструкции, долговечностью в нормальных условиях эксплуатации, технологичностью при изготовлении и нанесении на защищаемую стальную конструкцию. Этим требованиям соответствуют вспучивающиеся огнезащитные составы, огнезащитный эффект которых основан на образовании при тепловом воздействии пористой массы с низкой теплопроводностью, которая препятствует притоку тепла к защищаемой поверхности. Вспучивающиеся материалы являются одним из наиболее интенсивно развивающихся направлений в области разработки составов для огнезащиты стальных конструкций [4, 5].
Общие требования к огнезащитным материалам, предназначенным для несущих стальных строительных конструкций, изложены в ГОСТ Р 53295-2009 «Средства огнезащиты для стальных конструкций. Общие требования. Метод определения огнезащитной эффективности», однако данный документ не предлагает расчетных методов для предварительной оценки эффективности применяемых конструктивных решений.
При этом высокая эффективность, возможность использования индустриальных методов нанесения огнезащитных составов на поверхность конструкций обусловливает необходимость разработки расчетных методов, которые бы позволили оценивать эффективность огнезащиты стальных конструкций в условиях теплового излучения открытого пожара, а так же обосновывать рекомендации по выбору вида огнезащитного состава и его толщины в зависимости от характеристик пожара [6] на объектах нефтегазового комплекса.
Целью настоящей работы является разработка доступной и надежной методики проведения инженерных расчетов температурных полей для оценки эффективности огнезащиты стальных конструкций легковозводимых сооружений объектов нефтегазового комплекса в условиях пожара, возникающего при разгерметизации резервуара с перегретой ГЖ или газом под давлением, а также со сжиженными газами.
Для достижения поставленной цели необходимо:
- разработать тепловую модель для расчета температурных полей при прогреве незащищенных и защищенных стальных конструкций тепловым излучением с позиций надежности полученных решений и доступности их использования для оценки
эффективности огнезащиты в зависимости от теплофизических свойств применяемого огнезащитного материала, его толщины, условий развития пожара;
- провести программную реализацию алгоритма расчета теплового режима открытого пожара на объектах нефтегазового комплекса, который бы учитывал возможности как развития пожара при диффузионном режиме выгорания пожарной нагрузки, так и быстрого полного разрушения резервуара с последующим возникновением «огненного шара», оценкой продолжительности существования и характеристикой его теплового излучения.
Расчет металлических конструкций на огнестойкость в условиях интенсивного теплового излучения открытых пожаров на объектах нефтегазового комплекса сводится к решению трех совместных теплофизических задач:
- расчет теплового режима открытого пожара на объектах нефтегазового комплекса при диффузионном режиме выгорания пожарной нагрузки;
- расчет теплового режима открытого пожара на объектах нефтегазового комплекса при полном разрушении резервуара с последующим возникновением «огненного шара», оценкой продолжительности существования и характеристикой его теплового излучения;
- расчет распределения температуры по толщине защищенной металлической конструкции в зависимости от времени воздействия теплового излучения пожара.
Интенсивность теплового излучения q для пожара при проливе ГЖ или при горении твердых материалов вычисляется по формуле:
где Е/ - среднеповерхностная плотность теплового излучения пламени; уq - угловой коэффициент облученности; ф - коэффициент пропускания атмосферы.
Значение Е/ принимается на основе имеющихся экспериментальных данных, а при отсутствии таковых современные нормативы [6] допускают принимать величину Е/ равной: 100 кВт/м2 для сжиженных углеводородных газов и 40 кВт/м2 для нефтепродуктов.
Эффективный диаметр пролива с1 можно рассчитывать приведением формы пролива к кругу, при этом высота пламени Н, м, вычисляется по формуле:
Г
Н=42- а -
М
Л 0,61
V
Рв V 8
У
где М - удельная массовая скорость выгорания топлива, кг/(м2-с); рв - плотность окружающего воздуха, кг/м3; 8=9,81 м/-с2 - ускорение свободного падения.
Угловой коэффициент облученности определяется по формулам (1-3):
(1)
где ун - факторы облученности для вертикальной и горизонтальной площадок соответственно, определяемые с помощью уравнений (2) и (3):
¥ = — -
V п
г
8
А
\
V
Л2-
1У
агеН
8
)
А
Г
ТА2-
-агсгс
V'
'(А+1Х 8-1) (А-1)-(8+1)
л
У
(2)
1
Щ = -
н
л
(в -1/5) ( [(В+ТЦя^ТП (A -1/5) ( [(А+1)5-Т)
2
в -1
• агс^
V
(в -1)5 + О У л/
и2 -1
V
(а -1)5 +1
Л
У
где А=(^2+52+1)/(2 -5); в=(1+52)/(2- 5); 5=2г/й; И=2И/с1; г - расстояние от геометрического центра пролива до облучаемой конструкции.
Для случая теплового излучения пламени при горении разлившегося жидкого углеводородного топлива коэффициент пропускания атмосферы определяется по формуле:
ф=ехр [-7,010 4(г-0,5ф].
Интенсивность теплового излучения «огненного шара», который может сформироваться при быстром полном разрушении резервуара с горючей жидкостью, определяется на основе имеющихся экспериментальных данных, а при отсутствии таковых принимается значение Е/ равное 450 кВт/м2 [6].
Значение углового коэффициента облученности вычисляется по формуле:
И / Д + 0,5
ч
4-[(И / й, + 0,5)2 + (г / й, )2 ]
где Н - высота центра «огненного шара»; - эффективный диаметр «огненного шара»; г - расстояние от облучаемого объекта до точки на поверхности земли непосредственно под центром «огненного шара».
Эффективный диаметр «огненного шара» определяют по эмпирической формуле:
А=5,33да0,327,
где т - масса горючего вещества, кг, участвующего в образовании «огненного шара».
Величину высоты центра «огненного шара» принимаем равной половины его эффективного диаметра, а продолжительность (время его существования т5, с), определяется по эмпирической формуле:
г8=0,92т°,303.
Коэффициент пропускания атмосферы фш для случая «огненного шара» рассчитывается по формуле:
р = ехр[- 7,0-10-4 -(Vг2 + И2 - В5 /2)]. ш
Задача нахождения распределения температуры по толщине (х=0+5) защищенной металлической конструкции в зависимости от времени воздействия теплового потока при пожаре решается с помощью дифференциального уравнения теплопроводности с соответствующими граничными условиями второго или третьего рода.
Для пола, перекрытий, двутавровых балок и колонн при отсутствии в них внутренних источников теплоты может быть записано общее одномерное дифференциальное уравнение теплопроводности следующего вида:
дг ( д2г Ф дгл
— = а —2 +--
дт Vдx2 х дху
где х - координата тела, которой может быть и радиус цилиндра; Ф - постоянное число: для пластины Ф=0, для цилиндра (колонна) Ф=1.
В таком случае для пластины уравнение (4) принимает вид:
дг
= а-
д 2г
дт дх'
граничные условия для обогреваемой (х=0) и необогреваемой (х=^) поверхностей пластины:
,д г —л—
дх
Л —
с дх
х=0 =а/
1
—(ст)]
х=3 =а0
г ($,т)—г,
где г/, го - температура продуктов горения и температура воздуха со стороны
необогреваемой поверхности соответственно; а - коэффициент температуропроводности, зависящий от температуры материала ограждающей конструкции; X, - коэффициенты теплопроводности огнезащитного слоя и материала ограждающей конструкции; ау и ао - коэффициенты теплоотдачи на обогреваемой и необогреваемой поверхностях пластины соответственно.
Для цилиндра (колонны) уравнение теплопроводности (4) принимает вид:
дт
= а
гд 2г 1 дгл +
дг2 г дг
а граничное условие имеет вид:
лдг
дг
х=з=аг
у —кТ
Коэффициент теплоотдачи от продуктов горения к обогреваемой поверхности пластины вычисляется по формуле:
,, 0,0023 гг
а , = 11,63 е 1
Коэффициент теплоотдачи с необогреваемой поверхности к окружающему воздуху находится следующим образом:
- в случае если температура поверхности >60 0С, коэффициент теплоотдачи а0 вычисляют по формуле:
а0 = 11,63 е
0,0023 гг (З.т)
- в случае если температура охлаждаемой поверхности стенки воздухом <60 0С, то коэффициент теплоотдачи а0 вычисляют по формуле:
а0 = 4.07-^(¿тР^.
Теплофизические свойства материалов огнезащиты существенно отличаются от свойств конструкционной стали. Значения коэффициента теплопроводности различных огнезащитных покрытий в зависимости от температуры представлены на рис. 1 [7].
0 20 100 200 300 400
Температура, оС
Рис. 1. Зависимость коэффициента теплопроводности огнезащитной изоляции от температуры: 1 - каолиноперлитовая изоляция; 2 - асбестовермикулитовая изоляция; 3 - асбестоперлитовая изоляция; 4 - каолиновермикулитовая изоляция; 5 - пеноасбест
На рис. 2 представлен алгоритм расчета параметров теплового излучения открытого пожара при разгерметизации резервуара с перегретой ГЖ, а так же СПГ или горючим газом под давлением.
На основании представленного алгоритма была создана интегрированная диалоговая система, включающая в себя электронную таблицу Microsoft Office Excel, семейство макросов на Visual Basic for Applications и программный комплекс ELCUT 6.3.
Результаты расчета температурного поля защищенной металлической конструкции программный комплекс ELCUT 6.3 позволяет документировать в виде файла формата *.svg и может отображать в цветографической форме любым браузером (рис. 3).
Для подтверждения надежности результатов работы созданной специальной интегрированной диалоговой системы с ее помощью были выполнены расчеты по оценке огнестойкости стальных пластин с огнезащитным покрытием и проведено сравнение результатов расчетов с экспериментальными исследованиями [7].
Для исследования выбраны стальные пластины толщиной 3 мм, защищенные огнезащитным вермикулитовым составом, разработанным Ю.М. Тихоновым и регламентированным ВТУ-89 «Огнезащитное вермикулитовое покрытие ОВП-1 для стальных строительных конструкций»: вермикулит вспученный марки 150; фракций 0,6+5 мм - 80 %; вата каолиновая - 20 %; жидкое стекло плотностью 1,3 г/см3 - Т:Ж=1:2,5;
отвердитель (кремнефторид натрий) - 10 % от массы жидкого стекла; наполнитель - перлит вспученный - 15 % от массы жидкого стекла, фракций 0,16^0,3 мм.
Натурные испытания таких пластин проводились в отделе легких конструкций ЦНИИСК им. В.А. Кучеренко по температурно-временному режиму, рекомендованному ВНИИПО и регламентированному стандартом [8].
Для обеспечения аутентичности условий прогрева расчетная плотность теплового потока, падающего на обогреваемую поверхность, вызывала изменение температуры на ней, близкую к стандартной температурной кривой (рис. 4).
Разлив горючей жидкости
Исходные данные расчета: геометрические размеры резервуара, пожарная нагрузка, атмосферные условия
Огненный шар
Коэффициент пропускания атмосферы
Эффективный диаметр пролива
1
Высота пламени
1
Угловой коэффициент облученности пламени
1
Плотность теплового излучения пламени
Продолжительность выгорания розлива ГЖ
Массив свойств
воздуха
Л
Массив свойств продуктов горения нефтепродуктов
Угловой коэффициент облученности «огненного шара»
Массив характеристик обращаемых нефтепродуктов
Коэффициент пропускания атмосферы
Эффективный диаметр «огненного шара»
\
Высота «огненного шара»
\
Рис. 2. Алгоритм расчета параметров теплового излучения открытого пожара на объекте
нефтегазового комплекса
Температура Т[К]
Рис. 3. Температурное поле стальной колонны Я=160 мм через 20 мин с огнезащитным слоем (вермикулит) толщиной 10 мм
Рис. 4. Температура на поверхности огнезащищенной стальной пластины
Сравниваемые результаты расчета предела огнестойкости и результаты натурных испытаний стальных пластин, огнезащищенных вермикулитовым покрытием ОВП-1 [7], представлены в таблице.
Таблица. Результаты натурных испытаний на огнестойкость и расчета предела огнестойкости стальных пластин
Толщина, мм Средняя плотность покрытия до теплового воздействия, кг/м3 Предел огнестойкости, мин
натурные испытания результаты расчета
металла покрытия
3 12 480 43 47
3 15 400 48 51
3 25 505 93 87
3 30 483 86 78
Полученные расхождения в результатах расчета и натурных испытаний предела огнестойкости стальных пластин, защищенных вермикулитовым покрытием, не превышают 10 %.
Таким образом, в ходе проведенной работы были получены следующие результаты:
- предложена методика расчета теплового режима открытого пожара на объектах нефтегазового комплекса, которая учитывает возможности как развития пожара при диффузионном режиме выгорания пожарной нагрузки, так и быстрого полного разрушения резервуара с последующим возникновением «огненного шара», оценивает продолжительность существования и определяет характеристики его теплового излучения; разработан алгоритм расчета тепловых полей незащищенных и защищенных стальных конструкций, учитывающий зависимость температурного поля от изменения теплофизических свойств материала конструкций и применяемого огнезащитного материала по мере развития открытого пожара;
- сформирована интегрированная диалоговая система Microsoft Office Excel - макрос Visual Basic for Applications - программный комплекс ELCUT 6.3, реализующая этот алгоритм с интерфейсом, дружественным широкому кругу пользователей;
- исследована надежность работы сформированной интегрированной диалоговой системы путем сравнения результатов расчета предела огнестойкости и результатов натурных испытаний стальных пластин, огнезащищенных вермикулитовым покрытием ОВП-1.
Сравнительный анализ полученных результатов расчета и экспериментальных исследований подтвердил возможность использования программного продукта для предварительной оценки эффективности огнезащиты стальных конструкций на объектах нефтегазового комплекса.
Литература
1. Информация об авариях, произошедших на предприятиях, подконтрольных территориальным органам Федеральной службы по экологическому, технологическому, атомному надзору. URL: http://www.rostehnadzor.ru/chronicle.html. (дата обращения: 11.04.2017).
2. Novozhilov V. Non-Linear Dynamical Model of Compartment Fire Flashover // Journal of Engineering Mathematics. 2010. Vol. 67. No 4. P. 387-400.
3. Архитектурное материаловедение / Ю.М. Тихонов [и др.]: учеб. СПб.: Изд-во Академия, 2013. 288 с.
4. Страхов В.Л., Крутов А.М., Давыдкин Н.Ф. Огнезащита строительных конструкций / под ред. Ю.А. Кошмарова. М.: Информ.-изд. центр «ТИМР», 2000. 433 с.
5. Драйздел Д. Введение в динамику пожаров / под ред. Ю.А. Кошмарова, В.Е. Макарова; пер. с англ. М.: Стройиздат, 1990. 424 с.
6. НПБ 105-03. Определение категорий помещений, зданий и наружных установок по взрывопожарной и пожарной опасности (утв. Приказом МЧС РФ от 18 июля 2003 г. № 314). Доступ из справ.-правовой системы «Гарант».
7. Исследование поведения огне-теплозащитных материалов в условиях пожара / Ю.М. Тихонов [и др.]. Ч. I: Исследования поведения огнезащитных материалов на основе вермикулита, перлита и стеклофибры в условиях пожара: монография / под общ. ред. В.С. Артамонова. СПб.: Астерион, 2015. 235 с.
8. ГОСТ 30247.1-94. Межгосударственный стандарт. Конструкции строительные. Методы испытаний на огнестойкость. Несущие и ограждающие конструкции. М.: Изд-во стандартов, 1995.
References
1. Informaciya ob avariyah, proizoshedshih na predpriyatiyah, podkontrol'nyh territorial'nym organam Federal'noj sluzhby po ehkologicheskomu, tekhnologicheskomu, atomnomu nadzoru. URL: http://www.rostehnadzor.ru/chronicle.html. (data obrashcheniya: 11.04.2017).
2. Novozhilov V. Non-Linear Dynamical Model of Compartment Fire Flashover // Journal of Engineering Mathematics. 2010. Vol. 67. No 4. P. 387-400.
3. Arhitekturnoe materialovedenie / Yu.M. Tihonov [i dr.]: ucheb. SPb.: Izd-vo Akademiya, 2013. 288 s.
4. Strahov V.L., Krutov A.M., Davydkin N.F. Ognezashchita stroitel'nyh konstrukcij / pod red. Yu.A. Koshmarova. M.: Inform.-izd. centr «TIMR», 2000. 433 s.
5. Drajzdel D. Vvedenie v dinamiku pozharov / pod red. Yu.A. Koshmarova, V.E. Makarova; per. s angl. M.: Strojizdat, 1990. 424 s.
6. NPB 105-03. Opredelenie kategorij pomeshchenij, zdanij i naruzhnyh ustanovok po vzryvopozharnoj i pozharnoj opasnosti (utv. Prikazom MCHS RF ot 18 iyulya 2003 g. № 314). Dostup iz sprav.-pravovoj sistemy «Garant».
7. Issledovanie povedeniya ogne-teplozashchitnyh materialov v usloviyah pozhara / Yu.M. Tihonov [i dr.]. CH. I: Issledovaniya povedeniya ognezashchitnyh materialov na osnove vermikulita, perlita i steklofibry v usloviyah pozhara: monografiya / pod obshch. red. V.S. Artamonova. SPb.: Asterion, 2015. 235 s.
8. GOST 30247.1-94. Mezhgosudarstvennyj standart. Konstrukcii stroitel'nye. Metody ispytanij na ognestojkost'. Nesushchie i ograzhdayushchie konstrukcii. M.: Izd-vo standartov, 1995.