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

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

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

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

Добавлен: 10.04.2024

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

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

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

В обеих формулах I есть толщина равновесного слоя воды в сква­ жине, г — радиус последней, t — время, за которое уровень воды в скважине после откачки поднимется с отметки Z Y на отметку Z2. Коэффициент А определяют эмпирически. Судя по полученным дан­ ным, он в разных случаях несколько различен. Для однородной толщи его предельное значение считают равным rZ/0,19.

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

Ван Бэвел и Керкем [6] получили строгое решение уравнений течения воды в скважину, пронизывающую весь водоносный пласт до водоупора. Это решение имеет вид

dZ/dt = 1,6 (IK/r) 2,

где

СО

2 = 2 (—1)(л_1,/2 cos (nnZ/21)

n = l

п = 1, 3, 5, . . .;

(19.39)

Кі (nur/21) К 0 (nnr [21) ’

здесь К 0 и К і — функции Бесселя второго рода нулевого и первого порядков соответственно, имеющиеся в математических таблицах. Таким образом, наблюдения за ходом уровня воды в скважине поз­ воляют определить одновременно уровень и скорость подъема в любой заданный момент и затем вычислить К.

Для случая, когда скважина не достигает водоупора, Керкем [94] определил коэффициент 2, входящий в уравнение (19.39), теорети­ чески, в довольно сложной форме, доведенной до числа для некоторых определенных геометрий. Кроме того, Ван Бевел и Керкем [6] с по­ мощью электроаналогового метода получили кривые, связывающие 2 с различными отношениями между глубиной скважины и глубиной водоупора.

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


горизонтальной влагопроводности методом двух скважин, после чего из пьезометрических измерений можно вывести истинную вер­ тикальную влагопроводность. Лабораторная проверка метода была проведена Чайлдсом, Коулом и Эдвардсом [37], а полевая мето­ дика — Чайлдсом, Коллис-Джорджем и Холмсом [39].

В методе двух скважин подходящим буром проходят пару сква­ жин диаметром г и глубиной I каждая, на расстоянии 2d друг от друга; после этого ожидают, пока вода в скважинах не установится

на равновесных уровнях. Затем воду с постоянной скоростью Q перекачивают из одной скважины в другую, так что через некоторое время устанавливается стационарное состояние, при котором уро­ вень воды в одной скважине возрастает на величину, равную пониже­ нию уровня в другой. Разность этих уровней равна AZ. Если уровень грунтовых вод не слишком нарушен откачкой, течение воды между* скважинами происходит приблизительно по горизонтали, не считая окрестностей дна скважин. Концевой эффект можно учесть с помощью электроили гидравлических аналогий или же исключить, проводя опыты с различным заглублением скважин. Мы предположим, что

на отрезке I поток Q горизонтален.

Горизонтальный поток представляет собой хорошо известную двухмерную задачу для течения между точечным источником и точеч­ ным стоком. Поверхности скважин соответствуют паре круглых эквипотенциалей. Решение приводится в обычных руководствах (см., например, [146]) и имеет вид

О

п К „ 1 A Z

(19.40)

= ______

Ѵ

ch-l (d/r)

 

Поскольку все величины, входящие в это уравнение, можно изме­ рить, нетрудно вычислить Кн .

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

равной (КцКу) 2 , если бы все линейные размеры системы изменились j_

в отношении (Кн/Кѵ) 2 . Поскольку полость имеет нулевую длину, этот размер остается в преобразованной системе неизменным, изме­ нение же глубины в общем несильно влияет на коэффициент А в урав­ нении (19.32), когда его используют применительно к данному слу­

чаю.

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


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

л_

 

К = (КНКѴ)К

(19.41)

Найдя Кн из опыта с двумя скважинами, можно сразу же вычис-

лить Ку. Если при этом окажется, что полученное отношение (Кн ]КѴ) 2 настолько влияет на глубину полости в преобразованном простран­ стве, что это приводит к заметному изменению коэффициента А, это новое значение А подставляют в уравнение (19.32) и вычисляют следующее значение К ѵ, используя прежние опытные данные. Так, методом последовательных приближений получают наконец вели­ чину Ку, удовлетворяющую уравнение (19.32), в котором коэффи-

циент А соответствует величине (Кн Ку) 2 .

Обзор методов измерения коэффициента фильтрации в полевых условиях выполнен Доннаном [57]. Все методы характеризуются как сложные и подверженные значительным погрешностям. Некоторые недостатки связаны с реальными трудностями, коренящимися в при­ роде почв, например с нестабильностью, из-за которой невозможно Сохранять полость или скважину постоянных размеров. Другие воз­ ражения, однако, являются результатом недостаточного осознания статистической природы влагопроводности. Например, там, где влаго­ проводность определяется сложением почвы, а расстояние между границами структурных отдельностей имеет порядок десятка санти­ метров, размер образца или исследуемой зоны должен быть значи­ телен. Поэтому применение скважин или полостей диаметром 2— 3 см, что нередко практикуется, влечет за собой ошибки. Влияние вариабельности коэффициента фильтрации на вариабельность и асимметрию уровней воды при использовании метода двух скважин рассмотрено Чайлдсом, Коллис-Джорджем и Холмсом [39].

Д О П О Л Н Е Н И Я

Дополнение 44. Колебания жидкости в качающемся пермиметре.

Пусть в приборе, показанном на рис. 19.4, мениск в правом колене нахо­ дится на высоте Z R по отношению к фиксированной отметке О на том же колене.

Тогда, если а есть площадь поперечного сечения трубки, в которой движется

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

dQ/dt = adZR/dt.

(Д 44.1)

Высота отметки О относительно некоторого фиксированного нулевого

уровня, согласно уравнению (19.5), равна

ZH= Z + Z0 sin (2Ш/Т).

I ^ (19І5)

Благодаря симметричности движения относительно оси качаний высота соответствующей отметки на левом колене равна

zL = г —z0 sin (2ntjZ).

(Д 44.3)


Когда мениск в правом колене находится на высоте ZR над отметкой О, мениск в левом колене должен находиться на глубине ZR ниже отметки на этом

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

Фд = гй + Ѵ

44.4)

высота левого мениска над тем же уровнем равна

 

 

=

^4.5)

По определению потенциала [уравнение (9.3)], Фд и

есть гидравли­

ческие потенциалы на поверхностях материала соответственно в правом и левом коленах. Следовательно, потенциал в правом колене больше потенциала в левом колене на величину ЛФ, которую можно выразить, комбинируя уравнения (Д44.4)—(Д44.5) с уравнениями (Д44.2)—(Д44.3):

АФ = ФВ ~ Ф Ь = 2 {z0 sin (2ntfT)-\-ZRj.

(Д 44.6)

Если длина колонки

пористого

материала, измеренная вдоль

изгиба,

равна I, а площадь ее сечения А , то уравнение, выражающее закон Дарси, будет

иметь форму

 

 

 

 

dQJdt

ДФКА

( Д 44.7)

 

I

 

 

 

Подставив dQ/dt из уравнения (Д44.1) и ДФ из уравнения (Д44.6), получим

после небольших преобразований

 

 

dt

2КА

[z0 sin (2лt/T)-\-ZR],

( Д 44.8)

al

 

 

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

еще требуется найти. Поэтому решение уравнения (Д44.7) имеет вид

 

ZR = В s m ( 2 n t / T -\-а).

(Д 44-9)

Соответственно

 

 

 

 

 

dZр

2лВ

 

( Д 44.10)

 

- J L

^ ^

c o s (2m /f + a).

Подставив эти значения ZR и

dZR/dt в уравнение (Д44.8), получим

В

cos (2 яг/Г + а) =

2КА

[z0 sin (2ntJT)-\-B sin (2nt/T -fa)].

Т

al

 

Разложив тригонометрические функции аргумента {2nt/T + a), получим

уравнение

—ÿ - [cos (2яЦТ) cos a —sin (2яt/T) sin aj =

КА

=------ - [ ( z 0 - f i ? cos a) sin (2nt/T)-^-B sin a cos (2яt/T)].

Это уравнение должно удовлетворяться в любой момент времени, что воз­ можно, только если коэффициенты при cos (2 п t/T) в левой части равны таковым


в правой части; то же относится к коэффициентам при sin (2nt/T). Это условие

соответствует паре уравнений:

 

л cos ос

/L 4sina-

 

 

 

al

 

„ I .

КАТ

 

^ _

K A T z о

В ^sirm а ---------г- cos а

/'

лаі

 

лаі

 

Из уравнения (Д44.11) непосредственно получаем

tg a = —

лаі К А Т »

используя это выражение в уравнении (Д44.12), найдем

В = z0 C O S а

 

 

(l + tgü a ) 2

Избавимся от минусов при вычислениях, введя

Ѳ= —a,

при этом из уравнения (Д44.13) следует:

44.11)

44.12)

(Д 44.13)

44.14)

* Û

4

лаі

/

44.15)

tg Ѳ= .

tg a =

—g ^ j T

1

(19.7)

 

 

 

 

Из уравнений (Д44.9) и (Д44.14) найдем решение:

 

 

Zo

— sin (2лtJT — Ѳ),

 

 

 

 

 

( l + t g 2 Ѳ)

 

 

 

 

или, избавляясь от минуса,

 

 

 

 

 

 

 

Ч =-

Zo

sin (2лtJT -j- я — Ѳ).

 

44.16)

 

( l + t g 2 Ѳ)

 

 

 

 

 

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

нием

 

 

Zo

 

I

44.17)

 

 

Zn = -

 

 

 

 

 

I

 

(19-8)

 

 

 

 

 

 

 

 

d + tg » ѳ)

2

 

 

 

а уравнение движения (Д44.16) примет вид

 

 

 

 

г :

ZR = Z 0 sin (2я</7 +

я —Ѳ).

f

44.18)

I

 

(19.6)

 

 

 

 

 

 

Дополнение 45. Переходный режим в откачиваемой скважине.

 

 

 

Вначале

предположим,

что скважина, откачиваемая с постоянной ско­

ростью Q, полностью пронизывает водоносный пласт толщиной I, ограниченный

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

V= —К (ІФ/dr.