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

книги из ГПНТБ / Климентов П.П. Динамика подземных вод учеб. для геологоразведоч. техникумов

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

Суммируя почленно уравнения (V, 12), найдем общее падение напора при фильтрации воды через слоистую систему АЯ:

П

АН = 2

1

І А + А

+ , . . + А )

(ѴДЗ)

г—1

о

' ki ~ кг

kn >

 

 

 

 

 

Исходя из среднего для всей слоистой толщи коэффициента фильтрации при движении воды перпендикулярно напластованию kCVtJ_ разности напоров АН, длины пути фильтрации I и площади потока о, его расход можно определить по следующей формуле:

 

АЯ

 

 

(V,14)

 

Q kcp,± - ^ ' w,

 

следовательно, падение напора

 

 

 

 

QI

 

 

(Ѵ,15)

 

АН =

 

 

 

^cp,j_û)

 

 

 

Приравняв правые

части уравнений

(V, 13 и V, 14),

имеем:

QI

h

,

К

(V.16),

kcp,±(ù

kz +

' " +

ka )■

откуда средний коэффициент фильтрации слоистой неоднородной толщи, приводимой к условно однородной при движении в ней воды перпендикулярно напластованию, определится по формуле:

&ср,_ _

 

I

hi

Нг +

... + hn

 

 

 

Л2_

(V,17)

Г

+ Г + - - - + Г 1

 

hn

ki

k2

kQ

к1

«2

кп

При фильтрации нормально к напластованию также возможно приведение слоистой толщи к условно однородной методом вирту­ ального приведения мощности. Определение значения приведенной мощности, как это следует из формулы (Ѵ,17), можно проводить по следующей формуле:

^np,_L '— ko '

 

) = * » Ü "

mi

(V,18)

ki ' k2

ki

ka

1=1

73, б, при

Так, например, для

схемы,

изображенной

на рис.

фильтрации подземных вод перпендикулярно напластованию при­ веденная мощность слоистой системы, состоящей из четырех слоев, при условии приведения всей толщи к коэффициенту фильтрации первого слоя (k0 = ki), будет следующей:

ki

ki

ki

m np,J_ = mi 4“ T ~ m 2 +

— m 3 +

-T— m t-

k2

k3

kk

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

большим коэффициентом фильтрации по сравнению с другими слоями, приведенная мощность т пРіХ будет больше суммарной мощности системы т и, наоборот, — при наименьшем значении k0 приведенная мощность будет меньше суммарной мощности слоис­ той системы.

Сопоставление выражений (Ѵ ,9 и Ѵ ,1 7 ) показывает, что величи­ на среднего коэффициента фильтрации кср при движении воды по напластованию является максимальной, а величина среднего ко­ эффициента фильтрации при движении воды нормально к наплас­ тованию /гСр,х— минимальной.

Разберем для примера простейший случай залегания слоистых пластов пород, сложенных двумя слоями одинаковой мощности с коэффициентами фильтрации k x и k2. Найдем отношение а =

— kCplkcPix в зависимости от соотношения слоев с коэффициентами фильтрации ki и k2.

Для указанного сложения пластов пород, очевидно, будем иметь

 

,

__

ki + k2

 

__

 

2

 

 

 

 

 

^ср

-

&ср, д

—----------J- •

-

( Ѵ ,1 9 )

 

 

 

 

 

 

 

 

ki

k%

 

 

Отношение первого выражения ко второму дает:

 

 

ken

1 ,,

. . . /

1

1

\

I

/

ki

 

k2

2).

(V,20)

 

= 1 ( * . +

fe) ( I +

ф ) = J - ( ф - +

i î - +

«CP,_L 4

 

'

« î

k2>

4 '

k%

 

kil

1

 

ki

m найдем:

 

 

 

 

 

 

 

 

 

 

 

При — =

 

 

 

 

 

 

 

 

 

 

 

k2

 

 

 

 

 

 

 

 

 

 

 

 

 

 

<cp

 

( m + ± +

i )

-

±

<

s

± 1):

 

(V,21)

 

<-cp,J_

 

'

 

m

'

 

 

4

 

m

 

 

 

(OT+ l)2

 

>4.

Следовательно, kcv> k cv,x при

Отношение------------ всегда

 

m

 

 

 

 

 

 

 

 

 

 

 

любом значении m.

 

 

 

 

 

 

 

 

 

 

 

 

ДВИЖЕНИЕ ПОДЗЕМНЫХ ВОД В МНОГОСЛОЙНЫХ

 

 

ПЛАСТАХ

 

 

 

 

 

 

 

 

 

 

 

 

 

Количественная оценка условий движения грунтовых и

напорных потоков в слоистых

неоднородных

пластах

проводится

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

Так, в условиях равномерного движения подземных вод (при постоянстве мощности слоистого пласта по потоку) определение расхода потока проводится по общим формулам Дарси для единич­ ного расхода (Ѵ,7).

JLJ _

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

Так, например, для условий фильтрации безнапорного потока в четырехслойном пласте, ограниченном сечениями 1 и 2, в которых значение напора составляет Н х и Н2 (см. рис. 73), расход потока может быть определен по следующим формулам:

при использовании средневзвешенного значения коэффициента фильтрации kcp-

Я —«срп/—kcph —---- ,

где

 

 

 

Tl_2

 

kihi -J- k2h2-f- k3h3-f- &4/14

 

 

 

ki h2-j- h3-j- /Î4

 

при приведении всей толщи к коэффициенту фильтрации перво­

го пласта

(k0 = kl) и использовании значения

приведенной мощно­

сти /гпр:

 

 

 

 

 

 

Я =

kihjjp Я 1

- Я

2

 

 

 

Т 1 -2

 

 

где

hnp

ki ,

, k

3

 

h2

-f- — h3+

)■

 

 

kl

kl

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

Я (kihi -f- k2h2-f- k3h3-f- &4Й4)

---------- = Яі “Ь Я2 “Ь Яг Яі- (V,22)

 

T 1 _ 2

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

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

Движение грунтовых вод в слоистом пласте при расположении кривой депрессии в верхнем слое. Если пласт состоит из нескольких слоев с коэффициентами фильтрации k\, k2 k3 ... и мощностями —

m i, m2, m3,

a кривая депрессии расположена

в верхнем слое,

мощность которого изменяется от hn' в первом

сечении до hn" во

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

и которая рассматривается как

безнапорная. Кривая депрессии,

1

г

Рис. 75. Движение грунтовых вод в слоистой толще с кривой депрессии в пределах верхнего

слоя

располагающаяся в верхнем слое, является вместе с тем пьезомет­ рической кривой для нижней напорной части потока. Расход потока определяется как сумма расходов его верхней и нижней частей. Так, для условий горизонтального водоупорного ложа (г= 0) фор­ мула для определения расхода потока в слоистой толще имеет вид:

Q

-f- kzm.2 k3m3

 

m i )

Li—2

1Li-2

 

 

При

наклонном залегании водоупорного

ложа (і^О , НфК)

формула (Ѵ,23) видоизменяется:

^

_ JLj[

 

 

 

q = {kitiii -f- kztriz -|- k3m3-)- ...) -— ------- [-

 

 

 

Li-2

 

K + K ' H i - H z

 

 

+ kn-

 

( V , 2 4 )

 

2

 

 

где kn— коэффициент фильтрации верхнего слоя пласта.

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

только одного верхнего слоя. Решение задачи в таких условиях для горизонтальнослоистого пласта можно получить, используя функ­ цию Н. К. Гиринского G. Определив значение функции Гиринского для начального сечения потока Gx и для конечного— G2 , расход потока находят по формуле:

Я =

1

2

(V,25)

Значение функции Ги­ ринского G для сечения с мощностью потока h оп­ ределяется по следующей формуле:

G — kitrii (h — ei) +

+ ^2m2 {h-- Z2) + . .

. +

+ knmn(h — z » ) .

(V,26)

 

 

 

где Z\,

z 2, ...,

z n расстоя­

 

 

 

ние от середины соответ­

 

 

 

ствующего

слоя до водо­

Рис. 76. Движение грунтовых

вод

в слоис­

упора (рис. 76). Функция

той толще с кривой депрессии

в

пределах

Гиринского

G имеет раз­

нескольких слоев

 

 

мерность расхода. Напри­

 

 

 

мер,

для

рассматривае­

мой схемы слоистого пласта (рис. 76) значение функции Гирин­ ского для начального и конечного сечения потока определяется со­ ответственно следующими выражениями:

Gi = kimi (hi — 0,5mi) + k2m2 (hi — — 0,5m2) +

+ k3m3 (hi — mi — m2— 0,5m3) + k4m4 (hi mxm2 — m3 — 0,5m4) +

+ knmn(hi — mi— «2

tn3 — m4 — 0,5mn) ,

(V,27)

Gz — kitni (h2 — 0,5mi) +

k2m2 (h2 — nil — 0,5m2) +

 

+ k3m3(h2— mi — m2— 0,5m3).

(V,28)

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

Gi Gz

(V,29)

Gx = G i - — ---- - x

Gi—z

Зная значение функции Gx и ее выражение в соответствии с формулой (Ѵ,26), можно определить и соответствующую функции Gx мощность потока hx. Для удобства определения hx по Gx пред­ варительно строится график G= f(h), который используется затем для графического определения hx по значению Gx.

Движение напорных вод в слоистом пласте переменной мощно­ сти. Неравномерное движение напорного потока в слоистом пласте

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

Г и mi “b т1

Н2

т2-|- т2

Ні Н2

Ц= l ki -----

^---------

;--------

 

М г

2

7-1—2

 

2

7-1—2

 

 

... + ka

гпп + m

n"

Hi

Н2Л

 

,

----------------—

— J = [fei( m i +

m i ) +

 

2

 

 

7-I_2

 

 

 

k2(m2

m2 ) + -

 

+

ka(mn -)- отп

H i - H 2

 

)]•

(V,30)

 

 

 

 

 

 

 

 

27-1—2

Рис. 77. Движение напорных вод в слоистом пласте переменной мощности

Здесь от/ и от/'— мощность рассматриваемого слоя в крайних се­ чениях 7 и 2.

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

о т " — о т '

т ср = :------,— ------7 .

In от" In от'

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

я = н _______________________ 2дх _____ _______________

h (m' + от") + k2(от' + m p + ... + kn(от' + от")

(V,31)

где Ях — ордината пьезометрической кривой в сечении, отстоящем на расстоянии х от начального.

П р и м е р . Определить расход грунтового потока в песчано-га- лечниковых слоистых отложениях, вскрытых скважинами 17 и 18, расположенными по направлению движения потока на 577 м одна от другой, а также отметку уровня грунтовых вод на расстоянии 177 м от скважины 17 по следующим данным. Уровень воды в сква­ жине 17 # і 7= 118,16 м, в скважине 18 / / 18 — 115,16 м. Отметка по­ верхности плотных глин, являющихся водоупором, в обеих скважи­

 

 

 

 

 

 

 

 

 

нах равна 101,57 м. На гли­

Ш. /7

 

іу= о

 

 

 

Скв. 18

нах

снизу

вверх

залегают:

 

 

 

 

1)

тонкозернистые

глини­

11816\

 

 

 

 

 

 

стые пески мощностью Ш\ =

%=16\ • ' X

■ •

'

Х

ч

:

.115161

= 5

м,

с

коэффициентом

фильтрации

ki = 2

м/сут;

О ,

о

 

 

 

 

 

 

h ?

О

к °

0

о

0 —

0

0

2)

крупнозернистые

граве­

Л=зо

0 1

О

О

О

 

“ о

о

Г

листые

пески

мощностью

- А = у ' .

 

 

 

 

 

г 10157

 

 

 

 

 

т 2 = 6,4 м с коэффициентом

 

 

 

 

 

 

 

 

 

 

 

 

*ж А1

^

І-0

 

 

 

 

 

фильтрации

&2= 30

м/сут;

 

 

 

 

 

 

 

 

3)

мелкозернистые

пески,

 

 

 

 

 

 

 

 

 

Рис. 78. Движение грунтовых вод в сло­

средняя

мощность которых

 

истой толще

(к примеру)

/гСр= 10

м,

коэффициент

 

 

 

 

 

 

 

 

 

фильтрации

&3 = 3,6

м/сут

 

 

 

 

 

 

 

 

 

(рис. 78).

 

 

 

Р е ш е н и е . Построив разрез, видим, что кривая депрессии грун­ тового потока расположена в пределах верхнего песчаного слоя. При этом мощность водоносной части слоя в скважине 17 (прини­

мается за первое сечение)

составляет

/іі = # і — (ml+ m2) — z =

= 118,10—11,40—101,57 = 5,19

м, а в скважине 18 /і2= 115,16—

11,40—101,57 = 2,19 ж.

 

(Ѵ,23), рассматривая от­

Расход потока определяем по формуле

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

и верхнюю с переменной мощностью:

 

 

 

 

h

h2

, 2 ,2

q — (kitîly +

k2m2)

hi h%

hl—2

2LI-2

 

 

 

5,19 —

2,19

3,6 X

5,192— 2,192

= (2 X 5 + 30 X 6,4) X -

577

+

2X 577

 

 

 

 

=

1,1192 «

1,12 м3/сут.

 

Имея значение расхода потока и считая его неизменным по пути движения (№ =0), легко определить положение уровня в сече­ нии, отстоящем на расстоянии 177 м от скважины 17. Для этого до­ статочно рассмотреть грунтовый поток в пределах его верхнего слоя:

 

 

 

,2

,2

 

 

 

h\ h2

 

 

 

 

h

hi

hx

 

 

 

 

 

 

 

 

2x

=

q (kitrii + k2m2)

 

 

 

 

 

 

 

 

 

 

—2

 

 

 

 

 

 

 

 

,

,2

,2

 

 

 

 

 

 

 

 

 

Лі — /І2

 

 

 

 

 

 

 

 

 

«3—^7----- , откуда

 

 

 

 

 

 

 

 

 

Zbi_2

 

 

 

 

 

 

 

 

 

 

 

5,192—2,192 X 177=4,49 M.

 

 

 

 

 

 

 

 

 

577

 

 

С учетом величины hx отметка уровня в заданном сечении # х=

= Ш\ tn2-ЬzxТ йх=11,40 4-101,57 -Т 4,49 = 117,46 м.

 

 

 

ДВИЖ ЕНИЕ ПОДЗЕМНЫХ ВОД В ДВУХСЛОЙНОМ

 

 

ПЛАСТЕ

 

 

 

 

 

 

 

 

Двухслойный

пласт

является

частным

случаем

много­

слойного пласта

(при п 2), наиболее часто встречающейся схемой

фильтрации

подземных

вод, особенно грунтовых.

В предыдущем

параграфе

были

рассмотрены

некото­

 

 

 

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

 

 

 

которых

решения

для

двухслойного

 

 

 

пласта могут быть получены как част­

 

 

 

ный случай. Ниже дается более стро­

 

 

 

гое обоснование решений для двух­

 

 

 

слойного пласта на основе рассмотре­

 

 

 

ния дифференциальных уравнений для

 

 

 

единичного

расхода

потока.

случай

 

 

 

Разберем

наиболее

частый

 

 

 

двухслойного

строения

пласта, когда

 

 

 

верхний

слой

менее водопроницаем,

Рис. 79. Движение

грунто­

чем нижний. Нижний слой имеет мощ­

вых вод в двухслойном пла­

ность Ш\ и коэффициент фильтрации h

 

сте

 

(тх = const,

k\ const);

верхний

слой

 

 

 

имеет переменную мощность, изменяющуюся от hx в первом сече­ нии до h2 во втором сечении, и коэффициент фильтрации k2. Отсчет

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

(так

как

і = 0, то Hx= h x и Н2 = Іг2). Расчетная схема фильтрации

пока­

зана

на рис. 79.

 

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

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

ренциальным уравнением Дюпюи:

 

 

q =

dh

k2h

dh

klnil

(V,32)

 

dx

 

dx'

В уравнении (V,32) первый член правой части — расход напор­ ного потока в нижнем слое с постоянной мощностью ті и коэффи­ циентом фильтрации ku а второй член — расход грунтового потока в верхнем слое с переменной мощностью h и коэффициентом фильт­

рации k2 (см. рис. 79).

 

 

 

 

 

Из уравнения

(Ѵ,32), разделив переменные, найдем:

 

 

 

q dx = kimi dh — k2h dh.

 

(V,33)

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

от се­

чения 1 до сечения 2 будем иметь:

 

 

 

 

q {х2Xi) =

kiiîii (hi — h2) + kz

h2 — h 2

(ѵ,34)_

1

2 .

Из уравнения

(V,34),

приравнивая x2 — X\ = L\-2, получим рас-

четную формулу:

 

hi — hz , .

hi + hz

hi — h2

 

q =

,

(Ѵ,35)

kirrii— -------------1- kz------- ------------------------- .

 

 

-Ьі-2

^

'

-^1-2

 

Аналогично (Ѵ,35) расход потока в любом сечении, расположен­ ном на расстоянии х от сечения 1, определится выражением:

qi-х = kitni

hi hx

(V,36)

 

X

 

Учитывая постоянство расхода потока по всем его сечениям (по условию задачи, W = 0 и, следовательно, q\-2qi-*=q) , можно при­ равнять правые части уравнений (Ѵ,35 и Ѵ,36) и найти значение hx, необходимое для построения кривой депрессии. Либо при известной величине расхода потока, определенной по формуле (Ѵ,35), можно определить значение hx при различных х непосредственно из урав­ нения (Ѵ,36).

Полученные расчетные формулы являются строго обоснованны­ ми при сравнительно небольшом различии коэффициентов фильтра-

! ki

\

ции верхнего и нижнего слоев \ —- ^ 1 0

). Удовлетворительные ре-

kt

 

шения получаются и при величине—до 100. Однако если-г^> 100,

kz

kz

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

висимости от соотношения напоров в рассматриваемых слоях верхний слой будет либо питать, либо дренировать нижний слой (см. гл. IX—X ).

Напорно-безнапорное движение по схеме двухслойного пласта.

Основываясь на схеме двухслойного пласта, можно легко получить решение для напорно-безнапорного движения с использованием функции Н. К. Гиринского. В этом случае значение функции Гиринского для начального и конечного сечений потока (см. рис. 65) определяется как для двухслойного пласта: в качестве верхнего слоя рассматривается перекрывающий водоносные отложения водоупор с коэффициентом фильтрации k2 = 0 и мощностью потока, ко­ торая определяется положением пьезометрической кривой над кров­ лей нижнего водопроницаемого слоя, имеющего мощность т и ко­ эффициент фильтрации k \ —k.

Найдем значение функции Н. К. Гиринского G для начального

(х = 0) и конечного

(x — L і_2) сечений, учитывая приведенное выше

выражение (Ѵ,26) для ее определения:

 

 

 

 

П

 

 

 

 

G =

2 кіЩ(h

- Zi).

 

 

 

 

І—1

 

 

 

Для начального

сечения

при ky = k, Zi = 0,5m

и k2~0,

 

Ну — т

получим:

 

 

 

z2 = т-\------------,

 

 

 

Gy = km (Ну — 0,5т) -f k2 (Ну — т ) ^ Н у — т +

~ 2 ~ “ )]

=

 

= km(Hl — 0,5m)i

 

( Ѵ , 3 7 )

Для конечного сечения при ky = k и Zi = Q,bh2 имеем:

 

Gz =

kh2(ht -

0,5ht) =

h 2

(V,38)

0,5khi = k-^--

Далее, учитывая, что расход потока определяется выражением

Gy Gz

 

 

 

 

 

Q — —i------- находим:

 

 

 

 

•Ь1-2

 

 

 

 

 

т(Ну — 0,5т)-

2mH у — т2hz

 

 

 

=

 

 

Ly-z

k

 

 

 

 

2Ly-z

 

 

т(2Ну — т) — h2

(V,39)

2hZz

Выражение для расхода потока (Ѵ(39) идентично формуле (ІѴ.57), которая была получена на основе метода фрагментов (см. гл. IV, стр. 138).

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