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

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

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

vk.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

=10

 

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 =10

 

 

 

 

 

r

 

vk.com/club152685050 | vk.com/id446425943

318 Глава 5

С другой стороны, дифференцирование (5.24) дает

dH

 

 

= Pнас H1m .

dPr

 

P =1−0

 

 

 

 

 

 

r

 

Из двух соотношений следует

 

 

 

 

 

H1 =

fH (sг0 )

.

 

(5.46)

 

 

 

 

 

 

 

 

 

 

 

m

~

~ ~

 

 

 

 

 

 

 

Итак, Вогелю следовало искать зависимость Q

= Q(P) не в виде по-

~

~

~2

с двумя эмпирическими коэффициентами a и b , а

линома Q

= 1 aP bP

 

~

 

~m

с одним неизвестным параметром m. Зависи-

в виде функции Q

= 1 p

мость, предложенная Вогелем, не соответствует физическому содержанию

задачи. Следует, однако, отметить, что сами по себе эти две функции весь- ма близки друг к другу (см. рис. 5.17).

~

P

0,8

0,6

0,4

0,2

0

0,2

0,4

0,6

0,8

~

 

 

 

 

 

Q

Рис. 5.17. Графики универсальных индикаторных кривых, полученных по разным формулам:

кривая Вогеля, уравнение (5.33)

кривая, полученная по уравнению (5.43) при m = 1,8

Атеперь более внимательно проанализируем структуру композитной индикаторной кривой. Пусть PR > Pнас и r радиус зоны разгазирования

вблизи скважины. Движение в области r > 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

)Pm1 ,

 

(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

Глава 5

321

Как видим, основное отличие между этими зависимостями состоит в том, что аппроксимация Фетковича не учитывает скачка функции α (Pr ) в

точке Pr = 1. Вместо этого снижение фазовой проницаемости нефти в зоне

разгазирования предлагалось учитывать путем введения в уравнения при- тока псевоскина [22]. Предложенный нами подход является более естест- венным.

Уравнение (5.49) можно переписать в виде

Q = K(PR Pc1),

где

 

Pc ,

 

 

 

 

 

 

Pc1

=

 

1

 

Pнас

 

 

 

 

 

 

 

 

 

 

Pc > Pнас ,

 

 

fн (sг

0

)

 

P

1,8

 

 

1

 

c

 

 

 

 

 

1,8

 

 

 

 

 

 

 

 

Pнас

, Pc < Pнас .

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

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

5.4. Гиперболические законы распределения

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

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

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 г. Б. Мандельброт показал, что это