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

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

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

Используя этп разложения, урав ѳнпе семейства изобар представим в виде

рт-\- В (Ѳ) рт~1-)-. . . = С,

П

ігде т = 2 2 = —2L по условию теоремы отлична от пуля. £=і

Дальнейшие рассуждения аналогичны приведенным в предыдущем дока­ зательстве.

«Точки-скважпны» для системы (II.6. Г)' являются особыми точками типа «узел». Интегральные кривые системы (II.6.1) — траектории п совпадающие •с нпмн линии тока в этих точках мало отличаются от пучка прямых с центром

в точке хі, у,-. Справедливость этого

заключения вытекает пз того,

что лпнип

'тока исходят пз источника п

перпендикулярны изобарам, а последние вблизи

 

 

 

 

скважины

мало

отличаются

от окружностей.

 

 

 

 

Можно

показать

также,

что точки покоя, где

 

 

 

 

др/дх =

0, др/ду = 0, для системы (II.6.1) явля­

 

 

 

 

ются особыми точками типа «седло» и

через этп

 

 

 

 

точки

проходят

две

или несколько

эквипотен­

 

 

 

 

циален и интегральных кривых.

 

 

 

 

 

 

Линиями поворота первой и второй пазовем

 

 

 

 

кривые, представляемые

соответственно уравне­

 

 

 

 

ниями

Fx (х, у) =

О,

 

(х, у) =

0, где Fx и F„

 

 

 

 

заданы системой (II.6.3).

первой

(второй) линией

 

 

 

 

При пересечеішп

с

 

 

 

 

поворота

линии

тока,

определяемые системой

 

 

 

 

(II.6.1), имеют касательную, параллельную осп

Тис. 21. Линии тока

при

оу {ох), если только точка пересечения не яв­

ляется особой точкой. В точке пересечения имеем:

работе двух равнодебнтных

эксплуатационных

сква­

 

 

dx

= 0;

 

dy

 

= 0.

 

 

жин.

 

 

 

 

dt

 

dt

 

 

 

-Линия поворота:

1 — первая;

 

 

 

 

 

 

 

 

Следовательно, dy/dx =

сю.

 

 

 

 

2 — вторая.

 

 

 

 

 

 

 

 

 

 

Затем через все «точки-скважпны» проходит

 

 

 

 

первая

(вторая)

линия

поворота и в этих точ­

ках имеет касательную, параллельную

осп оу {ох).

Поместим начало коорди­

нат в первую

скважину (хг, у±) п уравнение F1 (х, у) =

0 представим в сле­

дующем виде:

 

п

 

 

 

 

 

 

 

 

 

 

 

Fx {x, I/) = М

(x?+ yf 2Ж£Х 2УіУ+ х2 + у2) +

 

 

у2) В (х, у),

2

(х2 +

 

 

1=2

 

 

 

 

 

 

 

 

 

 

 

где В {х, у) — полином от х,

 

 

 

 

 

 

 

 

 

 

Отсюда имеем:

 

 

дРг (0, 0)

 

dFx (0,

0)

 

 

 

 

 

Т і(0,

0 )= 0 ;

:0;

=

0,

 

 

ято и доказывает утверждение.

дх

 

 

ду

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Высказанные предложения могут помочь построить картину течения жидко­

сти в пластах горных пород, вскрытых скважинам. Для построения линий тока

.можно опираться на следующие положения.

1. Чтобы построить направление поля системы (II.6.1), данную точку М (х, у) «соединим со скважиной А х {х1Уу г) и по этой прямой отложим отрезок

 

 

_________ ____________

^«х

 

 

 

Ѵ { х - Х 1)2 + { у - у 1)2

R l ’

 

где

В 1 — расстояние

от точки М до скважины, по направлению к скважине

при

стоке и от нее

при источнике. Геометрическая

сумма векторов М А г,

М А 2, . . ., М А п дает

направление поля, определяемого

системой (II.6.1), что

непосредственно следует из принципа суперпозиции.

60

2. Будем называть сепаратрпссами лпшш тока, проходящие через точки покоя, среди них особый интерес представляют те, которые ограничивают область притяжения к данной скважине. Через точку покоя, представляющую простое «седло», проходят две сепаратрпссы. Бслп примем точку покоя за начало, а направление вектора скорости за ось ох, то уравнение сепаратрпссы вблизи точки покоя может быть найдено в виде

у = Cjx2-]- С2х3,

причем коэффициенты определяются из дифференциального уравнения (II.6.4). Вдали от точки покоя уравнение строится методом, изложенным в (II.6.1).

3. Линии тока исходят из нагнетательных скважин и кончаются в эксплу­ атационных. Бесконечно удаленную точку принимают также за скважину с ха­ рактеристикой

п

£ = - 2 ^ . (=і

При пересечении линий поворота (первой п второй) линии тока меняют направление и имеют касательную, параллельную оси оу (ох).

Липпи тока не пересекают друг друга (и сепаратрисе), кроме как в особых точках. Приведем примеры.

Пример 1. Скважины расположены в точках (а, о), (— а, о) и имеют харак­ теристики Х± = Х2 = —1. Система (II.6-1) запишется в следующем виде:

 

dx

 

Xа

 

 

27-f- а

 

 

 

dt

 

(x -a )2 -\ -iß

(х-)-а)2-|-у2 ’

 

d y

=

_______________У_________________________ У

 

 

откуда

dt

 

—а)'2-)-і/2

(x-j-a)a-|-y2 ’

 

 

 

dy

У(.r2 -fy 2 -)-a2 )

 

 

 

 

 

 

 

 

 

 

dx

X (х2 -{- 1/2— 0,2)'

 

 

Первая линия поворота состоит из оси оу и окружности

 

 

 

 

а-2 +

у2 =

а 2 .

 

 

 

Вторая линия поворота — из осп ох. Линии тока

показаны на рис. 21.

Общий интеграл уравнения будет

 

 

 

 

 

 

 

 

 

х2 +У 3 + Сху =

а2.

 

 

Пример 2.

Скважина

(а, о) — эксплуатационная

с характеристикой — Аі,

(—а, о) — нагнетательная с характеристикой Х2 (—а, о).

В бесконечно удален­

ной точке сток обильности Х2 — X.x. Система (а) имеет вид:

 

 

dx

 

Хг (х — а)

 

X2 (x + a )

 

 

 

dt

 

(x — a)2+

A2

1

(x-j-a)2 +

y2

>

 

dy

 

hV

 

1

X2y

2/2

 

откуда имеем

dt

 

(x —a)2+!/2

1

(X + O) 2 +

>

 

 

 

 

 

 

 

 

 

dx

—Xi (x — а) [(х + а)2 +

і/2 ]4

2 У [(x —a)2-

dt

 

 

[(* -a )2 + p 2 ][(* + a )2 + 0 ® ]

 

 

dy

ЯхУ[(а;+ а)2 +

у2] + Я2Р[(а: —a)2 + y2]_

dt

 

[(я —а)2 + р2] [(® +0)2+у2]

 

 

dy

________ Хгу [(х — а) 2

+ р2

] _ А 2г/((х — д) 2 4 -У2]______

dx

А,! (а —а) [(х +

я.)2

-)-і/2 ] _ Х 2 (х —а) [(а —а)2 + і/2 ]

61

Первая лшшя поворота будет:

 

 

*■1 (*г —а) [(-г*+

«)3 + І/3] — ^2 (®+<0 [(а; —о)2 +

1/2] = 0 ,

откуда

а)

 

у2 [Ѵ? (,г’ + а)

 

д2) [Хо (#

Я-і (а;+ о)] +

(а; —о)] = О,

плп

 

Ъ - х

 

 

 

 

 

 

где

 

6 + .г- (х'і- аі) ,

 

 

Ло “{-

 

 

 

 

а = Ь.

 

 

 

Х>2 — ^

 

Кривая состоит из

двух ветвей: одна — замкнутая, пересекающая ось оаг

в точках (а, о) и (Ь, о),

другая пересекает ось ох в точке х = —а и удаляется

в бесконечность при х =

•—Ь.

 

 

Рис. 23. Линии тока при работе двух равнодебптиых нагнетательных сква­

скважины. жин. Обозначение см. па рис. 21.

Вторая линия поворота будет

h y [О® + а)2 + уЦ Ъ2У [(а- — я)2 + т/2] = 0.

Она состоит из оси ох и кривой

[(а- + о)2 + 2/2]_Х27/ [(я —а)2 + і/2]=0.

Упрощая, получим:

 

 

(я — b)2 + J /2 = b 2 _ a2.

Это

окружность

с центром в точке Ъ и радиуса І^Ь2 — а2. Изображение

течения дано на рис.

22.

При

= Х2 уравнение интегрируется в квадратурах, общий интеграл его

*2 + (1/— С)2 = а2 + С2.

Интегральные кривые представлены на рис. 23.

Пример 3. Пусть три скважины расположены в точках (я, о), (—а, о), (о, о) и имеют равные характеристики — единицу. Первая линия поворота предста­ вляется следующим уравнением:

Я [Зу4 — 6я2і/2_|_3^4 — 4о2я2 -]- я4] — Q

62

II состоит из осп оу и двух замкнутых кривых, пересекающих ось ох в точках {± а , о), (± а /] /з , о). Сепаратрпссы, проходящие через точки покоя ( ± а / у 2 , о),

в первом приближении имеют следующее уравнение:

яЗ

1Г ’

Изображение течения дано па рис. 24.

Рис. 24. Линии тока при работе

Рис. 25. Линли тока при работе двух

трех эксплуатационных

скважин.

нагнетательных и одной эксплуатацион­

Обозначения см. на

рис.

21.

ной скважины.

 

 

 

Обозначения см. на рис. 21.

Две эксплуатационные

скважины

с характеристиками ^ расположены

в точках (о, а), (о, —а) и одна эксплуатационная скважина расположена в точке (о, о) и имеет характеристику —Х0. Первая линия поворота будет:

(2Х]. -

%о) г/1 + [2h (2*2 +

я2) - 2h (xi - Я2)] учд

2h ^ (*2 + я2) -

h

2 +

я2)2 =

о.

При 2Л.Х— %0

0 — кривая замкнутая,

пересекающая

ось

ох

в точке

(± а ,

о), а ось оу — в точке ( ±ß,

о), где

 

 

 

 

 

 

а

.

во

„о —(^о4-^і) + 1^?+4Ао?ч

 

 

 

 

а - V ü ^ h ' Р

 

2 h = T o

 

 

Вторая линия поворота также замкнута и пересекает ось ох в точках ( ± у ,

о),

где

 

 

 

_________

 

 

 

 

 

 

 

о _

^-о +

^і + 1^^і-Ь4Х0Х1

 

 

 

 

 

 

 

Ѵ '

 

2?ы-Хо

 

 

 

 

 

Изображение течения дано на рис. 25.

Г л а в а III

СОВМЕСТНАЯ ФИЛЬТРАЦИЯ ДВУХ НЕСМЕШИВАЮЩИХСЯ ЖИДКОСТЕЙ

ВПЛАСТАХ ГОРНЫХ ПОРОД

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

Например, на рис. 10, а наблюдается замещение нефти водой по всему разрезу пласта с образованием обводнеппых зон с высокой

Н | /

2

Рис. 26. Расчленение фильтрации жидкости по мощности пласта.

1 — нефть; 2 — вода.

степенью вытеснения нефти, причем здесь не видим подъема перво­ начального положения ВНК, находящегося в нижней части пласта. На рис. 10, б, в происходит постепенное замещение нефти в нижних горизонтах, которые можно рассматривать как подъем ВНК. На рис. 10, г наблюдается послойное движение нефти и водонефтяной смеси. Еще более расчлененная по мощности картина движения жидкостей наблюдается на рис. 26 для екв. 183, 184, 54 Покровского нефтяного месторождения.

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

64

в том числе сотрудниками лаборатории, Казанского университета. Приведем результаты опытов Б. И. Плещинского и др. [80, 81].

Для установления структуры течения был взят небольшой лоток с внутренними размерами 20 X 50 см, позволяющий набрать модель пласта мощностью до 2,2 см. Лоток изготовлялся из органического стекла и закрывался стеклянной крышкой. В качестве жидкости, моделирующей нефть, применялась смесь керосина и а-монобром- нафталина или масел с коэффициентом преломления света, равным коэффициенту преломления стекла. В качестве жидкости, модели­ рующей воду, бралась подкрашенная вода. Жидкость подводилась и отбиралась с помощью галерей, представляющих 10-мм трубки, рассверленные по стороне, обращенной в сторону рабочего поля модели, и прикрытые сеткой. Путем укладки стеклянной крошки разных размеров послойно создавалась одно-, двух- и трехслойная модель. Жидкостью модель заполнялась под вакуумом.

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

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

Сделаны следующие выводы.

1.Распределение скоростей движения частично на границе жидко­ стей близко к нормальному (гауссовскому).

2.При отсутствии остаточной воды и отношениях вязкостей

р,н/р,в =

1

-г 2 (см.

рис. 27, а) ВНК движется сплошным фронтом.

3. При Цн/р-в =

3-^-4

фронт вытеснения рвется (теряет устой­

чивость),

 

образуя

целые

необводненные участки (см. рис. 27, б).

При рн/ра

> 4

вода

движется отдельными струйками (см. рис. 27, в

е). Для

рн/рв

= 2 — 3 (в

промежутках интервала) опыт мог быть

отнесен к тому или к другому случаю.

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

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

5 Заказ 322

65

 

скорости движ ения воды по слоям размой проницаемости вы равни­

ваются (см.

рис. 27, е).

П еред фронтом

движ ения

нагнетаемой

п од ­

краш енной

воды

часто

наблю дается

движ ение

вала

погребенной

прозрачной

воды.

 

 

 

 

 

 

Перечисленные

обстоятельства движ ения границы

раздела

воды

и нефти ' необходим о учитывать при построении математической модели вы теснения нефти краевой водой.

Рис. 27. Результаты экспериментов на прозрачных моделях по вытеснению сверху, нижние — вид

Модели: а, г, 0, е — однослойные; б,

66

5*

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

пользуются несколькими схемами или моделями,

описывающими

это движение.

В зависимости

от

структуры потоков, характерных

особенностей строения

нефтяных

коллекторов,

свойств жидкостей и

целей

расчетов

эти модели должны

меняться

иногда

в

пределах

одной расчетной

площади.

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Построение

универсальной

модели,

описывающей все случаи

 

 

 

 

 

 

 

движения,

едва ли целесооб­

 

 

 

 

 

 

 

разно ввиду усложнений, ко­

 

 

 

 

 

 

 

торые вносятся

вместе с уче­

 

 

 

 

 

 

 

том

дополнительных

факто­

 

 

 

 

 

 

 

ров.

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Можно пойти другим пу­

 

 

 

 

 

 

 

тем: разработать унифициро­

 

 

 

 

 

 

 

ванный ряд моделей, в ко­

 

 

 

 

 

 

 

тором можно вести счет по

 

 

 

 

 

 

 

одной

из отработанных схем,

 

 

 

 

 

 

 

в

пределах

целого

пласта

 

 

 

 

 

 

 

(месторождения)

 

или

его ча­

 

 

 

 

 

 

 

сти,

и переходить

к

другой

 

 

 

 

 

 

 

при

изменении контрольных

 

 

 

 

 

 

 

условий.

 

Переход

 

должен

 

 

 

 

 

 

 

быть

 

выполним

 

без

суще­

 

 

 

 

 

 

 

ственного

нарушения

всей

 

 

 

 

 

 

 

модели и только путем за­

 

 

 

 

 

 

 

мены

 

ее

части

(блока),

при

 

 

 

 

 

 

 

этом

следует

иметь

возмож­

 

 

 

 

 

 

 

ность использовать

основной

 

 

 

 

 

 

 

массив параметров

и

исход­

 

 

 

 

 

 

 

ных

данных.

Для

выполне­

 

 

 

 

 

 

 

ния

последнего

условия

не­

 

 

 

 

 

 

 

обходимо,

чтобы методы из­

 

 

 

 

 

 

 

мерения

параметров и исход­

 

 

 

 

 

 

 

ных

данных

были

унифици­

 

 

 

 

 

 

 

рованы

в

данном

 

случае,

 

 

 

 

 

 

 

чтобы

входящие

величины

 

 

 

 

 

 

 

измерялись

одинаковым

ме­

 

 

 

 

 

 

 

тодом.

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Такой ряд вместе с зада­

 

 

 

 

 

 

 

нием

 

признаков

 

перехода

 

 

 

 

 

 

 

может играть

роль

универ­

 

 

 

 

 

 

 

сальной

схемы

 

счета

при

 

 

 

 

 

 

 

большой конкретизации и уп­

керосина,

водой

(левые

рисунки — впд

рощениях внутри каждой мо­

снизу).

 

 

 

 

 

 

дели ряда. Все это существен­

я — трехелойиые.

 

 

 

 

 

но

 

упрощает

составление

67

программ для ЭВМ, давая возможность придать им блочную структуру и комплектовать программы из отработанных блоков.

В настоящей главе предлагается ряд из 12 моделей, который удовлетворяет сформулированным требованиям. Классификация проведена в зависимости от структуры (слоистости) пластов и от вида функции Баклея — Леверетта (см. § 2).

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

темы,

отсылая читателя

к работе

[126].

Для

написания

§ 1, 5, 6

использованы работы

[26,

53,

78,

142,

143, 147], § 3 — [137,

138,

145],

§ 4

- [1 2

-1 4 ,

43,

44,

105,

112],

§

10 -

[37, 38, 41],

§

12 -

[1, 2,

26,

86, 87,

90, 120,

149]

и др.

 

 

 

 

 

Желающих ознакомиться с более общими в определенном смысле моделями отсылаем к работам [70, 50]. Отметим, что при исполь­ зовании этих моделей объем вычислений весьма велик, что заставляет искать пути его уменьшения, например упрощая схемы.

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

§ 1. Схема трехмерной фильтрации с пропиткой

Рассмотрим движение двухжидкостных систем.

Считается, что в каждой точке пласта существует два давления — нефти рі и воды рв. Разность этих давлений рн — pt = рк объяс­ няется наличием капиллярных сил и называется капиллярным скачком.

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

Vix =

kxk*i

dp\

.*

 

 

 

Pi

 

dx

 

 

 

Vl y =

-

kyk*

dPi

= н,

B)

(III.1.1)

Pi

dy

 

Viz =

k7k*

(

dp*

 

 

Pi

1

dz

 

f p i g )

 

Равенство і = н, в означает, что і принимает значения н или в. Введем приведенное пластовое давление pt для нефти и воды:

. .

Рі = Р*і + Pigz.

(III.1.2)

68

Тогда, подставляя р* в (III.1.1), получаем

 

dpi

 

I

Рі

дх

 

 

k y k t

dpi

А-тА-г* V pp,

(IIU .3)

Pi

Оу

Рі

 

A'zA'i

др;

 

 

Pi

dz

 

 

Для капиллярного скачка будем иметь следующее выражение:

 

 

Рк = Р н ~ Р„ — gz(pH— рв).

 

 

 

 

(III.1.4)

Уравнение неразрывности и состояния запишутся:

 

 

~

Pu) + ^

{ ? Ы « ) + 4 г

 

+

aiPötm)

= °;

 

(ІП -1-5)

 

р =

р(р); т — т (р);

к* = к*(р)-,

 

 

 

(III.1.6)

 

К = кЛр)'> ку=ку(р); кг =Ісг(рУ,

 

 

(III.1.7)

 

 

 

s„ + sB=

i.

 

 

 

 

 

(III.1.8)

Подставим теперь

в

(III.1.5)

значения

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

из (III.1.3) и продифференцируем последний член. Тогда

получим:

д

 

 

к у к І _ д р і _ _ \ , ± _

(

_ k ^ l

dPi \

дх

Рі

 

Pi

 

ду ) ~ г

dz

\

vt

Ц,.

dz

J

 

 

 

dm

1

dpi

 

 

dsi

 

 

(III.1.9)

 

 

 

dpi

J

dt

9im

 

dt

 

 

Разложим функции, определенные формулами (III.1.6)—(III.1.8),

в ряд Тейлора и оценим члены, входящие в (III.1.5),

аналогично

гл. 2, § 2. Тогда получим:

 

 

 

 

 

 

 

 

д_ /

кхк?

dpi

\

д

/

кук*

dpt \

. д

/

кгЩ

dpi

\ _

дх \

р(-

дх

)

ду

\

р,-

diI )

‘ dz

\

ці

dz

)

 

 

=

Si (nzß* +

ßc)

 

4 r

 

 

(III.1.10)

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

PK= PK(S); S = sB.

(III.1.11)

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

А* = **(«),

(III.1.12)

из (III. 1.4) получаем:

 

Рн — Рв = Рк (s) — gz (р„ — р„).

(III.1.13)

69

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