ВУЗ: Не указан
Категория: Не указан
Дисциплина: Не указана
Добавлен: 21.10.2024
Просмотров: 55
Скачиваний: 0
0,16 X / |
l--0,r |
|
|
|
0,1/ |
|
|
|
|
0,08 |
|
|
|
|
0.Ц |
|
|
|
|
0 |
OA |
|
|
|
0,Z |
1,2 |
1.0 |
1,6 1,6 X i |
|
- от • |
|
|
|
|
- 0,08 |
|
|
|
i f / |
- |
|
|
|
|
Ss
-Off
-
0,16
V*(X,'t)
Рис. 1. Графики значений функций N(x, i, n) и N*(x, t, n— 1).
13
Рис. 2. Графики значений функций М(х, t, п) и М*(х, t, п—1).
14
сей на дне канала; Но — наполнение канала; р — мутность воды; у — объемный вес скелета донных отло
жений; |
б — коэффициент водоотдачи |
(недостатка насы |
|||
щения); |
Ь\ — удаление ближайшей |
стороны |
канала |
от |
|
начала |
координат; |
62= b.i-j-2£>, / — ширина полосы, |
р = |
||
= +1 — параметр, |
учитывающий характер |
граничных |
|||
условий |
пласта (р = —1, когда обе |
границы |
проницае |
мы); erfu — функция ошибок, Ei(-«) — интегральный эк споненциал. Остальные обозначения приведены в табл .1.
Для того чтобы можно было производить прогнозы уровенного режима грунтовых вод в окрестностях ка нала, нами совместно с А. Е. Зайцевой с помощью ЭВМ «Минск-22» были протабулированы функции
N (х, /), N*(xyt), М(х, t) и М*(х, /), входящие в расчет ные зависимости (табл. 1). Графики этих функций при ведены на рис. 1 и 2.
Для иллюстрации методики расчетов рассмотрим следующий пример.
Пусть имеет место инфильтрация из канала при под держании постоянного наполнения в нем и следующих
исходных данных: |
м/сутки, |
&= 5 м/сутки, |
6= 0,15, |
|||||||
2& = 20 |
ж, |
&0 = 0,001 |
||||||||
hcр = 20 ж, |
а = |
kh.ср |
■15•20 |
=667 м3/сутки, # 0= 1 |
ж, |
|||||
б. |
= 0,15 |
|||||||||
ср |
’ |
|
|
q0 = 2 м/сутки. |
|
|||||
р = 0,225 кг/м3, |
у —1,45-103 кг/м3, |
|
||||||||
Найдем значения постоянных параметров |
|
|
||||||||
а = |
2р |
|
|
2-0,225 |
_ |
0,45 |
_0 311 |
|
- |
|
yk0H0 |
1,45-103 -1 - 10_3 ~ |
1,45 |
“ ’ |
ж2 |
’ |
|||||
|
||||||||||
р# 0 |
|
|
0,225-1- |
1 |
I 0,00Ы ,45Ю3\ |
|
||||
|
|
|
|
|
|
|
2-0,225 |
j |
|
|
|
|
|
|
|
|
1,45-10s |
|
|
||
|
|
|
|
= 6 - 104 ж; |
|
|
|
|
||
|
|
|
h |
6 - 104 |
|
0,36 |
сут |
|
|
|
|
|
|
k0Hо |
10-3- 1 |
|
ж2 |
|
|
15
Т а б л и ц а 1
Сводка формул для расчета инфильтрации из бассейнов при различных видах питания
Вид инфильт |
Режим |
Наименование вида |
|
|
рации |
работы |
инфильтрации |
|
|
канала |
|
|
|
|
Инфильтрация |
|
|
|
из канала в |
|
|
|
неограниченный |
|
|
|
пласт |
|
|
D |
|
|
|
Я |
|
|
|
Я |
|
|
|
0J |
|
|
|
Я |
|
|
|
К |
|
|
|
о |
|
|
|
я |
|
<7= |
Vat |
о |
|
|
я |
|
|
|
|
я |
|
|
|
к |
Инфильтрация |
|
|
о |
|
|
|
н |
из канала в |
|
|
о |
|
|
|
со |
пласт-полосу |
Расчетная зависимость |
Примечания и пояснения |
|
|
|
x-\-b |
х—b |
|
ь у |
erf — |
—erf — 7= |
|
||
а |
2 Vat |
2 Vat |
|||
х-\-Ь |
(x+b)‘ |
x—b |
(x—b)s |
||
7 = |
Ei |
|
|
p=-E1 |
4at |
2Van |
4at |
|
2 V aT |
||
s{xit)= |
—4 = 1 (1 + |
p) b V |
1 -+p [N(xa, t) |
||
|
о у aa { |
|
|
|
|
■N{Xl, t ) \ - [ N { Xi, t ) - N ( x 2,t)}.
N (x, f)=
К T 7C |
Я2 |
|
n= 1 |
1 |
sin (я—0,5)л x |
fV* (x, 0 = 71 ' |
X |
|
П=1 |
Xerf [(я |
0 15) то V 1 ]. |
2?
fk0H0
—дг—J— — х—
*i= —— .
х3- |
х-\-Ь2 - |
х—Ь2 |
|
|
д4= |
|
|
_ |
ЬО— |
— • |
|
6= |
|
, f=af//2 |
|
Инфильтрация |
s= —— ([—-^= -+v |
■(х+Ь)*— (x—b f X |
h |
|
JLl k |
||||||
|
из канала в |
3 / a " |
1L / |
е |
К " |
|
12а |
|
<?i= |
|||
|
неограниченный |
X erf V |
|
(х+6)3—(х—b)3 e—*t |
ir2 |
|
uk |
|||||
|
пласт |
|
l2aV~ |
V t )u |
|
|
|
|||||
|
|
|
|
|
|
|
|
|
||||
<v |
|
s= - ~ - f (1 + Р) Ь V t + p { M (Хз, 0 - |
|
|
|
|||||||
|
|
2oa 1. |
|
|
|
|
|
|
|
|
||
я |
|
|
|
|
|
|
|
|
|
|
|
|
я |
|
|
—М (хъ 7)] - |
[А4 (х4, 0 - |
М (**, <Т] }_ |
|
|
|
||||
я |
|
|
|
|
|
|||||||
Щ |
|
|
|
|
|
|
|
1г2 |
|
|
|
|
Он |
|
|
|
|
|
|
|
|
|
|
|
|
я |
Инфильтрация |
-----2 |
X I |
sin я т:х |
|
. , w |
|
|
i l l |
|
||
О |
из канала в |
|
*или |
|
|
|||||||
М (х, Л = |
— |
> |
•— ;-------- |
~ (к п |
|
|
||||||
|
пласт-полосу |
V |
' |
ТС |
|
Г,2/г2+Т |
|
|
|
a |
|
|
|
|
|
|
|
п= 1 |
|
|
|
|
|
|
|
|
|
|
— - |
2 |
sin (я—0,5) их |
X |
|
|
|
|||
|
|
|
Л4* (х, t)= |
— |
> . ■---- ----------------- |
|
|
|
||||
|
|
|
V ' |
т-. |
|
г .2 (я -0 ,5 )+ Т |
|
|
|
|
п=I
х е _ [Гь2(„_0,5)2+ Т] 7 .
Динамика формирования бугра грунтовых вод под каналом рассчитывается для рассматриваемого случая по зависимости
s= 12IV t |
erf |
2 at |
J |
|
( |
Х~ Ч |
У |
|
|
|
|
|
|
V 2 У at |
} _ |
|
|||
|
x + 10 Ei |
/ x + 10 у ! |
x — 10 |
|
|
||||
|
91,2 |
|
[ 2V~at |
) |
\ + |
91,2 |
|
|
|
|
|
X Ei |
x — 10 |
|
2T W |
|
|
|
|
|
|
2 V a t |
|
1,15 |
|
|
|||
|
|
|
|
|
|
||||
|
|
|
|
|
|
|
|
Т а б л и ц а 2 |
|
|
Формирование грунтового потока под |
|
|
||||||
|
|
инфильтрационным каналом |
|
|
|
||||
\ . x t м |
|
|
|
|
|
|
|
|
|
t, сут \ |
0 |
10 |
20 |
|
|
50 |
80 |
|
100 |
|
|
|
|
|
|
|
|
|
|
3,85 |
3,772 |
3,761 |
3,614 |
|
|
2,881 |
1,377 |
0,867 |
|
8,85 |
5,591 |
5,522 |
5,416 |
|
|
3,955 |
2,718 |
1,945 |
|
18,85 |
7,412 |
7,359 |
7,346 |
|
|
5,570 |
4,243 |
3,332 |
|
28,85 |
8,478 |
8,397 |
8,184 |
|
|
7,851 |
5,203 |
3,782 |
|
58,85 |
10,296 |
10,249 |
10,018 |
|
|
8,499 |
6,910 |
5,925 |
|
88,85 |
11,357 |
11,310 |
11,084 |
|
|
9,334 |
8,033 |
6,807 |
|
118,85 |
12,114 |
12,068 |
11,840 |
|
10,160 |
8,391 |
7,500 |
18
Результаты |
расчетов приведены в табл. 2. При этом |
|
для - |
|
/ х ± b у |
— ■' <Д,1 принималось, что—Ei |
||
. |
4 at |
[ 2У Н ) |
2,25 at |
|
|
а Д п |
---------- . |
_ _ |
|
(х ± bf |
Используя графики функций N (х, t) и М (х, t), ана логичные расчеты можно производить для всех случаев, рассмотренных в табл. 1.
Описанная методика прогноза уровенного режима грунтовых вод имеет значения при исследовании мелио ративного состояния орошаемых земель и расчетах искусственного регулирования запасов подземных вод.
В. И. ФОМЕНКО
ФИЛЬТРАЦИОННЫЙ РАСЧЕТ ДРЕНАЖНЫХ И ВОДОЗАБОРНЫХ СКВАЖИН
ПРИ ИНТЕНСИВНОМ ВОДООТБОРЕ
Основные зависимости динамики подземных вод по лучены в предположении, что движение воды в горных породах подчиняется линейному закону фильтрации (за кону Дарси). Данные натурных наблюдений [2, 8] и ла бораторных исследований в грунтовых лотках [4, 6] показывают, что при работе скважин возможно наруше ние линейного закона фильтрации.
Фильтрация с нарушением закона Дарси наиболее вероятна при заложении скважин в песчано-гравийных и гравийно-галечных отложениях речных долин, конусов выноса, древних ложбин стока, близком расположении естественных или искусственных источников питания, больших понижениях уровней воды. В рассматриваемых гидрогеологических условиях скважины могут эксплуа тироваться с высокими дебитами, при которых дополни тельные гидравлические сопротивления движению во ды в породе, обусловленные нарушением закона Дарси, становятся весьма существенными.
В результате экспериментальных исследований, про веденных различными авторами {1, 4, 5], установлено, что нарушение линейного закона фильтрации происходит при больших скоростях движения воды в горных поро
19
дах. Установлено также, что основной причиной нару шения закона Дарси является искривление струек воды в породе вследствие извилистости поровых каналов. Кроме этого, следует также учитывать возможность воз никновения локальной турбулентности в отдельных, наи более крупных порах.
В области нарушения закона Дарси фильтрация достаточно хорошо подчиняется двучленному закону сопротивления [5]:
|
|
= |
X |
, А |
У2 |
( 1) |
|
|
dr |
kt |
ki |
|
|
|
|
|
|
|||
где Н — напор; |
|
|
|
|
|
|
г |
— радиус; |
|
|
|
|
|
V — скорость фильтрации; |
|
|
||||
kt — коэффициент фильтрации; |
|
|||||
i |
— безразмерный |
параметр |
турбулентности; |
|||
— индекс |
( t= 1 |
относится к области пласта с |
||||
|
высокой |
водопроницаемостью, |
i = 2 относит |
|||
|
ся к пласту). |
|
решается в |
общей поста |
||
Рассматриваемая задача |
новке в предположении неоднородного скачкообразного уменьшения проницаемости пласта по радиусу. Частным случаем задачи является определение дебита скважины с гравийной обсыпкой.
Установившаяся фильтрация к скважине в такой по становке при сохранении закона Дарси рассматривалась М. Маскетом, К. Хилом, Л. Юреном [5], а при наруше нии закона Дарси задача решалась В. Н. Щелкачевым [8] и А. Ж. Муфтаховым [3]. Полученные зависимости отвечают частным схемам фильтрации воды к скважи нам, не доведены до практического использования, так как не указаны способы определения параметров, ха рактеризующих фильтрацию при нарушении закона Дарси.
Для выбора наиболее характерной расчетной схемы фильтрации воды к скважинам, отвечающей фактиче ским условиям их работы, проводились опыты в напор ном лотке, моделирующем секторный фрагмент области фильтрации. Как видно из рис. 1, 1, в зоне с высокой водопроницаемостью практически во всем интервале из менения дебита движение воды подчиняется двучленной зависимости (1). В породе (рис. 1, 2, 3) при малом де бите ,фильтрация подчиняется закону Дарси, а с увели-
20