На правах рукописи
КОСТАРЕВА Татьяна Викторовна
ПРОГНОЗ ЭКСТРЕМАЛЬНО ВЫСОКИХ УРОВНЕЙ ЗАГРЯЗНЕНИЯ ВОЗДУХА
ДЛЯ КРУПНЫХ ПРОМЫШЛЕННЫХ РЕГИОНОВ
(НА ПРИМЕРЕ УРАЛЬСКОГО РЕГИОНА)
Специальность 25.00.30 Метеорология,
климатология, агрометеорология
АВТОРЕФЕРАТ
диссертации на соискание ученой степени кандидата географических наук
Санкт – Петербург 2011 2
Работа выполнена в Главной геофизической обсерватории им. А. И. Воейкова
Научный руководитель: доктор географических наук Сонькин Лев Рахмилович
Официальные оппоненты: доктор физико-математических наук Е. Л. Генихович кандидат географических наук Н. С. Буренин
Ведущая организация: Российский государственный гидрометеорологический университет
Защита диссертации состоится « 15 » февраля 2012 г. в 14 часов на заседании диссертационного совета Д 327.005.01 по защите докторских и кандидатских диссертаций в ФГБУ «ГГО» по адресу: 194021, Санкт-Петербург, ул. Карбышева, д.
С диссертацией можно ознакомиться в библиотеке федерального государственного бюджетного учреждения «Главная геофизическая обсерватория им. А. И. Воейкова»
Автореферат разослан « 13 » января 2012 г.
Ученый секретарь совета по защите докторских и кандидатских диссертаций, доктор географических наук А. В. Мещерская
ОБЩАЯ ХАРАКТЕРИСТИКА РАБОТЫ
Актуальность темы Одной из острых экологических проблем в крупных промышленных регионах остается проблема загрязнения воздуха. Для её решения и соответственно защиты здоровья населения в периоды неблагоприятных метеорологических условий большой интерес представляет задача прогнозирования и предотвращения экстремально высоких уровней загрязнения воздуха в городах (ЭВУЗВ), с которыми связаны наиболее неблагоприятные для населения и окружающей среды эпизоды. Имеются в виду опасные случаи высоких концентраций загрязняющих веществ, которые наблюдаются редко и могут быть приравнены к стихийным бедствиям, таким как наводнения, цунами, ураганы и др.
Указанная задача представляет особый интерес и для Уральского региона, где расположено много городов с развитой промышленностью, много предприятий, выбрасывающих в атмосферу вредные вещества. В наиболее неблагоприятные дни загрязнение воздуха может достигать катастрофического уровня. Особенности городов региона заключаются в ориентации их промышленности на использование металла, нефти и газа. Уральский регион относится к регионам с высоким уровнем загрязнения воздуха.
Больше половины городов региона ежегодно фигурируют в Приоритетном списке с очень высоким уровнем загрязнения воздуха.
Цель и задачи исследования Целью данной диссертации является исследование закономерностей формирования ЭВУЗВ на территории крупных промышленных регионов, а также разработка эффективной схемы прогноза ЭВУЗВ с целью защиты атмосферы при неблагоприятных метеорологических условиях (НМУ) в масштабе региона и в каждом из городов данного региона, в частности и при отсутствии специально разработанных схем прогноза ЭВУЗВ.
В соответствии с этой целью в диссертации были поставлены следующие задачи:
1. Выявить и исследовать условия формирования ЭВУЗВ в городах Уральского региона 2. В дни с экстремально высоким уровнем загрязнения воздуха определить связь между значениями обобщенного интегрального показателя загрязнения воздуха (параметра Р) и значениями концентраций примесей.
3. Исследовать особенности формирования ЭВУЗВ в целом по региону 4. Разработать статистические схемы прогноза ЭВУЗВ с использованием метода графической регрессии и метода множественной линейной регрессии с предварительным исключением нелинейности связей для городов Уральского региона 5. Оценить эффективность схем прогноза ЭВУЗВ на зависимом и независимом материале.
1. Создана эффективная региональная система прогноза ЭВУЗВ в крупном промышленном регионе (на примере Уральского региона), а также в каждом из городов данного региона, с целью предотвращения высоких уровней загрязнения воздуха, формирующихся в периоды НМУ.
2. Показана возможность и целесообразность использования одного и того же синоптического предиктора Sn для промышленных городов в пределах региона, так как по материалам ряда городов установлена тенденция к одновременному формированию ЭВУЗВ в городах региона 3. Показано, что в дни с значениями параметра Р, соответствующими ЭВУЗВ, в городах региона реально имеют место опасные концентрации примесей в воздухе, значительно превышающие в ряде случаев предельно допустимые концентрации (ПДК).
4. Впервые по Уральскому региону проведена типизация синоптических процессов, способствующих накоплению вредных примесей в атмосфере с учетом географического происхождения, траектории движения барических образований, скорости ведущего потока по картам АТ 700 гПа или АТ 500 гПа, термического поля над регионом, а также длительности нахождения барического образования над рассматриваемым районом.
5. Даны рекомендации по прогнозу ЭВУЗВ с использованием единого синоптического предиктора в городах, для которых имеющихся материалов недостаточно для разработки методов прогноза загрязняющих веществ.
6. Получены общие по региону характеристики загрязнения воздуха с использованием метода разложения по естественным ортогональным функциям для оценки ожидаемого общего состояния загрязнения воздуха в регионе по материалам пяти Уральских городов.
Практическая ценность и реализация результатов работы Разработанные схемы прогноза ЭВУЗВ позволяют заблаговременно предсказывать в городах региона подавляющее большинство (79% - 92%) случаев экстремально высокого уровня загрязнения воздуха. Авторские испытания предложенного метода, проведенные по независимым материалам наблюдений пяти городов (Пермь, Екатеринбург, Соликамск, Березники, Губаха), подтвердили его эффективность. Схемы прогноза ЭВУЗВ в настоящее время используются в оперативной работе в городах Уральского региона. Они могут быть рекомендованы для оперативного использования во всех крупных промышленных регионах. Разработанные схемы прогноза ЭВУЗВ для городов Уральского региона размещены на сайте ГМЦ России в Методическом кабинете.
Основные результаты диссертационной работы вошли в Методические рекомендации, подготовленные для апробации в региональных управлениях Росгидромета, докладывались на международной научной конференции, проводившейся в Уральском отделении РАН (г. Екатеринбург, 2011г.), научно-методических курсах повышения квалификации для руководящих работников и специалистов Росгидромета "Современные задачи мониторинга загрязнения атмосферы" в Санкт-Петербурге (2006 г), на Ученом Совете ГГО, семинарах отдела исследования и мониторинга загрязнения атмосферы в ГГО; работа участвовала в региональном конкурсе исследовательских работ «Экология. Человек года», где была отмечена благодарностью; работа участвовала в конкурсе Росгидромета на лучшие научно-исследовательские и опытно-конструкторские работы, признана лучшей региональной работой, присуждена ведомственная премия.
Диссертация состоит из введения, трех глав, заключения, приложения.
Она содержит 234 страницы, включая 30 рисунков, 19 таблиц и библиографию из наименований.
СОДЕРЖАНИЕ РАБОТЫ
Во введении обосновывается актуальность темы, дается общая характеристика работы, формулируются цели и основные задачи исследования, отмечается его научная новизна и практическая ценность.
Первая глава посвящена обзору экологического состояния загрязнения воздушного бассейна в городах Уральского региона, рассматриваются основные закономерности распространения примесей в атмосфере городов, а также влияние метеорологических условий на создание высоких уровней загрязнения воздуха. Рассматриваются вопросы теоретических и экспериментальных исследований процесса распространения примесей в атмосфере городов, а также работы по краткосрочному прогнозированию загрязнения атмосферы, выполненные как в нашей стране, так и в ряде зарубежных стран.
В разделе 1.1 описано современное состояние загрязнения воздуха в городах Уральского региона. За последние годы по данным «Ежегодников состояния загрязнения атмосферы в городах на территории России» количество городов с высоким уровнем загрязнения атмосферы снизилось со 145 до 136. Количество городов Приоритетного списка с высоким уровнем загрязнения атмосферы с тридцати семи в 2007 г. снизилось до тридцати в 2008 г. Несмотря на это, в городах Уральского региона уровень загрязнения воздуха остается высоким. В городах Уральского региона наибольший вклад в выбросы загрязняющих веществ вносят предприятия химической, нефтехимической и топливной промышленности. За пятилетний период 2004–2008 гг. возросли средние за год концентрации в городе Перми взвешенных веществ, фторида водорода и бенз(а)пирена, в городе Екатеринбурге — аммиака, формальдегида и этилбензола, в городах Березники, Губаха, Соликамск — бензола и этилбензола. Среднегодовые концентрации отдельных примесей превышают ПДК во всех городах региона, Так в 2008 год в Перми средняя за год концентрация формальдегида составила 3 ПДК, фторида водорода — 1,4 ПДК, бенз(а)пирена — 2,6 ПДК; в г. Екатеринбурге средняя за год концентрация формальдегида составила 5 ПДК, диоксида азота — 2 ПДК, бенз(а)пирена — 1,6 ПДК.
Максимальные значения концентраций примесей, превышающие в 10 раз максимальные разовые ПДК (ПДКмр), в 2007 и 2008 годах наблюдались в большинстве городов Уральского региона. Число жителей, многократно испытывающих воздействие такого высокого уровня загрязнения воздуха, составляет в регионе более 5 млн. человек. В 2007 году города Пермь, Соликамск, Екатеринбург входили в Приоритетный список городов с очень высоким уровнем загрязнения воздуха в России, в 2008 году — город Екатеринбург.
В разделе 1.2 кратко изложены результаты теоретических исследований, выполненных в ГГО под руководством профессора М. Е. Берлянда. Приводятся описания основных закономерностей изменения характеристик загрязнения воздуха от источников выбросов различных типов под влиянием метеорологических условий. В этом же разделе излагаются результаты анализа влияния метеорологических условий на уровень загрязнения воздуха.
Изменения уровня загрязнения воздуха по городу в целом определяются множеством факторов, доминирующее значение среди которых принадлежит метеорологическим условиям: скорости и направлению ветра, устойчивости атмосферы, слоям инверсии и изотермии, температуре воздуха, облачности, туманам, инерционным факторам и др.
Результаты, представленные в данной главе, получены на основании выполненного автором анализа фактических материалов по городам Уральского региона с учетом имеющихся выводов ряда авторов в других городов России. Так, в городе Пермь при югозападном ветре жилая зона города, расположенная северо-восточнее промышленного комплекса, загрязняется более интенсивно, чем при других направлениях ветра. Для города Березники неблагоприятное направление ветра – северное и северо-западное. В целом в городах Уральского региона зависимость загрязнения воздуха от направления ветра соответствует характеру расположения источников и эффекту наложения выбросов.
В зависимости от скорости ветра, как и в ряде других городов, отмечаются два максимума загрязнения воздуха - при штиле и при скорости ветра 4 - 6 м/с. Первый максимум обусловлен действием низких источников и автотранспортом, создающих общий фон загрязнения воздуха. Второй максимум обусловлен мощными выбросами высоких источников, для которых опасная скорость ветра равна 4 - 6 м/с.
Связь между температурой воздуха и загрязнением приземного слоя атмосферы в городах обусловливается несколькими обстоятельствами. В городах Уральского региона отмечается тенденция к повышению уровня загрязнения воздуха при росте температуры воздуха, особенно в летний период (Рис.1). Это проявляется при рассмотрении температуры в слое ОТ-500/1000. На рис.1 показана зависимость параметра Р от температуры в слое 500/1000 летом в городе Соликамске.
Рис. 1 Зависимость параметра Р от температуры в слое 500/1000 летом, Зимой высокий уровень загрязнения воздуха в городах региона имеет место в первую очередь в условиях застоя воздуха. Обработка фактического материала наблюдений в нескольких городах Уральского региона показала, что повторяемость случаев высокого загрязнения воздуха при застоях воздуха в среднем составляет 74%.
Осадки оказывают значительное влияние на уровень загрязнения воздуха в городах.
По материалам наблюдений в Перми была изучена зависимость уровня загрязнения воздуха от осадков. В таблице 1 показаны значения параметра Р в дни с осадками и без осадков в феврале и июле.
В разделе 1.3 приводится обзор работ по краткосрочному прогнозированию загрязнения атмосферы.
В настоящее время разработаны разные методы прогноза загрязнения воздуха.
При этом отмечаются два подхода: детерминированный и статистический.
В основу детерминированных методов положены разработки теории атмосферной диффузии на основе математического описания распространения примесей с помощью уравнения турбулентной диффузии. Подробная библиография по данному вопросу приведена в монографии М. Е. Берлянда.
Статистические методы прогноза загрязнения воздуха в городах базируются на эмпирико-статистическом анализе материалов наблюдений за концентрациями примесей в воздухе и сопутствующими метеорологическими и синоптическими условиями. Методы краткосрочного прогнозирования загрязнения атмосферы получили в России широкое развитие. В работах группы специалистов Главной геофизической обсерватории, возглавляемой Л.Р. Сонькиным, разработана методология и успешно внедрены в оперативную практику Росгидромета методы прогноза от одиночных источников и интегральных показателей загрязнения по городу в целом. Такие методы активно используются и за рубежом. В основу этих разработок положен учет реального вида связей между загрязнением воздуха и метеорологическими и синоптическими условиями.
В качестве прогнозируемой величины обычно используется параметр Р, который является вероятностной характеристикой превышения концентраций, измеренных в течение дня на всех постах наблюдения в городе, по отношению к некоторому среднему уровню загрязнения. Он определяется по формуле: Р= m/n, где n – общее количество наблюдений за концентрацией примесей в городе в течение одного дня на всех стационарных постах;
m – количество наблюдений в течение этого дня с концентрациями q, которые превышают среднесезонные значение qср более чем в 1.5 раза (q>1.5qср ).
В разделе 1.4 рассматривается проблема прогноза экстремально высокого уровня загрязнения воздуха (ЭВУЗВ). Под ЭВУЗВ по городу в целом принимаются редко встречающиеся случаи (1-2%) с наибольшими значениями обобщенного показателя загрязнения воздуха. Параметр Р в этих случаях обычно равен или превышает 0,5.
Большую роль в формировании высоких уровней загрязнения воздуха играет инерционный фактор. Для установления синоптических процессов формирования ЭВУЗВ в каждом городе в отдельности по многолетнему ряду наблюдений за предшествующие годы отдельно по сезонам выбираются случаи, попадающие в 2% дней с наибольшими значениями параметра Р или другого интегрального показателя. Рассматриваются синоптические ситуации, соответствующие этим случаям, устанавливаются характерные ситуации. Затем такие ситуации отбираются по всему материалу, который используется при выполнении разработок. Если при заданной ситуации наряду с ЭВУЗВ имеют место случаи с относительно низким уровнем загрязнения, то ситуация уточняется. В первую очередь изучаются различия ситуации в случаях с экстремально высоким уровнем загрязнения воздуха и в случаях с низким уровнем загрязнения воздуха. С учетом характера синоптических ситуаций, особенностей развития синоптических процессов, а также инерционного фактора в ГГО им. Воейкова был сформулирован ряд правил для прогноза ЭВУЗВ. Правила относятся как к отдельному городу, для которого составляется прогноз, так и для всего региона, расположенного в зоне данной неблагоприятной ситуации. Сформулированные прогностические правила для прогноза ЭВУЗВ общие, поэтому в каждом регионе, с учетом его особенностей, необходимо их уточнять.
Во второй главе рассматриваются условия формирования ЭВУЗВ в городах Уральского региона, приводятся сведения о данных, использованных при разработке схем прогноза.
В разделе 2.1 рассмотрены результаты анализа экстремальных случаев загрязнения воздуха в городах Уральского региона. Предварительно для каждого из городов региона были рассчитаны свои критерии ЭВУЗВ. Для этого по всему ряду наблюдений, с учетом 2% - ной повторяемости, были выбраны случаи с наибольшими значениями параметра Р.
В г. Березники к ЭВУЗВ относятся случаи со значением параметра Р 0,42, как для летнего, так и зимнего периодов. В г. Соликамске для летнего периода это дни со значением параметра Р 0,38, для зимнего - Р 0,40. В г. Екатеринбурге для летнего периода это дни со значением параметра Р 0,40, для зимнего - Р 0,60. В г. Перми для летнего периода это дни со значением параметра Р 0,42, для зимнего - Р 0,45. В г.
Губаха для летнего периода это дни со значением параметра Р 0,45, для зимнего периода - Р 0,50. Далее в процессе анализа материалов пяти уральских городов достаточно определенно выявилось, что в дни с большими значениями параметра Р (в дни с ЭВУЗВ) в городах возникают опасные условия. Так, во всех городах региона в дни с экстремально высоким уровнем загрязнения воздуха на всех постах наблюдения за загрязнением атмосферного воздуха отмечались концентрации различных примесей, превышающие ПДК от 3 до 10 раз. Превышение концентраций примесей над 1,5q среднесезонными (используется при расчете ежедневного значения обобщенного интегрального показателя загрязнения воздуха) в летний период составило от 3 до 10,5 раз, в зимний - в 2-13 раз. В отдельные дни превышение концентраций вредных примесей над 1,5q среднесезонными было в 13-65 раз.
По данным для этих городов наибольшие концентрации диоксида азота, диоксида серы, оксида углерода превышают в дни с ЭВУЗВ 1,5qср в 7 – 9 раз и более, сероводорода, формальдегида, соляной кислоты, оксида азота – в 30 – 50 раз. Концентрации вредных веществ в воздухе в такие дни значительно превышают ПДК. В Перми и Соликамске при ЭВУЗВ концентрации сероводорода превышают ПДК в 8 – 10 раз. Концентрации других измеряемых примесей во всех городах региона превышают ПДК в такие дни в 3 – 5 раз.
Повторяемость случаев превышения ПДК по градациям для пяти городов Уральского региона представлена в таблице 2. В этой же таблице для сравнения приведена повторяемость случаев превышения ПДК по градациям для этих же городов в дни со значением параметра Р < 0,10, что соответствует пониженному уровню загрязнения воздуха.
Повторяемость (%) случаев превышения ПДК по градациям в дни с ЭВУЗВ (в числителе) и в дни с пониженным (Р < 0,10) уровнем загрязнения воздуха (в Анализ материалов показал, что во всех городах Уральского региона в дни с ЭВУЗВ одновременно на нескольких или на всех постах общегосударственной службы наблюдений и контроля над загрязнением атмосферы (ОГСНКА) сразу по нескольким примесям превышение ПДК составляет в зимний период 90% - 100%, в летний период – 83% – 100%; превышение 3 ПДК, составляет в зимний период от 21% до 64%, в летний — от 23% до 50 %, в г. Губаха 17%, в Перми 8%;. превышение 5 ПДК составляет от 8% до 29%. В дни с пониженным уровнем загрязнения воздуха превышение ПДК отмечается лишь на отдельных постах по отдельным примесям. В зимний период оно составляет от 14% до 23%, в летний период – 21% – 36%, лишь в г. Березники – 46%; превышение ПДК составляет в зимний период от 1% до 2%, в летний - от 2% до 4 %; превышение ПДК в отдельных городах региона составляет 1% - 2%.
Таким образом, проведенный анализ убедительно продемонстрировал, что во всех без исключения отобранных случаях с ЭВУЗВ концентрации примесей в воздухе во всех городах региона превышают ПДК.
В разделе 2.2 приводится классификация синоптических условий формирования ЭВУЗВ в городах Уральского региона.
По Уральскому региону такая классификация проведена впервые. Анализ выполнялся по приземным и высотным картам за 00 часов всемирного скоординированного времени (ВСВ) за шестилетний период (2000-2005гг). Были прослежены и изучены траектории движения барических образований, вычислены скорости их движения, по картам АТ 500гПа или АТ 700гПа были определены местоположение и ориентация ложбин и гребней. Термические условия проанализированы по картам АТ 850гПа и ОТ 500/1000. Учитывались географическое происхождение барических образований, а также длительность их нахождения над рассматриваемым регионом.
При установлении критерия ЭВУЗВ за шестилетний период было проанализировано около 1800 случаев. Учитывались все 12 месяцев года. При разработке схем прогноза ЭВУЗВ использовались данные только по летним и зимним месяцам.
За зимний период по данным города Перми было рассмотрено и проанализировано 426 случаев. Отмечалось двадцать два случая ЭВУЗВ, из них девять случаев отмечалось, когда над районом располагался центр стационарного антициклона, восемь случаев в малоградиентной северо-западной части гребня Сибирского антициклона и пять случаев в малоградиентном поле.
За летний период было проанализировано 429 случаев, ЭВУЗВ отмечался двенадцать раз: шесть случаев - в малоградиентном поле и пять случаев – в центре антициклона.
В результате анализа были выявлены синоптические ситуации, при которых формируется экстремально высокий уровень загрязнения воздуха в городах региона. Они приведены на рисунках 2-5.
Рис. 3 Малоградиентная северо-западная часть гребня Сибирского антициклона Рис. 5 малоградиентная южная периферия арктического Оказалось, что при выявленных ситуациях наблюдались и пониженные уровни загрязнения воздуха. Для выявления особенностей возникновения ЭВУЗВ в одних случаях и отсутствия в других ситуации уточнялись. Были установлены различия в близких ситуациях при ЭВУЗВ и при пониженных уровнях.
В процессе выполненного уточнения, получены следующие выводы:
В летний период при антициклонах, смещающихся из северо-западных районов ЕТР на территорию Уральского региона и антициклонах, сформировавшихся при ультраполярных вторжениях в тыловых частях циклонов, ЭВУЗВ никогда не отмечается.
ЭВУЗВ в летний период отмечается только при антициклонах, смещающихся с югозападных районов ЕТР на территорию Уральского региона или при распространении гребня Азорского антициклона на территорию региона. На высотах обычно отмечается малоградиентное поле повышенного давления. Над районом в эти дни расположена периферия очага тепла или очаг тепла, температура в слое 500/1000 выше 7о С.
В зимний период ЭВУЗВ при синоптической ситуации «центр антициклона»
отмечается при любом типе антициклонов. Характерными особенностями наличия ЭВУЗВ в центре антициклона в зимний период являются:
1. Антициклон сохраняется над районом два дня и более (стационарный малоподвижный антициклон).
2. Параметр Р' имеет повышенное значение (Р0,20).
3. Ветер слабый (0-2м/с), как приземный, так и градиентный 4. Наличие приземных инверсий мощностью 500-1700м, интенсивностью более 14о С.
5. Температура в слое 500-1000м от -10оС до -13оС.
6. На карте ОТ 500/1000 отмечается очаг тепла или его периферия.
Экстремально высокий уровень загрязнения воздуха в зимний период никогда не встречается, если:
1. Центр антициклона располагается над регионом первые сутки или быстро смещается.
2. Значения исходного уровня загрязнения воздуха низкие ( Р'0,20 ) 3. На карте ОТ 500/1000 над районом располагается очаг холода или его периферия.
4. Температура в слое 500-1000м ниже -13оС.
Экстремально высокий уровень загрязнения воздуха (ЭВУЗВ) при синоптической ситуации «северо-западная часть гребня Сибирского антициклона» в зимний период отмечается при следующих условиях:
1. Ветер у земли юго-восточный или переменный, не более двух метров в секунду.
2. Градиентный ветер – юго-восточный, южный не более 6м/с.
3. Наличие приземных инверсий мощностью 800-1200м и интенсивностью 2-14оС.
4. Синоптическая ситуация сохраняется над районом двое суток и более.
5. Ведущие потоки над районом юго-восточные (1400 – 1600) – не более 5м/с.
6. Наличие над территорией Урала мощного высотного антициклона ЭВЗУВ при этой же синоптической ситуации никогда не отмечается, если:
1. Градиентный ветер более 6м/с 2. Отсутствует приземная инверсия или отмечается приподнятая инверсия с нижней границей 0,4-0,6км.
3. Скорость ведущего потока по картам АТ 700гПа и АТ 500гПа 10 м/с и более.
4. Территория региона находится в северной или северо-западной части высотного гребня.
Экстремально высокий уровень загрязнения воздуха при синоптической ситуации «малоградиентное поле» в зимний период отмечается, если одновременно выполняются три условия:
1. Над рассматриваемым районом в вышележащих слоях расположен очаг тепла или его периферия.
2. Скорость ветра у земли и на высотах не более 5 м/с 3. В утренние часы отмечается приземная инверсия.
В летний период ЭВУЗВ отмечается только при синоптической ситуации «малоградиентная южная периферия арктического антициклона». Для этой синоптической ситуации характерно наличие обширного, меридионально расположенного антициклона с центром в районе Обской Губы. Территория Уральского региона находится под влиянием его малоградиентного гребня. На картах барической топографии (АТ 700гПа и АТ 500гПа) над территорией региона расположена ось высотного гребня или центр антициклона. ЭВУЗВ отмечается при наличии очага тепла или его периферии над рассматриваемым районом. Температура в слое 500/1000 в этот период от 7,50 до 11,50С.
Классификация синоптических процессов четко выявила синоптические ситуации, при которых формируется экстремально высокий уровень загрязнения воздуха в регионе. Полученные результаты будут использованы при разработке схем прогноза ЭВУЗВ для городов Уральского региона.
В третьей главе формулируется основная задача работы, которая заключается в прогнозировании ЭВУЗВ в крупном промышленном регионе, а также приводятся синоптико-статистические схемы прогноза ЭВУЗВ в городах Уральского региона.
В разделе 3.1 представлены основные положения синоптико-статистических схем прогноза ЭВУЗВ и подход к решению задачи прогноза.
Решение задачи прогноза ЭВУЗВ включает два этапа. Первый из них касается анализа метеорологических, в первую очередь, синоптических условий формирования высоких уровней загрязнения воздуха и установления на основе этого анализа прогностических правил, приведенных во второй главе. Второй этап включает разработку специальных статистических схем для данного вида прогнозов.
Для разработки схем прогноза ЭВУЗВ использовался ряд наблюдений, в который включены были только дни с повышенным уровнем загрязнения воздуха. Исходя из этого в схемы прогноза для городов Пермь, Березники, Соликамск и Губаха были включены все случаи со значением параметра Р>0,20, для города Екатеринбурга для зимнего периода это дни со значением параметра Р0,30.
В разделе 3.2 излагается способ установления синоптического предиктора, с целью его использования в схемах прогноза ЭВУЗВ.
Для учета зависимости уровня загрязнения воздуха от синоптических условий в прогностических схемах вводится количественный синоптический предиктор. Одной из возможностей получения такого предиктора является кодировка типов синоптических ситуаций баллами. Для разработки синоптико-статистических схем прогноза ЭВУЗВ существующий способ установления синоптического предиктора Sn неприменим, вследствие малого числа случаев ЭВУЗВ в отдельном городе. Для оценки факторов, влияющих на формирование ЭВУЗВ, был применен новый подход с использованием системы баллов. Комплексная оценка степени опасности синоптической ситуации Sn складывалась из суммы баллов ее отдельных составляющих: Sn= Бтс+Бпс+Бпр+Бос (1), где Бтс – балл текущей ситуации, Бпс – балл предшествующей ситуации, Бпр – балл преобразования синоптических процессов, Бос - балл неблагоприятных особенностей синоптических процессов.
При установлении балла текущей ситуации, все синоптические ситуации были разделены на группы в соответствии с их влиянием на создание высоких уровней загрязнения в городе. Была учтена классификация синоптических процессов, изложенная выше, а также повторяемость высокого и экстремально высокого уровня загрязнения воздуха при определенной синоптической ситуации (Таблицы 3 и 4).
№группы №группы К 1 группе были отнесены наиболее неблагоприятные ситуации и оценены в 4 балла.
Зимой это стационарный антициклон (Аст), западная и северо-западная части гребня Сибирского антициклона при скорости ветра Vо 6м/с (ПАСЗ, ПАЗ), летом – стационарный антициклон (Аст), малоградиентная южная часть арктического гребня (Га) и малоградиентное поле повышенного давления (М).
Неблагоприятные ситуации были включены во 2 группу и оценены в 3 балла. В эту группу вошли ситуации 1 группы в первый день их действия. В зимний период дополнительно малоградиентное поле (М).
В 3 группу зимой включены следующие синоптические ситуации – южная, юговосточная периферии циклона (ПЦЮ, ПЦЮВ), восточная периферия антициклона (ПАВ), летом – все периферии антициклона (ПАЮ, ПАЮЗ, ПАВ, ПАЗ, ПАСЗ), кроме северной и северо-восточной. К 4 группе в зимний период отнесены такие синоптические ситуации западная и северо-западная периферии антициклона (ПАЗ, ПАСЗ с V о > 6м/с), а в летний период - малоградиентное поле пониженного давления (Мц).
В 5 группу зимой включены восточная и северная части циклона (ПЦВ, ПЦС), центр циклона (Ц), летом – ложбина (Л), теплый сектор циклона (Цтс) и северная периферия антициклона (ПАС).
Все остальные синоптические ситуации, не вошедшие в предыдущие группы, были отнесены к 6 группе. Группы 3 и 4 были оценены в 2 балла, группа 5 в 1 балл и группа 6 в 0 баллов.
Обычно характер связей между параметром Р и синоптическими ситуациями текущего и предшествующего дней сохраняется. Поэтому балл предшествующей ситуации решено оценивать так же, как и текущей, но с коэффициентом 0,5: Бпс = 0,5Бтс.
Формирование ЭВУЗВ связано с преобразованием синоптических процессов, с переходом ситуаций из одной группы в другую или с сохранением ситуации в данной группе. Вклад такого преобразования в создание ЭВУЗВ оценивался баллами Бпр.
Наиболее неблагоприятные преобразования: переход ситуации из 4 группы к 1 или группам; переход любой неблагоприятной ситуации в группу 1 или 2 оцениваются в балла. Взаимные преобразования неблагоприятных ситуаций оцениваются баллами Бпр= или 2. Преобразование групп синоптических ситуаций при прогнозе ЭВУЗВ учитывалось в соответствии с таблицами 5 и 6.
Количественная оценка преобразования синоптических ситуаций Количественная оценка преобразования синоптических ситуаций Далее были установлены неблагоприятные особенности синоптических процессов Бос, которые способствуют значительному увеличению уровня загрязнения воздуха в городе. По всему ряду были выписаны предикторы, а также все особенности, отмечавшиеся в городах региона в дни с ЭВУЗВ. Прогноз каждой из перечисленных ниже особенностей добавляет к Sn от половины балла до двух.
В холодный период к ним относятся: приземная инверсия, очаг тепла или периферия очага тепла над городом по карте ОТ 500/1000, туман, дымка, наличие малоподвижного фронта окклюзии в окрестностях города и расположение северо-западнее города обширного глубокого циклона на расстоянии не менее 500 км. В летний период к таким особенностям относятся: приземная инверсия в сочетании со скоростью ветра 0-1м/с, туман, дымка, малоподвижный фронт окклюзии в окрестностях города.
В летний период также были выявлены особенности синоптических процессов, способствующие значительному снижению уровня загрязнения воздуха в городе при неблагоприятных ситуациях, это: расположение очага холода или периферии очага холода над городом на карте ОТ 500/1000, если прогнозируется прохождение холодного фронта и температура в слое 500/1000 менее 5оС. В этом случае баллу особенностей Бос присваивается обратный знак (таблицы 7 и 8) Количественная оценка особенностей синоптических В разделе 3.3 представлен сравнительный анализ случаев ЭВУЗВ в городах Уральского региона. Он был проведен отдельно для летнего и зимнего сезонов. При анализе наиболее опасных эпизодов загрязнения воздуха основной интерес представляют случаи, когда значительное увеличение концентраций вредных веществ в воздухе имеет место одновременно в разных городах. Известные катастрофические эпизоды характеризовались формированием ЭВУЗВ одновременно на всей территории Уральского региона. В г. Перми в зимний период отмечалось 23 случая ЭВУЗВ. В 18 случаях в одном или в нескольких городах региона также отмечался ЭВУЗВ, что составляет 78% от всех случаев. В летний период из 11 случаев ЭВУЗВ, отмечавшихся в г. Перми, в 8 случаях ЭВУЗВ отмечался и в других городах региона, что составляет 73% от всех случаев.
Для других городов Уральского региона был проведен такой же сравнительный анализ. Проведенный анализ показал, что в городах Уральского региона в 88% - 100% случаев одновременно отмечается высокий и экстремально высокий уровень загрязнения воздуха В разделе 3.4 приводятся схемы прогноза ЭВУЗВ, разработанные с использованием метода графической регрессии. В схемах были использованы два, но наиболее значимых предиктора: синоптический предиктор Sn и показатель исходного уровня загрязнения воздуха в городе параметр Р`. В результате анализа синоптических условий загрязнения воздуха в рассматриваемых городах Уральского региона было установлено, что зависимость параметра Р от синоптических условий в каждом из городов практически одинаковая. Это позволило решать задачу прогноза ЭВУЗВ в городах региона, принимая одинаковые значения предиктора Sn для каждого из городов, в том числе одинаковые неблагоприятные синоптические ситуации (Бт и Бп), неблагоприятные особенности (Бос) и неблагоприятные преобразования процессов (Бпр). Графики были построены для городов Пермь, Соликамск, Березники, Губаха и Екатеринбург, отдельно для зимнего и летнего периодов. Всего построено 10 графиков. На графиках красной линией отделена область, соответствующая ЭВУЗВ в городе. Схема включает только один график. Она очень удобна для использования в оперативной работе. По прогностическому значению Sn и рассчитанному на момент прогноза значению параметра Р по графикам находится ожидаемый уровень загрязнения воздуха в городе. На рисунках 6 и 7 представлены схемы прогноза для города Перми.
В летний период ЭВУЗВ в городах региона отмечается при значениях синоптического предиктора Sn = 7 - 9,5; в зимний период при Sn = 8 - 9,5.
В разделе 3.5 описаны результаты применения метода множественной линейной регрессии с предварительным исключением нелинейности связей. При разработке схем прогноза ЭВУЗВ использовался ряд наблюдений, в который включены были только дни с повышенным уровнем загрязнения воздуха, поэтому линейная связь между предиктантом и предикторами была нарушена. В связи с этим в данной схеме предикторы преобразованы. Разработка схемы прогноза ЭВУЗВ по методу множественной линейной регрессии с предварительным исключением нелинейности связей выполнена с использованием специально составленной программы на ПЭВМ, выполняющей следующие операции:
1. Анализ связей между предиктантом и предикторами, 2. Преобразование предикторов в соответствии с указаниями РД 52.04.306- 3. Расчет весовых коэффициентов и свободного члена уравнения регрессии.
Для включения в схему различных предикторов, сначала была проведена оценка значимости каждого из них. Для холодного периода, наиболее значимыми предикторами являются количественный синоптический предиктор Sn, значение параметра Р, предшествующее прогнозу и форма термического поля (очагу тепла на карте ОТ500/1000 придается =5, области холода - =1, промежуточные состояния – периферия очага тепла, прямолинейные изогипсы, периферия области холода). Для зимнего периода было получено прогностическое уравнение:
где Р(Sn), Р( ), Р(Р`) – преобразованные предикторы Для теплого периода, наиболее значимыми предикторами являются количественный синоптический предиктор Sn и значение параметра Р, предшествующее прогнозу. Было получено следующее прогностическое уравнение:
где Р(Sn), Р(Р`) – преобразованные предикторы Прогностические уравнения были разработаны по материалам города Перми. Для решения вопроса об использовании этих схем прогноза в других городах Уральского региона проводилась их апробация по материалам города Екатеринбурга за шестилетний период, т. к. анализ синоптического материала показал, что синоптические ситуации неблагоприятные для города Перми, в основном являются неблагоприятными и для других городов региона. Апробация схем прогноза по материалам города Екатеринбурга показала неплохие результаты. Оправдываемость схем по материалам г. Екатеринбурга для летнего и зимнего периода получилась 100%, предсказуемость случаев ЭВУЗВ для зимнего периода составила 67%, для летнего периода – 62%.
Таким образом, можно сделать вывод, что схемы прогноза ЭВУЗВ по методу множественной линейной регрессии с предварительным исключением нелинейности связей, разработанные для города Перми, могут быть использованы в оперативной работе и в других городах Уральского региона.
В разделе 3.6 приводится оценка эффективности схем прогноза ЭВЗУВ Испытание разработанных схем прогноза для городов Уральского региона проводилось на зависимом и независимом материале отдельно для зимнего и летнего периода.
Эффективность схем прогноза ЭВУЗВ оценивалась с помощью коэффициента корреляции R между фактическими и прогностическими значениями параметра Р, при помощи критерия Багрова Н* (Эффективной считается схема при Н* > 0,2), также рассчитывались оправдываемость прогнозов ЭВУЗВ и их предсказуемость. Критерий Багрова рассчитывался по формуле:
n11–число оправдавшихся прогнозов с явлением; n22–число оправдавшихся прогнозов без явления;
n21–не давалось, а было; n00–общее число прогнозов; n10–число прогнозов с явлением;
n20–число прогнозов без явления; n01–число дней с явлением; n02–число дней без явления; n12–давалось, а было без явления.
Для Перми за зимний период по зависимому материалу методом графической регрессии было составлено 24 прогноза ЭВУЗВ, оправдалось 23. Из 18 случаев ЭВУЗВ предсказаны все 18. Для летнего периода составлено 12 прогнозов ЭВУЗВ, из них оправдалось 8. Из 11 случаев ЭВУЗВ с Р0,42 предсказаны 10.
Эффективность прогностической схемы методом множественной линейной регрессии с предварительным исключением нелинейности связей также оказалась весьма высокой. Для зимнего периода коэффициент корреляции между прогностическими и рассчитанными по данным измерений значениями параметра Р составил 0,90. Для летнего периода множественный коэффициент корреляции составил 0,74. В зимний период было составлено 19 прогнозов ЭВУЗВ, из которых оправдались все 19. Всего было отмечено случаев ЭВУЗВ, предсказаны все 18. В летний период было составлено 8 прогнозов ЭВУЗВ, из которых оправдались 7. Из отмечавшихся 11 случаев ЭВУЗВ предсказаны были 10. В таблице 10 представлены результаты испытания прогностических схем для г. Перми на зависимом материале по обеим схемам.
Схемы для зимнего периода, разработанные по методу графической регрессии и по методу множественной линейной регрессии с предварительным исключением нелинейности связей, оказались более эффективными, чем для летнего периода.
Эффективность прогностических схем для летнего и зимнего периодов на зависимом регрессии В таблице 11 представлены результаты испытания прогностических схем ЭВУЗВ на зависимом материале для пяти городов Уральского региона. Оправдываемость прогнозов ЭВУЗВ по всем городам Уральского региона достаточно высокая и составляет в зимний период 89%, в летний период – 79%. В зимний период было составлено прогнозов ЭВУЗВ, из них оправдалось 85; в летний период составлено 82 прогноза ЭВУЗВ, оправдалось 65. Предсказуемость случаев ЭВУЗВ также высока: из отмечавшихся в зимний период 63 случаев ЭВУЗВ предсказаны 59; в летний период отмечалось 57 случаев ЭВУЗВ, предсказаны 49. Как видно из таблицы 11 значения критерия Багрова свидетельствуют об эффективности разработанных схем прогноза. В летний период они равны 0,73 - 0,84 и в зимний период 0,79 - 0,99.
Эффективность прогностических схем для летнего и зимнего периодов на зависимом Графической регрессии На независимом материале испытание схем прогноза ЭВУЗВ для городов Пермь, Березники, Соликамск и Губаха было проведено по данным за 2006 - 2008 годы, для города Екатеринбурга - по материалам 2009 -2010 годов.
Методом графической регрессии и методом множественной линейной регрессии с предварительным исключением нелинейности связей для Перми за летний период было составлено 78 прогнозов, за зимний период – 54 прогноза. Из таблицы 12 видно, что большую эффективность показала схема, разработанная методом графической регрессии.
Эффективность прогностических схем для летнего и зимнего периодов по В зимний период методом графической регрессии и методом множественной линейной регрессии с предварительным исключением нелинейности связей было составлено 8 прогнозов ЭВУЗВ, из которых все оправдались. В летний период методом графической регрессии было составлено 13 прогнозов ЭВУЗВ, из которых оправдались. Методом множественной линейной регрессии с предварительным исключением нелинейности связей было составлено 8 прогнозов ЭВУЗВ, из которых оправдались.
Предсказуемость случаев ЭВУЗВ для Перми по независимому материалу также достаточно хорошая. Из отмечавшихся в зимний период 8 случаев ЭВУЗВ, все были предсказаны. В летний период отмечалось 10 случаев ЭВУЗВ, методом множественной линейной регрессии были предсказаны 7 случаев, методом графической регрессии предсказаны 8 случаев. Значения критерия Багрова как для зимнего периода, так и для летнего по обеим схемам получились высокими (0,81 – 0,99). Это свидетельствует об эффективности разработанных схем.
В таблице 13 представлены результаты испытания прогностических схем для летнего и зимнего сезонов на независимом материале для всех городов Уральского региона.
Графической Всего по городам Уральского региона в зимний период было составлено прогнозов ЭВУЗВ, 34 прогноза оправдались; в летний период составлено 39 прогнозов ЭВУЗВ, из них оправдалось 36.
В зимний период отмечалось 35 случаев ЭВУЗВ, предсказаны были 32; в летний период из отмечавшихся 34 случаев ЭВУЗВ, предсказаны были 27 случаев.
Испытание схем прогноза ЭВУЗВ для городов Уральского региона на зависимом и независимом материале оказалось довольно успешным. Таким образом, можно сделать вывод о реальной возможности успешного прогнозирования ЭВУЗВ в городах региона.
Для определения успешности разработанных схем прогноза ЭВУЗВ по всему ряду наблюдений были составлены прогнозы уровней загрязнения воздуха с использованием прогностической схемы с учетом количественного синоптического предиктора (используется при ежедневном составлении прогнозов уровней загрязнения воздуха) и схемы прогноза ЭВУЗВ (рисунки 8 и 9). На графиках точки в виде ромбов - это прогнозы, составленные по схемам прогноза с учетом количественного синоптического предиктора, в виде квадратов – прогнозы, составленные с использованием схем прогноза ЭВУЗВ. При прогнозе пониженного уровня загрязнения воздуха прогнозы по обеим схемам совпадают.
На корреляционных графиках видно, что связь между прогностическими и фактическими значениями параметра Р линейная и в том и в другом случае, но корреляционная связь лучше в случае использования схем прогноза ЭВУЗВ. При использовании схемы прогноза с учетом количественного синоптического предиктора коэффициент корреляции составляет 0,71 и в летний и в зимний периоды; при использованиии схемы прогноза ЭВУЗВ - в летний период он равен 0,79, в зимний – 0,87.
Рис. 8 Корреляционная зависимость между фактическими и прогностическими Рис. 9 Корреляционная зависимость между фактическими и прогностическими Таким образом при прогнозе высоких и экстремально высоких уровней загрязнения воздуха более точное значение параметра Р дают схемы прогноза ЭВУЗВ.
В разделе 3.7 представлены рекомендации по прогнозу ЭВУЗВ в городах региона, при отсутствии разработанной прогностической схемы, а также дается оценка ожидаемого состояния загрязнения воздуха в регионе с использованием метода разложения по е. о. ф.
Рекомендации по прогнозу ЭВУЗВ в каждом из городов региона при отсутствии разработанной прогностической схемы заключается в следующем:
1. Прогнозируется синоптический предиктор Sn. Ожидаемое большое значение Sn уже является определенным предостережением. (Выше было показано, что в рассмотренных городах ЭВУЗВ формируется, когда Sn превышает 8 – 9 баллов).
2. Учитывается имеющаяся информация о состоянии воздушного бассейна в городах на момент прогноза.
3. Учитывается прогноз для близлежащего города.
На основании этого делается вывод о возможном формировании ЭВУЗВ.
С использованием метода разложения по естественным ортогональным функциям были получены общие по региону характеристики загрязнения воздуха. Использовались ряды параметра Р для пяти рассмотренных городов. Ряды для расчета r содержат более случаев. Ряды параметра Р в пяти городах региона(Р(х, t)) представляются в виде суммы членов разложения, каждый из которых характеризует определенное свойство процесса загрязнения воздуха в городах региона:
где n – число членов разложения по е. о. ф. (соответствует числу рассматриваемых городов в регионе и в нашем случае равно пяти);
- коэффициенты разложения, зависящие от времени t;
- естественные ортогональные функции, компоненты которых относятся к отдельным городам.
Важная особенность метода разложения по е. о. ф. заключается в том, что первые члены разложения дают наибольший вклад в общую изменчивость изучаемой характеристики. В нашем случае первая составляющая включает 45% изменчивости загрязнения воздуха в рассматриваемом регионе. Она характеризуется одинаковыми по знаку изменениями параметра Р во всех пяти городах. Компоненты 1 во всех городах положительны. Вклад второй составляющей в общую изменчивость параметра Р в регионе равен 21%. Для всех городов за исключением Екатеринбурга значения компонентов отрицательны. Вклады остальных трех составляющих незначительны.
Обобщенными характеристиками загрязнения воздуха в регионе являются зависящие от времени коэффициенты разложения Наибольший интерес представляют значения первого коэффициента разложения 1. Анализ показал, что при больших значениях 1 практически во всех городах отмечаются высокие значения Р, в основном достигающие ЭВУЗВ.
Рис. 10 Корреляционная зависимость между синоптическим предиктором Sn и Существует достаточно сильная корреляция между 1 и параметром Р в каждом их пяти городов. Значения r колеблются от 0,70 до 0,54. Большие значения 1 отмечаются и при больших значениях синоптического предиктора (Sn=8 и более), причем происходит одновременное повышение и понижение этих двух характеристик. Коэффициент корреляции между 1 и Sn составляет 0,74 (Рис. 10).
В результате выполненного исследования получены следующие основные результаты:
1 Создана эффективная региональная система прогноза экстремально высоких уровней загрязнения воздуха в крупном промышленном регионе (на примере Уральского региона), с целью предотвращения высоких уровней загрязнения воздуха, формирующихся в периоды неблагоприятных метеоусловий.
2. Разработаны и внедрены в практику схемы прогноза ЭВУЗВ для пяти городов Уральского региона.
3. Выявлена тенденция к одновременному формированию экстремально высоких уровней загрязнения воздуха (ЭВУЗВ) в близлежащих городах промышленного региона;
4. Показана возможность и целесообразность использования одного и того же синоптического предиктора Sn для промышленных городов в пределах региона;
5. Установлена эффективность синоптико-статистических схем прогноза ЭВУЗВ, включающих два наиболее значимых предиктора – синоптический предиктор Sn и параметр загрязнения воздуха в предшествующий день Р`.
6. Результаты испытания прогностических схем для г. Перми на зависимом материале показали, что схемы для зимнего периода, разработанные по методу графической регрессии и по методу множественной линейной регрессии с предварительным исключением нелинейности связей, оказались более эффективными, чем для летнего периода.
7. По материалам наблюдений в городах Уральского региона показано, что в дни с большими значениями параметра Р, соответствующими ЭВУЗВ, в городах реально имеют место опасные концентрации примесей в воздухе, значительно превышающие в ряде случаев ПДК.
8. По Уральскому региону впервые проведена типизация синоптических процессов, способствующих накоплению вредных примесей в атмосфере с учетом географического происхождения, траектории движения барических образований, скорости ведущего потока по картам АТ 700 гПа или АТ 500 гПа, термического поля над районом, а также длительности нахождения барического образования над рассматриваемым районом.
9. Даны рекомендации по прогнозированию ЭВУЗВ в промышленном регионе с использованием единого синоптического предиктора, в том числе в городах, для которых имеющихся материалов недостаточно для разработки прогностических положений.
10. По материалам пяти городов Уральского региона получены общие характеристики загрязнения воздуха с использованием метода разложения по естественным ортогональным функциям для оценки ожидаемого общего состояния загрязнения воздуха в регионе Результаты диссертации опубликованы в следующих работах:
1. Костарева Т.В. Схемы прогнозирования ЭВУЗВ для г. Пермь. //Информационный бюллетень за 2005 год. Состояние работ по прогнозу загрязнения воздуха в городах Российской Федерации. – М., 2006. с. 36-41.
2. Костарева Т.В., Сонькин Л.Р. Некоторые вопросы предупреждения и предотвращения высоких уровней загрязнения воздуха в периоды НМУ (на примере г.
Пермь и Пермской области). В сб.: Вопросы охраны атмосферы от загрязнения.
Информационный бюллетень НПК «Атмосфера» при ГГО, № 2 (32).
3. Костарева Т.В., Кириллова В.И., Николаев В.Д., Сонькин Л.Р. Прогноз экстремально высоких уровней загрязнения воздуха в городах промышленного региона. // Труды ГГО. – 2010. – Вып. 561. С. 92-103.
4. Костарева Т.В. Прогноз ЭВУЗВ в городах Уральского региона. В сборнике:
«Глубинное строение, геодинамика, тепловое поле земли, интерпретация геофизических полей». Материалы конференции. Екатеринбург: УрО РАН. 2011. С.205-208.