Разработка программного модуля для нахождения оптимальных предельно-допустимых выбросов в атмосферу от группы источников
1. Существующая система установления ПДВ для промышленных источников
1.1 Нормативы и показатели загрязнения атмосферы.
1.2 Данные об источниках загрязнения атмосферы.
1.3 Метеопараметры
1.4 Данные наблюдений за загрязнением атмосферы
1.5 Модели расчета загрязнения атмосферы
1.6 Унифицированные программы расчета загрязнения атмосферы
1.7 МетодЫ равного квотирования и МРН-87
2. Симплекс-метод
2.1 Общая характеристика симплекс метода
2.2 Алгоритм симплекс метода (первая симплекс таблица)
3. Формализация поставленной задачи
4. Программная реализация и пример практического применения
4.1 Выбор загрязняющих веществ
4.2 Обработка точек с повышенным загрязнением
4.3 Обработка источников
4.4 Обработка таблиц влияния источников на точки
4.5 Применении симплекс-метода
4.6 Вывод полученных результатов
4.7 Сравнение различных методов расчета ПДВ для реального предприятия
Заключение
Список литературы
Вспомогательные указатели
Перечень сокращений
Приложения
Unit1.pas
Simplex.pas
ВВЕДЕНИЕ
В нашей стране существует государственная система управления выбросами промышленных источников загрязнения атмосферы (ИЗА) (1-4). Целью системы является достижение и поддержание на нормативном уровне N загрязнения воздушного бассейна населенных мест, которое характеризуется приземной разовой (осредненной за 20 минут) концентрацией С.
Основными разделами данной системы является проведение воздухоохранной экспертизы вновь строящихся объектов и установление предельно допустимых норм (ПДВ) для существующих объектов.
В этих задачах расчетным путем устанавливается связь между выбросами Q предприятий и приземными концентрациями C загрязняющих веществ (ЗВ), выбрасываемыми промышленными источниками.
Для расчета концентраций С по известным значениям Q и характерным для данной территории метеопараметрам используется единая, утвержденная на государственном уровне методика ОНД-86 (2).
В настоящее время проектные расчеты осуществляются с использованием программных комплексов (ПК) сертифицированных главной геофизической обсерваторией им. А.И. Войекова (г. Санкт-Петербург).
При проведении атмосферноохранной экспертизы проекта нового предприятия, решается прямая задача моделирования загрязнения атмосферы.
При установлении ПДВ для существующего предприятия, решается обратная задача – подбираются такие значения выбросов Х, при которых в прилегающих к предприятию жилых районах приземные концентрации С не превышают санитарную норму N.
В соответствии с нормативными документами по нахождения ПДВ (1,4) в этом случае обратная задача решается путем подбора реальных (технологически применимых) для данного производства атмосфероохранных мероприятий (АОМ), таких как замена оборудования, установка газоочистки и т.д. Фактически, обратная задача решается при этом методом подбора. Однако при большом числе источников большую помощь в выборе первоочередных мероприятий дает решение задачи форматизированными методами. Пользуясь линейностью модели ОНД-86 по выбросам источников можно представить загрязнение атмосферы в контрольных точках жилой зоны в виде линейной формы:
,
где 
- концентрация в i-ой точке, 
- выброс j-го источника, 
- вклад j-го источника в i-й точке, который в дальнейшем будем называть коэффициентом влияния.
Санитарные требования приводят к системе линейных неравенств:
,
решение которой ищется на интервале
,
где 
- технологически обоснованный минимум выброса j-го источника.
На сегодняшний день в методической литературе описаны два метода нахождения решения поставленной задачи: МРН-87 (24) и метод равного квотирования (25). Оба метода дают частное решение поставленной системы неравенств из соображений удобства нахождения единственного решения. Однако любое предприятие заинтересовано в минимальных затратах, необходимых для установления нормативных выбросов 
. Для этой цели, к поставленной системе неравенств добавляем целевую функцию:
,
где в общем случае 
- стоимость снижения на единицу выброса для j-го источника. В данном виде решение Xj дает минимум затрат на достижение нормативного загрязнения атмосферы.
В предположении 
получаем
,
что эквивалентно 
.
Следовательно, задача сводится к поиску максимально допустимого по сумме сочетания выбросов Xj данного предприятия, позволяющего достичь нормативного загрязнения атмосферы. В связи с линейностью модели ОНД-86 по отношению к выбросам для поиска решения поставленной задачи может быть применён симплекс-метод.
Целью настоящей работы является разработка программного модуля для расчета оптимальных значений ПДВ с использованием симплекс-метода. Данный модуль должен не просто реализовать решение задачи линейного программирования, но и быть интегрирован в среду программного комплекса ЭРА-ВОЗДУХ, который предназначен для решения проектных задач, использует стандартизированные данные об источниках выбросов, метеорологических параметрах распространения примесей и выдает выходную информацию в установленной нормативными документами форме (подробности на сайте www.logos-plus.ru). После разработки и отладки модуля в работе предполагается решить задачу на основе реальных данных и сравнить оптимальное решение с методами, описанными в методической литературе.
1. Существующая система установления ПДВ для промышленных источников
1.1 Нормативы и показатели загрязнения атмосферы
Основными нормативами чистоты атмосферы, лимитирующими загрязнение (следовательно и выбросы ИЗА), являются предельно допустимые концентрации (ПДК) для населенных мест, которые подразделяются на среднесуточные (ПДКс - длительных периодов воздействия) и максимально-разовые (ПДКр - 20-ти минутного осреднения) (6,7,8). Такие ПДК отражают критерии качества воздуха для человека (гигиенические). В документах (1,5) при определении допустимых выбросов в атмосферу предусмотрен учет ПДК для растительности, животного мира, сельхозугодий и т.д. (экологические) с выбором в качестве определяющей наиболее жесткой величины. Однако пока утвержденные какими-либо нормативными документами экологические ПДК отсутствуют (1), а управление выбросами ведется с ориентацией только на гигиенические ПДК.
ПДКс имеет цель предупреждение хронического воздействия ЗВ при длительном вдыхании. Концентрация ЗВ в пределах ПДКс не должна оказывать на человека прямого или косвенного воздействия при неограниченно продолжительном вдыхании. Поэтому и гигиенистами (9) и специалистами по атмосферной диффузии (10) отмечается, что ПДКс, несмотря на название, можно использовать для таких времен осреднения, как год или сезон. Список ПДКсдополняют максимально разовые предельно допустимые концентрации (ПДКр), при установлении которых учитывается рефлекторное воздействие. В случае отсутствия какого-либо ПДК используется уравнение ПДКр=10ПДКс. Для ряда ЗВ, по которым пока не определены нормативы ПДК, расчетным путем установлены ориентировочно безопасные уровни воздействия (ОБУВ), используемые в качестве ПДКр. На сегодняшний день в перечне (7) около 3000 веществ имеющих какой-либо норматив ПДК. В нашей стране в качестве основной размерности ПДК принят мг/м3.
Рассмотрение списка ПДК показывает, что отношение ПДКр/ПДКслежит в интервале от 1 до 20.
Наличие ПДК, регламентирующих как длительное, так и кратковременное воздействие ЗВ на здоровье человека, свидетельствует, что показателями загрязнения атмосферы населенных мест должны являться концентрации, согласованные с ПДКс и ПДКр по времени осреднения. Таким образом, названными показателями будут среднегодовая Сг и разовая С концентрации, периодами осреднения которых являются год и 20 минут соответственно.
В настоящее время, согласно (1), считается, что чистота ВБ населенных мест для некоторого ЗВ соответствует нормативной, если на границах санитарно-защитных зон предприятий и в пределах жилой зоны выполняется неравенство
См < mПДКр (1.1)
где См– максимальная суммарная концентрация данного ЗВ; m = 0.8 для особо охраняемых территорий и m = 1.0 в остальных случаях. К особо охраняемым территориям относятся места массового отдыха населения (подробно см. (16)) и территории размещения лечебно-профилактических учреждений длительного пребывания больных и центров реабилитации. В качестве См берется вычисленная по (2) суммарная от всех источников концентрация, достигающаяся при неблагоприятном сочетании метеорологических условий (МУ) и параметрах выбросов ИЗА данного ЗВ. Однако, как отмечается в (10), “при планировании долгосрочных АОМ часто невозможно или экономически нецелесообразно предусматривать, чтобы абсолютно ни при каких МУ концентрации примесей не превышали ПДКр. Обычно оказывается достаточным учесть часто наблюдаемые неблагоприятные МУ, а на периоды особо опасных аномальных МУ вводить в действие оперативные АОМ”. Поэтому в неравенстве (1.1) под См подразумевается не абсолютный максимум, а разовая концентрация при так называемых “нормальных неблагоприятных метеоусловиях”. Частота появления МУ, способных вызвать большее загрязнение ВБ (называемых “аномальными, неблагоприятными”) при тех же уровнях выбросов зависит от метеорологического режима района расположения предприятий (10)
Отсутствие в нормативных указаниях (1) аналогичного (1.1) неравенства для годового осреднения
Сг < ПДКс (1.2)
объясняется тем, что если (1.1) выполняется, то (1.2) для того же ЗВ будет следовать автоматически. Действительно, отношение ПДКр/ПДКс за редким исключением не превосходит 10, в то врем как отношение См/Сг числом из интервала от 25 до 40 (10). Однако этот интервал относится прежде всего к случаю изолированного или группы близко расположенных ИЗА. В крупном промышленном городе, при наличии большого числа распределенных по территории и различных по типу источников, нижняя граница указанного интервала существенно снижается. Оценки отношения См/Сг на основе статистического анализа результатов стационарных наблюдений в различных городах России рассматриваются в работах Э.Ю.Безуглой (11,12). Фактические соотношения между См/Сг для серистого ангидрида, полученные по данным регулярных наблюдений на 863 стационарных постах, расположенных в городах России показывают, что нормирование выбросов по разовым концентрациям в 30% случаев не обеспечит соблюдение (1.2). Кроме того, следует подчеркнуть, что определение риска от загрязнения атмосферы для здоровья населения основано именно на долгопериодных концентрациях ЗВ. Отсюда следует, что в общем случае следует предусматривать проверку обеих условий (1.1) и (1.2).
С начала 90-х годов в состав нормативных данных по загрязняющим веществам добавились базовые удельные стоимостные показатели (19) по выбросам каждого ЗВ за превышение нормативов предельно допустимых выбросов (ПДВ). Под ПДВ для предприятия понимается такая совокупность выбросов его ИЗА (1,3,4,5), при которой (с учетом фонового загрязнения атмосферы от других предприятий) соблюдается условие (1.1). За выброс в пределах установленного проектом нормативов ПДВ лимита предприятие платит базовый норматив (руб/тонну в год), а за превышение – в пятикратном размере. Кроме того, начисляемая плата умножается на коэффициент экологической ситуации (территориальный показатель, от 1.0 до 1.44) и коэффициент индексации. Для примера можно сказать, что годовая плата бытовой котельной может составлять от нескольких тысяч до нескольких десятков тысяч рублей, а крупного химического или металлургического предприятия достигает нескольких миллионов.
Здесь хотелось бы еще раз подчеркнуть, что при разработке проекта норм ПДВ условие (1.1) проверяется расчетным путем на основе данных о выбросах ИЗА с использованием методики ОНД-86 (2). Следовательно, и плата за выброс зависит от локальной модели расчета загрязнения атмосферы, которая нормативно утверждена на федеральном уровне. Как и любая (а тем более - инженерная) модель столь сложного явления как атмосферная диффузия, она не описывает всех природных закономерностей и может допускать значительные погрешности. Тем не менее данная проблема в нормативных работах не обсуждается, поскольку в разрешении финансовых вопросов обязательно требуется ссылка на законодательно утвержденный документ, которым и является ОНД-86.
1.2 Данные об источниках загрязнения атмосферы
Естественно, важнейшей информацией для управления выбросами промышленных предприятий являются параметры их ИЗА. Существует ориентированная на данные цели система сбора и хранения такой информации на государственном уровне. При постановке многих научных и экспертных задач, не входящих в перечень нормативных проектных проработок, исследователями явно или неявно предполагается, что названная система обеспечивает (или обеспечит в ближайшем будущем) представление информации об ИЗА в необходимом пространственном и временном разрешении и с требуемой для поставленной задачи точностью. В этой связи представляется интересным рассмотреть информацию об ИЗА и те ограничения которые должны учитываться не только при постановке, но и при решения исследовательских задач применительно к реальным объектам промышленных городов или регионов.
Источник загрязнения атмосферы (ИЗА) представляет собой хозяйственное сооружение или естественный объект, из которого в атмосферу поступают ЗВ. К источнику загрязнения, как правило, относится один или несколько источников выделения, которыми называются промышленные установки, генерирующие ЗВ в процессе работы (например, водогрейный котел есть источник выделения, а труба – источник загрязнения).
Нормативная информация о параметрах ИЗА представляется предприятиями в природоохранные органы на основе проведения инвентаризации. В Федеральный Закон "Об охране атмосферного воздуха" (5) впервые (в сравнении с ранее действовавшим Законом) введена статья, касающаяся инвентаризации выбросов ЗВ в атмосферный воздух. В соответствии со статьей 22 Закона “юридические лица, имеющие источники выбросов ЗВ в атмосферный воздух, проводят инвентаризацию выбросов вредных (загрязняющих) веществ в атмосферный воздух и их источников в порядке, определенном специально уполномоченным федеральным органом исполнительной власти в области охраны атмосферного воздуха» (в настоящее время - Ростехнадзор РФ). Инвентаризацию проводят все действующие предприятия, организации, учреждения независимо от их организационно-правовых форм и форм собственности, производственная деятельность которых связана с выбросом ЗВ в атмосферу. Ответственность за полноту и достоверность данных инвентаризации несет предприятие (в лице руководителя). Инвентаризация выбросов (т.е. представление информации с точностью до каждого ИЗА) проводится 1 раз в 5 лет (1, 15). Кроме того, ежегодно предприятия отчитываются по форме государственного статистического наблюдения № 2-тп (воздух), которая содержит сведения только о годовых выбросах предприятия в сумме по всем источникам.
Определение параметров источников загрязнения атмосферы (ИЗА) для проекта нормативов ПДВ должно осуществляться при регламентных загрузке и условиях эксплуатации технологического и пыле-, газоочистного оборудования. В соответствие с новыми требованиями (1), параметры ИЗА следует фиксировать и на основных режимах работы технологического оборудования (установки) и различных стадиях технологических процессов. Однако в настоящее время для подавляющего большинства предприятий такая информация отсутствует, и данные о выбросах представлены только для максимальной регламентной нагрузки установок (максимальные разовые выбросы, г/сек) и в целом за год (годовые или валовые выбросы, т/год).
Для определения количественных и качественных характеристик выделений и выбросов загрязняющих веществ (3В) в атмосферу используются инструментальные и расчетные (балансовые, а также основанные на удельных технологических нормативах или закономерностях протекания физико-химических процессов) методы (1). К расчетным методам, как правило, относятся расчетно-аналитические методы, в которых в качестве параметров расчетных формул для определения величин выброса (г/с) используются значения многократно измеренных концентраций вредных веществ (мг/м3) в атмосферном воздухе для типовых источников выделения. Выбор методов определения количественных и качественных характеристик выделений и выбросов 3В в атмосферу зависит, в первую очередь, oт характера производства и типа источника.
Инструментальные методы являются превалирующими для ИЗА с организованным выбросов ЗВ в атмосферу (организованные источники), к которым, в основном, относятся:
· дымовые и вентиляционные трубы, вентиляционные шахты;
· аэрационные фонари;
· дефлекторы.
Расчетные методы применяются, в основном, для определения характеристик источников с неорганизованными выделениями (выбросами). К ним относятся:
· неплотности оборудования (в т.ч. работающего при избыточном давлении);
· погрузочно-разгрузочные работы;
· открытое хранение сырья, материалов и отходов;
· оборудование и технологические процессы как в производственных помещениях, не оснащенных вентиляционными установками, так и расположенные на открытом воздухе (например, передвижные сварочные посты, резервуары хранения нефтепродуктов и т.д.);
· пруды-отстойники и накопители, нефтеловушки, шлако- и хвостохранилища, открытые поверхности испарения и т.п.; взрывные работы; открытые стоянки автотранспорта;
· передвижные источники, эксплуатируемые на производственной территории (автотранспорт, тепловозы, дорожная и строительная техника, речные суда и т.п.).
При этом, как подчеркивается в (1), могут использоваться только методики, рекомендованные к применению в установленном порядке перечнем (13) или документами, дополняющими и корректирующими этот перечень.
В настоящее время данные инвентаризации содержат сведения о типе источника, координатах (x1,y1), (x2,y2) (м) его расположения на плоской карте (схеме) местности, высоте выброса H (м), геометрии выходного канала D (м). Кроме того, указывается выбрасываемый в единицу времени объем Vg (м3/сек) (или вертикальная скорость Wo(м/сек)) газовоздушной смеси, ее температура Tg (оС), показатели разового выброса Qi (г/сек) и суммарного за год (валового) выброса Qгi (т/год), i=1,…,I, где I - число выбрасываемых в атмосферу ЗВ.
Источники по типу подразделяются на точечные, линейные и площадные. Для точечного источника (x2,y2)=0. Линейные представляются прямолинейными отрезками конечной длины, при этом (x1,y1), (x2,y2) есть координаты концов отрезка. Площадные описываются прямоугольником с центром в точке (x1,y1) и длиной сторон (x2,y2). Для площадного ИЗА указывается острый угол, отсчитываемый против часовой стрелки от оси x до x2. В качестве показателей выбросов Qi и Qгi для линейных и площадных источников указываются суммарные значения, что предполагает равномерное выделение с единицы длины или площади. Источники произвольной формы приближаются набором отрезков или прямоугольников. Аналогично приближаются источники с неравномерным выбросом по протяженности или площади.
По величине погрешности данные об источниках можно разделить на две группы. Первая - характеристики источника как сооружения (высота H, диаметр D) и координаты его расположения на территории. Измерение таких данных осуществляется достаточно просто, в зависимости от потребностей задачи они могут быть легко уточнены и изменяются только при коренной перестройки промышленного предприятия. Поэтому при постановке задач можно предполагать, что перечисленные параметры ИЗА, взятые по результатам инвентаризации, заданы с малой погрешностью даже для локальных задач переноса аэрозолей. Вторая группа – характеристики условий и количества выхода ЗВ из источника: Tg, Vg (Wo), Q. Температура выброса является весьма устойчивой величиной для однотипных промышленных установок и погрешность ее задания на превосходит 2-3% для энергетических установок и 10-15% для остальных промышленных ИЗА. Большей изменчивостью обладает объем Vg (или скорость Wo газо-воздушной смеси), которые по имеющимся оценкам может изменяться для некоторых ИЗА в 2 раза. Такой важный параметр источника как выбрасываемое количество ЗВ в единицу времени Q также обладает значительной изменчивостью, и точность его получения весьма невысока. Реальные данные по конкретным предприятиям в большинстве случаев получают на основе расчетных методов, поэтому теоретически их погрешность неизвестна. На практике же специалисты оценивают ее величиной 50-100% для конкретных источников и 20% для предприятий в целом (20).
Из приведенного краткого рассмотрения следует, что нормативная информация о выбросах промышленных предприятий, представленная в инвентаризации, пригодна для моделирования максимального (причем без указания его наступления во времени) или среднего за достаточно длительный (сезон, год) период загрязнения атмосферы промышленными аэрозолями. Поскольку зависимость мощности разового выброса Q от времени не может быть в масштабе города (региона) напрямую получена по данным инвентаризации, моделирование временного хода загрязнения с учетом зависимости выбросов от времени может быть осуществлено лишь весьма приближенно, на основе косвенных данных по суточной или сезонной производительности промышленных объектов. Исключения составляют случаи специальных экспериментов, во время которых ИЗА оборудуются измерительной аппаратурой для контроля выбросов. А поскольку количество источников по многим ЗВ в промышленном городе измеряется сотнями, то очевидно, что такие эксперименты реальны только для уникальных ЗВ, допускающих надежные методы оперативного инструментального контроля и выбрасывающихся в атмосферу города незначительным числом ИЗА.
Можно также обоснованно полагать, что значения выбросов, полученные по данным инвентаризации в масштабе города (региона), если и содержат ошибку, то систематическую. Это связано с тем, что все выбросы получены на основе единообразных методов расчета, а случайные ошибки с определенной надежностью фильтруются контролирующими органами на этапе приемки и согласования инвентаризации. Тем не менее, даже при наличии систематической ошибки данные инвентаризации остаются пригодными для сравнительного анализа воздействия различных объектов на загрязнение атмосферы.
Что касается постановки и решения обратных задач по управлению выбросами в атмосферу, то тут следует учитывать как высокую чувствительность решения таких задач к погрешностям исходных данных, так и ограниченные возможности ИЗА реальных предприятий по оперативному изменению выбросов. Большинство промышленных установок не допускают отклонений от регламентного режима. Остановка многих из них либо вообще невозможна, либо сопровождается увеличением выбросов в атмосферу. Поэтому практическая ценность математических моделей (например - оптимизационных) оперативного управления разовыми выбросами весьма ограничена. А реальную пользу по планированию первоочередных АОМ могут принести модели оптимального планирования долгосрочных АОМ, где в качестве управляемых параметров являются суммарные годовые и максимальные за год разовые выбросы ИЗА.
1.3 Метеопараметры
Состояние локального загрязнения приземного слоя воздуха существенно зависит от метеорологических условий. Хорошо известно, что при одних и тех же параметрах выбросов ИЗА, в зависимости от метеоусловий, концентрация у земли может меняться на порядок и более.
С точки зрения распространения ЗВ в атмосфере метеоусловия подразделяются на нормальные и аномальные (10). Нормальные характеризуются, прежде всего, наличием ярко выраженного среднего направления ветра. Таковыми в крупных городах являются условия со скоростью ветра более 1-2 м/сек. При меньших скоростях (штиль или близкое к штилевому состояние) в результате рельефных особенностей и температурной неоднородности подстилающей поверхности могут образовываться локальные циркуляционные зоны, приводящие к накоплению ЗВ в слое дыхания (17,18). Ситуация становится особенно опасной при наличии вертикальной температурной инверсии, препятствующей уносу примеси в верхние слои атмосферы. Именно при таких метеоусловиях фиксируются максимальные уровни загрязнения при инструментальных наблюдениях.
Подразделение метеоусловий на нормальные и аномальные играет важную роль для осознания результатов инженерных расчетов загрязнения атмосферы. Дело в том, что все инженерные модели применимы только при нормальных метеоусловиях, поскольку единое направление ветра и его стационарность являются их непременным условием. Поэтому расчетная "максимальная" концентрация является не абсолютным максимумом загрязнения, а наибольшей из концентраций для нормальных условий. Даже если предположить, что методика расчета и параметры выбросов полностью соответствуют происходящему в природе, то превышение расчетного максимума все является равновозможным и зависит от частоты появления аномальных неблагоприятных метеоусловий.
В рамках любой локальной стационарной модели наиболее важными с точки зрения рассеяния примесей метеорологическими параметрами являются скорость и направление ветра, а также показатели диффузионной активности (устойчивости) атмосферы. Скорость и направление ветра измеряются непосредственно. Обзор методов измерений скорости ветра показывает, что относительная погрешность составляет от 15% (при скоростях порядка 5 м/с) до 55% (при скорости порядка 1м/с) (26). С точки зрения решения задач переноса аэрозолей локального масштаба представляет интерес, что погрешности при определении направления ветра могут привести к ошибке положения оси дымового факела на карте территории порядка 10-15%, в результате чего при достаточно устойчивой стратификации атмосферы факел просто не накроет расчетную точку и приведет к большой ошибке моделирования. Это следует учитывать при интерпретации понятия "опасное направление ветра" и определении основных виновников загрязнения заданной точки города. Сказанное еще раз подчеркивает, что на практике при решении краткосрочных задач нормирования выбросов в смысле неравенства (1) представляет интерес предсказание с разумной точностью максимума разовой концентрации даже без указания момента времени, когда это произойдет.
1.4 Данные наблюдений за загрязнением атмосферы
При обосновании системы мониторинга на территории России академик Ю.А. Израэль подчеркивает, что "только регулярные наблюдения в строго определенных местах и в строго установленные сроки являются источниками прямой и статистически обеспеченной информации о загрязнении окружающей среды"(23). Такого рода наблюдения применительно к загрязнению атмосферы осуществляются на сети стационарных постов Росгидромета в городах . Посты расположены только в крупных городах. Например, в Кемеровской области стационарные посты оборудованы в Новокузнецке (10), Кемерово (9) и Белово (1). Причем количество постов сокращается (в 90-х годах в Кемерово было 12 постов). На постах ежедневно (в 7, 13 и 19 часов местного времени) осуществляются отборы проб воздуха, которые доставляются в лабораторию местного подразделения Росгидромета, где анализируются стандартизированными методами. Контролируются только незначительное число наиболее распространенных ЗВ (10-20 примесей), в то время как в данных инвентаризации совокупности предприятий крупного промышленного города встречается на практике 100-200 веществ. Таким образом, большинство ЗВ на сети не контролируется ничем, кроме интуиции разработчиков проектов и согласующих эти проекты экспертов.
Обзор методов инструментального анализа воздуха (21,22) показывает, что количественные оценки погрешностей различных этапов лабораторных методов анализа полученных проб составляют 6-25% (с доверительной вероятностью 95%).
1.5 Модели расчета загрязнения атмосферы
Даже при наличие данных наблюдений за загрязнением воздуха на сети стационарных постов, одной из важнейших наукоемких задач охраны атмосферы является расчет загрязнения заданной территории по имеющимся данным о выбросах и условия распространения примесей. Действительно, оценить перспективный уровень загрязнения в зависимости от варианта промышленного развития можно только расчетными методами. Кроме того, методы инструментальных наблюдений в общем случае не могут указать вклад отдельного источника (предприятия) в измеряемую суммарную величину, что необходимо при определении основных виновников загрязнения, установлении ПДВ и начислении платы за выброс. Поэтому моделирование загрязнения атмосферы необходимо как для анализа, так и для прогноза.
В предисловии к (17) А.М. Яглом подчеркивает, что расчет диффузии примеси в атмосфере "не может быть сведен к какой-то задаче математической физики, а обязательно требует тех или иных нестрогих гипотез и приближенных допущений. По этой причине задача о распространении примеси в атмосфере не имеет одного общепринятого ‘правильного решения’, а характеризуется наличием целого ряда различных подходов к требуемому расчету, ни один из которых не может претендовать на полную строгость и точность".
Модели, естественно, можно классифицировать с различных точек зрения, и этому посвящена весьма обширная литература (10,17,18,20). С позиций, используемых для построения научных теорий, модели подразделяются на статистические и полуэмпирические. Статистические основаны на том предположении, что поступающее из источника ЗВ переносится вместе со средним потоком, а его распространение в поперечном потоку направлении происходит под воздействием вихрей, движение которых подчиняется определенным статистическим закономерностям. Полуэмпирические используют для получения результата те или иные решения уравнения турбулентного переноса в предположение об аналогии между турбулентной и молекулярной диффузией. Слово "полуэмпирические" подчеркивает, что для задания коэффициентов турбулентной диффузии необходимы эмпирические предположения, и только после этого можно начинать поиск точного или численного решения уравнения переноса.
В зависимости от времени осреднения рассчитываемой концентрации локальные модели можно подразделить на краткосрочные (разовые) и долговременные. Как статистические, так и полуэмпирические краткосрочные модели рассчитывают концентрацию атмосферной примеси, осредненную за 20-30 минут. Долгосрочные предназначены для оценки загрязнения, осредненного за большие промежутки времени (сезон, год) и основаны на осреднении разовых расчетов с использованием повторяемости характерных для данной территории условий распространения ЗВ (17,24,27,31).
С точки зрения простоты использования модели также можно разбить на два класса: научно-исследовательские и инженерные. Первые, будучи способными описывать достаточно тонкие особенности распространения ЗВ в атмосфере (сложный рельеф, штиль, особенности турбулентного режима, локальные циркуляции и т.д.) (17,18), являются безусловно более сложными, требуют высокой квалификации персонала как для применения моделей, так и для получения специфических и дорогостоящих исходных данных. Вторые, предназначенные для проектных расчетов, доведены до однозначно трактуемых числовых зависимостей, табличных и графических аппроксимаций и обеспечены системой сбора (расчета) исходной информации для возможности реализации необходимых количественных оценок в процессе выполнения проектных работ (1,2). Можно сказать, что именно наличие реальной на сегодняшний день системы обеспечения модели исходными данными является признаком того, что сама модель представляет практический интерес для целей управления (нормирования) промышленных выбросов.
Наиболее распространенной для инженерных приложений и принятия административных решений за рубежом является локальная модель Гауссовского факела (17,20), в основе которой лежит статистическая теория с рядом упрощающих предположений и большим количеством эмпирических таблиц для задания дисперсионных коэффициентов (20), определяющих процесс расширения факела при удалении от источника. Краткосрочная модель позволяет рассчитать в заданной точке территории разовую (среднюю за 20-30 минут) концентрацию в зависимости от трех параметров – скорости ветра, направления ветра и класса устойчивости (температурной стратификации) атмосферы. Гауссовская модель долговременного осреднения позволяет оценить среднюю за заданный период концентрацию по совместному распределению (за этот период) указанных трех параметров. Для практического использования модели достаточно иметь данные стандартных метеорологических наблюдений и параметров источников в объеме, рассмотренном в пункте 2.
Принятая на федеральном уровне Минприродой России (2) для нормативных расчетов модель атмосферного переноса ЗВ предназначена для расчета максимальных разовых концентраций См атмосферных примесей, фигурирующих в условии (1.1). Модель была создана в конце 60-х годов и в дальнейшем доведена до инженерных формул коллективом разработчиков под руководством проф. М.Е. Берлянда в Главной геофизической обсерватории им. А.И. Воейкова. В настоящее время она оформлена в виде руководящего документа ОНД-86 (2), вышедшего взамен ранее действующего СН-369-74 (с рядом уточнений, учетом коэффициента рельефа и застройки). Модель основана на численном решении стационарного полуэмпирического уравнения турбулентной диффузии для различных типов источников с последующей аналитической аппроксимацией. Численные расчеты проводились при атмосферных стратификациях, которым соответствуют наибольшие значения максимума приземной концентрации, достигаемые при “опасной” скорости ветра. Отсюда следует, что модель (2) изначально нацелена на расчет верхней оценки приземных концентраций при нормальных (соответствующих условиям стационарности) метеоусловиях. Поэтому аналитические аппроксимации, в отличие от Гауссовских моделей, позволяют оценивать концентрации ЗВ только в зависимости от двух параметров – скорости и направления ветра. А поскольку невозможно получить априорного аналитического решения для нахождения максимума (по скорости и направлению) суммарной концентрации для произвольного множества произвольно расположенных ИЗА, то указанная задача решается численно для каждого конкретного случая либо методом простого перебора, либо на основе анализа характерных особенностей заданного множества ИЗА. Для применения модели нужно знать, кроме параметров ИЗА, коэффициент рельефа территории, коэффициент осаждения исследуемой примеси, параметр характерной “опасной” стратификации атмосферы (все три этих параметра специально ориентированны на модель (2)) и климатические показатели, определяющие для заданной территории режим температуры окружающего воздуха и скорости ветра. Соответствующие количественные характеристики, вместе с фоновыми концентрациями, выдаются для проектных расчетов территориальными подразделениями Росгидромета.
1.6 Унифицированные программы расчета загрязнения атмосферы
С середины 70-х годов начали разрабатываться программы расчета загрязнения атмосферы (ПРЗА) для широко распространенных в то время ЭВМ типа БЭСМ и М-20. Наличие достаточно сложной вычислительной основы при расчете поля приземных концентраций и многовариантности методики (2) приводили к естественным ошибкам в алгоритмах и программных реализациях и, как следствие, к несопоставимости ре