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

книги / Физика и гидравлика нефтяного пласта

..pdf
Скачиваний:
1
Добавлен:
12.11.2023
Размер:
6.42 Mб
Скачать

ала

ф на цилиндрической поверхности радиуса R c• Решение та­

кой

задачи является сложным.

Зададимся

постоянной

интенсивностью q вдоль всей оси

скважины {2hq = Q o ). Тогда интеграл примет вид

 

Ф(г,

г) =

Qo

У ( г +

/г)2 + л2 +

(г +

/г)

(III. 140)

8лЛ ■In

У (2 -

й)2 4 - /-2 +

(2 -

Л)

С удалением от линии стоков (источни­ ков) потенциал непрерывно стремится

кнулю. В этом заключается отличительная особенность течения к линейному стоку (ис­ точнику) в бесконечном пространстве по сравнению с плоскопараллельным течением

ктакому же стоку (источнику). В послед­ нем случае, как известно, потенциал из­ меняется по логарифмическому закону,

стремясь

к ± о о в зависимости от знака

величины интенсивности q.

Если

предположить, что подлогариф­

мическое выражение равно постоянной ве­ личине, то получим уравнение эллипсоида вращения с фокусами в точках (0, 0, А) и (0, 0, —А). Подберем достаточно узкий эл­ липсоид, соответствующий / 3 >г, и заменим его цилиндром радиуса Ro и равновеликого

объема (рис. III.12).

Н. К. Гирин-

Тогда получим формулу

ского

 

2жро

Q =

 

In

I.6A ’

 

Яо

я I O=jQ0

///////.. у „////////

L

г

Рис. III.12. Схема к рас­ чету фильтрации жидко­ сти к скважине, вскрыв­ шей пласт бесконечной мощности

(III. 141)

где фо— приближенное значение разности потенциалов на по­ верхности фиктивной скважины радиуса R c и в бесконечности

( ф ( о о ) ) = 0 .

 

 

 

 

 

Фо = ф(/?0)

ф(°°,

оо) = ф (Rc),

(III. 142)

где фо— потенциал на этой поверхности.

на поверхности

Возьмем две точки

А и В соответственно

z = 0 и оси линии стоков г =

0 .

 

 

Для этих точек получим следующие значения потенциала из

(III.140):

 

 

 

 

 

Ф д =

<?о

,

л Л

* + ' * + А

 

Q ,

I n ------ 7— — ---------------

 

 

8лЛ

 

^ h 2+ r \ - h

(III. 143)

_

Qo

Z« + A

 

VB— W T ln- i ,—h

 

где га и гв — координаты точек А и В.

91

Пусть точки А и В настолько далеки от середины линии сто­ ков (z = 0, г = 0 ), что можно считать гд > /г и гв>/г.

Тогда из (III.142) имеем

 

 

 

(III. 144)

Фв =

Q

In

Q

8яЛ

4nzB *

 

 

 

(III. 145)

Но формулы вида Q/{4nrh) (III.144) и (III.145)— это выра­ жения для точечного стока (источника).

Следовательно, на достаточно большом расстоянии от линии стоков (источников) распределение потенциала сходно с рас­ пределением потенциала, возникающего от точечного стока (ис­ точника) .

Приведенное выше можно рассматривать как случай филь­ трации жидкости к скважине (или от нее), вскрывшей пласт очень большой толщины (теоретически бесконечно большой) на некоторую величину Л. При этом плоскость z = 0 (плоскость симметрии потока) служит кровлей пласта. Дебит этой сква­ жины, очевидно, равен половине дебита Qo линейного стока (ис­ точника) длины 2 А.

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

Пусть между двумя непроницаемыми параллельными плоско­ стями (кровлей и подошвой пласта) расположен в произвольной точке сток (источник). Если бы сток (источник) находился в бесконечном пространстве, то потенциал скорости имел бы вид <р — <7(4яг).

Для одной непроницаемой границы (например, кровли) до­ статочно было бы одного отображения этой скважины относи­ тельно кровли, т. е. такого же стока (источника) на расстоянии £ от последней (рис. III.13, а). При этом потенциал в любой точке полупространства получился бы путем суперпозиции обоих потоков (реального и фиктивного — «зеркального»).

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

92

жны быть такими же, как и у действительного стока (источ­ ника) (см. рис. III.13, б).

Только при такой системе отображений выполняется условие непроницаемости кровли и подошвы пласта. Таким образом, мы получим потенциал, возбужденный лишь одним точечным сто­ ком (источником) q между двумя непроницаемыми плоско­ стями. Реальную же скважину следует аппроксимировать неко­ торой воображаемой линией, каждый элемент которой погло-

Рис. III.13. Схема бесконечного числа отображений для определения дебита скважины, не полностью вскрывшей пласт

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

В результате этой операции определим потенциалы реаль­ ного и бесчисленного множества отображенных линейных сто­ ков (источников).

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

93

Задача, очевидно, не может считаться оконченной, если не задаться видом распределения интенсивности потока q{Z). Ведь ясно, что один и тот же дебит можно получить при любом рас­ пределении q(t,) по линии стоков, если соблюдается условие

ь

Q = \q(Qdt,

о

где Ь— вскрытая толщина пласта.

Однако не при всяком распределении интенсивности потока q(t,) выполняется необходимое условие постоянства потенциала на поверхности питания и в скважине.

Соответствующим подбором q{t,) Маскет получил следующую

приближенную

формулу для дебита скважины, не полностью

вскрывшей пласт:

 

 

 

 

п

 

 

 

 

h (Фк - Фс)

 

 

Ч ~

I

Го.

.

Г (0,8756) Г (0,1256)

1

4А ’

 

28

[

Яс

 

Г ( I -0,8756) Г ( I - 0,1258)

J

In RK

 

 

 

 

 

 

 

(III. 146)

где Фк и Фс — потенциалы цилиндрических поверхностей пита­ ния и скважины соответственно; RK и Яс — радиусы поверхности питания и скважины б = 6/А; Г — гамма-функция.

Функция

Ф (6) = In

Г (0,8758) Г (0,1258)

(III. 147)

Г(1 - 0,8758) Г (I - 0,1258)

Формуле (111.146) можно придать следующий вид:

 

Q =

h (Фк — Фс)

(III. 148)

 

In Як

 

 

Второй член в знаменателе называется фильтрационным со­ противлением, которое обусловлено несовершенством скважины по степени вскрытия:

 

4h

 

 

. - - И

2 In Яс -ф ( 8) ) —

In Яс

(III. 149)

При совершенной скважине 6 =

1, ф (6 ) = 0 , Ct =

0 мы при­

ходим к формуле Дюпюи.

 

 

 

 

Для малых значений аргумента имеем зависимости

 

l n r ( l + * ) « - L ln H i n ^ - Y * .

(Y =

0,577).

(III. 150)

Кроме того,

 

 

 

 

 

Г (*) Г ( 1 х) =

л/sln ллс,

Г (1 +

х) =

лсГ(лс).

(III. 151)

94

Используя эти

выражения,

преобразуем формулу Маскета

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

 

 

 

 

 

п

________________ 2 я Л ( Ф к — Ф с )

Лк

(III. 152)

 

I ,

I . 3 2

A 8

■ln

 

 

 

 

 

 

или

 

 

 

2

, 2 5 Л

 

 

 

 

 

 

 

I

I . I , 3 2

* 8

 

 

 

 

Фк — Фс = (

In75----

 

 

 

< I I U 5 3 )

2 я й

8

 

 

 

 

Из (111.153) следует, что фильтрационное сопротивление ме­ жду поверхностями радиусов RK и Rc распадается на два после­ довательно соединенных сопротивления: одно непосредственно в окрестности скважины (первый член правой части уравнения (II 1.153)); другое — во внешней по отношению к окрестности скважины (второй член). Это и есть формула эквивалентных фильтрационных сопротивлений для указанной пространствен­ ной задачи. Структура первого члена уравнения соответствует структуре формулы А. К. Гиринского (III.141), пригодной, как показано в [61], для скважины, которая вскрыла пласт беско­ нечной толщины. Только в этой формуле вместо коэффициента 1,32 стоит коэффициент 1,6.

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

радиусов RK и 2.25Л

 

(приблизительно). Формула (III.152) при

б = 1 не приводится

к формуле Дюпюи

 

(хотя мало от нее от­

личается).

 

 

 

 

 

 

 

 

Вместе с тем необходимо отметить, что выражение

<

Г

( 1 + 0 , 8 7 5 8 ) Г

( 1 +

0 , 1

2

5

8

)

Ш

Г

( 2 - 0 , 8 7 5 8 ) Г

( 2 -

0 , 1

2

5

8

)

близко к нулю.

Используя это обстоятельство, а также уравнение (111.152), можно получить из формулы Маскета выражение, справедливое для всего диапазона б:

Q =

2 я й ( Ф к - Ф

с )

(III. 154)

1 , 3 2 й

 

где

IRc

 

 

 

 

 

5 = д / 1 — 6 +

б2.

(III. 155)

При этом предельный переход б

0 дает для второго члена

знаменателя (111.154) значение, мало отличающееся от соответ­ ствующего члена выражения (II 1.152) (вместо коэффициента 2,25 имеем 2,42).

ГЛАВА IV

ФИЛЬТРАЦИЯ ОДНОРОДНОЙ СЛАБОСЖИМАЕМОЙ ЖИДКОСТИ В ПЛАСТЕ

ОСНОВНЫЕ ПОНЯТИЯ ОБ УПРУГОМ РЕЖИМЕ ФИЛЬТРАЦИИ

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

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

Упругие свойства воды и нефти

Коэффициент сжимаемости жидкости рж можно определить из следующего соотношения:

где рж — плотность жидкости.

Заменяя рж через M/Vm, где М — масса жидкости, объем ее при данном давлении р, из (IV. 1) получим

dVx/Vx = —К dp.

(IV. 1)

Кж —

(IV.2)

Согласно [8], при 20 °С и атмосферном давлении для чистой негазированной воды имеем рв = 4,69* 10-10 Па-1.

Этот коэффициент, как утверждается в [8], практически не зависит от температуры, давления и характера процесса (адиа­ батического или изотермического).

Однако это утверждение нельзя считать точным. Исследова­ ния в литературе [81] показывают явную зависимость р„ как от температуры, так и от давления. При очень высоких давле­ ниях ( р « 44 МПа и температурах от 20 до 100 °С) р„ составляет (3,9 — 4,6) •10-10 Па-1, т. е. величины его довольно близки к зна­

чению, указанному выше в

[8]. С понижением давления р„ за­

кономерно увеличивается,

достигая при р = 7 МПа (4,9 —

— 5,10) •10-10 Па-1 в том же интервале температур. Если рв мо­ нотонным образом изменяется с изменением давления, то влия­ ние температуры на этот показатель не является монотонным и в «окрестности» температуры 50— 60 °С имеет явно выражен­ ные минимумы 3,9- 10~10 Па-1 = 44 МПа) и 5,0- 10-10 Па-1 (р = 7 МПа).

96

Экстраполяция этих данных до р ~ 0,1 МПа дает основание считать, что при нормальных атмосферных условиях рв может находиться в пределах (4,8 — 5,1) •10-10 Па-1, т. е. несколько выше цифры из источника [8].

Важно отметить, что при температурах порядка 130— 150°С рв « 6 •10_4 МПа- 1 (рис. IV.1).

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

дены

соответствующие

эм­

 

пирические

формулы,

учи­

ygg.nna1

тывающие

влияние

 

этих

 

 

факторов.

 

 

 

анализ

 

Обстоятельный

 

результатов этих исследова­

 

ний выполнен В. Н. Щелка-

 

чевым

[81]. Отметим

 

лишь,

 

что степень

минерализации

 

воды

оказывает

более

су­

 

щественное

влияние

 

на

рв.

 

Из примера,

приведенного

 

В. Н. Щелкачевым, следует,

 

что добавление солей в ве­

 

совых

количествах

 

70—

 

270 мг/л

к

чистой

воде

 

уменьшает

рв

соответствен­

 

но на

12 и 40 %> в то время

 

как насыщение

воды

газом

Рис. IV. 1. Коэффициент сжимаемости чистой

увеличивает

этот

показа­

воды

тель всего лишь на 2 —4% .

О распространенности такой концентрации свидетельствуют

также данные по ряду месторождений мира,

опубликованные

в известном труде А. И. Леворсена [34].

по крайней мере

Анализ этих материалов показывает, что

в 30 % случаях это содержание превосходит то минимальное зна­ чение его, которое было принято В. Н. Щелкачевым в расчетах (70 мг/л). Рассмотренные в [81] случаи достаточно характерны и дают основание считаться с фактором минерализации при оп­ ределении упругих свойств пластовых вод.

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

следующего вида:

 

рвг = Рг(1 + аП ,

(IV.3)

где рвг— коэффициент сжимаемости воды, содержащей раство­ ренный газ; Г — количество растворенного газа в м3/м3. По одним

7 Заказ № 283

97

данным [81] аж0,05, по другим — значительно меньше. Но в этом или другом случае расчеты показывают, что раство­ ренный газ увеличивает сжимаемость воды на очень небольшую величину (необходимо учесть, что величина Г по порядку соот­ ветствует единице и в редких случаях превышает 2—3. В этом заключается причина незначительной разницы между рвг и рв). Изменение объема воды при растворений в ней газа (метана, азота, углекислого газа) настолько незначительно, что его в рас­ четах можно не учитывать [81].

В работе Ш. Г. Гиматудинова Д10] показано, что с раствори­ мостью углеводородных газов, несмотря на ее небольшую вели­ чину, необходимо считаться в случаях газовых залежей с под­ стилающей водой при значительности их площади контакта и больших давлениях в пласте.

Коэффициент сжимаемости рн нефти зависит от состава пла­ стовой нефти, количества растворенного газа, абсолютного дав­ ления и температуры. Чем больше растворенного в нефти газа, тем больше коэффициент сжимаемости. Для дегазированной нефти рн колеблется от 4 •10-4 до 7 •10-4 МПа-1, что в общем соответствует коэффициенту сжимаемости для воды. При рас­ творении в нефтях (особенно легких) больших количеств газа происходит увеличение коэффициента сжимаемости нефти до 140 •10-4 МПа-1, что уже на порядок выше значения рн дегази­ рованной нефти.

Для

расчетов коэффициентов сжимаемости нефти с раство­

ренным

в ней газом существуют различные методы

[68, 81].

Перейдем к интегрированию уравнения (IV.2). С учетом

всего

вышесказанного выражение

(IV.2) следует

переписать

так:

 

 

 

 

dVJVx = —р(р,

Т) dp

(IV. 4)

подчеркнув тем самым, что коэффициент сжимаемости смеси (нефть и вода) — это функция давления и температуры. Коли­ чество минерального вещества и растворенного газа мы в эту функциональную зависимость не включаем, так как и количе­ ство растворенных минеральных веществ и газа зависит от дав­ ления. Что же касается температуры, то мы с самого начала оговоримся, что процесс в пласте является изотермическим. Как справедливо отмечает В. Н. Щелкачев [81], в большинстве за­ дач, связанных с исследованием разработки нефтяных пластов в условиях упругого режима, такое предположение вполне спра­ ведливо.

Примем, что при давлении р = р0 Кж = Кн<о. Тогда, инте­ грируя (IV.4), получим

Р

 

SР ж ф ^

(IV.5)

98

или

У* = У* о exp 1^- S! рж (6) rfsJ.

(IV-6)

Для раскрытия интеграла необходимо знать аналитический или графический вид зависимости р®(р). Задачу можно сущест­ венно упростить, если учесть, во-первых, малость величины Р® и, во-вторых, диапазон изменения давлений |р — ро|. Действи­ тельно, возьмем р® = 10-10-4 МПа- 1 и |р— ро| = 4 -10 МПа, тогда

Рж I Р — Ро! = 4 •10-2.

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

V* = V*o(l - SРж(6)л).

(IV.7)

А так как величину Р®(р) в рассмотренном диапазоне изме­ нения давления практически можно считать постоянной [рав­ ной, например, Р®(ро)], то предыдущее равенство перепишется так:

о =

________ 1

У ж

- У

ж

о

________ 1

\ & У Ж \

-

Рж

У ж о

р Ра

 

Уж о

I Ар I

 

 

1

_ Рж

Рж О

,

,

1 #

Арж

(IV.8)

 

Рж

Р

РО

 

 

Рж о

Др

 

Здесь мы уже убрали символ (р) при Р®, подчеркнув тем самым постоянство (конечно, в приближенном понимании) ко­ эффициента сжимаемости.

Упругие свойства коллекторов нефти

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

Если р — пластовое давление, то сжимающее скелет пласта давление можно в первом приближении принять равным

рск = о — р.

(IV. 9)

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

7*

99

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

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

Остановимся на модели В. Н. Щелкачева [81], которая пред­ ставляет собой идеализированную породу, состоящую из оди­ наковых упругих шаров диаметра D в кубической упаковке, за­ жатой двумя горизонтальными плоскостями.

Пусть к горизонтальной плоскости сверху приложено равно­ мерно распределенное давление (горное давление). Усилие, при­

 

ходящееся на одну вертикальную ко­

 

лонну шаров, равно F. Площадь прило­

 

жения этой силы,

очевидно,

равна D2.

 

Обозначим через р

давление

жидкости,

 

заполняющей пористую среду.

 

 

 

При этом происходит деформация ша­

 

ров (рис. IV.2). Вычислим величину си­

 

лы, действующей на одну вертикальную

 

колонну шаров. Если бы давление в по­

 

ристой среде равнялось нулю, то усилие

 

на

колонну равнялось бы F.

Давление

 

р в пористой среде должно противодей­

 

ствовать сжимающему усилию и дейст­

Рис. IV.2. Колонна недефор-

вует снизу не на всю сжимающую гори­

мированных (а) н сжатых

зонтальную плоскость, а только на часть

(б) упругих шаров

ее,

а именно

на

площадь

D2nd2l4

 

(d — диаметр

отпечатка — контура дав­

ления, переменная величина, зависящая от результирующего усилия).

Однако в области упругих деформаций d<^D, поэтому мо­ жно считать, что давление действует на всю площадку D2 и ре­ зультирующее сжимающее усилие будет приближенно равно

F0 = F - pD2 = D2( o - р).

(IV. 10)

Согласно разработанной Г. Гертцем теории сжатия двух со­ прикасающихся упругих шаров, сближение а двух соседних ша­ ров определяется из следующего уравнения:

_а_

-|3/

9 ^ ( 1 - v2)2

(IV. 11)

D

| /

£>4£2

 

где D — диаметр шара; F0— усилие; Е — модуль упругости ма­ териала шаров; v — коэффициент Пуассона для материала шаров.

100