
книги из ГПНТБ / Бызова, Н. Л. Рассеяние примеси в пограничном слое атмосферы
.pdf5.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 и |
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а |
|
|
|
|
|
- |
|
|
|
|
|
|
ИЭМ, |
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 и коэффициента диффузии К