Б.2 Метод расчета размеров зон, ограниченных НКПР газов и паров, при аварийном поступлении горючих газов и паров ненагретых легковоспламеняющихся жидкостей в помещение
Нижеприведенные расчетные формулы применяют для случая 100 m / (rг,пVсв)< 0,5 СНКПР [СНКПР - нижний концентрационный предел распространения пламени горючего газа или пара, % (об.)] и помещений в форме прямоугольного параллелепипеда с отношением длины к ширине не более 5.
Б.2.1 Расстояния XНКПР , YНКПР и ZНКПР рассчитывают по формулам
, (Б.5)
, (Б.6)
, (Б.7)
где К1 - коэффициент, принимаемый равным 1,1314 для горючих газов и 1,1958 для легковоспламеняющихся жидкостей;
К2 - коэффициент, равный 1 для горючих газов;
для легковоспламеняющихся жидкостей;
К - коэффициент, принимаемый равным 0,0253 для горючих газов при отсутствии подвижности воздушной среды; 0,02828 для горючих газов при подвижности воздушной среды; 0,04714 для легковоспламеняющихся жидкостей при отсутствии подвижности воздушной среды и 0,3536 для легковоспламеняющихся жидкостей при подвижности воздушной среды;
h — высота помещения, м.
d, l, b и C0 приведены в А.2.3.
При отрицательных значениях логарифмов расстояния XНКПР , YНКПР и ZНКПР принимают равными 0.
Б.2.2 Радиус Rб и высоту Zб, м, зоны, ограниченной НКПР газов и паров, вычисляют исходя из значений XНКПР , YНКПР и ZНКПР для заданного уровня значимости Q.
При этом Rб > XНКПР , Rб > YНКПР и Zб > h + Rб для ГГ и Zб > ZНКПР для ЛВЖ (h — высота источника поступления газа от пола помещения для ГГ тяжелее воздуха и от потолка помещения для ГГ легче воздуха, м).
Для ГГ геометрически зона, ограниченная НКПР газов, будет представлять цилиндр с основанием радиусом Rб и высотой hб = 2 Rб при Rб £ h, hб = h + Rб при Rб > h, внутри которого расположен источник возможного выделения ГГ. Для ЛВЖ геометрически зона, ограниченная НКПР паров, будет представлять цилиндр с основанием радиусом Rб и высотой Zб = ZНКПР высоте источника паров ЛВЖ h < ZНКПР и Zб = h + ZНКПР при h ³ ZНКПР. За начало отсчета принимают внешние габаритные размеры аппаратов, установок, трубопроводов и т. п.
Б.2.3 Во всех случаях значения расстояний XНКПР , YНКПР и ZНКПР должны быть не менее 0,3 м для ГГ и ЛВЖ.
Примеры
1. Определить размеры зоны, ограниченной НКПР паров, образующейся при аварийной разгерметизации аппарата с ацетоном, при работающей и неработающей общеобменной вентиляции.
Данные для расчета
В центре помещения размером 40 х 40 м и высотой hп = 3 м установлен аппарат с ацетоном. Аппарат представляет собой цилиндр с основанием диаметром da = 0,5 м и высотой ha = 1 м, в котором содержится 25 кг ацетона. Расчетная температура в помещении tp = 30 °С. Плотность паров ацетона rа при tр равна 2,33 кг/м3. Давление насыщенных паров ацетона pн при tр равно 37,73 кПа. Нижний концентрационный предел распространения пламени СНКПР = 2,7 % (об.). В результате разгерметизации аппарата в помещение поступит 25 кг паров ацетона за время испарения Т = 208 с. При работающей общеобменной вентиляции подвижность воздушной среды в помещении u = 0,1 м/с.
Расчет
Допустимые значения отклонений концентраций d при уровне значимости Q = 0,05 будут равны: 1,27 — при работающей вентиляции; 1,25 — при неработающей вентиляции (u = 0). Предэкспоненциальный множитель С0 будет равен:
при работающей вентиляции
% (об.),
Сн = 100рн /р0 = 100 · 37,73/101 = 37,36 % (об.),
Vсв = 0,8 Vп = 0,8 · 40 · 40 · 3 = 3840 м3;
при неработающей вентиляции
% (об.).
Расстояния XНКПР , YНКПР и ZНКПР составят:
при работающей вентиляции
м,
м,
м;
при неработающей вентиляции
м,
м,
м.
Таким образом, для ацетона геометрически зона, ограниченная НКПР паров, будет представлять собой цилиндр с основанием радиусом Rб и высотой Zб = hа+ZНКПР, так как hа > ZHKHP, при работающей вентиляции
Zб = 1 + 0,2 = 1,2 м, Rб = 9,01 м;
при неработающей вентиляции
Zб = 1 + 0,03 = 1,03 м, Rб = 10,56 м.
За начало отсчета принимают внешние габаритные размеры аппарата.
2. Определить размеры зоны, ограниченной НКПР газов, образующейся при аварийной разгерметизации газового баллона с метаном, при работающей и неработающей вентиляции.
Данные для расчета
На полу помещения размером 13 х 13 м и высотой Hп = 3 м находится баллон с 0,28 кг метана. Газовый баллон имеет высоту hб = 1,5 м. Расчетная температура в помещении tр = 30 °С. Плотность метана rм при tр равна 0,645 кг/м3. Нижний концентрационный предел распространения пламени метана СНКПР = 5,28 % (об.). При работающей общеобменной вентиляции подвижность воздушной среды в помещении u = 0,1 м/с.
Расчет
Допустимые отклонения концентраций при уровне значимости Q = 0,05 будут равны: 1,37 при работаюшей вентиляции; 1,38 при неработающей вентиляции (u = 0).
Предэкспоненциальный множитель С0 будет равен:
при работающей вентиляции
% (об.);
при неработающей вентиляции
% (об.);
Расстояния XНКПР , YНКПР и ZНКПР составят:
при работающей вентиляции
,
,
.
следовательно XНКПР , YНКПР и ZНКПР = 0;
при неработающей вентиляции
м,
м,
м.
Таким образом, для метана при неработающей вентиляции геометрически зона, ограниченная НКПР газов, будет представлять собой цилиндр с основанием радиусом Rб = 3,34 м и высотой hб = h + Rб = 3 + 3,34 = 6,34 м. Ввиду того, что hб расчетное больше высоты помещения hп = 3 м, за высоту зоны, ограниченной НКПР газов, принимаем высоту помещения hб = 3 м.
ПРИЛОЖЕНИЕ В (рекомендуемое) МЕТОД РАСЧЕТА ИНТЕНСИВНОСТИ ТЕПЛОВОГО ИЗЛУЧЕНИЯ ПРИ ПОЖАРАХ ПРОЛИВОВ ЛВЖ И ГЖ
В. 1 Интенсивность теплового излучения q, кВт/м2, рассчитывают по формуле
q = Ef · Fq · t, (B.1)
где Ef — среднеповерхностная плотность теплового излучения пламени, кВт/м2;
Fq — угловой коэффициент облученности;
t — коэффициент пропускания атмосферы.
В.2 Ef принимают на основе имеющихся экспериментальных данных. Для некоторых жидких углеводородных топлив указанные данные приведены в таблице В. 1.
Таблица B.1— Среднеповерхностная плотность теплового излучения пламени в зависимости от диаметра очага и удельная массовая скорость выгорания для некоторых жидких углеводородных топлив
Топливо
| Ef, кВт/м2, при d, м
| т, кг/(м2 · с)
|
|
|
|
|
|
|
| СПГ (метан)
|
|
|
|
|
| 0,08
| СУГ (пропан-бутан)
|
|
|
|
|
| 0,1
| Бензин
|
|
|
|
|
| 0,06
| Дизельное топливо
|
|
|
|
|
| 0,04
| Нефть
|
|
|
|
|
| 0,04
| Примечание— Для диаметров очага менее 10 м или более 50 м следует принимать Ef такой же, как и для очагов диаметром 10 м и 50 м соответственно
|
При отсутствии данных допускается Ef принимать равной 100 кВт/м2 для СУГ, 40 кВт/м2 для нефтепродуктов.
8.3 Рассчитывают эффективный диаметр пролива d, м, по формуле
, (В.2)
где S — площадь пролива, м2.
8.4 Рассчитывают высоту пламени Н, м, по формуле
, (В.3)
где т — удельная массовая скорость выгорания топлива, кг/(м · с);
rв — плотность окружающего воздуха, кг/м3;
g— ускорение свободного падения, равное 9,81 м/с2.
8.5 Определяют угловой коэффициент облученности Fq по формуле
, (В.4)
где ,(В.5)
где А = (h2 + + 1) / 2S1 , (В.6)
Sl = 2r/d (r— расстояние от геометрического центра пролива до облучаемого объекта), (В. 7)
h = 2H/d; (B.8)
, (В.9)
B = ( 1+S2 ) / ( 2S ), (B.10)
B.6 Определяют коэффициент пропускания атмосферы t по формуле
t = exp[ -7,0 · 10 -4 ( r - 0,5 d)] (B.11)
Пример — Расчет теплового излучения от пожара пролива бензина площадью 300 м2 на расстоянии 40 м от центра пролива.
Расчет
Определяем эффективный диаметр пролива d по формуле (В. 2)
м.
Находим высоту пламени по формуле (В.3), принимая
т = 0,06 кг / (м2 · с), g = 9,81 м/с2 и rв = 1,2 кг/м3:
Находим угловой коэффициент облученности Fq по формулам (В.4) — (В. 10), принимая r = 40 м:
h = 2 · 26,5 / 19,5 = 2,72,
S1 =2 · 40 / 19,5= 4,10,
А = (2,722 + 4,102 + 1) / (2 · 4,1) = 3,08,
B = (1 + 4,12) / (2 · 4,1) =2,17,
Определяем коэффициент пропускания атмосферы т по формуле (В. 11)
t = exp [ - 7,0 · 10 -4 (40 - 0,5 · 19,5 )] = 0,979.
Находим интенсивность теплового излучения q по формуле (В.1), принимая Еf= 47 кВт/м2 в соответствии с таблицей В. 1:
q = 47 · 0,0324 · 0,979 = 1,5 кВт/м2.
ПРИЛОЖЕНИЕ Г (рекомендуемое) МЕТОД РАСЧЕТА РАЗМЕРОВ ЗОН РАСПРОСТРАНЕНИЯ ОБЛАКА ГОРЮЧИХ ГАЗОВ И ПАРОВ ПРИ АВАРИИ
Г.1 Сущность метода
В настоящем приложении установлен порядок расчета изменения во времени концентрации газа в облаке при мгновенном выбросе и непрерывном истечении сжиженного углеводородного газа (СУГ), плотность которого больше плотности воздуха.
Г.1.1 Мгновенный выброс СУГ
Г. 1.1.1 Мгновенный выброс СУГ может происходить при повреждении резервуара или иного аппарата, в котором СУГ находится под давлением.
За счет внутренней энергии СУГ его массовая доля d мгновенно испаряется, образуя с капельками жидкости облако аэрозоля. За счет больших скоростей вихревых потоков происходит быстрое вовлечение в облако воздуха и быстрое испарение оставшейся части СУГ.
Массу воздуха Ма0, кг, мгновенно вовлекающуюся в облако для такого испарения, рассчитывают по формуле
Ма0 = ( l - d) MgLg / ( Cp.a ( Ta - Tg ) + XwLw), (Г.1)
где Мg — масса выброшенного СУГ, кг;
Ср.a — удельная теплоемкость воздуха, Дж/(кг·К);
Lg — удельная теплота парообразования СУГ, Дж/кг;
Ta — температура окружающего воздуха, К;
Тg — температура кипения СУГ при атмосферном давлении, К;
Хw — массовая доля водяных паров в воздухе;
Lw — удельная теплота парообразования воды, Дж/кг.
d определяют из соотношения
d = 1 - ехр (- Сp.g (Ta - Tg) / Lg), (Г.2)
где Cp.g — удельная теплоемкость СУГ, Дж/(кг · К).
Г. 1.1.2 Принимают, что образовавшееся облако дрейфует по ветру со скоростью vd = 0,6nв (nв — скорость ветра) и имеет в начальный момент форму цилиндра, высота которого равна его радиусу. С течением времени высота облака уменьшается, а радиус растет.
Изменение во времени радиуса, высоты облака и концентрации газа в нем в этой фазе (называемой фазой падения) определяется путем решения методом Рунге-Кутта (реализованным в виде стандартной программы на ЭВМ) системы обыкновенных дифференциальных уравнений:
dMa / dt = ra p r2 a2 a3 nв Ri-1 + 2 ra а1 (dr / dt) p r h,
dT / dt =((dMa / dt) Cp.a ( Ta - T ) + p r2 ( Tgr - T )1,333 ) / ( Ma Cp.a + Mg Cp.g ) , (Г.З)
dr / dt = a4 (gh (rg.a - ra) / rg.a) 0,5,
где Ma — масса воздуха в облаке, кг;
ra — плотность воздуха, кг/м3;
r — радиус облака, м;
а1, a2, a3, a4 — коэффициенты ( а1 = 0,7, а2 = 0,5, a4 = 1,07, a3 = 0,3 для классов устойчивости А—В (классы устойчивости даны по Паскуиллу, таблица Г.1); 0,24 — для С—В; 0,16 — для E—F);
Ri — число Ричардсона, определяемое из соотношения
Ri = (5,88 h0,48 g / ( a32 nв2 )) ( rg.a - ra ) / ra ;
h — высота облака, м;
Т— температура облака. К;
Тgr — температура земной поверхности. К;
rg.a — плотность паровоздушного облака, кг/м3.
Таблица Г.1— Классы устойчивости атмосферы по Паскуиллу
Класс по Паскуиллу
| Типичная скорость ветра, м/с
| Описание погоды
| Вертикальный градиент температуры, К/м
| А
|
| Безоблачно
| >>>0,01
| В
|
| Солнечно и тепло
| >>0,01
| С
|
| Переменная облачность в течение дня
| >0,01
| D
|
| Облачный день или облачная ночь
| »0,01
| Е
|
| Переменная облачность в течение ночи
| <0,01
| F
|
| Ясная ночь
| Инверсия (отрицательный градиент)
|
Решением системы вышеуказанных уравнений являются зависимости
Ma = Ma(t), Т= Т(t), r= r(t).
Для решения системы уравнений необходимы дополнительные соотношения
rg.a = (Ma + Mg ) / (Ma / ra + Mg / rg ) ( Ta / T ). (Г.4)
В качестве критерия окончания фазы падения принимают выполнение условия:
(rg.a - ra ) / rg.a < 10-3.(Г.5)
Зависимость h = h(t) находим из соотношения
h(t)=(Ma / ra + Mg /rg ) (T /Ta )( 1/(p r(t)2). (Г.6)
Г. 1.1.3 Когда плотность паровоздушного облака незначительно отличается от плотности воздуха (т. е. после окончания фазы падения), его движение определяется как фаза пассивной дисперсии и описывается процессами турбулентной диффузии.
Концентрацию газа в точке с координатами (х, у, z) в фазе пассивной дисперсии определяют из формулы
(Г.7)
где sy, sz — среднеквадратичные отклонения, зависящие от величины xc - x0 ;
хc — координата центра облака в направлении ветра, м
x0 — координата точки окончания фазы падения, м;
sy (xc - x0); sz (xc - x0) зависят от класса устойчивости по Паскуиллу.
При xc = x0 принимается sy0 = r / 2,14, sz0 = h / 2,14;
при xc > x0
Г.1.2 Непрерывное истечение СУГ
Для описания непрерывного истечения СУГ из резервуаров или иных аппаратов предполагается, что результирующая концентрация газа в паровоздушном облаке является суммой концентраций от отдельных элементарных газовых объемов и рассчитывается по формуле
, (Г.8)
где Q = т· t j,— масса СУГ в j-м элементарном объеме, кг;
т — массовая скорость истечения СУГ, кг/с;
xj— координата центра j-го элементарного объема, м;
— среднеквадратичные отклонения распределения концентраций в j-м элементарном объеме, м.
- определяют аналогично в Г. 1.1.3.
Пример — Расчет динамики паровоздушного облака в открытом пространстве
Для расчета динамики паровоздушного облака (движения в пространстве границы облака, определяемой НКПВ) принимается, что в некоторый момент времени t0 начинается истечение пропана с массовой скоростью 1,3 кг/с, скорость ветра составляет 1 м/с, градиент температуры составляет 0,667 К/м.
Процедура расчета, реализованная на ПЭВМ, представлена на блок-схеме (рисунок Г.1). Результаты расчета границы облака для двух значений времени t0 + 10 с и t0 + 300 с представлены на рисунке Г.2.
Рисунок Г. 1 — Алгоритм расчета параметров паровоздушного облака
t0 — время начала истечения
Рисунок Г. 2 — Границы паровоздушного облака по НКПВ на различные моменты времени от начала истечения
ПРИЛОЖЕНИЕ Д (рекомендуемое) МЕТОД РАСЧЕТА ИНТЕНСИВНОСТИ ТЕПЛОВОГО ИЗЛУЧЕНИЯ И ВРЕМЕНИ СУЩЕСТВОВАНИЯ «ОГНЕННОГО ШАРА»
Д. 1 Расчет интенсивности теплового излучения «огненного шара» q, кВт/м2, проводят по формуле
q = Ef · Fq · t, (Д.1)
где Ef — среднеповерхностная плотность теплового излучения пламени, кВт/м2;
— угловой коэффициент облученности;
t коэффициент пропускания атмосферы.
Д.2 Ef определяют на основе имеющихся экспериментальных данных. Допускается принимать Ef равным 450 кВт/м2.
Д.3 Fq рассчитывают по формуле
, (Д.2)
где Н— высота центра «огненного шара», м;
Ds — эффективный диаметр «огненного шара», м;
r — расстояние от облучаемого объекта до точки на поверхности земли непосредственно под центром «огненного шара», м.
Д.4 Эффективный диаметр «огненного шара» Ds рассчитывают по формуле
Ds =5,33 m 0,327, (Д.3)
где т — масса горючего вещества, кг.
Д.5 H определяют в ходе специальных исследований. Допускается принимать H равной Ds/2.
Д.6 Время существования «огненного шара» ts, с, рассчитывают по формуле
ts = 0,92 m 0,303. (Д.4)
Д.7 Коэффициент пропускания атмосферы т рассчитывают по формуле
t = ехр [-7,0 · 10-4 ( - Ds / 2)]. (Д.5)
Пример — Определить время существования «огненного шара» и интенсивность теплового излучения от него на расстоянии 500 м при разрыве сферической емкости с пропаном объемом 600 м3 в очаге пожара.
Данные для расчета
Объем сферической емкости 600 м3. Плотность жидкой фазы 530 кг/м3. Степень заполнения резервуара жидкой фазы 80 %. Расстояние от облучаемого объекта до точки на поверхности земли непосредственно под центром «огненного шара» 500 м.
Расчет
Находим массу горючего т в «огненном шаре» по формуле
т = Vr a = 600 · 530 · 0,8 = 2,54 · 105 кг,
где V— объем резервуара, м3 (V = 600 м3);
r — плотность жидкой фазы, кг/м3 (r = 530 кг/м3);
a — степень заполнения резервуара (a = 0,8).
По формуле (Д.3) определяем эффективный диаметр «огненного шара» Ds
Ds = 5,33 (2,54 · 105)0,327 = 312 м.
По формуле (Д.2), принимая H = Ds /2 = 156 м, находим угловой коэффициент облученности Fq
По формуле (Д.5) находим коэффициент пропускания атмосферы t:
t = ехр [-7,0 · 10-4 ( - 312/2)] = 0,77.
По формуле (Д.1), принимая Ef = 450 кВт/м2, находим интенсивность теплового излучения q
q = 450 · 0,037 · 0,77 = 12,9 кВт/м2.
По формуле (Д.4) определяем время существования «огненного шара» ts
ts = 0,92 (2,54 · 105)0,303 = 40 с.
ПРИЛОЖЕНИЕ Е (рекомендуемое) МЕТОД РАСЧЕТА ПАРАМЕТРОВ ВОЛНЫ ДАВЛЕНИЯ ПРИ СГОРАНИИ ГАЗОПАРОВОЗДУШНЫХ СМЕСЕЙ В ОТКРЫТОМ ПРОСТРАНСТВЕ
Е. 1 Исходя из рассматриваемого сценария аварии, определяют массу т, кг, горючих газов и (или) паров, вышедших в атмосферу из технологического аппарата (приложение А).
Е.2 Избыточное давление Dp, кПа, развиваемое при сгорании газопаровоздушных смесей, рассчитывают по формуле
, (Е. 1)
где р0 — атмосферное давление, кПа (допускается принимать равным 101 кПа);
r — расстояние от геометрического центра газопаровоздушного облака, м;
mпp — приведенная масса газа или пара, кг, рассчитанная по формуле
mпр = (Qсг / Q0)mг,п Z,
где Qсг — удельная теплота сгорания газа или пара, Дж/кг;
Z— коэффициент участия, который допускается принимать равным 0,1;
Q0— константа, равная 4,52 · 106 Дж/кг;
mг,п — масса горючих газов и (или) паров, поступивших в результате аварии в окружающее пространство, кг.
Е.3 Импульс волны давления i, Па · с, рассчитывают по формуле
. (Е.3)
Пример — Рассчитать избыточное давление и импульс волны давления при выходе в атмосферу пропана, хранящегося в сферической емкости объемом 600 м3, на расстоянии 500 м от нее.
Данные для расчета
Объем емкости 600 м3. Температура 20 °С. Плотность сжиженного пропана 530 кг/м3. Степень заполнения емкости 80 % (по объему). Удельная теплота сгорания пропана 4,6 · 107 Дж/кг. Принимается, что в течение времени, необходимого для выхода сжиженного газа из емкости, весь пропан испаряется.
Расчет
Находим приведенную массу mпр по формуле (Е.2):
mпр = 4,6 · 107 / 4,52 · 106 ·(0,8· 530 · 600) · 0,1 = 2,59 · 105 кг.
Находим избыточное давление Dp по формуле (Е.1)
Dp = 101 [0,8 (2,59 · 105)0,33 / 500 + 3 (2,59 · 105) 0,66 / 5002 + 5 (2,59 · 105) / 5003] = 16,2 кПа.
Находим импульс волны давления i по формуле (Е.3):
i = 123 (2,59 · 105)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)
где p0 — атмосферное давление, кПа (допускается принимать равным 101 кПа);
r — расстояние до разрушающегося технологического оборудования, м;
mпр — приведенная масса, кг, рассчитанная по формуле
mпр = Eиэ / Q0. (Ж.4)
где Eиэ — энергия, выделяющаяся при изэнтропическом расширении среды, находящейся в резервуаре, Дж;
q0 — константа, равная 4,52 · 106 Дж/кг.
Ж.4 Eиэ , Дж, рассчитывают по формуле
Eиэ = Сэфф m ( Т - Ткип). (Ж.5)
где m — масса вещества в резервуаре, кг;
Сэфф — константа, равная 500 Дж/(кг·К);
Т — температура вещества в резервуаре в момент его взрыва, К;
Ткип — температура кипения вещества при атмосферном давлении, К.
При наличии в резервуаре предохранительного клапана Т, К, допускается рассчитывать по формуле
, (Ж.6)
где А, В, Са — константы Антуана вещества;
рк — давление срабатывания предохранительного клапана, кПа. Константа А должна соответствовать давлению, выраженному в килопаскалях.
Пример — Расчет параметров ударной волны при BLEVE
Данные для расчета
Рассчитать параметры положительной фазы волны давления на расстоянии 750 м от эпицентра аварии, связанной с развитием BLEVE на железнодорожной цистерне вместимостью 50 м3 с 10 т жидкого пропана. Цистерна имеет предохранительный клапан на давление срабатывания 2,0 МПа.
Расчет Энергию, выделившуюся при расширении среды в резервуаре, рассчитывают по формуле (Ж.5).
Eиэ = Сэфф m ( Т - Ткип),
где m = 4 · 104 кг — масса пропана в цистерне;
Сэфф — константа, равная 500 Дж/кг·К);
Ткип = - 43 + 273 = 230 К — температура кипения пропана при постоянном давлении.
Т, К, находим по формуле (Ж.6)
где рк = 2,000 кПа, А = 5,949, В = 812,648, Са = 247,55.
Отсюда
Получим Eиэ
Eиэ = 4 · 104(332-230)500 = 2,06 · 109 Дж.
Находим приведенную массу mпр, кг, по формуле (Ж.4)
mпр = 2,06 · 109 / (4,52 · 106 ) = 456 кг.
Вычислим Dp и i по формулам (Ж.2) и (Ж.3)
Dр = 101 (0,8 · 4560,33 / 750 + 3 · 4560,66 / 7502 + 5 · 4563 / 750 ) = 0,86кПа,
i = 123 · 4560,66 / 750 = 9,7 Па · с.
ПРИЛОЖЕНИЕ И (рекомендуемое) МЕТОД РАСЧЕТА ПАРАМЕТРОВ ИСПАРЕНИЯ ГОРЮЧИХ НЕНАГРЕТЫХ ЖИДКОСТЕЙ И СЖИЖЕННЫХ УГЛЕВОДОРОДНЫХ ГАЗОВ
И. 1 Интенсивность испарения W, кг/(с·м2), определяют по справочным и экспериментальным данным. Для ненагретых выше температуры окружающей среды ЛВЖ, при отсутствии данных допускается рассчитывать W по формуле1)
W = 10-6 h pн, (И.1)
_______
1) Формула применима при температуре подстилающей поверхности от минус 50 до плюс 40 °С.
где h — коэффициент, принимаемый по таблице И.1 в зависимости от скорости и температуры воздушного потока над поверхностью испарения;
М — молярная масса, г/моль;
pн — давление насыщенного пара при расчетной температуре жидкости tр, определяемое по справочным данным, кПа.
Таблица И.1
Скорость воздушного потока в помещении,
| Значение коэффициента hпри температуре t, °С, воздуха в помещении
| м/с
|
|
|
|
|
| 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, по формуле1)
, (И.2)
_______
1) Формула применима при температуре подстилающей поверхности от минус 50 до плюс 40 °С.
где М — молярная масса СУГ, кг/моль;
Lисп — мольная теплота испарения СУГ при начальной температуре СУГ Тж, Дж/моль;
Т0 — начальная температура материала, на поверхность которого разливается СУГ, соответствующая расчетной температуре tp, К;
Тж — начальная температура СУГ, К;
lтв — коэффициент теплопроводности материала, на поверхность которого разливается СУГ, Вт/(м · К);
а — эффективный коэффициент температуропроводности материала, на поверхность которого разливается СУГ, равный 8,4·10-8 м2/с;
t — текущее время, с, принимаемое равным времени полного испарения СУГ, но не более 3600 с;
число Рейнольдса (n— скорость воздушного потока, м/с; d — характерный размер пролива СУГ, м;
uв — кинематическая вязкость воздуха при расчетной температуре tр, м2/с);
lв — коэффициент теплопроводности воздуха при расчетной температуре tр , Вт/(м · К).
Примеры — Расчет параметров испарения горючих ненагретых жидкостей и сжиженных углеводородных газов
1 Определить массу паров ацетона, поступающих в объем помещения в результате аварийной разгерметизации аппарата.
Данные для расчета
В помещении с площадью пола 50 м2 установлен аппарат с ацетоном максимальным объемом Vaп = 3 м3. Ацетон поступает в аппарат самотеком по трубопроводу диаметром d = 0,05 м с расходом q, равным 2 · 10-3 м3/с. Длина участка напорного трубопровода от емкости до ручной задвижки l1 = 2 м. Длина участка отводящего трубопровода диаметром d = 0,05 м от емкости до ручной задвижки L2 равна 1 м. Скорость воздушного потока и в помещении при работающей общеобменной вентиляции равна 0,2 м/с. Температура воздуха в помещении tр=20 °С. Плотность r ацетона при данной температуре равна 792 кг/м3. Давление насыщенных паров ацетона рa при tр равно 24,54 кПа.
Расчет
Объем ацетона, вышедшего из напорного трубопровода, Vн.т составляет
м3,
где t — расчетное время отключения трубопровода, равное 300 с (при ручном отключении).
Объем ацетона, вышедшего из отводящего трубопровода Vот составляет
Объем ацетона, поступившего в помещение
Va = Vап + Vн.т + Vот = 3 + 6,04 ·10-1 + 1,96 · 10-3 = 6,600 м3.
Исходя из того, что 1 л ацетона разливается на 1 м2 площади пола, расчетная площадь испарения Sр = 3600 м2 ацетона превысит площадь пола помещения. Следовательно, за площадь испарения ацетона принимается площадь пола помещения, равная 50 м2.
Интенсивность испарения равна:
Wисп = 10-6 · 3,5 · 24,54 = 0,655 · 10-3 кг/(с · м2).
Масса паров ацетона, образующихся при аварийной разгерметизации аппарата т, кг, будет равна
т = 0,655 · 10-3 · 50 · 3600 = 117,9 кг.
2 Определить массу газообразного этилена, образующегося при испарении пролива сжиженного этилена в условиях аварийной разгерметизации резервуара.
Данные для расчета
Изотермический резервуар сжиженного этилена объемом Vи.р.э = 10000 м3 установлен в бетонном обваловании свободной площадью Sоб = 5184 м2 и высотой отбортовки Ноб = 2,2 м. Степень заполнения резервуара a = 0,95.
Ввод трубопровода подачи сжиженного этилена в резервуар выполнен сверху, а вывод отводящего трубопровода снизу.
Диаметр отводящего трубопровода dтp = 0,25 м. Длина участка трубопровода от резервуара до автоматической задвижки, вероятность отказа которой превышает 10-6 в год и не обеспечено резервирование ее элементов, L= 1 м. Максимальный расход сжиженного этилена в режиме выдачи Gж.э = 3,1944 кг/с. Плотность сжиженного этилена rж.э при температуре эксплуатации Тэк = 169,5 К равна 568 кг/м3. Плотность газообразного этилена rг.э при Тэк равна 2,0204 кг/м3. Молярная масса сжиженного этилена Мж.э = 28 · 10-3 кг/моль. Мольная теплота испарения сжиженного этилена Lиcn при Тэк равна 1,344 · 104 Дж/моль. Температура бетона равна максимально возможной температуре воздуха в соответствующей климатической зоне Tб = 309 К. Коэффициент теплопроводности бетона lб=1,5Вт/(м·К). Коэффициент температуропроводности бетона а = 8,4 · 10-8 м2/с. Минимальная скорость воздушного потока umin = 0 м/с, а максимальная для данной климатической зоны umax = 5 м/с. Кинематическая вязкость воздуха nв при расчетной температуре воздуха для данной климатической зоны tр = 36 °С равна 1,64 · 10-5 м2/с. Коэффициент теплопроводности воздуха lв при tр равен 2,74 · 10-2 Вт/(м · К).
Расчет
При разрушении изотермического резервуара объем сжиженного этилена составит
м3.
Свободный объем обвалования Vоб = 5184 · 2,2 = 11404,8 м3.
Ввиду того, что Vж.э < Vоб примем за площадь испарения Sисп свободную площадь обвалования Sоб, равную 5184 м2.
Тогда массу испарившегося этилена mи.э с площади пролива при скорости воздушного потока u = 5 м/с рассчитывают по формуле (И.2)
Масса mи.э при u = 0 м/с составит 528039 кг.
ПРИЛОЖЕНИЕ К (рекомендуемое) МЕТОДЫ РАСЧЕТА ТЕМПЕРАТУРНОГО РЕЖИМА ПОЖАРА В ПОМЕЩЕНИЯХ ЗДАНИЙ РАЗЛИЧНОГО НАЗНАЧЕНИЯ
К.1 Условные обозначения
V— объем помещения, м3;
S— площадь пола помещения, м2;
Аi — площадь i-го проема помещения, м2;
hi — высота i-го проема помещения, м;
— суммарная площадь проемов помещения, м2;
— приведенная высота проемов помещения,м;
П— проемность помещения, рассчитывается по формуле (К.1) или (К.2), м0,5;
Рi — общее количество пожарной нагрузки i-го компонента твердых горючих и трудногорючих материалов, кг;
q — количество пожарной нагрузки, отнесенное к площади пола, кг/м;
qкр.к — удельное критическое количество пожарной нагрузки, кг/м2;
qк — количество пожарной нагрузки, отнесенное к площади тепловоспринимающих поверхностей помещения, кг/м2;
Пср — средняя скорость выгорания древесины, кг/(м2 · мин);
Псрi — средняя скорость выгорания i-го компонента твердого горючего или трудногорючего материала, кг/м2 · мин);
— низшая теплота сгорания древесины, МДж/кг;
Не нашли, что искали? Воспользуйтесь поиском по сайту:
©2015 - 2024 stydopedia.ru Все материалы защищены законодательством РФ.
|