книги из ГПНТБ / Климентов П.П. Динамика подземных вод учеб. для геологоразведоч. техникумов
.pdfв рассматриваемом сечении. С учетом этого соотношения дифферен циальное уравнение (IV,15) будет иметь вид:
( dh |
dz |
(IV,16) |
q = — khy |
dS ) . |
|
dS |
|
В уравнении (ІѴД6) переменными являются величины h, z и S поэтому его непосредственное решение вызывает значительные за труднения.
г
Рис. 55. Схема к определению расхода грунтового потока с наклонным водоупором
Строгое решение приведенного уравнения (IV, 16) с учетом из менения всех входящих в него величин получено H. Н. Павловским
185] для случаев прямого (/>0) и обратного |
(КО) |
уклонов |
водо |
упорного ложа (см. рис. 55 и рис. 56). |
простое |
решение |
для |
Приближенное, но достаточно точное и |
|||
оценки условий движения потока грунтовых |
вод в пластах |
с на |
|
клонным |
водоупором получено |
||
Г. Н. Каменским [56]. |
|
Рис. 56. Расчетная схема к построе нию кривой депрессии грунтового по тока с наклонным водоупором
Определение расхода подзем ных вод по Г. Н. Каменскому.
Для определения расхода подзем ных вод при фильтрации их в ус ловиях однородного пласта с на клонным водоупором Г. Н. Камен ский считает возможным восполь зоваться известной формулой Ж. Дюпюи (ІѴ,9), выведенной для определения расхода потока грунтовых вод при горизонталь ном залегании водоупорного ло-
жа, учитывая при этом, что значения напоров потока должны от считываться от плоскости сравнения. Тогда формула (IV,9) при обретает вид:
q = |
,h i + h2 |
Ih + Uz |
(IV,17) |
k -------- - |
LI-2 |
||
|
|
|
|
hi -f- h2 |
|
|
безнапорного |
где --------- представляет собой среднюю мощность |
|||
водоносного горизонта в пределах |
рассматриваемого |
участка по |
|
тока. |
(IV.17) |
может быть обоснована более |
|
Приведенная формула |
строго, если для решения задачи воспользоваться гидравлическим методом. Для этого рассмотрим участок грунтового потока с на
клонным водоупором (і> 0), ограниченный двумя |
вертикальными |
|
сечениями 1 и 2, в которых мощность потока |
имеет |
значения hi и |
h2 а пьезометрический напор соответственно |
Нх и Я2. Значения |
пьезометрических напоров измеряются от горизонтальной плоско сти сравнения 0—0, совпадающей с осью х принятой системы ко ординат (см. рис. 55). Питание потока в пределах участка отсутст вует и его расход является постоянным по всем сечениям (<7= const). Принимается, что при малых значениях уклонов потока и водоупорного ложа в качестве длины пути фильтрации рассмат ривается его горизонтальная проекция (по оси х).
В этих условиях единичный расход грунтового потока определя ется следующим дифференциальным уравнением Дюпюи:
q ^ - k h ~ . |
(IV,18) |
ах |
|
Разделяем переменные (в данном случае х и |
Н) и интегрируем |
это уравнение (IV,18) с учетом изменения |
переменных величин в |
|
пределах ограничивающих сечений: |
|
|
Хі |
Н2 |
|
\ ß - d x = |
— \ dH. |
(IV,19) |
Jkh
Влевой части уравнения величина h также переменная и зави сит от значения х. Применяя теорему о средней, будем считать, что
1 |
„ |
u |
hi -f- h2 |
мощность потока h равна его средней величине |
«ср = |
------------ ; |
|
в этих условиях интеграл левой части уравнения |
(IV, 19) определит |
||
ся следующим образом: |
|
|
|
Я |
-(xz — xi). |
|
(IV,20) |
khcp |
|
||
kh cp |
|
|
После интегрирования обеих частей уравнения (IV,19) получим:
!L— [xz - x |
1) = |
Hl - H |
z. |
(IV,21) |
Учитывая, что х2—Х і=Д _2, |
а |
и |
hi~\~ hz |
формулы |
«ср —>---- ----- , из |
(IV,21) найдем выражение для определения расхода потока с на клонным водоупором:
q = |
, hi |
hz Hi - Я2 |
(IV,22) |
k --------- |
|||
|
|
h 1 - 2 |
|
Формула (IV,22) является |
приближенной |
вследствие замены |
при интегрировании переменной мощности потока h ее средним зна чением hCp. Однако во многих, практически важных примерах, где движение подземных вод происходит со свободной поверхностью при наклонном залегании водоупорного ложа, выполненные расче ты по приближенной формуле Г. Н. Каменского не показали суще ственных отклонений по сравнению с расчетами по строгим форму лам H. Н. Павловского.
Построение кривой депрессии. Положение уровня грунтовых вод при наклонном водоупоре в любом сечении потока определяется величиной напора, отсчитываемого от плоскости сравнения. Следо вательно, для построения кривой депрессии надо располагать вы ражением для определения напора Нх в искомом сечении потока. Это выражение может быть получено на основе выведенной форму лы для расхода потока подземных вод (IV,22). Для этого достаточ но составить два уравнения для определения расхода потока по двум различным парам сечений, в одну из которых входит и сече ние с искомой величиной напора Нх и, приравняв эти уравнения (в силу равенства расхода), получить формулу для определения искомой величины Нх.
Пусть, например, требуется определить значение пьезометриче ского напора Нх в сечении, расположенном на расстоянии х от лево го ограничивающего поток сечения. Значения мощности потока и пьезометрического напора в ограничивающих поток сечениях 1 и 2 считаются известными и соответственно равными hu Hi и /г2, Н2 (см. рис. 56).
В соответствии с выше приведенной формулой (IV,22) напишем выражения единичного расхода для пар сечений 1—2 и 1—х:
q1—2= |
hi -f- hz |
Hi — Hz. |
k -------------- -------- и |
||
Qi- |
hi -f- hx |
Hi — Hx |
k |
(IV,23) |
Учитывая, что q\ - 2 = q\~x, приравниваем правые части |
уравне |
|
ний и после сокращений найдем: |
|
|
(hi + hi) H l ~ H- 2 = |
(hi + hx)Hi ~ H x . |
(IV,24) |
4 . 1 - 2 |
X |
|
В уравнении (IV,24) два неизвестных: hx и Ях. Одно неизвестное можно исключить, для чего мощность пласта кх можно заменить разностью отметок уровня воды Нх и поверхности водоупорного пласта zx, т. е. hx= Hx—zx (см. рис. 56).
Подстановка найденного значения hx в формулу (IV,24) дает:
( h + ht) Н і~ Нг = ( h + Нх - zx)^—^ . |
(IV,25) |
|
4 - 1 - 2 |
X |
|
Решая уравнение (IV,25) |
относительно Я х, можно определить |
|
отметку уровня воды в любом сечении потока. Для этого |
помимо |
значений напора и мощности потока в двух известных сечениях не обходимо еще знать превышение водоупора в искомом сечении над плоскостью сравнения (zx=ix).
При построении кривой депрессии потока можно для определе
ния ее ординат пользоваться следующим уравнением [29]: |
|
hx = ~\lh?— іх (hi — 0,25ûc)--- ^-x — 0,5tx. |
(IV,256) |
Это уравнение получается из формулы (IV,23), в которой вели чина Я х представлена через разность Hx= hx+ ix (і — уклон водо упора). При этом плоскость сравнения проводится через водоупор в
первом сечении, где х = 0 (см. рис. 56). |
|
||||
Для определения ординат по выше |
|
||||
приведенной формуле |
(IV,256) |
пред |
|
||
варительно находится значение |
расхо |
|
|||
да потока q по формуле |
(IV, 22). |
|
|||
При наклонном |
залегании |
водо |
|
||
упорного ложа форма кривой депрес |
|
||||
сии зависит от изменения мощности по |
|
||||
тока. Принято различать потоки грун |
|
||||
товых вод с прямым и обратным укло |
|
||||
ном водоупорного ложа. При прямом |
Рис. 57. Кривая подпора при |
||||
уклоне водоупорного |
ложа направле |
||||
прямом уклоне водоупорно |
|||||
ние движения воды и уклон ложа сов |
го ложа |
||||
падают (г> 0), при обратном — проти |
|
воположны ( і< 0). По Н. Н. Павловскому при прямом уклоне во доупорного ложа возможны две формы кривой депрессии: кривая спада и кривая подпора (см. рис. 55 и 57). При обратном уклоне (і<0) возможна лишь кривая спада (см. рис. 56). При кривой спа да мощность водоносного горизонта по направлению движения грунтовых вод уменьшается, а при кривой подпора — увеличивается (рис. 57).
О характере влияния мощности потока на форму кривой деп рессии можно судить по изменению напорного градиента. Из об щей формулы закона Дарси (IV,3) выражение для напорного гра диента имеет вид
|
/ = j —. |
|
(IV,26) |
|
При неизменном расходе |
потока |
(q = const) и |
постоянном |
ко |
эффициенте фильтрации |
(&=const) |
напорный |
градиент, |
как |
это видно из формулы (IV,26), зависит от изменения мощности потока. При уменьшении мощности потока по направлению движе ния напорный градиент увеличивается, и кривая депрессии приоб-
Рис. 58. Расчетная схема грунтового потока с перемен ным уклоном водоупорного ложа
ретает выпуклый характер, что типично для кривой спада (см. рис. 56). При увеличении мощности потока по пути движения на блюдается обратная картина, и кривая депрессии приобретает во гнутый характер, что свойственно для кривой подпора (см. рис. 57).
Движение грунтовых вод при переменном уклоне водоупорного ложа. Если водоупорное ложе на различных участках потока име ет различные уклоны и направления падения, то для решения зада чи по построению депрессионной кривой и определению расхода потока необходимо рассмотреть последовательно все участки от дельно с учетом граничных условий в точках перегиба водоупорно го ложа. Пусть, например, водоупорное ложе имеет три участка с различным уклоном и направлением падения водоупора (рис. 58). Известными являются значения мощности потока и напора в край них, ограничивающих поток сечениях {h\, Н\ и hit Я4) и отметки водоупора в точках его перегиба (z2 гъ) . Требуется определить отметки уровня грунтовых вод в промежуточных сечениях потока (Я2, Я 3) и величину его расхода.
Через точки перегиба водоупорного ложа проводим вертикаль ные сечения 2 и 3. Расстояния между всеми сечениями считаются известными (Т-і_2, 7-2_3, Т-з-Д. Теперь, используя известную форму лу для расхода потока с наклонным водоупором (IV,22), напишем выражения для единичных расходов на каждом из его участков с одинаковым уклоном водоупорного ложа с учетом значений мощ ности потока и напора на границах участков 1—2, 2—3, 3—4.
hi + hz H i - H z
2 7-1-2
h2+ h3 Hz - H 3
2 7-2-3
h3~Г /14 H3- H ,
(IV,27)
27-3-4
Вуравнения (IV,27) входят неизвестные в сечениях 2 и 3 мощ
ности (h2 hz) и напоры потока (Я2, Я 3), связанные следующими соотношениями (см. рис. 58): H2 = h2 + z2\ H3 = hz + z3.
Учитывая это соотношение, исключим из уравнений (IV,27) не известные величины h2 и hz, выразив их как разность отметок уров
ня воды и водоупора (h2=H2—z2; h3—H3—z3): |
|
|||
|
Яі —(—Я2— z2 |
Hi — Н3 |
|
|
ai_о = k ------------------------ ------------ |
|
|||
H2— Zz -f- Я3— z3 H2- H 3 |
|
|||
q2-3 = |
|
2 |
7-2-3 |
|
|
|
|
||
^3-4 = |
H3— z3-f- hi |
Я3- Я 4 |
(IV,28) |
|
k |
2 |
L3-4 |
||
|
|
|
||
Уравнения (IV,28) имеют одинаковую левую часть в силу неиз |
||||
менности расхода потока |
по |
пути движения (qi-2 = q2-z = qz- 4 — |
||
= <7= const). Последовательное |
решение уравнений этой системы |
позволяет получить искомые значения ординат депресеионной кри вой Я 2и Я3.
Для приближенного решения уравнений (IV,28) можно предва
рительно определить величину |
расхода потока q |
по формуле |
(IV, 17) с подстановкой значений |
мощности и напора |
потока в его |
крайних граничных сечениях и использовать это значение при по следовательном решении системы уравнений (IV,28).
Пр и м е р . Грунтовый поток, движущийся в аллювиальных пес чаных отложениях с коэффициентом фильтрации k=\3,4 м/сут, дренируется рекой. В скважине, заложенной на урезе реки (скв. 2), отметка уровня грунтовых вод составляет 128,4 м, отметка водо упора 114,4 м. В скважине, заложенной на расстоянии 1000 м от уреза реки вверх по потоку (скв. 1), отметки уровня воды и водо-
упора соответственно равны 140,04 и 108,04 м (рис. 59). Питание грунтовых вод на участке скважин 1—2 отсутствует (Ф' = 0). Тре буется определить расход потока и положение его свободной по верхности в сечении х, расположенном на расстоянии 563,5 м от скважины /, если отметка водоупора в этом сечении равна 111,62 м.
Скб
Рис. 59. Расчетная схема к определению расхо да и ординат кривой депрессии грунтового по тока
Как видно из исходных данных, водоупорное ложе имеет обрат 108,04— 114,4
ный уклон, величина которого |
— 0,00636. |
КЮ0 |
|
Решение задачи проводим с использованием формулы |
(IV,22). |
Для облегчения решения примем, что ось абсцисс (ох) совпада ет с горизонтальной плоскостью сравнения и с направлением дви
жения подземных вод. Плоскость сравнения |
принимаем |
располо |
|||
женной на нулевой отметке, тогда отметки |
уровня |
грунтовых вод |
|||
в сечениях будут равны |
пьезометрическим |
напорам |
( # t = 140,04 и |
||
Н2 —128,4 м), а отметки |
водоупора — высоте |
положения |
водоупо |
ра над плоскостью сравнения [z\ = 108,04; z2 = 114,4 и zx= lll,6 2 м).
Мощность потока на его границах |
(в сечениях 1 и 2) равна: h\ — |
|
= //і— = 140,04—108,04 = 32 л; h2 = H2—z2 = 128,4— 114,4= 14 м. |
||
Единичный расход потока определяем по формуле (IV,22): |
||
|
hi~\-h2 |
# і — Н2 |
Я1 -2 |
k ----------X |
|
|
2 |
El—2 |
3 2 + 14 |
140,04 |
128,4 |
= 13,4 X —тг--- X |
3,587 м3/сут. |
|
|
1000 |
Отметку уровня воды в сечении х можно определить на основе приведенного выше выражения (IV,25), которое получается если составить и приравнять выражения для расхода потока на участках 1— X и X—2. Искомой величиной здесь является Нх, значение мощ ности потока Лз вводится в формулу в виде hx= Hx—zx:
u h + h i H i — H 2 |
h i + # x — zx H i — tfX |
||
k -----------------------= k ------------------------------- |
|||
2 |
LI-2 |
2 |
X |
Подставляем численные значения и решаем полученное уравне
ние относительно Нх: |
|
|
|
|
|
3,587 = ,3,4 X 31 |
± |
П |
х У 40'04 - |
||
|
|
|
2 |
563,5 |
|
3,587 X 2 X 563,5 |
(Нх — 79,62) X (140,04 — Ях) . |
||||
Гзд |
“ |
||||
|
|
|
Откуда: # х2—219,66 Нх+ 11451,67=0 и Я х= 134,55 ж.
С учетом отметки водоупора в сечении х мощность потока соот ветственно равна /іх = Я х—гх= 134,55—111,62 = 22,93 м.
УСТАНОВИВШЕЕСЯ ДВИЖЕНИЕ РАДИАЛЬНОГО ПОТОКА ПОДЗЕМНЫХ ВОД
Естественные потоки подземных вод, как уже отмеча лось, могут быть сведены к одномерным плоско-параллельным и к
двухмерным в плоскости или в |
|
||||||
разрезе. |
Во |
многих |
случаях |
|
|||
двухмерные в плане потоки рас |
|
||||||
сматривают |
как |
радиальные, |
|
||||
линии |
тока |
которых |
помеща |
|
|||
ются в вертикальных |
плоскос |
|
|||||
тях, расположенных под углом |
|
||||||
и сходящихся или расходящих |
|
||||||
ся |
по |
направлению |
движения |
|
|||
подземных вод. Типичным при |
|
||||||
мером сходящегося радиально |
|
||||||
го |
потока |
является |
движение |
|
|||
подземных вод к скважине или |
|
||||||
колодцу (см. рис. 39). В при |
|
||||||
родных |
условиях |
радиальные |
|
||||
потоки наблюдаются в излучи |
|
||||||
нах |
и |
на |
непрямолинейных |
|
|||
участках речных долин, где их |
|
||||||
воды дренируются руслами рек. |
|
||||||
При этом |
линии |
токов будут |
|
||||
сходиться |
или расходиться и |
|
|||||
поток будет иметь переменную |
|
||||||
ширину. В таких условиях от |
|
||||||
несение расхода потока на еди |
|
||||||
ницу его ширины уже не допу |
Рис. 60. Схема радиального потока |
||||||
стимо (так как В переменная). |
грунтовых вод: |
Вначале рассмотрим ради альный поток грунтовых вод,
а —разрез |
(кривая |
депрессии показана |
условно), |
б — план |
(радиальный расхо |
дящийся поток), в — план (радиальный сходящийся поток)
основанием для которого служит горизонтально залегающий водоупор (рис. 60). При этом примем две возможных схемы движения потока. Первая схема соответствует расширению потока в плане по пути его движения (рис. 60, б), т. е. отвечает схеме радиально-рас- ходящегося потока. Вторая схема соответствует сужению потока в плане и отвечает характеру радиально сходящегося потока (рис. 60, в ) . Изменение ширины потока происходит по линейному закону. Тогда, считая направление оси ох совпадающим с направлением движения и принимая первое ограничивающее поток сечение за исходное, ширину потока В в любом сечении на расстоянии % от исходного (сечение 1) можно определить по формуле:
В = B i В1-2Bi X, (IV,29)
где В{ и В2— ширина потока в ограничивающих его сечениях 1 и 2, расположенных на расстоянии LI_2 одно от другого.
Решение, как и ранее (см. гл. IV, стр. 114), может быть полу чено гидравлическим методом на основе интегрирования дифферен циального уравнения для расхода потока с той лишь разницей, что здесь вводится в рассмотрение переменная ширина потока В и уравнение записывается не для единичного, а для общего расхода потока. Размещение координатных осей и принятые обозначения показаны на рис. 60.
Общее выражение уравнения Дюпюи для расхода радиального потока в любом его сечении в дифференциальной форме имеет вид:
Q = - k h B d^ ~ . |
|
(IV,30) |
|
С учетом линейного изменения |
ширины потока |
В уравнение |
|
(IV,30) видоизменяется на: |
|
|
|
Bz — Ві |
, dh |
(IV,31) |
|
Q = |
|
х) h — . |
|
f |
U - Z |
dx |
|
После разделения переменных и интегрирования этого уравне ния (IV,31) в пределах от сечения 1 до сечения 2 получим выраже ние для определения общего расхода грунтового радиального по тока при горизонтальном залегании водоупорного ложа:
= к B z - B i |
f h - h t |
ІпВ2— ln ß i |
(IV,32) |
2LI_2 |
Формула для определения ординаты кривой депрессии в любом сечении, расположенном на расстоянии х от начального сечения 1, может быть получена на основе приравнивания выражений для расхода потока на участках 1—2 и 1 — х и имеет следующий вид:
,2 |
hl2 ~ В 1 ln В х— In В х |
h \ - h \ |
hi — |
/?х— В г ln В 2— In В х |
X, ( I V , 33) |
|
В\~2 |
где в х — ширина потока в сечении, отстоящем на расстоянии х от сечения 1, и определяемая по приведенной формуле (IV,29).
При наклонном залегании водоупорного ложа для определения расхода радиального потока может быть применена приближенная формула Каменского (IV,22), полученная для условий одномерного грунтового потока с наклонным водоупором. Для этого она должна быть переписана применительно к определению полного расхода потока с учетом переменной его ширины:
„ 7 ^hi-\-hz |
Bz— В1 |
Hi — Hz |
(IV,34^ |
|
2 |
l nß2— lnßi |
L i - z |
||
|
Эта формула (IV,34) применима для определения расхода ра диального потока с линейным изменением ширины потока В по пу ти его движения. При нечетко выраженном характере изменения ширины потока иногда используется следующая приближенная формула:
hiBi hzBz |
H i - H z = / і |
+ К 2 |
H i - H z |
Q = k ---------------- |
L\-z |
2 |
(IV,35) |
|
Bi-z |
В этой формуле используется среднее значение площади сече
ния потока, определяемое как среднеарифметическое из |
значений |
||
/ |
, |
Ві + |
Fz \ |
площади в ограничивающих поток сечениях I |
с ср = |
— ------ |
. |
Для определения ординат кривой депрессии радиального пото ка используется прием рассмотрения отдельных фрагментов с по следующим приравниванием расходов аналогично тому, как это делается для одномерного плоскопараллельного потока при на клонном залегании водоупорного ложа (см. гл. IV, стр. 122).
Решения для оценки условий фильтрации радиальных потоков к скважинам рассмотрены ниже (гл. IX и X).
ДВИЖЕНИЕ НАПОРНЫХ ВОД В ПЛАСТАХ ПОСТОЯННОЙ И ПЕРЕМЕННОЙ МОЩНОСТИ
Движение напорных вод в пластах постоянной мощности.
В условиях изолированного однородного пласта постоянной мощно сти напорный поток характеризуется постоянным по всем его сече ниям значением скорости фильтрации, т. е. имеет место равномер ное движение подземных вод.
Расчетные формулы для определения расхода подземных вод напорного потока и построения депрессионной кривой могут быть получены как на основе интегрирования дифференциального урав нения для единичного расхода (гидравлический метод), так и с по мощью интегрирования основного дифференциального уравнения фильтрации. Рассмотрим последовательно оба эти метода.
1. Выделим фрагмент напорного потока в пласте постоянной мощности (т = const), ограниченный двумя вертикальными сече