CC BY
303
УДК 519.6
Ю. Л. ЗАЕЦ, В. В. БЕЛЯЕВА (ДИИТ)
РАСЧЕТ ОБРАЗОВАНИЯ ВЗРЫВООПАСНОЙ КОНЦЕНТРАЦИИ В ПОМЕЩЕНИИ ПРИ АВАРИЙНОЙ УТЕЧКЕ ГАЗА
Розроблено тривимiрну чисельну модель розрахунку процесу розповсюдження забруднюючих речовин у виробничому примщенш при аварiйному викидi токсичного газу. Модель базуеться на чисельному штегру-ваннi рiвняння конвективно-дифузiйного переносу домшки та моделi течи нестисло! рвдини. Модель дозво-ляе розрахувати час утворення вибухонебезпечно! концентрацп у примщенш при аваршному викидi газу. Наводяться результати обчислювального експерименту по прогнозу рiвня забруднення повiтряного середо-вища у виробничому примiщеннi.
Разработана трехмерная численная модель расчета процесса распространения загрязняющих веществ в производственном помещении при аварийном выбросе токсичного газа. Модель основывается на численном интегрировании уравнения конвективно-диффузионного переноса примеси и на модели невязкой несжимаемой жидкости. Модель позволяет рассчитывать время образования взрывоопасной концентрации в помещении при аварийной утечке газа. Приводятся результаты вычислительного эксперимента по прогнозу уровня загазованности воздушной среды производственного помещения.
The 3D numerical model to simulate the toxic gas dispersion in industrial room after accident ejections was developed. The model is based on the K-gradient transport model and equation of potential flow. The model permits to calculate the time of appearance of concentration when the blast can take place. The results of numerical experiment of air pollution in industrial room are presented.
Введение
Одной из важных задач в области экологической безопасности является повышение качества прогнозной информации при распространении токсичных веществ в производственных помещениях после аварии с целью оценки масштабов возможных последствий аварий, выявления в производственном помещении зон взрывоопасных концентраций. Используемые на практике инженерные методики прогноза уровня загрязнения воздушной среды в производственных помещениях основываются, как правило, на аналитическом решении одномерного уравнения переноса примеси либо на применении балансового соотношения «приток -отток - накопление массы» для помещения (нульмерная модель).
На практике расчет динамики изменения уровня загазованности помещения при работе аварийной вентиляции осуществляют на базе зависимости [6]
С = Спр + — + Kp
с — с
пр K
exp(-tKp)
, (1)
где Спр - концентрация загрязнителя в приточном воздухе, мг/м ; Kp - кратность воздухообмена;
С — концентрация загрязнителя, мг/м3;
t - время, ч; _ —
— = _ — удельное, часовое выделение загряз-V
нителя;
V - объем помещения, м3;
G - интенсивность выделения загрязнителя в воздушную среду помещения, мг/ч; С - начальная концентрация загрязнителя, мг/м3.
Данная зависимость представляет собой аналитическое решение уравнения материального баланса загрязнителя в помещении (нульмерная модель) [6]
г§ = щ„+ с - 1С,
где Ь — воздухообмен, м3/ч.
Относительно приведенной выше расчетной зависимости, следует отметить следующее:
• модель определяет концентрацию загрязнителя С (токсичного газа) не в помещении, где произошла авария, а в удаляемом воздухе. Это не дает возможности проектировщику оценить степень загазованности воздуха в различных местах помещения, особенно в плохо проветриваемых местах (за технологическим оборудованием) и тем самым дезориентирует проектировщика относительно возможности возникновения вторичной аварии;
• модель не позволяет учесть влияние положения приточных и вытяжных отверстий на организацию движения воздушного потока внутри помещения и, следовательно, на динамику формирования зоны загрязнения в помещении, где произошла авария;
• модель не дает возможности учесть влияние технологического оборудования на процесс рассеивания токсичного газа в помещении.
Отметим, что указанная расчетная зависимость используется на практике как для прогноза уровня загазованности производственных помещений при аварии (при отсутствии работы аварийной вентиляции), так и при работе аварийной вентиляции.
На практике одним из важнейших вопросов является, в частности, расчет времени, когда будет достигнута взрывоопасная концентрация в помещении при аварийной утечке. Для решения этой задачи обычно используются эмпирические зависимости и расчетные формулы, основанные на балансовых соотношениях. Например, время образования взрывоопасной концентрации в помещении при неработающей вентиляции определяется по зависимости [3]
т = -
0.5CnV £0
где §о — величина утечки газа; V - объем помещения; Сп — нижний концентрационный
предел воспламенения.
Очевидно, что данная зависимость представляет собой балансовое соотношение и дает оценку времени, когда величина нижнего предела взрывоопасной концентрации будет достигнута в половине объема помещения. Данная зависимость не позволяет выявить зоны в помещении, где может быть достигнута взрывоопасная концентрация (например, в застойных зонах за технологическим оборудованием). При наличии вентиляции время образования взрывоопасной концентрации в помещении рассчитывается по зависимости [3]:
T=3600ln
£0
£0-
CnnV 72-104
1. Математическая модель. Для моделирования процесса миграции загрязняющего вещества в производственном помещении будем использовать трехмерное уравнение переноса примеси [3, 4]
дС диС дС ди—^) С д С — +-+--=— —)+
д дх ду дг дх дх
§>тг-*. (2)
где С - концентрация загрязняющего вещества в помещении; и, V, ^ - компоненты вектора скорости воздушной среды; wS - скорость оседания примеси; р = (рх, ру, - коэффициент турбулентной диффузии; Q - интенсивность выброса токсичного вещества в помещении; д(г — г ) - дельта-функция Дирака; г,= (х,, у, - координаты источника выброса.
Для расчета поля скорости воздушного потока в помещении, индуцированного работой вентиляции, делается допущение, что движение воздушной среды в помещении - потенциальное. В этом случае компоненты скорости воздушной среды определяются соотношениями
u =
dP_ dx
v = ■
cp_
dy
w = ■
cp
где Р - потенциал.
Уравнение для определения потенциала имеет вид
д P д P д P п —т + —т + —г = 0.
dx2 dy2 dz
(3)
где п - кратность воздухообмена.
Данная зависимость также дает лишь приближенную оценку для определения такой важной, с точки зрения экологической безопасности, величины, как т .
В этой связи особую актуальность приобретает разработка теоретических методов прогноза распространения токсичных веществ в производственных помещениях химически опасных объектов, позволяющих учесть те существенные факторы, которые влияют на данный процесс и многочисленные сценарии аварийных выбросов (утечки, разливы и т.д.). Применение таких методов на практике позволит существенно повысить качество прогнозов и может служить научным обоснованием для принятия инженерных решений, направленных на минимизацию экологического ущерба и последствий аварий в производственных помещениях химически опасных объектов.
Постановка краевых условий для уравнения (1) рассмотрена в работе [3]. Для уравнения (3) ставятся следующие граничные условия:
дР
• на твердых стенках -= 0, где п -
дп
единичный вектор внешней нормали;
• на входной границе (границы втекания воздушного потока в помещение)
дР
— = Vn, где Vn - известное значение дп
скорости;
• на выходной границе Р = Р * (x = const, y) + const. (условия Дирихле).
2. Метод решения. Численное интегрирование уравнения (2) осуществляется с использованием неявной попеременно-треугольной разностной схемы расщепления [3]. Для численного интегрирования уравнения (3) используется идея установления решения по времени, т. е. интегрируется уравнение вида
п
дР
__дР+д2р + д2 р
дт дх2 дУ2 дг1
где т - фиктивное время.
Численное интегрирование данного уравнения проводится с использованием неявной разностной схемы условной аппроксимации. В разработанном пакете программ осуществляется также численное интегрирование уравнения (3) с помощью метода Либмана. Расчет поля потенциала проводится одновременно с использованием указанным методов с целью контроля результатов расчета.
3. Практическая реализация. Как уже отмечалось, одним из опасных последствий аварий может быть вторичная авария - взрыв, пожар в помещении, если в нем образовались зоны, где концентрации паров находится в определенных пределах и есть источник воспламенения [2]. Рассмотрим применение разработанной трехмерной численной модели для решения этой задачи. Параметры задачи: размеры помещения 10 м х 10 м х 10 м; воздухообмен - 2.77 м3/с; коэффициент диффузии - 0.1 м3/с; интенсивность выброса аммиака -G _ 97 г/с, выброс осуществляется на временном интервале [0...5 с]; координаты источника выброса: х = 4.5 м; у = 4.5 м; г = 2.5 м; шаг по времени йг = 0.005 с; подача воздуха и отвод его из помещения осуществляется через одно приточное и одно вытяжное отверстие, расположенное на противоположных стенках; концентрация аммиака в приточном воздухе Спр _ 0. Расчет выполняется на «грубой» (реальной) разностной сетке 12х12х12 узлов.
Результаты решения данной задачи представлены в табл. 1. Здесь же показаны результаты расчета величины концентрации аммиака в помещении, полученные по аналитической модели (нульмерная модель) (1).
Таблица 1
Динамика загрязнения воздушной среды в помещении при аварийной утечке токсичного газа
г, с Стах, г/м3 (расчет по разработанной трехмерной численной модели) С, г/м3 (расчет по аналитической модели)
2 98.97 0.19
3 107.30 0.29
4 112.90 0.38
5 115.57 0.48
Как видно из анализа данных, представленных в табл. 1, динамика загрязнения воздушной среды в помещении при аварийной утечке токсичного газа характеризуется ростом концентрации аммиака на рассматриваемом временном интервале. Обращает на себя внимание существенное различие значения концентра-
ции, определенной на базе аналитической модели и с помощью трехмерной численной модели. Этот факт объясняется тем, что в применяемой на практике аналитической модели, расчет выполняется не по реальной мощности источника выделения токсичного газа G, а по усредненной по объему помещения величине выброса. Таким образом, применение аналитической модели дает заниженное прогнозное значение концентрации токсичного газа в помещении при аварии. Отметим, например, что взрывоопасная концентрация для аммиака находится в диапазоне 107.200 г/м3, таким образом, этот предел уже будет достигнут через 3 с после аварии. Как видно, аналитическая модель не дает этого результата, что может дезориентировать проектировщика.
Выводы. В работе предложена трехмерная численная модель, на основе которой построен метод расчета динамики загрязнения воздушной среды в помещении при аварийной утечке токсичного газа.
Данная модель обладает рядом преимуществ, крайне важных для проектировщика: возможность учета влияния расположения технологического оборудования в помещении, положения приточных и вытяжных отверстий вентиляции, возможность моделирования различной формы облака токсичного газа, которое формируется на месте аварии; возможность моделирования выброса, утечки токсичного газа в любом месте помещения. Предложенная модель позволила впервые решить важную задачу в области экологической безопасности -прогнозирование процесса распространения загрязняющих веществ в производственных помещениях при авариях. Дальнейшее развитие данной модели должно быть направлено на создание численной модели для расчета рассеивания тяжелых газов в производственных помещениях.
БИБЛИОГРАФИЧЕСКИЙ СПИСОК
1. Аварии и катастрофы. Предупреждение и ликвидация последствий. Уч. пос. в 5 кн. / Под ред. В. А. Котляревского и А. В. Забегаева. - М.: Изд-во АСВ, 2001. - 200 с.
2. Згуровский М. З. Численное моделирование распространения загрязнения в окружающей среде / М. З. Згуровский, В. В. Скопецкий, В. К. Хрущ, Н. Н. Беляев. - К.: Наук. думка, 1997. - 368 с.
3. Крикунов Г. Н. Безопасность жизнедеятельности / Г. Н. Крикунов, А. С. Беликов, В. Ф. Залу-нин, В. Н. Довгаль. - Д.: УкО ИМА-пресс, 1995. - 196 с.
4. Марчук Г. И. Математическое моделирование в проблеме окружающей среды. - М.: Наука, 1982. - 320 с.
5. Самарский А. А. Теория разностных схем / 2-е изд., испр. - М.: Наука, 1983. - 616 с.
6. Эльтерман В. М. Вентиляция химических производств / 3-е изд., перераб. - М.: Химия, 1980. -288 с.
Поступила в редколлегию 05.02.2008.
