Файл: Аполлов, Б. А. Курс гидрологических прогнозов учебник.pdf

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

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

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

Добавлен: 15.10.2024

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

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

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

Совместное решение уравнений (2.XI) и (З.ХІ) дает

 

А //= - і -

(В At — А А/ —АС),

(4.XI)

где Ѵо —п — ц — объемная влажность грунта

над

капиллярной

зо­

ной горизонта грунтовых вод.

 

 

 

 

 

 

 

Перепишем (4.ХІ) в виде

 

 

 

 

 

 

 

АНѴ$

і

/I

-|-

АС

 

ѵ п

В At

1

 

В At

 

(р .л і;

Обозначим величину, стоящую в правой части

(5.XI), через ß и по­

ложим, что она меняется мало. Тогда

 

 

 

 

 

Д Я = в 4 ^ .

 

(6.XI)

 

 

 

 

 

 

 

Величина ß может быть определена по данным о средних вели­ чинах АЯ, Ѵо и В для рассматриваемой площади. В гидрогеологии величина B A t вычисляется по эмпирическим формулам. Например, для условий гидрогеологического стационара в Щемнлове (под Мо­ сквой) для весны получено выражение

В М = л 'Х|_|ц — £ Хі-ііі ~Ъ (-Hv — ^іѵ)>

где х'хі- ііі и £ хі- ііі — осадки и испарение за зимние месяцы с но­ ября по март; лдѵ и В \ \ — то же за апрель.

Для стационара среднее значение Ѵ0 = 0,137 и

ß = 0,4. Отсюда

АН—2,92 ( лГхі- ііі — ^ х і - н і ~ г Л‘іѵ — ^ iv ) -

( 8 . XI)

Пользуясь (8.XI), можно давать долгосрочный

прогноз весен­

него повышения уровня грунтовых вод в районе названного стацио­ нара. Аналогично многим прогнозам стока неизвестные величины осадков и испарения, в данном случае за апрель, приходится при­ равнивать к их средним многолетним значениям. Как мы уже ви­ дели в гл. ѴП, при прогнозах весеннего стока потерн талых вод определяются более точно, чем это можно сделать с помощью вы­ ражения (7.ХІ).

При прогнозе высоты подъема грунтовых вод весной можно ис­ ходить из того, что

АH {= ^ L P i ,

(9.XI)

Р

 

где АН і и pi — подъем уровня и потери талых вод в данном году (потери вычисляются на основе ожидаемой величины весеннего

стока и данных, использовавшихся при прогнозе этого стока); АН

и р — то же средние многолетние.

Рассмотренные способы прогноза величины АН, конечно, при­ ближенны и причины этого достаточно ясны из изложенного.

Более точного решения задачи долгосрочного прогноза весен­ него подъема уровня грунтовых вод можно достичь путем отыска­

334


ния эмпирических зависимостей на основе уравнения (2.XI) п со­ ответствующих данных многолетних наблюдений. Остановимся на этом вопросе.

Первый член правой части уравнения (2.XI) зависит при посто­ янном значении A t от начального уровня грунтовых вод Я0, вели­ чины потерь талых вод р и влажности почво-груита в зоне аэрации \Ѵа. Отсюда

ДH = f ( H 0, р , W a) .

При переменной продолжительности времени подъема уровня грун­ товых вод A t имеем

A H — f ( Н 0, р , W a, A t).

Способы расчета и прогноза величины р были изложены в гл. VII. На значительных территориях после наступления весеннего мак­ симума уровня грунтовых вод начинается плавное закономерное его понижение. Это дает большие возможности для долгосрочного прогноза хода уровня грунтовых вод на протяжении лета и даже осени, а также зимы. То же относится к подземному стоку рек;

об этом говорилось в гл. VI.

§3. ПРОГНОЗЫ МЕТОДАМИ ДИНАМИКИ ПОДЗЕМНЫХ ВОД

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

Когда питание грунтовых вод за счет инфильтрации отсутствует, уравнение неустановившегося движения их можно записать в виде

ж=--Ѣ[*<"+"> S +4 - <*+"> Щ] • <іо-х|)

где т — недостаток насыщения грунта в слое над капиллярной кай­ мой; k — коэффициент фильтрации; — глубина водоносного слоя; Я — уровень грунтовых вод (в месте выхода источника Я = 0); t — время; X и у — расстояние соответственно вдоль и поперек грунто­ вого потока.

Если Я мало по сравнению с /г, то (10.Х1) примет вид

дН

д

дН

д_

 

 

 

(И.XI)

т —тт“= —^— kh

дх

ду

 

 

 

dt

дх

 

 

 

 

Уравнение (11.XI) можно переписать в виде

 

 

 

 

дН ии ( д2Н

J д2Н \ ,

д (kh)

dH ,

д {kh)

 

dH

(12.XI)

!П~ д Г ~ ПП[~ дЖ

~~df~)~x

д х ~

' 1 F " 1“

ду

ду

 

Величины kh, d(k h)/dx, d ( k h ) / d y зависят от координат х и у, а не от времени. Поэтому когда все Я меняются пропорционально, то и дН/дх, дН/ду, д2Н /дх2 и д2Н / д у 2 меняются в одинаковом отношении;

335


значит, и m d H / d t меняется в таком же отношении. Отсюда следует, что т д Н /dt пропорционально Я и

Н = Н йе ~ ч ,

(13.XI)

где Но — начальный уровень грунтовых вод.

ему пропорцио­

Так как уровню Я пропорционален уклон, а

нальна скорость, то при отсутствии пополнения

грунтовых вод за

счет инфильтрации расход источника будет, по Бусспнеску и Майе,

Q = Q 0e - “'f

(14.XI)

где Q — расход воды источника в момент времени t, отсчитывае­

мого от начала бездождного периода;

Q0 — расход воды в начале

этого периода; а — коэффициент истощения.

Кривую, являющуюся графиком (14.XI), часто называют кривой

истощения. Формулы (13.XI) и (14.XI)

справедливы, когда мощ­

ность водоносного слоя /г, из которого

вытекают источники, на­

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

Буссинеск и Майе получили приближенное решение также для случая, когда мощность водоносного слоя Іг невелика по сравнению с изменениями уровня грунтовых вод Я и поэтому не может быть

принята постоянной. Тогда

 

 

 

 

 

 

 

Н

(1+Щ/)2

 

 

 

(15.ХІ)

 

 

 

 

 

 

 

 

Q— Qo(l +Яі02 .

 

 

 

(16.XI)

где аі — тоже коэффициент

истощения, но

имеющий

уже

другое

значение.

 

 

 

 

по формуле

Значение коэффициента истощения а определяется

 

 

 

а — *2kh

 

 

 

(17.XI)

 

 

 

■ірР

 

 

 

 

а коэффициента аі — по формуле

 

 

 

 

 

 

 

kV

 

 

 

(18.XI)

 

 

а і — 6 , 7 7 4 а /3 ’

 

 

 

где k

и /г— коэффициент фильтрации и средняя мощность водонос­

ного

слоя; р — водоотдача

грунта; I — расстояние

от

водораздела

грунтовых вод до

выхода источника; V — объем водоносного слоя

на протяжении /.

Для потока единичной

ширины

V = h l .

Тогда

из (17.XI) и (18.ХІ) имеем

 

 

 

 

 

 

 

 

а= 1,71а1.

 

 

 

(19.XI)

В табл. 14 приведены ординаты кривых истощения, отвечающих формулам (14.ХІ) и (16.ХІ), выраженные в долях величины Qо при аі = 0,01. Из нее следует, что обе формулы дают практически одн-

336


наковые величины до момента снижения расхода воды в 2,0— 2,5 раза по сравнению с начальным Qo. Лишь в дальнейшем появ­ ляются различия, но и то не очень большие, причем формула (14.XI) дает меньшие величины.

Т а б л и ц а 14

Ординаты кривых истощения

t

е—0,0171/

1

t

 

 

(1+0,01 <)2

 

1

0,98

098

100

10

0,84

083

200

50

0,42

044

500

70

0,30

035

 

1

о

"

 

О

0,18

0,03

0,00

1

(1 -і-о.оі о2

025

011

003

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

Q = Qo [(1 - 7 ) е~’’71а<+ Т ( і Л р я ] .

(20.XI)

где параметр у может меняться от 0 до 1.

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

Q=

(Qo —<7г) [ ( 1 - т ) е - 1,71аг+ 7

+

(21 .XI)

Заметим,

что формулы (14.XI), (16.XI) и

(21 .XI)

приводились

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

спада половодья

(и паводков) и в

течение

лета для

прогноза

стока рек.

а 1г у и qr

можно определять

следующими спосо­

Параметры а,

бами.

(14.XI),

получим

lnQ = lnQo— at.

Следова­

Логарифмируя

тельно, а найдется как тангенс угла наклона прямой на графике связи величины ln Q с t.

Так как кривые, отвечающие формулам (14.XI) и (16.XI), прак-

Q

тически совпадают при -——-3=0,7, то параметр а можем наити по Ѵо

формуле (19.XI). Молено определять его и с помощью графика связи величин ) ИІ.

Параметры формулы (21.XI) находятся следующим образом. Вначале по данным многолетних наблюдений за стоком данной реки находится величина qr. Приближенно она может быть опреде­ лена также с помощью графика на рис. 132. Затем, пользуясь

22 Зак. № 708

337


рассмотренными приемами нахождения а н см, определяются их числовые значения в формулах Q qr= (QoЯг)е~а1 и Q qr=

Qо qr

и строятся вместе графики этих формул при найден­

 

(1+cCiO2

ных значениях а и см. Здесь же вычерчивается график Q (t) по фактическим данным. После этого величина у находится как отно­ шение отрезка ab к ас (рис. 133).

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

а

Q-q

щий способ определения пара­

Рис. 134.

Схема пополнения п ис­

 

 

 

метра

у.

 

1, 2

и

3

— кривые истощения

соответ­

тощения в бассейне запасов под­

расходом

в замыкающем створе

 

 

(14.X I),

(16.XI) и

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

ственно

по формулам

 

 

фактическая.

Пока мы рассматривали только процесс истощения запасов под­ земных вод. Однако кривые истощения могут быть с успехом ис­ пользованы и для анализа всего процесса пополнения и истощения запасов подземных вод. Исходной посылкой для этого является тот факт, что пополнение запасов подземных вод происходит сравни­ тельно быстро, а истощение подземных вод является безусловно преобладающим во времени процессом (рис. 134). Пополнение мо­ жет рассматриваться как воздействие на динамическую систему — бассейн, а истощение — как ее ответная реакция. На рис. 134 пери­ оды пополнения (h-ti) и (А— /з), а периоды истощения от 0 до tu {ts— k) и ( h ti). Пользуясь кривыми истощения и исследовав связь между величинами поверхностного стока и пополнения под­ земных вод, Т. С. Абальян и Г. П. Калининым был разработан спо­ соб выделения на гидрографе реки подземной составляющей стока.

338