Файл: Вопросы водного хозяйства [сборник]..pdf

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

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

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

Добавлен: 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-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