B = (1 + s2) / (2s), (b.10)

ПРИЛОЖЕНИЕ Д

(рекомендуемое)

 

МЕТОД РАСЧЕТА ИНТЕНСИВНОСТИ ТЕПЛОВОГО ИЗЛУЧЕНИЯ И ВРЕМЕНИ СУЩЕСТВОВАНИЯ «ОГНЕННОГО ШАРА»

 

Д.1 Расчет интенсивности теплового излучения «огненного шара» q, кВт/м², проводят по формуле

q = E_f · F_q ·t, (Д.1)

где E_f - среднеповерхностная плотность теплового излучения пламени, кВт/м²;

F_q - угловой коэффициент облученности;

t - коэффициент пропускания атмосферы.

Д.2 E_f определяют на основе имеющихся экспериментальных данных. Допускается принимать E_f равным 450 кВт/м².

Д.3 F_q рассчитывают по формуле

, (Д.2)

где Н - высота центра «огненного шара», м;

D_s - эффективный диаметр «огненного шара», м;

r - расстояние от облучаемого объекта до точки на поверхности земли непосредственно под центром «огненного шара», м.

Д.4 Эффективный диаметр «огненного шара» D_s рассчитывают по формуле

D_s =5,33 m^0,327, (Д.3)

где т - масса горючего вещества, кг.

Д.5 H определяют в ходе специальных исследований. Допускается принимать H равной D_s/2.

Д.6 Время существования «огненного шара» t_s, с, рассчитывают по формуле

t_s = 0,92 m^0,303. (Д.4)

Д.7 Коэффициент пропускания атмосферы т рассчитывают по формуле

. (Д.5)

Пример - Определить время существования «огненного шара» и интенсивность теплового излучения от него на расстоянии 500 м при разрыве сферической емкости с пропаном объемом 600 м³ в очаге пожара.

Данные для расчета

Объем сферической емкости 600 м³. Плотность жидкой фазы 530 кг/м³. Степень заполнения резервуара жидкой фазы 80 %. Расстояние от облучаемого объекта до точки на поверхности земли непосредственно под центром «огненного шара» 500 м.

Расчет

Находим массу горючего т в «огненном шаре» по формуле

т = V · r · a = 600 · 530 · 0,8 = 2,54 · 10⁵ кг,

где V - объем резервуара, м³ (V = 600 м³);

r - плотность жидкой фазы, кг/м³ (r = 530 кг/м³);

a - степень заполнения резервуара (a = 0,8).

По формуле (Д.3) определяем эффективный диаметр «огненного шара» D_s

D_s = 5,33 (2,54 · 10⁵)^0,327 = 312 м.

По формуле (Д.2), принимая H = D_s /2 = 156 м, находим угловой коэффициент облученности F_q

По формуле (Д.5) находим коэффициент пропускания атмосферы t:

.

По формуле (Д.1), принимая E_f = 450 кВт/м², находим интенсивность теплового излучения q

q = 450 · 0,037 · 0,77 = 12,9 кВт/м².

По формуле (Д.4) определяем время существования «огненного шара» t_s

t_s = 0,92 (2,54 · 10⁵)^0,303 = 40 с.

 

 

ПРИЛОЖЕНИЕ Е

(рекомендуемое)

 

МЕТОД РАСЧЕТА ПАРАМЕТРОВ ВОЛНЫ ДАВЛЕНИЯ ПРИ СГОРАНИИ ГАЗОПАРОВОЗДУШНЫХ СМЕСЕЙ В ОТКРЫТОМ ПРОСТРАНСТВЕ

 

Е.1 Исходя из рассматриваемого сценария аварии, определяют массу т, кг, горючих газов и (или) паров, вышедших в атмосферу из технологического аппарата (приложение А).

Е.2 Избыточное давление Dp, кПа, развиваемое при сгорании газопаровоздушных смесей, рассчитывают по формуле

, (Е.1)

где р₀ - атмосферное давление, кПа (допускается принимать равным 101 кПа);

r - расстояние от геометрического центра газопаровоздушного облака, м;

m_пp - приведенная масса газа или пара, кг, рассчитанная по формуле

m_пр = (Q_сг / Q₀) m_г,п Z,

где Q_сг - удельная теплота сгорания газа или пара, Дж/кг;

Z - коэффициент участия, который допускается принимать равным 0,1;

Q₀ - константа, равная 4,52 · 10⁶ Дж/кг;

m_г,п - масса горючих газов и (или) паров, поступивших в результате аварии в окружающее пространство, кг.

Е.3 Импульс волны давления i, Па · с, рассчитывают по формуле

. (Е.3)

Пример - Рассчитать избыточное давление и импульс волны давления при выходе в атмосферу пропана, хранящегося в сферической емкости объемом 600 м³, на расстоянии 500 м от нее.

Данные для расчета

Объем емкости 600 м³. Температура 20 °С. Плотность сжиженного пропана 530 кг/м³. Степень заполнения емкости 80 % (по объему). Удельная теплота сгорания пропана 4,6 · 10⁷ Дж/кг. Принимается, что в течение времени, необходимого для выхода сжиженного газа из емкости, весь пропан испаряется.

Расчет

Находим приведенную массу m_пр по формуле (Е.2):

m_пр = 4,6 · 10⁷ / 4,52 · 10⁶ ·(0,8· 530 · 600) · 0,1 = 2,59 · 10⁵ кг.

Находим избыточное давление Dp по формуле (Е.1)

Dp = 101 [0,8 (2,59 · 10⁵)^0,33 / 500 + 3 (2,59 · 10⁵)^0,66 / 500² + 5 (2,59 · 10⁵) / 500³] = 16,2 кПа.

Находим импульс волны давления i по формуле (Е.3):

i = 123 (2,59 · 10⁵)^0,66 / 500 = 1000 Па · с.

 

 

ПРИЛОЖЕНИЕ Ж

(рекомендуемое)

 

МЕТОД РАСЧЕТА ПАРАМЕТРОВ ВОЛНЫ ДАВЛЕНИЯ ПРИ ВЗРЫВЕ РЕЗЕРВУАРА С ПЕРЕГРЕТОЙ ЖИДКОСТЬЮ ИЛИ СЖИЖЕННЫМ ГАЗОМ ПРИ ВОЗДЕЙСТВИИ НА НЕГО ОЧАГА ПОЖАРА

 

Ж.1 При попадании замкнутого резервуара со сжиженным газом или жидкостью в очаг пожара может происходить нагрев содержимого резервуара до температуры, существенно превышающей нормальную температуру кипения, с соответствующим повышением давления. За счет нагрева несмоченных стенок сосуда уменьшается предел прочности их материала, в результате чего при определенных условиях оказывается возможным разрыв резервуара с возникновением волн давления и образованием «огненного шара». Расчет параметров «огненного шара» изложен в приложении Д. Порядок расчета параметров волн давления изложен ниже. Разрыв резервуара в очаге пожара с образованием волн давления получил название BLEVE (Boiling Liquid Expanding Vapour Explosion - взрыв расширяющихся паров вскипающей жидкости).

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

Ж.2.1 Рассчитывают d по формуле

d = С_р (T - T_кип) / L, (Ж.1)

где С_р - удельная теплоемкость жидкой фазы, Дж/кг;

Т - температура жидкой фазы, соответствующая температуре насыщенного пара при давлении срабатывания предохранительного клапана, К;

T_кип - температура кипения вещества при нормальном давлении, К;

L - удельная теплота испарения при нормальной температуре кипения Т_кип, Дж/кг.

Ж.2.2 Если d < 0,35, BLEVE не происходит. При d ³ 0,35 вероятность возникновения данного явления велика.

Ж.3 Параметрами волны давления, образующейся при BLEVE, являются избыточное давление в положительной фазе волны Dp и безразмерный импульс положительной фазы волны i.

Dp, кПа, и i, Па·с, рассчитывают по формулам:

, (Ж.2)

, (Ж.3)

где p₀ - атмосферное давление, кПа (допускается принимать равным 101 кПа);

r - расстояние до разрушающегося технологического оборудования, м;

m_пр - приведенная масса, кг, рассчитанная по формуле

m_пр = E_иэ / Q₀, (Ж.4)

где E_иэ - энергия, выделяющаяся при изэнтропическом расширении среды, находящейся в резервуаре, Дж;

Q₀ - константа, равная 4,52 · 10⁶ Дж/кг.

Ж.4 E_иэ, Дж, рассчитывают по формуле

E_иэ = С_эфф m (Т - Т_кип). (Ж.5)

где m - масса вещества в резервуаре, кг;

С_эфф - константа, равная 500 Дж/(кг·К);

Т - температура вещества в резервуаре в момент его взрыва, К;

Т_кип - температура кипения вещества при атмосферном давлении, К.

При наличии в резервуаре предохранительного клапана Т, К, допускается рассчитывать по формуле

, (Ж.6)

где А, В, С_а - константы Антуана вещества;

р_к - давление срабатывания предохранительного клапана, кПа.

Константа А должна соответствовать давлению, выраженному в килопаскалях.

Пример - Расчет параметров ударной волны при BLEVE

Данные для расчета

Рассчитать параметры положительной фазы волны давления на расстоянии 750 м от эпицентра аварии, связанной с развитием BLEVE на железнодорожной цистерне вместимостью 50 м³ с 10 т жидкого пропана. Цистерна имеет предохранительный клапан на давление срабатывания 2,0 МПа.

Расчет

Энергию, выделившуюся при расширении среды в резервуаре, рассчитывают по формуле (Ж.5).

E_иэ = С_эфф m (Т - Т_кип),

где m = 4 · 10⁴ кг - масса пропана в цистерне;

С_эфф - константа, равная 500 Дж/кг·К);

Т_кип = -43 + 273 = 230 К - температура кипения пропана при постоянном давлении.

Т, К, находим по формуле (Ж.6)

где р_к = 2,000 кПа, А = 5,949, В = 812,648, С_а = 247,55.

Отсюда

_━━━

Получим E_иэ

E_иэ = 4 · 10⁴ (332 - 230) 500 = 2,06 · 10⁹ Дж.

Находим приведенную массу m_пр, кг, по формуле (Ж.4)

m_пр = 2,06 · 10⁹ / (4,52 · 10⁶) = 456 кг.

Вычислим Dp и i по формулам (Ж.2) и (Ж.3)

Dр = 101 (0,8 · 456^0,33 / 750 + 3 · 456^0,66 / 750² + 5 · 456³ / 750) = 0,86 кПа,

i = 123 · 456^0,66 / 750 = 9,7 Па · с.

 

 

ПРИЛОЖЕНИЕ И

(рекомендуемое)

 

МЕТОД РАСЧЕТА ПАРАМЕТРОВ ИСПАРЕНИЯ ГОРЮЧИХ НЕНАГРЕТЫХ ЖИДКОСТЕЙ И СЖИЖЕННЫХ УГЛЕВОДОРОДНЫХ ГАЗОВ

 

И.1 Интенсивность испарения W, кг/(с·м²), определяют по справочным и экспериментальным данным. Для ненагретых выше температуры окружающей среды ЛВЖ, при отсутствии данных допускается рассчитывать W по формуле¹⁾

, (И.1)

______________

¹⁾ Формула применима при температуре подстилающей поверхности от минус 50 до плюс 40 °С.

 

где h - коэффициент, принимаемый по таблице И.1 в зависимости от скорости и температуры воздушного потока над поверхностью испарения;

М - молярная масса, г/моль;

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

Таблица И.1

 

Скорость воздушного потока в помещении, м/с	Значение коэффициента hпри температуре t, °С, воздуха в помещении
	10	15	20	30	35
0,0	1,0	1,0	1,0	1,0	1,0
0,1	3,0	2,6	2,4	1,8	1,6
0,2	4,6	3,8	3,5	2,4	2,3
0,5	6,6	5,7	5,4	3,6	3,2
1,0	10,0	8,7	7,7	5,6	4,6

 

И.2 Для сжиженных углеводородных газов (СУГ) при отсутствии данных допускается рассчитывать удельную массу паров испарившегося СУГ m_СУГ, кг/м², по формуле¹)

, (И.2)

___________

¹⁾ Формула применима при температуре подстилающей поверхности от минус 50 до плюс 40 °С.

 

где М - молярная масса СУГ, кг/моль;

L_исп - мольная теплота испарения СУГ при начальной температуре СУГ Т_ж, Дж/моль;

Т₀ - начальная температура материала, на поверхность которого разливается СУГ, соответствующая расчетной температуре t_p, К;

Т_ж - начальная температура СУГ, К;

l_тв - коэффициент теплопроводности материала, на поверхность которого разливается СУГ, Вт/(м · К);

а - эффективный коэффициент температуропроводности материала, на поверхность которого разливается СУГ, равный 8,4·10^-8 м²/с;

t - текущее время, с, принимаемое равным времени полного испарения СУГ, но не более 3600 с;

число Рейнольдса (v - скорость воздушного потока, м/с; d - характерный размер пролива СУГ, м;

u_в - кинематическая вязкость воздуха при расчетной температуре t_р, м²/с);

l_в - коэффициент теплопроводности воздуха при расчетной температуре t_р, Вт/(м · К).

Примеры - Расчет параметров испарения горючих ненагретых жидкостей и сжиженных углеводородных газов

1 Определить массу паров ацетона, поступающих в объем помещения в результате аварийной разгерметизации аппарата.

Данные для расчета

В помещении с площадью пола 50 м² установлен аппарат с ацетоном максимальным объемом V_aп = 3 м³. Ацетон поступает в аппарат самотеком по трубопроводу диаметром d = 0,05 м с расходом q, равным 2 · 10^-3 м³/с. Длина участка напорного трубопровода от емкости до ручной задвижки L₁ = 2 м. Длина участка отводящего трубопровода диаметром d = 0,05 м от емкости до ручной задвижки L₂ равна 1 м. Скорость воздушного потока u в помещении при работающей общеобменной вентиляции равна 0,2 м/с. Температура воздуха в помещении t_р = 20 °С. Плотность r ацетона при данной температуре равна 792 кг/м³. Давление насыщенных паров ацетона р_a при t_р равно 24,54 кПа.

Расчет

Объем ацетона, вышедшего из напорного трубопровода, V_н.т составляет

м³,

где t - расчетное время отключения трубопровода, равное 300 с (при ручном отключении).

Объем ацетона, вышедшего из отводящего трубопровода V_от составляет

м³.

Объем ацетона, поступившего в помещение

V_a = V_ап + V_н.т + V_от = 3 + 6,04 · 10^-1 + 1,96 · 10^-3 = 6,600 м³.

Исходя из того, что 1 л ацетона разливается на 1 м² площади пола, расчетная площадь испарения S_р = 3600 м² ацетона превысит площадь пола помещения. Следовательно, за площадь испарения ацетона принимается площадь пола помещения, равная 50 м².

Интенсивность испарения равна:

кг/(с · м²).

Масса паров ацетона, образующихся при аварийной разгерметизации аппарата т, кг, будет равна

т = 0,655 · 10^-3 · 50 · 3600 = 117,9 кг.

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

Данные для расчета

Изотермический резервуар сжиженного этилена объемом V_и.р.э = 10000 м³ установлен в бетонном обваловании свободной площадью S_об = 5184 м² и высотой отбортовки Н_об ⁼ 2,2 м. Степень заполнения резервуара a = 0,95.

Ввод трубопровода подачи сжиженного этилена в резервуар выполнен сверху, а вывод отводящего трубопровода снизу.

Диаметр отводящего трубопровода d_тp = 0,25 м. Длина участка трубопровода от резервуара до автоматической задвижки, вероятность отказа которой превышает 10^-6 в год и не обеспечено резервирование ее элементов, L = 1 м. Максимальный расход сжиженного этилена в режиме выдачи G_ж.э = 3,1944 кг/с. Плотность сжиженного этилена r_ж.э при температуре эксплуатации Т_эк= 169,5 К равна 568 кг/м³. Плотность газообразного этилена r_г.э при Т_эк равна 2,0204 кг/м³. Молярная масса сжиженного этилена М_ж.э = 28 · 10^-3 кг/моль. Мольная теплота испарения сжиженного этилена L_иcn при Т_эк равна 1,344 · 10⁴ Дж/моль. Температура бетона равна максимально возможной температуре воздуха в соответствующей климатической зоне T_б = 309 К. Коэффициент теплопроводности бетона l_б=1,5Вт/(м·К). Коэффициент температуропроводности бетона а = 8,4 · 10^-8 м²/с. Минимальная скорость воздушного потока u_min = 0 м/с, а максимальная для данной климатической зоны u_max = 5 м/с. Кинематическая вязкость воздуха v_в при расчетной температуре воздуха для данной климатической зоны t_р = 36 °С равна 1,64 · 10^-5 м²/с. Коэффициент теплопроводности воздуха l_в при t_р равен 2,74 · 10^-2 Вт/(м·К).

Расчет

При разрушении изотермического резервуара объем сжиженного этилена составит

м³.

Свободный объем обвалования V_об = 5184 · 2,2 = 11404,8 м³.

Ввиду того, что V_ж.э < V_об примем за площадь испарения S_исп свободную площадь обвалования S_об, равную 5184 м².

Тогда массу испарившегося этилена m_и.э с площади пролива при скорости воздушного потока u = 5 м/с рассчитывают по формуле (И.2)

Масса m_и.э при u = 0 м/с составит 528039 кг.