
книги из ГПНТБ / Гвахария, В. К. Испарение с водной поверхности водоемов Кавказа
.pdfчае настолько мал в сравнении с объемом притекающей и оттекающей воды, что полностью может уложиться в допус
тимую |
величину ошибки измерения притока |
и |
оттока. |
|
|||||||
Например, по расчетам Г. С. Метревели, |
|
испарение |
с |
||||||||
Сионского |
|
водохранилища |
(Грузинская ССР) |
|
равно |
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]. Однако, доводы М. П. Тимофеева против приводимого выше пояснения врвд-ли являются убедительными. Ссылка на несовпадение показаний скорости ветра по анемометру и
флюгеру, несостоятельна, поскольку при построении |
формул |
||||
|
|
ЕМА |
+ В\Ѵг)(в0-ё2) |
|
(20)' |
используются данные только лишь анемометров, |
точность |
||||
которых зависит |
не |
от |
географического пункта, где |
измеря |
|
ют скорость ветра, |
а от |
технических свойств самого |
прибора, |
||
и остается на любых широтах и высотах одинаковой. |
|||||
Как показала практика, коэффициенты |
А и В меняются |
||||
в пространстве, Что говорит, очевидно, о их |
зависимости от |
||||
географических |
условии. |
|
|
||
Очевидно, не является случайностью |
то, что |
формулы |
Б. Д. Зайкова, Браславского—Викулиной и др., составленные по материалам «равнинных» испарителей, дают весьма близ кие между собой значения параметров А и В, тогда как зна чение тех же параметров в формулах, составленных по гор ным странам, группируются вместе (формулы В. К. Давы дова, О. Й. Халатяна, В. К- Гвахария). Для объяснения этото явления следует, очевидно, провести более углубленное исследование.
Существуют эмпирические формулы, в которых разность ео—е20о заменена дефицитом влажности воздуха D (формулы В. К- Давыдова, Б. Д. Зайшва и др.). Однако, использова
ние дефицита влажности воздуха вместо |
разности eQ |
— |
е200 |
||||
нам |
представляется |
не |
совсем 'оправданным. |
Дело |
в том, |
||
что |
дефицит влажности |
Q формируется |
за счет макропро- |
||||
Цесоов, протекающих |
в |
атмосфере, тогда |
как |
разность е0 — |
|||
—е20д. |
определяется термическими условиями |
воды того |
или |
29