Гомель 2004

Таблица 2

Примечание:

склад 1 - кирпичное здание;

склад 2 - здание из сборного железобетона;

административное - многоэтажное здание с металлическим каркасом;

ВНБ - водонапорная башня;

К-150 - наземный кольцевой трубопровод;

ЛЭП - воздушная линия электропередачи низкого напряжения.

1. Расчёт степени воздействия ударной волны на объекты и человека при детонационном взрыве газо-паровоздушного облака

Участок А

Описание горючего вещества:

Диэтиловый эфир этиловый эфир, этоксиэтан ()

Физико-химические свойства: жидкость. Молярная масса 74,12 г/моль плотность 713,5кг/м³ при 20^оС; плотность пара по воздуху 2,6; температура кипения +34,5^оС; коэффициент диффузии пара в воздухе 0,0772 см²/с (расч); теплота образования -1252,2 кДж/моль; теплота сгорания-2531 кДж/моль; в воде растворим.

Пожароопасные свойства: Легковоспламеняющаяся жидкость. Температура вспышки -41^оС; температура самовоспламенения 180^оС; концентрационные пределы распространения пламени: 1,7-49 % (об.); температурные пределы воспламенения: НТПВ -44^оС, ВТПВ 16^оС; минимальное содержание кислорода для диффуз. горения 15,4 %, скорость выгорания 10,83*; максимальная нормальная скорость распространения пламени 0,49 м/с; максимальное давление взрыва 720МПа; минимальная энергия зажигания 0,2 мДж;

Средства тушения: Вода в виде компактных или распылённых струй, пены, порошок ПСБ-3.

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

Масса паров диэтилового эфира в облаке складывается из массы вещества, испарившегося с поверхности разлива m_исп и массы насыщенных паров диэтилового эфира m_н.п., содержащихся изначально в свободном объёме резервуара до момента аварии:

m_г = m_исп + m_н.п.

1.1 Расчёт массы насыщенных паров горючего в резервуаре

Объём насыщенного диэтилового эфира в резервуаре:

, м³ 1.2

где V_рез - объём резервуара, м³4

V_ж - объём горючей жидкости в резервуаре, м³;

- объёмная доля насыщенного пара горючей жидкости (в долях от единицы).

Объём горючей жидкости в резервуаре:

, м³

где m_ж - масса жидкости в резервуаре, кг;

_ж - плотность диэтилового эфира при температуре 10^оС (по справочной таблице приложения 4 находим плотность диэтилового эфира а при 10^оС, равную 713,5 кг/м³).

Парциальное давление насыщенного пара диэтилового эфира при температуре 10оС:

,

где А, В, С_А - константы уравнения Антуана, равные соответственно 3,6875, 903,588 и -66,69 (приложение 1).

Объёмная доля насыщенного пара диэтилового эфира при 10С:

где P_о - атмосферное давление (по условию равно 99,5кПа);

Следовательно, в данных условиях:

Молярный объём при температуре 10С (283) и давлении 102,5кПа

По формуле Менделеева-Клапейрона:

Масса насыщенных паров горючей жидкости в резервуаре:

где V_н.п. - объём насыщенных паров жидкости, м³;

V_м - молярный объём при заданных условиях, м³/кмоль.

1.2 Расчёт массы жидкости, испарившейся с поверхности разлива

Интенсивность испарения диэтилового эфира в неподвижную среду:

где Р_н.п. - парциальное давление насыщенных паров жидкости, кПа (рассчитывается по формуле 1.5);

М - молярная масса жидкости, кг/кмоль.

Масса диэтилового эфира , испарившегося с поверхности разлива

Так как, согласно исходным данным, разлив диэтилового эфира происходит в пределах обвалования, площадь которого меньше площади разлива на неограниченную поверхность, то принимаем, что площадь разлива равна площади обвалования: Fпр = Fобв. Следовательно:

где W_исп - интенсивность испарения;

_исп - продолжительность испарения до момента воспламенения облака;

F_обв - площадь обвалования.

Следовательно, по формуле 1.1 масса диэтилового эфира в облаке ГПВС равна:

1.3 Расчёт тротилового эквивалента при детонационном взрыве облака ГПВС

Стехиометрическая концентрация паров диэтилового эфира в облаке ГПВС:

Уравнение полного горения паров диэтилового эфира в воздухе:

+6(О₂ + 3,76N₂) 4СО₂ + 5Н₂О +63,76N₂

Из уравнения реакции видно, что газопаровоздушной смесью в данном случае является смесь + 6(О₂ + 3,76N₂) - т. е. левая часть уравнения. Следовательно:

где n_г - стехиометрический коэффициент перед горючим веществом;

- сумма стехиометрических коэффициентов компонентов исходной ГПВС.

Объём облака ГПВС:

где V_м - молярный объём при заданных условиях;

m_г - масса диэтилового эфира в облаке ГПВС.

Плотность стехиометрического облака ГПВС:

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

М - молярная масса горючего вещества, кг/кмоль.

Низшая теплота сгорания:

диэтиловый эфир является химическим веществом, имеющим определённую химическую формулу

Для индивидуальных химических веществ низшую теплоту горения определяется через нахождение энтальпии реакции их сгорания. Энтальпия сгорания вещества есть изменение энтальпии реакции его горения в расчёте на 1 моль (или 1 кмоль).

Энтальпия сгорания:

Для реакции полного горения паров диэтилового эфира

где , и - стандартные энтальпии образования углекислого газа, паров воды и паров диэтилового эфира соответственно.

Согласно справочным данным: = -393,5 кДж/моль,

= -242,5 кДж/моль, = -276,96 кДж/моль.

Следовательно:

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

Массовая теплота сгорания стехиометрической ГПВС:

где - плотность стехиометрической смеси.

Тротиловый эквивалент:

где V_ГПВС - объём образовавшегося облака ГПВС, м³;

4184 - теплота взрывчатого разложения 1 кг тротила, кДж/кг.

1.4 Расчёт степени воздействия ударной волны на различные объекты

1) Склад 1 - кирпичное здание (r₁ = 100 м):

Приведённый радиус (по формуле 1.19):

Приведённое давление во фронте ударной волны (по формуле 1.20):

Следовательно:

Избыточное давление во фронте ударной волны (по формуле 1.21):

;

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

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

Избыточное давление отражения (по формуле 1.25):

где Р_о - атмосферное давление, равное по условию 102,5кПа,

Согласно справочной таблице приложения 8, для складского кирпичного здания:

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

2) Склад 2 - здание из сборного железобетона (r₂ =105 м):

Приведённый радиус:

Приведённое давление во фронте ударной волны:

Следовательно:

Избыточное давление во фронте ударной волны:

;

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

Избыточное давление отражения :

Согласно справочной таблице приложения 8, для здания из сборного железобетона:

Следовательно, на расстоянии 105 м от эпицентра взрыва облака ГПВС помещение склада 2 подвергнется слабым разрушениям: будет разрушена кровля, фонари, оконные переплёты, незначительно повреждены внутренние малопрочные перегородки, разрушено 100% остекления. Значительного физического ущерба людям, находящимся в помещении склада, причинено не будет.

3) Административное здание - многоэтажное здание с металлическим каркасом (r₄ = 75 м):

Приведённый радиус:

Приведённое давление во фронте ударной волны:

Следовательно:

Избыточное давление во фронте ударной волны:

;

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

Избыточное давление отражения :

Согласно справочной таблице приложения 8, для здания с металлическим каркасом:

Следовательно, на расстоянии 75м от эпицентра взрыва облака ГПВС помещение административного здания подвергнется средним разрушениям: повреждены оконные проёмы, разрушено 100% остекления, происходит изгибание металла. Значительного физического ущерба людям, находящимся в помещении административного здания, причинено не будет.

4) Водонапорная башня (r₅ =60м):

Приведённый радиус:

Приведённое давление во фронте ударной волны:

Следовательно:

Избыточное давление во фронте ударной волны:

;

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

Давление скоростного напора (по формуле 1.23):

,

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

Лобовое давление (формула 1.24):

Согласно справочной таблице приложения 8, для металлической водонапорной башни:

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

5) Наземный кольцевой трубопровод К-150 (r₆ = 50 м):

Приведённый радиус:

Приведённое давление во фронте ударной волны:

Следовательно:

Избыточное давление во фронте ударной волны:

;

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

Согласно справочной таблице приложения 8, для наземного трубопровода:

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

6) Грузовой автопарк (r₇ = 75 м):

Приведённый радиус:

Приведённое давление во фронте ударной волны:

Следовательно:

Избыточное давление во фронте ударной волны:

;

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

Следовательно, на расстоянии 75 м от эпицентра взрыва облака ГПВС грузовая автотехника подвергнется очень слабым разрушениям: частичное разрушение остекления. Опрокидывания техники не будет. Текущего (среднего) ремонта практически не требуется.

7) Линии электропередачи: (r = 100-90 м)

Для воздушной линии электропередачи низкого напряжения, согласно справочным данным (приложение), слабые разрушения наблюдаются при величине Р_ф 20-60 кПа. На расстоянии 90 м избыточное давление во фронте ударной волны составляет 60 кПа (согласно расчётам в пункте 4). Следовательно, при взрыве ГПВС в данной обстановке может произойти обрыв ЛЭП

Зоны разрушений:

В результате бризантного и фугасного воздействия ударной волны условно можно выделить четыре зоны разрушений: зону полных R_I (ДР_ф > 50 кПа), зону сильных R_II (ДР_ф = 30 - 50 кПа), зону средних R_III (ДР_ф = 22 - 30 кПа) и зону слабых R_IV (ДР_ф = 10 - 22 кПа) разрушений. При значении ДР_ф < 10 кПа можно выделить зону лёгких разрушений R_V (ДР_ф = 1 - 10 кПа).

Коэффициент разрушений К_i:

Радиусы зон разрушений (по формуле 1.26):

;

;

;

;

.

Радиус облака ГПВС (по формуле 1.27):

где F_пр - площадь пролива жидкости, м²

Из проведённых расчётов видно, что радиусы зон разрушений R_I, R_II, меньше радиуса облака ГПВС.

1.5 Расчёт вероятности поражения человека ударной волной.

Группа людей (r₃ = 50 м):

Приведённый радиус:

Приведённое давление во фронте ударной волны:

Следовательно:

Избыточное давление во фронте ударной волны:

;

Основными поражающими факторами при воздействии ударной волны на человека является избыточное давление во фронте ударной волны ДР_ф, импульс I и осколочное действие. Величина ДР_ф обуславливает барическое воздействие на человека (т. е. баротравму), а I - метательное воздействие (т. е. опрокидывание или отброс).

Приведённый импульс во фронте ударной волны (формула 1.28):

Следовательно:

Импульс во фронте ударной волны (формула 1.30):

«Пробит» - функция (формула 1.31):

где:

Следовательно:

Согласно справочной таблице приложения 10, рассчитанное значение «пробит»-функции соответствует вероятности поражения человека Q_вп 91,5%.

1.6 Расчёт минимального безопасного расстояния человека от эпицентра взрыва

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

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

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

Участок А

2.1 Расчёт массовой скорости выгорания диэтилового эфира. Линейная скорость выгорания (по формуле 2.1):

Горение горючей жидкости в пределах обвалования можно представить как горение в резервуаре большого размера. Для резервуаров достаточно крупных размеров (d > ?):

где М - молярная масса горючей жидкости, кг/кмоль;

с_ж - плотность горючей жидкости при температуре окружающей среды, кг/м³;

л - коэффициент теплопроводности жидкости, кал/(м•с•град);

g - ускорение свободного падения, равное 9,81 м/с²;

Т_кип - температура кипения жидкости, К;

с_р - теплоёмкость жидкости, кал/(моль•К);

з - динамическая вязкости жидкости, П (Пуаз);

ДН_исп - теплота испарения жидкости, кал/моль

Q_н - молярная теплота сгорания жидкости, кал/моль.

Учитывая, что 1 кал = 4,1865 Дж, для диэтилового эфира

Значение Q_н в кДж/моль рассчитано в разделе 1.3.

Остальные данные берутся из справочных таблиц приложения 1-3.

Следовательно, для диэтилового эфира при t = 10^оС:

Массовая скорость выгорания диэтилового эфира (по формуле 2.2):

где с_ж - плотность жидкости, кг/м³.

2.2 Расчёт плотности теплового потока на различных расстояниях от эпицентра горения

Массовая низшая теплота сгорания диэтилового эфира (по формуле 2.4):

,

где - молярная низшая теплота сгорания диэтилового эфира, рассчитанная в разделе 1.3.

Интенсивность тепловыделения факела пламени (по формуле 2.3):

,

где х_m - массовая скорость выгорания диэтилового эфира , кг/(м²?с);

F_пр - площадь поверхности разлива диэтилового эфира , м².

Интенсивность излучения факела пламени (по формуле 2.5):

Доля тепла, расходуемого на излучение факела пламени б_изл, для диэтилового эфира, можно принять равной 0,25:

Условная высота пламени при горении разлива диэтилового эфира:

+6(О₂ + 3,76N₂) 4СО₂ + 5Н₂О +63,76N₂

Стехиометрическое отношение воздуха к летучим продуктам горения:

По формуле 2.8:

где - коэффициент перед кислородом в уравнении реакции полного горения диэтилового эфира.

- коэффициент перед i-м продуктом горения в уравнении реакции.

Диаметр очага горения

Удельная теплоёмкость воздуха:

Плотность воздуха при t = 10^оС:

Выражение безразмерного параметра (по формуле 2.7):

где - удельная теплоёмкость воздуха, Дж/(кг?К);

Т_в и с_в - температура и плотность окружающего воздуха, К и кг/м³ соответственно;

- массовая теплота сгорания диэтилового эфира , кДж/кг;

g - ускорение свободного падения (9,81 м/с²);

q - интенсивность тепловыделения факела пламени, кВт;

d - диаметр очага горения, м;

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

Условная высота пламени (по формуле 2.6):

Коэффициент пропускания атмосферы для расстояния r = 50 м (по формуле 2.10):

где ф_атм - коэффициент пропускания атмосферы для теплового излучения;

r - расстояние от группы людей до эпицентра очага горения, м.

Коэффициент облучённости людей факелом пламени

Из рисунка видно, что величина a соответствует радиусу пролива жидкости (r_пр= 21,85 м), b - расстояние r от заданного объекта или человека до эпицентра горения (r₃ = 50м)

Следовательно:

;

Плотность теплового потока на расстоянии 50м от эпицентра горения (по формуле 2.13):

где - интенсивность излучения факела пламени, кВт (формула 2.5);

2.2 Воздействие теплового излучения на объекты и человека

Согласно таблицам приложения 14-17, критические плотности теплового потока для большинства строительных материалов и конструкций, а также для техники, намного превышают рассчитанную величину для r = 50м. Следовательно, можно сделать вывод, что в данной обстановке помещение склада 1, склада 2, административного здания, водонапорная башня и автотехника находятся на безопасном удалении от очага горения.

Безопасная плотность теплового потока для людей без специальной защиты составляет 1,39 кВт/м². Следовательно, в данном случае группа людей находится в зоне безопасного теплового воздействия, так как < 1,39 кВт/м².

Минимальное безопасное расстояние человека от эпицентра очага горения:

1) Приблизительное минимальное безопасное расстояние (по формуле 2.14):

2) Коэффициент облучённости для расстояния :

3) Коэффициент пропускания атмосферы при = 48м;

4) Плотности теплового потока на расстоянии :

5) Минимальное безопасное расстояние (методом линейной интерполяции по формуле 2.15):

Так как согласно исходным данным, люди находятся на расстоянии, большем чем , то для них вероятность теплового поражения человека Q_вп 0%..

2.3 Минимальное безопасное расстояние для прибывших подразделений

Для людей в специальной защитной одежде, к которой относится в том числе и боевая одежда спасателей, критическая плотность теплового потока равна 4,2 кВт/м². Расстояние от эпицентра горения, которому соответствует данная плотность теплового потока, рассчитывается по формулам 2.14 - 2.15 методом последовательных приближений с учётом того, что в формулу 2.14 подставляется 4,2 кВт/м²:

1) Приблизительное минимальное безопасное расстояние (по формуле 2.14):

2) Коэффициент облучённости для расстояния :

3) Коэффициент пропускания атмосферы при = 30м:

4) Плотности теплового потока на расстоянии :

5) Минимальное безопасное расстояние (методом линейной интерполяции по формуле 2.15):

Согласно расчётам в разделе 2.2, при r = 48м = 1,5кВт/м².

Следовательно:

ГЛАВА II. УЧАСТОК «В»

Таблица 3. Исходные данные

2. Расчёт степени теплового воздействия на объекты и человека при горении огненного шара

Участок В