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

ПРОЦЕССЫ НЕФТЕГАЗОДОБЫЧИ

.pdf
Скачиваний:
24
Добавлен:
24.08.2019
Размер:
8.3 Mб
Скачать

vk.com/club152685050 | vk.com/id446425943

Глава 5

307

Пример. Обработка данных исследования кернов одного из пластов Нефтеюганского региона по описанным выше алгоритмам привела к зави- симости

sно

=

 

 

sнн

,

(5.29)

 

+ k

0,15 1/ 2

 

1

sнн

 

 

где проницаемость измеряется в миллидарси. В этом случае, действитель-

но, σ = 0,15 < 1.

На рис. 5.10 приведены графики функций β = β (k), полученные с использованием (5.29), при sнн = 0,5 , sнн = 0,65 и sнн = 0,85 (кривые 1, 2 и 3 соответственно). Как видим, большой разброс значений коэффициента вытеснения при фиксированных значениях проницаемости, наблюдаю- щийся при проведении лабораторных исследований, можно объяснить влиянием начальной нефтенасыщенности.

5.3. Метод асимптотических координат

Положим, что имеется некоторая величина F , зависящая от двух па- раметров p и q . Пусть в условиях эксперимента задавались определенные значения параметра q = q1, q2 , q3,..., qn и определялась зависимость F от p при фиксированных q . В том случае, когда вид полученных кривых в плоскости (p, F ) носит качественно сходный характер, часто удается по- добрать специальные координаты, с помощью которых исследуемую сложную двумерную поверхность F = F(p, q) удается описать с помощью нескольких более простых плоских кривых (при этом семейство кривых в плоскости (р F ), соответствующих различным значениям q , сжимается

водну универсальную кривую).

Вкачестве примера рассмотрим экспериментальные зависимости де- бита жидкости Q от расхода газа V и диаметра подъемника d , получен-

ные [15] на лабораторной установке, моделирующей работу газлифтной скважины (см. рис. 5.11).

Перейдем к переменным

 

 

 

 

 

 

 

 

 

 

x =

V V0 (d )

 

 

 

 

 

 

 

 

 

 

,

(5.30)

 

 

 

V

(d )V

 

(d )

 

 

 

 

m

0

 

 

 

 

 

 

y =

 

Q

,

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Q (d )

 

 

где V0 (d )

расход

газа,

 

 

 

m

 

 

 

 

при котором

 

начинается подъем

жидкости

(Q 0 при

V V0 ),

Vm (d )

расход газа, соответствующий максималь-

ному дебиту жидкости Qm (d ).

 

 

 

 

 

 

 

Тогда все пять кривых на рис. 5.11 можно представить в виде одной универсальной зависимости y = f (x), показанной на рис. 5.12. Характер-

vk.com/club152685050 | vk.com/id446425943

 

 

 

 

 

308

 

 

Глава 5

 

 

 

 

 

ные величины V0 (d ), Vm (d ), Qm (d ) в зависимости от диаметра подъемника

приведены на рис. 5.13. Кривая

y = f (x)

может быть аппроксимирована

аналитической зависимостью

 

 

 

 

0 ≤ x ≤ 1;

 

 

 

y = xα exp (1 xα ),

0,9;

.

(5.31)

 

α =

 

x > 1

 

 

 

 

 

0,5;

 

 

 

Q, 10–3 м3/с

 

 

 

 

 

 

 

 

10

 

 

 

 

 

 

 

 

 

 

 

 

 

1

 

 

 

 

 

5

 

 

 

 

 

 

 

 

 

0

25

50

 

75

 

100

V, 10–3 м3/с

 

Рис. 5.11. Зависимости Q = Q (ν ) для подъемников различного диаметра ● d = 0,100 м, d = 0,075 м, d = 0,063 м, d = 0,050 м,

d = 0,038 м, 1 – восстановленная по трем точкам зависимость

y

 

 

 

 

0,75

 

 

 

 

0,50

 

 

 

 

0,25

 

 

 

 

0

1

2

3

x

Рис. 5.12. Нормированная зависимость дебита жидкости от расхода газа

d = 0,100 м, d = 0,075 м, d = 0,063 м, d = 0,050 м, d = 0,038 м

vk.com/club152685050 | vk.com/id446425943

Глава 5

309

/с

 

 

 

 

 

/с

3 м

 

 

 

 

 

 

 

 

 

3 м

–3

 

 

 

 

 

 

 

 

–3

, 10

 

 

 

 

 

 

 

 

 

, 10

m

 

 

 

 

 

 

 

0

V

 

 

 

 

 

 

 

 

V

40

 

 

 

 

 

 

 

 

 

 

10

 

 

 

 

 

 

 

2

 

 

 

 

 

 

 

 

 

 

 

 

 

 

20

 

 

 

 

 

 

 

 

 

 

5

 

 

 

 

 

 

 

 

1

 

 

 

 

 

 

 

 

 

 

 

 

 

0

 

 

 

 

 

 

 

 

 

 

0

 

 

 

 

 

 

 

 

 

 

0,025

0,050

0,075 d, м

 

 

 

 

 

 

а)

 

 

 

 

 

10

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

м3/с

 

 

 

 

 

 

 

–3

 

 

 

 

 

 

 

10

 

 

 

 

 

 

 

,

 

 

 

 

 

 

 

 

 

 

 

 

m

 

 

 

 

 

 

 

Q

 

 

 

 

 

 

 

5

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

3

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

0

0,025

0,050

0,075

d, м

 

б)

 

 

Рис. 5.13. Зависимости характеристических величин от диаметра подъемника

1 – V0 (d ), 2 – Vm (d ), 3 – Qm (d )

Теперь для любого подъемника нам достаточно знать значения трех величин V0 (d ), Vm (d ), Qm (d ), чтобы рассчитать зависимость Q(V ) по формуле

 

 

 

V0

α

 

 

 

 

V0

α

 

Q = Q

V

 

 

V

 

 

 

 

 

 

 

.

 

 

 

 

m V

V

 

 

 

V

V

 

 

 

 

 

m

0

 

 

 

 

m

0

 

 

 

vk.com/club152685050 | vk.com/id446425943

310

Глава 5

Описанный способ нормирования кривых родственен простому обезразмериванию и широко применяется при обработке эксперименталь-

ных данных. В работе [16] он назван методом асимптотических коорди-

нат, поскольку вид нормирующих преобразований устанавливается путем изучения поведения кривых в некоторых предельных случаях ( p 0 ,

p → ∞ ,

F

0

и т. д.). Так, в нашем случае характеристические точки

 

p

 

 

 

 

определяются условиями Q = 0 и

dQ

= 0.

 

 

 

 

dV

 

Представление исследуемой двумерной поверхности с помощью

плоских кривых (в нашем примере представление зависимо- сти Q = Q(V , d ) набором кривых, приведенных на рис. 5.12–5.13) облегчает

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

5.3.1. Восстановление характеристик газлифтных скважин

Предположим, что у нас имеется некоторое множество газлифтных скважин, работающих в примерно одинаковых условиях, часть из которых была тщательно исследована на различных режимах работы с получением регулировочных кривых Q = Q(V ). Представив эти кривые в асимптотиче- ских координатах вида (5.30), можно определить вид зависимости y = f (x)

и считать, что регулировочные кривые остальных, неисследованных, сква- жин в координатах (x, y) имеют такой же вид.

Это позволяет существенно упростить исследование второй группы скважин: вместо того чтобы проводить полномасштабные эксперименты, на каждый из них достаточно сделать замеры Q и V при трех различных режимах закачки газа. Тогда неизвестные значения характеристических величин V0 , Vm , Qm определяются путем решения относительно них сис-

темы из трех уравнений

 

Vi

V0

 

=

Qi

 

 

 

 

,

(5.32)

 

V

Q

f V

 

 

m

0

 

 

m

 

 

где Qi , Vi дебит жидкости и расход газа на i -м режиме работы скважины

(i = 1, 2, 3 ).

vk.com/club152685050 | vk.com/id446425943

Глава 5

 

311

Определив V0 , Vm , Qm , можно

восстановить всю зависи-

мость Q = Q(V ) по формуле

 

 

 

V0

 

 

V

 

 

V

Q = Qm f V

.

m

0

 

 

В частности, оптимальный расход газа, Vоп , определяемый условием

 

 

 

d Q

 

 

dQ

(Vоп )

 

Qоп

 

 

 

 

 

= 0 или

 

 

 

 

=

 

,

 

 

 

 

 

 

 

 

 

 

 

 

dV V

 

 

dV

 

Vоп

 

находится из выражения

(Vm Vo )xоп ,

 

 

 

 

 

 

Vоп = V0 +

 

 

где xоп есть корень уравнения

 

 

 

 

 

 

 

 

 

 

 

 

 

f

=

~

1

,

 

 

 

(5.33)

 

 

 

 

 

 

 

+ x

 

 

 

 

 

 

 

 

f

V0

 

 

 

 

 

 

~

V0

 

 

 

 

 

 

 

 

= Vоп .

 

 

V0 =

 

, Qоп дебит жидкости при V

 

 

Vm V0

 

 

 

 

 

 

 

 

 

 

 

(x) немонотонна, то единствен-

 

Заметим, что поскольку функция f

ное решение системы (5.32) можно получить только при условии, что про каждую точку (Vi , Qi ) заранее известно, на какой ветви зависимо- сти Q(V ) левой или правой она находится.

Из вида кривых Q = Q(V ) (см., например, рис. 5.12) ясно, что без ап- риорной информации, представленной в виде универсальной функции, восстановить зависимость Q(V ) всего лишь по трем точкам невозможно.

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

Универсальная кривая y = f (x) дает формализованное представле- ние априорной информации о результатах исследований, проведенных ра- нее на других скважинах. Надежное восстановление зависимости Q = Q(V ) по малому числу замеров возможно потому, что при учете априорной ин- формации происходит «обогащение» информации о данной скважине она пополняется и уточняется за счет предыдущего опыта исследования подобных скважин.

Для апробации предложенного алгоритма попытаемся восстановить зависимость Q = Q(V ) для d = 0,100 м по трем экспериментальным точкам

(27,1; 4,67), (79,5; 7,83),и (100,7; 7,67) обведенным на рис. 5.11 кружками.

Численное решение системы в этом случае дает Qm = 7,8 103 м3/с,

V = 12 103

м3/с. Зависимость Q = Q(V ), определяемая при найденных

0

 

vk.com/club152685050 | vk.com/id446425943

 

 

 

312

 

 

 

Глава 5

 

 

 

значениях параметров, представлена на рис. 5.11 штриховой линией. Оп-

тимальный режим, получаемый по уравнению 5.33, характеризуется вели-

чинами V

 

= 29,6 103

м3/с, Q

= 5 103

м3/с.

 

 

оп

 

оn

 

 

 

 

5.3.2. Расчет притока нефти к скважине с забойным давлением

ниже давления насыщения

 

 

 

В настоящее время при расчете индикаторных кривых (зависимостей

дебита нефти Q от забойного давления Pc ) для скважин, работающих при

забойном давлении ниже давления насыщения, широко используются ре-

зультаты исследования Вогеля [17], который путем численного решения

уравнений движения газированной нефти при разных значениях параметра

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

которых представлен на рис. 5.14. Эти кривые соответствуют различным

стадиям истощения пласта и характеризуются двумя параметрами пла-

стовым давлением PR i

(определяемым по значению Pc при Q = 0 ) и мак-

симальным дебитом Qm i , достигаемым при Pc = 0 ( i

номер кривой в се-

мействе). При расчете каждой серии кривых начальное пластовое давление

принималось равным давлению насыщения ( PR1 = Pнас ).

 

 

Pс

 

 

 

 

 

 

PR1

 

 

 

 

 

 

PR2

 

 

 

 

 

 

PR3

 

 

 

 

 

 

PR4

 

 

 

 

 

 

PR5

 

 

 

 

 

 

 

Qm5 Qm4

Qm3

Qm2

Qm1

Q

Рис. 5.14. Индикаторные кривые при различных значениях пластового давления

vk.com/club152685050 | vk.com/id446425943

Глава 5

313

Переходя к асимптотическим координатам

~

 

Q

 

~

 

P

 

Q

=

 

,

P

=

c

,

Q

P

 

 

 

 

 

 

 

 

m

 

 

 

R

 

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

эту зависимость в виде полинома второй степени и пришел к уравнению

~

~

~2

.

(5.33)

Q

= 1 0,2P

0,8P

Путем многочисленных расчетов им было показано, что уравне- ние (5.33) действительно универсально: оно применимо для пластов с са- мыми различными фильтрационными характеристиками и PVT-свойствами флюидов. Ошибка, допускаемая при применении уравнения Вогеля, в среднем не превышает 10%. Поскольку это уравнение не содержит в явном

виде значения газового фактора, оно применимо и для обводненных сква- жин, если под Q понимать дебит жидкости [18].

При разработке месторождений методом заводнения пластовое дав- ление, как правило, поддерживается выше давления насыщения, т. е. сква-

жины работают (при Pc < Pнас ) в режиме локального разгазирования,

когда газ в свободном виде выделяется только в некоторой области вблизи скважины (размеры этой области обычно не превышают несколько десят-

ков сантиметров).

Pc

PR

0,8

0,6

0,4

0,2

0

0,2

0,4

0,6

0,8

Q

 

 

 

 

 

 

 

 

 

 

 

Qm

Рис. 5.15. Индикаторные кривые в асимптотических координатах

vk.com/club152685050 | vk.com/id446425943

 

 

314

Глава 5

 

 

 

PС

 

 

 

A

 

 

 

PR

 

 

 

B

 

 

 

Pнас

 

 

 

 

C

 

 

Qнас

Qm

Q

 

Рис. 5.16. Композитная индикаторная кривая

 

Для определения дебита скважины в условиях локального разгазиро- вания предложено [18] использовать композитную индикаторную кривую (рис. 5.16), при построении которой используют следующие предположе-

ния:

 

- при Pc < Pнас (участок AB на рис. 5.16) зависимость Q от Pc

прямоли-

нейна:

 

Q = K(PR Pc ),

(5.34)

где K коэффициент продуктивности скважины в отсутствие газа;

-при 0 < Pc < Pнас отрезок индикаторной кривой (участок BC ) подобен

кривой Вогеля, т. е. описывается уравнением (5.33) с

 

 

~

 

Q Qнас

 

 

~

 

Pc

 

 

 

 

 

 

Q =

 

 

 

,

P

=

 

 

 

 

 

 

 

 

Q Qнас

 

Pнас

 

 

 

 

 

 

 

 

m

 

 

 

 

 

 

 

 

 

 

 

или

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Q Qнас

 

 

 

Pc

 

 

 

 

Pc

2

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

= 1 0,2

 

 

0,8

 

 

 

,

(5.35)

 

Q

Qнас

 

Pнас

 

Pнас

 

m

 

 

 

 

 

 

 

 

 

 

 

 

 

где

Qнас = K(PR Pнас );

vk.com/club152685050 | vk.com/id446425943

Глава 5

315

- кривая BC касается прямой AB, т. е. углы их наклона в точке B равны:

dQ

(Pнас + 0) =

dQ

(Pнас 0)

 

dP

 

 

 

dP

 

 

c

 

 

 

c

 

 

или

 

 

 

 

 

 

K = − 1,8

(Q

Qнас ).

(5.36)

 

 

Pнас

 

m

 

 

 

 

 

 

 

С учетом (5.34)–(5.36)

K(PR Pc ),

 

 

 

Q =

+

KPнас

Qнас

 

 

 

 

1,8

 

 

 

 

получим окончательно

Pc > Pнас

 

 

 

 

2

 

 

 

 

 

 

 

 

 

Pc

Pc

 

, Pc < Pнас . (5.37)

 

 

 

 

 

1

0,2

 

 

0,8

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Pнас

 

Pнас

 

 

 

 

 

 

 

 

 

 

Результаты Вогеля с аналитической точки зрения

Итак, Вогель обнаружил замечательный факт существования уни- версальной формы представления индикаторных кривых. Из предыдущего ясно, что при дальнейшем использовании и обобщении этого результата на случай PR > Pнас были приняты следующие предположения:

1.Индикаторная кривая при Pc < Pнас описывается полиномом второй степени.

2.Углы наклона прямолинейного и криволинейного участков индикатор- ной кривой в точке P = Pнас равны.

Рассмотрим с позиций теории фильтрации газированной жидкости, насколько обоснованы эти предположения.

Уравнение стационарной совместной фильтрации нефти и газа мож- но записать в виде [19, 20]

 

 

P

 

ρг (P) fг (Sг ) µн (P)

 

 

 

 

 

=

 

 

 

 

 

B(P),

(5.38)

G 1

 

 

 

 

 

 

 

 

ρ гс fн (Sг ) µ

г (P)

 

 

Pнас

 

 

 

где G газовый фактор, м3/м3; ρг плотность газа при данных условиях,

кг/м3; ρгс плотность газа в нормальных условиях, кг/м3; fн , fг относи-

тельные фазовые проницаемости нефти и газа; Sг насыщенность газом,

µн (P), µ г (P) вязкость нефти и газа при пластовой температуре и давле-

нии P , Па с; B(P) объемный коэффициент нефти. Соотношение (5.38) можно переписать в безразмерном виде

g(sг ) = η(Pr ),

vk.com/club152685050 | vk.com/id446425943

316 Глава 5

где

 

 

 

 

 

 

 

 

 

f (sг ) µно

 

 

 

 

 

 

 

 

g(sг ) =

 

 

 

 

 

 

 

 

 

 

r

µ

го ,

 

 

 

 

 

 

 

fн (sг )

 

 

 

 

 

 

η (P ) =

 

 

G(1 Pr )

 

ϕ r

(Pr )

,

 

 

 

 

 

B0

 

 

 

 

 

 

 

 

r

 

β (Pr )λ (Pr ) ϕ н

(Pн )

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

P = P

 

, ϕ г (P )

=

г (Pr ),

 

ϕн (P ) = µ

н (Pr )

,

r

Pнас

 

 

r

 

 

 

µ го

 

 

r

 

µно

 

 

(Pr )

 

 

 

 

 

 

 

 

 

 

λ =

ρ г

,

 

β (P ) = B(Pr ),

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

ρ гс

 

 

 

r

 

Bо

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

µно , µ го вязкости нефти и газа при давлении насыщения, Па с, Bо объ-

емный фактор нефти при p = pнас .

 

 

 

 

 

 

 

 

 

Обращая функцию g(sг ), получим из (6.39)

 

 

 

 

 

 

sг = ψ (P ) = g1[η(P )].

 

 

 

 

 

 

 

 

 

 

 

r

 

 

 

 

r

 

 

 

 

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

в режиме стационарной фильтрации газированной

жидкости насыщенность газом sг в каждой точке пласта однозначным об- разом связана с давлением. Это, в свою очередь, позволяет связать значе- ния фазовых проницаемостей с давлением:

 

 

 

 

fг

(sг ) = fг

[ψ (Pr )],

 

 

 

 

fн

(sг ) = fн

[ψ (Pr )],

и линеаризировать уравнения фильтрации путем введения функции

 

 

 

 

H (Pr ) = Pнас Prα (Pr )dPr ,

 

 

fн [ψ (Pr )]

 

 

0

где α (Pr ) =

 

.

 

 

ϕ

н (Pr )β (Pr )

 

 

 

 

 

 

Эта функция называется псевдодавлением или (в отечественной ли-

тературе) функцией Христиановича.

 

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

быть получены из уравнений однофазной фильтрации путем простой заме-

ны давления P на псевдодавление H

(P). В частности, формула Дюпюи

для радиального притока жидкости к скважине примет вид [19]:

 

 

 

 

 

 

 

 

 

Q = K[H (PRr )H (Pc r )] ,

(5.40)

где

K =

 

 

2π k h

 

 

, PRr

приведенное пластовое

давле-

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

R

 

+ S

 

 

 

 

 

 

 

 

 

 

 

 

 

 

B0µно ln r

 

 

 

 

 

 

 

 

 

c

 

 

 

 

 

 

 

 

PR

 

 

 

 

 

 

 

 

 

 

 

 

 

 

приведенное значение забойного давления, k и h

 

 

 

 

 

ние PRr =

 

, Pсr

 

 

 

Pнас