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

Физика вод суши by Винников С.Д., Викторова Н.В. (z-lib.org)

.pdf
Скачиваний:
1
Добавлен:
19.02.2024
Размер:
9.8 Mб
Скачать

схождения, предполагая, что они ограничивают ячейки конвекции.

Выполненные на Ладожском озере подробные исследования пока­ зали, что при глубине воды 8 м расстояние между линиями схож ­ дения d a 3 м, а при глубине 60 м d a 35 м, т. е. расстояние d уве­ личивается с глубиной водоема. Глубина же проникновения цир­ куляции растет с увеличением скорости ветра. По имени ученого,

впервые описавшего этот вид конвективного течения, в литературе закрепился термин «циркуляция Ленгмюра».

I Таким образом, циркуляция Ленгмюра - это результат плот­ ностной неустойчивости, возникающей при охлаждении поверхно-

I стного слоя воды под действием ветра.

Плотностная конвекция и ветровое перемешивание в стоя­ чих водоемах являются причинами образования на некоторой глу­ бине слоя температурного скачка и расслоения их водных масс на три зоны (рис. 7.11): эпилимнион (верхняя зона), металимнион

(средняя зона, или слой температурного скачка) и гиполимнион

(нижняя застойная зона).

Рис. 7.11. Схема ветрового перемеш ивания воды.

1 - распределение температуры воды до воздействия ветра, 2 - распределение температуры воды после ветрового воздействия, 3 - распределение плотности воды до воздействия ветра, 4 - распределение плотности воды после ветрового воздействия.

Описанный процесс конвекции в чистом виде наблюдается в водоемах больших размеров в плане при относительно постоян­ ной глубине. Реальные ж е водоемы ограничены в плане, а глубина

I их уменьшается до нуля у берегов. В этих водоемах при развитии конвекции возникают конвективные течения, схематически пока­ занные на рис. 7.12.

Рис. 7.12. Схема конвективных течений:

а) при охлаждении водоема; б) при нагревании водоема.

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

Итак, установление условий возникновения упорядоченной конвекции Бенара, Ленгмюра и др. позволило определить ее зна­ чение в числе других факторов, формирующих деятельный слой водоема (эпилимнион). При этом необходимо также отметить важность учета формирования в водоеме и водотоке зон с повы­ шенной концентрацией растворенных солей и присутствие раз­ личных взвесей (минеральных, органических и др.), обусловлен­ ных конвективным водообменном, при отборе проб воды на хими­ ческий анализ, мутность и на наличие загрязнителей. Например,

вслучае циркуляции Ленгмюра наиболее загрязненная вода будет

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

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

7.4. Конвективный водообмен в устье реки

Явление проникновения морских вод в устьевую зону реки по ее дну (так называемый клин соленых вод) издавна привлекало внимание многих исследователей. Оно характеризуется конвек­ тивными движениями значительных масс воды в этой зоне реки.

В отдельных случаях длина этого клина может достигать сотни

222

километров. Сразу же отметим, что основными факторами, опре­ деляющими наличие в устье реки клина соленых вод, являются течение пресной воды в реке и соленость воды моря. Следователь­ но, увеличение стока реки (расхода воды) приведет к уменьшению его длины и, наоборот, уменьшение стока - к увеличению длины.

Ярким примером последнего может служить р. Днепр, сток кото­ рой уменьшился после ее зарегулирования водохранилищами, что привело к более глубокому проникновению соленых Черномор­ ских вод в устье реки, вплоть до г. Херсона. При этом из-за увели­ чения солености воды ухудшились ее технические свойства и воз­ можность использования в хозяйственных целях, а также погибла вся речная биомасса в этой зоне.

Следует также отметить, что от проникновения морской во­ ды в устье реки (лимана, эстуария) зависит ее скоростной, терми­ ческий, ледовый и биологический режимы. Однако механизм об­ разования клина соленых вод полностью еще не раскрыт. Поэтому задача по-прежнему считается актуальной.

С целью решения задачи о клине соленых вод предпринима­ лись экспериментальные и теоретические исследования. П одроб­ ные экспериментальные исследования в лабораторных и натурных условиях (на р. Миссисипи) выполнены Д.Г. Кейлеганом. Теоре­ тическим путем решить эту задачу пытались многие исследовате­ ли. В их разработках исходными данными всегда служили законо­ мерности для неустановившегося движения жидкости, что предпо­ лагает нахождение решения для нестационарного клина соленых вод. Однако для многих рек с достаточной степенью точности его можно считать стационарным. Поэтому введение в рассмотрение неустановившегося движения жидкости часто является неоправдан­ ным. Тем более, если учесть, что объем счетных работ и сложность расчетов характеристик потока при неустановившемся движении жидкости во много раз больше, чем при установившемся.

7.4.1. Расчет стационарного клина соленых вод устьевой зоны реки при открытой водной поверхности

Будем рассматривать двухслойное установившееся движе­ ние воды в устьевой зоне широкой реки, предполагая, что течение происходит между прямолинейными и параллельными берегами.

223

Это условие позволяет нам считать его характеристики по ширине реки неменяющимися, т.е. рассматривать плоскую задачу. В верх­ нем слое реки течение под действием силы тяжести направлено в сторону моря, а в нижнем, обусловленное разностью плотностей соленой морской и пресной речной воды, - в обратную сторону

(рис. 7.13). В результате турбулентности и диффузии в потоке мас­ сы более тяжелой (соленой) воды нижнего слоя при своем движе­ нии переходят в верхний слой потока. Следовательно, верхний слой состоит из масс воды реки и моря ( q v + q E03E).

Рис. 7.13. Схема проникновения морских вод в устье реки: 1 - линия раздела по­ токов речной и морской воды (линия нулевых скоростей); 2 - эпюра скорости.

При названных условиях для любого сечения рассматривае­ мого участка реки можем записать уравнение баланса расхода воды:

н

 

Jv d z = qp

(7.29)

О

 

или

 

9н+<7возв+?р=?р

С7 -30)

и уравнение баланса количества солей:

 

н

 

\ v S d z = 0

(7.31)

о

 

или

 

+ (<7р + ?возв)^в = 0 ,

(7.32)

где Я - глубина; v - скорость течения; q n , qB03Bи q

- элементар­

ные расходы водысоответственно нижнегослоя,перешедшего из

224

нижнего в верхний (расход возвратного потока) и реки; z — ось прямоугольных координат, направленная вверх; S n и S B - средняя

соленость воды (кг/м3) нижнего и верхнего слоев.

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

зований.

Тогда уравнение (7.32) примет вид

v HhSH + (qp - v tth )S B = 0 ,

(7.33)

где vH и h - средняя скорость течения и толщина нижнего слоя.

Исходя из предпосылки о распределении скорости течения по глубине по прямой, можно записать

v = - ^ ----- ( г - А ) .

Л

(7.34)

Н- 2 h H

h

Если в (7.34) подставим z = — , то найдем выражение для

определения средней скорости нижнего слоя:

Qnh

V « = - r r Z

( 7 -3 5 )

А после совместного решения (7.33) и (7.35) найдем выра­ жение для определения средней солености верхнего слоя:

к 2

н .

(7.36)

S a = — — ^

в ( H - h ) 2

 

Составим уравнения равновесия сил для верхнего слоя потока

^ - Р т - Р Ув= 0

(7.37)

и для поверхности раздела потоков (эта поверхность проходит че­ рез точки на вертикали, где скорости равны нулю)

Р,

- Р у

= 0 ,

(7.38)

гг=Л

Гг=А

 

4

'

I

 

 

 

225

а также уравнение диффузии, описывающее поступление солей из нижнего слоя в верхний

= 0 ,

(7.39)

z = h

где Ад - коэффициент турбулентной диффузии солей.

Вуравнении (7.37) слагаемые определим по формулам:

-сила тяжести

PK = y J { H - h ) ,

(7.40)

- сила трения между нижним и верхним слоями

Р, = А,

dv

(7.41)

'dz z = h

-сила давления, обусловленная изменением солености во вдоль оси х, совпадающей с направлением течения реки

(7.42)

Рт = Ь Н ~ Н ) Л ’

где ув —удельный вес воды верхнего слоя; I - уклон реки; Аг -

коэффициент турбулентного обмена; 1 - носит линейную размер­ ность.

После подстановки выражений (7.40) - (7.42) в (7.37)

полу­

чим:

 

 

y sI ( H - h ) - A z ^ -

бУв ( Я -A ) 1 = 0 .

(7.43)

dz z = h

dx

 

Раскрывая слагаемые в (7.38) по аналогии с (7.40) и (7.42), по­

лучим

 

 

дУг.Р

1 = 0,

(7.44)

У г / ~

dx

 

 

где угр - удельный вес воды на границе между слоями.

226

дУв

^Г.р

Выразим в производных —А и

------ удельный вес через

йх

ас

соленость воды S по формуле (при этом влиянием температуры воды можно пренебречь)

y = y 0 + k s S ,

(7.45)

где Уо _ удельный вес пресной воды; ks = 0,8. Одновременно про­ изводную от скорости в (7.43) заменим конечными разностями,

считая, что |vH| = |vB03B| . Тогда уравнения (7.43) и (7.44) примут вид

(здесь опустим единицу в последних слагаемых с целью простоты записи этих уравнений):

dS.

y J

+ 2 A I ------

ы----------

Ь — ®- = 0 ,

(7.46)

U

z h ( H - h )

s дх

 

 

 

 

dST„

(7.47)

 

 

 

= о .

z=h

Если выполним интерполяцию между средними значениями солености воды S H и S B, учитывая, что ее распределение по глу­ бине принято по прямой, то найдем значение солености на границе между слоями по формуле

s r.p = | s . + 1-

s„.

(7.48)

 

н

 

Подставим это выражение в (7.47), тогда получим:

Угр/ Н

*- +

дх

 

(7.49)

{ H - h f

• = 0.

h - ( S H- S B)

2hHS„

дх

227

Анализ показал, что в уравнениях (7.46) и (7.49) в первом

приближении можно принять у в = у

= у . Тогда, исключая в этих

уравнениях уI , получим:

 

 

 

 

[ i h H S ^ i S ^ - S J h 2] ^ -

 

 

 

 

 

дх

 

 

 

 

- [2hHSu + (SH- S B) ( H - h ) 2

l d x

- - A*V-aH

.

1

н вА

J

k s ( H _ h y

 

С целью определения градиентов солености в (7.50) обра­ тимся к уравнению диффузии (7.39). Выразим в нем расход воды в нижнем слое через скорость и выполним дифференцирование:

hSH

SvH

dh

dSH

dS

(7.51)

+v,A .

 

+ vnh - j l = ^

 

ox

ox

dx

oz z=h

 

Чтобы определить в (7.51) производную — - , воспользуемся dx

зависимостью (7.33), тогда получим:

(7.52)

 

 

h S

-

S

 

 

 

 

 

s ^

_

s zs „

dvn _

q p

S E dh

q p

"

&

Ub &

Эх

h 2

S H- S B dx

h

 

 

(7.53)

 

{SH- S 3f

Подставим (7.53) в (7.51), одновременно производную

a s

запишем в конечных разностях, получим:

dx z= h

dSH

c 2 dSB _

2 Aj S

H~ S j,3

(7.54)

S l

=

q p

.

ox

ox

Н

 

Итак, получили систему из двух линейных уравнений (7.50)

и (7.54). Решением ее будут следующие выражения:

228

 

 

- S .) j[2AKS„ + (s„ - S , \ H - h f ] -

3

8SK ^

q pH

2 ~ h

B/L

J

k s ( H - h )

dx

h [2 H S H+ {Sa - S B)h]S2 -

\2 h H S H+ { S H- S B\ H - h ) 2] S 2B

(7.55)

2A^ { S n - S j h [ 2 H S H+ {S«

dSB _ qvH ________________________________ ks ( H - h )

&h [2 H S H + (SK - S B)h]S2 - [ 2 hHSH + (SH- S B\ H - h ) 2] S 2 '

(7.56)

Таким образом, имея зависимость (7.35), (7.36), (7.55) и

(7.56), можем рассчитать характеристики потока на устьевом уча­ стке реки и, в частности, дальность проникновения клйна соленых

вод. При этом должны быть заданы удельный расход реки q p , глу­

бина потока Н - / ( х ) , соленость морской воды в нижнем слое и

толщина нижнего слоя потока \ в створе на выходе реки в море

(в устье реки). Если последняя величина не известна, то она может быть найдена подбором в процессе решения задачи. Что касается

коэффициента турбулентного обмена Az и турбулентной диффу­ зии Ад , то их рекомендуется определять по зависимости, предло­

женной В.М . Маккавеевым для случая разноплотностных потоков:

 

 

 

 

А

Ъ Т Я .

 

А

=^ _ А

(7.57)

 

 

 

 

 

М С

 

 

д

уср

 

где

у ср

-

удельный

вес воды,

средний для обоих слоев;

У ,к 1 = кн! + 1?возв1 + -k l ~

суммарный

элементарный

расход воды

обоих слоев; М =

1440 -

коэффициент, получен обратным пересче­

том

по

данным

р. М иссисипи

для

разноплотностного потока;

!

J

1/6

коэффициент Шези;

п -

коэффициент шероховатости

С = —Н

 

-

 

п

 

 

 

 

 

 

 

 

 

русла; g -

ускорение свободного падения.

 

 

В качестве примера на рис.

7.14 приведены результаты рас-

j чета длины клина соленых вод в реке Оби при удельном расходе

I

229

воды реки q p = 0 ,5 м2/с (другие данные: глубина реки Н = 12 м,

коэффициент шероховатости дна « = 0,015, соленость воды Кар­ ского моря S = 0,03 т/м3).

JlM

Рис. 7.14. Кривые характеристик потока в зоне клина соленых вод р. Оби.

В результате этого расчета получена толщина клина в на­ чальном створе h0 = 5,01 м, а его длина L = 170 км.

На этом же рисунке приведена кривая изменения уклона вод­ ной поверхности I на рассматриваемом участке. Анализируя эту кривую и кривую изменения глубины клина h = / ( х ) , приходим к выводу, что в зоне распространения клина соленых вод имеет ме­ сто неравномерное движение жидкости и возникает кривая спада водной поверхности по типу перелива через широкий порог. Вы­ полненный анализ показал, что принятие глубины потока Н = const

на точности расчетов не сказывается.

На рис. 7.15 приведена кривая изменения длины клина соле­ ных вод в зависимости от удельного расхода воды. Для сравнения на этом же рисунке приведена и кривая, построенная для тех же условий по эмпирической зависимости Кейлегана:

230