Оценка загрязнения атмосферного воздуха пылью по данным снегосъёмки на основе реконструкции полей выпадений

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

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

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

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

• ^* Работа выполнена при поддержке Программы РАН 4.9 и отраслевой научно-исследовательской программы «Гигиеническое обоснование минимизации рисков для здоровья населения России» (на 2011–2015 гг.) программы Роспотребнадзора по гигиенической оценке и минимизации риска для здоровья населения.

длительного загрязнения территорий города, и в этом плане перспективно использование данных мониторинга загрязнения снежного покрова [2].

Сопряжёнными исследованиями загрязнения пылью атмосферного воздуха и снежного покрова на стационарных постах Росгидромета ряда крупных городов юга Западной Сибири, включая Новосибирск, Кемерово, Барнаул, Томск, установлены коэффициенты пересчёта уровней загрязнения. Это даёт возможность по результатам изучения снежного покрова проводить ориентировочную гигиеническую оценку загрязнения воздушного бассейна города [4].

Геохимическими исследованиями выявлены количественные закономерности выпадений пыли в окрестностях промышленных площадок, что фиксируется в виде аномалий в снежном покрове, депонирующем загрязнение, и легкодоступном для изучения по любой заранее заданной сети точек отбора проб [2]. Учёт же дополнительной вполне доступной информации о параметрах источников, характеристиках дисперсного состава выбросов пыли, текущих метеорологических условиях позволяет оптимизировать систему точек отбора проб снега, повысить точность восстановления полей выпадений, оценить по данным внешнего мониторинга суммарные выбросы пыли в атмосферу [5].

Численная реконструкция поля длительных выпадений пыли. При расчете средней концентрации в приземном слое атмосферы определяющее значение имеют часто встречающиеся метеорологические условия. К ним относятся так называемые нормальные метеоусловия, для которых применима степенная аппроксимация скорости ветра и коэффициента вертикального турбулентного обмена [1]. Использование этих аппроксимаций, асимптотик полуэм-пирического уравнения турбулентной диффузии и свойств статистических характеристик распределения скорости ветра и вертикального турбулентного обмена в приземном слое атмосферы позволяет выразить плотность выпадений полидисперс-ной примеси за длительный промежуток времени в виде следующей регрессионной зависимости [5]

*

а ( r , ф , е > = 0 1 G ( r , 6 2 , 6 3 ) P ( ф+ 180 ° ), (1)

G ( r , θ 2 , θ 3 ) =

0->                                     to fw2exp(-0,w)(c) ,

™ I-I dto, 0   Г(1 + to)   v r / где r , ϕ – полярные координаты, P(ϕ) – приземная роза ветров; Г(m) – гамма-функция Эйлера, c – константа, зависящая от высоты источника, θ1,θ2,θ3 – неизвестные параметры, определяемые по данным наблюдений.

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

o(r,Ф,—) = 51 r52exp| -s^ IP(Ф+180°), (2) v r J где неизвестные параметры s1,s2,s3 также оцениваются по данным наблюдений.

Параметр s1 линейно зависит от ин тенсивности источника, а параметры s2,s3 – от характеристик дисперсного состава пыли.

Повторяемость направлений ветра P ( ϕ ) обычно задаётся в табличном виде по 8 или 16 румбам. Для её непрерывного описания удобно использовать между румбами следующую линейную интерполяцию по углу φ.

ϕ ∈

P_i ( ϕ ) = p_i + ^{pi + 1 - pi} π N

π i π ( i + 1) N ^, N

где pi – повторяемость i -го направления ветра, i=1,...,N.

Замечание 2. При отсутствии надёжной информации о розе ветров для данной точки местности или переменной интенсивности эмиссии источника пыли восстановление поля плотности осадка можно выполнить поэтапно. При фиксированном ϕ ₀ , т.е. на радиально расположенном маршруте пробоотбора, провести оценивание величин 0 1 ( ф ₀) = 0 1 P ( ф ₀ + 180 ° ), 0 2 , 0 3 . Затем, учитывая, что параметры θ ₂, θ ₃ практически от φ не зависят, провести оценивание величины θ 1 ( ϕ ) для других углов φ.

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

Qсум = ∫∫σ(ξ,η)dξdη, (4) S где S – территория вокруг источника, на которую происходит выпадение пыли, σ(ξ,η) – плотность выпадений пыли, представленная в декартовых координатах.

Если область S является кольцеобразной относительно источника, то с учётом (1) и замечания 2 соотношение (4) представляется в более удобном виде

2 π R 2

Q сум = ∫ θ 1 ( ϕ ) d ϕ ⋅ ∫ G ( r , θ 2 , θ 3 ) rdr . (5) 0 R ₁

Соотношение (5) позволяет существенно упростить вычисление суммарных выпадений пыли в различных областях и оптимизировать количество точек отбора проб. Например, с учётом интерполяционной формулы (3) в случае 8-румбового задания частей света и без использования розы ветров, оценку величины Q _сум можно провести по 10 опорным точкам наблюдений.

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

Город Искитим расположен в юговосточной части Новосибирской области в 55 км от областного центра – города Новосибирска. Искитимский цементный завод находится в северной части города на одной промплощадке. С северной и восточной сторон промплощадка предприятия примыкает к р. Бердь. С южной и западной сторон от промплощадки расположен жилой сектор. Ближайшее расстояние до жилых кварталов составляет 30–50 м.

Материалы и методы. Объектами исследования служили выбросы загрязняющих веществ от стационарных источников ОАО «Искитимцемент», снежный покров на территории г. Искитим и за его пределами. Материалами исследований являлись отчёты ОАО «Искитимцемент» о выбросах в атмосферу загрязняющих веществ стационарными источниками в период 2012 и 2013 г., результаты визуального обследования и физико-химического анализа состава проб снеговой воды. Маршруты отбора проб располагались по 8 румбам относительно основных источников выброса неорганической пыли – двух близко расположенных 80-метровых труб. Точки наблюдений находились в диапазоне расстояний от 0,4 до 3 км.

Отбор проб снега проводили с помощью пластмассовой трубы диаметром 10 см. В каждой точке отбора осуществлялась выемка от 2 до 10 кернов снега. Пробы снега таяли при комнатной температуре, талая вода фильтровалась через фильтр (синяя лента), измерялось значение рН. В пределах 1,5 км от основных источников величина рН варьировалась от 9 до 12. Исследова- ние химического состава снеговой воды и выделенного осадка проводили в аккредитованных лабораториях Новосибирского НИИ гигиены, Института неорганической химии СО РАН. Статистическая обработка и математическое моделирование выполнялись в Институте вычислительной математики и математической геофизики СО РАН.

Результаты и их обсуждение. Отбор проб снега был проведён по 8 радиальным маршрутам и более чем в 40 точках местности. Это позволило осуществить детальный численный анализ процессов выпадений пыли от основных источников предприятия, установить количественные закономерности содержания осадка в снеге на различных направлениях выноса. На рис . 1–3 представлены результаты численного восстановления по данным наблюдений поля плотности осадка неорганической пыли.

Рис. 1. Восстановленные на основе зависимости (1) выпадения неорганической пыли в северовосточном ( а ) и северном ( б ) направлении от цементного завода; ○, • – опорные и контрольные точки измерений

б

Из анализа рис. 1 вытекает, что согласие расчётов с данными наблюдений в контрольных точках вполне удовлетворительное. Максимум выпадений неорганической пыли находится примерно в 450 метрах от основных источников, что указывает на достаточно разнородный дисперсный состав оседающих частиц. Вынос пыли в зимнее время в северо-западном направлении является преобладающим и обусловлен, по-видимому, орографическими особенностями местности.

Рис. 2. Измеренные (■) и восстановленные (□) на основе зависимости (1) выпадения неорганической пыли (г/м2) в точках отбора проб снега

Полученные закономерности за рассматриваемый зимний период позволили определить величину суммарного выпадения неорганической пыли от выбросов в атмосферу на разном расстоянии от цементного производства. Суммарное содержание неорганической пыли в снеге на территории в радиусе 1 км от основных источников выброса ОАО «Искитимцемент» составило 626 тонн, в радиусе 2 км – 875 тонн, в радиусе 3 км – 942 тонны, в радиусе 4 км – 969 тонн. Полученные оценки суммарных выпадений существенно расходятся с данными инвентаризации валовых выбросов пыли цементным заводом, проведённой в 2012 г. Согласно этим данным, валовый выброс пыли должен был составить в зимнем сезоне 2012/2013 гг. порядка 100 тонн.

Рис. 3. Восстановленное поле плотности выпадений неорганической пыли (г/м2) в окрестностях Искитимского цементного завода по данным маршрутных снегосъёмок в конце зимнего сезона 2012/2013 гг.

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