CC BY
155
УДК 504.3.054
Ю. А. Тунакова, Р. А. Шагидуллина, И. Г. Григорьева УЧЕТ МЕТЕОРОЛОГИЧЕСКИХ ФАКТОРОВ, ФОРМИРУЮЩИХ УРОВЕНЬ ЗАГРЯЗНЕНИЯ ПРИЗЕМНОГО СЛОЯ АТМОСФЕРЫ В ЗОНАХ ДЕЙСТВИЯ ПОЛИМЕРНЫХ ПРОИЗВОДСТВ
Ключевые слова: метеорологические условия, уровень загрязнения, атмосферный воздух.
Рассматриваются способы учета вклада метеорологических условий в формирование уровня приземного слоя атмосферного воздуха.
Key words: meteorological conditions, level of pollution, atmospheric air.
Ways of the account of the contribution of meteorological conditions in formation of a level of a ground layer of atmospheric air are considered.
Введение
В зонах действия полимерных производств приземная концентрация атмосферных примесей определяется воздействием множества источников выбросов, что приводит к высокой динамичности содержания и соотношения химических элементов в воздухе [1].
Известно, что уровень загрязнения атмосферного воздуха промышленно развитых территорий формируется под влиянием ряда техногенных, метеорологических и физико-географических факторов. Поскольку техногенные и физико-географические факторы характеризуются малой изменчивостью, в основном динамичность поля приземных концентраций примесей определяется набором действующих
метеорологических факторов.
Если туманы, высокая влажность, штили, инверсионное распределение температуры создают условия для накопления примесей в атмосфере, то высокие скорости ветра и сильные осадки, вымывающие примеси, способствуют ее очищению. Большую опасность для состояния атмосферы представляют приземные температурные инверсии. Они затрудняют обмен воздуха по вертикали и, следовательно, уменьшают самоочищающую способность атмосферы. Создаются условия, ограничивающие подъем выбросов и способствующие накоплению примесей в приземном слое атмосферы [2].
Результаты исследования
Анализ материалов наблюдений в городах РФ с учетом тесноты и реального вида связей между метеорологическими факторами и концентрациями примесей в воздухе показал, что к определяющим метеопараметрам относятся направление переноса примесей, скорость их переноса, атмосферная устойчивость и связанная с ней степень вертикального перемешивания примесей, термическое состояние воздушной массы, от которого зависит начальный подъем выбросов, вымывание примесей осадками, их аккумуляция в туманах.
В целом возможными предикторами [3] могут быть: у0 - скорость ветра на высоте флюгера, м/с; d — направление ветра; v1 — скорость ветра на уровне изобарической поверхности АТ925, м/с; V —
средняя скорость в слое перемешивания, м/с; ДТ — разность температур на уровне земли и АТ925, °С; Т0
— температура воздуха у земли, °С; Нп — высота слоя перемешивания, м (иногда вместо Нп используется величина Нтд — слой термодинамического перемешивания, которая характеризует процесс вертикального обмена, определяемый как термической стратификацией, так и вертикальным сдвигом ветра; Н2—нижняя граница приподнятой инверсии, м (если инверсия приземная, то Н2 = 0); Н1 — мощность приземной инверсии, м (при отсутствии приземной инверсии Н1 = 0); ДР/Дп
— градиент приземного давления, гПа/градус меридиана);
Для количественной оценки
самоочищающей способности атмосферы используются различные показатели: потенциал загрязнения атмосферы (ПЗА), климатический (метеорологический) потенциал самоочищения атмосферы (КПСА) и коэффициент Селегей. При расчетах ПЗА учитываются метеорологические факторы, неблагоприятные для рассеивания примесей - туманы, штили и застои воздуха, инверсии температуры воздуха.
Величина ПЗА показывает во сколько раз средний уровень загрязнения атмосферного воздуха в конкретном районе, с определенной повторяемостью неблагоприятных для рассеивания примесей метеорологических условий, будет выше или ниже чем в условном районе, принятом за эталон. Так, для территорий, где преобладают низкие источники загрязнения, а так же высокие источники с холодными выбросами ПЗА может быть определен по формуле: ПЗА = 2,5ехр[(0,04/(72-71)2)-(0,471/(72-71))], где и 72 - аргументы вероятности Ф^), а Ф связано с соответствующими повторяемостями Р1 и Р2 соотношениями: Ф^) = 1-2Р1, Ф(72) = 1-2Р2. В свою очередь,
Р1^п) = РИн + Рсл - Рз - Рт;
Р2^>1^п) = Рз + Рт,
где Рш - повторяемость приземных инверсий, Рсл -повторяемость скорости ветра в интервале 0-1 м/с, Рз - повторяемость застоя воздуха (приземные инверсии со слабыми ветрами), Рт - повторяемость туманов [4].
В отличие от потенциала загрязнения атмосферы, климатический потенциал самоочищения
атмосферы и коэффициент Селегей учитывают повторяемость климатических процессов способствующих самоочищению атмосферы,- Так климатический потенциал самоочищения КМ определяется по формуле:
Рш + Р{
Км =-—
Ро + Ру
где РШ - число дней со штилем, скорость ветра V = 0 - 1 м/с; Р(- - число дней с влажностью > 80 %; РО - число дней с осадками > 0,5 мм; РV - число дней со скоростью ветра > 6 м/с [5],
В значительной степени схожим показателем, характеризующим климатический потенциал загрязнения атмосферы, является коэффициент Селегей [6,7]:
K =
Рш + Рт Pv + Ро
где РШ - число дней со штилем, скорость ветра V = 0 - 1 м/с; Рт - количество дней с туманом; Ро - количество дней с осадками > 0,5мм; Ру - число дней с сильным ветром > 6м/с,
Согласно [6,7], чем больше по абсолютной величине Кс, тех хуже условия рассеяния примеси в атмосфере, Если Кс < 1, то повторяемость процессов способствующих самоочищению атмосферы преобладает над повторяемостью процессов, способствующих накоплению вредных примесей в ней, Но недостатком используемых показателей является тот факт, что получаемая оценка по данным 1-2 метеопостов, функционирующих в городах, имеет общий фоновый характер без пространственно-временной дифференциации, что затрудняет ее использование в эколого-геохимических исследованиях, В этой связи на наш взгляд целесообразно учет влияния метеорологических параметров проводить по показателям загрязнения воздуха, получаемым по данным постов наблюдений за загрязнением атмосферы, которых в крупном городе функционирует 5 и более, Уровень загрязнения воздуха, изменяющийся на всей территории города под влиянием метеорологических условий, характеризуется обобщенными (интегральными) показателями, которые рассчитываются по материалам фактических наблюдений для всего дня и для отдельных частей дня, К таким обобщенным показателям относятся:
- осредненная по всему городу и по всем срокам наблюдений данного дня или части дня концентрация примеси в воздухе (0),
0 = Ч /?ср
где 0 — нормированная средняя концентрация в городе; дср — среднесезонная концентрация примесей в городе, Ч — средняя по городу концентрация примеси за данный день или часть дня
Ч = N £ Ч,
где ч — концентрация примесей на /-м посту; N — количество измерений по одной примеси и всем постам за все сроки наблюдений данного дня; зави-
а,
сящие от времени коэффициенты (а,) при членах разложения данных наблюдений за концентрациями примесей на стационарных постах в городе по естественным ортогональным функциям,
N
=Еф, j-q j,
j=i
где фу — компоненты ,-й естественной функции, относящиеся j-му стационарному посту; qj — отклонение от среднесезонной концентрации примеси на j-м посту в заданный момент времени; N — количество стационарных постов в городе.
Коэффициент при первом члене разложения а1 приближенно характеризует процесс одновременного роста или снижения концентраций примесей в воздухе в целом по городу. Он содержит наибольшую информацию об изменчивости концентраций. Коэффициент а2 чаще всего характеризует процесс переноса загрязненного воздуха из одной части города в другую. Последующие коэффициенты описывают другие особенности загрязнения воздуха в городе. Вследствие ортогональности естественных функций коэффициенты разложения а, статистически независимы, и прогноз каждого из них можно осуществлять независимо друг от друга. Установлено, что для удовлетворительного описания поля концентраций примеси достаточно использовать первые 2—3 члена разложения. Таким образом, при использовании а,- в качестве обобщенных показателей наряду с интегральным загрязнением воздуха в городе описываются другие, более детальные особенности распространения примесей на городской территории. Коэффициент aj близок к параметру Р и к средней по городу концентрации примесей.
Параметр Р, который рассчитывается по формуле:
Р = m/n,
где n — общее количество наблюдений за концентрацией примесей в городе в течение одного дня на всех стационарных постах; m — количество наблюдений в течение этого же дня с концентрациями q, которые превышают среднесезонное значение q^ более чем в 1,5 раза (q > 1^ср).
Таким образом, параметр Р представляет собой отношение количества существенно повышенных концентраций (относительно среднего значения) к общему числу измерений в течение дня. Для получения ежедневных значений параметра Р предварительно рассчитываются среднесезонные значения концентраций примесей для каждого стационарного поста. Расчеты средних значений концентраций проводятся отдельно для каждого года. По всему используемому ряду наблюдений отмечаются единичные концентрации, которые превышают j,5q^. Для каждого дня определяются величины m, n и параметр Р.
Параметры Р, Q и aj являются относительными характеристиками и не зависят от среднего уровня загрязнения воздуха. Следовательно, их значения в основном определяются метеорологически-
ми условиями. В связи с этим оказывается возможным использовать для анализа значения указанных параметров за несколько лет вместе, даже если в течение этого периода наблюдались изменения суммарного количества выбросов и уровня концентраций.
Параметр Р рассчитывается для совокупности примесей, а при достаточном количестве измерений в течение дня — для отдельных примесей, а параметры Q и а!—только для отдельных примесей. Таким образом, параметр Р позволяет характеризовать общее состояние загрязнения воздуха в городе и в то же время при наличии достаточного количества измерений он может быть показателем городского загрязнения отдельными примесями. Нами проводился корреляционный анализ, который показал наличие тесных корреляционных связей между параметром Р, рассчитанным по совокупности примесей и легко измеряемыми метеопараметрами: направление ветра (п), скорость ветра (у), температура (Т), давление (р), влажность воздуха ; Данные метеопараметры были выбраны в качестве значимых предикторов. Исследования проводились в зоне действия полимерных производств - на территории Нижнекамского промышленного узла. Р = -0,2207-0,1028 хТ-0,0049хТ2; Р=0,029хр-21,4; Рч =0,0558+0,2401 ху-0,043 ху2; Рч =0,1085+0,0011 хп-2,0326х10-6хп2; Р = -86,6768+2,1845 ху-0,0137ху2; коэффициент множественной корреляции Я=0,64; коэффициент детерминации Я2=0,41; критерий Фишера Б=21,1; уровень значимости р=0,00015.
При исследовании многофакторных процессов возникает необходимость установления зависимости между несколькими переменными. Методом множественной регрессии все предикторы были сведены в несколько моделей, линейно характеризующих их взаимосвязь с параметром Р в пределах одного из сезонов (зима, весна, лето, осень). Так, для весны рассчитано следующее уравнение: Р=8,33-0.0014хТ-0,0098хр-0,058ху+0,00025хп+0,0103х;Г
Параметры модели: коэффициент множественной корреляции Я=0.67, коэффициент детерминации Я2=0,45, критерий Фишера Б=18,5, уровень значимости р=0.0034.
Для лета, уравнение регрессии выглядит следующим образом:
Р=0,002852х; + 0,000033 хп-0,02462ху+0,000473 хр-0,006634хТ - 0,159924;
Параметры модели: коэффициент множественной корреляции Я=0,73, коэффициент детерминации Я2=0,53, критерий Фишера Б=24,0, уровень значимости р=0,011.
Соответствующее регрессионное уравнение для осени:
Р=2,577650+0.004827хТ - 0,002603 хр-
0.03777.у+0.0008 хп-0,006604х;
Параметры модели: коэффициент множественной корреляции Я=0,82, коэффициент детерминации Я2=0,67, критерий Фишера Б=14,0, уровень значимости р=0,0083.
Таким образом, использование параметра Р для характеристики уровня загрязнения приземного слоя атмосферы позволяет учитывать влияние метеорологических факторов на рассеивание примесей.
Литература
1. Тунакова Ю.А., Новикова С.В., Шагидуллина Р.А., Шмакова Ю.А. Вестник Казанского технологического университета, 12, 71-74 (2012).
2. Методические указания по прогнозированию загрязнения воздуха в городах с учетом метеорологических условий. РД 52.04.78-86. - М.: Госкомгидромет. СССР, 1986.
3. Безуглая Э.Ю. Метеорологический потенциал и климатические особенности загрязнения воздуха городов. Л.: Гидрометеоиздат, 1980. -184с.
4. Безуглая Э.Ю. К определению потенциала загрязнения воздуха. Л.: Труды ГГО, 1968, вып.234, С.69-80.
5. Линевич Н.Л., Сорокина Л.П. География и природные ресурсы. 4. (1992).
6. Селегей Т.С., Юрченко И.П. География и природные ресурсы. 2. .132-137 (1990).
7. Селегей Т.С. География и природные ресурсы. 1, 44-48 (1994).
© Ю. А. Тунакова - д-р хим. наук, проф. каф. технологии пластических масс КНИТУ, ]иНарго1"@таП.ги; Р. А. Шагидуллина - канд. хим. наук, нач. управления государственной экологической экспертизы и нормирования воздействия на окружающую среду Мин-ва экологии и природных ресурсов РТ, Raisa.Shagidullina@tatar.ru; И. Г. Григорьева - доц. каф. общей химии и экологии КНИТУ-КАИ.
