Добавил:
Опубликованный материал нарушает ваши авторские права? Сообщите нам.
Вуз: Предмет: Файл:

книги из ГПНТБ / Бызова, Н. Л. Рассеяние примеси в пограничном слое атмосферы

.pdf
Скачиваний:
27
Добавлен:
19.10.2023
Размер:
7.53 Mб
Скачать

5.3.4. Среднесуточные

значения

ау

Учет влияния

крупномасштабных

колебаний

направления вет­

ра при времени

действия источника

порядка суток рассматривал­

ся Артемовой (1967). Она использовала предположение о том, что для приближенной оценки среднесуточные значения поперечной

дисперсии

примеси

°\уг(х)

можно

рассматривать

как результат

наложения

средних 20-минутных зон загрязнения, оси которых на­

правлены

вдоль среднего за 20 мин направления ветра. Рассеяние

внутри 20-минутных

промежутков

определяется

мелкомасштаб­

ной турбулентностью, а изменения

положения оси—мезомасштаб-

ной изменчивостью направления ветра. На материалах измерений, проведенных на высотной мачте ИЭМа в разное время года, Ар­

темовой были выделены пять типов изменения направления

ветра

в течение суток, характеристика которых приведена в табл.

5.19.

Примеры экспериментально полученных распределений повторяе­ мости среднего за 20-минутные промежутки направления ветра и аппроксимирующие их нормальные законы показаны на рис. 5.8.

 

 

 

Таблица 5.19

Типы (I V) изменения

направления ветра

и сопутствующие им барические

образования

 

(« — повторяемость)

 

Характеристика

Характеристика

баричес­

Л

направления вет­

ев .

ра в течение

ких образований

суток

 

 

V S

Рис. 5.8. Примеры экспе­ риментально полученных распределений повторяе­ мости направления ветра и аппроксимирующие их

нормальные за.ионы:

I, II, III — типы

изменений

направления

ветра

I Очень ус­

Устойчивые

бари

0,2131 12,2

26

тойчивое

ческие центры

 

22

II Средне-

Слегка подвижные1 0,335 19,2

устойчивое

барические центры

0,537 30,8

12

• III Малоус­

Формирование,

тойчивое

разрушение,

смеще-,

 

 

IV Резкий

иие барических

 

 

поворот (один центров

 

 

31

или более)

Прохождение

 

в течение

фронтов

 

 

 

суток

 

 

 

 

V Крайне

Штиль

 

 

 

неустойчивое

 

 

 

 

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

 

°2 сут(*) = 4>(*) +

°2 ос«(-*)-

 

(5.25)

В свою очередь,

а\й{х) определяется

выражением

(2.2) или

(5.24),

а колебания оси

 

 

 

 

 

°1п(*) = Вух,

 

 

(5.26)

где Ву =:"1Л<в2 >

— среднеквадратичное значение

закона

распре­

деления среднего 20-минутного направления ветра в течение суток. Отсюда

 

 

°cyr{x)

=

V B l

+ b2yX,

 

 

(5.27)

причем Ву > Ьу

почти

при

всех

комбинациях

условий,

поэтому

практически

а'1

определяется выражением

(5.26).

 

Величина

B v ,

составляющий

в условиях

средней полосы Рос­

сии от 0,20 до 0,54 в зависимости от характера

погоды,

очевидно,

зависит от климатических особенностей местности. Приведем для сравнения аналогичные, хотя не идентичные с этими данные, по­ лученные в США 1[МАЭ, 1968]. Результаты опытов с последова­ тельными выпусками тетронов в течение большого промежутка времени — от 8 до 48 ч показали, что среднеквадратичное попе­ речное отклонение траекторий тетронов от среднего их положения за период выпуска описывается выражением вида (5.26), причем значения Ву в окрестности Лос-Анжелеса, где резко выражен су­ точный ход ветра, составляли от 0,8 до 1,1, а для местности с ме­ нее выраженным суточным ходом — порядка 0,5.

В табл.

5.20 приведены рассчитанные

Артемовой отношения

° 2 о / ° с у т на

разных расстояниях для трех

типов устойчивости на­

правления ветра, а также значения показателя степени сп, полу­

ченные

при

аппроксимации

зависимости

o-c y T (jr)

от времени

из­

мерения

Т выражением

типа

(5.19). При отсутствии

резких изме­

нений направления ветра на

расстояниях

до 5 км

cti составляет

0,16—0,25. Сравнивая этот результате выражением

(5.19), можно

отметить, что как приближенное эмпирическое правило при

кли­

матических условиях, близких к условиям средней

 

полосы, . оно

пригодно в диапазоне Т от 1 мин до 24 ч.

 

 

Таблица

5.20

 

 

 

 

 

 

 

 

 

 

Значения

Ого/Осут на разных

расстояниях

от источника

 

 

 

 

и значения показателя степени

а\

 

 

 

 

 

Тип изменения направления

 

 

 

X км

I

 

II

 

 

 

in

 

 

 

 

 

 

 

 

 

 

 

а-.,(|с у т

 

а.„/сг

 

сут

 

 

 

 

 

 

сут

 

 

 

 

 

 

о д

0,58

0,12

0,42

0,20

0,27

 

0,28

 

 

0,5

0,52

0,15

0,36

0,24

0,23

 

0,34

 

 

1,0

0,50

0,16

0,34

0,25

0,22

 

0,36

 

 

5,0

0,44

0,19

0,30

0,27

0,19

 

0,39

 

!0,0

0,41

0,20

0,28

0,28

0,17

 

0,41

 

5.4. Максимальная плотность осадка. Длина следа

5.4.1. Зависимость осевой концентрации от времени действия источника

Если

два источника выпускают одинаковую примесь с одной

и той же

высоты

в течение

разного времени Т и Т

соответствен­

но, то, согласно

(2.1),

отношение осевых

плотностей осадка р и

р' имеет

вид

 

 

 

 

 

 

,

 

р(х)

_

°T'{x)qi(x)

<

1

(5.28)

 

 

р'(х)

 

о у (х) <7i (х)

 

 

 

 

 

 

при Т'<Т. Функции qi(x) и ц\ {х) определяются гравитационным оседанием и вертикальной диффузией, которая в рамках полуэм­ пирической теории не зависит от времени действия источника, и выражение (5.28) при q'(x)=q(x) может служить оценкой р{х)1р'{х) в этом предположении.

Другая оценка может быть получена из (1.50). В случае изо­ тропной турбулентности на оси струи имеет место соотношение

q' (х)

 

(5.29)

о\ {х)

вт (х)

Вне оси, на фиксированном уровне

qt(x. z)

- - 5

°\Лх)

ехр

 

 

(5.30)

•Ц-^

 

2а-г (х)

2ау

Ц\ (A", Z)

 

ат(х)

 

( J C )

где z0 и z'0 — координаты оси струи, осредненной за время Т и Т' соответственно. Экспоненциальный множитель в (5.30) из-за меандрирования струи в вертикальном направлении может прини­ мать различные значения очень большого диапазона. Поскольку, однако, временной масштаб диффузии в вертикальном направле­ нии существенно меньше, чем в горизонтальном, на тех расстоя­ ниях от источника, где наблюдается наземный след, влияние меандрирования при 2 = 0 может быть большим только при очень сильной неустойчивости и малом Т (менее I мин).

Обратимся теперь снова к

результатам

опытов, проведенных

при прочих равных условиях,

но с разной

продолжительностью

действия источника, с точки зрения параметра вертикальной диф­ фузии. Эти результаты показывают, что положение максимума плотности осадка, действительно, очень мало зависит от Т и опре­ деляется в основном метеорологическими условиями. Тем не менее

оказалось, что при меньшем времени

действия

источника макси­

мум располагается, как правило, немного

дальше (в 17

случаях

из 23), хотя в

среднем этот

эффект

лежит

в пределах

одного и

того же класса

устойчивости.

Таким

образом,

при Т'<Т

эффек­

тивный коэффициент диффузии оказывается как бы несколько за-

1 п<\

ниженным. Количественная характеристика этого занижения

Го =

(5.31)

В' '

где В и В' — значения В, соответствующие Т и Т', приведена в табл. 5.21, построенной по результатам опытов табл. 5.9.

 

 

 

 

 

Таблица 5.21

Средние отношения

осевых концентраций

(р),

поперечных

размеров

следа

( а г ), их квадратов и характеристик

вертикального рассеяния Г\ и г2 во время парных опытов

при разном

времени действия

источника Т'<Т

для серий

 

 

табл.

5.9

 

 

Серия

Г' ' Г

pip'

Т' т

г.

 

а

0,42

0,56

0,64

0,88

1,17

б

0,32

0,50

0,57

0,88

0,79

в

0,50

0,58

0,71

0,82

0,86

г

0,59

0,56

0,74

0,73

0,83

Среднее

0,41

0,53

0,64

0,84

0,85

Выражение (5.28) может служить ориентировочной оценкой от­ ношения р(х)1р'(х) с избытком, а (5.29) с недостатком. Как показано в работе Осипова (1969), при Т<Т соотношение

4 '

^ Р(х) ^

^_

< 1

(5.32)

ат

р' (х)

от

 

 

в среднем достаточно хорошо подтверждается, что видно также из данных табл. 5.21. Приведенные в ней средние значения ° г ' О т ­ служат оценкой влияния времени действия источника на горизон­ тальное поперечное рассеяние частиц. Для оценки зависимости вертикального рассеяния от Т приведены две независимые вели­ чины: г2 по (5.31) и полученная из (5.28)

р (X) а г (X)

<7, (X)

(5.33)

г, = р'(х)сг{х)

q l ( x )

 

которые оказались в среднем равными, чего и следовало ожидать, поскольку они характеризуют один и тот же эффект. Отметим, что если в среднем при изменении Т в 10—15 раз (как было во время опытов) количественная характеристика горизонтальной диффузии меняется в пределах 35%", то количественная характеристика вер­ тикальной диффузии при этом меняется значительно меньше — около 15%.

5.4.2. Максимальная

плотность осадка

Поскольку максимальная плотность осадка связана с попереч­ ной шириной следа, которая зависит от времени действия источ-

нйка, Целесообразно рассматривать

значения

_ р0Нзу

(л-0 )

Рп =

(5.34)

 

QB

в функции от v = - — . Такая нормировка позволяет обобщить все

BU

результаты, независимо от класса устойчивости и других пара­ метров.

Измеренные значения максимальной плотности осадка обнару­ живают большой разброс, при близких значениях v величина рп может меняться на порядок, а отдельные значения (как правило, полученные при особых метеорологических условиях, в частности, во время очень слабого ветра) даже выходят за эти пределы. На рис. 5.9 приведены сгруппированные и осредненные значения рп,

каждая точка представляет результаты 8—12 опытов. Значения В для нормировки выбирались в зависимости от высоты источника и класса устойчивости.

На рис. 5.9 даны также аналогичные результаты из работы Стюарта (1968), где приведены максимальные значения проинте­ грированной по оси у плотности осадка и теоретические кривые, полученные по схеме с линейно растущим k(z) (для сравнения с данными ИЭМ — кривая а, с данными Стюарта — кривая б). Начиная с v = l эти кривые практически совпадают. Отметим, что результаты ИЭМа имеют тенденцию при больших v ложиться не­ сколько ниже теоретической кривой, а результаты Стюарта — не­ сколько выше. Самая нижняя точка рисунка получена по данным

работы Осипова и др. (1970), для опытов была

использована

прак­

тически

однородная

примесь,

состоящая

из частиц

размером

16—20 мкм. Самая

левая,

напротив,

может

содержать

ошибки,

связанные с неоднородностью

составов

1 и

(см. табл.

4.6).

Выражения (2.2),

(2.67)

и

(5.24)

позволяют

построить

зави­

симость

величин PoHoy{xQ)/.Q

и p0H2IQ

от w/U,

если принять для

каждого класса устойчивости среднее значение В и во втором слу­ чае by. Соответствующие кривые представлены на рис. 4.14—4.17. Оказалось, что хотя значения v при полном диапазоне изменения устойчивости различаются более чем на порядок, полученные за­ висимости poH2/Q и poHoy/Q от v очень близки, что позволяет для практических расчетов рекомендовать простые формулы вида

5.35)

и

(5.36)

Согласно формулам (2,69) — (2.72) и (2.80) — (2.83) в пре­ деле при v со показатели у\ и у2 принимают значения 2,5 и 1,5,

i

а при v-*0 оба они стремятся к единице, если граничное условие имеет вид (2.61), и к нулю при условии полного поглощения на границе и &(2)=const. Экспериментальные результаты, представ­ ленные на рис. 5.9, не позволяют выявить различие, возникающее

Рис. 5.9. Зависимость

рп 0Наy(x0)/Q

от

v=w/UB:

 

/ — по результатам опытов

ИЭМа, 2— по результатам Стюарта

(1968);

о —расчетная кривая к точкам /, б —то же для точек 2

 

при малых v, но противоположный

предел

при больших v здесь

достигается.

 

Аи А2,

yi

 

 

В табл. 5.22 приведены значения

и уг, полученные

по экспериментальным данным как регрессия

po#2 /Q и

poKay/Q

на w/U. Приведены также значения

этих

параметров,

соответст-

гующие эмпирической зависимости

(Петрова и Мирошкина, 1967)

для #>300 м. Значения А2 и у2 получены с учетом эмпирической зависимости из этой же работы

Рис. 5.10. Зависимость Х/х0 от v

по

результатам опытов ИЭМа; расчет- .

пая кривая

по

(5.9) н рис. 5.3

Экспериментальные значения Л2 и у2 для классов устойчивости 4—6 и данных Петровой и Мирошкиной очень плохо согласуют­ ся с расчетными при v-со, а для классов 1—3 они меньше, так

 

 

Значения А, 71, А2 и \2

 

 

 

Таблица

5.22

 

 

в формулах (5.35) и (5.36)

 

 

 

Автор, класс

устойчивости

 

1,

л.

 

 

 

 

-

 

 

 

 

 

 

ИЭМ,

1 — 3

 

0,28

1,18

0,20

1,12

 

 

ИЭМ,

4—6

 

1,7

1,63

1,4

1,51

 

 

Петрова и Мирошкмна, [1967]

200

2,8

1,6

1,50

 

 

Предельная при

 

60

2,5

2,3

1,50

 

 

-1 -*• сю , класс 6

 

 

То же, класс 5

 

1,6

1,50

 

 

То же, класс 4

 

 

1,1

1,50

 

как

в этих случаях и значения v в среднем

меньше. Значения

А\

и у\

согласуются

с расчетными

значительно

хуже из-за

искажаю­

щего влияния поперечного рассеяния, которое, кроме класса устой­ чивости, зависит от времени действия источника.

5.4.3. Длина следа. Общее количество осевшей примеси

Характеристикой спадания плотности осадка за зоной его мак­ симального значения служит величина Аз (2.68). Эксперименталь­ но полученные значения Ao.t =Я, т. е. расстояния до точки, где плотность осадка составляла 0,1 ро, обнаруживают большой раз­ брос. Согласно (5.9) эта величина линейно связана с х0 и зависит от параметра v.

На рис. 5.10 представлены после группировки и осреднения экспериментальные значения Х/х0 в зависимости от v, а также те­ оретическая зависимость, полученная с помощью уравнения (2.68). Из рисунка видно, что в среднем характер зависимости Х/х0 от v совпадает с теоретическим, хотя отклонения довольно велики. Они могут быть, в частности, связаны с влиянием неоднородности со­ става.

Рис. 5.11. Связь теоретически

(т)т ) и эксперименталь­

но (г|э) рассчитанных

количеств уловленной

на пло­

щадь полигона

примеси

ло отношению к выпущенной:

/ — регрессия

т)э нат)т ,

2 — ее

доверительные

пределы; 3 —

 

•г)э=0,9Б т)т

 

 

Общее количество осевшей на площадь полигона примеси рас­ считывалось по картам-схемам с помощью планиметрирования. Отношение этой величины к общему количеству выпущенного аэро­ золя т)э сравнивалось с г\г— этой же величиной, рассчитанной в соответствии с (5.8), в предположении, что в поперечном к оси следа направлении вне площади полигона неучтенной примеси практически не содержится. На рис. 5.11 представлено соответ­ ствие величин т]э и т]г, коэффициент корреляции между ними составляет 0,65, в среднем имеет местот)э = 0,96 -цт

5.5.Дополнительные результаты

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

дущих разделах, «о позволяют рассмотреть некоторые частные во­ просы, каждый из которых освещается тон или иной выборкой из массива опытов.

5.5.1. Особые метеорологические условия

В качестве особых метеорологических условий при проведении опытов были выделены слабые ветры (25 опытов), приподнятые инверсии (семь опытов) и слабые осадки (четыре опыта). Малое количество результатов с приподнятой инверсией и слабыми осад­ ками не позволило сделать существенных выводов по этим груп­ пам особых условий. Можно только отметить, что во время при­ поднятых инверсий параметр вертикальной диффузии близок к зна­ чениям, типичным для безразличной стратификации (табл. 5.5 и 5.8), а параметр горизонтальной диффузии оказывается промежу­ точным между значениями, типичными для безразличной страти­ фикации и слабой устойчивости (классы 4 и 5). При слабых же осадках параметр горизонтальной диффузии соответствует безраз­ личной стратификации, а параметр вертикальной диффузии — не­ устойчивой (классы 2 или 3); иначе говоря, максимум осадка фор­ мируется ближе, чем это соответствует состоянию устойчивости, которая при осадках бывает безразличной. Что касается значений Ро и л, они не обнаружили никаких систематических отклонений от общего массива данных.

Более подробно рассмотрены результаты опытов, проведенных при слабых ветрах. В качестве критерия слабого ветра были при­

няты

условия

UH<

1,0

м/с

при Н=2

м; I)п

<

2 м/с при 25 < Н <-

< 97

м; UH<

3,0

м/с

при

97^.Н

< . 300

м.

Их распределение в

зависимости от класса устойчивости и высоты источника во время опытов приведено в табл. 5.23. Звездочками отмечены случаи, удовлетворяющие более жестким условиям по скорости ветра, а именно: [/я<-0,5, 1,0 и 2,0 м/с соответственно тем же градациям высоты.

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

ных величин, связанных с ее значением, для

этих случаев оказал­

ся больше. Тем не менее удалось отметить

некоторые признаки,

по которым результаты опытов заметно и систематически отли­

чаются

от средних результатов данного класса

устойчивости.

В табл. 5.24 представлены средние значения

Н/хт

и Ьу для

этих опытов. Сравнение с данными табл. 5.5. и

5.18

показывает,

что они

превышают аналогичные характеристики соответствую­

щих классов; менее всего это заметно для класса

1, потому что из

16 опытов, соответствующих этому классу, 12 удовлетворяют кри­ териям слабого ветра, а в остальных скорость ветра также не­ велика. Это превышение, иногда достаточно большое, наблюдается, однако, не во всех случаях. Для классов 5 и 6, в частности, оно отмечалось только в переходное время суток.

Таблица 5.23

Распределение условий слабого ветра в зависимости от класса устойчивости и высоты источника

И м

Класс

устойчивости

 

2

3

5

6

1

2

9*

1 *

25

 

 

 

 

49

5

1

1

1

 

73

2*

2

2*

 

 

97

1

 

 

 

1 "

1 ОТ

 

1

 

 

1 1

1

 

3

169

1*

1*

193

1

 

 

1*

301

 

 

 

 

Всего

12

5

3

9

о

Таблица 5.24

Средние значения HjxT (числитель) и Ь у (знаменатель) для условий слабых ветров

11 м

 

Класс

устойчивости

 

1

2

3

5

6

 

 

0,176

_

_

 

 

2

0,26

25-49

0,22

0,15

0,154

0,056

0,24

0,31

0,17

0,040

 

73

0,25

0,16

0,21

0,14

0,11

0,14

 

 

 

97—121

0.20

0,090

 

 

0.50

0,096

 

67058

 

 

 

169-301

0,21

0,082

0,067

0,13

0,12

0,067

0,057

 

Что касается

значений

максимальной

концентрации, то эта ве­

личина

испытывает разброс, значительно

больший,

чем в случаях

с более

сильным ветром. При этом

можно

отметить, что при сла­

бых

ветрах

и неустойчивой стратификации

преобладают аномаль­

но

низкие

значения концентрации (на рис. 4.16 соответствующие

значения

°^Х °^

отмечены):

из

20

таких

результатов

 

 

 

 

Q

 

 

 

 

 

 

14 ниже средних расчетных значений, причем

восемь из них — су­

щественно ниже (на порядок и более), и только три результата

выше, но в пределах общего разброса. В

случаях

слабого

ветра

при инверсиях (классы устойчивости 5 и 6)

плотность

осадка ока­

зывается близкой к расчетной.

 

 

 

 

Таким образом, слабые ветры при неустойчивой

стратифика­

ции или при переходном времени суток имеют тенденцию

увели­

чивать параметры диффузии по сравнению с теми, которые

харак­

терны для соответствующих классов устойчивости,

и

уменьшать

плотность осадка по сравнению с расчетной. При инверсиях

какой-

либо тенденции к изменению этих величин по сравнению со сред­ ними для соответствующего класса устойчивости не замечено.

5.5.2. Влияние гравитационного оседания на параметры рассеяния

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

Соседние файлы в папке книги из ГПНТБ