CC BY
63
Серия «Науки о Земле» И З В Е С Т И Я
2013. Т. 6, № 2. С. 47-56 Иркутского
Онлайн-доступ к журналу: государственного
http://isu.ru/izvestia университета
УДК 551.510.42(571.53)
Оценка загрязнения воздушной среды г. Иркутска автотранспортом
A. В. Аргучинцева ([email protected])
B. К. Аргучинцев ([email protected])
C. А. Новикова ([email protected])
Аннотация. Рассматривается проблема загрязнения атмосферного воздуха г. Иркутска автотранспортными средствами. На основе математической модели строятся области распространения загрязняющих веществ, поступающих в атмосферный воздух от автомобилей на участке дороги. Оценивается уровень загрязнения атмосферного воздуха на территории, прилегающей к предприятию «Иркутский аэропорт».
Ключевые слова: город, загрязнение, атмосфера, автотранспорт, моделирование.
Введение
В городах автотранспорт является одним из наиболее значимых загрязнителей окружающей среды. Не является исключением и г. Иркутск, в котором количество автотранспорта постоянно возрастает, доля старых машин с эксплуатацией свыше десяти лет велика. Плохое качество дорог и дорожного покрытия, постоянно разрушающегося из-за больших перепадов температур; отсутствие развязок; узкие улицы, загруженные с утра до вечера с обеих сторон припаркованным автотранспортом; отсутствие подземных парковок - все это приводит к недостаточной пропускной способности, возникновению автомобильных пробок, часто повторяющимся циклам торможения и трогания с места (во время работы на холостом ходу выброс ингредиентов увеличивается примерно в шесть раз по сравнению с выбросом движущегося автомобилея). В выхлопных газах содержится свыше 200 различных загрязняющих ингредиентов, наиболее токсичными из которых являются оксид углерода, сажа, формальдегид. Попадая в атмосферу из низко расположенных труб глушителей, они рассеиваются на уровне дыхания детей и взрослого человека.
На рассеяние примесей в условиях города существенно влияют планировка улиц, их ширина и направление, высота зданий, наличие зеленых массивов - все, что образует как бы разные формы наземных препятствий воздушному потоку и приводит к возникновению локальных метеорологических условий в городе.
Наблюдения показывают, что даже при постоянных объемах и составах транспортных выбросов в результате влияния метеорологических условий уровни загрязнения воздуха могут различаться в несколько раз. Учет этого влияния важен при подготовке документов о качестве атмосферного воздуха, разработке воздухоохранных мероприятий, планировании размещения городов и промышленных объектов, прогнозировании уровня загрязнения.
В связи с этим при оценке эффективности выполнения мероприятий по охране атмосферного воздуха недостаточно иметь только сведения о сокращении выбросов. Требуется надежная информация за длительный период о содержании примесей в атмосфере и климатических условиях их распространения в атмосфере [3; 6].
Объекты исследования и подготовка данных
На первом этапе были проведены расчеты по методике [5], которая используется для оценки загрязнения воздушной среды движущимся автотранспортным потоком на участке автодороги. Данная методика устанавливает порядок расчета и определения выбросов автотранспорта для их использования при проведении сводных расчетов загрязнения атмосферы городов; может быть применена ко всем категориям автотранспортных средств при эксплуатации в городских условиях. В качестве исходных данных для расчета выбросов автотранспорта в атмосферу на действующих автодорогах используются результаты натурных обследований структуры и интенсивности автотранспортных потоков с подразделением по основным группам автотранспортных средств.
Осенью 2013 г. нами проведены натурные наблюдения по изучению автотранспортных потоков (интенсивности движения) на участке автодороги ул. Ширямова (в районе аэропорта г. Иркутска). Ситуационный план исследуемой территории представлен на рис. 1. Наблюдения проводились за автомобилями трех категорий: грузовые, легковые, автобусы (маршрутные такси). Подсчет проходящих по данному участку автодороги транспортных средств проводится в течение 20 мин в часыпик: 8.30-8.50, 13.2013.40, 18.20-18.40. Подсчет проходящих по обследуемому участку транспортных средств проводился с использованием видеозаписывающего устройства (видеокамеры) (табл. 1).
Рис. 1. Ситуационная карта-схема исследуемой территории
В табл. 1 представлены данные видеонаблюдения в будние дни в часыпик на участке дороги ул. Ширямова.
Таблица 1
Данные наблюдений в будние дни на участке дороги ул. Ширямова
Тип автотранспортных средств Со стороны аэропорта В сторону аэропорта Сумма
Время: 8.30 - 8.50
Легковые автомашины 388 448 836
Автобусы, маршрутные такси 12 18 30
Грузовые автомашины 20 25 45
Время: 13.20 - 13.40
Легковые автомашины 380 454 834
Автобусы, маршрутные такси 13 16 29
Грузовые автомашины 21 27 48
Время: 18.20 - 18.40
Легковые автомашины 398 444 842
Автобусы, маршрутные такси 19 18 37
Грузовые автомашины 18 41 59
Как видно из табл. 1, в наблюдаемые сроки (утро, день, вечер) количество всех видов автотранспорта почти одинаковое (разброс обусловлен некоторой случайностью появления машин и возможными погрешностями счета).
По методике [5] были проведены расчеты выбросов СО в атмосферный воздух от автотранспортных средств трех категорий (легковые, грузовые, автобусы и маршрутные такси).
Выброс /-го загрязняющего вещества (г/с) движущимся автотранспортным потоком на автодороге (или ее участке) с фиксированной протяженностью Ь (км) рассчитывался по формуле:
т к
МЬ =----------XМЬг • • Гу г/с,
Ь/ 1200 1 к ’г к Ук *
где МЬк,/ (г/км) - удельный пробеговый выброс >го вредного вещества автомобилями к-й группы, определяемый по табл. 2; к - количество групп автомобилей; Ок - фактическая наибольшая интенсивность движения, т. е. количество автомобилей каждой из к групп, проходящих через фиксированное сечение выбранного участка автодороги в единицу времени (20 мин) в обоих направлениях по всем полосам движения; гукл _ поправочный коэффициент, учитывающий среднюю скорость движения транспортного потока (Ук/ , км/ч на выбранной автодороге (или ее участке), определяемый по табл. 3); Ь (км) - протяженность автодороги (или ее участка), из которого исключена протяженность очереди автомобилей перед запрещающим сигналом светофора, включающая длину соответствующей зоны перекрестка (для перекрестков, на которых проводились дополнительные обследования).
Таблица 2
Удельные пробеговые выбросы загрязняющих веществ МЬк/ (г/км) для различных групп автомобилей
Наименование группы автомобилей № груп- пы В ы б р о с , г/км
О N02 СН Сажа 802 Формаль- дегид Бенз(а) - пирен
Легковые I 3,5 0,9 0,8 0,7-10-2 1,510-2 3,210-3 0,3 10-6
Автофургоны и микроавтобусы, до 3,5 т II 8,4 2,1 2,4 3,810-2 2,8-10-2 8,4-10-3 0,810-6
Грузовые, от 3,5 до 12 т III 6,8 6,9 5,2 0,4 5,1-10-2 2,2-10-2 2,110-6
Грузовые, свыше 12 т IV 7,3 8,5 6,5 0,5 7,3 10-2 2,5 10-2 2,6-10-6
Автобусы, свыше 3,5 т V 5,2 6,1 4,5 0,3 4,2-10-2 1,8. 10-2 1,810-6
Таблица 3
Значения коэффициентов г^, учитывающих изменения количества выбрасываемых вредных веществ в зависимости от средней скорости движения
К Скорость движения V км/ч
V 5 10 15 20 25 30 35 40 45 50 60 70 80 100 110
Гук.1 1,4 1,35 1,30 1,20 1,10 1,00 0,90 0,75 0,65 0,50 0,30 0,40 0,50 0,65 0,75
Гук,1 да2 1,0 1,0 1,0 1,0 1,0 1,0 1,0 1,0 1,0 1,0 1,0 1,0 1,0 1,0 1,2
Удельные выбросы загрязняющих веществ представляют собой усредненные для рассматриваемых групп автомобилей значения удельных выбросов загрязняющих веществ при их движении на участках городских автомагистралей, а также в условиях их пребывания на пересечениях городских автомагистралей. При назначении удельных выбросов загрязняющих веществ для автомагистралей принималось во внимание то, что в городских условиях автомобиль на конкретном участке автодороги совершает непрерывные разгоны и торможения, перемещаясь с некоторой средней скоростью, определяемой дорожными условиями.
В табл. 4 приведены данные по выбросам СО движущимися автотранспортными средствами на участке дороги ул. Ширямова, рассчитанные по формулам, представленным выше [5]. Из таблицы видно, что максимальное количество выбросов приходится на легковые автомобили, что прямо пропорционально количеству автотранспортных средств, зафиксированных на данном участке дороге при натурных наблюдениях.
Таблица 4
Выбросы СО движущимися автотранспортными средствами на участке дороги ул. Ширямова
Тип автотранспортного средства Выброс СО, мг/с
Легковые 490
Автобусы 40
Грузовые 60
Моделирование загрязнения
Для описания ветровых потоков в городе предлагается численная модель мезомасштабных процессов, базирующаяся на основе системы гидротермодинамических уравнений, выражающих основные законы сохранения энергии, импульса и массы [1; 2]. Это уравнение движения:
dV 1 , „ г г г
— =------gradp - 2а х у + g + Dv , (1)
dt р
уравнение неразрывности:
—+р divv = 0, (2)
dt
уравнение притока тепла:
&Т_____1_&р _
& Срр &
уравнение состояния, записанное в общем виде:
р _ ЯрТ,
уравнение баланса кинетической энергии турбулентности:
ЛЬ_
Л
ди У ( ду
+ 1 -
+ с,
§ др р ді
уравнение скорости диссипации турбулентной энергии:
Лв , в
— = с1в к — Л Ь
соотношение Колмогорова: соотношение Смагоринского:
ді, к„ = сЬ2 /в
+ с,
§ д£.
р ді
(3)
(4)
— в + DЬ, (5)
2
в
— С2в — + Dв, (6)
Ь
- + -
+
дх ду) 1 дх ду
Л д д д д
где — =----+ и---+ V-----+ н—;
Л дt дх ду ді
_ д , ду д , ду д , ду Dy = — кх — +—ку — +—су кг —. дх дх ду ду ді ді
Здесь ¥ - любая из функций рассматриваемой задачи; t - время; х, у, і -оси декартовой прямоугольной системы координат; и, V - горизонтальные и н - вертикальная компоненты вектора V скорости движения воздуха; р - плотность атмосферы; р - давление; Т - температура; ш - вектор угловой скорости вращения Земли, направленный параллельно оси Земли к Северному полюсу; § - сила тяжести; ср - удельная теплоемкость при
постоянном давлении; Я - газовая постоянная сухого воздуха; Ь - кинетическая энергия турбулентности; в - скорость диссипации кинетической энергии турбулентности; кх, ку, кг - коэффициенты турбулентного обмена
по горизонтали и вертикали; с, су - эмпирические константы; си = су = сн = 1 ; Ь = к0 • А^ / 2 ; Ь - масштаб, пропорциональный шагу горизонтальной сетки А^ ; к0 - безразмерный параметр, аналогичный постоянной Кармана.
2
2
Система уравнений рассматривается в параллелепипеде ^ {х,у,ъ: 0 < х < X, 0 < у < У, 0 < ъ < И}, где х = 0, х = X, у = 0, у = У, г = 0, г = Z - границы области счета. Краевые условия задаются по данным наблюдений.
Уравнения модели (1)-(6) интегрируются в декартовой прямоугольной системе координат с применением метода фиктивных областей, который позволяет проводить расчеты в условиях городской застройки. Численный алгоритм решения строится на основе метода расщепления по физическим процессам и геометрическим переменным [4].
Найденные на основе системы (1)-(6) скорости движения и турбулентные характеристики используются для решения дифференциального уравнения переноса и турбулентной диффузии примесей:
^ + ёгу^у) - — _ Бз -аз + р, (7)
д ё дъ
где 5 - концентрация любой субстанции; а - коэффициент неконсерватив-ности примеси; Q(х, у, г, 0 - функция, описывающая распределение и мощности источников рассматриваемой субстанции; _ скорость гравита-
ционного осаждения аэрозолей, определяемая по формуле Стокса [14]:
^ _ (2Рпёг2)/(9ру), рп и г - соответственно плотность и радиус частицы взвеси; g - ускорение свободного падения; р и V - плотность и вязкость воздуха.
Из-за отсутствия детальной информации из наблюдений в качестве начальных условий принимается 5 = 0. На границах области счета по горизонтальным осям и при г = 100 м задаются условия Неймана, на уровне подстилающей поверхности - учитывается отражение и поглощение примеси в зависимости от свойств подстилающей поверхности.
Уравнение (7) решается численно методом расщепления [4].
Расчеты проводились для крупномасштабных потоков различных направлений при следующих значениях параметров: шаг по времени выбирался с учетом критерия Куранта; вертикальный градиент фоновой температуры - (-0,65 °С/100 м); верхняя граница области счета задавалась на высоте 100 м; шаги по горизонтали и вертикали 2 м; горизонтальное разрешение 100x50 точек.
В качестве иллюстрации приведем фрагменты расчетов загрязнения атмосферы над выбранным участком автомагистрали при общем глобальном переносе, обеспечиваемом ветрами направлений: восточного (рис. 2), западного (рис. 3) и наиболее повторяющегося - северо-западного (рис. 4). Из рисунков видно, как ориентация улиц и зданий изменяют направления ветров.
Рис. 2. Векторное поле скорости при восточном переносе и концентрации СО. Область сгущения краски - концентрации более 3 мг/м3
Рис. 3. Векторное поле скорости при западном переносе и концентрации СО. Область сгущения краски - концентрации более 3 мг/м3
Рис. 4. Векторное поле скорости при северо-западном переносе и концентрации СО. Область сгущения краски - концентрации более 3 мг/м3
Выводы
Используемая в работе численная модель [1, 2] является более универсальной по сравнению со стандартными методиками, поскольку она основана на наиболее полных уравнениях гидротермодинамики атмосферы с учетом динамических уравнений турбулентности. Кроме того, модель может быть применена и для исследования загрязнения атмосферного воздуха в особо опасных ситуациях, близких к штилевым.
Возможности модели и получаемые результаты могут быть использованы для создания экологических паспортов и для принятия управленческих решений по уменьшению негативного влияния транспорта на экологическую ситуацию промышленных городов. Среди таких мероприятий -это прежде всего приведение в надлежащий порядок городских магистралей (расширение улиц; асфальтирование; хорошее качество покрытия; устранение многочисленных «лежачих полицейских», выполненных в грубой форме с нарушениями техники безопасности; многоярусное озеленение), контроль за выбросами отработанных газов автомобилей в соответствии с действующими государственными стандартами, выпуск на линию автотранспортных средств в соответствии с действующими стандартами.
Список литературы
1. Аргучинцев В. К. Моделирование мезомасштабных гидротермодинамических процессов и переноса антропогенных примесей в атмосфере и гидросфере региона оз. Байкал / В. К. Аргучинцев, А. В. Аргучинцева. - Иркутск : РИО ИГУ, 2007. - 255 с.
2. Моделирование и управление процессами регионального развития / А. В. Аргучинцева [и др.]. - М. : Наука, Физматлит, 2001. - 432 с.
3. Базаров И. П. Термодинамика / И. П. Базаров. - М. : Высш. шк., 1991. -
375 с.
4. Марчук Г. И. Приоритеты глобальной экономики / Г. И. Марчук, К. Я. Кондратьев. - М. : Наука, 1992. - 263 с.
5. Методика определения выбросов автотранспорта для проведения сводных расчетов загрязнения атмосферы городов. - СПб., 2010. - 15 с.
6. Яблокова Е. П. Экологические проблемы автомобильного транспорта : аналит. обзор / Е. П. Яблокова ; ГПНТБ СО РАН. - Новосибирск : Север, 2005. -113 с.
Estimating of Air Pollution Due to Emissions by Transport in Irkutsk
A. V. Arguchintseva, V. K. Arguchintsev, S. A. Novikova
Annotation. The article is devoted to the problem of air pollution due to emissions of pollutants from transport of Irkutsk. Of the basis of mathematical model we are building area of pollution is distribution from transport of Irkutsk. Also we are assessing the level of air pollution in area, which located near Irkutsk airport.
Key words: city, pollution, atmosphere, transport, modeling.
Аргучинцева Алла Вячеславовна доктор технических наук, профессор, Иркутский государственный университет 664003, Иркутск, ул. К. Маркса, 1 декан географического факультета тел.: (395-2) 42-56-84
Arguchintseva Alla Vyacheslavovna Doctor of Technical Sciences, professor Irkutsk State University 1, K. Marx st., Irkutsk, 664003 Dean, Faculty of Geography tel.: (395-2) 42-56-84
Аргучинцев Валерий Куприянович
доктор технических наук, профессор,
Иркутский государственный университет
664003, Иркутск, ул. К. Маркса, 1
зав. кафедрой метеорологии
и охраны атмосферы
тел.: (395-2) 52-10-94, 52-10-88
Новикова Светлана Александровна аспирант
Иркутский государственный университет 664003, Иркутск, ул. К. Маркса, 1 тел: 8(9086) 61-58-58
Arguchintsev Valery Kupriyanovich Doctor of Technical Sciences, professor Irkutsk State University 1, K. Marx st., Irkutsk, 664003 Head of the Department of Meteorology and protection of the atmosphere tel.: (395-2) 52-10-94, 52-10-88
Novikova Svetlana Alexsandrovna Post-Graduate Student Irkutsk State University 1, K. Marx st., Irkutsk, 664003 tel.: 8(9086) 61-58-58
