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

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

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

 

 

Сложность зависимостей

 

 

(10.23) и (10.24) ограничивает

 

 

возможность

интегрирования

 

 

уравнения (10.21) в аналитиче­

 

 

ском виде. Поэтому часто ис­

 

 

пользуют

численные

методы

 

 

решения этого уравнения: ме­

 

 

тод конечных разностей и ме­

 

 

тод конечных элементов. Од­

 

 

нако и в этом случае встреча­

 

ет

ются трудности, обусловлен­

 

 

ные большим объемом подго-

Рис. 10.2. Схема гистерезиса

товительных работ и вычисле-

зависимости потенциала влаги

ш й на ЭВМ;

особенно при пе-

от объемной влажности.

ременных

начальных

и гра-

, ,

,

1 - обезвоживание, 2 - увлажнение.

с

 

 

г

 

 

ничных условиях и рассмотре­

 

 

нии двух-

и трехмерных про­

цессов переноса влаги.

Поэтому поиски приближенных аналитических выражений, отражающих основные закономерности движения почвенной вла­ ги, весьма актуальны. Для ряда характерных случаев движения влаги в почве такие выражения получены [53,56].

При решении уравнения (10.21) чаще всего назначаются

следующие начальные и граничные условия:

 

1)

начальные условия:

 

 

 

при т = 0

Wz = f x {z)\

(10.25)

2)

граничные условия первого рода (влажность):

 

а) на поверхности почвы (z = 0)

 

 

п р и т > 0

Wa = const

(10.26)

или

 

 

 

 

Wu = f 2{т),

(10.27)

б) на глубине почвы (z =

Г)

 

 

п р и т > 0

й^/=/з('с).

(10.28)

361

Глубина в условии (10.28) определяется точкой, начиная

скоторой значение влажности с ростом z не меняется;

3)граничные условия второго рода (поток влаги):

а) на поверхности почвы (z = 0) скорость переноса влаги v

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

 

 

при т > 0

vn = const

(10.29)

или

 

 

vn = - D dW /dz - kz ,

(10.30)

б) на глубине почвы (z = I)

 

п р и т > 0

v , = const.

(10.31)

Влагоперенос в почве при установившемся режиме. Рас­ смотрим случай, когда увлажнение почвы происходит с постоян­ ной скоростью, что соответствует граничным условиям (10.29) и (10.31). В этом случае будет наблюдаться перенос влаги по глуби­ не, но при установившемся режиме в каждой ееточке; количество же влаги в этих точках не меняется (dW /dx = 0).Тогда уравнение

(10.15) при нисходящем потоке влаги (вертикальная ордината от­ считывается снизу) примет вид:

± к

dqV -1 = 0 .

(10.32)

dz

v dz

 

Проинтегрировав уравнение (10.32), получим:

kz (dq>w / d z - l) = v ,

(10.33)

где v - постоянная скорость потока влаги (увлажнения).

Представим уравнение (10.33) в виде

d y w /d z = v/kz + 1.

(10.34)

Решая это уравнение относительно z и

интегрируя его от

уровня грунтовых вод, где г = 0 и ф ^ = 0 , довысоты z, где водный потенциал равен ф^, получаем:

362

ф(}/

Z= \ { d ( - y w ) / { l - 4 K ) l

(10-35)

о

Выражая коэффициент влагопроводности аналитической за­ висимостью kz = /( ф ж), получаем решение интеграла (10.35).

В случае установившегося движения воды от уровня грунто­ вых вод к поверхности, получим следующее решение уравнения

(10.15):

 

4V

 

* = - JVqv/O + v/ kz)]-

(Ю.36)

о

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

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

цию A(z, х) = 0 . При решении этого уравнения почву по глубине

принимают бесконечной или полубесконечной. Многие исследо­ ватели его решение осуществляют с применением преобразования Больцмана, заключающегося во введении новой переменной:

t! =z/ ( 2 ^ ) ,

(10.37)

где D q - коэффициент диффузии при влажности Wn . Такой прием позволяет превратить дифференциальное уравнение диффузии в частных производных в обыкновенное дифференциальное уравне­ ние.

Если решение уравнения (10.22) осуществляется для относи­ тельно малых значений влажности, то слагаемым dkz/ dz можно пренебречь. Тогда уравнение примет вид

d W

d f ^ d W r\

дх

, D -

(10.38)

dz I

dz

Решением этого уравнения при начальных и граничных ус­ ловиях (10.25) и (10.26) будет следующее выражение:

363

V = \ 4 d J z{W, - W j d w j d ^ о,

(10.39)

где

 

Wr = ( W - W n)/{W2 - W n).

(10.40)

Движение грунтовых вод в зону аэрации по капиллярам. При наличии грунтовых вод происходит подъем их по капилляр­ ным порам почвогрунта под действием капиллярных сил. Высота подъема в зависимости от почвогрунтов (диаметра пор) может достигать 4 - 5 м. Поскольку поры неодинакового диаметра, то подъем воды в них происходит на разные уровни. Скорость подъ­ ема воды в разных порах различная: в тонких капиллярах она меньше, хотя подъемная сила в них больше, чем в капиллярах большего диаметра.

Следуя работам А.П. Летунова и А.А. Роде, получим зави­ симость высоты и скорости капиллярного подъема воды от времени.

Так как рассматривается только капиллярное движение вла­ ги в почве, то уравнение (10.14) с учетом уравнения (10.6) примет вид

^ =

(10.41)

дт dz dz

Здесь градиент капиллярного (гидравлического) потенциа­ ла может быть выражен следующим образом:

0ф* / & = (йшах - К )/К .

(10.42)

где /гтах - предельная возможная высота капиллярного подъема в

данной почве; hT - высота подъема в момент времени т.

Переходя к количеству воды, поднимающемуся по капилля­ ру через единицу поперечного сечения, уравнение (10.41) с учетом (10.42) можно переписать в следующем виде:

dh/dx = kz (/zmax - hx)/h%.

( 10.43)

Интегрируя это выражение, получаем:

3 6 4

* = 7 - t a » '’’fc a .A 'W - * ,) ] - * ,) ■

(10.44)

к

В начале капиллярного подъема воды, когда /гт мало по сравнению с /гтах, уравнение (10.43) можно записать в виде:

dh/dx

= kz hmax/h x .

(10.45)

Выполнив интегрирование уравнения (10.45), получим:

 

h2

1

 

Т =

-5Г “

Г

z

(10'46)

 

max

 

 

или, решая относительно высоты подъема,

 

К =

 

 

(10.47)

Введем обозначение /г0

= yj2kzhmsK , тогда

выражение

(10.47) примет вид

 

 

 

 

hx

= h 0x ° \

(10.48)

а в более общем виде

 

 

 

 

К = \ х \

 

(10.49)

или

 

 

 

 

т =

(йт/А о Г .

(10-5°)

где (3 - показатель степени, зависит от вида почвогрунтов.

Дифференцирование выражения (10.48) приводит к формуле

vh = dh /dx = fih0^ - 1.

(10.51)

Решая совместно уравнения (10.50) и (10.51), получаем:

vA= p ^ /p^ ’“1/p).

(10.52)

365

Подставив в выражение (10.52) среднее значение показателя Э = 0,3, соответствующее, например, суглинистым почвам, полу­ чим:

v, * 0,3 Л3’33/ / ^ 33 = A /h 2* ,

(10.53)

где А = 0,3/?о'33 - постоянная.

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

10.4. Мерзлотное пучение некоторых почвогрунтов

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

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

Такие линзы прослеживаются, например, повсеместно вдоль берегов северных рек, в обнажениях которых они достигают раз­ меров порядка 20 - 30 м. Образование линз льда и связанное с ним пучение грунта необходимо учитывать, поскольку мерзлотное пу­ чение почво-грунта приносит существенный ущерб хозяйственной деятельности человека. Например, образование линзы льда в грун­ товой дороге приводит к ее вспучиванию. Весной же при оттаива­ нии, грунт дороги сползает на стороны, а после расстаивания са­ мой линзы образуется провалина. Таким образом, уменьшается прочность и надежность дороги. Давление при морозном пучении грунта может достигать порядка 20 МПа.

366

В рассматриваемом частном примере методом борьбы с мерзлотным пучением грунта является пропитывание его хими­ ческими растворами в целях цементации (упрочения) или засыпка под полотно дороги щебенки или песка на глубину промерзания грунта. В этих материалах вода не задерживается.

Чем же объясняется образование и рост с нижней поверхно­ сти линзы чистого льда в промерзающем грунте.

В настоящее время трудами ученых и, в частности Б.Ф. Рельтова, установлено, что причиной перемещения воды из области с малым давлением в область с более высоким давлением (снизу вверх), определяемым по формуле

 

Р = РмА .г + Р А >

 

 

(10.54)

где рмг и рл

- плотность мерзлого грунта и льда;

hWT и /гл - тол­

 

 

щина слоя мерзлого грунта и льда, яв­

 

 

ляется химический и термический ос­

 

 

мос.

 

 

 

 

5

 

Химический осмос

-

это про-

5

h

никновение молекул

растворителя

 

 

через

полупроницаемую

 

мембрану,

 

 

разделяющую чистый растворитель и

 

1

раствор или два раствора различной

 

2

концентрации.

 

 

 

 

 

 

Развивающееся

в .

результате

 

 

осмоса избыточное гидростатическое

 

 

давление измеряется столбом раство­

 

 

ра высотой h (рис. 10.3),

при котором

Рис. 10.3. Схема осмометра.

устанавливается осмотическое равно­

1 - растворитель,2-раствор,

весие. Это давление получило назва-

ние о см о т и ч е с к о го . 0 Н0

МОЖвТ быть

3 - п о л у п р о н и ц а е м а я м ем б р ан а .

вычислено для разбавленных раство­

 

 

ров неэлектролитов по закону Вант-Гоффа:

 

 

 

 

л

= C R T ,

 

 

(10.55)

где п —осмотическое давление;

С = n/V — молярная концентра­

ция раствора; п - число молей неэлектролита; V -

объем раствора;

367

R коэффициент пропорциональности, численно равный универ­ сальной газовой постоянной; Г - абсолютная температура.

Согласно этому закону, одномолярный раствор неэлектроли­ та при О °С обладает осмотическим давлением 22,6 • 105 Па, а де-

цимолярный - 2,26 • 105 Па. Таким образом, если рассматривать почвогрунт как полупроницаемую мембрану с большой концен­ трацией солей в растворе у его поверхности по сравнению с ниже­ расположенными слоями (т.е. имеет место градиент концентрации раствора), то можно обнаружить подъем воды к поверхности и об­ разование некоторого давления.

Рассмотрим пример.

Пусть почвогрунт, содержащий раствор концентрацией ст0,

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

Итак, при замерзании поч-

 

вогрунта с раствором соль будет

 

выпадать в осадок, что приведет к

 

увеличению концентрации раство­

 

ра на границе мерзлого и талого

 

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

 

ст0 . Ниже этой границы концен­

 

трация раствора будет изменяться

 

от максимальной сттах до естест­

 

венной по кривой 3, показанной на

Рис. 10.4. Схема распределения

рис. 10.4.

 

Наличие

градиента концен­

концентрации раствора в почво-

трунте при его промерзании.

трации д а / д г

приводит к явлению

 

химического осмоса, т. е. вода из

1- зона чистого льда в грунте (<з = 0),

нижних слоев грунта будет под­

2 - нижняя граница льда (ст= иш ),

3 - ход концентрации в талом грунте

ниматься к верхним - к поверхно­

О Т СУщах Д О С7о.

сти льда и там замерзать. Поступ­

 

ление воды вверх осуществляется со скоростью

368

vocu=~kOd°/dz,

(10.56)

где k0 — коэффициент пропорциональности, зависящий от почвог-

рунта; z - ордината, направленная вниз.

При промерзании грунта образование линз льда будет про­

исходить при выполнении следующего условия:

 

v „ _ > 5 V r/ * .

(Ю.57)

т.е. когда скорость промерзания грунта меньше скорости подтока влаги к границе (к линзе).

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

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

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

к гидрологическим наблюдениям.

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

10.5. Таяние снежного покрова

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

369

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

Основоположником учения о снежном покрове: его структу­ ре, физических свойствах и таянии, следует считать нашего сооте­ чественника А.И. Воейкова (1871). В дальнейшем этим важным во­ просом занималось очень большое число ученых. Не будем их здесь перечислять, отметим только, что существенный вклад в изучение снежного покрова и исследования процессов формирования талых вод как экспериментальный, так и теоретический внес сотрудник Государственного гидрологического института (ГГИ) П.П. Кузь­ мин. Результатом его исследований явились три монографии, по­ священные снежному покрову ( Физические свойства снежного по­ крова (1957), Формирование снежного покрова и методы определе­ ния снегозапасов (1960), Процесс таяния снежного покрова (1961)). Эти исследования явились итогом разработки методов количест­ венной оценки характеристик снеготаяния. В основе этих методов лежат уравнения теплового и водного балансов. П.П. Кузьминым изучалось снеготаяние на открытых полях, на полях, защищенных лесными полосами, и в лесу в зависимости от его загущенности.

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

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

3 7 0