СПОСОБ ПРОСТРАНСТВЕННОЙ КОЛИЧЕСТВЕННОЙ ОЦЕНКИ УРОВНЯ ЗАГРЯЗНЕНИЯ АТМОСФЕРНОГО ВОЗДУХА Российский патент 2013 года по МПК G01W1/00 

Изобретение относится к охране среды обитания, в частности к определению уровня загрязнения атмосферного воздуха, и может быть использовано при контроле чистоты воздуха населенных мест.

Важнейшей проблемой санитарно-эпидемиологического контроля территорий является получение достоверной документальной количественной информации, привязанной к координатам местности.

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

Например, известны способы оценки степени антропогенного загрязнения воздушной среды с использованием фитоиндикации, в частности, по проценту пораженной растительной ткани (1), по количеству содержащихся в ней биохимических индикаторов (2). Главным существенным недостатком этих способов является низкая оперативность, т.к. морфологические и структурные перестройки проявляются в растениях под действием загрязнения не сразу, а через длительный промежуток времени.

Также известны способы экологического зонирования территории путем проведения космической съемки и последующей программной обработки полученного изображения по спектральной яркости (3, 4). Однако указанные известные способы имеют следующие недостатки: неоперативны, отсутствуют количественные характеристики качества атмосферного воздуха, обладают низкой точностью при изучении небольших по площади территорий, имеющих значимые различия уровней поля приземных концентраций атмосферного воздуха.

Из уровня техники известен способ выявления зон химического загрязнения атмосферного воздуха (5), согласно которому исследуемую территорию разбивают на условную сетку квадратов, размер которых определяется в зависимости от того, насколько подробно следует получить границу загрязнения. Проводят отбор пробы в центре квадрата. С использованием биосенсора на основе суспензии подвижных микроорганизмов определяют величину индекса токсичности пробы и при ее величине меньше минимально допустимого или больше максимально допустимого значения естественного природного фона область, характеризуемую пробой, считают принадлежащей загрязненной зоне.

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

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

Известен способ определения предельно допустимого выброса (ПДВ) токсических веществ для каждого промышленного предприятия, который устанавливают на уровне, обеспечивающем соблюдение гигиенических нормативов в воздухе населенных мест при наиболее неблагоприятных для рассеивания метеоусловиях (7). Подобные расчетные методики, реализованные в программе "Эфир-5", требуют учета большого числа параметров, которые, как правило, по своей природе нестационарны, ориентированы на средние значения. В результате существенно снижается точность оценки фактического поля приземных концентраций. Внесение дополнительных конкретизирующих данных значительно увеличивает трудоемкость расчетов, а сложности учета всех влияний являются источником систематической ошибки при вычислениях. Таким образом, объективно оценить уровень загрязнения воздуха можно только в местах отбора проб.

Наиболее близким к предлагаемому изобретению является способ контроля качества атмосферного воздуха (8) путем измерения концентрации загрязняющих веществ в точках контроля, сравнение измеренных концентраций с предельно допустимыми значениями (ПДК) и осуществление контрольных операций по анализу влияния на точки контроля источников выброса в атмосферу загрязняющих веществ, концентрации которых превышают величину ПДК, при этом дополнительно измеряют азимут и угол рассеивания загрязнений, а источники выброса, на которых осуществляют контрольные операции, выбирают по измеренным азимуту ветра и угла рассеивания загрязнений.

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

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

Указанный технический результат достигается предлагаемым способом пространственной количественной оценки уровня загрязнения атмосферного воздуха, согласно которому проводят выбор территории, которую необходимо исследовать на предмет состояния уровня загрязнения атмосферного воздуха, на исследуемой территории в точках натурных замеров проводят натурные инструментальные замеры концентраций загрязняющих веществ в атмосферном воздухе с фиксированием даты, времени отбора пробы воздуха и метеорологических характеристик атмосферы в момент отбора указанной пробы, при этом новым является то, что карту исследуемой территории покрывают регулярной сеткой, выделяют узловые точки в местах пересечения линий сетки и на указанной карте отмечают расположение точек, в которых были проведены натурные замеры концентраций загрязняющих веществ, выявляют все источники загрязнения атмосферного воздуха на исследуемой территории и проводят сбор данных о параметрах выбросов от указанных источников загрязнения атмосферного воздуха, выполняют расчет приземных концентраций загрязняющих веществ в точках натурных замеров и в узловых точках заданной регулярной сетки от указанных источников загрязнения атмосферного воздуха с применением стандартных математических моделей и программных средств, для каждой точки, где имеются данные и расчета рассеивания, и натурных замеров, определяют коэффициент соответствия как отношение измеренной концентрации к рассчитанной, далее точки натурных замеров концентраций загрязняющих веществ объединяют на карте непересекающимися отрезками в треугольники, образуя систему треугольников с вершинами в точках натурных замеров, для каждого треугольника решают уравнение плоскости с установлением коэффициентов уравнения, зависящих от координат x и y вершин треугольника - точек натурных замеров, и значений коэффициентов соответствия в них, далее относят каждую узловую точку расчетной сетки с координатами x_i и y_i к какому-либо треугольнику образованной системы треугольников или устанавливают, что она лежит вне указанной системы треугольников, для каждой узловой точки с координатами x_i и y_i, лежащей внутри системы треугольников, рассчитывают коэффициент соответствия по уравнению плоскости соответствующего треугольника, для узловых точек, лежащих вне системы треугольников, расчет коэффициента соответствия выполняют методом экстраполяции, для этого значение коэффициента соответствия в узловой точке принимают равными коэффициентам соответствия в ближайшей точке, лежащей на внешней границе системы треугольников, ранее рассчитанные приземные концентрации загрязняющих веществ в узловых точках заданной регулярной сетки умножают на полученные коэффициенты соответствия с получением уточненной концентрации загрязняющих веществ в узловых точках сетки, затем строят карту пространственного распределения уточненных концентраций загрязняющих веществ, по которой количественно оценивают уровень загрязнения атмосферного воздуха на исследуемой территории.

Достижение технического результата обеспечивается за счет следующего.

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

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

Благодаря последующему проведению расчета рассеивания загрязняющих веществ от указанных источников загрязнения атмосферного воздуха с применением стандартных математических моделей и программных средств (например, модель «Методики расчета концентраций в атмосферном воздухе вредных веществ, содержащихся в выбросах предприятий» - ОНД 86, модель рассеивания Гаусса, программные продукты УПРЗА «Эколог», «Gaussian Dispersion Model Calculator" и др.) обеспечивается получение расчетных концентраций загрязняющих веществ в узловых точках заданной регулярной сетки (каждая ячейка сетки характеризуется одной узловой точкой) и в точках непосредственного измерения качества воздуха. В результате формируется информационная база с пространственно дифференцированными характеристиками уровня загрязнения атмосферы в каждой расчетной точке, в том числе в ячейках регулярной сетки.

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

Затем производится интерполяция значений коэффициентов соответствия, установленных для реперных точек, на узлы сетки.

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

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

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

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

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

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

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

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

Другими словами, моделирование пространственного распределения загрязнения по данным натурных инструментальных исследований методами интер- и экстраполяции с использованием результатов расчетов рассеивания позволяет минимизировать неопределенности каждого метода (только метода натурных измерений или только метода расчетов рассеивания) в отдельности и получить наиболее точные результаты при условии корректной аппроксимации данных при реализации предлагаемого способа.

Предлагаемый способ иллюстрируется рядом чертежей, где на Рис.1 отображено разбиение пространства исследуемой территории с помощью триангуляции Делоне (один из наиболее корректных способов триангуляции), причем жирными точками отмечены координаты точек натурных измерений; на Рис.2 - система треугольников с указанием номеров точек натурных замеров; Рис.3 - иллюстрация к установлению принадлежности каждой узловой точки, расположенной внутри системы треугольников, к одному из ряда получившихся треугольников; Рис.4 - приведен результат интерполяции коэффициентов соответствия для узловых точек, расположенных внутри системы треугольников; Рис.5 - показаны аппроксимированные значения коэффициента соответствия во всех узлах регулярной сетки; Рис.6 - карта-схема пространственного распределения расчетных максимально-разовых концентраций загрязняющих веществ (азота диоксида) на исследуемой территории и коэффициенты соответствия расчетных и натурных данных до проведения процедуры аппроксимации; Рис.7 - карта-схема пространственного распределения расчетных максимально-разовых концентраций загрязняющих веществ (азота диоксида) на исследуемой территории и коэффициенты соответствия расчетных и натурных данных после проведения процедуры аппроксимации; Рис.8 - представлены коэффициенты соответствия в пяти контрольных точках, выбранных для подтверждения точности и достоверности предлагаемого способа, до проведения процедуры аппроксимации (в данном случае - на примере полей распределения максимально-разовых концентраций азота диоксида); на Рис.9 представлены коэффициенты соответствия в пяти контрольных точках, выбранных для подтверждения точности и достоверности предлагаемого способа, после проведения процедуры аппроксимации.

Предлагаемый способ осуществляют следующим образом, реализуя его на конкретном примере.

1. Например, была выбрана территория N площадью почти 800 кв. км, на которой размещены крупные отраслевые комплексы тяжелой промышленности, электроэнергетики, нефтегазопереработки, машиностроения, химии и нефтехимии, деревообработки, полиграфии и прочие, а также ряд крупных автомагистралей.

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

3. Карту исследуемой территории покрывают регулярной сеткой (Рис.1) с шагом 100×100 м (шаг сетки определяется необходимой конечной точностью результатов анализа). Сетка размером 50,5 км на 34,8 км состоит из 176 594 ячеек и 7 реперных точек.

4. На указанной карте отмечают расположение семи реперных точек, т.е. тех, в которых выполнены инструментальные натурные замеры качества атмосферного воздуха по 23 приоритетным химически опасным веществам. Замеры проводились ежедневно дискретно через равные промежутки времени 3-4 раза в день в 1, 7, 13, 19 часов по местному времени. Необходимым условием корректной оценки уровня среднегодового загрязнения атмосферы является проведение не менее 200 разовых натурных замеров концентрации каждого из установленных приоритетных химически опасных веществ в течение не менее 50 дней в каждой из точек расположения стационарных или передвижных постов наблюдения. В таблице 1 приведен фрагмент измеренных максимально-разовых концентраций в точках расположения постов замеров.

5. На исследуемой территории выявляют все источники загрязнения атмосферного воздуха. Например, на данной территории расположено более 11500 источников, выбрасывающих в атмосферу более 450 химических веществ, и 25 линейных участков улично-дорожной сети.

6. Затем осуществляют сбор данных о метеорологических характеристиках атмосферы мест расположения источников выбросов и параметрах источников выбросов загрязняющих веществ, которые должны содержать следующие сведения (согласно ГОСТ): наименование предприятия, номер или наименование цеха, название источника выброса, номер источника загрязнения атмосферы (ИЗА), тип, высота источниках, диаметр устья, скорость и объем истекающей газовоздушной смеси, температура выброса, координаты X и Y, ширина (для площадного источника), масса выброса каждого загрязняющего вещества по каждому источнику (г/сут и т/год). Например, для источника N были сформированы данные: цех №1, источник слесарно-расточное отделение, номер источника 1, тип источника 1, высота 10,50 м, диаметр устья 0,40 м, скорость 3,00 м/с, объем 0,377 м³/с, температура 18°С, координаты Х - 1795, Y - 3109 (в системе координат изучаемой территории), ширина источника 0,00 м, железа оксид (код 0123) - 0,003 г/с, 8,52202 т/год, масло минеральное нефтяное (код 2735) - 0,005 г/с, 12,83675 т/год, эмульсол (код 2868) - 0,00001 г/с, 0,03534 т/год, корунд белый (код 2930) - 0,0006 г/с, 1,80676 т/год. Участки улично-дорожной сети рассматриваются как линейные стационарные источники.

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

8. В 7 реперных точках (точках натурных замеров) определяют коэффициент соответствия K_i как отношение измеренной (фактической) концентрации к рассчитанной:

$$K_i=\frac{C_i^f}{C_i^r}$$

где i - номер поста;

$$C_i^f$$ - фактические концентрации загрязняющего вещества на i-м посту наблюдений;

$$C_i^r$$ - расчетные концентрации загрязняющего вещества на i-м посту наблюдений;

В таблице 3 приведен фрагмент рассчитанных коэффициентов соответствия в точках расположения постов замеров.

9. Для интерполяции коэффициента соответствия на основании данных в точках постов наблюдения проводят процедуру триангуляции. Для этого точки натурных замеров концентраций загрязняющих веществ объединяют на карте с регулярной сеткой в треугольники, на базе которых, используя, например, метод триангуляции Делоне (или другой адекватный метод триангуляции), определяют коэффициенты соответствия во всех узлах регулярной расчетной сетки. Более подробно это определение осуществляют следующим образом.

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

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

Для каждого треугольника решают уравнение плоскости с установлением коэффициентов уравнения, зависящих от координат x и y вершин треугольника (точек натурных замеров) и значений коэффициентов соответствия в них. Уравнение плоскости представляет собой непрерывную линейную функцию двух переменных, которая может быть записана в следующем виде (2):

$$K\left(x,y\right)=a_0+a_{1}x+a_{2}y,\left(2\right)$$

где

-a ₀, a ₁, a ₂ - произвольные постоянные коэффициенты;

- K(x,y) - коэффициент соответствия для точки с координатами x, y;

- значения функции в вершинах треугольника, соответствующие значению коэффициента соответствия на постах образующих этот треугольник, обозначают, как k₁, k₂, k₃;

- с помощью координат вершин треугольника и значений k₁, k₂, k₃ вычисляются неизвестные постоянные коэффициенты a ₀, a ₁, a ₂.

- получают систему трех линейных алгебраических уравнений относительно неизвестных коэффициентов a ₀, a ₁, a ₂ (3):

$$k_i\equiv K\left(x_i,y_i\right)=a_0+a_1{x}_i+a_2{y}_i,i=\overline{\mathrm{1,3,}}\left(3\right)$$

где i - номер вершины треугольника.

- решив систему (3), получают однозначное выражение функции (2) через ее узловые значения (значения в узлах сетки) для каждого треугольника.

Определяют принадлежность каждой узловой точки, расположенной внутри указанной системы треугольников, к одному из ряда семи получившихся треугольников по следующему алгоритму:

Указанная точка, расположенная внутри указанной системы треугольников, соединяется отрезками с вершинами каждого из треугольников с образованием в свою очередь трех треугольников с площадью S₁, S₂, S₃ (Рис.3).

При этом, если площадь исходного треугольника S равна сумме площадей образовавшихся трех треугольников S=S₁+S₂+S₃, где S₁ - площадь треугольника 1; S₂ - площадь треугольника 2; S₃ - площадь треугольника 3, то считается, что точка принадлежит данному треугольнику. Если же S<(S₁+S₂+S₃), то данная точка не принадлежит данному треугольнику.

Площадь треугольников вычислялась по следующим формулам:

$$S=\frac{\left|\left(x_2-x_1\right)\left(y_3-y_1\right)-\left(x_3-x_1\right)\left(y_2-y_1\right)\right|}{2}$$ ;

$$S_1=\frac{\left|\left(x_t-x_1\right)\left(y_3-y_1\right)-\left(x_3-x_1\right)\left(y_t-y_1\right)\right|}{2}$$ ;

$$S_2=\frac{\left|\left(x_2-x_t\right)\left(y_3-y_t\right)-\left(x_3-x_t\right)\left(y_t-y_t\right)\right|}{2}$$ ;

$$S_3=\frac{\left|\left(x_2-x_1\right)\left(y_t-y_1\right)-\left(x_t-x_1\right)\left(y_2-y_1\right)\right|}{2}$$ .

Затем рассчитывают значения коэффициента соответствия во всех точках (узлах регулярной сетки), внутри каждого треугольника по уравнению (2) с использованием полученных коэффициентов для соответствующего треугольника.

На Рис.4 приведен результат интерполяции коэффициентов соответствия для узловых точек, расположенных внутри системы треугольников (более интенсивный цвет отражает более высокие значения коэффициентов соответствия).

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

11. В результате этих действий получают аппроксимированные значения коэффициента соответствия во всех узлах регулярной сетки (Рис.5).

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

12. Производят расчет уточненных концентраций загрязняющих веществ в каждой узловой точке расчетной сетки на исследуемой территории согласно формуле (4):

$$C^u\left(x,y\right)=K\left(x,y\right)\cdot C^r\left(x,y\right),где\left(4\right)$$

C^u - уточненные концентрации загрязняющего вещества в расчетной точке (x,y),

K - коэффициент соответствия в расчетной точке (x,y);

C^r - расчетные концентрации загрязняющего вещества в расчетной точке (x,y).

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

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

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

На Рис.6 приведены карта-схема пространственного распределения расчетных максимально-разовых концентраций загрязняющих веществ (азота диоксида) на исследуемой территории и коэффициенты соответствия расчетных и натурных данных до проведения процедуры аппроксимации.

На Рис.7 приведены карта-схема пространственного распределения расчетных максимально-разовых концентраций загрязняющих веществ (азота диоксида) на исследуемой территории и коэффициенты соответствия расчетных и натурных данных после проведения процедуры аппроксимации.

На Рис.6 видно, что до аппроксимации при анализе соответствия расчетных и натурных данных в точках расположения постов лабораторного мониторинга выявлены различия результатов, о чем свидетельствуют вариации коэффициентов соответствия для анализируемого химического вещества (в данном случае азота диоксида) в интервале 0,85-10,48.

После проведения процедуры аппроксимации, как видно на Рис.7, картина распределения полей концентраций загрязняющих веществ изменяется, при этом в точках натурного мониторинга (постах) коэффициенты соответствия становятся равными 1, а в ряде контрольных точек, которые были выбраны на исследуемой территории для подтверждения достоверности полученных результатов, сходимость расчетных и натурных данных возросла до 70-95% при ранее отмечаемых 8-50% (Рис.8, 9).

На рисунке 8 представлены коэффициенты соответствия в 5 контрольных точках, выбранных для подтверждения точности и достоверности предлагаемого способа, до проведения процедуры аппроксимации (в данном случае - на примере полей распределения максимально-разовых концентраций азота диоксида), на рисунке 9 представлены коэффициенты соответствия в 5 контрольных точках, выбранных для подтверждения точности и достоверности предлагаемого способа, после проведения процедуры аппроксимации.

Другими словами, доказательством того, что предлагаемый способ пространственной количественной оценки уровня загрязнения атмосферного воздуха является точным и достоверным, является факт, что после всей поэтапной процедуры аппроксимации полученные результаты в контрольных точках стали близкими по значению, а коэффициенты соответствия с 0,76-10,58 достигли значений 0,82-1,12, т.е. значительно выросла сходимость натурных и полученных в результате аппроксимации данных. Следовательно, верифицирующие измерения подтвердили корректность оценки качества приземного слоя атмосферного воздуха в узловых точках сетки.

Таким образом, заявляемый способ позволяет:

- максимально учесть все источники загрязнения атмосферного воздуха при интер- и экстраполяции данных сети постов натурного наблюдения на всю исследуемую территорию;

- максимально нивелируя недостатки и используя достоинства расчетного и натурного методов, добиться повышения соответствия данных о загрязнении атмосферного воздуха, полученных с помощью расчетов рассеивания по методике ОНД-86 и данных по реальным концентрациям загрязняющих веществ, полученных с помощью лабораторных замеров на местности;

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

Таблица 1 Фрагмент измеренных максимально-разовых концентраций в точках расположения постов замеров   Координаты поста, м Концентрации веществ, мг/м³ N поста X Y Азота диоксид Аммиак Углерода оксид Фенол Формальдегид Взвешенные вещества 1 8781 9801 0,120 0,120 12,000 0,008 0,065 0,400 2 3274 1080 0,390 0,470 6,000 0,012 0,162 0,400 3 1439 -3189 0,610 0,160 32,000 0,020 0,137 0,500 4 -1022 -1671 0,380 0,200 6,000 0,009 0,113 0,600 5 -3919 -6228 0,310 0,200 5,000 0,020 0,086 0,500 6 -18425 -764 0,410 0,190 4,000 0,012 0,164 0,400 7 3578 -2062 0,130 0,310 14,000 0,015 0,109 0,400

Таблица 2 Фрагмент результатов расчетов концентрации распространения ряда химических загрязняющих веществ на исследуемой территории Координаты узлов сетки Концентрации веществ, мг/м³ Х Y Азота диоксид Аммиак Углерода оксид Фенол Формальдегид Взвешенные вещества -25540 -16543 0,01586119 0,00323080 0,15707225 0,00007355 0,00165417 0,00135560 -25440 -16543 0,01593381 0,00325589 0,15720700 0,00007394 0,00166386 0,00136710 -25340 -16543 0,01600856 0,00328130 0,15788955 0,00007432 0,00167335 0,00137880 -25240 -16543 0,01608345 0,00330702 0,15779110 0,00007469 0,00168264 0,00139069 -25140 -16543 0,01615838 0,00333307 0,15859545 0,00007505 0,00169182 0,00140277 -25040 -16543 0,01624157 0,00335944 0,15870280 0,00007553 0,00170350 0,00141505 -24940 -16543 0,01632229 0,00338614 0,15919865 0,00007603 0,00171500 0,00142754 -24840 -16543 0,01640275 0,00341317 0,15950485 0,00007652 0,00172629 0,00144024 -24740 -16543 0,01648380 0,00344054 0,15970675 0,00007702 0,00173773 0,00145316 -24640 -16543 0,01656397 0,00346825 0,16020445 0,00007750 0,00174863 0,00146693 -24540 -16543 0,01664390 0,00349629 0,16012680 0,00007798 0,00175929 0,00148099 -24440 -16543 0,01672359 0,00352468 0,16080810 0,00007844 0,00176972 0,00149530 -24340 -16543 0,01680300 0,00355343 0,16063125 0,00007889 0,00177989 0,00150984 -24240 -16543 0,01688208 0,00358254 0,16132170 0,00007933 0,00178978 0,00152465 -24140 -16543 0,01696462 0,00361202 0,16130405 0,00007977 0,00180111 0,00153970 -24040 -16543 0,01705017 0,00364186 0,16175015 0,00008019 0,00181368 0,00155503 -23940 -16543 0,01713546 0,00367208 0,16188460 0,00008060 0,00182602 0,00157062 -23840 -16543 0,01722200 0,00370268 0,16209800 0,00008116 0,00183811 0,00158649 -23740 -16543 0,01730993 0,00373367 0,16237730 0,00008173 0,00184993 0,00160264 -23640 -16543 0,01739744 0,00376504 0,16236900 0,00008230 0,00186145 0,00161907 -23540 -16543 0,01748454 0,00379681 0,16278600 0,00008285 0,00187269 0,00163580

Таблица 3 Фрагмент рассчитанных коэффициентов соответствия в точках расположения постов замеров   Координаты поста, м Коэффициент соответствия N поста X Y Азота диоксид Аммиак Углерода оксид Фенол Формальдегид Взвешенные вещества 1 8781 9801 2,02 26,63 19,17 23,57 6,29 18,26 2 3274 1080 2,50 76,86 19,62 35,47 3,42 22,43 3 1439 -3189 2,58 7,84 107,11 15,28 2,49 34,20 4 -1022 -1671 2,26 11,62 10,03 18,98 1,83 15,63 5 -3919 -6228 2,74 16,61 16,93 42,96 2,73 19,11 6 -18425 -764 10,48 20,51 6,44 18,64 23,51 89,58 7 3578 -2062 0,85 32,42 75,22 25,74 8,25 23,07

Таблица 4 Результат интерполяции коэффициентов соответствия для узловых точек, расположенных внутри системы треугольников Координаты узлов сетки Коэффициент соответствия Х Y Азота диоксид Аммиак Углерод а оксид Фенол Формальдегид Взвешенные вещества -25540 -16543 10,47967 16,60703 6,435936 42,95037 23,51028 89,57657 -25440 -16543 10,47967 16,60713 6,435935 42,95375 23,51039 89,57648 -25340 -16543 10,47964 16,60714 6,435936 42,96286 23,51032 89,57644 -25240 -16543 10,47965 16,6071 6,43594 42,95086 23,51008 89,5764 -25140 -16543 10,4797 16,60721 6,435935 42,95803 23,51018 89,57634 -25040 -16543 10,47965 16,60723 6,435936 42,96306 23,51042 89,57634 -24940 -16543 10,47966 16,60711 6,435934 42,95673 23,5102 89,57647 -24840 -16543 10,47965 16,60714 6,435936 42,95609 23,50995 89,57604 -24740 -16543 10,47968 16,60699 6,435933 42,96287 23,50998 89,57651 -24640 -16543 10,47967 16,60722 6,435939 42,95484 23,51041 89,57619 -24540 -16543 10,4797 16,60703 6,435937 42,95973 23,51005 89,57657 -24440 -16543 10,47969 16,60718 6,435938 42,96277 23,50993 89,576 -24340 -16543 10,37071 16,60705 6,583725 42,95855 23,21773 89,57638 -24240 -16543 10,3707 16,60721 6,583727 42,95979 23,21738 89,57662 -24140 -16543 10,26171 16,60705 6,731517 42,96101 22,92475 89,57654 -24040 -16543 10,26172 16,60717 6,731518 42,96047 22,92466 89,57641 -23940 -16543 10,2617 16,60721 6,731517 42,95285 22,92472 89,57609 -23840 -16543 10,15271 16,60716 6,879308 42,95219 22,63249 89,57636 -23740 -16543 10,15273 16,60698 6,879305 42,95852 22,63221 89,57657 -23640 -16543 10,04372 16,60699 7,027099 42,95261 22,33957 89,57611 -23540 -16543 10,04373 16,6071 7,027097 42,95715 22,33952 89,57635

Таблица 5 Результаты аппроксимированных полей концентраций ряда химических загрязняющих веществ на исследуемой территории Координаты узлов сетки Концентрации веществ, мг/м³ Х Y Азота диоксид Аммиак Углерода оксид Фенол Формальдегид Взвешенные вещества -25540 -16543 0,166220 0,053654 1,010907 0,003159 0,038890 0,121430 -25440 -16543 0,166981 0,054071 1,011774 0,003176 0,039118 0,122460 -25340 -16543 0,167764 0,054493 1,016167 0,003193 0,039341 0,123508 -25240 -16543 0,168549 0,054920 1,015534 0,003208 0,039559 0,124573 -25140 -16543 0,169335 0,055353 1,020710 0,003224 0,039775 0,125655 -25040 -16543 0,170206 0,055791 1,021401 0,003245 0,040050 0,126755 -24940 -16543 0,171052 0,056234 1,024592 0,003266 0,040320 0,127874 -24840 -16543 0,171895 0,056683 1,026563 0,003287 0,040585 0,129011 -24740 -16543 0,172745 0,057137 1,027862 0,003309 0,040854 0,130169 -24640 -16543 0,173585 0,057598 1,031066 0,003329 0,041111 0,131402 -24540 -16543 0,174423 0,058063 1,030566 0,003350 0,041361 0,132662 -24440 -16543 0,175258 0,058535 1,034951 0,003370 0,041606 0,133943 -24340 -16543 0,174259 0,059012 1,057552 0,003389 0,041325 0,135246 -24240 -16543 0,175079 0,059496 1,062098 0,003408 0,041554 0,136573 -24140 -16543 0,174086 0,059985 1,085821 0,003427 0,041290 0,137921 -24040 -16543 0,174964 0,060481 1,088824 0,003445 0,041578 0,139294 -23940 -16543 0,175839 0,060983 1,089729 0,003462 0,041861 0,140690 -23840 -16543 0,174850 0,061491 1,115122 0,003486 0,041601 0,142112 -23740 -16543 0,175743 0,062005 1,117043 0,003511 0,041868 0,143559 -23640 -16543 0,174735 0,062526 1,140983 0,003535 0,041584 0,145030 -23540 -16543 0,175610 0,063054 1,143913 0,003559 0,041835 0,146529

Источники информации

1. Надежкина Е.В. и др. Практикум по экологии и химии окружающей среды. Пенза, 2003. - с.193.

2. Патент РФ №2346270.

3. Патент РФ №2132606.

4. Патент РФ №2018156.

5. Патент РФ №1804482.

6. Указания по расчету рассеивания в атмосфере вредных веществ, содержащихся в выбросах предприятий. СН 369-74. - М.: Стройиздат, 1974, с.41.

7. Программа «Эфир-5» и инструкция по автоматизированному расчету загрязнения атмосферы вредными выбросами промышленных предприятий с целью установления ПДВ (ВСВ). - М., 1983, с.143

8. Авторское свидетельство РФ №1008685.

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

Способ пространственной количественной оценки уровня загрязнения атмосферного воздуха, согласно которому проводят выбор территории, которую необходимо исследовать на предмет состояния уровня загрязнения атмосферного воздуха, на исследуемой территории в точках натурных замеров проводят натурные инструментальные замеры концентраций загрязняющих веществ в атмосферном воздухе с фиксированием даты, времени отбора пробы воздуха и метеорологических характеристик атмосферы в момент отбора указанной пробы, отличающийся тем, что карту исследуемой территории покрывают регулярной сеткой, выделяют узловые точки в местах пересечения линий сетки и на указанной карте отмечают расположение точек, в которых были проведены натурные замеры концентраций загрязняющих веществ, выявляют все источники загрязнения атмосферного воздуха на исследуемой территории и проводят сбор данных о параметрах выбросов от указанных источников загрязнения атмосферного воздуха, выполняют расчет приземных концентраций загрязняющих веществ в точках натурных замеров и в узловых точках заданной регулярной сетки от указанных источников загрязнения атмосферного воздуха с применением стандартных математических моделей и программных средств, для каждой точки, где имеются данные и расчета рассеивания, и натурных замеров, определяют коэффициент соответствия как отношение измеренной концентрации к рассчитанной, далее точки натурных замеров концентраций загрязняющих веществ объединяют на карте непересекающимися отрезками в треугольники, образуя систему треугольников с вершинами в точках натурных замеров, для каждого треугольника решают уравнение плоскости с установлением коэффициентов уравнения, зависящих от координат x и y вершин треугольника - точек натурных замеров, и значений коэффициентов соответствия в них, далее относят каждую узловую точку расчетной сетки с координатами x_i и y_i к какому-либо треугольнику образованной системы треугольников или устанавливают, что она лежит вне указанной системы треугольников, для каждой узловой точки с координатами x_i и y_i, лежащей внутри системы треугольников, рассчитывают коэффициент соответствия по уравнению плоскости соответствующего треугольника, для узловых точек, лежащих вне системы треугольников, расчет коэффициента соответствия выполняют методом экстраполяции, для этого значения коэффициента соответствия в узловой точке принимают равными коэффициентам соответствия в ближайшей точке, лежащей на внешней границе системы треугольников, ранее рассчитанные приземные концентрации загрязняющих веществ в узловых точках заданной регулярной сетки умножают на полученные коэффициенты соответствия с получением уточненной концентрации загрязняющих веществ в узловых точках сетки, затем строят карту пространственного распределения уточненных концентраций загрязняющих веществ, по которой количественно оценивают уровень загрязнения атмосферного воздуха на исследуемой территории.