Контакты
Подписка
МЕНЮ
Контакты
Подписка

Оценка мощности взрывов газопаровоздушных смесей при аварийных проливах сжиженного природного газа

Оценка мощности взрывов газопаровоздушных смесей при аварийных проливах сжиженного природного газа

В рубрику "Пожарная безопасность" | К списку рубрик  |  К списку авторов  |  К списку публикаций

Оценка мощности взрывов газопаровоздушных смесей при аварийных проливах сжиженного природного газа

Перспективы развития топливо-энергетического комплекса (ТЭК) Российской Федерации, выход добывающих компаний России в шельфовую зону Арктики, стратегическая задача упрочения позиций России как энергетической сверхдержавы предъявляют повышенные требования к обеспечению безопасности и надежности функционирования объектов ТЭК. Основная опасность, приводящая к чрезвычайным ситуациям на объектах ТЭК, связана со взрывами газопаровоздушных смесей (ГПВС), возникающими вследствие аварийного пролива СУГ или СПГ. На основе прикидочных расчетов можно оценить возможную мощность взрыва при внезапном разрушении резервуара, содержащего до 100 тыс. м3 СУГ или СПГ. По мнению авторов, мощность такого взрыва эквивалентна десятку взрывов, осуществленных США в Хиросиме и Нагасаки [1]. Учитывая важность поднятой проблемы и необходимость разработки обоснованных решений повышенной пожаровзрывобезопасности объектов ТЭК, разработка методики оценки мощности взрывов ГПВС в открытом пространстве особенно актуальна
Анатолий
Гришкевич
Старший научный сотрудник ФГБУ "27 НЦ" МО РФ, к.т.н.
Владислав
Филин
Доцент ФГБУ "27 НЦ" МО РФ, к.т.н.
Виктор
Ушаков
Старший научный сотрудник ФГБУ "27 НЦ" МО РФ, к.т.н.
Геннадий Маньковский
Профессор, ФГБУ "27 НЦ" МО РФ, д.т.н.

Общие положения

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

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

Сферическая детонационная волна может возникнуть и непосредственно в ГПВС от слабого энергетического источника, например от искры, если размер облака превышает некоторое критическое значение и в пределах концентраций, указанных в табл. 1 [2].


Характеристики взрыва газопаровоздушных смесей приведены в табл. 2.


Для расчета поражающего действия взрывов ГПВС определяют параметры детонационной волны расширяющихся продуктов детонации (ПД) и воздушной ударной волны (ВУВ).

Начальный радиус r0, м, полусферического облака ГПВС в зависимости от его объема V0, м, определяют по формуле

Поскольку согласно закону Авогадро киломоль µГ идеального газа занимает объем Vа=22,4 м3, для газовой смеси стехиометрического состава концентрацией ССТХ и молекулярной массой µГ горючей компоненты массы Св запишем приближенную формулу для объема облака

Значение коэффициента ϑ принимают в зависимости от способа хранения продукта: ϑ= 1 - для газов при атмосферном давлении; ϑ=0,5 - для газов, сжиженных под давлением; ϑ=0,1 - для газов, сжиженных охлаждением; ϑ=0,02… 0,07 - при растекании легковоспламеняющихся жидкостей. При взрыве детонационная волна распространяется со скоростью D, м/с, (Qm - удельная теплота взрыва, Дж/кг)

то есть время полной детонации облака tв, с, равно

Максимум избыточного давления на фронте детонационной волны (химический пик)

Через весьма малый интервал времени τ давление в детонационной волне уменьшается в два раза (так называемая точка Жуге)

Характеристики распространенных газо- и паровоздушных смесей приведены в табл. 2, где приняты обозначения: ΔP2 - избыточное давление детонационной волны (эффективное давление), γСТХ - показатель адиабаты продуктов детонации, - плотность, Qm.СТХ и Qv.СТХ – удельная теплота взрыва по массе и объему смеси соответственно, а индексом СТХ помечен стехиометрический состав. Параметры ударной волны на расстояниях, превышающих размер облака R0 ≥ r0, м, определим по формулам, аппроксимирующим численное решение задачи о детонации пропановоздушной смеси (Б.Е. Гельфанд и др., 1985 [2, 8, 9]). Полученные универсальные зависимости максимума эффективного избыточного давления ΔPm, Па, и удельного импульса I, Па.с, в ударных волнах от расстояния R от центра взрыва в энергетических координатах удовлетворительно согласуются с экспериментальными данными для горючих смесей углеводородов с воздухом или кислородом, независимо от состава смеси:

Здесь mT - тротиловый эквивалент наземного взрыва полусферического облака ГПВС, кг,

где:

m=ρСТХ•V0 - масса горючего облака, кг;
QT - удельная теплота взрыва тротила (принята 4,184 МДж/кг);
P0 - атмосферное давление, Па, (P0=1,013 x 105 Па).

Для сшивания параметров на границе облака величина r0 заменяется на близкое значение R0, определяемое из условия непрерывности функции ΔPm(R) в точке:

Эффективное время θ, с, действия фазы сжатия ударной волны с давлением, аппроксимируемым треугольником, определяется по формуле

В качестве критерия мощности поражающего действия (МПД) объемно-детонирующего взрыва выбрана площадь приведенной зоны поражения человека с учетом степени его защищенности и времени потери работоспособности. К объемно-детонирующим взрывам относятся взрывы ГПВС или облака топливо-воздушной смеси (ТВС), поражающее действие которых обусловлено действием воздушной ударной волны (ВУВ), генерируемой при детонации облака ТВС, образующегося при диспергировании горючего (как правило, углеводородного) в аэрозольное или газообразное состояние. За рубежом такие взрывчатые системы носят название FAE (Fuel Air Explosive). Параметры ВУВ определяются энергетическими и физико-химическими характеристиками горючего, условиями формирования топливо-воздушного облака у цели и возбуждения в нем взрывного процесса в детонационном режиме. Реализация энергии аэрозольного заряда может осуществляться по двум основным технологическим схемам:

  • по двухтактной схеме, где сначала формируется топливо-воздушное облако, а затем через некоторое время в нем возбуждается с помощью инициирующего заряда реакция химического превращения топлива в режиме детонации (объемно-детонирующий взрыв - ОДВ);
  • по однотактной схеме, в которой реакция химического превращения топлива протекает под действием ударной волны (УВ), образованной диспергирующим зарядом (объемный взрыв).

Объемно-детонирующие взрывы наиболее часто происходят при авариях с веществами, относящимся к различным классам органических соединений: непредельным углеводородам (1,3 - пентадиен); органическим окисям (окись пропилена); циклическим соединениям; рецептуры пастообразных композиций, в основу которых входят углеводороды с добавками порошков металлического горючего и промотирующие добавки. Проведенные исследования позволили выявить характерные черты газового взрыва, основными из которых являются: малая доля энергии, уходящая в воздушную УВ (примерно половина от полной выделившейся энергии при детонационном взрыве); очень резкий спад давления вблизи заряда при взрыве в режиме детонации; формирование волны сжатия при сгорании зарядов, которая по мере распространения трансформируется в ударную волну.

Обработка полученных опытных данных позволила авторам [3] предложить эмпирические зависимости для основных параметров воздушной УВ при детонации сферических зарядов топливо-воздушных смесей на основе метода энергетического подобия, который согласуется с вышеизложенным подходом к оценке параметров ВУВ (Б.Е. Гельфанд и др. [8, 9])

Здесь:

E0 - полная энергия взрыва в ккал;
r - расстояние от центра заряда в м;
ΔPm - максимальное избыточное давление на фронте УВ в атм.;
i - удельный импульс избыточного давления фазы сжатия в кгс-с/м2;
τ - длительность фазы сжатия в мс.

Время задержки подрыва облака ТВС определяется временем испарения основной массы аэрозоля горючего, для больших емкостей, например баллонов, содержащих более 100-200 кг горючего, это время может составлять 80-120 мс.

Отличительной особенностью объемного взрыва ТВС является высокое поражающее действие по людям в укрытиях и сооружениях за счет высокого избыточного давления в облаке ОДС и его сплошности, а также способности его "затекать" внутрь негерметичных сооружений

Рассмотрим основные закономерности поражения человека в результате воздействия ПФ ВУВ при взрыве объемно-детонирующих систем.

Закономерности поражения человека ВУВ при взрыве объемно-детонирующих систем

Закон поражения, выражающий зависимость вероятности поражения цели от удаления точки разрыва от цели, называют координатным законом поражения (КЗП). Универсальный подход к оценке поражающего действия ВУВ, основанный на результатах исследований последних лет, был предложен специалистами ВАА им. М.И. Калинина и ВМедА им. С.М. Кирова. Данный подход заключается в том, что в качестве основной характеристики ВУВ, определяющей эффект первичного поражения человека, принимается эффективная удельная энергия ВУВ, отнесенная к единице поверхности и представляющая собой суммарную удельную избыточную механическую энергию воздушных потоков, пришедших к цели вслед за фронтом ВУВ [4, 5].

Значение среднепоражающего давления на фронте ВУВ зависит от длительности фазы сжатия, особенно сильно при длительностях до 10 мс. Поэтому наиболее общей и универсальной характеристикой ВУВ является удельная энергия взрыва, представляющая собой суммарную удельную избыточную механическую энергию воздушных потоков, пришедших вслед за фронтом ВУВ. Удельная энергия определяется по формуле [4]

где:

Eуд - энергия взрыва, Дж/м2;
P(t) - избыточное давление, Па;
ρ(t) - плотность воздуха в зависимости от времени, кг/м3;
u(t) - массовая скорость воздуха, м/с;
t - длительность фазы сжатия ВУВ, с;
k - показатель ударной адиабаты воздуха (k=1,4).

Учитывая отсутствие поражения при действии подпороговых значений удельной энергии ВУВ и зависимость среднепоражающего значения удельной энергии ВУВ от длительности фазы сжатия ВУВ, был доработан параметрический закон поражения. Полученные зависимости позволяют оценивать вероятность поражения человека по 5 степеням тяжести в зависимости от величины воздействующей эффективной удельной энергии и длительности фазы сжатия ВУВ. Данные зависимости имеют следующий вид [11]:

где:

Pусл – условная вероятность поражения по i-типу поражения;
kэ - эмпирический коэффициент ориентации человека;
ko – коэффициент ослабления энергии ударной волны;
a, b – коэффициенты, зависящие от степени тяжести поражения;
t+ – длительность фазы сжатия, мс;
Eэф.уд – эффективная удельная энергия ВУВ, кДж/м2;
Ei0,5 – эффективная удельная энергия ВУВ, вызывающая поражение i-й степени тяжести, с вероятностью, равной 0,5, кДж/м2;
Ei0,5=fi при 2 ≤ t+ ≤ 5 мс ;
Ei0,5=ci+di•t+ при 2 ≤ t+ ≤5 мс .

График трехмерной функции поражения живой силы в зависимости от величины удельной эффективной энергии ВУВ и длительности фазы сжатия приведен на рис. 1.


Коэффициенты для расчета различных степеней тяжести поражения приведены в табл. 3.


При однократном воздействии ВУВ потери пострадавших определяются удельной энергией этой ВУВ. При одновременном и многократном воздействии ударных волн различной интенсивности нанесение ущерба людям находят через эффективную суммарную удельную энергию всех воздействующих волн.

Максимальный интервал времени учета воздействия ПФ ВУВ не должен превышать периода морфофизиологического, психофизиологического и психического восстановления организма.

Параметры воздушной ударной волны объемно-детонирующего взрыва

Основными параметрами ВУВ, определяющими ее поражающее действие на человека и разрушающее по материальным объектам, являются: избыточное давление на фронте ВУВ, импульс положительной фазы сжатия, длительность фазы сжатия и эффективная удельная энергия ВУВ. Профили изменения основных параметров ВУВ при взрыве объемно-детонирующих систем отличаются от аналогичных зависимостей для случая конденсированных взрывчатых веществ (КВВ), описываемых по формулам М.А. Садовского и ГОСТ В-25.801-83 [7]. Зависимости параметров ВУВ для случая наземного взрыва объемно-детонирующих систем были получены на основе метода энергетического подобия взрыва и теории размерностей [1, 7, 8], численные значения коэффициентов в уравнениях получены методом регрессионного анализа многочисленных экспериментальных данных отечественных и зарубежных исследователей, в которых масса заряда горючего изменялась в пределах от 3 до 1000 кг. Получены зависимости основных параметров ВУВ (избыточного давления на фронте, импульса давления фазы сжатия и длительности фазы сжатия) от расстояния. Данные выражения построены в виде преобразований Сахса и имеют следующий вид:

где:

ΔPф - избыточное давление на фронте ВУВ, Па;
P0 - атмосферное давление, Па;
a0 - скорость звука в невозмущенной атмосфере (a0 = 340 м/с), м/с;
Е - полная энергия взрыва облака ТВС (Е=МгQвзр), Дж;
Мг - масса горючего, кг;
Qвзр - удельная теплота взрыва горючего, Дж/кг;
τ+ - длительность фазы сжатия, с;
I+ - импульс давления фазы сжатия, Па.с;
R - расстояние от центра облака, м.

Значения эффективной удельной энергии ВУВ (в кДж/м2) рассчитываются по уравнению (5) по значению избыточного давления на фронте ВУВ (в изб. атм) и длительности фазы сжатия (в мс) со средней относительной ошибкой, не превышающей 1% (по сравнению с результатами численного интегрирования):

Поражение человека рассчитывается через приведенную площадь поражения. Принимая допущение, что в точках, равноудаленных от центра взрыва, поражающие характеристики ВУВ имеют одинаковые параметры, а приведенная площадь поражения имеет форму круга, можно рассчитать ее значение по формуле

где:

R0 - радиус от центра взрыва, где вероятность поражения цели равна 1 м;
n - количество точек, разбивающих оставшуюся площадь поражения, шт.;
Ri+1 - радиус от центра взрыва до очередной точки, м;
Рi+1, Рi - граничные вероятности поражения.

Значения параметров ВУВ (избыточное давление на фронте ВУВ, импульс фазы сжатия, длительность фазы сжатия и удельная эффективная энергия) при взрыве облака ГПВС (1000 тонн пропана) для различных расстояний приведены в табл. 4.


Графики зависимостей параметров ВУВ (избыточное давление на фронте ВУВ, импульс фазы сжатия, длительность фазы сжатия и удельная эффективная энергия) при взрыве облака ГПВС (1000 тонн пропана) от расстояния до центра облака приведены на рис. 2.


Рассчитываем площади приведенных зон поражения объемного взрыва по человеку для различных типов поражения и степени укрытости (табл. 5).


С другой стороны в учебном пособии [2] приведена классическая схема сравнения мощности взрыва полусферических зарядов на поверхности земли: КВВ (тротил) и полусферического облака ГПВС (пропан), приводящие к практическому удвоению мощности за счет отражения ВУВ. Учитывая, что пропан при молекулярной массе 44 у.е., по сравнению с воздухом (не менее 29 у.е.), при естественном испарении не будет, как правило, подниматься на высоту более 150 м, как приведено в решении [2], то можно принять, что высота облака пропана не будет превышать 5-6 м, что составит площадь распространения для объема ГПВС 6,316 - 106 м3 около 1,148 - 106 м2.


Таким образом, на площади одного км2 при инициировании КВВ или случайной искре с необходимой энергией инициирования может произойти объемно-детонирующий взрыв с величиной избыточного давления на фронте ВУВ не менее 12-15 атм, что приведет к 100% поражению людей со степенью не ниже смертельной и тяжелой, что может составить по мощности взрыва не менее 50 килотонн в тротиловом эквиваленте (по ударной волне эффект сопоставим с взрывом ядерного заряда).

Выводы

1. Систематизированы результаты исследований по поражающему действию воздушных ударных волн при взрывах ГПВС. На основе анализа отечественных и зарубежных данных разработана универсальная параметрическая функция поражения человека в зависимости от удельной энергии, воздействующей на него ВУВ для всех общепринятых типов поражения.

2. Получены эмпирические зависимости основных параметров ВУВ (избыточного давления во фронте ВУВ, импульса фазы сжатия, длительности фазы сжатия и удельной энергии) от расстояния в едином приведенном виде на основе универсального метода энергетического подобия.

3. Разработана методика расчета мощности поражающего действия по человеку различной степени укрытости для всех общепринятых типов поражения. Основой методики является разработанный универсальный алгоритм на основе метода энергетического подобия, обеспечивающего широкий интервал использования методики в зависимости от характеристик заряда ТВС. Приведен пример численного расчета параметров объемного взрыва ГПВС, показывающий высокую сходимость с экспериментальными результатами и обоснованность оценки мощности взрыва СПГ, приведенной д.т.н. И.М. Абдурагимовым и к.т.н. Г.Н. Куприным. Мощность взрыва облака ГПВС более 1000 тонн сжиженного природного газа по воздушной ударной волне будет сопоставима с мощностью взрыва ядерного оружия около 10 килотонн.

Список использованных источников

  1. Абдурагимов И.М., Куприн Г.Н. Нерешенные проблемы пожаровзрывобезопасности энергоресурсов (СУГ и СПГ) как оборотная сторона успехов энергетической стратегии Российской Федерации // Пожаровзрывобезопасность. Издательство ПОЖ-НАУКА-2014. – Том 23. № 4. – С. 42–49.
  2. Аварии и катастрофы. Предупреждение и ликвидация последствий. Учебное пособие в 4-х книгах. Книга 1 / Под ред. В.А. Котляревского. – М.: Издательство АСВ, 1995. – 320 с. с ил.
  3. Физика взрыва / Под ред. Л.П. Орленко. – Изд. 3-е, переработанное. – В 2 т. Т. 1 – М.: ФИЗМАТЛИТ, 2002. – 832 с.
  4. Кудрин И.Д., Филатов А.И., Крюков А.П. Патоморфологическая классификация поражений биологических объектов и прогностическая военно-медицинская оценка поражений живой силы от воздушной ударной волны взрыва неядерных боеприпасов. – Л.: ВМедА, 1985. – 47 с.
  5. White C.S. The Scope of Blast and Shock Biology and Problem Areas in Relating Physical and Biological Parameters. – Annals of the New York Academy of Sciences, 152, Article 1, p.p. 89–102. Lovelace Foundation for Medical Education and Research, Albuquerque, New Mex., 1968.
  6. Взрывные явления. Оценка и последствия: в 2-х книгах. Пер. с англ. / Бейкер У., Кокс П. Под редакцией Зельдовича Я.Б. – М.: Мир, 1986. – 319 с.
  7. Садовский М.А. Механическое действие воздушных ударных волн взрыва по данным экспериментальных исследований / Физика взрыва. – № 1 – М.: АН СССР, 1952.
  8. Гельфанд Б.Е., Сильников М.В. Баротермическое действие взрывов. – СПб.: Астерион, 2006. – 658 с., ил.
  9. Гельфанд Б.Е., Сильников М.В. Объемные взрывы: монография. – СПб.: Астерион, 2008. – 374 с., илл.
  10. Средства поражения и боеприпасы: Учебник для втузов / А.В. Бабкин, В.А. Велданов, Е.Ф. Грязнов и др. – М.: Изд-во МГТУ им. Н.Э. Баумана, 2008. – 984 с. ил.
  11. Маньковский Г.И., Сальников А.Л., Ушаков В.С. Оценка параметров воздушной ударной волны при взрывах пиротехнических изделий и ее воздействие на человека // Современные проблемы пиротехники (Материалы III Всероссийской конференции). – Сергиев Посад: ФНПЦ "НИИ прикладной химии", 2005, с. 117–119.

Опубликовано: Каталог "Пожарная безопасность"-2017
Посещений: 7239


  Автор
Геннадий Маньковский

Геннадий Маньковский

Профессор, ФГБУ "27 НЦ" МО РФ, д.т.н.

Всего статей:  1


  Автор
Анатолий Гришкевич

Анатолий Гришкевич

Старший научный сотрудник ФГБУ "27 НЦ" МО РФ, к.т.н.

Всего статей:  1


  Автор
Владислав Филин

Владислав Филин

Доцент ФГБУ "27 НЦ" МО РФ, к.т.н.

Всего статей:  1


  Автор
Виктор Ушаков

Виктор Ушаков

Старший научный сотрудник ФГБУ "27 НЦ" МО РФ, к.т.н.

Всего статей:  1

В рубрику "Пожарная безопасность" | К списку рубрик  |  К списку авторов  |  К списку публикаций