![](/user_photo/_userpic.png)
книги из ГПНТБ / Климентов П.П. Динамика подземных вод учеб. для геологоразведоч. техникумов
.pdfВ случае двухслойного строения ложа, когда основной водонос ный пласт отделен от реки слабо-проницаемым -слоем мощностью 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 — на |
|
|
пор на его границе).
Значение функции Fо (т’ т )
определяется по графику, изобра женному на рис. 104, причем ве личина % вычисляется по фор муле:
|
khçpt |
(VII,46) |
|
|
|
(lL2 ’ |
|
||
|
|
|
|
|
где |
/гСр — средняя мощность грунтового потока, |
|
||
|
|
, |
Уі + hi + hn |
(VII,47) |
|
|
«cp = |
--------- ------------• |
|
. |
Формула (VII,45) используется для расчетов |
при значении |
||
2hi |
|
|
|
|
1 < |
-J— (см. подробно работы [5, 27 и др.]). |
|
Неустановившийся подпор грунтовых вод в ограниченном пото ке с постоянным расходом на верховой границе. Такая схема отве чает условиям прислонения водоносных отложений речной доли ны к коренным слабопроницаемым отложениям берега или цо кольной террасы (рис. 105). Расход воды, поступающей через такую границу, является постоянным (в частном случае Q = const = 0). Расчетная формула для определения ординаты кривой депрессии при мгновенном заполнении водохранилища имеет вид:
= £ + |
т )], |
(VII,48) |