ПРОЦЕССЫ НЕФТЕГАЗОДОБЫЧИ
.pdfvk.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
Глава 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 10−3 м3/с,
V = 12 10−3 |
м3/с. Зависимость Q = Q(V ), определяемая при найденных |
0 |
|
vk.com/club152685050 | vk.com/id446425943 |
|
|
|
||||
312 |
|
|
|
Глава 5 |
|
|
|
значениях параметров, представлена на рис. 5.11 штриховой линией. Оп- |
|||||||
тимальный режим, получаемый по уравнению 5.33, характеризуется вели- |
|||||||
чинами V |
|
= 29,6 10−3 |
м3/с, Q |
= 5 10−3 |
м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 |
|
|
|
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
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 ) = g−1[η(P )]. |
|
|
|
|||||||||||
|
|
|
|
|
|
|
|
r |
|
|
|
|
r |
|
|
|
|
Следовательно, |
в режиме стационарной фильтрации газированной |
жидкости насыщенность газом sг в каждой точке пласта однозначным об- разом связана с давлением. Это, в свою очередь, позволяет связать значе- ния фазовых проницаемостей с давлением:
|
|
|
|
fг |
(sг ) = fг |
[ψ (Pr )], |
|
|
|
|
fн |
(sг ) = fн |
[ψ (Pr )], |
и линеаризировать уравнения фильтрации путем введения функции |
||||||
|
|
|
|
H (Pr ) = Pнас P∫rα (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нас |
|
|
|
|
|
|
|