Файл: Климентов П.П. Динамика подземных вод учеб. для геологоразведоч. техникумов.pdf

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

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

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

Добавлен: 10.04.2024

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

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

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

Для получения уравнения состояния сжимаемой жидкости свя­ жем изменение объема dVm с изменением плотности жидкости dp, учитывая при этом, что масса жидкости Мжпри изменении объема жидкости остается неизменной (закон сохранения массы). В со­ ответствии с этим изменение объема жидкости при колебаниях дав­ ления пропорционально изменению ее плотности и вместо выраже­ ния (11,50) можно пользоваться следующей формулой:

 

 

 

ßm

1 dp

 

 

 

(И,51)

 

 

 

p ~ 'd F '

 

 

 

 

 

 

 

 

 

Более строго переход от одной

формулы

(11,50)

к другой

(11,51) может быть обоснован

следующим образом. Учитывая, что

масса жидкости

д,

 

 

 

т/

=

Мж

 

МЖ=Ѵ жр и соответственно

---- , заменим в

 

 

 

Гж и dVw.,

 

 

Р

 

формуле (11,50)

значения

выразив

их через

значения

Мт и р:

Мт

 

 

 

 

 

 

 

VЖ

> а

dVm =

 

 

 

. ( 1 1 , 5 2 ) ,

Р

 

 

 

тогда:

 

 

 

 

 

 

 

 

 

 

dp

 

 

 

 

 

 

 

 

 

 

 

 

 

р

■Мя

 

dp

 

 

 

ß>K--

 

 

 

 

(H,53)

 

Мж

dP

~dP '

 

 

 

 

 

Это выражение (11,53) идентично приведенному выше выраже­ нию (11,51) и представляет собой уравнение состояния сжимаемой жидкости при упруговодонапорном режиме ее фильтрации в диф­ ференциальной форме. Интегрированием уравнения (11,51), с уче­ том изменения давления в пределах от Ро до Р и плотности — от ро до р, найдем уравнение состояния в конечном виде:

р = р 0 - е М р - р о),

( Н , 5 1 )

которое после разложения в ряд и некоторого упрощения может быть записано в следующем приближенном виде:

р =

Po +

Poßffi (Р — Ро) ■

(И,55)

В уравнении (11,55) р0 — плотность жидкости

при атмосферном

давлении (для воды ро =

1 2

\

 

1

)■

 

Аналогичным образом при исследовании фильтрации в услови­ ях упруговодонапорного режима учитывается и объемная упру­ гость пласта (пористой среды). При этом также принимается до­ пущение, что водоносный пласт ведет себя как упругое тело, под­ чиняясь закону Гука (необратимые процессы деформации горных пород не учитываются). Изменение объема пласта под влиянием гидростатического давления на скелет породы связывают с изме­


нением объема порового пространства. Характер изменения объе­ ма порового пространства породы в процессе фильтрации можно уяснить из следующего рассмотрения.

Давление на скелет породы, слагающей водоносный пласт РСКг определяется как разность между горным давлением Ргор (давле­ ние на водоносный пласт, которое создается весом вышележащих горных пород: РГоР=іЛуг.п) и пластовым гидростатическим давле­ нием жидкости Рпл:

Рек — Prop Рпл.

(11,56)

Из формулы (11,56) видно, что при уменьшении пластового дав­ ления Рпл давление на скелет породы Р ск возрастает и объем по­ рового пространства должен уменьшаться. Изменение объема по­ рового пространства упругой пористой среды происходит в соот­ ветствии с законом Гука:

dVпор

 

ßnn dPск-

 

 

(И,57)

~ к Г

 

 

 

 

 

 

 

 

 

Знак минус в формуле (11,57)

показывает,

что объем

 

порового

пространства увеличивается с уменьшением

давления

на скелет

породы Р ск.

 

 

 

 

 

 

Способность пласта уменьшать или увеличивать объем пор при

изменении давления на скелет

породы

характеризуется

коэффи­

циентом объемной упругости пласта рпл, размерность

и

физиче­

ский смысл которого видны из формулы

(11,57), откуда:

 

 

Рпл

1

dVпор

 

 

(П,58)

dPск

 

 

 

 

 

 

Как и коэффициент объемной упругости жидкости, коэффи­ циент объемной упругости пласта рпл показывает, на какую часть первоначального объема пласта Ко изменяется объем порового пространства при увеличении или уменьшении давления на едини­

цу. Размерность коэффициента

рПл зависит от размерности давле-

1

1

ния и выражается в --д

- или “

• Экспериментально установлено,

что значение коэффициента объемной упругости горных пород ßnjl

колеблется в пределах (0,3 2,0) X Ю~6 — .

м

Для решения практических задач выражение (11,58) удобнее переписать, используя значение пластового давления Рпл. Значе­ ние пластового давления Р пл находим из соотношения (11,56), от­ куда следует, что dPCK= d(Prov—Рпл). Принимая горное давление Prop величиной неизменной (что практически справедливо для мно­ гих задач), найдем, что dPCK= —dPnJI и соответственно выражение (11,58) изменится на:

1 dVпор


Формула (11,59) характеризует уравнение состояния пористой

среды при упругом режиме фильтрации.

пространства

по

Учитывая, что изменение объема пористого

 

dVпор

это

отношению к первоначальному объему горных пород-— —

 

уо

 

dVпор

d(nV о)

\

и есть изменение пористости dn (так как——— =

——------= dn],

Ко

'

выражение (11,59) может быть переписано в отношении dn следую­ щим образом:

dn = ßmi dPпл*

(И,60)

Уравнение (11,60) также является уравнением состояния упру­

гой пористой среды.

подземный поток дви­

Уравнения неразрывности потока. Если

жется без образования в нем пустот и разрыва сплошности, то он

 

подчиняется уравнению неразрывно­

 

сти, которое математически выража­

 

ет закон сохранения массы движу­

 

щейся воды. Уравнение неразрывно­

 

сти может быть получено на основе

 

рассмотрения

баланса

массы

жид­

 

кости в пределах

постоянного

эле­

 

ментарного

объема,

выделенного

 

внутри пористой среды. Для приме­

 

ра рассмотрим вывод уравнения не­

 

разрывности потока в условиях не­

 

установившейся

фильтрации

упру­

 

гой жидкости в упругой пористой

Рис. 18. Схема к выводу уравне­

среде. С этой целью выделим

в по­

ния неразрывности

ристой среде

элементарный парал­

 

лелепипед с ребрами dx, dy и dz, па­

раллельными координатным осям (рис. 18). Объем выделенного

элемента пористой среды dV является

постоянным и определяется

выражением (11,61):

const.

(И,61)

dV — dx-dy-dz =

Объем порового пространства выделенного элемента с учетом

пористости п составляет ДУдор:

 

 

dVnov п dV = п dx-dy-dz.

(11,62)

Масса жидкости, заполняющая поровое пространство в объеме

выделенного элемента М с учетом ее плотности р, равна:

 

М — р (/Упор = p-tidV = pndx-dy-dz.

(И,63)

Рассмотрим изменение массы жидкости в объеме выделенного параллелепипеда за промежуток времени dt за счет поступления ее через грани параллелепипеда. В направлении оси X через грань ЪЬ'с'с за промежуток времени dt в выделенный элемент поступает


масса жидкости, равная:

 

 

 

Mi = pvx dy-dz-dt,

(П,64)

где

ѵх — составляющая скорости

фильтрации

по оси X (рѵх —

массовая скорость фильтрации

по оси X); dy

dz — площадь гра­

ни

bb'c'c.

 

 

 

В то же время через противоположную грань параллелепипеда

dd'a'a, отстоящую от первой на расстоянии dx, вытекает масса во­ ды М2 , равная:

М2 — pux dy-dz-dt -f ^ ^ V*- d x d y dz dt.

(11,65)

dx

 

Накопленная в объеме параллелепипеда масса воды за время dt за счет ее фильтрации в направлении оси х соответственно со­ ставляет:

дМх dt — Mi — М2

d(pvx) dx-dy-dz-dt

 

dt

dx

 

 

d(pvx) dVdt.

(И,66)

 

dx

 

Аналогичным образом накопление массы воды в объеме парал­ лелепипеда за время dt за счет ее фильтрации в направлении осей у и z соответственно составляет:

дМу

d(pay) dV dt и

(11,67)

dt

dt

ày

 

dMz

d(pvz)

(П,68)

dt

dV dt.

àz

âz

 

Общее накопление массы воды в объеме выделенного элемента пористой среды за время dt определяется суммированием выраже­ ний (11,66) — (11,68):

дМ

âMx

dMу

дМг

 

dt

dx

dy

dt =

 

dz

 

= _

[ d(pPx)

a(poy)

a (po,) 1 ■dV dt.

ЛП,69)

 

*- dx

dy

dz J

 

Так как объем выделенного параллелепипеда во времени не изменяется (условие 11,61), то изменение массы воды в объеме выделенного элемента может произойти только за счет изменения плотности воды р и объема пористого пространства. Для определе­ ния изменения массы воды в объеме параллелепипеда за счет из­ менения плотности и пористости необходимо продифференциро­ вать выражение (11,63):

дМ

d(flp) dV dt.

(11,70)

dt =

dt

dt

 

Приравняв правые части уравнений (11,69 и 11,70) и сократив

их на dVdt, получим уравнение неразрывности,

характеризующее


неустановившуюся фильтрацию упругой жидкости в упругой по­ ристой среде:

д(рах)

д(рру)

д{9ѵг) =

д (пр)

 

дх

ду +

дг

dt '

'

При жестком режиме фильтрации, когда упругие свойства фильтрующейся жидкости и пористой среды не учитываются, т. е. плотность воды р и пористость пород п считаются в каждой дан­ ной точке пористой среды не зависящими от времени, уравнение неразрывности пространственного потока имеет вид:

djpvx) а(риу) d(poz)

(11,72)

дх "

ду

dz

При отсутствии той или иной составляющей скорости фильтра­ ции частная производная от массовой скорости по соответствую­ щей оси в уравнение неразрывности потока не входит. Так, напри­ мер, при установившейся фильтрации двухмерного потока в пло­ скости ху уравнение неразрывности принимает вид:

д(р^х) , д(рцу) =

дх ду

ОСНОВНЫЕ ДИФФЕРЕНЦИАЛЬНЫЕ УРАВНЕНИЯ ФИЛЬТРАЦИИ ПОДЗЕМНЫХ ВОД

Для получения основных дифференциальных уравнений фильтрации в динамике подземных вод используются два метода: а) метод синтеза трех видов уравнений, определяющих условия фіильтрации — уравнений движения подземных вод, уравнения не­ разрывности потока и уравнений состояния жидкости и пористой среды; б) балансовый метод, в основу которого положено изучение изменения водного баланса элемента потока подземных вод под влиянием факторов, предопределяющих условия их движения.

Для примера рассмотрим вывод дифференциальных уравнений фильтрации методом синтеза трех уравнений (на примере жестко­ го режима фильтрации подземных вод) и балансовым методом (на примере неустановившейся фильтрации грунтовых вод).

При жестком режиме фильтрации подземных вод уравнения движения пространственного потока при соблюдении закона Дар­

си могут быть записаны согласно формуле

(11,42):

dH

dH

- d H

Vx kx

Vy ky —-

;

vz k%— . (11,74)

dx

dy

 

dz

Фильтрующаяся жидкость и

пористая среда

принимаются не­

сжимаемыми, т. е. уравнения состояния имеют вид:

 

р = const,

п = const.

(И,75)