Файл: Справочное руководство по проектированию разработки и эксплуатации нефтяных месторождений. Проектирование разработки.pdf

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

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

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

Добавлен: 29.02.2024

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

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

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

 

Обратимся

 

к

численному

 

примеру. В численных расчетах

 

приняты

следующие

значения

 

исходных

данных.

Отношение

 

вязкостей воды и нефти |ы0 = 0,4;

 

содержание

погребенной

воды

 

sn =

0,135.

Функции

относи­

 

тельных фазовых

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

 

тей, полученные

в лаборатории

 

подземной

 

 

гидродинамики

 

ВНИИ,

взяты

из работы [26].

 

По

ним

рассчитаны

функции

Рис. XII.9. Изменение дебита жидкости в за-

/ (s)

и

Ф (s)

при

jul0 =

0,4 и

висимости от времени

найдены

значения

постоянных

 

параметров

smax =

0,9;

S(j>= 0,702; Фф=1,54; snp = 0,784.

 

Приближенный расчет дебита

Q проводился при постоянном перепаде давле­

ния (Др =

1) по формуле (XI1.54) до прорыва воды и по формуле (XI1.56) после ее

прорыва. В последнем случае предварительно было рассчитано изменение во времени средней водонасыщенности по формуле (XI1.59).

Результаты выполненных расчетов по приближенным формулам (XII.54), (XI1.56) и (XI 1.59) и расчетов £0 точным формулам (XI 1.42), (XI1.47) и (XI1.50)

приведены в виде зависимости Q* и Q от т (рис. XI1.9).

На первой стадии вытеснения результаты приближенных и точных расчетов практически совпадают. Существенные расхождения начинаются после прорыва воды.

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

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

Вместо обычных относительных фазовых проницаемостей рассмотрим некото­ рые эффективные нелокальные характеристики.

р.Qb (01

в““ Qkbp(t) 9

F* —

Qh (0 l

H

 

Qk Дp (t) >

s

-

s

4-

V" (<)

S-

 

Sn+

Siml •

(XI 1.60)

(XII.61)

(XII.62)

В правые части уравнений (XI 1.60)—(XI1.62)

входят легко измеряемые

в эксперименте величины. Исключая из них время

можно установить зависи­

мости FI (s) и F* (s). Если эти характеристики окажутся достаточно универсаль­

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

т


Рис. XII. 10. Относительные фазовые проницаемости для нефти и воды

Для решения вопроса о влиянии отношения вязкостей на указанные эффек­ тивные характеристики воспользуемся приведенными выше формулами, получен­ ными в точной постановке задачи о вытеснении нефти водой.

Представим выражения (XI 1.60)—(XI 1.62) в безразмерном виде в соответ­ ствии с обозначениями (XII.31)

f• _ QF (se)

вАР

S = Sn + УцЗдесь

Т7 '- М У vи —

«. =

Q H -/(« /)]

' Н

ДрРо

 

 

. .

.

t S

« -

s„.

(XI 1.63)

(XII.64)

Подставляя в формулы (XII.63) значения Q из (XII.50) и VH из (ХП-64), получим окончательные расчетные формулы

FT =

Фeft

F* =

Ф « р - Ц

в

 

г(Ф/) ’

г(Ф,) ро

s =

sc +

\ ~ ! е

 

(XII.65)

Ф/

 

 

где 7 (se) = ft.

г?о


Отсюда видно, что FJ, F* и s зависят от времени т неявно через насыщенность в выходном сечении se.

Для расчетов /, Ф и г (Ф) использовали те же относительные фазовые про­ ницаемости, что и в приведенном выше примере. Расчеты были проведены для трех

значений |д0: 0,05; 0,1; 0,4.

), F* (s) представ­

Полученные в результате расчетов графики функций

лены на рис. XI 1.10. Здесь же для сравнения показаны обычные Fn(s) и FH(s). Как видно из рис. XII. 10, нелокальные F*, Z7*, как и следовало ожидать,

значительно отличаются от локальных FHи Fu, особенно для нефти. Для воды

при 5 > 0,7 Fв практически совпадает с Fв. Важной особенностью вводимых фильтрационных характеристик является то, что они практически не зависят от отношения вязкостей. Поэтому для их определения в широком диапазоне измене­ ния насыщенности можно использовать в экспериментах высоковязкие вытесняе­ мые жидкости. В тех же случаях, когда заданы FH, Fb, они вычисляются по изло­ женной методике.

§ 6. ПРИТОК НЕФТИ к СКВАЖИНАМ В НЕФТЕГАЗОВЫХ ЗАЛЕЖАХ С ПОДОШВЕННОЙ ВОДОЙ

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

Рассмотрим задачу о максимально возможном отборе нефти без допуска

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

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

Ниже излагается более общий расчетный метод определения предельного

дебита и соответствующего ему положения интервала вскрытия [22]. Предположим, что нефть из горизонтального пласта притекает в скважину

с интервалом вскрытия b (рис. XII. 11). До начала отбора нефти водонефтяной и газонефтяной контакты являлись горизонтальными плоскостями. При дебите^, не превышающем некоторого предельного, образовавшиеся конусы подошвенной воды и газа будут неподвижны и устойчивы. Положение интервала вскрытия, при котором безгазовый и безводный дебит достигает максимального значения, назо­ вем оптимальным. При таком положении интервала вскрытия конусы газа и воды 6УДУТ одновременно находиться в предельном устойчивом состоянии.

т


Рис. XII.11. Водяной и газовый конусы

Используя условие устойчивости конуса [46] и гидростатический закон распределения давления в газе и воде, можно написать систему уравнений:

ГдФ(г,гП

_

MVh — Vr)

 

(XI1.66)

L

dz J r=o« ?=yr

p

 

 

 

ГдФ (r,z) 1

_

kr (Yb — Yh)

 

(XII.67)

L dz

Jr=0, z=k—yB

|x

 

Ф (0.yr) = Ф (tfo.0) -

kr (v"r ~ Yr) Ur,

 

(VI1.68)

 

 

ф (0.ft -

yB) = ф (Яо.0) -

^ ^ В - Т я )

Ув

(XI1.69)

 

 

 

г

 

 

где Ф (г, г) — потенциал скорости фильтрации; уГ, ув — соответственно высоты газового и водяного конусов; kr — проницаемость пласта в горизонтальном на­ правлении; уг, ун, Ув — соответственно удельные веса газа, нефти и воды; р — вязкость нефти; R0 — расстояние от оси скважины, где распределение давления считается гидростатическим; h — толщина нефтенасыщенной части пласта на расстоянии R0.

Приняв допущение приближенной теории конуса о том, что отклонение поверхности раздела жидкостей от начальной плоской формы не влияет на рас­ пределение потенциала в нефтяной зоне пласта, для определения функции Ф (г, г), которая представляет собой решение уравнения Лапласа, воспользуемся потен­ циалом напорного течения к несовершенной скважине в однородно-анизотропном пласте постоянной толщины h.

Потенциал напорного течения при г = 0 можно представить в виде

ф(о,|) = Фо--£г(1,л,ЛтБ),

где

Ф0 = Ф(/?, 0); Q — дебит;

hr — расстояние от начальной плоскости раздела между нефтью и газом до середи ны интервала вскрытия; k* — проницаемость пдаста в вертикальном на­ правлении.

Ж


Безразмерная функция F при R > 1 имеет вид:

 

 

 

“ т ;

(ХП.?0)

 

 

 

 

(£ + Q,5й + йг) (£ + 0|5Я — йг)

 

H l.fi,Л,-) = In

(£ — 0,5ft — йг) (| — 0,5й + йг)

 

_.

Г[ 1 +0,5(6 + 0,5й + йг)]П1 - 0 ,5 ( |- 0 ,5 й + йг)1

 

 

Г[1 +0,5(5 —0,5й — йг))Г [1 —0,5(| г 0,5й + йг)]

 

 

Г[1 +0,5 (| + 0,5й — йг)] Г [ 1— 0,5 (| — 0,5й + йг)]

(XII.71)

*

П1 + 0,5 (5 — 0,5й + йг)] Г [ 1— 0,5 (5 + 0,5й — й,.)| ’

 

где Г — гамма-функция.

Эту задачу проще всего решать графически, для чего удобнее пользоваться

безразмерной формой уравнений (XII.66)—(XII.69):

 

Гд/ (£,/*,/*г"[

 

_ 4лЯу

 

 

 

L

l.n.Hr]

_ ;

= Кй

 

(XI 1.72)

 

J

\=yr

q

 

 

 

Ъ=уг

 

4л/г

 

 

_

г д/(5,й,йг) 1

 

= /С2;

(XI 1.73)

 

L

dl

J 1='—ув

 

 

 

/(£г.Мг) =Ki9r-

2h In

 

 

(XI 1.74)

/ (1 — £в,Мг) —

— 2/г 1п^ ,

 

(XI 1.75)

Здесь

 

 

 

 

___ :

 

 

 

 

 

 

 

 

 

 

 

 

4

kr (Ув — Ун) Ла

(XI 1.76)

 

Ун — Уг

 

 

 

 

V=

 

 

 

 

 

Ув — Ун ’ Уг ~ h ’

“ /I *

 

При заданных значениях у, /г и Я задача сводится к нахождению из системы четырех уравнений (XI 1.72)—(XI 1.75) четырех неизвестных величин: qt £п, ЛГ.

Построив графики функции / (£> Л, /гг) для различных значений h и /гг и проведя

касательные к ним из точек ^0, —2h X

Х1п

* -) и

^ 1,—2/г 1 п - ^ ,

находим

 

угловые коэффициенты этих касательных

 

К\ и /С2, а также уг и ув. Затем вычисляем

 

-

я

В

-

 

 

 

 

 

 

 

 

результате получаем зависимости

 

 

 

Уг = Уг(У),

Ув = Ув(У),

я —я (V)•

 

 

 

 

 

 

 

 

(XI 1.77)

 

 

Пример графического решения пока­

 

зан на рис. XII. 12, где изображена

кривая

 

/ (Л, 6Г, Б) для Л = 0,2, Лг = 0,4

и опре-

 

деляется £в, Рг*

Ki и

^2

при

2/i In — =

 

=

0,37.

расчетов

представлены

рис> \II.12. Распределение потенциала

на

Результаты

рис. XII. 13, а—е

В виде

графиков

на стенке скважины

223