CC BY
101
УДК 551.510.42
ПРОГРАММА РАСЧЕТА ПОЛЯ ПРИЗЕМНОЙ КОНЦЕНТРАЦИИ ВРЕДНЫХ ВЕЩЕСТВ В АТМОСФЕРЕ
БОРЗЕНКОВ Б.И., МАЛЮК В.Г._____________
Описываются методика, программная реализация и результаты практических расчетов поля приземной концентрации газовоздушных выбросов техногенного характера в атмосферу.
Освоение природных ресурсов, строительство и эксплуатация промышленных комплексов, возникновение и рост городов сопровождается выбросами в атмосферу вредных примесей, способных переноситься на значительные расстояния. Антропогенное загрязнение атмосферы примесями (в том числе химически и биологически активными) связано с деятельностью мощных промышленных предприятий, разнообразием выпускаемых ими материалов, несовершенством технологических процессов. Выбросы загрязняющих веществ могут носить в том числе и аварийный характер [1,2]. Загрязнение приземного слоя атмосферы в значительной степени зависит от метеоусловий. Инвентаризация выбросов состоит в определении их характеристик и систематизации сведений о распределении источников на местности. По результатам мониторинга производится прогнозирование возможного вредного влияния загрязнения, что выполняется в целях экспертизы проектов промышленных сооружений, определения границ санитарно-защитных зон и зон загрязнения в случаях аварийных выбросов. Такое прогнозирование невозможно без использования математических моделей.
1. Методика расчета
Существует нормативный документ, регламентирующий методику расчета приземных концентраций вредных веществ, предназначенный для обязательного использования при проектировании и строительстве предприятий и жилых зданий, нормировании вредных выбросов в атмосферу, экспертизе и согласованию экологоохранных мероприятий [3 ]. Методика, используемая в нем, основана на расчете наибольшей величины приземной концентрации, которая может устанавливаться на определенном расстоянии от места выброса. Учет влияния группы источников на результирующее загрязнение производится на основе расчета процессов рассеивания выбросов от одиночного точечного источника, для которого величина максимальной приземной концентрации вредных веществ См, мг/ м3 определяется по формуле:
Cm = A • M • F • m • и/(и2 • ^Vi -Дт),
где А—коэффициент, зависящий от температурной стратификации атмосферы, которая определяет условия вертикального и горизонтального рассеивания вредных веществ в атмосферном воздухе; M — количество вредного вещества, выбрасываемого в
атмосферу; F — безразмерный коэффициент, учитывающий скорость оседания вредных веществ в атмосферном воздухе; m и n — безразмерные коэффициенты, учитывающие условия выхода газовоздушной смеси из устья источника выброса; H —высота источника выброса над уровнем земли; W0 — средняя скорость выхода газовоздушной смеси из устья источника выброса; ДТ — разность между температурой выбрасываемой газовоздушной смеси и температурой окружающего атмосферного воздуха; V1 — обьем газовоздушной смеси, выбрасываемой в атмосферу, который определяется по
формуле V1 =
п- D2W(
0
где D — диаметр устья
источника выброса.
Величина безразмерного коэффициента m определяется по формуле
m =
1
0.67 + 0.1# + 0.34f
3 W0 • D
т2
где f = 10
H -ДТ
Величина безразмерного коэффициента n опреде-
\удт и
ляется в зависимости от параметра Vm = 0.653,
следующим образом: при Vm<0.5 — n=4.4Vm; при 0.5 < Vm < 2 - n = 3.13 - 2.13Vm + 0.532V2 ; при Vm>2 - n=l.
Скорости ветра вводятся в расчет, как правило, по направлениям розы ветров.
Для определения составляющей скорости ветра U в направлении от источника выброса к точке, в которой производит- Ui
ся вычисление теку- Рис. 1. К определению щей концентрации, скорости ветра U, направлен-мы использовали ли- ного к расчетной точке
нейную аппроксимацию (рис. 1).
Полагая скорости U1 и U2 на краях сектора розы ветров известными, можно определить составляющую скорости ветра в направлении на точку:
U =
________U_________
cos a-sin а(1 - U1V2/U2) .
Величина опасной скорости ветра Um, при которой имеет место наибольшая приземная концентрация вредных веществ в атмосферном воздухе, определяется в зависимости от параметра Vm следующим образом: при Vm<0.5 — Um =0.5; при 0.5 < Vm < 2 - Um =Vm; при Vm > 2 — Um = Vm(1 + 0.12#).
Величина максимальной приземной концентрации Сми с учетом показателей местной средней скорости ветра U, отличающейся от опасной скорости ветра Um, определяется по формуле
C
ми
r • C
m
(1)
152
РИ, 2000, № 2
где r — безразмерная величина, определяемая в зависимости от соотношения q=U/Um: r = 0.67q + \.61cq - 1.34q3 при q < 1; r = 3q/(2q2 -q + 2) при q >1.
Расстояние Xm от источника, на котором достигается максимальная приземная концентрация, определяется по формуле: Xm = d • H при f < 2 ; Xm = (5 -f)• d• H/4 при f > 2, где d — безразмерная величина, определяемая в зависимости от параметра Vm: при Vm<0.5 — d = 2.4б(і + 0.28^f); при
0.5 < Vm < 2 - d = 4.95 • Vm(l + 0.28^6); при Vm > 2
d = VVm(l + 0.28^f).
Расстояние Хми, на котором достигается максимальная предельная концентрация Сми, с учетом местной средней скорости ветра Um определяется по формуле:
Xми = Р • Xm , (2)
здесь p — безразмерная величина, определяемая в зависимости от соотношения q=U/Um: при q < 0.25 — p=3; при 0.25 < q < 1 — p = 8.4^1 - q)5 +1; при q > 1 — p = 0.32q + 0.68 .
Текущая концентрация С вещества на расстоянии X от источника определяется по эмпирическому нормированному графику одномерного распределения иммиссии с подветренной стороны источника газовоздушных выбросов. Этот график представляет собой зависимость величины f=C/C™ от соотношения х=Х/Хми [3]. Там же приводятся формулы аппроксимации данного графика. Нами была использована более простая и более точная формула
f(x) = a1 |a3 - exp (- xa2 |exp(- a4x), (3)
где a!, a2, a3, a4 — параметры аппроксимации. Результаты тестирования данной формулы изложе -ны в работе [4].
Рассматривая концентрацию примеси от произвольной системы источников как суперпозицию вкладов от отдельных источников, можно рассчитать поле концентраций при текущих метеопараметрах с учетом розы ветров.
2. Программная реализация
Описанная выше методика используется в программе “Pole”, предназначенной для экологических расчетов, анализа экстремальных ситуаций, а также для использования в системах оперативного мониторинга воздушного бассейна в промышленных районах [5,6]. Все расчеты проводятся в привязке к крупномасштабной карте местности, которая вводится в компьютер в виде bmp-файла. Исходные данные для расчета хранятся в файле проекта, включающем информацию о расположении и характеристиках источников выбросов, а также текущих значениях метеопараметров и скорости ветра: коэффициент стратификации; температура атмосферного воздуха, оС; скорость ветра по направлениям розы ветров, м/с.
Для каждого источника выбросов указываются: номер источника; имя источника; название вредного вещества; координаты источника; поток вредного вещества, г/с; коэффициент оседания вещества;
скорость выхода вещества, м/с; диаметр источника выброса, м; высота источника выброса, м; температура выбрасываемой смеси, оС.
Результаты расчета двумерного поля концентраций представляются в графическом виде в привязке к карте местности. Предусмотрено графическое выделение зон, в пределах которых рассчитываемые значения превосходят предельно допустимые концентрации конкретного вредного вещества. Допускается изменение масштаба и скроллинг изображения.
Результаты расчета поля концентраций могут быть представлены в трехмерном изображении. Допускается вращение изображения вокруг любой оси.
На рис. 2,3 приведены диаграммы распределения примесей (вид сверху).
Рис. 2. Распределение Рис.3.Распределение
примесей от одного Примесей от трех
источника источников
Исходный код программы занимает около 0,05 Мб, объем загрузочного модуля составляет 0,105 Мб.
Литература: 1. Внуков А.К. Защита атмосферы от выбросов энергообъектов. М.:Энергоатомиздат, 1992. 176 с. 2. Павленко Ю. С., НиконовА.М. Прогнозирование распространения загрязнений в приземном слое атмосферы/ /Гигиена и санитария. 1993. №6. С.67-69. 3. Методика расчета концентраций в атмосферном воздухе вредных веществ, содержащихся в выбросах предприятий (ОНД-86). Л.: Гидрометеоиздат, 1987. 93 с. 4. Борзенков Б.И., Малюк В.Г., Леонидов В.И., Ключка Р.В. Аппроксимация функции распределения приземной концентрации при-меси//Труды междунар. научн.- техн. конф. “Вопросы проектирования, эксплуатации технических систем в металлургии, машиностроении, строительстве”, Ст. Оскол, 1999. С. 12-13. 5. Борзенков Б.И., Малюк В.Г., Леонидов В.И., Ключка Р.В. Прогнозирование распространения загрязнений в атмосфере//Тез. докл. на междунар. научн.- техн. конф. “Вопросы проектирования, эксплуатации технических систем в металлургии, машиностроении, строительстве”, Ст. Оскол, 1999. С.25.
6. Борзенков Б.И., Малюк В.Г., Леонидов В.И. Система оперативного контроля экологии воздушного бассейна промышленного района//Тез. докл. на 4-й междунар. конф. “Теория и техника передачи, приема и обработки информации”, Харьков, 1998. С.80.
Поступила в редколлегию 20.11.99
Рецензент: д-р техн. наук, проф. Прохач Э.Е.
Борзенков Борис Иванович, канд. техн. наук, доцент кафедры биомедицинских электронных устройств и систем ХТУРЭ. Научные интересы: разработка электронных систем, моделирование физических процессов. Хобби: шахматы. Адрес: Украина, 61166, Харьков, пр. Ленина, 14, тел. 40-93-64, 37-25-74.
Малюк Виктор Григорьевич, канд. техн. наук, доцент кафедры биомедицинских электронных устройств и систем ХТУРЭ. Научные интересы: моделирование процессов в биосистемах, компьютерные технологии в медицине. Хобби: художественная литература, водный туризм. Адрес: Украина, 61166, Харьков, пр. Ленина, 14, тел. 40-93-64, 94-11-71.
РИ, 2000, № 2
153
