Распространение загрязняющих веществ в атмосфере в условиях застройки

Размещено на http://www.allbest.ru/

Распространение загрязняющих веществ в атмосфере в условиях застройки

Н.Н. Беляев

Е.Ю. Гунько

Т.М. Пустоварова

В последнее время особый интерес приобретают задачи, связанные с прогнозом токсичного поражения людей на промплощадках в случае аварий, сопровождающихся выбросом в атмосферу токсикантов. Для решения этих задач в настоящее время не предоставляется возможным использование моделей гидродинамики типа модель, LES модель и т. д. Это связано с тем, что практическое применение таких моделей требует мелкой расчетной сетки, что приводит к значительным затратам машинного времени при проведении вычислительного эксперимента. В этой связи для прогноза уровня загрязнения воздушной среды в условиях застройки на практике обычно используется эмпирические модели. Их существенным недостатком является то, что модели не позволяют учесть влияние зданий на процесс распространения токсичных веществ.

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

Математическая модель

Для расчета процесса переноса токсичного газа в условиях застройки будем использовать трехмерное уравнение миграции примеси (модель градиентного типа) [2-4]

, (1)

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

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

,

где - потенциал.

Уравнение для определения потенциала имеет вид

. (2)

Постановка краевых условий для уравнения (1) рассмотрена в работе [3]. Для уравнения (2) ставятся следующие граничные условия.

- на твердых стенках

,

где - единичный вектор внешней нормали;

- на входной границе (границы втекания воздушного потока в помещение)

,

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

- на выходной границе

.

(условия Дирихле).

Метод решения

Численное интегрирование уравнения (1) осуществляется с использованием неявной попеременно-треугольной разностной схемы расщепления [3]. Особенностью данной разностной схемы является то, что при решении уравнения (1) расчет разбивается на четыре шага. На каждом шаге расщепления строится неявная разностная схема. Однако, определение неизвестного значения концентрации токсичного вещества определяется по явной формуле бегущего счета.

Для численного интегрирования уравнения (2) используется идея установления решения по времени, т.е. интегрируется уравнение вида

,

где - фиктивное время.

Численное интегрирование данного уравнения проводится с использованием неявной разностной схемы условной аппроксимации. В разработанном пакете программ осуществляется также численное интегрирование уравнения (2) с помощью метода Либмана. Расчет поля потенциала проводится одновременно с использованием указанным методов с целью контроля результатов расчета.

Практическая реализация

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

Размещено на http://www.allbest.ru/

Рисунок 1. Схема размещения зданий на промплощадке: 1, 2 - здания; 3 - грибовидное облако токсичного газа

Рассматривается следующая постановка задачи. На промплощадке в результате аварии происходит залповый выброс аммиака. В результате выброса на месте аварии образуется грибовидное облако (см. рисунок 1), имеющее размеры: высота облака -15 м, ширина «ножки» - 10 м, ширина «шляпки» 30 Ч 30 м. Размеры расчетной области 120 Ч 120 Ч 60 м; скорость ветра 3,5 м/с; коэффициент турбулентной диффузии 2 м2/с. Профиль ветра - равномерный. Средняя концентрация аммиака в облаке составляет 1272 г/м3, т.е. в облаке находится около 7 тонн токсичного газа.

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

,

где С - концентрация токсичного газа в точке расположения рецептора (человека) на маршруте эвакуации; t - время.

Размещено на http://www.allbest.ru/

Рисунок 2. Схема размещения объектов (план): 1, 2 - здания; 3 - облако токсичного газа; ----- -маршрут эвакуации людей

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

.

Перейдем к рассмотрению результатов проведенного вычислительного эксперимента. Рассмотрим сначала вопрос о токсичном поражении людей при их эвакуации. Будем считать, что персонал передвигается по маршруту эвакуации со скоростью 3 м/с. Рассматриваются два сценария: в первом сценарии полагается, что эвакуация людей по маршруту (от торца первого здания) началась сразу после залпового выброса, т.е. в момент времени t=0.

Динамика изменения величины токсодозы для первого сценария выглядит следующим образом:

t = 5 c…………………TD = 2,33 мг· мин/л;

t = 7,5 с……………….TD = 2,33 мг· мин/л;

t = 10 с………………..TD = 7,88 мг· мин/л;

t = 15 с………………..TD = 13,3 мг· мин/л;

t = 15 с………………..TD = 15,64 мг· мин/л;

t = 20 с………………..TD= 17,69 мг· мин/л;

t = 22,5 с……………...TD = 18,62 мг· мин/л.

Во втором сценарии полагается, что эвакуация людей из первого здания началась в момент времени t = 8 c, т.е. это время понадобилось на то, чтобы люди покинули здание. В этом случае величина токсодозы изменяется следующим образом:

t = 10 c………………..TD = 1,11 мг· мин/л;

t = 12,5 с……………...TD = 7,26 мг· мин/л;

t = 15 с………………..TD = 12,55 мг· мин/л;

t = 17,5 с……………...TD = 17,06 мг· мин/л;

t = 20 с………………..TD = 21,10 мг· мин/л;

t = 22,5 с……………...TD = 25,00 мг· мин/л.

Если учесть, что токсичное поражение людей может наступить при величине токсодозы TD = 10 мг. мин/л, то становится очевидным, что на данном маршруте эвакуации уже через 15 с после аварии будет иметь место токсичное поражение эвакуируемых людей.

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

Рисунок 3. Зона загрязнения для момента времени t=0

На рисунке 4 показано как токсичное облако вплотную подошло к первому зданию и начинает формироваться зона загрязнения за этим зданием. Из рисунков 5 - 6 хорошо видно, как интенсивно происходит загрязнение воздушной среды возле всех зданий на проплощадке после аварии. Как видно из рисунка 7 спустя 80 с после аварии загрязненной останется область за вторым зданием, что создаст угрозу токсичного поражения людей, покидающих это здание.

Рисунок 4. Зона загрязнения для момента времени c

Рисунок 5. Зона загрязнения для момента времени t=8 c

Рисунок 6. Зона загрязнения для момента времени t=23 c

Рисунок 7. Зона загрязнения для момента времени t=80 c

Таким образом, результаты проведенного вычислительного эксперимента позволяет сделать вывод о том, что аварии аналогичного типа на промплощадках способны вызвать крайне негативные последствия вплоть до смертельного поражения людей. Хорошо видно, что уровень загрязнения воздушной среды в помещении начинает возрастать через 14 с после аварии. На рисунках 3-7 показана динамика формирования зоны загрязнения на промплощадке для различных моментов времени, соответствующих временному интервалу выделения аммиака. Хорошо видно, что формирующаяся зона загрязнения может быть условно разбита на две части. Первая часть - это шлейф, простирающийся над зданиями. Вторая часть - это подзоны, которые образуются между зданиями и создающие опасность токсичного поражения людей, покидающих здания.

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

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

Перечень ссылок

загрязняющий вещество токсичный аварийный

1. Аварии и катастрофы. Предупреждение и ликвидация последствий. Учеб. пособ. в 5-ти книгах / Под редакцией В.А. Котляревского и В. Забегаева). - М.: Из-во АСВ, 2001 - 200 с.

2. Згуровский М.З. Численное моделирование распространения загрязнения в окружающей среде / М.З. Згуровский, В.В. Скопецкий, В.К. Хрущ, Н.Н. Беляев. - К.: Наукова думка, 1997. - 368 с.

3. Марчук Г.И. Математическое моделирование в проблеме окружающей среды. - М.: Наука, 1982. - 320 с.

4. Самарский А.А. Теория разностных схем / 2-е изд., испр. - М.: Наука, 1983. - 616 с.

5. Эльтерман В.М. Вентиляция химических производств / 3-е изд., перераб. - М.: Химия, 1980.- 288 с.

6. ОНД-86. Методика расчета концентраций в атмосферном воздухе вредных веществ, содержащихся в выбросах предприятий.-Л: Гидрометеоиздат, 1987. - 94 с.

Размещено на Allbest.ru