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

Возгорание облака топливно-воздушной смеси приводит к образованию фронта пламени, распространяемого через его взрывоопасную часть. В зависимости от скорости распространения может образоваться взрывная волна. Это явление известно как взрыв парового облака.

Расчет характеристик взрывных волн и зон разрушений при возникновении аварийных ситуаций вследствие утечки взрывоопасных газов представлен при моделировании процессов на установке подготовки сырья в производстве полипропилена, на газо-фракционирующей установке (ГФУ), на установке риформинга ЛЧ-35-11/1000 и на установке каталитического крекинга Г-43-107.

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

146

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

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

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

Поскольку с ГРС производится отгрузка сжиженных газов в железнодорожных и автомобильных цистернах, возможны аварии, связанные с их заполнением и выбросом газов в атмосферу.

Газофракционирующая установка с блоком очистки головных фракций с установок предназначены для получения пропан-бутановой фракции и компонентов автобензина.

На ГФУ имеются емкости для хранения СГ. При аварийной ситуации возможен выброс пропан-бутановой фракции в атмосферу.

147

Схематично сценарии развития аварии с выбросом газа из емкостей на территории ГФУ представлены на рис. 2.13.

Установка риформинга ЛЧ-35-11 /1000, предназначенная для предварительной очистки сырья и каталитического реформирования бензиновых фракций, включает ряд блоков, представляющих наибольшую опасность возникновения пожара и взрыва (табл. 2.29, 3.30).

Таблица 2.29

Характеристики блоков установки риформинга ЛЧ-35-11/1000

Блок

Объем, м3

Объем, занимаемый водородом, м3

Масса водорода в блоке, кг

Доля жидкости, испаряющаяся при разрушении аппарата, %

Реакторы*
Р101	86,4	0,45	0,76	96
Р102	23	1,48	0,71	100
PI 03	54	3,95	1,90	100
PI 04	103,5	5,42	2,60	100
Сепаратор С101	32	16,00	53,6	0
Отпарная колонна К101	217	-	-	82
Колонна К102	161	-	-	82
Теплообменники
Т101-1...4	9,15	0,1	0,174	96

* Объем без катализатора.

148

Таблица 2.30

Основные характеристики реакторов

Тип реактора	Давление, атм	Температура, °С	Расход, м3/ч
			Бензин	ВСГ

Р-101	47	400	170000	40000
Р-102	14-16	520	170000	200000
Р-103	14-16	520	170000	200000
Р-104	14-16	520	170000	200000

Максимальный диаметр трубопроводов, проходящих по территории установки, - 350 мм. По ним перекачивается взрывоопасный водородсодержащий газ (ВСГ), содержащий до 80-90% Н2. Пожаровзрывоопасность указанных блоков установки в наибольшей степени определяется наличием в них бензина в перегретом состоянии, поскольку при возможной их разгерметизации могут наблюдаться явления, характерные для сжиженных углеводородных газов.

В качестве примера аварии с перегретым веществом можно привести аварию, произошедшую в Фликсборо (Великобритания) (табл. 2.31). В результате утечки из реакторов около 140 т циклогексана (находившегося в реакторах при температуре 155°С и давлении 0,9 МПа), его частичного испарения и загорания произошел взрыв, о силе которого можно судить по данным, приведенным ниже:

Расстояние, м	Избыточное давление взрыва, кПа
70	69
120	34
200	28
240	18
300	14
350	9

Установлено, что в случае, когда доля жидкости, испаряемая за счет перегрева, превышает 35% об., при разрушении аппарата возможно образование огневого шара по следующим сценариям:

Сценарий А: разгерметизация арматуры у одного из блоков; образование паровоздушной струи; быстрое загазование площадки установки; сгорание паровоздушного облака со взрывом, с разрушениями; последующее струйное горение бензина вследствие разгерметизации других узлов оборудования; воздействие струйных факелов на аппараты с их последующим разрывом по типу BLEVE.

149

Таблица 2.31

Характеристики разрушений при аварии в Фликсбор

Населенный пункт

Расстояние, от эпицентра, мили

Количество домов

Количество разрушенных домов

Amcotts	0,5	77	73
Flixborough	0,7	79	72
Normanby	1,8	38	35
Burton upon Stather	2,0	756	644
Gunness	2,4	333	26
Luddington	2,6	131	7
Keadby	2,75	495	5
Thealby	3,0	63	29
Alkboro	4<5	158	62
Garthorpe	4,6	111	56
Winterton	4,9	1503	11
West Halton	6,0	48	3
Whitton	6,0	48	3
Winteringhman	6,25	275	3

Сценарий В: разгерметизация арматуры у одного из блоков с мгновенным загоранием, образование струйного факела, воздействие струйного факела на аппараты с их последующим разрывом по типу DLEVE.

Сценарий С: разгерметизация арматуры у одного из компрессоров в помещении, образование паровоздушной струи, быстрое загазование помещения компрессорной, взрывное сгорание паровоздушного облака в помещении.

Сценарий D: разгерметизация арматуры на трубопроводе либо разрыв трубопровода подачи водородсодержащих газов в блоки, загорание паровоздушной струи, воздействие струйного факела на оборудование.

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

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

150

значения (10-30 м/с). В иностранной литературе огневые шары получили наименование - явление BLEVE. Характеризуются большим тепловыделением и могут вызвать пожары на больших площадях.

Одним из основных факторов, влияющих на параметры ударных волн при взрывах газовоздушных смесей, является масса газа (М), содержащегося в облаке. Расход жидкости при истечении из емкости (трубопровода) в атмосферу определяется по формуле:

G = ρ_ж[2ΔР/ρ_ж]^1/2, (2.17)

где ρ_ж - плотность жидкой фазы; ΔР - избыточное давление (ΔР = Р_нас - Р_атм); Р_нас - давление насыщенных паров.

Масса жидкости при истечении из емкости в атмосферу определяется по формуле:

M = GtS, (2.18)

где t - время истечения; S - площадь сечения отверстия.

Концентрация газа, находящегося в смеси с воздухом в пределах воспламеняемости, рассчитывается по формуле гауссовой модели, описанной ранее (в части 2.5):

(2.19)

где С - текущая концентрация вещества в облаке (зависит от времени, прошедшего с момента аварии (t) и расстояния от центра облака (r)); С₀ - исходная концентрация вещества в облаке; k = k₀(a/r₀) - коэффициент турбулентной диффузии, м²/с; а - исходный радиус облака, м; k₀ - коэффициент турбулентной диффузии для облака с радиусом r₀= 6,24 м, зависящий от состояния атмосферы (табл. 2.32) и определяемый в соответствии с устойчивостью атмосферы.

Атмосферная устойчивость характеризуется адиабатическим градиентом температуры в градусах на 100 м высоты над поверхностью (AT/AZ/100 м) и естественным вертикальным градиентом температур в зависимости от изменения давления с высотой (чем выше высота, тем ниже температура). Если вертикальный градиент температур меньше адиабатического градиента, то атмосфера

151

Таблица 2,32

Классификация и характеристика состояния атмосферы

Класс устойчивости

Характеристика устойчивости

Описание

Скорость ветра, м/с

Температурный градиент ΔТ/ΔZ, °С/100м

А	Очень неустойчивое состояние, сильно развитая конвекция	Очень солнечная летняя тихая погода	1	менее -1,9
В	Неустойчивое состояние, умеренная конвекция	Солнечно и тепло	2	-1,9...1,7
С	Слегка неустойчивое состояние, слабая конвекция	Переменная облачность в течение дня	5	-1,7...1,5
D	Безразличное нейтральное состояние	Облачный день или облачная ночь	5	-1,5...0,5
Е	Слегка устойчивое состояние, слабая инверсия	Переменная облачность в течение ночи	3	-0,5...1,5
F	Устойчивое состояние, умеренная инверсия	Ясная ночь	2	1,5... 4

152

считается неустойчивой и по коэффициенту турбулентной диффузии разделяется на шесть классов, приведенных в табл. 2.33.

Таблица 2.33

Коэффициент турбулентной диффузии ko, м2/с

Класс устойчивости атмосферы
А	В	С	D	Е	F
2,17	1,39	1,04	0,65	0,43	0,26

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

Со временем облако расплывается, образуя на границе пропан-воздух взрывоопасную смесь, которая первоначально увеличивается, затем уменьшается до полного рассеивания. Это видно по изменениям во времени расстояний от центра облака (рис. 2.15), на которых пропановоздушная смесь имеет нижний концентрационный

153

предел распространения пламени (НКПР), стехиометрический и верхний концентрационные пределы распространения пламени (ВКПР) соответственно.

Взрывоопасное облако (для наиболее устойчивого состояния атмосферы F) имеет максимальный размер через 150 с после аварийного выброса пропана. В дальнейшем его размеры (в том числе, и объем, рис. 2.15, 2.16) уменьшаются.

154

Максимальная масса пропана, участвующего во взрыве, составляет 12% от общей массы. Изменение во времени массы пропана, находящегося во взрывоопасной концентрации с воздухом, показано на рис. 2.17.

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

Таблица 2.34

Максимальные массы газа (М), участвующего во взрыве, при аварийных выбросах СУГ в атмосферу

Установка	М, т
	Z1* = l,0	Z1 = 0,5	Z1 = 0,25

Газораздаточная станция	9,6	4,8	2,4
Установка подготовки сырья при производстве полипропилена	4,8	2,4	1,2
Газофракционирующая установка	4,8	2,4	1,2

* Z₁ - часть испарившегося газа, перешедшего во взрывоопасное состояние.

155

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

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

Результаты расчетов давления ударной волны при детонационном горении ТВС, образуемой при аварийных ситуациях от технологических блоков установки Г-43-107, приведенные в табл. 2.35, показывают, что наиболее опасными в данном случае являются блоки Р-101/2 и К-201.

Таблица 2.35

Результаты расчета давления (Р) ударной волны при детонационном горении ТВС от расстояния (R) технологических блоков установки Г-43-107

Технологические блоки установки
П-101		Р-101/2		К-201		Р-201		К-304		К-306
Параметры
R, м	Р, кПа	R, м	Р, кПа	R, м	Р, кПа	R, м	Р, кПа	R, м	Р, кПа	R, м	Р, кПа

68	100	102	100	226	100	51	100	62	100	67	100
101	70	150	70	333	70	76	70	91	70	99	70
173	28	257	28	570	28	130	28	157	28	170	28
503	14	749	14	1663	14	379	14	457	14	495	14
1007	2	1497	2	3326	2	758	2	915	2	990	2

156

Поскольку невозможно заранее предсказать, произойдет ли детонация парогазовой смеси после инициирования или смесь сгорит, рассмотрено дефлаграционное (взрывное) горение парогазовой смеси.

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

• - в процессе взрывного горения фронт взрывной волны движется со сверхзвуковой скоростью, а фронт пламени распространяется с дозвуковой скоростью, зона между этими фронтами, в которой находится слой сжатого воздуха, возрастает;
• - избыточное давление и другие параметры (скорости потока воздуха и звука в потоке, плотность воздуха, скоростной напор) в слое сжатого воздуха возрастают от фронта взрывной волны к фронту пламени;
• - после окончания взрывного горения на расстоянии Rr (равном радиусу огненного шара) фронт максимального давления отрывается от фронта пламени и распространяется в сторону фронта взрывной волны со скоростью, равной сумме скоростей звука и потока газа, и на некотором расстоянии К_дог от центра взрыва догоняет фронт взрывной волны.

При дефлаграционном взрыве условно выделяются три зоны действия взрывной волны.

В 1-й зоне, расположенной в радиусе 0 < R < R_r от центра взрыва, происходит взрывное горение горючей смеси, воспринимаемое как развитие огненного шара. В этой зоне на сооружение последовательно действуют взрывная волна, движущаяся перед фронтом пламени, а затем раскаленные до 1600-2000°С продукты взрыва, избыточное давление в которых на 1-3% меньше, чем перед фронтом пламени до окончания взрывного горения, после которого на здания и сооружения, расположенные в 1-й зоне, воздействует давление разрежения, распространяющееся от конечного положения фронта пламени к центру взрыва.

Во 2-й зоне, расположенной в интервале R_r < R < R_дог, на здание последовательно воздействуют фронт взрывной волны и избыточное давление, повышающееся до максимального ΔР_r. После

157

прохождения максимального давления ΔР_m на здание воздействует понижающееся избыточное давление, а затем "давление разрежения".

На всем протяжении 1-й и 2-й зон избыточное давление на фронте взрывной волны ΔР_ф и другие параметры фронта остаются постоянными. Снижение значений указанных параметров с расстоянием, пройденным взрывной волной, начинается на границе 2-й и 3-й зон при R = К_дог и продолжается при R > R_дог.

В то же время на протяжении всей 2-й зоны происходит быстрое снижение величин максимального давления АР_т и скоростного напора q_m соответственно от значений АР_mах и q_max на границе с 1-й зоной до значений ΔР_ф и q_ф на границе с 3-й зоной.

В 3-й зоне, находящейся на расстоянии R > R_дог, взрывная волна перестает подпитываться энергией от фронта пламени и превращается в типичную сферическую ударную волну, в которой максимальное давление и другие параметры (скорость потока, плотность, скоростной напор) приходятся на фронт ударной волны. После прохождения фронта ударной волны все параметры в волне снижаются и затем наступает фаза разрежения.

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

Для определения параметров взрывной волны необходимо знать значение видимой скорости распространения пламени дефлаграционного взрыва V. Строгое определение значений видимой скорости пламени невозможно из-за влияния очень большого количества факторов и оценивается по формуле:

V = К₁К₂{220К₃ [1 - ехр(-(ЗК₄ + lg(MZ))/12,9)] + 30U_H sin(π/2 С_ср/С_стх)} м/с, (2.20)

где К₁ = 1,0-1,5 - коэффициент для учета интенсивности инициирования взрывного горения; К₂ - коэффициент для учета места воспламенения облака ГВС, равный 1 при воспламенении в центральной части облака и 0,5 при воспламенении на периферии; К₃ - коэффициент для учета плавучести (оседаемости) облаков в зависимости от соотношения плотностей и температур газопаровоздушной смеси (ГПВС) и окружающего воздуха (для нормальных условий К₃ = 0,98-1,01); К₄ - коэффициент для учета вида препятствий, вызывающих турбулизацию ГПВС перед фронтом пламени (К₄ = 1-3,5); U_H - скорость горения.

158

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

R_r = 0,782(М8_ср/С_ср)^1/3, (2.21)

где М - масса газа, участвующего во взрыве; Бср - среднее значение коэффициента расширения, определяется через полусумму величин:

ε = (Е_стех + ε_прод )/2,

где ε_стех - степень расширения продуктов сгорания стехиометрической смеси; ε_прод - степень расширения продуктов сгорания при нижнем концентрационном пределе распространения пламени; С_ср - средняя концентрация горючего вещества в смеси:

С_ср = (C_нкпр + C_стех)

где С_нкпр - нижний концентрационный предел распространения пламени; С_стех - стехиометрическая концентрация смеси.

При дефлаграционном взрыве на фронтах воздушной волны и пламени независимо от количества горючего вещества при постоянной скорости распространения пламени на одинаковых приведенных расстояниях R₀ = R / R_r в пределах 1-й и 2-й зон имеются одинаковые значения величин: избыточного давления, скоростного напора, плотности и массовой скорости потока газов.

Максимальное избыточное давление ΔР_mах определяется как

ΔР_mах = 24,22 + [(V - 140)/4,745]^1,232. (2.22)

Во 2-й зоне между фронтами взрывной волны и пламени в области 1 ≤ R₀ ≤ R_дог максимальное избыточное давление ΔР_m составит:

ΔP_m = ΔP_max/[l + B(R₀-l)^c], (2.23)

где R₀ = R/R_r.; В и С - коэффициенты, определяемые по табл. 2.36 в зависимости от величины расчетной скорости распространения пламени U_p.

Максимальное значение скоростного напора q_m определяется по аппроксимационной формуле:

q_m = 10^Y, (2.24)

где Y = А[1 - ln(R₀)/lnR_q⁰]N + C_q; A, N, C_q - величины, принимаемые по данным табл. 2.36 в зависимости от Up; Rq - приведенное расстояние, далее которого величина скоростного напора < 0,25 кПа; Rijпринимается по данным табл. 2.36 (0 < R₀ < 1; q_m = 0).

159

Таблица 2.36

Коэффициенты и показатели степени аппроксимационных формул

Коэффициент, показатель	Расчетная скорость распространения пламени Up, м/с
	100	140	170	200	240

В	0,573	0,553	0,525	0,502	0,491
С	1,055	1,069	1,077	1,098	1,165
N	1,091	1,106	1,121	1,112	1,059
F	1,630	1,804	1,791	2,005	2,190
Rq0	3,913	5,450	5,987	8,491	11,000
Cq	0,922	0,774	0,568	0,614	0,602

Расчеты максимальных уровней давления при дефлаграционном взрыве облака СУГ даны на рис. 2.18.

Результаты расчета изменения во времени массы газа, находящегося во взрывоопасной концентрации, для различных классов устойчивости атмосферы показывают (рис. 2.19), что, например, при аварийном выбросе 80 т СУГ при классе устойчивости А время рассеивания облака tup равно 200 с, а при классе

160

устойчивости F - t_дp равно 1400 с. В первом случае значительный дрейф облака практически исключается, во втором - наиболее вероятен.

Дальность дрейфа облака t_др (до расстояния, на котором еще сохраняются его взрывоопасные свойства) определяется по формуле:

L_др bu_вt_др (2.25)

где u_в - скорость ветра; b - константа (в зависимости от местности b = 0,6-1,0); t_др - время дрейфа.

Используя расчетные методики для любого класса устойчивости атмосферы, можно определить время дрейфа t_др, а по формуле (2.25) - дальность дрейфа b_др, при которых в облаке содержится определенная масса М во взрывоопасной концентрации. Расчетные данные приведены в табл. 2.37-2.39.

С учетом того, что при классе устойчивости атмосферы F время рассеивания облака имеет самое большое значение (следовательно, возможен его дрейф на большие расстояния), в табл. 2.40 приведены: время дрейфа облака t_др, его дальность b_др, масса газа в облаке М во взрывоопасной концентрации при этих условиях.

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

161

Таблица 2.37

Аварийный выброс 80 т СУГ

Класс устойчивости атмосферы	Значения tдр, с и Lдр, м при содержании массы газа во взрывоопасной концентрации М, т
	2		4		6		8
	tдр, с	Lдр, м	tдр, с	Lдр, м	tдр, с	Lдр, м	tдр, с	Lдр, м

Фаза формирования взрывоопасного облака
А	41,7	25	58,3	35	66,7	40	75,0	45
В	66,7	40	91,7	55	108,3	65	116,7	70
С	91,7	55	125,0	75	133,3	80	141,7	85
D	141,7	85	208,3	125	233,3	140	241,7	145
Е	216,7	130	308,3	185	350,0	210	366,7	220
F	375,0	225	508,3	305	575,0	345	600,0	360
Фаза рассеивания взрывоопасного облака
А	133,3	80	116,3	70	100,0	60	83,3	50
В	208,3	125	183,3	110	158,3	95	133,3	80
С	283,3	170	241,7	145	208,3	125	175,0	105
D	450,0	270	391,7	235	333,3	200	283,3	170
Е	675,0	405	591,7	355	508,3	305	433,3	260
F	1133,	680	983,3	590	841,7	505	725,0	435

Таблица 2.38

Аварийный выброс 40 т СУГ

Класс устойчивости атмосферы	tдр, с и Lдр, м при содержании массы газа во взрывоопасной концентрации М, т
	1		2		3		4
	tдр, с	Lдр, м	tдр, с	Lдр, м	tдр, с	Lдр, м	tдр, с	Lдр, м

Фаза формирования взрывоопасного облака
А	34,3	20,6	48,6	29,1	51,4	30,9	57,1	34,3
В	57,1	34,3	77,1	46,3	80,0	48	85,77	51,4
С	74,3	44,6	105,7	63,4	111,4	66,8	120,0	72,0
D	120,0	72	165,7	99,4	182,9	109,7	194,3	116,6
Е	182,9	109,7	251,4	150,9	280,0	168,0	291,4	174,9
F	297,1	178,3	417,1	250,3	468,6	281,1	491,4	294,9
Фаза рассеивания взрывоопасного облака
А	108,6	65,1	91,4	54,9	77,1	46,3	68,6	41,1
В	171,4	102,9	147,7	87,4	125,7	75,4	108,6	65,1
С	228,6	137,1	194,3	116,6	168,6	101,4	142,9	85,7
D	365,7	219,4	317,1	190,3	271,4	162,9	234,3	140,6
Е	554,3	332,6	480,0	288,0	414,3	248,6	357,1	214,3
F	920,0	552,0	794,3	476,6	691,4	414,9	591,4	354,8

162

Таблица 2.39

Аварийный выброс 20 т СУГ

Класс устойчивости атмосферы	tдр, м и Lдр, м, при содержании в облаке массы газа во взрывоопасной концентрации М, т
	1		2,4		4
	tдр, м	Lдр, м	tдр, м	Lдр, м	tдр, м	Lдр, м
	Фаза формирования облака		Фаза максимального облака		Фаза рассеивания облака

А	34,3	20,6	51,4	30,9	74,3	44,6
В	57,1	34,3	74,3	44,6	114,3	68,6
С	80,0	48,0	105,7	63,4	154,3	92,6
D	131,4	78,9	165,7	99,4	251,4	150,9
Е	200,0	120,0	240,0	144,0	382,9	229,7
F	331,4	198,9	400,0	240,0	634,3	380,6

Таблица 2.40

Параметры дрейфа облака при классе устойчивости атмосферы F

Выброс 80 т СУГ			Выброс 40 т СУГ			Выброс 20 т СУГ
tдр, м	Lдр, м	tдр, м	Lдр, м	tдр, м	Lдр, м	tдр, м	Lдр, м	tдр, м

300	150	1,12	200	120	0,56	300	180	0,76
400	300	3,44	400	240	1,76	400	240	2,40
500	375	9,60	500	300	4,80	500	300	1,76
600	450	7,76	600	360	3,84	600	360	1,14
700	600	3,92	800	480	1,80	700	420	0,60

в атмосферу) и наименее вероятным. Во втором случае (рис. 2.206) радиус зон поражения с сильной степенью разрушения составляет 200-300 м.

Далее приведены результаты расчета давления ударных волн при дефлаграционном горении ТВ С, образующихся при авариях технологических блоков установки каталитического крекинга Г-43-107. Расчет выполнен для наиболее неблагоприятного варианта при инициировании горения в центральной части облака. Рассмотрены ситуации при различных коэффициентах турбулизации (1,5-3,5), характеризующих минимальное и максимальное загромождение пространства в месте воспламенения смеси.

163

Исходные данные для расчетов скорости распространения пламени для блоков установки Г-43-107 приведены в табл. 2.41 и 2.42. Параметры и характеристики распространения пламени были рассчитаны исходя из состава ПГФ и физико-химических данных для индивидуальных веществ.

Таблица 2.41

Исходные данные для расчета скорости распространения пламени

Показатель	Технологические блоки установки
	П-101	Р-101/2	К-201	Р-201	К-304	К-306

Средняя концентрация горючего, г/м3	17,32	23,70	47,20	42,68	57,43	57,78
Средний коэффициент расширения	4,41	4,53	4,72	5,22	5,97	6,61
Максимальная скорость распространения пламени, м/с	2,22	1,82	0,38	0,35	0,37	0,46

164

Таблица 2.42

Свойства индивидуальных веществ

Вещество	Свойства
	Сстех, г/м3	СНКПР, г/м3	ζстех	ζнкпр	Uн, м/с	Сср	ζcp

Н2	24,7	3,5	6,9	2,1	2,670	14,1	4,5
СН4	63,5	31,7	7,4	5,0	0,338	47,6	6,2
С2Н4	67,6	234,4	7,8	5,0	0,740	45,5	6,4
С2Н6	72,6	36,3	7,9	5,1	0,401	54,5	6,5
C3H6	78,2	39,1	8,2	5,0	0,583	58,7	6,6
С3Н8	74,1	37,1	8,1	5,1	0,455	55,6	6,6
С4Н8	79,0	38,0	8,2	5,1	0,375	58,5	6,7
С4Н10	75,8	38,0	8,0	5,2	0,379	56,9	6,6
С5+	78,4	55,9	8,1	5,0	0,385	67,2	6,6
С6+	79,1	66,0	8,2	5,1	0,335	72,6	6,7

В табл. 2.43 для технологических блоков установки Г-43-107 приведены результаты расчета давления ударной волны в зависимости от расстояния при дефлаграционном горении ТВС при различных скоростях горения смесей. Анализ полученных результатов показывает, что авария на одном из блоков установки может привести к разрушению коммуникаций и технологического оборудования смежных блоков, что способно вызвать образование и взрыв топливно-воздушной смеси, т.е. развитие сценария по принципу "домино". Наиболее опасными в этом отношении являются реактор гидроочистки Р-101/2 и ректификационная колонна К-102.

Далее, в качестве примера, приведены результаты расчета максимального давления в волне сжатия при дефлаграционном взрыве при возможных аварийных ситуациях на установке ЛЧ-35-11/1000. Характеристики блоков установки даны в табл. 2.29.

На рис. 2.21 приведены максимальные значения давления при взрыве массы Мг = 15 т перегретого бензина. При этом скорость пламени изменялась в пределах: 103,4-158,0 м/с, что соответствует минимальной и максимальной загроможденное™ пространства в месте воспламенения смеси. Взрыв такого количества перегретого бензина (1-й тип аварии по сценарию А) возможен при холодном разрушении резервуаров К-101 или К-102. Частота подобного события составляет 1,3 · 10^-7 год^-1, поэтому оно маловероятно.

При разгерметизации межблочных трубопроводов и арматуры установки ЛЧ-35-11/1000 с загоранием паровоздушной струи бензина возможно образование струйных факелов в соответствии

165

Таблица 2.43

Результаты расчета
давления ударной волны при дефлаграционном горении ТВС
в зависимости от расстояния (R)
для технологических блоков установки Г-43-107

Технологические блоки установки
П-101		Р-101/2		К-201		Р-201		К-304		К-306
V1 = 166,3 м/с		V1 = 165,5 м/с		V1 = 162,7 м/с		V1 = 152,9 м/с		V1 = 154,7 м/с		V1= 150,8 м/с
R, м	Р, кПа	R, м	Р, кПа	R, м	Р, кПа	R, м	Р, кПа	R, м	Р, кПа	R, м	Р, кПа

30	12,5	30	16,9	30	23,7	30	6,4	30	7,6	30	7,6
50	7,7	50	10,8	50	16,0	50	3,9	50	4,6	50	4,6
75	5,1	75	7,3	75	11,2	75	2,5	75	3,0	75	3,1
100	3,8	100	5,5	100	8,5		1,9	100	2,3	100	2,3
150	2,5	150	3,6	150	5,7			150	1,5	150	1,5
200	1,9	200	2,7	200	4,3
		300	1,8	300	2,8
				500	1,7

V2= 108,1 м/с		V2= 109,4 м/с		V2 = 110,9 м/с		V2=93,4 м/с		V2=96,l м/с		V2 = 93,3 м/с
R, м	P, кПа	R, м	P, кПа	R, м	P, кПа	R, м	P, кПа	R, м	P, кПа	R, м	P, кПа

30	5,7	30	7,9	30	11,7	30	2,6	30	3,2	30	3,1
50	3,5	50	5,0	50	7,9	50	1,6	50	1,9	50	1,9
75	2,3	75	3,4	75	5,5			75	1,3
100	1,7	100	2,6	100	4,2
		150	1,7	150	2,8
				200	2,1
				300	1,4

* V₁ и V₂ - максимальные и минимальные скорости горения соответственно со сценариями А, В. Длина струйных факелов при авариях на различных блоках приведена ниже:

Блоки установки	Длина струйного факела, м
	осесимметричного	веерообразного

Реактор Р101	55	20
Реактор Р102	30	8
Реактор Р103	30	8
Колонна К101	30	8
Колонна К 102	50	16
Реактор Р104	50	16

166

Струйные факелы могут образовываться при разгерметизации запорной арматуры или разрыве трубопроводов для транспортировки водородсодержащего газа (см. сценарий D). При расчетах параметров струйных факелов при истечении водородсодержащего сырья диаметр аварийного отверстия принимался равным 1/5 максимального диаметра трубопровода установки. Результаты расчетов параметров струйного факела приведены в табл. 2.44.

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

Таблица 2.44

Параметры струйного факела

Давление ВСГ в трубопроводе, МПа

Длина факела, м

Малое основание факела, м

Большое основание факела, м

Угол расхождения факела, °

1,4	22,4	0,10	3,6	5,5
4,7	37,3	0,16	5,9	4,4

167

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

Установлено, что время нагрева, при котором наступает снижение прочности, для всех блоков составляет от 1 до 3 мин. Вследствие этого, рекомендуется на блоках Р101, Р102, Р103, Р104, К101, К102 установить систему водяного орошения с автоматическим пуском с интенсивностью подачи воды не менее 0,1 лм^-2с^-1 для поверхностей без арматуры и не менее 0,5 лм^-2с^-1 - для поверхностей с арматурой.

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

Расчет плотности теплового потока излучения и зон поражения от огневого шара. При образовании огневого шара (в результате полного разрушения холодного резервуара) его радиус R_хол (м), а также длительность существования t_хол (с) определяются по формулам:

R = 29 · М^1/3; t_хол = {

4,5 · М1/3 при М < 37т;

8,2 · M1/3 при М > 37 т;

  (2.26)

где М - масса сжиженных углеводородных газов или перегретой ЛВЖ, т.

Диаметр огневого шара dbl (м), а также длительность его существования tbl (с), при взрыве резервуара по двум вариантам разрушения, что возможно при потере прочности и разрыве резервуара в условиях пожара, определяются по зависимостям:

d_b1 = 3,44 [44,8М/μ]^1/3;   t_b1 = 0,31 [44.8М/μ]^1/3, (2.27)

где М - масса СУГ или перегретой ЛВЖ, кг; μ - молекулярная масса.

168

Плотность потока теплового излучения от огневого шара q рассчитывается по зависимости:

q = qt_fbFτ, (2.28)

где F - коэффициент облученности, рассчитываемый как

F = R²x/(R² + х²)³/², (2.29)

R - радиус огневого шара, м; х - расстояние до эпицентра, м; q_fb - средняя поверхностная плотность потока теплового излучения, принимаемая равной 225 кВт, м"2; τ - прозрачность атмосферы, определяемая по формуле:

τ = 1-0,0565 lnх. (2.30)

Образование огневого шара в случае полного разрушения резервуара в холодном состоянии является одним из наиболее опасных явлений. Это связано с высокой скоростью процесса образования шара, а также с большими радиусами теплового поражения. Вместе с тем указанное событие является маловероятным (частота возникновения - 1,3 · 10^-7год^-1), что значительно ниже верхней границы допустимого уровня воздействия опасных факторов пожара на одного человека - 10^-6год^-1 (согласно ГОСТ). Таким образом, с помощью современных методов расчета можно эффективно прогнозировать последствия пожаров и взрывов углеводородных выбросов и топливно-воздушных смесей.

Вопросы:

• Опишите возможные сценарии аварий на установках ГФУ, риформинга и каталитического крекинга.
• Опишите алгоритм расчета радиуса и времени распространения давления взрывных волн при взрывах газовоздушных смесей.
• Охарактеризуйте три зоны разрушения при взрывах технологических емкостей.
• Чем отличается по воздействию на окружающую среду и физическим характеристикам детонационный взрыв от дефлаграционного?
• Назовите особенности взрыва на открытых и закрытых технологических установках.
• Как вы понимаете эффект домино при взрывах на технологических установках?
• В чем заключаются особенности возникновения и развития зон поражения от огневого шара?

169