ПРОЦЕССЫ НЕФТЕГАЗОДОБЫЧИ
.pdfvk.com/club152685050 | vk.com/id446425943
Глава 5 |
317 |
проницаемость и мощность пласта, R – радиус зоны дренирования сква- жины, rc – радиус скважины, S – скин-фактор (предполагается, что P = PR при r = R).
Прямое вычисление функции H (Pr ) представляет собой нелегкую задачу, основные трудности которой связаны с заданием относительных фазовых проницаемостей нефти и газа. Однако вид этой функции можно определить и без вычислений, используя только отмеченный Вогелем факт существования универсальной зависимости.
Действительно, из (5.40) следует |
|
H (Pcr ) |
|
|||||||||||
|
Q |
|
|
|
= 1 − |
|
|
|||||||
|
|
|
|
|
|
|
|
|
. |
(5.41) |
||||
|
KH (R |
Rr |
) |
H (P |
|
) |
||||||||
|
|
|
|
|
|
|
|
Rr |
|
|
|
|||
Как показал Вогель, правая часть (5.41) должна зависеть не от Pc |
||||||||||||||
и PR в отдельности, а только от их отношения Pc PR . Это возможно, если |
||||||||||||||
функция H (Pr ) представляет собой степенную функцию: |
|
|||||||||||||
|
H (P ) = Pнас H |
|
Pm |
, |
|
|
(5.42) |
|||||||
|
r |
|
|
|
|
|
|
|
1 r |
|
|
|
|
|
где H1 – некоторая постоянная. |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
Тогда (5.41) можно переписать в виде, подобном (5.33): |
|
|||||||||||||
~ |
= 1 |
~m |
, |
|
|
|
(5.43) |
|||||||
|
Q |
− P |
|
|
|
|
|
|
||||||
если принять |
|
|
|
|
|
|
|
Pm . |
|
|
|
|
||
|
Q = KPнас H |
|
|
|
|
|
(5.44) |
|||||||
|
m |
|
|
|
|
|
1 |
Rr |
|
|
|
|
||
Значение показателя степени |
m можно определить из условия ра- |
|||||||||||||
венства производных правых частей (5.33) и (5.43) в точке |
Pr = 1. Легко |
|||||||||||||
показать, что из этого условия следует m = 1,8. |
|
|
|
|||||||||||
Коэффициент H1 в уравнении (5.42) |
связан со значением относи- |
тельной фазовой проницаемости нефти на границе зоны разгазирования. |
||||||||||
Действительно, из определения функции H (Pr ) |
как интеграла с пе- |
|||||||||
ременным верхним пределом Pr |
следует |
|
|
|
|
|
|
|||
|
|
|
fн (sг ) |
|
|
|
|
|
|
|
dH |
|
|
|
|
|
|
|
|||
dP |
|
= Pнас |
|
|
|
|
|
. |
(5.45) |
|
|
|
|
|
|
||||||
|
|
|
β (P )ϕн (P ) |
|
|
|
|
|||
r |
P |
=1−0 |
|
r |
r |
|
P |
=1 |
|
|
|
r |
|
|
|
|
|
|
|
||
|
|
|
|
|
|
|
|
r |
|
|
Поскольку свободный газ становится подвижным только при дости- |
жении некоторой критической газонасыщенности sг0 , то на границе разга-
зирования происходит скачок газонасыщенности от нуля до sг0 . При этом
относительная фазовая проницаемость нефти скачком уменьшается от 1 до |
||||
значения fн |
(sг0 ). Поскольку β |
(1) = 1, ϕн (1) = 1, то из (6.45) следует |
||
|
dH |
|
|
= Pнас fн (sго ). |
|
dPr |
|
|
|
|
|
P =1−0 |
||
|
|
|||
|
|
|
r |
|
vk.com/club152685050 | vk.com/id446425943
Глава 5 |
|
|
|
319 |
|
нофазную фильтрацию к «укрупненной» скважине с радиусом r |
и забой- |
||||
ным давлением Pнас , поэтому здесь применима формула Дюпюи в виде |
|||||
|
B0µн0Q |
|
|
|
|
|
|
R |
|
|
|
PR − Pнас = 2π k h |
|
(5.47) |
|||
ln r |
. |
||||
|
|
|
|
|
(Предполагается, что изменением произведения объемного фактора нефти
на ее вязкость при P > Pнас можно пренебречь, считая β (Pr )ϕн (Pr ) = 1
при Pr > 1 [21].)
Для области разгазирования rc < r < r по аналогии с (5.40) имеем
H (1)− H (Pcr ) = |
B µ |
н |
0 |
Q |
r |
|
|
|
|
|
|
|
|||||||||||||
0 |
|
|
|
|
|
|
|
|
|
|
+ S . |
|
|
(5.48) |
|||||||||||
2π k h |
|
|
|
|
|
|
|||||||||||||||||||
|
ln r |
|
|
|
|
||||||||||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
c |
|
|
|
|
|
|||||
Сложив уравнения (5.47) и 5.48) и учитывая соотношение |
|
||||||||||||||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||||||
|
|
|
R |
|
|
|
r |
|
|
|
|
|
R |
|
, |
|
|
|
|
||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||||||||||||
|
|
ln r |
|
+ ln r |
|
|
= ln r |
|
|
|
|
|
|||||||||||||
|
|
|
|
|
|
c |
|
|
|
|
c |
|
|
|
|
|
|
||||||||
получим |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
Q = K(PR − Pнас )+ K[H (1)− H (Pcr )]. |
|
|
|
||||||||||||||||||||||
Принимая во внимание уравнения (5.42), (5.46) и обобщая на слу- |
|||||||||||||||||||||||||
чай Pc > Pнас , имеем |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
K(PR − Pc ), |
|
Pc > Pнас , |
|
|
|
|
|
|
|
|
|
|
|
||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
1,8 |
|
|
|
|
|||
Q = |
+ |
KPнас fн |
(sг0 ) |
|
− |
Pc |
|
|
|
|
|
< Pнас . |
(5.49) |
||||||||||||
|
|
|
|
|
|
||||||||||||||||||||
Qнас |
|
|
|
|
|
|
|
|
|
1 |
|
|
|
|
|
|
|
, P |
|
||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|||||||||||
|
|
|
|
1,8 |
|
|
|
|
|
|
|
|
|
|
|
|
|
c |
|
|
|||||
|
|
|
|
|
|
|
|
|
|
Pнас |
|
|
|
|
|||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
Легко видеть, что угол наклона индикаторной кривой, описываемой уравнением (5.49), скачком меняется в точке P = Pнас . Действительно,
dQ |
|
|
= −K , |
||
dPc |
|
|
|||
|
P = P +0 |
||||
|
|||||
|
|
|
|
c |
нас |
dQ |
|
|
|
|
= −Kfн (sг0 ). |
dPc |
|
|
|
|
|
|
P = P |
−0 |
|||
|
|||||
|
|
c |
|
нас |
|
Этот слом индикаторной кривой связан со скачком газонасыщенно- сти, который, как уже отмечалось, имеет место на границе зоны разгазиро- вания.
Следовательно, предположение о гладкости индикаторной кривой, которое используется при построении композитной кривой Вогеля [18], неверно.
Отметим, что это предположение ведет к завышению прироста деби- та нефти, достигаемого при снижении забойного давления от Pнас до нуля, в 1 fн (sг0 ) раз. Поскольку относительная фазовая проницаемость нефти в области малых газонасыщенностей меняется достаточно рез-
vk.com/club152685050 | vk.com/id446425943 |
|
|
|
|
|
|
|
320 |
Глава 5 |
|
|
|
|
|
|
ко |
fн (sг0 ) ≈ 0,8, то оценка прироста дебита нефти может быть завышена |
||||||
на 20%–30%. |
|
|
|
|
|
|
|
|
Полученные выше соотношения являются, по существу, уточнением |
||||||
формулы Фетковича [22], который предложил аппроксимировать функцию |
|||||||
|
α (P ) = |
fн |
|
|
|
|
|
|
r |
ϕ н β |
|
|
|
||
|
|
|
|
|
|||
на интервале 0 < Pr < 1 ( 0 < P < Pнас ) линейной функцией, проходящей че- |
|||||||
рез начало координат (см. рис. 5.18, прямая 1), что, как легко видеть, также |
|||||||
приводит к степенной зависимости H от P вида (5.42) с m = 2 . |
|
||||||
|
Согласно нашему подходу функцияr |
α (Pr ) |
аппроксимируется зави- |
||||
симостью |
|
|
|
|
|
|
|
|
α (P ) = fн |
(sг |
0 |
)Pm−1 , |
|
(5.50) |
|
|
r |
|
|
r |
|
|
|
которая может быть получена путем дифференцирования (5.42) с уче- |
|||||||
том (5.46). Кривая 2 на рис. 5.18 представляет собой график функции (5.50) |
|||||||
при fн (sг0 ) = 0,8 и m = 1,8. |
|
|
|
|
|
|
|
|
α(Pr) |
|
|
|
|
|
|
|
1 |
|
|
|
|
|
|
|
1 |
|
|
|
|
|
|
|
|
|
|
|
2 |
|
|
|
0 |
|
|
|
|
1 |
Pr |
|
Рис. 5.18. Аппроксимации функции α (Pr ): |
|
|||||
|
1 – аппроксимация Фетковича, 2 – формула (5.50) с |
fн (sг0 ) = 0,8 и m = 1,8 |
vk.com/club152685050 | vk.com/id446425943
322 |
Глава 5 |
нормальный и логнормальный законы. Некоторое теоретическое обосно- вание этому убеждению дает центральная предельная теорема, но часто нормальная функция распределения используется только в силу удобства и привычности.
Б. Мандельброт показал, что не менее универсальным законом рас- пределения является гиперболический [25]. Гиперболические (степенные) законы распределения являются ближайшими «родственниками» фракта- лов – с этим и связана их широкая распространенность (мы уже упоминали об этом в главе 1).
5.4.1. Формы представления и свойства гиперболических зависимостей
Случайная величина X называется гиберполически распределенной,
если
Р(Х ≥ х) = |
А |
, 0 < x < ∞ , |
(5.51) |
хα
где Р(Х ≥ х) – вероятность того, что X ≥ x , A, α – некоторые постоянные положительные величины.
Если случайная величина дискретна, то вместо (5.51) используется
соотношение |
|
|
|
|
|
|
|
Nr (Х ≥ х) = |
|
|
А |
, |
(5.52) |
||
|
|
|
|||||
|
|
|
|
хα |
|
||
где Nr (Х ≥ х) – общее число случаев, в которых (Х ≥ х). |
|||||||
Предположив, что дискретные значения X |
ранжированы в порядке |
||||||
убывания, получим |
|
|
|
|
|
|
|
Nr (Х ≥ х(r))= r , |
|
||||||
где x(r) – значения x , имеющие ранг r . |
|
|
|
|
|||
Тогда, обращая (5.52), получим |
|
|
|
|
|
|
|
x |
= |
B |
|
, |
|
(5.53) |
|
|
|
|
|||||
(r) |
|
r β |
|
|
|
|
|
|
|
|
|
|
|
где B = A1/α , β = 1/α .
Это соотношение позволяет оценить значение случайной величины по ее рангу, т. е. оказывается еще одним примером применения порядко- вых статистик (см. раздел 5.1).
Гиперболическое распределение, представленное в виде (5.51)–(5.52) называют законом Парето, в честь итальянского экономиста Вильфредо Парето, обнаружившего, что количество людей с доходом, превышающим некоторую величину x , уменьшается с ростом x гиперболически.
Дискретный закон распределения в форме (5.53) впервые ввел Джордж Ципф для описания частоты употребления в текстах слов различ-
vk.com/club152685050 | vk.com/id446425943
Глава 5 |
323 |
ной величины. Й. Корчак показал, что закону (5.53) подчиняется также распределение числа озер [25].
Согласно (5.51) при х → 0 Р(Х ≥ х)→ ∞ . Но эта расходимость не
существенна, поскольку реальные величины (доходы граждан, длина слов, размеры островов) всегда имеют ограничения снизу (и сверху). В этом смысле соотношения (5.51)–(5.53) имеют промежуточно-асимптотический характер.
Как уже отмечалось в главе 1, гиперболические зависимости мас- штабно-инвариантны. Для уточнения этого утверждения рассмотрим ус-
ловное распределение |
|
|
|
|
|
|
|||
g(x, x0 ) = P(X ≥ x |
|
X ≥ x0 ), |
|
||||||
|
|
||||||||
определяющее вероятность того, что X ≥ x |
при условии X ≥ x0 . Посколь- |
||||||||
ку, по правилу умножения вероятностей, |
|
|
|||||||
P(X ≥ x) = P(X ≥ x |
|
X ≥ x0 ) P(X ≥ x0 ), |
|||||||
|
|||||||||
то |
λ (x) |
|
|
|
|||||
g(x, x ) = |
, |
(5.54) |
|||||||
λ (x0 ) |
|||||||||
0 |
|
|
|
|
|||||
|
|
|
|
|
|
|
|
||
где λ(x) = P(X ≥ x). |
|
|
|
|
|
|
|||
Функция g(x, x0 ) при различных значениях |
x0 соответствует раз- |
личным уровням рассмотрения исследуемой системы. Так, если речь идет о доходах, то функция g(x, x0 ) определяет распределение доходов среди
населения, уровень жизни которого выше предела, определяемого величи- ной x0 .
Распределение доходов будет масштабно-инвариантным, если функ- ции g(x, x0 ) при разных x0 подобны друг другу, т. е. если g(x, x0 ) зависит
не от x и x0 в отдельности, а только от их безразмерной комбина- ции х / х0 :
|
|
|
|
|
|
|
|
|||
|
|
|
|
|
x |
|
|
|
|
|
|
|
|
|
|
|
|
||||
g(x, x0 ) = g1 |
. |
|
||||||||
Подставив сюда x0 = 1, имеем |
|
x0 |
|
|
|
|||||
|
λ(x) |
|
|
|
|
|
||||
|
g (x) = |
, |
|
|
|
|||||
|
|
|
|
|||||||
1 |
|
|
λ(1) |
|
|
|
||||
|
|
|
|
|
|
|
|
|
||
откуда |
λ (x / x |
) |
|
|
||||||
|
λ (x) |
|
|
|||||||
|
|
= |
|
λ (1)0 |
|
. |
(5.55) |
|||
|
λ (x0 ) |
|
|
|||||||
Прологарифмировав (5.55) и осуществив преобразование |
|
|||||||||
ln λ = U (z), z = ln x , |
|
|||||||||
получим |
|
|
|
|
|
|
|
|
||
U (z)− U (zo ) = U (z − z0 )− U (0), |
(5.65) |
где z0 = ln x0 .
vk.com/club152685050 | vk.com/id446425943
324 Глава 5
Переходя к пределу z0 → z , из (5.56) можно получить
U ′(z) = U ′(0) = const ,
что возможно, если только функция U (z) линейна:
U (z) = a + bz
или
ln λ(x) = a + b ln x ,
откуда и следует (5.51) с A = ea , α = −b .
Таким образом, гиперболическое распределение (и только оно)
удовлетворяет условию масштабной инвариантности (5.55). |
|
|||||||||||||||
Легко видеть, что условное распределение g(x, x0 ) имеет вид |
|
|||||||||||||||
|
|
|
|
|
|
|
|
|
|
1, |
x < x0 |
|
||||
g(x, x ) = Р(X |
≥ x |
|
x ≥ x ) |
= |
x |
−α |
. |
|
||||||||
|
|
|||||||||||||||
0 |
|
|
|
|
|
|
0 |
|
|
|
|
, x ≥ x0 |
|
|||
|
|
|
|
|
|
|
|
|
|
|
|
|
||||
|
|
|
|
|
|
|
|
|
|
|
|
|||||
|
|
|
|
|
|
|
|
|
|
x |
|
|
||||
|
|
|
|
|
|
|
|
|
|
|
0 |
|
|
|
||
Соответствующая функция плотности распределения определяется |
||||||||||||||||
как |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||
|
dg(x, x0 ) |
0, |
x < x0 |
|
|
|
|
|||||||||
f (x) ≡ − |
= |
α |
|
|
−(α +1) |
. |
(5.57) |
|||||||||
|
||||||||||||||||
|
dx |
|
|
|
|
|
x |
|
|
|
, |
x ≥ x0 |
|
|||
|
|
x |
|
|
|
|
||||||||||
|
|
|
x |
|
|
|
|
|||||||||
|
|
|
|
0 |
0 |
|
|
|
|
|
|
В какой-то мере универсализм гиперболических распределений можно объяснить тем, что они характерны для систем, образование кото- рых контролируется процессами кластеризации. А этот механизм доста- точно широко распространен в природе [26]. Так, в социологии давно из- вестен феномен, который можно выразить словами «успех порождает ус- пех». Часто употребляемые слова становятся все более употребительными, крупные города разрастаются быстрее, часто цитируемые статьи все чаще цитируются и т. д. Все это – примеры социальной кластеризации, феноме- на, которому Р. Мертон дал название «эффект Матфея», имея в виду биб- лейское изречение «имущему дается».
Гиперболическая функция описывает резко неоднородное, асиммет- ричное распределение. Покажем это, воспользовавшись координатами Ципфа.
Пусть (х(1), x(2), ...x(N ) ) – ранжированная в порядке убывания выбор-
ка значений аддитивной величины x , |
распределенной по закону (5.53), |
|
N – объем выборки, σ (n) – сумма первых n значений х(r) |
(n = 1, 2, ..., N ); |
|
σ 0 – сумма всех значений x : |
|
|
n |
σ 0 = σ (N ). |
|
σ (n) = ∑ X(r) , |
|
r =1
vk.com/club152685050 | vk.com/id446425943
|
|
Глава 5 |
325 |
|||
Введем безразмерные переменные |
|
|
|
|||
µ = |
|
σ (n) |
ν = |
n |
||
|
|
, |
|
. |
||
|
σ (N ) |
N |
Легко видеть, что величина µ представляет собой долю σ 0 , «накоп-
ленную» в результате n реализаций, а ν – соответствующая доля реализа- ций. Так, если говорить о доходах, то µ есть доля совокупных доходов,
принадлежащих ν -й части населения; если объем выборки велик, то сумму можно заменить интегралом и считать n >> 1. Тогда, в предположе- нии β < 1, получим
n |
|
|
|
B |
(n1− β − 1)≈ |
B |
|
|
σ = ∫ |
B |
dr = |
|
|
n1− β , |
|||
|
|
|
|
|
||||
1 r β |
|
1 |
− β |
|
1 − β |
|
||
откуда |
|
|
|
|
|
|
|
|
|
|
|
|
|
µ = ν 1− β . |
|
|
|
График этой зависимости при β = 0,85 |
(значение, характерное для |
распределения доходов) представлен на рис. 5.19 (кривая ACB ). Этот гра- фик (так называемая кривая Лоренца) наглядно показывает неравномер- ность распределения, описываемого гиперболическим законом: уже при малых значениях ν величина µ близка к единице. Так, при ν = 0,2
µ = 0,8 . На примере распределения доходов это означает, что всего лишь
20% населения получают 80% доходов, в то время как остальные 80% имеют всего лишь 20% доходов. В более общем виде это правило, назы- ваемое принципом Парето, формулируется так: «В больших системах 80% случаев вызываются 20% причин и наоборот». Следует отметить, что граничные значения 80% и 20% достаточно условные, поскольку при дру- гих значениях β это может быть 90% и 10% или 70% и 30% и т. д.
Легко видеть, что |
абсолютно равномерное распределение доходов |
описывается прямой АВ |
(β = 0). Чем сильнее кривая Лоренца отклоняет- |
ся от прямой АВ, тем больше неравномерность распределения, поэтому мерой неоднородности может служить величина
|
L = SАСВА , |
|
|
|
|
|
SАВЕ |
|
|
|
|
представляющая собой отношение площади криволинейной фигуры ACBA |
|||||
к площади треугольника ABE . Эта величина называется коэффициентом |
|||||
Лоренца (или Джини [27]). Поскольку SАВЕ = 0,5 и |
|
||||
1 |
ν 1− β dν − SАВЕ = |
|
1 |
− 1 , |
|
SАСВА = ∫ |
|
||||
2 |
− β |
||||
0 |
|
2 |
|||
|
|
|
|
vk.com/club152685050 | vk.com/id446425943 |
|
|
|||
326 |
|
Глава 5 |
|
|
|
то |
|
|
β |
|
|
|
|
L = |
|
|
|
|
|
2 − β . |
|
|
|
Очевидно, что 0 < L < 1. |
|
|
|
||
µ |
|
|
|
|
B |
0,8 |
C |
|
|
|
|
|
|
|
|
|
|
0,6 |
|
|
|
|
|
0,4 |
|
|
|
|
|
0,2 |
|
|
|
|
|
A |
|
|
|
|
E |
|
|
|
|
|
|
0 |
0,2 |
0,4 |
0,6 |
0,8 |
ν |
|
|
Рис. 5.19. Кривая Лоренца |
|
|
Следует подчеркнуть, что коэффициент Лоренца определяется по интегрированным данным, поэтому эта величина устойчива к погрешности данных, что является общим свойством всех интегральных методов (см.
раздел 2.1.3).
Ограничение β < 1, принятое нами, не является обязательным. Соот- ношения, подобные полученным выше, могут быть выведены и при β > 1, если считать выборку ограниченной (N < ∞). Правда, полученные при этом формулы будут несколько сложнее.
5.4.2. Закон Парето в оценке запасов углеводородов
Начиная с работы Крига [28], было принято считать, что распределе- ние запасов минеральных богатств, включая нефть и газ, подчиняется лог- нормальному закону. Однако в 1962 г. Б. Мандельброт показал, что это