Файл: Чайлдс Э. Физические основы гидрологии почв.pdf

ВУЗ: Не указан

Категория: Не указан

Дисциплина: Не указана

Добавлен: 10.04.2024

Просмотров: 240

Скачиваний: 5

ВНИМАНИЕ! Если данный файл нарушает Ваши авторские права, то обязательно сообщите нам.

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

dQ/dt = K 0+ j

 

 

 

что можно записать также в форме

 

 

 

 

- JL

f

31.4)

d Q l d t = 4иr A t

2 + £ .

I

(13.12)

Интегрируя, получаем

 

 

 

_і_

 

f

(Д 31.5)

Q = At 2

+B t,

1

(13.11)

где

 

 

 

 

С

 

 

 

A = C Xdc,

 

 

Со

С

 

 

 

 

 

В ^ К о +

I \idc.

 

 

со

 

 

Постоянные А и В являются характеристиками почвы, зависящими от пре­ дельных влажностей сд и С.

Дополнение 32. Вычисление скорости впитывания в слоистый профиль.

Рассмотрим профиль, который состоит из верхнего слоя толщиной L с влаго­ проводностью К jj ж нижнего слоя с влагопроводностью KL. В соответствии

с гипотезой Грина и Эмта гидростатический напор на фронте увлажнения ра­ вен fHjj, когда фронт находится в пределах верхнего слоя, и fHL, когда он

находится в нижнем слое.

Тогда, как показывают уравнения (13.2)—(13.3),

dQ\dt = А ц -{-BjjjiljL)

,

32.1)

для 0 < (1/L) < 1,

I

(13.14а)

где

Аи = Ки>

В и — К ц ( Но — f H j j ^ l L .

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

Когда фронт находится в нижнем слое, длина эффективно проводящей верхней колонки почвы равна L, а длина нижней колонки равна I L, где I

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

dQldt = K u (H0- H J + L ) / L ,

dQJdt = К L ( H j - f H L + l - L ) ! ( l - L ) ,

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

E 9- f H L + l = {àQiat)[LIKv + { l - L ) l K L \,


A U {B L I A L + 4 L )

I

32.2)

dQ/dt ~ 1 + { А и 1Аь ) ( 1 1 Ь - \ ) ’

1

(13.146)

где

A L = K l .

В L — K L (Но- f H L)lL;

A L и В J можно найти из опыта по впитыванию в материал нижнего слоя.

Поскольку уравнение (Д32.1) по существу то же, что и уравнение (13.2)» и отличается оно только значением постоянной В, его можно интегрировать

точно так же (Дополнение 30), после чего получим уравнение, аналогичное урав­ нению (13.5), а именно

t = (L /A u ) A f u [ l / L - ( B u / A u ) l n { l + ( A u / B u) l / L } ] , {

где Д/ ц — прирост влажности в верхнем слое при прохождении фронта увлаж­

нения.

Так же можно интегрировать и уравнение (Д32.2). Запишем его в форме

Д / ь ^ і + ^ Ѵ ^ / ь - і ) ]

AU (BL/AL +1/L)

32.4)

 

где

AfL — прирост влажности при прохождении фронта увлажнения в ниж­

нем

слое.

Преобразуем теперь переменную I по уравнению

 

 

B J A L + 1 /L = h

 

32.5)

так что

 

 

32.6)

dtjL = dX.

 

 

Используя эти подстановки в уравнении (Д32.4), получим

 

 

dt = (L ДfL dXJAj^j [(AL I A JJ

1

BL / A L )/K-\- 1 ].

 

 

Интегрируя от L до I, что соответствует моментам tL и t, найдем

 

 

t - t L = (LIAl ) д/ г [(Аь / А п -

1 -

B J A l ) ІпХ + К].

 

 

Переходя обратно от X к I по соотношению (Д32.5), получим

 

 

г

 

llL~\-BLj A T

 

■]

t - t L = (L /A l ) A}L { ( A J A t j - 1 - B J A l ) ln

+ ll L - 1J .

 

 

j

32.7)

 

 

1

(13.165)


Движение грунтовых вод

14Л. Природа течения грунтовых вод

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

1) изменения профиля влажности являются прежде всего след­ ствием изменений скорости выпадения осадков или скорости испа­ рения с поверхности;

2) величина подобных изменений сравнительно одинакова на больших площадях;

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

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

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

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

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

Уравнение неразрывности, относящееся к рассматриваемым ус­ ловиям, в своем наиболее общем виде представляет собой уравнение Лапласа, т. е. уравнение (11.22):

д*ф/дх* + д*Ф/ду* + д2Ф/дг2 = 0.

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

14.2. Условия, в которых осуществляется движение грунтовых вод

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

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

Ф = z ^\- H — Z — d -r d = Z.

269



Движение грунтовых вод

14Л. Природа течения грунтовых вод

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

1) изменения профиля влажности являются прежде всего след­ ствием изменений скорости выпадения осадков или скорости испа­ рения с поверхности;

2) величина подобных изменений сравнительно одинакова на больших площадях;

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

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

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

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

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

Уравнение неразрывности, относящееся к рассматриваемым ус­ ловиям, в своем наиболее общем виде представляет собой уравнение Лапласа, т. е. уравнение (11.22):

д2Ф/дх* + д^Ф/ду* -}-д*ФJdsfl = 0.

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

14.2. Условия, в которых осуществляется движение грунтовых вод

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

Когда вода поступает в почву или уходит из нее через ложе заполненного водой канала, поверхность, отделяющая почву от воды в канале, является эквипотенциалью. Например, если поверхность свободной воды находится на высоте Z, то высота любой другой точки в водозаполненном канале на глубине d под поверхностью есть Z — d, а гидростатический напор Н на этой глубине равен d. Отсюда

Ф = г + Н = Z — d + d — Z.