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

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

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

В случае двухслойного строения ложа, когда основной водонос­ ный пласт отделен от реки слабо-проницаемым -слоем мощностью h\ с коэффициентом фильтрации k u значение AL определяется по сле­ дующей формуле:

 

(VII,31)

где cth {В 1/

\ представляет собой гиперболический котан-

1 '

k2h2h j

гене, определяемый по таблицам; k2 и h2— коэффициент фильтрации

и мощность основного водоносного пласта

(рис. 99, б).

/

В

\

Для крупных рек и водохранилищу при — >

20 — 30 J значе­

ние AL определяется выражением:

(VII,32)

При слоистом строении ложа рек и водохранилищ (рис. 99, в) значение AL можно рассчитывать по формуле:

(VII,33)

где kr и kB— средневзвешенные значения коэффициента фильтра­ ции при движении воды по напластованию и перпендикулярно на­ пластованию (при горизонтальной и вертикальной фильтрации кР и kB определяются по формулам (Ѵ,9 и Ѵ,17), рассмотренным в гл. V.

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

Пусть в створе из двух скважин, заложенных по потоку в непо­ средственной близости от русла реки с неоднородным строением ложа, замерено положение уровня грунтовых вод в период отсут­ ствия инфильтрации (№ =0). Расстояние первой скважины от уре­ за реки — L\, расстояние между скважинами 1 и 2 L2 положение уровня в скважинах соответственно обозначим Нх и Н2 в реке Н0

(рис. 100). Тогда расход потока

q на участке между скважинами

1 и 2 может быть записан так:

 

Ы,

 

Я

2 -

(VII,34)

Ч = Т

 

 

Для определения расхода потока на участке от скважины I до реки напишем выражение, аналогичное приведенному выше (VII,34), введя в формулу параметр AL, учитывающий неоднород-

Рис. 100. Схема грунтового потока к определе­ нию AL

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

Q

т Ні -

Но

(VII,35)

U + AL

 

 

Учитывая неизменность расхода потока по пути его движения (№ = 0), приравняем правые части выражений (VII,34 и VII,35) и получим формулу для определения значения AL:

AL = HJ ~ HJ U - U

(VII,36)

п 2 — Т7і

 

Из приведенного ясно, что если ложе реки обладает дополни­ тельным фильтрационным сопротивлением вследствие неоднород­ ности его строения, то значение АL, учитывающее это сопротивле­ ние эквивалентным увеличением длины пути фильтрации, будет больше нуля. При несущественной неоднородности ложа реки его дополнительное фильтрационное сопротивление будет незначитель­ ным и величина AL, определяемая по формуле (VI 1,36), будет близкой к нулю. В гидравлическом отношении наличие дополни­ тельного фильтрационного сопротивления ложа реки приводит к до­ полнительной потере напора потока, которая затрачивается на пре­ одоление этого сопротивления. О величине потерь напора при этом можно судить по значению напорного градиента на участке выхода потока в реку. Существенное его увеличение на этом участке отве­

чает наличию значительного дополнительного фильтрационного со­ противления ложа реки. Отсюда ясно, что учет неоднородности ло­ жа водохранилища в расчетных формулах можно осуществлять либо путем соответствующего увеличения длины потока L (на AL), либо путем учета дополнительной потери напора АН. Значение па­ раметра АН может быть определено так же, как и значение пара­ метра ЛL, по данным о положении уровня в трех точках с исполь­ зованием приведенных выражений (VII,34 и VII,35).

Обозначая в формулах (VII,34 и VII,35) падение напора на участке Ь2 через АНХ- 2 а на участке Ьх через АЯ0-і и вводя пара­ метр АН, получим после исключения q следующее выражение:

ЛЯ = АЯо_1- Д Я '1_ 2^-.

(VI 1,37)

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

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

ПОДПОР ГРУНТОВЫХ вод в условиях

НЕУСТАНОВИВШЕЙСЯ ФИЛЬТРАЦИИ

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

Процесс фильтрации подземных вод при развитии подпора носит явно неустановившийея характер и описывается дифференциальным уравнением (II.88) Буссинеска. Прогноз развития подпора осуще­ ствляется на основе аналитических решений дифференциальных уравнений либо методом моделирования.

Аналитические методы основываются на линеаризации диффе­ ренциальных уравнений и их последующем решении при соответст­

вующих начальных и граничных условиях. К настоящему времени получены аналитические решения для одномерных линейных полуограниченных и ограниченных потоков при мгновенном, постепен­ ном и других видах изменения уровня воды на их границах. Имеют­ ся отдельные решения для двухмерных плановых потоков и слож­ ных условий инфильтрационного питания [5, 13, 22, 32, 33, 62, 70,

106 и др.].

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

расчетов. Предварительно

отметим,

Ф(Л)

что точность прогноза развития под­

 

пора на основе аналитических реше­

 

ний значительно менее высокая,чем

 

точность определения стационарно­

 

го подпора,

вследствие

недостаточ­

 

ной достоверности параметров (k и

 

ц), используемых в расчетах, и зна­

 

чительной схематизации

природных

 

условий при получении решений.

 

Неустановившийся подпор грун­

 

товых

вод в условиях

полуограни-

 

ченного потока. Решение для схемы

 

однородного полуограниченного по­

 

тока с горизонтальным

водоупором

 

(і = 0)

при отсутствии инфильтраци­

 

онного

питания

(117= 0)

и мгновен­

 

ном изменении

горизонта

воды на

 

урезе

водохранилища

 

получено

 

H. Н. Веригиным [5, 33]. Принимая

 

мощность

потока

в зоне

подпора

Рис. 101. График функции

равной

некоторой

средней величине

Ф(Я)

/іСр и решая линеаризованное таким

образом дифференциальное уравнение (VI,2), H. Н. Веригин полу­ чил следующую расчетную формулу для определения ординаты кривой депрессии в процессе развития подпора во времени:

У* = УЬ\ + ( У \ - І Ч ) [ \ - Ф ( Х ) ] ,

(VII,38)

где ух— искомая ордината кривой депрессии в сечении, располо­ женном на расстоянии х от уреза водохранилища через время t, считая от момента заполнения водохранилища; hx— мощность по­

тока в расчетном сечении до подпора; hi и у\ — мощность

потока

на урезе водохранилища до и после подпора (см. рис. 92);

Ф(к)

специальная функция (интеграл вероятности Гаусса), значение ко­ торой определяется в зависимости от величины безразмерного ар­ гумента і по графику (рис. 101).

Величина аргумента К определяется выражением:

X

X

 

X =

(VI 1,39)

 

2 Уat

 

где k — коэффициент фильтрации;

ц — недостаток насыщения

по-

,

khcv

род; пСр — средняя мощность потока в зоне подпора;« = ------

 

ц

 

коэффициент уровнепроводности.

Таким образом, аргумент X учитывает положение сечения, в ко­ тором определяется подпор, время его развития и параметры пото­ ка. Значение средней мощности потока hcv можно определять по формуле:

/іср

2j/i + hi

(VII,40)

3

 

 

Для периода прекращения фильтрации воды из водохранилища

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

.

Уі-^-ha

 

находится по отношению «ср «

— - — (здесь пп— мощность по­

тока грунтовых вод в сечении, где влияние подпора практически не ощущается).

Как видно из графика функции

Ф(Л), представленного на

рис. 101, при ^~оо аргумент

Ф (Я =0)=0.

При этом

формула

(VII,38) обращается в формулу

(VI 1,6)

стационарного подпора.

Расход грунтового потока в любом сечении на расстоянии х от

уреза водохранилища в момент

времени t

определяется

выраже­

нием:

 

 

 

 

k (угhl)

—j — ,

1Ѵ7

1

(VII,41)

 

 

 

 

2 У nat

Из формулы (VII,41) следует, что расход потока на урезе водохра­ нилища (при х — 0) можно определять по уравнению (VII,41а):

Ц уі

ht)

(VII,41а)

Ях=0

<7о>

2 фnat

где (ft — расход потока на урезе водохранилища до подпора (опре­ деляется по формулам гл. IV и V ).

В работах [5, 29, 106, 108] приведены решения для определения неуетановившегося подпора в полуограниченном пласте при линей­ ном и ступенчатом законах изменения уровня в водохранилище, а также для условий неоднородного .-троения потока.

Пользуясь приведенными расчетными формулами (VII,38 и VII,41), можно построить кривые депрессии и определить расход потока на разные периоды времени t от начала развития подпора.

П р и м е р (по H. Н. Биндеману). Определить положение кри­ вой депрессии в процессе развития подпора. Участок сложен мел­ козернистыми песками, имеющими коэффициент фильтрации k — = 4,77 м/сут, недостаток насыщения р = 0,20. Водоупорное ложе по­ тока горизонтально и 'имеет отметку 0. Мощность водоносного пласта у реки до подпора Ді = 5 м, после наполнения водохранили­ ща у і= 12 м.

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

отстоящих на 50, 100, 250 и 500 м от уреза

водохранилища

через

50,

100 и 250 суток от момента наполнения

водохранилища

и при

t =

оо («стационарная» кривая).

 

 

Порядок расчета дается на примере определения мощности по­ тока в сечении, отстоящем на 100 м от водохранилища, через 250 суток от начала его наполнения. Мощность водоносного пласта до подпора в этом сечении определена /іюо= 6,98 м.

По формулам

(VI1,40 и VII,39) определяем значения /гср и Я: '

/іср 2Х 12 +

5 =

9,67 м и

Я =

 

100

 

0,208.

 

'

0,20 X 9,67 X 250

 

: 1 / І Ѵ

По графику (см. рис. 101)

находим значение Ф(Я) —0,235. Наконец,

вычисляем ух по формуле

(VI,38): ух2 = 6,982+ (122—52) (1—0,235) =

= 139,75; г/х = У 139,75= 11,82 м. В табл. 5 приводятся данные о мощности водоносного пласта для различных сечений и моментов времени, полученные аналогичным расчетом. Изменения кривой депрессии грунтовых вод во времени показаны на рис. 102.

регаводохрани­Расстояниеотбе­ лища,м

X

 

 

50

ос

100

ос

250

оо

500

оо

 

 

 

Время от началнаполненияводохранилища /, суткиТ а б л и ц а

5

 

0

 

50

 

100

 

 

250

 

ео

(стационар­

 

 

 

 

 

 

 

 

 

 

ныйподпор)

г<

у

X

Ф(А)

У

X Ф(Л)

У

X

Ф(\)

У

X

(>4

V

%

 

 

 

 

 

ф

 

1 6,070,2350,25511,200,1650,18511,570,1040,11511,920 012,48

1 6,980,4700,49510,430,3290,36211,170,2080,23511,820 0 12,95

1 9,101,1750,8809,850,8240,75010,610,5210,53711.740 0 14,26

112,002,3500,99912,011,6281,98412,041,0420,85312.750 016,22

Неустановившийся подпор грунтовых вод в ограниченном пото­ ке. Для ограниченного по длине потока получены решения при различных граничных условиях на второй (правой) границе потока:

а) для случая открытой второй границы и выполнения на ней граничного условия первого рода H = f(t) ;

б) для случая закрытой второй границы и выполнения на ней граничного условия второго рода (Q = const).

Ограниченный поток с постоянным уровнем на верховой грани­

це. Решение для условий

однородного грунтового потока

с гори­

зонтальным водоупором

(і = 0) и наличием

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

(ІѴ>0)

при мгновенном изменении уровня воды в

водохранилище

(мгно­

венный подпор) получено H. Н. Веригиным [5, 33]. Расчетная фор­ мула для определения ординаты кривой депрессии в процессе раз­ вития подпора в любом произвольном сечении на расстоянии х от

уреза водохранилища имеет вид:

 

 

 

г

U-

-X

 

(VII,42)

У:■= Ÿ Лх+ ({/ï- ■AÏ)

 

L\—2

Z-1—2

 

 

где Z.J—2— ширина междуречья

 

 

 

X

 

(длина

потока) А

%

 

Li-2

специальная функция (ряд Фурье), значения которой определяют-

 

 

X

и величи-

ся в зависимости от положения расчетного сечения -,----

khCpt

 

I. 1—2

 

at

 

 

ны параметра t =

L2

по графику (рис. 103).

Значения

 

1-2

 

Рис. 102. Схема формирования уровня грунто­ вых вод во времени при подпоре

Six, —---- ) отсчитываются по оси ординат для соответствующей

по значению т кривой.

При £ = оо значение S (т, —---- ) — 0 и приведенная выше фор-

'Х-1-2 '

мула (VII,42) становится аналогичной формуле (VII,11) для опре­ деления стационарного подпора.

Величина фильтрационного расхода потока на урезе водохрани­ лища (х = 0) в любой момент времени t определяется по формуле (VI 1,43):

 

 

 

2

^2

 

 

 

 

 

9х=о = * - ^ — î-[Sq(T)+l]+<7o.

(VI М3)

 

Значения функции 5q(t) приведены ниже [5]

 

 

1.

. .

О

0 , 0 2 0 , 0 3 0 , 0 4 0 , 0 5 0 , 0 7 5 0 , 1 0 0 , 1 5 0 , 2 0 , 3 0 , 4

0,5 0,7

оо

S

(т)

оо

3 2, 251, 82 1 , 5 2 1 , 0 8

0 , 7 8 0 , 4 6 0 , 2 8 0 , 1 0 , 0 3 9 0 , 0 1 5 0 , 0

4 5 0

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

Так, например, при мгновенном изменении уровня воды на ле­ вой и правой границах ордината кривой депрессии ух в расчетном сечении определится следующей формулой, полученной из выраже­ ния (VII,42):

 

 

 

-------------\--------

+ (У2

---------s (

L i - 2 — X VI

т,

(ѴІІ,44)^

L ^1-2

'

L i - 2

" ‘

Значение функции 5 (т,

- —

) определяется

с учетом начала

 

L i - 2

 

 

 

развития подпора на каждой из границ отдельно. Все обозначения формулы (VII,44) ясны из рис. 92.

Рассмотрим пример расчета по формуле (VII,42).

П р и м е р (по Н. Н. Биндеману).

Требуется

определить мощ­

ность грунтового потока при подпоре

в сечении,

находящемся на

расстоянии %= 500 м от уреза водохранилища, через 100 суток пос­ ле наполнения водохранилища. Мощность водоносного пласта в этом сечении до подпора была h50о= 13,75 м, в начальном сечении /іі = 12 м. Превышение НПГ водохранилища над водоупором у\ — = 18 м. Водоупор горизонтальный. Коэффициент фильтрации плас­ та, сложенного среднезернистыми песками, равен £=10 м/сут, не­ достаток насыщения ц = 0,20. Расстояние уреза водохранилища до соседней реки Т[_2 = 200 м.

Определим /гср по формуле (VI 1,40): 2 X 1 8 + 1 2

hep

16 м.

"

3

Вычисляем значение т по формуле (VII,42), откуда:

=

khcvî

= 10 X 16 X 100

 

т

(XL ^ 2

0 ,2 0 X 2 0 0 0 2

’ '

По графику (рис.

103)

при L1—2

 

500

=

0,25 находим значе­

 

2000

ние функции:

 

 

X

\

 

 

 

 

S

/

0,54.

 

 

 

т,

------) =

 

 

 

 

 

4-1-2

 

 

 

 

Определяем по формуле

(VI 1,42)

мощность водоносного пласта

при подпоре в расчетном сечении ух:

 

 

 

у і = 1 3,752 + (

1 8 * -

1 2 *

) X ( ' ~

^ Ö

5 Ö T ~

-

° - 5 4 ) = 2 2 6 ’ 8 ;

 

Ух =

У 226,8 =

15,06 м.

 

Следовательно, в указанном сечении за 100 суток уровень грун­ товых вод повысится на: ÿsoo—^5оо= 15,06—13,75=1,31 м.

Неустановившийся подпор грунтовых вод в ограниченном пото­ ке с наклонным водоупором. При наклонном залегании однородно­ го водоносного пласта и постоянстве расхода грунтового потока

до и после подпора мощность водоносного пласта ух в процессе развития подпора можно вычислять по формуле В. М. Шестакова (см. расчетную схему на рис. 95):

ух =

hx -\- іх (/гх

0,25іх) - f - (i/i

hi) Fo{^ T , —j

0,5ix,

 

 

 

 

 

(VII,45)

где I — уклон водоупора,

принимаемый

положительным при пони­

жении его к водохранилищу и отрицательным при

понижении от

»водохранилища;

L — расстояние

 

 

от уреза водохранилища до уда­

 

 

ленной

границы,

где

подпор

 

 

практически затухает и мощность

 

 

потока hn является неизменной

 

 

(т. е. L — длина потока,

hn — на­

 

 

пор на его границе).

Значение функции (т’ т )

определяется по графику, изобра­ женному на рис. 104, причем ве­ личина % вычисляется по фор­ муле:

 

khçpt

(VII,46)

 

 

(lL2 ’

 

 

 

 

 

где

/гСр — средняя мощность грунтового потока,

 

 

 

,

Уі + hi + hn

(VII,47)

 

 

«cp =

--------- ------------•

.

Формула (VII,45) используется для расчетов

при значении

2hi

 

 

 

1 <

-J— (см. подробно работы [5, 27 и др.]).

 

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

= £ +

т )],

(VII,48)

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