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

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

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

Интегрируя уравнение (III.(5.22') методом характеристик, полунаем

 

d.v

dt

ds

 

 

 

~ T

Т Г '

 

 

Система его первых интегралов

 

 

 

S =

CjJ

X = Wut +

Со.

 

Пусть при t — О, X =

X (s,

0).

Тогда

с.г — х (s, 0)

и решение

задачи будет иметь вид

 

 

 

 

 

x = x{s,

0) + wJ.

(III.6.24)

Для анализа фильтрации заметим, что в формуле (III.6.23) ѵ, т

постоянны, следовательно, фильтрация зависпт только от вида кри­ вых относительной проницаемости. Примерные кривые / (s), / ' (s) представлены на рис. 30. Из рисунка видно, что кривая / ' (s) воз­ растает не монотонно, а существует максимум функции /(, следова­

 

тельно, рассматриваемый

случай от­

 

носится к схеме фильтрации с обра­

 

зованием фронта насыщенности.

 

 

Для установления положения фронта

 

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

воспользуемся

уравне­

 

нием

(III.6.21).

Считая поток

единич­

 

ной ширины, получаем для суммар­

 

ного

количества воды

 

 

 

 

Qt {t) = mH j

(s — cB0)dx,

 

Ппс. 30. Построение

фронта

 

о

 

 

 

 

 

фронта насыщен­

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

где дгф — координаты

 

ности.

 

 

 

 

Интегрируя по частям, будем иметь:

 

 

 

 

 

 

>

xds,

 

 

 

Qi = тНхф (5ф — сво)— mH J

 

 

 

 

 

1~сво

 

 

 

где Хф — нефтенасыщенность на фронте; св0 — остаточная

водонасы­

щенность.

0 Q (0) = 0,

Хф = Хф (s, 0), тогда и з -(III.6.24)

Пусть при t =

«следует

 

 

 

 

 

 

 

wtt ds,

 

 

Qt = mHxip{si(l— cls0) — mH J

 

 

 

 

1_св0

 

 

 

откуда, используя (III.6.23), получаем:

 

 

 

 

Qi = mHXф{Sф — cBo)Hvt f

df.

 

 

 

 

 

l “CBO

 

 

 

*90

Однако H v t =

Q i есть суммарное количество жидкости, протека-

емое через сечение единичной ширины.

 

Учитывая, что / (я) == 1 при я =

1 — св0, получим

 

Qi = тНхф (яф — сво)— Qt [/ (яф)— 1],

 

т е.

т _

/ W

Qt

 

 

(III.6.25)

 

Ф

(*ф~св0)тН

 

 

Из (III.6.23) и

(III.6.24)

при

условии х (я, 0) = 0

следует, что

 

 

_ ѵі

di

(III.6.26)

 

 

т

ds

 

 

 

Сравнивая (III.6.25) и (III.6.26), получим:

 

 

 

df _ _

/Сф)

 

откуда

v Ht

ds

Яф

Сво Qt,

 

/, = / (Яф)

 

 

 

 

(III.6.27)

 

 

С вО

 

 

 

Насыщенность на фронте определяется как значение, при котором

касательная к кривой /' (я) проходит через точку я =

св0, / (я) = 0'

Распределение

насыщен-

 

s

 

 

ностей приводится

на рис. 31.

 

 

 

 

Используя формулу(ІІІ.6.27),

Рпс. 31. Построение фронта на-

Рис. 32.

Зависимость величины скачка

сыщенностн.

насыщенности от отношения вязкостей ш

 

вогнутости

кривых Баклея — Леверетта

построим график значений Яф в зависимости от величины р,н/р в, фильтрующихся жидкостей и вогнутости кривых относительных про­ ницаемостей. Относительные проницаемости записывали, используя: (1.2.8) в виде:

&B = (S — св)[1 — а ( і — я)];

А£ = (1 — s —сн)(1 — as),

где а — параметр, характеризующий выпуклость кривых.

91

Рассчитанные кривые приведены на рис. 32. Из рисунков видно, что величина фронта насыщенности падает при увеличении отно­ шения р.н/|.ів и значения параметра а, т. е. вогнутости кривых.

§ 7. Схема «струй»

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

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

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

щенность равномерно распределена по

разрезу пласта. Тогда к*,

е,- — постоянные по разрезу пласта.

 

Средние скорости фильтрации и движения частиц жидкости

определятся соотношениями (III.5.7)

и (III.5.8), откуда имеем:

і

(III.7.1)

 

(III.7.2)

Будем различать следующие случаи фильтрации по схеме струй. Предположим, что строение пластов таково, что: 1) вода и нефть движутся по пластам с лучшей и худшей проницаемостью с равной вероятностью; 2) вода (вытесняющая жидкость), попав в пропластки

случшей проницаемостью, имеет тенденцию двигаться по ним. Первому положению отвечает строение пластов, когда пропластки

схорошей проницаемостью чередуются с пропластками с худшей проницаемостью так, что жидкость при движении встречает пласты

схорошей и плохой проницаемостью с равной вероятностью. Второму положению 1 отвечают пласты с выдержанными харак­

теристиками по их мощности.

1. Перейдем к рассмотрению первого случая [22]. Получим усло вия, при которых реализуется схема струй.

Предположим, что вытеснение нефти водой происходит таким образом, что вода целиком замещает подвижную нефть в отдельных

1 Этот случай будет рассмотрен в § 9 данной главы.

<92

пропластках. Осредняя скорость фильтрации нефти и воды по мощ­ ности пласта, получим:

где Hi д — часть мощности пласта, занятая подвижной нефтью или водой; кі0 — относительные проницаемости при полном замещении подвижной нефти или воды.

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

н„

 

 

Введем следующие обозначения:

 

 

Ка==~ і Ы к (Z) dZ'

К в= я Ь \ * (z) dz■

(HI.7.4)

о

о

 

В ычислив значения средних вытесняемостей, имеем:

Вторые интегралы,

стоящие в

скобках,

равны нулю,

так

как

в первом случае для

этих слоев

ен = 0, а

во втором ев =

0.

При

сделанных предположениях относительно схемы вытеснения е,- должны быть постоянными по разрезу пласта; обозначая их еі0,

получим

 

 

Е =

-

(III.7.5)

Для средней скорости фильтрации і-той жидкости получаем:

..

>

 

 

ѵ‘= — і й - ѵ р )

ѴР,

(Ш.7.6)

93

где

о

Однако имеем к„ = ено; А:* = еио. Тогда получим

V: =

КіЕі

(III.7.7)

 

w -

Определив составляющие средних скоростей движения частиц, жидкости W;, буДем иметь

 

 

нід

 

»' = т Ы

“" * = -

Ж J r S r V P * = - - é r V P .

(ПІ.7.8)

На основании

(III.7.7)

уравнения (III.5.2) могут быть

записаны

в следующем виде:

 

 

 

д

іЕіН dp

д

KjEjH

N t

— scHB*

дх

р,- дх

*7

Иі

ду

t

i t

Суммируя (III.7.9)

при разных индексах,

dp

I

m IT dsj

(III.7.9).

dl

+

m-“ dt

получаем:

д

КЕН

др

.

д

КЕН

dp ^

Ң р * (>Р

(111.7.10)

дх

рнп

дх

~т~

ду

р„в

ду

dt

 

где

 

 

N„ + N B= N;

s„ß* +

sBßS = S*;

(111.7.11)

Е Н+ Е В= Е;

 

 

 

 

h .,=

PnA \\Eи -j- \ittK BEB

 

(III.7.12)

 

 

 

 

 

 

 

При расчетах функции давления одно из уравнений (III.7.9) целесообразно заменить на (III.7.10), для которого формулируется краевая задача.

2. Распределение подвижной жидкости в нефтяном пласте

Рассмотрим вопрос об аппроксимации функций

вытесняемости

в нефтяном пласте. Предположим, что переходный

слой занимает

часть пласта и имеет мощность Н. В переходном слое

насыщенность

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

z = 0

с но>

(III.7.13)

1

SB =

SH= 1 Сн0,

=

1

^воi

z = H

1

^во) ” ^во*

SB =

94

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

= /н (2)і SB= /B(z).

Эти функции должны быть выбраны таким образом, чтобы выпол­ нялись краевые условия (III.7.13), т. е. должно быть:

/н (0) =

сно; /„

( В )

•=- 1 — сп

/в(0) =

1 — сн0;

(Н)

(III.7.14)

= св

Чтобы заданные кривые определяли распределение насыщен­ ности единственным образом, в качестве кривых зависимостей насы-

Рис. 33. Распределешіе подвижной жпдкостп в пласте.

Слой: I — переходный; II — совместной фильтрации.

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

Н

Н

II

Е я= -ff j s„ßHdz;

E B= -jj-j sBßBdz; E =

j (sHßH+ sBßB) dz.

0

0

0

Если примем, что количество связанного компонента остается постоянным, т. е. s„ßH= сн0; sBßB = sB — св0^ (первое предположе­ ние), то получаем:

 

н

Ен — ң

J [/н (z) сно] dz,

 

о

 

н

Е ^ =

(/в (z) — сво] dz.

 

о

95

Если же принять, что связь между эффективностью вытеснения нефти и воды и содержанием вытесняемого агента задается в виде

ßn «и ' ^но» ßi3 SB

то получаем

Ii

Е н ~ ~JT § f" (z) t/ч (z)— CHO1dz;

0

H

E B= -jT$ U(z) [/в (z) — cB0] dz. 0

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

Зададим зависимость насыщения в виде линейных функций мощ­ ности пласта:

S H — Q.]Z -f*

S B z==z Cl2Z - | - Ъ2.

Краевые условия (III.7.13) дают

при z = О

6і сІІ0, Ъп:=1 ен0,

при z — Н

1 св0 ■ ахН + сн0, св сіъН- 1 ^но*

Из последних условий получаем-

а, =

1 СНО С»0 .

CLty

1 — <н

 

Н

 

н

Подставляя аг, Ьг, а.2, Ь.г в исходные уравнения, имеем:

s„= — -C|^ ~ CBOZ + C„0; « ,= - 1~ с" ° ~ с»°г + і - е Ію. (III.7.15)

В первом предположении получаем для нефти

ня

Е а = ~ ^

sHßHdz = -|г j

1 C„ o CBO Z dZ.

о

о

H

 

Интегрируя последнее выражение, имеем

 

_

1

СНО

^вО

(III.7.16)

 

 

 

 

 

Для воды получим

 

 

 

 

 

я

я

 

 

 

 

Е в = ± - 1 * вРв = ±

^ - (

1-

с™ - Св° 2+ 1 - с но- с во)с12=:2

 

 

 

я

 

 

 

 

1 — 5но

с ВО

(III.7.17)

 

 

 

2

 

 

 

 

 

 

96

Для величины Е получим

Е = Е а + Е л= 1 - с т - с „ .

(III.7.18)

Таким образом, при первом предположении вытесняемость нефти и воды в переходном слое будет одинаковой и равна половине вытесняемости для нефтяной и водной зон.

Решим ту же задачу, принимая второе предположение о зависи­ мости эффективности вытеснения компонента от насыщенности, предполагая, что распределение насыщенности от мощности пере­ ходного слоя удовлетворяет соотношениям (III.7.18):

ни

 

о

 

о

 

 

Подставляя sHи интегрируя, получаем

 

 

н •

н

 

 

 

(2aH + Зсно).

-ff- ^ (az + сно) azdz =-jt-§ (a2z2+ сно az) dz =

о’

о

 

 

 

 

Подставляя значения а, имеем:

 

 

 

 

( 1 ---СНО----Сво) ( 2 - р СНО---- 2Сво)

(III.7.19)

 

Еп

6

 

 

Аналогично для воды можно написать:

 

 

 

я

 

 

 

 

Е в=

-JJ- ^ ( oz -J-1

сн0) (

яz —}—1

[- сно

св0) dz.

 

о

 

 

 

 

Интегрируя,

получаем

 

 

 

 

Еа =

d (2 — 2Сно — ^во) Н "Ь (1

^но) (1

^НО ^во)'

Подставляя значение я, будем иметь

 

 

 

( 1 — с н — с в ) (2 — 2 с н о - ) - с в о )

(ПІ.7.20)

 

 

6

 

 

и для величины Е получаем

 

 

 

g HO

С В о ) (4

^H O СВо)

 

Ң

J j J ____ ('l

 

При втором предположении вытесняемость нефти и воды не будет одинаковой. Так как связанной с породой воды обычно больше, чем нефти, то количество подвижной нефти в сечении переходного слоя будет больше, чем воды, что следует из (III.7.19) и (III.7.20).

Перейдем к рассмотрению разреза нефтяного пласта. Рассмотрим два случая.

Случай 1. Предположим, что нефть и вода фильтруются при на­ личии переходного слоя от нефти к воде (см. рис. 33). Будем считать,

7 Заказ 322

97

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

ции нефти в нефтяной ловушке. Примем, что

здесь sH= 1 — сво.

В переходном слое нефтенасыщенность изменяется от величины

5н = сно иа нижней границе переходного слоя

до величины s„ —

= 1 — сВо на верхней границе. Будем рассматривать слой совмест­ ной фильтрации воды и нефти за нижнюю границу. Примем поверх­ ность, ниже которой фильтрация отсут­

 

 

ствует, т. е.

ниже

этого слоя, ß„ = 0.

 

 

За верхнюю границу слоя совместной

 

 

фильтрации

будем считать поверхность,

 

 

выше

которой отсутствует фильтрация

 

 

воды,

т. е.

ßB =

0.

нефтенасыщенпость

 

 

Таким образом,

 

 

изменяется

(снизу

вверх) (s„= сно) от

Pnc.

34. Распределение мощ­

подошвы пласта

до

нижней границы

ности

пласта прп вычислении

переходного

слоя,

изменяется в пере­

подвижной нефти и воды.

ходном слое

и

снова

постоянна (sH —

 

 

= 1 — сво) на верхней границе переход­

ного слоя. Водоиасыщенность постоянна

(sB= 1 — сно) от подошвы

пласта до нижней границы

переходного

слоя,

изменяется в пере­

ходном слое и снова постоянна (sB = сво) на верхней границе пере­ ходного слоя. На нижней границе слоя совместной фильтрации неф­ тенасыщенность пусть будет равна s„ = сн0 и на верхней s„ = 1 —

— св0, а водонасыщенность соответственно sB = 1 — сн0 и sB — св0. Будем считать, что вне слоя совместной фильтрации эффектив­ ность вытеснения нефти и воды постоянна (ßB= ßB0; ß„ = 0) для нижней части пласта и (ßn = ßH0; ß„ = 0) для верхней. Тогда, внося

обозначения мощности:

Н ' — слоя

совместной фильтрации;

И в

слоя

фильтрации воды,

Н„ — слоя фильтрации

нефти,

имеем

для

зон

I, II, I I I (рис. 34):

I — е =

sBßB; I I — е =

s„ß„;

I I I

е =

®нРн ~Ь SBßB.

Определим вытесняемость нефти, воды и смеси нефть — вода при их совместной фильтрации. Считая, что проницаемость постоянна, имеем:

Н

Н н

Нд

Н'

E = -Jf \

&dz=z-Jf \<Eadz + - ~ j

eBdz + - ^ j e dz;

0

0

0

0

E = 7Г Ь

dz+ I r ] s»dz+ Ж ] &dz-

0

0

0

98

Обозначая средине значения насыщенности

 

и.

SB dz,

н

получаем

 

ң — J sHdz,

 

 

 

 

н

Е:

5в Р в Я в

I ® Н рнЯ „

I

1

я

я

+4-1edz-

 

Введя среднее значение вытесняемости в переходном слое, полу­ чим

О

Тогда можем записать:

TP _

SB ß ß ^B

I 5n ß lI# H I Н п

Ѣ '

- ~ Н

I—

Среднее значение вытесняемости ё было определено в предыдущем пункте в двух предположениях относительно связи между sH и ßH в переходном слое. Из последней формулы имеем:

ZP

внРцЯц

Я

-

Tj,

5ВРВЯВ

!

Я

- _

“ н

ң

> ң

Ьн’

-“ в—

н

Н

Ьв’

 

 

Е„

Е в= Е.

 

 

 

 

Если можно

считать,

что

вне

зоны

совместной

фильтрации sH,

sB постоянны по мощности пласта и равны соответственно sH= 1 —

— своі sB= 1 — сно, принимая предположение

первое, получим:

— 1

*"но

‘'воі

SBß B•

1

 

Сно

Св

и далее

 

 

 

 

 

 

 

 

 

 

(1

сно

сво)

( и

1

Я

\

 

 

--------я --------Iя “+

Т

)

Е „=■

( 1

с Н О

с В о )

(г

г

I

Н \

 

 

--------я -------- \ н * +

^

г ) ’

 

Е

1

сно

св0

 

 

 

(III .7.21)

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

чения множителей КНЕН, КВЕВ, которые необходимо знать для

99

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