Добавил:
Опубликованный материал нарушает ваши авторские права? Сообщите нам.
Вуз: Предмет: Файл:

книги из ГПНТБ / Булыгин В.Я. Гидромеханика нефтяного пласта

.pdf
Скачиваний:
17
Добавлен:
24.10.2023
Размер:
11.6 Mб
Скачать

решения задачи при предположениях, сделанных для рассматрива­ емого случая. Получаем

я

K ,E n -= if

К Е Н=

ноЯнд

кк*0Н'

Я в/і~в0Я вд

(III.7.22)

 

я

Я

 

2 Я

2 Я

второе выражение — по аналогии с первым.

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

В этом случае мы

будем иметь К я — К в = К. Тогда

из (III.7.22)

следует:

 

 

 

 

 

 

Для средней вязкости получаем

 

 

 

 

 

М'нИ'В(£ң~ТЯВ)

 

(III.7.23)

 

 

HB* H

Рн-^в

 

 

 

 

 

Отношение средней скорости частиц воды к средней скорости

всей жидкости будет равно:

 

 

 

 

 

wB

Рн (Я н + Я в)

 

Е п + Е в

 

 

W l\B

И в^н + Ц н ^ в

 

р I Цв

р

 

 

 

 

 

Ев-Т— Ец

 

Последнее выражение можно представить

в виде:

 

 

 

 

 

 

 

(III.7.24)

где

 

а = Ин/Цв=1/Ц0-

 

 

 

 

 

 

Будем

считать,

что величина

ЕВ]ЕН мала.

Тогда,

разложив в

(III.7.24)

последнюю скобку в ряд, получим приближенно:

 

 

Wв

 

а Ев ).

(III.7.24)

или грубо

 

WHв

 

 

Ц>в

H_

 

 

 

 

 

 

 

 

^HD

I^B

 

 

Отношение дебитов воды и нефти в эксплуатационных скважинах

будет равно:

Яв

Е в

р.ң

 

 

 

 

 

(III.7.25)

 

 

Ян

Ен

Цв

 

 

 

 

 

100

Пример. Определить отношение скорости движения воды к сред­ ней скорости движения жидкости в разрезе пласта при следующих данных. Движение воды и нефти происходит послойно и можно при­ нять схему 3. Величины ßHи ßB считаются постоянными, пусть они

равны ß„ =

0,8, ßB=

0,7;

насыщенность изменяется в пределах 0 ^

=? s ^

1; я = 2,5. По формуле (III.7.24) получаем данные, сведенные

в табл. 1.

 

 

 

 

 

 

Т а б л и ц а 1

 

 

 

 

 

 

 

 

5Н

 

S B

 

W B

SH

 

S B

 

WB

ßH

ßB

W HB

ßH

ßB

V 'H B

1

0,8

0

0,7

2,5

0,3

 

0,7

 

1,27

0,9

0,1

2,22

0,2

0,8

1,15

0,8

0,2

1,95

0,1

0,9

 

1,07

 

 

 

 

 

 

0,7

0,3

1,78

0

 

1,0

1,00

0,6

0,4

1,61

 

 

0,01*

0,246 *

 

0,5

0,5

1,47

0,05*

0,232 *

0,4

0,6

1,37

0,10*

0,222*

 

' ’

 

 

 

 

 

* Данные получены по формуле (III. 7.24) при тех же данных для малых водонасы­ щенностей.

§ 8. Упрощенная схема «струй»

При расчетах вытеснения нефти водой из пластов иногда доста­ точно использовать простейшую модель фильтрации. Для этой цели можно, например, пласт и породу считать несжимаемыми и аппрок­ симировать кривые относительной проницаемости двумя прямыми. Согласно (1.2.8) это будет соответствовать значению а = 0, Ъ = 0. Тогда получаем

/св — ев == s св0, hu = ен = 1 s сн0.

(III.8.1)

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

^

= ~ і і г ѵ р ;

(іи .8.2)

^

=

(III.8.2*)

Заметим, что при постоянных значениях проницаемости в раз­ резе пласта (III.8.2) и (III.8.2') совпадают с (III.7.7) и (III.7.8), если за е,- принять среднюю по разрезу пласта вытесняемость. Тогда рассматриваемый случай фильтрации будет соответствовать схеме струй при постоянных значениях проницаемости в разрезе пласта (чтобы перейти к нему, достаточно заменить е,- на Е;).

101

Запишем уравнение материального баланса:

д

Щ

др \ .

д / К к * о р

= тН

дъі

(111.8.3)

дх

Р£

дх ) '

ду \ ре ду

 

 

dt

 

Подставляя (III.8.2)

в (III.8.3), получаем

 

 

 

 

~^mIlEiwu +

-^- mHBlWty = mH

дЕ(

(111.8.4)

 

 

 

 

 

dt

 

 

Суммируя почленно уравнения (III.8.4) при разных индексах і , получим

- ^ 7 7JlHewm + -Jjj-tnHeWaa = 0,

(III.8.5)

где средняя скорость воды н нефти равна:

wH„= -

w„

(III.8.6)

Суммируя уравнения (II 1.8.3) почленно, получим

Для уравнения (III.8.7) может быть сформулирована краевая задача, решая которую, получим значение давления р. Предположим, что оно известно. Тогда из (III.8.2) могут быть найдены средние скорости движения частиц нефти и воды, а используя (III.8.6), будем иметь среднюю скорость совместной фильтрации нефти и воды и со­ ставляющие скорости:

 

н:іл-'

К /_Ён_ I _ б в _ \ др_ .

 

ет

V Ин ‘

P D /

д х ’

 

 

д ' нв у

К

< е „

.

g n

\

др

(Ш .8.8)

ет

\ Р н

'

Рв

/

ду

Сравнивая с (III.8.2), можно определить скорость каждой из фаз,

например воды, в виде:

 

 

 

р

 

/ (^в) ^НВ)

(III.8.9)

fcB

 

 

f t p

 

P H ^ B

(III.8.10)

I( b>

рнев + р в8н

/ — функция Баклея — Леверетта для рассматриваемого случая. Подставив в (III.8.9), (III.8.4), получим:

d m H e f w Ha

 

d m H e f w HB

_ т ң дев

д х

'

ду

dt ’

откуда аналогично § 6 получаем

д] .

д}

5ев

+

ду

g f -

102

Учитывая, что

d f ( е в)

. _

д г ъ

д х

' Ев

дх

и обозначая

 

e2|X„LlB

 

 

wu- & U w m - w HB (адів+ 8ві11н)

получим окончательно

дев Зев

Ö Ë B .

w“* - t e + w“y ~ b = - a r

ИЛИ

(III.8.11)

(III.8.12)

ГО

— это уравнение Кельвина,

II

(III.8.12')

контролирующее движение фронта

жидких частиц ев = const, в котором wa — скорость движения

фронта. wu может быть представлено в виде:

"

е вр „ + 8 „ ц в

(III.8.13)

т

Сравнивая (III.8.2), (III.8.6) и (III.8.13), получаем

wuWuв — wnu\.

(Ш.8.14)

— это условие согласования

скоростей.

Учитывая (Ш .8.1), видим,

что для этого случая скорость линии равной вытесняемости, совпа­ дают со скоростью линий равной насыщенности.

Из (III.6.16) следует, что wu = —еве/|лв (ln f)'\7p-

Следует отметить, что упрощенная схема струй удобна для вы­

числения движения нефти и воды по

картам

изобар. В

этом слу­

чае по картам вычисляются dpjdn,

затем

определяется скорость

движения изогидры (линии равной

насыщенности) и

графически

строятся перемещения Апс линий равной насыщенности

 

Апе— wu At

 

(III.8.15)

за отрезок времени At. Для этой цели ввиду погрешностей в пара­ метрах задачи применять более сложную схему нецелесообразно. Примеры, рассчитанные с использованием схемы, приводятся в сле­ дующей главе.

§ 9. Вытеснение нефти водой из слоистых пластов

Слоистые неразобщенные пласты

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

103

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

Введем функцию усеченного нормального распределения прони­ цаемости. Если обозначить x Q— предел усечения (рис. 35), то со­ гласно работам [24, 25] функция распределения запишется

XX \

- Ф ( .

х0 — х \

 

_ ф ( а )

 

0

)

 

 

' Ж0

 

х \

 

 

1 —Ф (Ч

 

 

 

 

что

 

 

 

 

 

 

 

 

Ф' Г

Т

)

/

А.2= /(*) =

 

\

я

 

 

1 -Ф

 

 

 

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

г

М г —X-{- /о; М* = х-іг Іа (х0— х)-\-а2,

здесь

X

 

Ф : _J_

1 Т Т

Ге 2 dt; Ф

' ___ е 2 .

У2л

J

У2п

(III.9.1)

(III.9.Т)

Обозначим х х = і, тогда для значения х 0 будем иметь

a:0 = a:^l + i - |j .

Т а б л и ц а 2

SB

е

к н

Кв

SH

е

і

0,999

А

1 со

0,001

- 3

0,977 0,841

Я + 0,06 А + 0,29

А'—2,37 А - 1,52

0,023 0,159

—2 - 1

0,500

СО О0 о . О +

0,500

0

0,159

А+ 1,52

А0,29

0,841

1

0,023

А + 2,37

О 1

О

СО

0,977

2

0,001

А + 3,14

А

0,999

3

104

Рассмотрим проницаемость пла­

стов и обозначим X = К, а = /с„. Тогда из последней формулы по­ лучим

а:0= £ ( і + і-^ ).(Ш .9 .2 )

Величину і будем называть индексом обводненности. Если ис­ пользовать вероятное отклонение проницаемости кг, то будем иметь

ftp — х

а ВРіі

где рі — 0,6745.

Рис. 35. Усеченное нормальное рас­ пределение.

X — среднее значение; х0 — предел усечения.

В результате получаем

или

*0 = Я ( і + * ! - ! :) , (НІ.9.3)

где кг = рхка.

Рассчитывая проницаемость для различных усечений, сведем по­ лученные данные в табл. 2 и 3.

Средняя проницаемость для различных значений усечения (со­ отношений подвижной жидкости в пласте) будем записывать в сле­ дующем виде:

K ^ K ^ l + I ^ y , (III.9.4)

к ѣ= к ( \ - і ^ -

О

0,02

05

О

о

іЛ

см

о

ю

о

о-

о

т-< 05

о

со

05

о

о

о

тН

St«

 

о

со

+

 

о

 

 

 

ю

00

 

 

ю

о

0,98

 

со

о

со

+

1

 

 

 

со

st«

 

 

со

чгН

 

со

см

см

см

см

05

+

1

О

 

о

st«

 

 

со

ю

 

05

со

 

+

]

г—

 

о

 

 

Ь<

 

 

0505

тн

1

 

ю

о

1

о

 

 

st«

о

ю

 

со

05

 

о

1

см

1

+

о

Ssd

 

 

 

 

 

о

2,68-

о

см

со

 

05

 

СМ

 

 

+

1

о

1

 

і*:

 

 

00

ю

 

 

о

ю

см

со

о

со

о

+

1

о

1

 

st«

 

 

 

S?

о

sf«

 

 

1

і-Н

сз

со 1а

•У

 

Sjj-t

105

к и = к ( l + Л - ^ ) ;

£ в= я ( і - / 2-^-),

где значения индексов I х, 12, которые будем называть индексами проницаемости, могут быть взяты из таблиц или подсчитаны по фор­ мулам (III.9.1).

Расчеты движения нефти н воды в слоистых пластах следует вести по схемам «струй», при этом будем опираться на формулы (III.7.7)—(III.7 .20) и л и ' предыдущего параграфа. Средние абсолют­ ные проницаемости, определенные их индексом насыщенности і, іг ,

рассчитываются по формулам (III.9.4), которые и

являются

исход­

ными.

 

Ен = еІЮ, Е а =^епо.

 

 

 

 

 

 

 

Будем

иметь

получаем

 

 

 

 

 

Тогда

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

 

 

 

 

 

 

 

 

 

vt = - ^ E

V

p

 

 

 

(III.9.5)

н уравнение (III.7.9) перейдет в следующее:

 

 

 

 

д

KjSi0H

dp

д

КіЕі0Н

dp

■Nt

 

д І>

I ... и

d S

i

(III.9.6)

dx

іі£-

дх

ду

Pi

ду

 

dt

~rmH

at

где Ki определяются по формулам (III.9.4) и являются функциями

насыщенности. Таким образом, переменилась роль К £ и е/0: теперь Кі является функцией насыщенности, а е0 — постоянная.

При использовании упрощенной схемы «струй» вместо формулы

(III.8.1)

следует

взять

(III.9.4).

 

имеем

 

Тогда

вместо

уравнения

(III.8.3)

 

 

о

A'te,o

др_ I _д_

A',et0

dp

=

d s j

 

dx

щ

dx

' dy

щ

dy

 

dt

Аналог функции Баклея — Леверетта

будет:

 

 

/ ( в в ) = ~ >

Ці,* вБвѴ

..... ■

 

 

 

 

 

Р н А

В6 В0 -Г

Р ВА ЦЁНО

 

Видим, что /'(ев) — монотонная функция, т. е. фильтрация должна происходить без скачков насыщенности.

§10. Фильтрация в пачках одновременно вскрытых пластов

спроницаемыми перемычками

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

Юб

клиниваться и замещаться непроницаемыми породами. Будем счи­ тать, что в пластах находятся жидкости с одинаковыми физическими свойствами (разноцветные жидкости).

Ограничим рассматриваемую область поверхностью цилиндра, образующие которого перпендикулярны плоскости приведения давле­ ния с поверхностью продуктивной части Основания этой пло­ скости представляют верхнюю и нижнюю границы пачки продуктив­ ных пластов с проекцией F на плоскости приведения. Фильтрацион­ ный поток Q через боковую поверхность цилиндра F, и поверхности скважин F(j будет выражаться следующей формулой:

где п — число

скважин (/); I — число пластов і.

С учетом дебитов скважин

можно

записать:

I

I

п

п I

Этот поток должен компенсироваться изменением объема вследствие сжатия жидкости и горных пород:

где р* — упругоемкость t-того пласта; т, — объем продуктивной части пластов; t — время.

Сравнивая, получаем

П

(Ш.10.1)

В левой части равенства (III.10.1) каждое слагаемое сумм рас­ пространяется только на один пласт, что дает возможность для каж­ дого из пластов перейти от поверхностного интеграла к объемному. С учетом линейного закона фильтрации это дает следующее:

 

п

 

2 Л Ж - й - ) ѵ

4 r > < = - 2 , i

( I I U 0 -2)

1=1 X

 

 

И Л И

n

107

где

А-, (kxi,

kyi,

kzi) — тензор

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

пластов;

р,- — вяз­

кость жидкостей.

 

 

 

 

 

(/сЛ./р),- и

упругоем-

 

Введем средине значения гидропроводиостп

костп ß* пластов

формулами:

 

 

 

 

 

 

 

 

 

 

 

 

Н і

 

 

 

 

( і

г ) -

- Г І Т f

d z '

 

=

 

( I I U 0 '4)

 

 

 

 

О

 

 

о

 

 

где

Hi — мощность

і-того пласта.

 

 

 

 

Так как

принято

кгі]\хі dpt]dz

vzi, TO

 

 

 

 

11i

 

 

 

ih

 

 

 

 

 

0

 

=

0j

J 7 Vzi dz =

ѴІІ — VzV

(ШЛ0.5)

 

 

 

 

T . e. последний интеграл представляет потерю жидкости через по­ дошву и кровлю г'-того пласта.

Вводя (III.10.4) п (III.10.5) в (III.10.3), получаем

П 2 £[(-П*'£]+т[(П*‘т]-

F 4=1

п

)dxdy = 2 ?/•

J =I

Стягивая площадь F к некоторой малой площади, для которой можно составить надежное среднее плотности отбора N (х, у; t), будем иметь:

Ж [(тг)^] + і [(Н я- w

\

-

£=1

 

 

 

- m i H

i ^ r + v°zi-vzta} = N (x, у

t),

(III. 10.6)

где

 

 

 

J

fN{x, у, t) dF —2 Qj.

 

 

F J

3 = 1

 

 

В качестве закона фильтрации через перемычку примем

vh = — а,о (А —Рі-і).

 

(III.10.7)

Значение аі0 будем относить к подошве пласта, а аік — к его кровле.

108

Далее считаем, что значения дебитов скважин известны раздельио по всем пластам и обозначим их д( -. Тогда уравнение (III.10.6) рас­ падается на систему уравнений:

= Ф * ) і И і ~ 9І

у ; «) + «fo(P<P t - i )

-o*ft(P<+i—А);

(H I.10.8)

I

j ІѴ) (я, у ,

t ) d x d y =

Y i

q t j .

(III.10.9)

F

 

 

J =

1

 

Значения a l0 могут не совпадать со значениями a /_1 к, например,

при различных знаках разности

(р{рі_г).

 

 

§ 11. Схема напорной фильтрации в пачках одновременно вскрытых пластов

При выводе уравнений (III.10.8) для определения давления в пач­ ках одновременно вскрытых пластов с проницаемыми перемычками считалось, что переток жидкости из одного пласта в другой незна­ чителен, так что при этом вязкость фильтрующей смеси не изме­ няется. Теперь мы освободимся от этого условия [25]. Будем счи­ тать, что в каждом из пластов течение происходит по схеме струй. Примем, что для каждого из пластов проницаемость и эффективность вытеснения постоянны по мощности, это естественно предположить для тонких пластов, разделенных перемычками. Тогда на основании (III.7.9) для і-того пласта имеем:

д

kniffI

dpi

д

kniffі

dpi

дг

 

тт о*

dpi

I

тт

dsHj .

 

 

 

я ,, ....

 

-^ні

5нЬ“

іРн(

at

~ г Щ -О і

dsBiо, ,

дх

fi„

дх

ду

dt

 

 

ду

Цн

*,н‘

“»‘-

‘г-1

 

 

dt

д

kBjfft

dpi

д

квіНі

дрс

дг

IT ft*

dpi

тпіН і

(II I .ll.l)

dt

дх

(.1/

дх

 

 

Щ

ві- * віРві

dt

где ки1, кві — эффективные проницаемости для нефти и воды і-того пласта.

Будем нумеровать подошву пласта индексом і — 1, а кровлю г. Обозначим их проводимости, для нефти и воды соответственно:

ан і}{-1 —

Kt-г

КI

 

Пнн и

ѵі

где Щі і - і , . . ., К іі Hi- i , Hi, t — соответственно эффективные про­ ницаемости для нефти и воды и мощности подошвы и кровли.

Расход нефти и воды через единицу площади для подошвы будет выражаться следующим образом:

tf«H

@4,1, i - l ( P i

P i - l)l H n в = ÖB, i, i - 1 { P i

P i - 1)> (III.1 1 .2 )

109

Соседние файлы в папке книги из ГПНТБ