Анализ последствий аварийного разрушения сосудов с газами под давлением

АНАЛИЗ ПОСЛЕДСТВИЙ АВАРИЙНОГО РАЗРУШЕНИЯ СОСУДОВ С

ГАЗАМИ ПОД ДАВЛЕНИЕМ

И.Н. Пантелеев, доцент, к.ф.-м.н, ВГТУ, г. Воронеж

Рассмотрена задача об аварийном разрушении различных ёмкостей, предназначенных для хранения и перевозки сжатых или сжиженных газов. Проведенный анализ может использоваться для оценки опасности химических производств, построении карт риска, экспертных оценках последствий аварий на трубопроводах и хранилищах сжиженных газов.

В различных емкостях могут находиться инертные газы под давлением, сжатые углеводородные газы, перегретые водяные пары и т.д. Основной причиной разрушения емкостей, содержащих газы под давлением, является превышение внутренним давлением допустимого значения, вследствие нарушения нормальных условий эксплуатации емкостей, либо отказа устройств стравливания избыточного давления.

Возможны и иные причины взрыва сосуда с газом, например, истечение части газа вследствие внешнего пожара, образование небольшого облака газовоздушных смесей, взрыв или возгорание облака, разогрев сосуда с газом, повышение давления в сосуде и его полномасштабный взрыв [1-4].

При оценке последствий взрывов газовоздушных смесей в емкостях под давлением соотношения между параметрами смеси описываются на основе уравнения идеального газа:

р = Р-Т ,

где Р - давление в сосуде, Па, р - плотность газа, кг/м3, Т - абсолютная температура газа в сосуде, К, Яг - газовая постоянная, Дж/(г-К).

Я

* = мм.

Я = 8.314 Дж/(моль-К) - универсальная газовая постоянная, М - молярная масса газа, г/моль.

При использовании уравнения (2) для описания реальных процессов следует учитывать ограничение газообразного состояния среды: плотность р газа в сосуде не должна превышать значения плотности насыщенных паров при данной температуре Т.

Процессы, сопровождающие взрыв сосудов, содержащих газ под давлением, относятся к адиабатическим. Как известно, адиабатическим изменением состояния системы называется такое изменение, которое протекает без обмена теплом между системой и окружающими телами. При адиабатическом процессе система не получает извне тепла и не отдает его окружающей среде. Для адиабатического протекания процесса система должна быть окружена совершенно не теплопроводной средой. Поскольку подобное невозможно, то

156

модель любого реального процесса взрыва сосуда с газом под давлением может рассматриваться лишь как более или менее точное приближение к модели адиабатического процесса.

При взрыве металлической емкости, содержащей газ под давлением, образуются осколки, поражающее действие которых при таких взрывах зачастую бывает определяющим, а также образуется воздушная ударная волна и, возможно, интенсивное тепловое излучение. Энергия взрыва Ее при этом распределяется, в основном, между кинетической энергией Ек осколков, энергией ЕВ формирования воздушной ударной волны и тепловой энергией ЕТ:

е.- Ек - Ев - Ет = 0,

где Ее - суммарный эффект взрывного превращения и расширения энергоносителя (источника взрыва);

ЕВ - энергия, расходуемая на формирование ударной волны;

ЕТ - энергия, расходуемая на тепловое излучение;

ЕК - кинетическая энергия разрушения оболочки емкости и ее движения (сообщение осколкам начальной скорости полета).

Пренебрегая расходом энергии на разрушение оболочки и энергией, передаваемой осколками окружающему воздуху (его нагреванию осколками), для оценки сверху примем предположение, что вся доля энергии взрыва Ек расходуется только на разлет осколков:

—-К2 ^

= ек ,

где т - масса оболочки сосуда, содержащего газ, кг;

Уо - начальная скорость разлета осколков (принимается равной максимальной скорости в момент разлета), совпадающая со скоростью движения оболочки, м/с.

Исходя из уравнения (4), можно записать выражение для начальной скорости разлета осколков при взрыве сосуда, содержащего газ под давлением:

Ко .

V —

Поскольку скорость детонации при взрывах газовых смесей в резервуарах для продуктов детонации малой плотности равна:

р

О = у1 2 -

Ра

где Qв - удельная теплота взрыва единицы массы ГВС, то на величину начальной скорости полета осколков в случае взрыва сосуда, содержащего горючий газ, накладывается ограничение в виде:

- п У-1

- = Р<—Л ,

- у +1

где у - показатель адиабаты газа;

О - масса газа в резервуаре в момент аварии, кг.

у-1

При р > -— начальная скорость осколков

Л +1

у-1

При р < -— начальная скорость осколков

Л +1

к *

V. < К

4 Р

¡у2 -1

Б

у +1' полета

осколков

При QB = 0 (негорючие газы) начальная скорость определяется формулой (5).

При взрыве емкости под внутренним давлением Р инертного (негорючего) газа с плотностью в момент аварии р начальная скорость разлета осколков вычисляется по формуле (5), где значение энергии

Ек = Ее • ^ представляется в виде работы адиабатического расширения газа, сжатого избыточным давлением АР = Р - РА :

РК

Ек = К-

у-1

1 -

V Р у

, кДж,

где РА - атмосферное давление, кПа;

Р - давление газа в сосуде в момент его разрушения, кПа; ^ - доля потенциальной энергии, расходуемая на кинетическую энергию осколков, кДж;

у - показатель адиабаты газа;

V - внутренний объем резервуара в момент аварии, м3. При аварийном вскрытии резервуара со сжатым взрывоопасным газом значение энергии Ек = Ее • ^ , затрачиваемой на формирование кинетической энергии осколков определяется зависимостью:

Е = К

- у-1' \

Р Кв 1 - Г р) у

+— а

(у -1) V Р, у

,кДж,

где Gг = G •у • z - масса газа, участвующего во взрыве, кг;

QB - удельная теплота взрыва единицы массы ГВС, кДж/кг; V, z - коэффициенты, определяемые согласно (2) и (8). При расчетах параметров воздушной ударной волны используется часть суммарной энергии, определяемой по формуле:

определяется

Тротиловый эквивалент по соотношением:

Ев = Еъ •кв, кДж воздушной ударной волне

Ес

ш.

ТУВ

кг

о,9-а о,9-а где Qт - удельная теплота взрыва тротила, кДж/кг.

Максимальная дальность разлета осколков, при рассматриваемых авариях, может составить для резервуаров 1000.. .1500 м, для технологических установок и газгольдеров - до 3000 м.

Поражающее действие осколка определяется, в частности, его способностью пробить преграду. При высокоскоростном соударении

металлических тел принято считать, что сопротивление внедрению осколка в металл достаточно постоянно и зависимость глубины h внедрения осколка в металлическую преграду определяется функцией от квадрата скорости V встречи осколка с преградой.

Для оценки способности осколка вызвать воспламенение жидких углеводородных топлив, используется показатель Х в виде удельного импульса:

X = —^ = о-Ро-V.

50 7Г-Г0

Ударное действие осколка по конденсированному взрывчатому веществу, инициирующее детонацию, определяется критическим давлением и длительностью ударного импульса, точнее - притоком энергии в единицу времени, т.е. мощностью воздействия.

Для оценки этой величины её следует рассматривать как волновой процесс изменения сжимающих напряжений в веществе во времени. Экспериментальные граничные значения мощности воздействия осколка для возбуждения детонации прессованного тротила составляют безопасные мощности - не более 50 МВт; мощности, обеспечивающие возбуждение взрыва - не менее 300 МВт.

Список использованной литературы

1. Пантелеев И.Н., Пантелеев А.И. Анализ возгорания газовых облаков при аварийных выбросах дисперсных веществ в атмосферу. Материалы III Всероссийской научно-практической конференции. Воронеж: ВИ ГПС МЧС России, 2012. С.125-128.

2. Пантелеев И.Н., Пантелеев А.И. Моделирование распространения облаков тяжелых газов при техногенных авариях. Сб.ст. VIII Междунар. научн.-техн. конф. «Аналитические и численные методы моделирования естественнонаучных и социальных проблем».- Пенза, Изд-во ПГУ / 2013. С.184-188.

3. Пантелеев И.Н., Пантелеев А.И. Моделирование аварийного разрушения сосудов с газом под давлением. В сб. «Материалы XIII международного семинара ФММС», Воронеж, 2015. Ч.2. с. 123-128.

4. Пантелеев И.Н. Математическое моделирование экологического ущерба при выбросе загрязняющих веществ в атмосферу. Современные методы прикладной математики, теории управления и компьютерных технологий: сб. тр. IX междунар. конф. «ПМТУКТ-2016» / Воронеж: Изд-во «Научная книга», 2016. с. 261-264.