![](/user_photo/_userpic.png)
книги из ГПНТБ / Булыгин В.Я. Гидромеханика нефтяного пласта
.pdfИнтегрируя уравнение (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 = |
||||||
«следует |
|
5Ф |
|
|
|
|
|
|
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) |
|||
|
|
5ф |
С вО |
||
|
|
|
|||
Насыщенность на фронте определяется как значение, при котором |
|||||
касательная к кривой /' (я) проходит через точку я = |
св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