книги из ГПНТБ / Воронцов, П. А. Аэрологические исследования пограничного слоя атмосферы
.pdfпри ветрах с водоема 0,66—0,70. С высотой наблюдается вна-
чале рост ——, а затем уменьшение, что связано с профилями
ветра и температуры воздуха в этих пунктах. Суточный ход скорости ветра на маяке Сухо, по наземным наблюдениям за длительный период, приведенным А. И. Григорьевым, также имеет наибольшие скорости ветра в утренние сроки и мини мальные днем.
Джонс (Jones, 1953) приводит анализ двухмесячных шаро пилотных наблюдений над ветром на судовых станциях в Се
верной Атлантике. |
На основании подсчета он |
делает |
вывод, |
что повторяемость |
различных типов изменения |
скорости |
ветра |
с высотой, в частности, убывания скорости ветра, связана со стратификацией атмосферы. В качестве индекса стратифика ции им принята разность температур воды и воздуха над мо
рем. Повторяемость случаев с уменьшением |
скорости ветра |
с высотой над морем увеличивается с ростом |
этого индекса, |
т. е. с ростом неустойчивости. |
на профиль v |
Рассмотрим влияние бризовой циркуляции |
ив береговой зоне Черного моря по материалам шаропи
лотных |
|
наблюдений, проведенных |
автором |
летом |
в |
1936— |
|||||||||||
1937 гг. |
(табл. |
23) |
|
|
|
|
|
|
|
|
|
|
|
|
|
||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
Таблица 23 |
|||
|
|
Величина — при морских бризах, |
Черное море |
|
|
|
|||||||||||
|
|
|
|
|
|
|
|
|
Высота, |
км |
|
|
|
|
|
||
Пункт |
|
|
0,0 |
0,1 |
0,2 |
0,3 |
0,5 |
0,7 |
1,0 |
1,25 |
1,5 |
2,0 |
2,5 |
3,0 |
|||
|
|
|
|
||||||||||||||
Сочи . |
. |
. |
5,3 |
0,45 0,66 |
0,84 |
0,64 |
0,57 |
0,53 0,51 |
0,62 0,87 0,92 1,00 |
|
|||||||
Гагра . |
. |
. |
4,7 |
0,60 |
0,66 |
0,82 |
0,70 |
0,56 0,47 0,58 |
0,88 0,96 1,00 1,00 1,25 |
||||||||
Батуми |
. |
. |
4,2 |
0,62 |
0,68 |
0,89 0,63 |
0,52 |
0,55 0,64 |
0,78 0,90 1,00 1,20 1,38 |
||||||||
Поти . |
. |
. |
4,0 |
0,37 |
0,63 0,88 |
0,78 |
0,63 |
0,75 1,20|1,12 0,83 1,00 1,37 |
|
||||||||
Первые максимумы v и ~~~ отмечены на всех пунктах на |
|||||||||||||||||
уровне 0,2 |
км, |
минимумы тех |
же величин |
совпадают |
с зоной |
||||||||||||
смены |
направления |
ветра на |
противоположное, |
этот |
уровень |
вкаждом пункте определяется мощностью бризового потока.
Если на уровне 2 м отношение —— зависит от местных
условий закрытости горизонта, то с 0,1—0,2 км эта величина на всех пунктах примерно одинакова. Уровень начала геостро фического ветра, судя по профилю ветра, был в Сочи около 2,5 км, в Гагре и Батуми примерно 2 км и в Поти примерно 0,8 км. Число наблюдений было невелико, всего 15—20 по каждому
59
пункту, но они все же позволяют получить типичную для бри зовых ветров картину.
Профили v и ——— над Пахта-Аралом и Голодной Степью.
Средние скорости ветра над полупустыней (Голодная Степь) и орошаемым оазисом (Пахта-Арал), несмотря на близкое рас положение пунктов наблюдений (25 км), резко различны. Вы звано это особенностями термического режима над рассматтриваемыми пунктами, а также влиянием местных ветров Ангренской долины, особенно в ночное время.
Суточный ход v выражен сравнительно хорошо до уровня
1,5 км. Высота обращения ветра около 50 м. Быстрый рост v
отмечен в нижнем слое 50—100 м, |
особенно |
при |
сильных вет |
|||
рах, |
суточный ход |
Vz— приведен в |
табл. 24. |
|
|
|
Величины ог и Яп были подсчитаны ориентировочно по |
||||||
профилю ветра. Приведенные значения |
— дают |
хороший |
||||
пример суточного хода этих величин. |
|
|
|
|||
Профили v и - |
в пересеченной местности. В пересеченной |
|||||
местности нарушается как профиль |
ветра, |
так |
и |
отношение |
||
v* |
. Рассмотрим |
распределение |
v и |
по |
наблюдениям |
в двух пунктах: Кольском заливе и Мурманске, расположен ных на расстоянии 7 км друг от друга (табл. 25). Пункт
наблюдений на берегу Кольского залива находится в относи
тельно узком |
фиорде с |
грядами холмов высотой 120—150 м. |
||||
Над заливом |
в зимнее |
время |
ветер |
направлен вдоль фиорда, |
||
т. е. |
имеет примерно северное |
или |
южное |
направление, при |
||
этом |
происходит сжатие воздушного потока |
с ростом скорости |
в нижних слоях, это хорошо подтверждают данные табл. 25 по
распределению и и максимумам —— в нижнем слое.
Пункт наблюдений в Мурманске расположен на относи тельно ровном месте и имеет нормальный для данного района профиль ветра.
Другой пример местных орографических искажений про
филя ветра приведен в табл. 26.
При боре в Новороссийске сжатие потока воздуха с ростом скорости ветра происходит на Мархотском перевале, однако и в Новороссийске профиль ветра еще сохраняет типичные черты орографического искажения. Если при ветре с моря про
филь v и —— сходен с другими береговыми районами, то при боре максимум v и --^г— будет в самом нижнем слое, причем
скорости ветра у земли могут быть в 2—3 раза больше
60
Таблица 24
V,
Величины — Голодная Степь и Пахта-Арал. Июль 1952 г.
|
|
"п |
|
|
|
|
|
Высота, |
км |
|
|
|
|
|
Часы |
V |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
км |
0,0 |
0,05 |
0,1 |
0,2 |
0,3 |
0,5 |
0,75 |
1,0 |
1,25 |
1,5 |
2,0 |
2,5 |
|
|
|
||||||||||||
|
|
|
|
|
|
Го л одна я |
Степ ь |
|
|
|
|
|
|
|
0 |
9,2 |
0,3 |
0,13 |
0=72 |
0,88 |
0,96 |
1,00 |
0,88 |
0,76 |
0,51 |
0,34 |
— |
— |
— |
6 |
6,8 |
0,5 |
0,25 |
0,58 |
0,67 |
0,92 |
0,98 |
1,00 |
0,82 |
0,72 |
0,64 |
— |
— |
— |
12 |
4,4 |
2,5 |
0,66 |
0,91 |
1,13 |
1,11 |
0,98 |
0,93 |
0,88 |
0,76 |
0,61 |
0,70 |
0,93 |
1,00 |
16 |
5,3 |
2,0 |
0,44 |
0,75 |
0,83 |
0,89 |
0,82 |
0,78 |
0,77 |
0,72 |
0,70 |
0,78 |
1,0 |
— |
20 |
9,4 |
0,3 |
0,10 |
0,70 |
0,81 |
0,92 |
1,00 |
0,91 |
0,87 |
0,78 |
0,63 |
0,50 |
0,43 |
— |
|
|
|
|
|
|
|
П ахта-Арал |
|
|
|
|
|
|
|
0 |
6,0 |
0,5 |
0,10 |
0,81 |
0,88 |
0,92 |
0,96 |
1,00 |
1,02 |
1,06 |
0,95 |
0,83 |
0,67 |
— |
6 |
4,8 |
0,5 |
0,21 |
0,58 |
0,65 |
0,78 |
0,98 |
1,00 |
1,00 |
0,85 |
0,75 |
0,61 |
0,86 |
1,23 |
12 |
3,7 |
2,0 |
0,33 |
0,65 |
0,67 |
0,77 |
0,74 |
0,67 |
0,65 |
0,59 |
0,65 |
0,78 |
1,00 |
1,45 |
16 |
4,8 |
1,5 |
0,21 |
0,35 |
0,38 |
0,41 |
0,42 |
0,55 |
0,64 |
0,75 |
0,92 |
1,00 |
1,02 |
1,29 |
20 |
4,6 |
0,5 |
0,11 |
0,80 |
0,84 |
0,94 |
0,98 |
1,00 |
1,03 |
1,06 |
1,19 |
— |
— |
— |
|
|
Vz |
|
|
|
|
|
|
|
|
Таблица 25 |
|
|
Распределение v и |
Кольский залив и |
Мурманск. |
Январь—февраль 1953 и 1954 гг. |
|
|||||||
|
—. |
|
||||||||||
|
Элемен |
|
|
|
|
|
Высота, |
м |
|
|
|
|
Место |
k |
|
|
|
|
|
|
|
|
п |
||
ты |
2 |
25 |
50 |
100 |
150 |
200 |
300 |
400 |
||||
|
|
500 |
Кольский залив с —
Мурманск
Ветер
V
Vr V vz
S С *
450 |
6,5 |
3,7 |
5,0 |
4,8 |
4,4 |
4,0 |
4,4 |
4,7 |
5,7 |
6,4 |
6,8 |
30 |
— |
— |
— |
0,77 |
0,74 |
0,68 |
0,61 |
0,68 |
0,72 |
0,86 |
0,98 |
1,02 |
— |
840 |
9 |
— |
— |
— |
3,0 |
4,2 |
5,0 |
6,5 |
7,6 |
8,3 |
8,8 |
21 |
— |
— |
— |
|
— |
0,33 |
0,46 |
0,55 |
0,72 |
0,84 |
0,92 |
0,98 |
— |
|
|
|
|
|
V2 |
|
|
|
|
|
Таблица 26 |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
Величины v (первая строка) и — (вторая строка) для Новороссийска |
|
|
|
|||||||||
|
|
|
|
|
|
|
Высота, |
км |
|
|
|
|
vr |
k |
0,0 |
0,1 |
|
0,2 |
0,3 |
0,4 |
0,5 |
0,6 |
0,8 |
1,0 |
1,5 |
|
|
|
Бора (СВ) |
0,8 |
4,6 |
1,9 |
13,2 |
11,6 |
10,0 |
9,3 |
8,4 |
7,3 |
6,6 |
4,6 |
3,2 |
2,0 |
|
— |
— |
— |
2,86 |
2,52 |
2,17 |
2,01 |
1,82 |
1,56 |
1,42 |
1,00 |
0,69 |
0,43 |
С моря(ЮЗ) |
0,8 |
3,9 |
— |
2,5 |
2,8 |
3,0 |
3,2 |
3,4 |
3,6 |
3,8 |
3,9 |
4,0 |
4,2 |
|
— |
— |
— |
0,64 |
0,72 |
0,77 |
0,82 |
0,87 |
0,93 |
0,98 |
1,00 |
1,02 |
1,03 |
скорости геострофического 'ветра. С высотой эти величины очень быстро уменьшаются.
Блэккейдор (Blackadar, 1957) отмечает наличие максимума ветра на высотах ниже 1000 м. В таких случаях максимум v
обычно превышает vr и появление пмакс не сопровождается
существенными изменениями во времени градиента давления. Наиболее отчетливо цмакс наблюдается в ночные часы. В боль шинстве случаев, особенно летом, при значительных суточных
колебаниях температуры высота слоя с wMaKC примерно совпа дает с верхней границей приземной инверсии. Наличие макси мума v оказывает существенное влияние на развитие ночных инверсий; турбулентный обмен имеет место в инверсионном слое и отсутствует выше этого слоя.
Профиль коэффициента турбулентного обмена в пограничном слое
Процесс турбулентного перемешивания имеет весьма боль шое значение в формировании пограничного слоя и предста вляет интерес как для решения прикладных задач, так и для
общих теоретических исследований.
Лучше всего процесс турбулентного перемешивания харак теризуется коэффициентом турбулентного обмена k, изученным
весьма подробно в приземном слое. |
Работы М. И. Будыко, |
Д. Л. Лайхтмана, М. П. Тимофеева |
(1953), А. С. Монина, |
А. М. Обухова (1954) и др. позволяют |
с достаточной точностью |
и сравнительно просто определять k с помощью градиентных наблюдений над распределением v и t° в нижнем слое 1,5—
2 м. По наземным наблюдениям вычисляют k на высоте одного
метра, а затем на других уровнях с помощью различного рода соотношений, например
kz~ kzx~* , |
(30) |
где 8 — параметр, характеризующий |
устойчивость атмосферы |
в нижнем слое 1,5—2 м; z— высота |
в метрах. Такая эстрапо- |
ляция обычно проводится до уровня высоты приземного под слоя h, выше которого величина k принимается постоянной.
Применение этой методики для расчета k в пограничном слое далеко не всегда возможно в связи с рядом особенностей в распределении v и t° с высотой, а также резким изменением спектра вихрей по высоте.
Методика подсчета коэффициента k для пограничного слоя.
Методов непосредственного подсчета k в пограничном слое
значительно меньше. Кроме метода подсчета k по данным ве
личины геострофического ветра и значений температуры воз
духа на земле и на высоте пограничного слоя по Д. Л. Лайхт-
ману [формулы (10) и (12)1, существует другая группа методов подсчета k по профилю ветра.
63
В работах К. И. Кашина (1939) и И. В. Ханевской (19486) для вычисления k было использовано распределение ветра по
высотам.
Д. Л. Лайхтманом (1944, 19476) было выведено несколько формул для вычисления среднего в пограничном слое значения k по профилю ветра. При этом были взяты следующие предпо
ложения: 1) движение стационарно по времени; 2) траектория частиц имеет малую кривизну, поэтому центробежной силой можно пренебречь; За) барический градиент линейно меняется с высотой, 36) барический градиент не меняется с высотой.
Условие За введено для случая, когда горизонтальные гра диенты температуры велики (массы неоднородны), вследствие
чего направление барического градиента сильно меняется с вы сотой и наблюдается левое вращение ветра с высотой.
Условие 36 введено для случаев, когда горизонтальные
градиенты температуры малы (однородные массы) при правом
повороте ветра и росте его скорости с высотой.
При условии За формулу Д. Л. Лайхтмана можно записать в таком виде:
k= ----------------------------------- 2/(lge)2------------------------------------- . |
(31) |
|||||||
| |
1 |
дх р,)1 |
ду Т01) |
' |
дх lp0 |
дх Tol ' J I |
|
|
\ |
|
|
|
dz |
|
J |
|
|
При условии 36 |
£ —2Z(lge)2 |
|
.$2\ |
|||||
|
|
|
|
|||||
|
|
|
\ |
|
dz |
/ |
|
|
Здесь |
/ = 2со |
sinq)—параметр Кориолиса; |
e—основание на |
|||||
турального |
логарифма; |
и и |
v — составляющие вектора |
ветра |
||||
по оси ох и оу; |
др |
др |
дТ |
дТ |
|
|
||
~q~ , —, |
-§■— и |
-----горизонтальные гра |
диенты давления и температуры; R-—газовая постоянная; g — ускорение силы тяжести; То и ро — температура и давление на
данной станции на |
уровне |
будки; ug — составляющая геостро |
||
фического ветра пг |
по оси ох. |
только от широты места, |
||
Числитель этих |
формул |
зависит |
||
а поэтому постоянен |
для |
данной |
станции. Чтобы вычислить |
|
знаменатель формулы |
(31), нужно, |
кроме распределения ветра |
по высоте, знать наземные горизонтальные градиенты темпера туры и давления.
На практике далеко не всегда удается вычислить величину k по формулам Д. Л. Лайхтмана (31) и (32) из-за того, что не все наблюдения дают четко выраженный с высотой правый или левый повороты и рост скорости ветра, а также не всегда уда
ется надежно вычислить горизонтальные градиенты давления и
температуры по синоптическим картам. Кроме того, для подсче
64
та k по этим методам приходится пользоваться сглаженными кривыми распределения ветра по высоте.
Третий способ подсчета коэффициента турбулентного обмена k в пограничном слое основан на использовании пуль
саций различных физических величин с помощью так называе мых структурных формул.
Структурные методы определения k изложены в работах
Гессельберга (Hesselberg, 1929), Эртеля (Ertel, 1930),
Е. С. Ляпина (1948), А. С. Дубова (1949) и др.
Метод Гессельберга основан на использовании регистраций
вертикальной скорости в движущемся воздушном потоке.
Формула Гессельберга имеет вид
«z— 2 ’
где w—средние значения вертикальной скорости на уровне z,
т—среднее время существования w одного знака.
Этот метод был использован при обработке данных полетов Леттау и Швердфегера на аэростатах. Методика получения исходных данных оказалась весьма сложной, и метод не полу чил широкого распространения.
Метод Эртеля более удобен. Для определения k нужно знать пульсации и средний вертикальный градиент в данной
точке любого метеорологического элемента, который можно рассматривать как субстанцию, пассивно переносимую пото ками воздуха. Однако метод Эртеля также обладает рядом не достатков, которые приводят к резкому расхождению величин k. Более подробно критическое изложение метода дано в рабо
тах Л. Ф. Щербаковой (1949) и И. С. Борушко (1949). |
|
||||
По методу |
Е. С. Ляпина, необходимо знать пульсации вер |
||||
тикальной и |
горизонтальной |
составляющих |
скорости |
ветра |
|
в фиксированной точке пространства на любом уровне. |
|
||||
Формула Ляпина для расчета k |
имеет вид |
|
|||
|
А __ |
W-tV |
|
|
(34) |
|
~' ”2AtF ’ |
|
|||
|
|
|
|||
где w и До — средние значения пульсации вертикальной |
и во |
||||
ризонтальной составляющих скорости ветра; |
т — среднее |
время |
|||
существования |
w одного знака; v —■ средняя |
скорость воздуш |
|||
ного потока. |
|
|
|
|
|
Обобщение метода Ляпина для случая подвижного прием |
|||||
ника скорости |
ветра было сделано |
А. С. Дубовым. Для |
этой |
цели были использованы записи вертикальных ускорений центра тяжести самолета, который являлся и приемником ско рости ветра.
В настоящее время по ускорению самолета можно относи тельно точно определять только величины w и почти невоз можно измерить Ди.
5 Заказ № 345 |
65 |
А. С. Дубов делает допущение о изотропности структуры
ветра, т. е. допускает —— = 1 и получает формулу для k
|
|
kz |
~2~ |
|
(35) |
|
Принимая |
w = ЬДп, где |
Ап— перегрузки |
самолета |
в до |
||
лях, b — некоторый |
коэффициент, |
т — время |
сохранения |
пере |
||
грузок одного |
знака, |
формула (35) |
для расчета k приобретает |
|||
вид |
|
, |
b^nvz |
|
|
|
|
|
|
|
(36) |
||
|
|
к'=— |
|
|
Для приземного слоя воздуха влияние земной поверхности обычно приводит к резкой анизотропности структуры воздуш-
кого потока и отношение W ..- значительно меньше единицы.
Согласно ряду исследований влияние земной поверхности про стирается до высоты 30—40 м, выше которой при изотермиче ских условиях приближенно достигается изотропность и выпол
няется условие — 1. По исследованиям автора, отношение
■- = 1 будет при у = 1°- Устойчивость стратификации резко
сказывается на величине w и продолжительности существова
ния w одного знака..
Величина k, по формуле А. С. Дубова, в инверсиях должна получаться несколько завышенной, а при неустойчивой страти фикации заниженной. Величина k будет также зависеть от масштаба вихрей турбулентных движений. В соответствии с тем, какой участок спектра турбулентных движений рассмат ривается, будут получаться и значения k. В дальнейших под
счетах k (будем пользоваться короткопериодными участками спектра вихрей, размеры которых имеют порядок около 100 м. Расчет k по полной эпюре вихрей дал бы очень большие зна чения, что связано с влиянием процессов большого масштаба. Автором структурная, формула Ляпина иногда применялась для расчета k по данным регистрации структуры воздушного потока на привязном аэростате.
Расчет коэффициента турбулентности для некоторого слоя, постоянного по высоте и времени, по амплитудам температуры воздуха, известным для нескольких высот этого слоя, был пред ложен Л. Р. Орленко (1955). Суточный ход температуры воз духа на высоте z можно представить в виде синусоиды, и если
ограничиться только одной гармоникой, то |
|
А^Аое~(г~гв)^^, |
(37) |
где Аг и Ао — амплитуды температуры воздуха |
на высотах z |
66
и Zo; k — коэффициент турбулентности; т — период колебаний,
для суток т = 24 часа = 86 400 сек.
Логарифмируя и дифференцируя это выражение, получим
где —представляет собой тангенс угла наклона прямой
tg а и системе координат (In Az, z), тогда
Величины коэффициента k в пограничном слое. Подсчетом ве
личин |
коэффициента |
турбулентного обмена в |
пограничном |
слое, |
по методике |
Д. Л. Лайхтмана, много |
занималась |
М. П. Чуринова (1954, 1953 а, б). Некоторые эксперименталь ные определения k приведены также в работах В. А. Девятовой
(1957), Л. Р. Орленко (1955), автора (1952 г, 1956 в) и др.
Было установлено, что средние для всего пограничного слоя величины k растут от утра к полудню обычно в 2—4 раза
и к вечеру убывают. Величины k в полуденные часы в среднем
находятся в пределах 20—40 м2/сек. Согласйо исследованиям М. П. Тимофеева (1950) и ряда других авторов, установлено, что в дневные часы значения на уровне 1 м имеют величины 0,2—0,4 м2/сек. Эти величины примерно в 100 раз меньше k для
пограничного слоя и, следовательно, слой, где k растет, можно ориентировочно считать в 100 м.
По данным Л. Р. Орленко (1955), средние за сутки значе ния k в слое 100—400 м в Пахта-Арале в июле 1952 г. соста вили 16,8 м2/сек.
Летом k примерно в 2—4 раза больше, чем осенью и зимой.
Например, для Долгопрудной, по работе В. А. Девятовой (1957), наблюдались такие величины k за отдельные сезоны года:
Сезон |
Зима |
Весна |
Лето |
Осень |
|
k м2/сек. |
14,1 |
12,9 |
23,2 |
|
9,5 |
Значения kz, |
подсчитанные по |
структурным |
формулам, |
дают возможность получить распределение kz по высоте, а не
только средние значения k для всего пограничного слоя.
Для подсчета величины kz автор использовал регистрации перегрузок Ап, испытываемых в неуправляемом полете самоле том ПО-2 над некоторыми подстилающими поверхностями.
Расчет производился |
по |
структурной формуле |
Дубова (36) |
|
(рис. 9). |
|
|
высоте |
150— |
Максимальные величины kz наблюдались на |
||||
200 м в Пахта-Арале |
и |
Каменной Степи и 100 |
м в |
районе |
5* |
67 |
■пос. Воейково. Считая, что рост k с высотой при монотонном из менении t и v будет происходить только в приземном слое, полу
чается, что высота слоя может быть не 30—50 м, как это при
нимается рядом исследователей, а |
значительно |
больше — до |
|||
0,1—0,2 км. Интересно также наличие |
больших |
значений kz |
|||
в зимнее время над пос. Воейково. |
|
(36), |
получились |
явно |
|
Величины kz, по формуле Дубова |
|||||
заниженными, что вызвано влиянием коэффициента Ь. |
мала, |
||||
Принятая в подсчетах величина |
Ь = 3,62, |
очевидно, |
еенеобходимо увеличить.
Вычисление величины k было сделано также по материалам
Рис. 9. Распределение k над некоторыми
подстилающими поверхностями.
1 — орошаемый оазис, июль; 2 — Голодная Степь, июль; 3—массив лесных полос в Каменной Сте пи, июнь—июль; 4 — поле, март—апрель, пос. Воей ково, 5 — то же.
подъемов на привязном аэростате прибора с регистрацией структуры ветра. Для этой цели была использована формула
Ляпина |
(34). |
|
|
|
|
Таблица 27 |
|
|
|
Величины kz, пос. Воейково |
|||
|
|
|
|
|||
|
|
|
|
Высота, |
м |
|
Сезон |
|
2 |
100 |
200 |
300 |
п |
|
|
|||||
Лето . . . |
1,0 |
27,2 |
20,2 |
20,6 |
31 |
|
Осень . . |
. |
26,1 |
24,0 |
18,8 |
48 |
|
Зима . . |
. |
0,88 |
33,1 |
29,5 |
29,9 |
22 |
Весна . . |
. |
1,2 |
31,8 |
28,0 |
31,0 |
27 |
68