Разработка программного модуля для нахождения оптимальных предельно-допустимых выбросов в атмосф (стр. 7 из 10)

Санитарные требования приводят к системе линейных неравенств:

решение которой ищется на интервале

где

- технологически обоснованный минимум выброса j-го источника.

На сегодняшний день в методической литературе описаны два метода нахождения решения поставленной задачи: МРН-87 [24] и метод равного квотирования [26]. Оба метода кратко рассмотрены в параграфе . Они дают частное решение поставленной системы неравенств из соображений удобства нахождения единственного решения. Однако любое предприятие заинтересовано в минимальных затратах, необходимых для установления нормативных выбросов

. Для этой цели, к поставленной системе неравенств добавляем целевую функцию:

Где в общем случае

- стоимость снижения на единицу выброса для j-го источника. В данном виде решение X_j дает минимум затрат на достижение нормативного загрязнения атмосферы.

В предположении

что эквивалентно

Следовательно, задача сводится к поиску максимально допустимого по сумме сочетания выбросов X_j данного предприятия, позволяющего достичь нормативного загрязнения атмосферы. В связи с линейностью модели ОНД-86 по отношению к выбросам для поиска решения поставленной задачи может быть применён Симплекс-метод.

4. Программная реализация и пример практического применения

Для достижения поставленной задачи по разработке интегрированного в ПК ЭРА-ВОЗДУХ программного модуля расчета оптимальных ПДВ проведено изучение структуры файлов, в которых головной модуль передают данные расчетного блоку УПРЗА ЭРА. Разработаны соответствующие процедуры для автоматизированного чтения всех необходимых файлов.

4.1 Выбор загрязняющих веществ

После указания директории с данными для расчета, программа сканирует файлы в папке WORK по маске «htop*.ppp», таким образом выбирая загрязняющие вещества, для которых имеются начальные данные и возможно провести расчет.

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

4.2 Обработка точек с повышенным загрязнением

procedure get_point (s:string;var countPoint:integer;var point_pdk:tExtArray);

Процедура извлекает значения N_i (ПДК), а так же количество точек I (i=1,..,I), в которых выбросы превышают N_i, из файла вида «htop*.ppp» в директории /WORK/.

Внутренняя структура htop*.ppp представляет собой текстовый файл, содержащий таблицу с данными о контрольных точках. Значение N_i содержатся в 8-ом столбце. Чаще всего N_i в точках равно 1 или 0.8 для особо охраняемых территорий (санатории, зоны отдыха).

4.3 Обработка источников

procedure get_funnel(s:string; var countFunnel : integer ; var funnel_name : tsArray; var funnel_m:tExtArray;var funnel_min:tExtArray);

Процедура извлекает данные из файла вида «ist_*.txt», которые находятся в директории /DAT/. «ist_*.txt» - это текстовый файл, в котором в табличном виде представлена информация об источниках выбросов, в том числе:

· количество источников J;

· уникальный код источника;

· существующий выброс источника Qj;

· минимально возможный выброс Qjmin (не всегда указывается).

4.4 Обработка таблиц влияния источников на точки

procedure get_pointfunnel ( s : string; countPoint : integer; countfunnel : integer; funnel_name : tsArray; funnel_m : tExtArray; var pointfunnelx2 : tExtArrayx2; var point_cf : tExtArray);

Процедура извлекает коэффициенты влияния a_ij из файлов вида «10pd*.ppp», где «10pd*.ppp» - текстовый файл, содержащий отчет о результатах работы программы.

4.5 Применении симплекс-метода

procedure get_simplexsolve ( countPoint : integer; countFunnel : integer ; point_pdk : tExtArray; point_cf : tExtArray; funnel_m : tExtArray; funnel_min : tExtArray; pointfunnelx2 : tExtArrayx2; var x : tExtArray; var s_temp : string);

Процедура, используя данные расчетов программы «ЭРА-воздух», при помощи симплекс-метода рассчитывает оптимальные выбросы Xj для источников при заданных условиях.

4.6 Вывод полученных результатов

Результаты полученных вычислений выводятся на форму программы, а так же в файлы вида: «h_pd*.gpv», где * - это код вещества, для которого производился расчет.

4.7 Сравнение различных методов расчета ПДВ для реального предприятия

В расчете загрязнения атмосферы диоксидом азота на расчетном прямоугольнике 9 на 9 км с шагом 500 метров проведен расчет максимальных разовых концентраций, создаваемых 70-ю источниками выброса различного типа и высоты. В итоге получено поле максимальных концентраций, в которых есть области превышения норматива N, который в данном случае равен ПДК.

Картина загрязнения представлена на рисунке 1:

Рис 1. Загрязнение атмосферы в окрестности исследуемого предприятия.

Программный комплекс Эра-Воздух производит автоматический выбор точек превышения норматива N и позволяет найти расчетные значения ПДВ на основе методов МРН-87 и Метода равного квотирования.

При выполнении работы в программный комплекс добавлен новый модуль, в соответствии с алгоритмом описанным в пункте 5. Этот модуль позволяет рассчитать оптимальное значение ПДВ с использованием симплексного метода. Результаты расчетов ПДВ различными методами представлено ниже в таблицах 4.1, 4.2 и 4.3.

Таблица 4.1. Результаты расчета ПДВ (методом равного квотирования)

ПРИМЕСЬ=0301 Азот(IV) оксид (Азота диоксид)

Город :001 Кемерово, Объект : 0025

-----------------------------------------------------------------------------

п/п| выброса | м | г/с | выброс |рования | П Д В |выброса |

---|-----------|-------|----------|---г/с----|--------|----г/с----|---------|

15 00250010354 40.0 8.2390 0.0 0.485 3.9961 2.062

53 00250010944 40.0 8.2370 0.0 0.461 3.7941 2.171

Остальные источники не подлежат нормированию.

-----------------------------------------------------------------------------

В сумме по 0301 51.8545 0.0 0.832 43.1687 1.201

-----------------------------------------------------------------------------

Таблица 4.2. Результаты расчета ПДВ (метод МРН-87)

ПРИМЕСЬ=0301 Азот(IV) оксид (Азота диоксид)

Город :001 Кемерово, Объект :0025

-----------------------------------------------------------------------------

п/п| выброса | м | г/с | выброс |рования | П Д В |выброса |

---|-----------|-------|----------|---г/с----|--------|----г/с----|---------|

14 00250010353 33.4 1.0750 0.0 0.739 0.7942 1.354

15 00250010354 40.0 8.2390 0.0 0.739 6.0867 1.354

16 00250010356 60.0 2.1500 0.0 0.739 1.5884 1.354

24 00250010656 45.0 1.5000 0.0 0.739 1.1082 1.354

41 00250010892 60.0 9.0430 0.0 0.739 6.6807 1.354

52 00250010943 33.4 0.5280 0.0 0.739 0.3901 1.354

53 00250010944 40.0 8.2370 0.0 0.739 6.0852 1.354

54 00250010946 45.0 2.1470 0.0 0.739 1.5861 1.354

63 00250011198 24.4 0.5110 0.0 0.739 0.3775 1.354

66 00250011225 24.4 0.5890 0.0 0.739 0.4351 1.354

Остальные источники не подлежат нормированию.

----------------------------------------------------------------------------

В сумме по 0301 51.8545 0.0 0.829 42.9676 1.207

----------------------------------------------------------------------------

Таблица 4.3. Результаты расчета ПДВ (симплекс метод)

ПРИМЕСЬ=0301 Азот(IV) оксид (Азота диоксид)

Город :001 Кемерово, Объект :0025

---------------------------------------------------------

|-----------|----------|---г/с----|----г/с----|---------|

|00250010353| 1.075000 | 0.000000 | 0.0196222 | 0.01825 |

|00250010943| 0.528000 | 0.000000 | 0.0000000 | 0.00000 |

|00250011225| 0.589000 | 0.000000 | 0.0000000 | 0.00000 |

- - - - - - - - - - - - - - - - - - - - - - - - - - - - -

| в сумме: 51.854470 49.682092 | 0.95811 |

---------------------------------------------------------

Решение задачи линейного программирования показывает, что максимально возможный выброс по заводу в целом можно обеспечить при достижении нормы загрязнения, если закрыть источники 1225 и 0943, а на источнике 0353 снизить выброс примерно в 50 раз. О реальности такого решения могут судить технологические службы предприятия.

Итоговая таблица сравнения трех методов:

Метод	Существующий выброс	Расчетное значение	Процент снижения
Равного квотирования	51.8545	43.1687	16,75
МРН-87	51.8545	42.9676	17,14
Симплексный	51.8545	49.6820	4,19

Заключение

· Таким образом, в процессе выполнения дипломной работы проведен обзор существующей системы установления ПДВ для источников загрязнения атмосферы промышленных предприятий. При этом рассмотрена система обеспечения нормативных задач управления выбросами в атмосферу исходными данными. Исследована структура исходных данных их качество. Рассмотрены нормативные требования к загрязнению атмосферы населенных мест, которые накладывают ограничения на выбросы промышленных предприятий в атмосферу. Показано, что эти требования и свойство линейности нормативной модели расчета загрязнения атмосферы ОНД-86 по отношению к выбросам ИЗА позволяют представить процедуру нахождения расчетных значений ПДВ в виде задачи линейного программирования.