litceysel.ru
добавить свой файл
1
Мониторинг пограничного слоя атмосферного воздуха


(в условиях промышленного города)

Н.Н. Безуглова, Г.С. Зинченко, К.Ю. Суковатов

Институт водных и экологических проблем СО РАН, г. Барнаул


Monitoring of atmospheric boundary layer (in conditions of industrial town)

N.N. Bezuglova, G.S. Zinchenko, K.Yu. Sukovatov

Institute for Water and Environmental Problems SB RAS, Barnaul


Atmospheric boundary layer (ABL) is characterized by wide spatio-temporal variability. Using up-to-date data instrument for local and remote observation of different surface air characteristics lets amplify modeling and monitoring resources. By the instrumentality of ultrasonic weather complex AMC-03 located at IWEP SB RAS (Barnaul) the continuous series of meteorological and turbulent parameters for surface air by target average time were got. To estimate the turbulent regime in town surface air the meteorological and turbulent parameters for different synoptic processes were analyzed. Limit of range of certain turbulent parameters for different synoptic processes was defined.


Введение

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

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

Использование современных информационно-измерительных средств локальных и дистанционных измерений различных характеристик пограничного слоя атмосферы позволяет расширять и совершенствовать возможности моделирования и мониторинга атмосферных процессов. Применение акустического метеорологического комплекса АМК-03 позволяет использовать пульсационный метод для определения турбулентных характеристик атмосферного воздуха. В акустическом термоанемометре датчика метеопараметров используется методика измерений трех ортогональных компонент вектора скорости ветра и температуры воздуха с помощью ультразвука, основанная на зависимости его групповой скорости в атмосфере от температуры и скорости ветра. Применение ультразвукового метода измерения температуры и скорости ветра и обеспечивает малую инерционность измерений (не более 10-3 сек), очень высокую чувствительность к турбулентным изменениям указанных метеорологических параметров, отсутствие влияния солнечной радиации на результаты измерений. Вычисление трех ортогональных компонент скорости ветра осуществляется с учетом углов отклонения Gx, Gy термоанемометра от вертикального положения, что позволяет снизить требования по точности установки термоанемометра по отношению к вертикали, а также скомпенсировать ошибки, возникающие при раскачивании термоанемометра сильными порывами ветра. С помощью акустического метеорологического комплекса АМК-03 были получены мгновенные значения некоторых метеорологических параметров. Из мгновенных значений компонент скорости ветра Vx, Vy, w и температуры воздуха T вычислялись временные ряды турбулентных пульсаций температуры и компонент скорости ветра - продольной и поперечной относительно направления среднего вектора ветра <Vm>, а также вертикальной составляющей скорости ветра .


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


Основные положения теории Монина-Обухова

В настоящее время для изучения атмосферной турбулентности применяется теория Монина-Обухова. В нижней части атмосферного пограничного слоя, где можно пренебречь силой Кориолиса (высота несколько десятков метров [1]) уравнение баланса турбулентной энергии в предположении стационарности и горизонтальной однородности потока ( и ) будет иметь вид:

, (1)

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


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

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

. (2)


Градиент скорости определяется следующим образом

, (3)

Где , и (4)

Где – масштаб Монина-Обухова, -динамическая скорость. Наиболее надежно определены функции для скорости ветра и температуры, которые в общепринятой форме в настоящее время имеют вид [2,3]:


, ,

, , , (5)

где величины имеют небольшой разброс по данным разных источников =5-7 =15-16.

По результатам многочисленных исследований была установлена применимость теории подобия для описания характеристик атмосферной турбулентности над однородной подстилающей поверхностью, по крайней мере, на некотором расстоянии от поверхности, в широком диапазоне условий () [2-4].


Анализ рассчитанных величин параметров турбулентности

Описанные выше положения являются основой для вычисления параметров атмосферной турбулентности по данным ультразвукового анемометра-термометра. Измеряемыми величинами являются: температура воздуха, скорость и направление горизонтального ветра, скорость вертикального ветра, относительная влажность воздуха, атмосферное давления. По измеренным величинам мгновенных значений метеопараметров и их пульсаций рассчитываются 15 параметров турбулентности. Для оценивания режима турбулентности в приземном слое атмосферы промышленного города были проанализированы имеющиеся данные метеорологических параметров и турбулентных характеристик для различных типов синоптических процессов: антициклонального и процессов, связанных с прохождением атмосферных фронтов. Поскольку речь идет о данных, получаемых от метеорологического комплекса, установленного в одной точке (здание ИВЭП СО РАН, высота 25 м) в 2009 г., то можно говорить только об анализе временных изменений метеорологических параметров и турбулентных характеристик атмосферного воздуха города, полученных за период 2009-2010 гг.


В нашей работе выполнен анализ двух параметров турбулентности: вертикального потока тепла () и полной энергии турбулентности ():

,
(6)

где
,,- дисперсии турбулентных пульсаций компонент скорости ветра , , ;

,
(7)

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

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




Рис.1. Суточный ход полной энергии турбулентности

Результаты и обсуждение

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


Как видно из рис.1 в условиях антициклона наблюдается выраженный суточный ход полной энергии турбулентности, в дневные часы при увеличении флуктуаций скорости ветра увеличивается энергия турбулентности, величина ее меняется от 0,5 до 4,5 м22. При прохождении атмосферного фронта (19-20 часов) величина полной энергии турбулентности увеличивается на два порядка и достигает м22.

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

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


Список литературы


  1. А.С. Монин, А.М. Яглом. Статистическая механика. Москва: Наука, 1965,- 639 с.

  2. J.A. Businger, J.C. Wyngard, Y. Isumi and E.F. Breadley. Flux-profile relationships ithe atmospheric surface layer. J.Atmos.Sci., vol.28, 1971.- P. 181-189.

  3. A.J. Dyer. A review of flux-profile relatiinships. Boundary-Layer Meteorol., vol. 7, 1974.- P.363-372.

  4. А.В. Гармашов. Исследования приводного слоя атмосферы диссипативным

методом. /Системы контроля окружающей среды: Сб. науч. Тр. НАН Украины.

МГИ:-Севастополь. 2005.- 376с.