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

книги из ГПНТБ / Воронцов, П. А. Аэрологические исследования пограничного слоя атмосферы

.pdf
Скачиваний:
17
Добавлен:
30.10.2023
Размер:
23.11 Mб
Скачать

при ветрах с водоема 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 и в нижнем слое 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

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