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

книги из ГПНТБ / Гвахария, В. К. Испарение с водной поверхности водоемов Кавказа

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

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

тимую

величину ошибки измерения притока

и

оттока.

 

Например, по расчетам Г. С. Метревели,

 

испарение

с

Сионского

 

водохранилища

(Грузинская ССР)

 

равно

4,73

млн. м3

за год; годовой объем суммарного

притока в этомвоДо.

хранилище равен 377,4 млн. м3 , т. е. испарение

 

"от

притока

составляет всего около 1,5%, остальное стекает

[41].

 

•Как известно, допустимая ошибка

определения

стока

горных

рек,

характеризующихся

чрезвычайно

 

изменчивым

режимом и неустойчивым

руслом

(живым сечением), лежит

в пределах

± 5 % , а иногда

и более того. Конечно,

определять

испарение

в

этих условиях

как

остаточный

член

баланса,

недопустимо. В этом случае Г. С. Метревели испарение было подсчитано по известной формуле А. П. Браславского и 3. А.

Викулиной

[1].

 

 

 

 

 

 

 

 

 

 

2.

Метод

испарителей

 

 

 

Метод испарителей, т. е. метод определения

испарения с

©одной поверхности посредством

специальных

 

сосудов-испа­

рителей

основан

'на принципе

применения /метода водного

баланса,

при условии,

когда

все

члены этого

баланса

опре­

деляются

с

высокой

точностью,

чего можно

достичь

лишь

в искусственных

условиях.

 

 

 

 

 

Уравнение водного баланса

испарительного

сосуда

чрез­

вычайно примитивно: оно состоит лишь из трех членов для

любого отрезка

времени

 

 

 

 

 

 

 

 

 

 

 

E = âU-{-X,

 

 

 

 

 

(4)

где

Е—испарение,

Ш—изменение

уровня

(объема) воды

в

ис­

 

парителе

между

сроками

измерения

и

К—количество

 

осадков,

выпавших на поверхность

испарителя

в

том

 

же промежутке времени.

 

 

 

 

 

 

 

 

Очевидно, что в 'случае отсутствия осадков

высота

слоя

испарившейся воды будет равна высоте слоя

убывшей

из

испарителя воды, т. е. вся убывшая из

испарителя

вода

бу­

дет

израсходована лишь

на .испарение. По

знаку

испарение

(Е)

и изменение объема воды в сосуде ( A U) противооолож-

20

ны: убыль будет соответствовать испарению, а прибыль, как правило, конденсации1 .

Первые попытки определения испарения с помощью сосудов-испарителей, относятся к последним десятилетиям прошлого века. В 1880—1890 годах начинается более или менее широкое применение чаш-испарителей (эвапорометров) Внльда. В силу малого количества воды в них, а также ис­ кусственных условий, создаваемых вокруг испарителей ме­ теорологическими будками (затененность, 'Непродуваемость и т. д.), показания испарителей Вильда были весьма дале­ кими от фактических.

Лишь в самом конце прошлого века

появились

более

совершенные приборы для измерения испарения

с

водной

поверхности.

Это испарители

Лермонтова-Лгобослаівского,

конструкция

которых в какой-то мере послужила

прототи­

пом современных испарителей, широко известных

под наз­

ванием ГГИ—3000. Правда,

испарители

Лермонтова-Л юбо-

славокого были плавучими,

но

от этого,

конечно,

 

суть их

не меняется и принцип их работы вполне современен и пов­ торяется в ГГИ—3000, будь они установлены на суше либо на воде.

Сеть испарителей, разбросанных на территории Совет­ ского Союза, довольно обширна: около 300 малых испаромеров (ГГИ—3000) и свыше 60 испарительных бассейнов с площадью 20-кваідратных метров и более (имеется несколь­ ко испарителей с S = 100 м2 ) составляют основную сеть, на­ ходящуюся в ведении ГУГМС. Кроме этого, существует зна­ чительное число испарительных установок, находящихся в ведении других организаций, заинтересованных в определе­ нии величин испарения на тех или иных конкретных объек­ тах. Большинство испарителей ведомственного подчинения

1 В случае, когда выпадают

осадки, изменение объема между срока­

ми измерения, также может оказаться положительным,

т. е.

будет от­

мечена

не убыль воды в бассейне, а прибыль.

Однако,

это

не значит,

что в этом промежутке времени

испарения

не

было:

при оштенсиадых

ливнях

за час, например, может

выпасть столько осадков, что слой, об­

разованный ими в испарителе, не

попарится

полностью за «ее

последую­

щие II

часов (при двухсрочном

измерении

испарения).

 

 

21

•Страдает одним общим недостатком — они нестандартны, т. е. испаряющая площадь, глубина, высота бортика над во­ дой и над землей этих испарителей не соответствуют габари­ там стандартных 20-метровых бассейнов2 .

Наряду с наземными .испарителями в последнее время расширяется масштаб применения плавучих испарителей, ос­ нащенных в основном ГГИ3000, хотя предпочтительнее, ко­

нечно, устанавливать на плотах большие

бассейны, показа­

ния которых можно считать идентичными

с величиной испа­

рения непосредственно с водоема3 .

 

 

Преимущество испарителей, установленных посреди вод­

ного

объекта

(т. и. плавучий иопаритель)

перед установлен­

ными

на суше

(на берегу водоема, либо в

отдалении от

не­

го), очевидно:

первые испытывают на себе воздействие

тех

же метеорологических, аэродинамических и водно-термичес­

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

объектом,

тог­

да как вторые находятся под воздействием

гидро-метеоро-

логическнх элементов, отличающихся от первых

как коли­

чественно, так

и качественно. Поэтому и величина испарения

по наземным

испарителям

не будет совпадать

с

испарением

с водоема.

 

 

 

 

 

Изучение

испарения с

водной поверхности

методом

ис­

парителей широко развито и за пределами СССР. В зару­

бежных

странах

в

качестве иопарнтелей используют

такие

же

сосуды

(испарители) как

и в

Советском

Союзе, однако,

их

габариты

и способ установки

несколько

отличаются от

наших.

Принцип

действия

идентичен.

 

 

 

Следует

подчеркнуть,

что

ряд

технических причин

дела­

ет

несовершенным

показания

даже стандартных испарите­

лей (возвышение бортика над водной поверхностью, метод

измерения слоя

испарившейся ©оды с

помощью

бюретки

и т. д.).

 

 

 

2 Подробно на

вопросы о «нестандартных»

бассейнах

и о возмож­

ности использования результатов наблюдений на них говорится в сле­

дующей

главе.

3 На

некотерых водоемах Советского Союза установлены плавучие

иопарители с большой площадью испаряющей поверхности. Среди них следует отметать Валдайский и .Кингирокий плавучие бассейны с S = 20 м2. 22

Как говорилось выше, большое различие в метеорологи­ ческих, аэродинамических и водно-термических условиях в районе берегового испарителя и непосредственно над водое­ мом, не дает возможности использовать показания берегово­ го испарителя для подсчета 'испарения с водоема непосред­ ственно. В показания береговых установок следует вводить соответствующие коррективы в виде поправочных коэффи­ циентов Къ, Ks, указанных в главе VI настоящего труда.

Сложнее обстоит дело с испаромерамн ГПТ—3000, в по­ казания которых приходится вводить редукционные коэф­ фициенты, величина которых сильно колеблется как во вре­ мени, так и в пространстве и который зависит от многих фак­ торов. Главными из «их являются климатические условия, широтное и высотное положение испарителя, а также время года.

3. Метод теплового баланса

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

 

Тепловой баланс «деятельного слоя» водоема описывает­

ся

следующим уравнением:

 

 

 

 

 

 

 

R = LE+P+B,

 

 

 

(5)

где

R — радиационный

баланс; LE—затраты

тепла на

испарение

 

(£—испарение,

L—теплота

парообразования);

Р—тепло­

 

обмен

с атмосферой и В—теплообмен с

нижележащими

 

слоями

воды.

 

 

 

 

 

 

В уравнении (5)

неизвестны

Е, Р, В

и R.

Ясно,

что для

его решения следует определить радиационный баланс, ве­ личину теплообмена с атмосферой и с нижележащими слоями.

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

23

[29, 51 и др.], с помощью которых R можно определить с довольно высокой точностью. Величина 'радиационного ба­

ланса может быть определена как расчетным

путем,

так и

измерена

непосредственно.

 

 

 

 

 

 

 

 

 

Сложнее обстоит

дело

с определением остальных

двух

членов уравнения

и В).

 

 

 

 

 

 

 

 

 

Теплообмен с атмосферой (Р) можно рассчитать через

испарение

по формуле:

 

 

 

 

 

 

 

 

 

 

 

£ = - 0 r 6 2 2

Çïî-

 

i l l ;

Р = - С

p / С т — ,

 

(6)

 

 

 

Р

 

дг

 

 

 

 

dz

 

 

где в, р, р , Г—влажность,

плотность,, давление

и температура

воздуха, Ср —теплоемкость

воздуха

при пост, давлении, Ке "

Кт — коэффициенты

обмена

влаги

и тепла,

Z—имеет

направле­

ние вверх.

 

 

 

 

 

 

 

 

 

 

 

 

 

По уравнению

(6) легко

определить

LE+P,

оно равняется

LE+P^LE

 

Р

\

 

I

'

 

К

т С »

р Т"

 

 

( 1 + — )

= LE

1

 

 

ö

z

 

 

 

 

L E

I

 

V

 

'

L .

0,622

P

dz

 

 

 

 

 

 

 

 

 

 

 

 

 

 

=LE

Ii+-1*1—

 

. — ] .

 

 

(7>

 

 

 

V

 

0,622 L

' - 1

 

 

 

 

При решении этого уравнения принималось, что коэффициенты обмена влаги и тепла, согласно Боуену, равны е = Кт), а вместо дТ и де приняты их конечные значения ДГ и Ае.

Подставляя значение LE+P из (7) в уравнение теплового

баланса (5)

 

 

R = LE ( 1 +

_ ? " Л _

. ^L) +В\

V

0,622 L

Ае)

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

Ь Е = , . %р

д г '

<8>

0,622 L Ае

Касаясь сути формулы (8), Ä. Р. Константинов- [3*2] указывает на ее неточность, и уточняет ее введением' коэф­ фициентов ат и ае после чего она примят следующий вид:

24

1 +

CpP

aT

AT

( У )

0,622 L

ae

Ae

 

Для окончательного решения данной формулы следует рассчитать величину теплообмена с нижележащими слоями воды (В).

Расчет величины В по зависимости

 

 

 

 

 

 

В=-К*р*

С* —

,

 

(10)

 

 

 

 

 

 

 

 

 

дг

 

 

 

где

К*,

р*,

С*

и

Т* — коэффициент

обмена,

плотность,

удельная

теплоемкость

и-

температура

воды,

соответст­

венно,

связан с трудностями

определения

коэффициента об­

мена

(К*) и градиента температуры у поверхностных слоев.

 

М. П. Тимофеев

рекомендует

для большого

интервала

времени

принять

BœO,

т. к. средняя

температура вводы за

большой

период

меняется

мало.

Коль

скоро

принимаем

В ça О

то, очевидно,

что выражение

(9)

можно

переписать

следующим

образом:

 

 

 

 

 

 

 

 

0,622 L

ае

Ае

 

 

 

 

 

Дальнейшие расчеты зависят уже от- наличия

материа­

лов наблюдения над AT к Ае,

а также

величины

радиаци­

онного

баланса..

 

 

 

 

 

 

 

Касаясь преимущества расчета испарения методом теп­

лового баланса перед определением по

методу

водного

ба­

ланса,

А. П. Браславский и К. Б". Шергина [2] пишут,

что

«определение

величины испарения

как

остаточного

члена

уравнения теплового баланса, деленного

на окрытую

тепло­

ту испарения,

рационально потому,

что удельный

вес

затрат

тепла на испарение в этом уравнении1 во много раз (в сред­ нем в 50) больше, за счет высокого значения величины скры­ той теплоты испарения», тогда, как абсолютная величина ис­ парения может оказаться соизмеримой- с полрешностью из ­ мерения элементов водного баланса, при: определении испа­ рения этим последним методом.

25-

•А. Метод турбулентной диффузии

Метод исходит из предположения, что вся влага, испа­ рившаяся с подстилающей поверхности, отводится от послед­ ней лишь вертикально, за счет турбулентной диффузии па­

ров в

воздухе.

 

 

Перенос водяных паров через .какую-то плоскость, рас­

положенную над испаряющей поверхностью,

в этом случае

можно

описать

следующим уравнением:

 

 

 

Е = Р К е ^ - ,

(12)

 

 

dz

 

где -р—плотность

воздуха, Кекоэффициент турбулентной диф-

фузии,

dq

 

влажности зоз-

и ——средний вертикальный градиент

 

dz

 

 

духа.

А. Р. Константинов [32] считает, что данная формула не точно описывает процесс испарения при турбулентной диф­ фузии и предлагает ввести в формулу (12) дополнительный •множитель ае , который, по автору, учитывает различие «эф­ фективного» и среднего значения градиента влажности и влияние архимедовых сил на интенсивность вертикального обмена воздушных масс. Тогда формула (12) принимает

.следующий вид:

E^cKaäd-l.

(13)

 

dz

.Для расчета а е предлагается зависимость следующего вида

а е = 1 + 0 г 7 2 ( К 1 - 2 8 З Д - 1 ) ,

(14)

где Ri—число Ричардсона, а Z0—.шероховатость подстилающей поверхности.

Для водной поверхности Z0

принимается равным

0,01 см, а

Ri рассчитывается по формуле

 

 

. ^ , 0 » = - 0 , 0 7 8

Т ° ' 2 ~ Т ^ ,

(15)

 

(W2,0-WQ,2f

 

если

значение

ветра измерялось в двух точках (в точках 0,2

и 2,0

метров

над испаряющей поверхностью-), 'и по формуле

.26

 

 

 

# i 1 > 0

= -0,025 I n 3 ' —

l^-Ull,

 

 

(16)

 

 

 

 

 

^0

 

l V 1,0

 

 

 

когда ветер измерялся

в одной

точке, на высоте

1,0

метра

над

испаряющей поверхностью.

 

 

Т обозначена

 

 

В формулах

(15)

и

(16)

через

температура

воздуха, а через W ветер.

 

 

 

 

 

 

 

Для

того,

чтобы рассчитать

испарение по

формуле

(13),

 

 

 

 

 

 

 

do

 

 

 

 

следует

иметь

натурные

данные

— ,

т. е. знать градиент влаж-

 

 

 

 

 

 

 

dz

 

 

 

(К)

ности воздуха, а также коэффициент

турбулентного

обмена

и коэффициент

ае.

 

 

 

 

 

 

 

 

 

Для

расчета

'коэффициента

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

построенную на основе теории Прандтля-Каірімана, которая имеет следующий вид:

/C = x'-Z2

,

(17)

 

dz

где х —• безразмерная постоянная н, по данным Никурадзе, равна 0,38. По исследованиям А. Р. Константинова, для рав­ новесных условий атмосферы значение х совпадает с вели­ чиной, полученной Никурадзе для гладких труб.

Складывается впечатление, что формула (13) указывает наиболее простой путь для определения величины испарения. Однако при близком рассмотрении становится ясным, что современные методы определения К и ас , вернее методы измерения метеоэлементов, входящих в расчетные формулы (14) и (17), we являются настолько совершенными, чтобы обеспечить приемлемую точность рассчитываемого испаре­ ния.

5. Метод эмпирических формул

Попытки рассчитать испарение с помощью формул, в которые входят гидрометеорологические факторы, обусло­ вливающие испарение в природе, предпринимались неодно­ кратно.

В настоящее еремя формулы именно этого вида явля­ ются наиболее распространенными и расчеты испарения с водных объектов ведутся, в основном с их помощью.

27

Большинство формул

берет

начало с данного еще в на­

чале XIX в. Дальтоном

выражения

 

 

 

Е = С(е0г),

 

 

 

(18)

где С — .коэффициент, численно

равный

скорости

испарения

•при разности упругости

водяного

пара,

равной еди­

нице.

 

 

 

 

 

Зависимость (18) характеризует испарение

в природе с

водной поверхности при спокойной атмосфере.

 

 

Однако в природе условия

испарения редко

являются

столь идеальными и расчет испарения по

(18) не всегда да­

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

чета

испарения с водной

поверхности

имеют

следующий об­

щий

вид:

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

£

=

(

«

-

О

9)

где f(Wz) — A-\-B\Vz

является

т. н.

«ветровым

фактором»,

т. е.

показывает

зависимость испарения от движения

воздушных

масс,

от ветра. А и В являются

параметрами,

которые

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

эмпирически на основании

данных

испарительных

бассейнов.

 

В зависимости от того, в каких

физико-географических

условиях

располагаются

испарительные

бассейны, по

дан­

ным которых найдены значения А и В, числовые

значения

этих

параметров

различны. Например,

числовое

значение

параметра

А колеблется

в

пределах

0,13—0,3,

а

параметра

В — в пределах

0,23—0,72 у разных

авторов.

 

 

 

 

Следует отметить, что возраста.™е параметра

А сопро­

вождается

уменьшением

значений

параметра

В, что приво­

дит к такому положению, что почти

все формулы

этого

типа

в диапазоне скоростей ветра 3—5 .м/еек дают близкие вели- ІЧИ.НЫ рассчитанных испарений, тогда как за пределами этих

скоростей

расхождения между

ними увеличиваются.

Таким

образам,

становится

ясным, что для территории со скоростя­

ми ветра

менее 3 м/сѳк, или более 5 м/сек, следует находить

свои, региональные

значения (коэффициентов

А и В.

 

Формулы этого

типа, «при

всем своем

несовершенстве,

отвечают

ясной механической

модели: Коэффициент

А при

28

 

 

 

 

 

.отсутствий ветра й разности упругости водяного пара равной •единице, соответствует тому переносу водяного пара с ис­ паряющей поверхности, который происходит в результате подъема вверх масс воздуха, принявших -испарившуюся вла­ гу и оказавшихся поэтому более легкими, чем вышерасполо­ женные более сухие массы. Коэффициент В учитывает влия­ ние динамической турбулентности на испарение» [24].

Ряд авторов не соглашается с таким объяснением физи­ ческой сути параметров А и В, в особенности параметра А [64]. Однако, доводы М. П. Тимофеева против приводимого выше пояснения врвд-ли являются убедительными. Ссылка на несовпадение показаний скорости ветра по анемометру и

флюгеру, несостоятельна, поскольку при построении

формул

 

 

ЕМА

+ В\Ѵг)(в02)

 

(20)'

используются данные только лишь анемометров,

точность

которых зависит

не

от

географического пункта, где

измеря­

ют скорость ветра,

а от

технических свойств самого

прибора,

и остается на любых широтах и высотах одинаковой.

Как показала практика, коэффициенты

А и В меняются

в пространстве, Что говорит, очевидно, о их

зависимости от

географических

условии.

 

 

Очевидно, не является случайностью

то, что

формулы

Б. Д. Зайкова, Браславского—Викулиной и др., составленные по материалам «равнинных» испарителей, дают весьма близ­ кие между собой значения параметров А и В, тогда как зна­ чение тех же параметров в формулах, составленных по гор­ ным странам, группируются вместе (формулы В. К. Давы­ дова, О. Й. Халатяна, В. К- Гвахария). Для объяснения этото явления следует, очевидно, провести более углубленное исследование.

Существуют эмпирические формулы, в которых разность ео—е20о заменена дефицитом влажности воздуха D (формулы В. К- Давыдова, Б. Д. Зайшва и др.). Однако, использова­

ние дефицита влажности воздуха вместо

разности eQ

е200

нам

представляется

не

совсем 'оправданным.

Дело

в том,

что

дефицит влажности

Q формируется

за счет макропро-

Цесоов, протекающих

в

атмосфере, тогда

как

разность е0

—е20д.

определяется термическими условиями

воды того

или

29

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