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

1464

.pdf
Скачиваний:
14
Добавлен:
15.11.2022
Размер:
27.36 Mб
Скачать

Исходя из уравнения (13, XI). имеем:

Рк ~ Qr — Рат^

1 “I"' Р { Р^к

^Рат) "Ь 2 0Л (Рк

Рат)

1 + /?(рг Рат) "Ь 2 Р (Рг — Рат)

(14, XI)

^ =

=

Рат/^(Рк

Р г)(1 “Ь P p ')i

 

где р' — среднеарифметическая величина между избыточными давле­ ниями на контуре питания и в галлерее.

~/ Рк ~ЬРг

Рат*

Р = “ 1 —

 

Из формул (9, XI) и (5, XI) легко найти объемный расход жидко-

сти Q = ^

Q

kF Рк

Рг

1

(15, XI)

(За

LK

Q'

 

 

Как видно из формулы (15, XI), объемный расход сжимаемой жид­ кости в случае установившегося движения есть величина переменная, поскольку плотность жидкости д есть функция х (см. формулу 11, XI).

Разделив объемный расход жидкости на площадь F, найдем ско­ рость фильтрации v, являющуюся также переменной величиной:

_к_ Рк ~ Рг 1

(16, XI)

Рр LK Q*

Знак минус в правой части формулы (16, XI) указывает, что ско­ рость фильтрации направлена противоположно направлению оси х.

Подставляя в формулу (15, XI) вместо д величину £ат, полу­ чим приведенный к атмосферному давлению объемный расход жид­ кости фат

Фат

 

kF

(вк ~ вг)

(17,

XI)

PfLQ&T

LK

 

 

 

Подставляя сюда вместо

(дК — дг) его значение

из уравне­

ния (14, XI), имеем:

 

 

 

 

 

Фат

kF

(QK -

рг) (1 + (}р').

(18,

XI)

 

v-

ь,к

 

 

 

В случае же установившейся одномерной фильтрации несжимае­ мой жидкости объемный расход жидкости QH (см. формулу 14', IX) равен:

Q» = ^ЕГ

(19’ Х1)

Сравнение формул (18, XI) и (19, XI) показывает, что при одинако­ вых граничных условиях приведенный к атмосферному давлению объ­ емный расход сжимаемой жидкости QaT отличается от расхода несжи­ маемой жидкости QHл и ш ь на множитель (1 +/?р'). Но величина /Зр' ма­

ла по сравнению с 1. Так при /3= 5-10“ 5^ и р' = 100 am, (Зр' = 5 •Ю_3. Поэтому при решении большинства практических задач с достаточной точностью [учитывая точность определения других коэффициентов, входящих в уравнение (16, XI) и в уравнение (17, XI)] можно считать, что расход сжимаемой и несжимаемой жидкости при одинаковых гра­ ничных и прочих условиях одинаков и для определения расхода сжи­ маемой жидкости надо пользоваться формулой (17, XI).

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

приближенное равенство (14, XI), имеем:

 

 

dQ _ Рк ~ Qr rsj £ат/3(Рк Рг)(1 Н~ (Зрг)

 

 

dec

LK

LK

 

 

Подставляя это значение градиента плотности в формулу (3', XI),

получим:

 

 

 

 

d p ^ e ат (Рк ~Рг)(1 + /V )

(20,

XI)

dx

Р

Ьк

 

 

В случае же одномерной фильтрации несжимаемой жидкости гра­

диент давления

согласно формуле (15', IX) равен:

 

 

 

( dp\

_ Рк ~Рг

(21,

XI)

 

\ d x )H

LK

 

 

 

Разделив (20, XI) на (21, XI) и учитывая малую сжимаемость жид-

кости, имеем:

dp

 

 

dx

J L ( I+ f& )sn .

(22, XI)

 

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

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

§ 2. Радиальная установившаяся фильтрация сжимаемой жидкости

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

dp

(23, XI)

dr'

Из формулы (3, XI) градиент давления будет:

ф

= л 1

(24, XI)

dr

Р Q dr

 

CL D

Подставляя это значение ^ в уравнение (23, XI), получим:

к Ф

(25, XI)

Рр dr'

 

Массовый расход жидкости QMполучим, умножив модуль массо­ вой скорости gv на площадь сечения F = 27гг6, где Ь— мощность пла­ ста:

QM= |gv\ F =

2nkb r dg

 

Рр dr

Разделяя переменные g и г, имеем:

P^QM dr

(26, XI)

2тгкЬ г

Граничные условия формулируются следующим образом: при

r = R c

в = в с ,

г = Дк

(27, XI)

Q = QKI

где дс — плотность жидкости на контуре скважины радиуса Д^

gK— плотность жидкости на круговом контуре области питания радиуса Д^

Интегрируя уравнение (26, XI) в пределах от дс до дКи от Дс до Дк и решая полученное уравнение относительно QM, найдем формулу мас­

сового расхода жидкости в виде:

 

 

 

 

 

 

27Гkb {QK — Qc)

(28,

XI)

Q M

=

 

 

 

 

Pt*

In

 

 

 

 

 

 

ь нГ

 

 

Проинтегрируем уравнение (26, XI) в пределах от дс до д и от Дс

до г:

 

 

 

 

 

 

f

j

PH'QM f

dr

 

 

J

d e =

b M

j

- •

 

 

Qc

 

 

R c

 

 

 

отсюда

 

 

 

 

 

 

 

 

 

 

 

(2S'

XI)

Но из формулы (28, XI) имеем:

 

 

 

 

PI^QM _ QK ~ Qc

27Гkb Дк

1пж

Подставляя это выражение в уравнение (29, XI), получим:

е = г, + ^ г 1 тг

(зо. XI)

14 Подземная гидравлика

Формулы (29, XI) и (30, XI) дают изменение плотности жидкости в пласте в направлении г. Для определения давления в различных точ­ ках пласта нужно найденные по формулам (29, XI) и (30, XI) значения д подставить в уравнение (3, XI).

Найдем объемный расход жидкости Q:

2 п k b (вк

в с )

1

(31, XI)

Рг

R

Q'

 

Разделив объемный расход Q на площадь F = 2пгЬ, найдем ско­ рость радиальной фильтрации жидкости, направленную противопо­ ложно направлению оси г:

v — — к

( в к

в с ) 1

(32, XI)

Рг

R

£г *

 

1п7Ггьс

Подставляя в формулу (31, XI) д = £ат, получим величину приве­ денного к атмосферному давлению расхода сжимаемой жидкости Q&T:

Q&T

2 п k b

QK — Qc

(33,

XI)

Рцвъ?

,

RK

 

 

'"Ж

 

 

По аналогии с формулой (14, XI) имеем:

 

 

(0к - &с) = Рат/?(Рк -

Рг)(1 + № ) ,

(34,

XI)

где

 

 

 

 

 

Р' = Рк +Рс

 

Рат*

 

 

 

2

 

 

 

 

Подставляя это значение (дК— дс) в уравнение (33, XI), получим:

2тгk b

(#к

Qc)

(35,

XI)

Q&T

1 Дк

(1 + № ■

*

 

 

 

1пж

Формула (35, XI) отличается от известной формулы (21, IX) для дебита несжимаемой жидкости на постоянный множитель (1 + /?р').

Поскольку, как указывалось в § 1, величина /Зр' мала по сравнению

сединицей, при решении практических задач дебит (объемный рас­ ход Q) скважины при радиальной фильтрации сжимаемой жидкости

сдостаточной точностью можно определить по формуле Дюпюи

Q =

p i , щ

1пж

Произведя далее вычисления, аналогичные приведенным в § 1, по­ лучим формулу (22, XI), из которой следует, что кривые распределения давления в пласте при радиальной установившейся фильтрации сжи­ маемой и несжимаемой жидкости практически совпадают.

В связи с изложенным очевидно, что все соображения, приведен­ ные в § 2 и 3 главы IX, о форме изобар и линий тока, о зависимости дебита скважины от ее радиуса и расстояния до контура питания, о «во­ ронке депрессии» и пр. справедливы и для установившейся фильтрации сжимаемой жидкости в аналогичных условиях.

§3. Приближенное решение задачи

онеустановившейся радиальной фильтрации сжимаемой жидкости

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

Выясним, что происходит в пласте при вскрытии его скважиной в условиях упругого режима. Схема пласта изображена на рис. 54.

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

жину.

Назовем этот процесс увеличения «радиуса действия» RKскважи­ ны или, иными словами, расширения «воронки депрессии» первой фа­ зой неустановившегося движения1. Характерной особенностью ее яв-

1 Величина RK является условной и определяется по методу последовательной

смены стационарных состояний (см. дальше).

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

С того момента, когда радиус воронки депрессии гк достигнет кон­ тура области питания, на котором в условиях водонапорного режима давление рк можно во многих случаях считать постоянным, движение жидкости станет установившимся и при решении практических задач можно (как указывалось в § 1 и 2 настоящей главы) считать жидкость несжимаемой.

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

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

1. Случай, когда дебит скважины является постоянным

Пусть ко времени t (время, истекшее с момента ввода скважины в эксплуатацию) падение давления распространилось на расстояние RK от скважины. Тогда масса жидкости, извлеченной из пласта за время t, равна QM£, где массовый расход QM= const. С другой стороны, мас­ са извлеченной из пласта жидкости равна изменению массы жидкости, находящейся в пласте.

Следовательно,

QMt = I 2 ( Q K - ?),

(36, XI)

где i? — объем порового пространства области пласта, ограниченной окружностью радиуса Дк,

где дк — плотность жидкости в точках пласта, для которых координа­ та г ^ RK; QK = const;

Q — средневзвешенная по объему Q плотность жидкости

Q

(38, XI)

h j * * 1'

Остальные обозначения прежние.

Подставляя в уравнение (36, XI) вместо массового расхода жид­ кости QM и объема порового пространства ft их значения из фор­

мул (28, XI) и (37, XI), получим:

 

 

 

 

 

Ь£Ь

= n{R2 _ R2)bn{gK _ g).

(39,

XI)

^

In #*

 

 

 

 

 

 

 

it с

 

 

 

 

 

 

Обозначим

 

 

 

 

_

 

 

 

е = £±-

R* =

ё* .

?

£ =

(40,

XI)

 

9к’

к

Д .’

вк

v

'

Тогда уравнение (39, XI) после сокращения на 7г6^к и вынесения

за скобки Щ можно переписать в виде:

 

 

 

 

 

J 1

= 1п Дк( К 2

-

1)(1 - О-

(41,

XI)

Для нахождения величины £ определим средневзвешенную по объ­ ему плотность жидкости д. Выделим на расстоянии г от скважины кольцевой элемент пласта шириною dr (рис. 76). Объем порового про­ странства этого элемента равен:

dft = 2nbmrdr.

42, XI

Подставляя значения ft и dft, приведенные в формулах (37, XI) и (42, XI), в уравнение (38, XI), получим:

2

JД к g(r)rdr.

Д к - Rc Rc

Подставляя вместо д(г) его выражение из формулы (30, XI), раз­ делив уравнение (43, XI) на дк и вводя переменные e, R* и f , значения которых даются формулами (40, XI), получим:

i ‘ l ^ i J ( e + h k ,,' R) R'JR'-

<44>Х1)

к1

где R* = -щ.

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

e = i - ( i - £ )

1

 

1

(45,

XI)

2 Inд*

R :2 - l

 

г

 

что дает

 

 

 

 

 

1 - £ = (1 — е)

 

 

1

(46,

XI)

21пД*

R

f - 1

 

 

 

Подставляя полученное значение (1 - е )

в уравнение (41, XI), име­

ем:

 

 

 

 

 

2А(1 - е)

 

 

1

 

 

* = 1 п Я Ж 2 - 1 ) ( 1 - е )

 

Д*2 - 1

 

 

 

21пД*

 

И Л И

Rt2 -

1

 

 

 

 

 

 

 

2k ,t = ^K

- I n Rt.

 

 

mPfiR2

 

 

 

 

 

 

 

 

Откуда

 

 

 

 

 

mPuRZ ( Д*2 -

1

 

(47,

XI)

t =

 

- I n Rt

2k)

 

 

 

 

 

Вводя безразмерное время

 

 

 

 

 

т =

 

t

 

(48,

XI)

 

mPuR2

 

 

 

получим:

Рис. 66. Зависимость безразмерного радиуса воронки депрессии R* от без­ размерного времени т в случае эксплуатации скважины при постоянном де­ бите (Q = const).

На рис. 66 приведена кривая зависимости безразмерного радиу­ са воронки депрессии R* от безразмерного времени т, построенная по формуле (49, XI).

Зная значения R* в различные моменты времени, легко най­ ти соответствующие им величины забойного давления рс. Из форму­ лы (21, IX) имеем:

(50, XI)

Рс ~ Рк 2тгkb

где Q — объемный расход жидкости.

Распределение давления определяется по формуле (23, IX)2.

2. Случай эксплуатации скважины при постоянном противодавлении

2По формуле (23, IX) Дюпюи можно определять фиктивное распределение дав­ ления, а не истинное.

Соседние файлы в предмете [НЕСОРТИРОВАННОЕ]