Scientific journal

Fundamental research

ISSN 1812-7339

"Перечень" ВАК

ИФ РИНЦ = 1,674

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

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

Приведем краткое изложение метода идентификации, описанного в [5, 6]. Входной информацией являются прямые измерения приземных концентраций загрязняющих веществ на территории города вне санитарно-защитных зон предприятий и оперативная информация о температуре воздуха, скорости и направлении ветра. Основными выходными данными являются характеристики выбросов единичного источника в атмосферу, значения которых могут изменяться в реальном масштабе времени (согласно [3, 5, 6]): Mi[t] – масса выбросов загрязняющего вещества в единицу времени от i-го источника, г/с и Vi[t] – расход газовоздушной смеси из i-го источника, м3/с, где i – индекс источника выброса загрязняющего вещества; I [t] – множество источников, оказывающих влияние на значения концентраций в некоторой точке; t – момент текущего времени; g∈{1,…,G} – индекс загрязняющего вещества, G – общее число веществ, по которым осуществляется мониторинг (далее изложение ведется применительно к одному из веществ, поэтому индекс вещества g не будет указываться).

Остальные параметры стационарных источников выбросов в атмосферу (в частности, высота трубы над поверхностью земли Hi, м; диаметр устья трубы Di, м; координаты размещения i-го источника (xi, yi), м) задаются для всех i∈I[t] и рассматриваются в задаче идентификации как константы.

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

Метод идентификации был разработан в работах [5, 6] и основан на поиске минимума среднеквадратического критерия рассогласования между вычисленными и измеренными концентрациями загрязняющих веществ

где Cf ( j) – измеренная концентрация загрязняющего вещества в точке j, c координатами xj, yj; Cm(xj, yj, M, V) – вычисленные значения по модели ОНД-86 [], M, V – вектора масс и объемов выбросов, элементами которых являются массы и объемы выбросов отдельных источников; N – количество точек измерения.

Поиск минимума основан на простейшем градиентном методе наискорейшего спуска [1]

(1)

где является нормированным градиентом критерия по искомым переменным, рассчитанным для точек (xj, yj), j = 1,…, N; λ(v) – шаг алгоритма, v – итерация алгоритма.

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

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

Эксперименты проводились для оценки параметров выбросов диоксида серы из труб пяти различных предприятий г. Камышина Волгоградской обл. для следующих значений метеорологических условий и экологических параметров (такие же, как предложенные в [4–6]): экологические константы: A = 200, F = 1, η = 1; зимнее время, ∆T = 135 °C; угол между направлением ветра и абсциссой равен 0,25π; абсцисса направлена на север; угол отсчитывается против часовой стрелки; шаги сетки на зоне рассеивания – 100 м по осям абсцисс и ординат; опасная скорость ветра um = 4,46 м/с.

Параметры источников выбросов приведены в таблице.

Пример результата вычислительных экспериментов первой группы представлен на рис. 1.

Рис. 1. Смещенность в оценке объемов выбросов источников № 1 (точное значение 165 м3/с)

По результатам экспериментов этой группы можно сделать вывод, что, с одной стороны, отсутствие одной точки измерения (четыре точки на пять источников) приводит к смещенности в оценке параметров источников (рис. 1). А с другой стороны, отсутствие измерения одного источника может не привести к смещенности оценок объемов выбросов ни одного источника. Это объясняется тем, что подфакельные оси некоторых источников совпадают между собой и с направлением ветра или точки измерения находятся достаточно близко (расстояние до 200 м).

Пример результатов экспериментов второй группы приведен на рис. 2.

Рис. 2. Ошибка идентификации объема выбросов источника № 1 в зависимости от отклонения концентраций в точке измерения

Представляет определенный интерес поведение алгоритма при больших величинах отклонений концентрации в точке измерения. Так, было установлено, что при превышении концентрации в 3 раза ошибка алгоритма идентификации по определению объемов выбросов первого источника составит порядка 73 %.

Результаты экспериментов этой группы были получены и на других источниках выбросов.

Таким образом, из проведенных экспериментов можно сделать следующие выводы:

1. При определенных условиях расположения точек контроля, источников выбросов и направления ветра можно уменьшить количество исходных данных для алгоритма идентификации.

2. При наличии неучтенного источника выбросов ошибка идентифицируемых параметров линейно зависит от мощности источника и появляется смещенность в оценке параметров стационарных источников.

Работа выполнена при финансовой поддержке РФФИ, грант 14-07-97011.

Фоменков С.А., д.т.н., профессор, профессор кафедры «Системы автоматизированного проектирования и поискового конструирования» Волгоградского государственного технического университета, г. Волгоград;

Томашевский Ю.Б., д.т.н., профессор, заведующий кафедрой «Системотехника» Института электронной техники и машиностроения Саратовского государственного технического университета имени Гагарина Ю.А., г. Саратов.

Работа поступила в редакцию 05.12.2014.