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

книги из ГПНТБ / Чайлдс Э. Физические основы гидрологии почв

.pdf
Скачиваний:
12
Добавлен:
27.10.2023
Размер:
21.58 Mб
Скачать

Через бесконечно большое время экспоненциальный сомножитель стано­ вится равным нулю, так что Z сводится к постоянному Z u , определяемому фор­

мулой

Zu / L = A + quB J K . {

Тем самым подтверждается уже известный факт, что Z u есть окончательная

стационарная высота уровня грунтовых вод при стационарной интенсивности осадков qu и, следовательно, определяется уравнением (18.7).

Поэтому из уравнения (Д42.6)

Z

Z U

 

- К t ! Y BL

(Д42.8)

L

L ~

L

 

 

 

Высота уровня грунтовых вод Z в начальный момент, когда величине t

произвольно приписывается нулевое значение, равна Z0, так что приведенное

выше уравнение принимает вид

 

 

 

 

 

 

 

Z Q — Z u = ос,

 

 

 

(Д42.9)

Подставив в уравнение (Д42.8) значение а из уравнения (Д42.9), получим

окончательную форму решения:

 

 

 

 

 

Z - Z . - f t - Z . ) . - * " ' ™ .

 

{

M g .io)

Дополнение 43. Колеблющийся уровень грунтовых вод.

 

 

 

Дифференциальное уравнение движения уровня (18.16) имеет вид

 

„ d Z , К

f Z

. \ —

. .

п

/

(Д43.1)

у_йГ+ _ж

( т _ Л ) _ 9 “ 9о8ш(юг)_0-

I

(18.16)

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

грунтовых вод Z и амплитуду колебаний Z 0.

Такое движение описывается сле­

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

 

 

Z = Z-\-ZQ sin (coi +Ѳ).

{

В соответствии с ним

 

 

d Z / d t =(oZ0 cos (ші + Ѳ).

(Д43.3)

Подставив эти значения Z и d Z / d t в уравнение (Д43.1), получим

 

ycüZ0 (cos соt cos Ѳ— sin шг sin Ѳ) +

 

(sin coi cos Ѳ+ cos a t sin

— g—g0 sin cqf = 0.

(Д43.4)

Уравнение это должно оставаться справедливым для всех моментов времени, поэтому постоянные члены и коэффициенты при sin coi и cos ші должны быть по отдельности равны нулю. Благодаря этому условию получаем три уравнения:

 

К 7

(Д43.5)

Y(ÙZ Q C O S Ѳ+ -р2Д sin Ѳ= 0»

—YCOZQsin Ѳ

K Z 0

(Д43.6)

cos Ѳ— g0= 0,

B L

K_

(Д43.7)

В

Из уравнения (Д43.5) сразу же имеем

tg в = — Y a > B L /K .

 

Из уравнения (Д43.8) следует, что

 

sin Ѳ=

Y a B L

i_ •

 

 

(К 2 4- У2©2В2£ 2) 2

cos Ѳ=

-----------------------

.

f (Д43.8) t (18.20)

( ^ 2 + У 2 Ш 2 Д 2 Ь 2 ) 2

Подставив эти значения sin Ѳи cos Ѳв уравнение (Д43.6), получим выра­ жение для Zÿ.

Z0

B q o / K

 

(Д43.9)

L

J. '

 

 

 

 

(1 + У 2 со2Я2£,2/£2)2

 

 

Из уравнения (Д43.7) непосредственно следует:

 

 

 

Z / L = A + B q / K . .

[

(Д43.10)

 

t

(18-18)

 

 

ГЛАВА 19

Методы измерения коэффициента влагопроводности

19.1. Общий обзор проблемы

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

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

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

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

точно характеризует ту величину, которая относится к полевым усло­ виям. Это нетрудно показать, рассмотрев какую-либо гипотетиче­ скую структуру. Так, на рис. 19.1 изображена толща с регулярным расположением вертикальных и горизонтальных трещин, расстоя­ ния между которыми равны соответственно W и L. Горизонтальные трещины непрерывны, а вертикальные смещены на отрезок WJ2 по отношению друг к другу, так что через слой находятся одна над другой. Теперь понятно, что образец шириной менее W]2 и длиной превосходящий L никак не сможет обеспечить транзитную вертикаль­ ную влагопроводность, и измеренная величина последней будет устойчиво равна нулю, тогда как в поле или в достаточно большом образце она будет существенно отличаться от нуля.

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

тока.

1

1

1

1

1

1

1

1

1___ 1

1

J

 

•л

г —ч

1

1

1

1

1

!

п п

1____1

L . . J

L

 

. w

1

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

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

19.2. Лабораторные пермиметры1

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

1 Пермиметр — фильтрационный прибор. — П р и м , перев.

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

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

 

 

фильтр помещают и на верхнюю

 

 

поверхность. Таким образом, вла­

 

 

гопроводность колонки фактически

 

 

не постоянна

по высоте, и общая

 

 

разность потенциалов, измеряемая

 

 

по разности уровней воды в верх­

 

 

нем и нижнем резервуарах, харак­

 

 

теризует

лишь средний

градиент

 

 

потенциала в

составной

колонке,

 

 

а не

истинный градиент

потен­

 

 

циала

в

образце. Поэтому

в ко­

 

 

лонку вблизи ее концов вводят

 

 

манометры. Разность уровней воды

 

 

в их открытых коленах есть мера

 

 

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

 

 

между точками, в которых они

 

 

установлены. Так

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

 

 

тинный

градиент

потенциала.

 

 

Поток измеряют, наблюдая за

 

 

поступлением воды в нижний при­

 

 

емник, а коэффициент влагопровод­

 

 

ности (в данном случае коэффи­

Рис. 19.2.

Пермиметр с постоянным

циент фильтрации)

вычисляют не­

напором.

Напор в данном случае

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

по

закону

Дарси

поддерживается с помощью мариот­

* = <гг)(ж ).

това устройства.

 

где QJt есть измеренная величина потока, А — площадь поперечного сечения образца в колонке, АФ — разность потенциалов, измеря­ емая манометрами, D — расстояние между манометрами.

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

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

ʜ 2W 2SZ v5n =
Пермиметр с падающим напором.

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

ции, так и

градиент

потенциала

в подобном

приборе

непрерывно

меняются.

Поэтому

необходимо

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

і=дф

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

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

величине потока, а высота ме­ Рис. 19.3. ниска в тот же момент представ­ ляет разность потенциалов.

Зависимость между общим количеством профильтровавшейся воды и временем выводится следующим образом. Пусть Q есть рас­ ход в данный момент. Если а — площадь сечения верхней трубки, а I — высота мениска над постоянным уровнем воды в нижнем мано­

метре, то

(19.1)

Q-- —а (dl/dt).

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

Q ^ K A l / D ,

(19.2)

где А — площадь поперечного сечения колонки. Следовательно, из уравнений (19.1) и (19.2) имеем

—a (dljdt) = KAljD.

(19.3)

Если в моменты £0 и £ наблюдаемые уровни мениска равны соот­ ветственно 10 и I, то, интегрируя уравнение (19.3) от t 0 до t, полу­ чаем

1Q (WO — (А/а) (K/D) (t —10).

(19.4)

19.3. КачающиЗся пермиметр

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

2 Такой раствор приходит

Рис. 19.4.

Качающийся пермиметр.

в равновесие

с

почвой,

1 — нейтральное

положение, 2 — уровень отсчета.

ПО ЧТИ НѲ М ѲН ЯЯ

Ѳв

И ОННЫ Й

состав, и остается в равновесии столько, сколько необходимо для определения влагопроводности. Пермиметр, удовлетворяющий этим требованиям, был разработан Чайлдсом и Пуловассилисом [40]. Одна из его конструкций показана на рис. 19.4. Исследуемый образец помещают в U-образную трубку, так что обе поверхности его остаются свободными. Фильтров или мембран не требуется. Насыщающий раствор образует небольшие столбики жидкости над поверхностями образца. Благодаря этим столбикам существует составляющая потенциала, связанная с гидро­ статическим напором. Свободная жидкость стоит в капиллярах меньшего диаметра, чем трубка с образцом, что повышает точность измерения, как в пермиметре с падающим напором. Прибор устана­ вливается на раме, которая с помощью мотора и эксцентрика может качаться вокруг горизонтальной оси, участвуя в простом гармониче­ ском колебании с любым желаемым периодом Т.

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

О, соответствующую равновесному положению мениска. Эта отметка

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

zR = z-\-z0sin(2nt/T).

(19.5)

Рис. 19.5. Типичные кривые, получаемые при работе с кача­ ющимся пермиметром.

Через достаточный промежуток времени в пористой колонке установится движение жидкости, которое по частоте будет соответ­ ствовать колебательным движениям с периодом Т ; мениск в каждом колене будет подниматься и опускаться относительно отметки О

сэтой же частотой. Однако такое согласованное движение мениска

вобщем случае не совпадает по фазе с возбуждающей частотой; амплитуда его также неизвестна.

Если высота мениска в правом колене относительно О равна ZR, то, как показано в Дополнении 44, колебание жидкости описывается уравнением

ZB= Z osin(2jt*/7, + j t - 0 ) ,

(19.6)

где Z0 — амплитуда колебаний, (я — Ѳ) — разность фаз по отноше­ нию к возбуждающему колебанию. Ѳ определяется формулой

tg Ѳ=

anl

(19.7)

 

T A T

 

где I и А — соответственно длина и площадь сечения образца пори­ стого материала в колонке, а — площадь поперечного сечения капилляров, в которых движутся мениски. Амплитуда Z 0 опреде­ ляется соотношением

------ ZS~- T ‘

(19.8)

(1 + tg2 Ѳ) 2

 

На рис. 19.5 изображены результаты опыта

в форме графика,

на который нанесены колебания колена прибора zR и колебания ме­ ниска в колене ZR. На этом графике можно непосредственно измерить разность фаз (л — Ѳ) и таким образом найти Ѳ. Можно также изме­ рить амплитуду Z 0 и возбуждающую амплитуду z0, найдя затем Ѳ косвенно, по уравнению (19.8). В любом случае из уравнения (19.7) можно вычислить К, поскольку все остальные входящие в него вели­ чины известны из опыта.

Точность измерений максимальна, когда Ѳблизко к я/4 или к 45°, что можно обеспечить, подбирая режим и размеры установки. При таком Ѳ

tg0 = —(anl)/(KAT) = —1.

Период Т и отношение плошадей сечения а]А произвольны, по­ этому, сделав грубые предварительные измерения К, можно подо­ брать параметры, необходимые для точных измерений.

19.4. Измерение коэффициента влагопроводности ненасыщенных почв

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

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

основанные на этой закономерности (см. Г л о б у с А.

М. Экспериментальная

гидрофизика почв. Гидрометеоиздат, 1969, Л.) — П р и м ,

перев.

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

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

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

Равномерно увлажненные колонки с образцом можно получить посредством инфильтрации. В параграфе 12.2 было показано, что если происходит инфильтрация с постоянной скоростью в достаточно длинную колонку, заканчивающуюся уровнем грунтовых вод, то. верхняя часть колонки имеет постоянную влажность, и влагопровод-, ность этой части равна скорости инфильтрации. Таким образом, для нахождения влагопроводности достаточно только поддерживать постоянной и измерять скорость инфильтрации, а также измерять влажность почвы после установления стационарного состояния. Применяя этот метод, Чайлдс и Коллис-Джордж [38] устанавливали колонку под разными углами к вертикали, за счет чего можно было получать различные градиенты потенциала. Таким способом им

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