Выбор управлений, наиболее эффективных с точки зрения «природоохранных» и «производственных» критериев, является непростой задачей, решение которой вряд ли возможно без применения метода математического моделирования. В данной работе рассматривается математическая модель основанная на численном решении управления переноса и диффузии загрязняющих примесей.
Проблемы охраны и управление качеством окружающей среды порождает класс задач, связанных с поиском оптимальных решений при подготовке народно-хозяйственных проектов, осуществление которых сопряжено с воздействием на природную среду, а также при планировании природоохранных мероприятий, требующих управления выбросами действующих промышленных объектов с учетом особенностей гидрометеорологического режима и ограничений санитарно и социально-экономического характера.
В связи с этим в практике хозяйствования все значение приобретают методы улучшения качества окружающей среды. К этим методам можно отнести:
– реконструкция и усовершенствование действующих технологических процессов, обеспечивающих снижение выбросов примесей и вредных отходов;
– разработка и внедрение малоотходных (замкнутых) технологических процессов, обеспечивающих комплексное использование всех компонентов и минимальное поступление выбросов в окружающую среду.
В настоящее время имеются несколько типов моделей, отражающих те или иные аспекты взаимодействия общества и среды с учетом загрязнения окружающей среды и его социально-экономических последствий.
В работах [1–5] сформулирован ряд математических моделей для решения такого рода задач. В данной работе в качестве целевой функции выступает функционал стоимости ущерба от отдельных источников и затрат на их оптимизацию. Эти функции зависят от концентрации примесей и могут зависеть от входных параметров модели.
Пусть рассматриваемый регион расположен в ограниченной трехмерной области и на его территории имеется nпромышленных предприятий, производящих выбросы вредных веществ в атмосферу.
Не ограничивая общности, будем считать источники выбросов точечными и расположенными внутри области. Для описания процесса распространения примеси от указанных источников воспользуемся линейным уравнением турбулентной диффузии (4, 5)
(1)
С краевым и начальными условиями
(2)
(3)
(4)
(5)
Где – концентрация примеси в точке
в момент времени t;
– коэффициенты горизонтальной и вертикальной диффузии;
– функция, описывающая скорость изменения примеси вследствие химических превращении;
β – коэффициент, характеризующий взаимодействие примеси с подстилающей поверхностью;
– значение концентрации примеси S – боковой поверхности области D и в начальный момент времени;
– известная функция, описывающая неорганизованные источники выбросов;
Qm – начальная мощность выброса источника, m=1, .
Будем считать заданными начальные и краевые условия (2)-(5), скорость ветра , коэффициенты турбулентного обмена, и ограничимся рассмотрением модели «интегральной» концентрации [5, 6, 7]. Предположим, кроме того, что коэффициенты (1)-(5) и оператор δ не зависят от искомого решения.
Введем следующие обозначения: – вектор размерности компонента, которого описывает относительное уменьшение интенсивности выброса от i-го источника,
Ei – максимум относительного уменьшения интенсивности выброса от i-го источника, i=1,
.
В силу линейности модели (1)-(5), концентрация прмеси в точке после уменьшения интенсивности выбросов определяется по формуле [4, 7]:
(6)
Где – решение набора задач
при z = 0, H (7)
m=1, n
а функция удовлетворяет задаче
при z = 0, H
(8)
При выработке долговременной политики улучшения качества атмосферы используются, как правило, значения осредненных концентраций примеси приземного слоя [5, 7].
Интервал осреднения T может быть различной величины в зависимости от конкретных обстоятельств, например: год, сезон.
Обозначим осредненные концентрации примеси через и определим их как
(9)
Где
Решения задач (6)-(8) находим с помощью метода конечных элементов [4, 7].
Пусть для численного интегрирования этих задач в D введена сеточная область
Где – сетка в направлении координат x,y,z с шагом
, i=0, I,
;
Поскольку нас будут интересовать приземные концентрации, то введем обозначения
(10)
где
Таким образом, множество векторов описывает возможные последствия загрязнения в точках приземного слоя в зависимости от варианта управления
.
Мы рассмотрим модель основанный на понятие функции стоимости регулирования источников [3, 4, 5].
Введем, следуя [8], следующие обозначения: Gm(lm) – функция, характеризующая стоимость уменьшения интенсивности выбросов на m – M предприятий на величину lm, m=1, n; общая стоимость регулирования источников в пределах данного региона.
Пусть S – стоимость всех средств, используемых для улучшения качества атмосферы. Тогда множество E можно считать заданным в виде
(11)
Для построения зависимостей Gm(lm), m=1, n, имеется несколько путей. Например, когда уменьшение выбросов производится за счет предварительной очистки исходного сырья или топлива, то стоимость регулирования может быть определена как функция от веса элементов, отдельных в результате очистки. Гораздо сложнее обстоит дело, когда для уменьшения выбросов используются такие методы, как модернизация и реконструкция существующего производства. В этом случае затраты на эти мероприятия могут дать одновременно и положительный производственный эффект. Даже такая сугубо атмосфероохранная мера, как повышение высоты трубы, увеличивает интенсивность процесса горения. Что повышает эффективность производства за счет более полной утилизации сырья и топлива [1, 6].
Следуя [8], будем понимать под стоимостью противозагрязняющих мероприятий на предприятии ((Gm(lm))) сумму всех издержек, которые несет данное предприятие, при уменьшении объема вредных выбросов на величину em, m=1, n и неизменном объеме выпускаемой продукции. Основные статьи этих издержек связаны с закупкой других, более дорогих видов сырья и материалов, дополнительными капиталовложениями эксплуатационными затратами на основания новых малоотходных технологий, увеличением себестоимости производимой продукции и, следовательно, уменьшением прибыли от ее реализации.
Для расчета этих составляющих элементов стоимости предотвращения загрязнения может быть успешно применен метод экономико-математического моделирования [3, 4].
Рассмотрим далее метод построения функций Gm(lm), основанный на использовании одной из простейших и наиболее употребительных моделей линейной производственной модели [1, 9]. Будем исходить прежде всего из того, что рассматриваемой производственной единице (m-му предприятию) установлено плановое задание по выпуску Jm видов продукции в объемах за время Т. Для выполнения этого задания предприятие распологает технологическими способами .
Обозначим через hm1 интенсивность использования 1-ой технологии на предприятии, – вектор интенсивного функционирования m – го предприятия. Представим экономико-математическую модель работы данного предприятия следующим образом:
(прибыль) (12)
При ограничениях
(плановое задание) (13)
(предельно допустимые выбросы) (14)
(себестоимость) (15)
(капиталовложения) (16)
(экспулатационные затраты) (17)
Здесь используются следующие обозначения:
– прибыль m-го предприятия при использовании е-й технологии с единичной мощностью;
– объем выпуска продукции вида j на m-м предприятии по способу производства;
– мощность выброса примеси на m-м предприятии по е-й технологии с единичной мощностью;
– себестоимость производства единицы продукции е-го вида для m-го предприятия по технологическому способу;
– капитальное вложение на предотвращения загрязнения атмосферы в е-ю технологию на m-м предприятии (
для вновь осваиваемых и реконструируемых технологий);
– затраты на эксплуатацию и содержание газоочистных установок и других очистных сооружении в связи с применением е-го технологического способа на m-м предприятии;
– лимит капиталовложений для е-го предприятия;
– лимит эксплуатационных затрат на m-м предприятии;
– предельная величина себестоимости выпускаемой m-м предприятием продукции.
Оптимальный план функционирования m-го предприятия определяется с помощью решения задачи линейного программирования (12)-(17). Переменная em участвует в этой задаче как параметр. Величина Em в рамках принятой модели может быть определена из решения задачи линейного программирования:
При ограничениях (12)-(17). Если – оптимальное решение этой задачи, то
.
Пусть – вектор оптимального решения задачи оптимизации (12)-(17), зависящий от параметра. Используя введенные обозначения вычислим следующие величины:
– прибыль m-го предприятия при оптимальном плане
– себестоимость продукции на m-м предприятии при оптимальном плане
– объем капталовложений на оздоровление атмосферы, требуемый m-му предприятию при данном режиме работы
– стоимость всех эксплуатационных затрат, связанных с работой очистных сооружений и установок при плане работы m-го предприятия.
Определим суммарные издержки m-го предприятия (функцию Gm(lm)), возникающие вследствие уменьшения выбросов на lm. Тогда
(19)
Отметим, Gm(lm) что функция lm является кусочно-линейной функцией на [0, Em ]. Это следует из общего свойства решения задач линейного программирования, вектор-функция «склеена» из кусков линейных отрезков в Rm. Поэтому для построения функции стоимости Gm(lm) достаточно иметь решения нескольких задач линейного программирования вида (12)-(17),которые последовательно принимают значения, равные точкам излома графика вектор-функции
.
Таким образом, построенная экономико-математическая модель управления используется для описания процессов распространения загрязнителей в численных моделях. Это позволяет получить оценку уровней загрязнения в точках рассматриваемого региона, которые далее могут быть использованы для формирования критерия качества воздушного бассейна области. Целевая функция представлена в виде свертки кусочно-линейной функции.