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

книги / Справочное руководство по проектированию разработки и эксплуатации нефтяных месторождений. Проектирование разработки

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

в результате решения интегрального уравнения

относительно функции

dx/dt:

 

t

 

| ф ( / - т ) ^ Л = М > ,

(XVI 11.23)

О

 

полученного из уравнения (XVIII. 16). После определения функции dx/dt в за­ висимости от t функцию х (/) находим при условии х (0) = 0.

Продолжительность безводного периода tx вычисляют как корень трансцен­ дентного уравнения

* (*i) =

(XVIII.24)

Долю нефти в добываемой жидкости определяют как отношение расхода

нефти к суммарному расходу воды и нефти. Используя закон Дарси в трещинах в виде

Цщ —---- Fв (si)

др_

 

дх 9

 

 

Н-в

 

 

 

kx

г

др

(XVI 11.25)

и т ~ ~

^

Рн (5х)ж

 

и уравнение

неразрывности

 

ивгТ Ищ

 

Я it)

 

(XVI 11.26)

 

hS

9

 

 

 

получаем

 

 

 

 

Ф = _________ q(t) |iB______

(XVI11.27)

дх

hSkx [Fв (s{) + ixqFh (si)]

 

Здесь Wbi» tiH — скорости фильтрации воды и нефти в системе трещин; kx — про­ ницаемость трещин. Тогда из (XVI 11.25) и (XVIII.27) доля нефти в добываемой жидкости составит

f(Н)

“hi

\

= (1 - F ( Sl))x=Li.

 

 

ит + Пв1 / x=L

 

В результате согласно

уравнению (XVIII. 16) получим

 

 

1

 

при t < tx,

(XVII 1.28)

 

 

 

 

Коэффициент использования запасов нефти определяется как отношение объема нефти, извлеченной к рассматриваемому моменту времени из пористых блоков пласта за счет капиллярной пропитки, ко всему извлекаемому объему нефти. Это дает

 

J*<T>dx

 

(XVI11.29)

 

hSm2s2Ll

t < i i.

 

 

 

 

 

11

 

 

 

j Лб

 

x

 

Li

403

Фильтрационное сопротивление получим из (XVI 11.27):

 

 

 

J£в_ Г

dx

 

(XVI11.30)

 

 

 

hSki J

Fв (Si) + jliq/^h (si)

 

 

 

 

о

 

 

 

 

Для простоты положим, что кривые относительных проницаемостей в тре­

щинах

прямолинейны,

т е.

 

 

FB(s1) =

s1,

 

 

 

(XVIII.31)

Fh (si) =

1— Si.

 

 

 

Тогда

 

формула

(XVI11.30)

упростится:

 

Q —

Цн

(

L-i — (1 — Mo) J

(XVI11.32)

 

-

Skxh \

 

 

 

 

где

 

 

 

 

t ^ tj ,

 

 

 

 

( x(th

 

(XVIII.33)

H

 

t

, .

 

t > t V

 

 

 

 

 

Из уравнения

(XVIII. 16)

имеем

 

F(Sl) = 1 -

-у- j

(p(t — T(x))dx.

(XVI11.34)

 

 

 

 

 

и

 

 

 

Подставив

(XVIII.34)

в формулу (XVIII.32) и

применив

ь

х

 

 

ь

 

 

 

j

j

/ {У)dy dx = J

(b — у) f (у) dy,

 

a a

 

 

 

a

 

 

 

получим

 

 

 

 

 

 

Ax

 

 

 

 

 

 

 

 

 

 

Skxh

L1 — (• — Mo) Ml (I — /(H)) + -y - J X(f (t — x (*)) dx) ,

 

 

 

 

 

 

 

 

0

 

 

 

 

 

 

 

 

(XVIII.35)

где Аг

задано

формулой (XVIII.33), а /(Н) — формулой (XVIII.28).

Конкретные расчетные формулы будут получены из соотношений (XVI 11.24), (XVIII.25), (XVIII.29), (XVIII.30), (XVIII.35) после подстановки в них зави­

симостей ср (/ — т) и г|о =

(/ — т).

Будем считать скорость закачки q (/) постоянной и равной q. Перейдем к без­

размерным

величинам:

 

B=

 

q

(XVIII.36)

H

a^hSm.ySoL^ 9

 

0 =

a2/,

0! = a2/lt

(XVIII.37)

s =

(XVIII.38)

ё

PZ-x

 

Здесь Xz — характерное время капиллярной пропитки блока; р — коэффициент,

характеризующий отношение времени капиллярной пропитки к времени разработки для рассматриваемой аппроксимации функции капиллярной пропитки.

404

Тогда уравнение (XVIII.23) примет вйД

I (в) е_ (о—е (D)

= Кя ИЛИ

J ^ е - 6 (?))

о

о

6 d% d%- 1Лгсе°.

о

Зависимость £ (0), полученная в работе [29] как решение этого интеграль­ ного уравнения при условии £ (0) = 0, в наших обозначениях может быть за­ писана так:

1(0) = Ф(Кб) (е + 4 " )

+ ] /- ^ - е ~ ° .

 

(XVI11.39)

Тогда

безразмерная

продолжительность

безводного

периода 9Х согласно

(XVI 11.24)

определяется

как

корень уравнения

 

ф (/§ 3 (0! + -§-)

+ У

• ^ е_0‘ = - у

(xvi 11.40)

Все остальные показатели определяют после подстановки в формулы (XVIII.28), (XVIII.29), (XVIII.34) зависимостей г|о (t — т) и ф (t — т) из фор­ мул (XVIII.21) и (XVIII.22) и I (0) из формулы (XVIII.39). Для удобства введем следующие обозначения:

Si (б) = ф ( /б )

(б + 4 “)

+ Y i t

е~6

 

 

 

 

(XVIII.41)

 

 

 

 

 

 

 

 

 

 

г б—v

 

dz\

h

(V. 6) =

1— ф (Vу) Ф (Кб — у) —

 

е

 

v J

*2 +

 

 

 

 

 

 

 

 

 

 

 

У

 

 

 

 

 

 

 

 

 

 

 

 

(XVIII.42)

h

(У. б) =

/ v r

 

 

 

 

 

 

 

2 j

 

I 6Ф (г) (Кя*Ф (г) ег2- 6 + 2e~6) Кб - г2 -

 

----- (26

Г 1—22) arctg •

_____

е

z!_ e] d2 +

 

 

K S '" U

, ‘

‘ 7 “'~‘ь

Кб + К

б ^

 

Г

 

 

+

е

(46 +

1)Ф (/у) +

2е у Уу

arctg ■

Vv

 

 

 

 

 

 

 

 

 

 

 

 

Кб + К б- v

 

 

ev К я

Кб -

уФ2(Ку)(2б+ 1)j ;

 

 

 

 

(XVIII.43)

 

 

 

 

Уу

 

 

 

 

 

 

(XVIII.44)

J3 (У. б) = 2

j

Ф (К в - 2 * )(г Ф (г )4 --^ -^ * .

 

Функции Ух и У3 удовлетворяют условиям J1 (6x6) = 1, У3 (6хб) = б.

405

Тогда все основные формулы можно представить в виде

ll (01) =

1

 

 

 

 

(XVIII.45)

Р ’

 

 

 

 

 

 

 

 

 

i (0) =

ix (0),

 

0< 0!,

 

(XVI11.46)

- j- , 0>Oi.

 

 

 

 

 

 

/(Н, =

{■

1,

 

0 < 0 t,

(XVI11.47)

(0i,

0)*

0

0i»

 

 

 

 

,

г

ре,

 

о <

0j

 

(XVI 11.48)

11

1

РА (0х.

0).

e > 0 i.

 

 

 

•М9. 0),

 

 

 

 

е<01,

^

 

J2 (01> 0) Н

2р”

— ^l(01»0))»

(XVI 11.49)

 

0 > 01-

Фильтрационное сопротивление Q через безразмерную функциюф выражается

следующим образом:

 

 

 

 

 

 

 

 

 

 

 

 

 

(XVIII.50)

Итак,

 

последовательность

расчетов

следующая:

1) вычисляют

коэффициент р

по формуле

(XVI 11.36);

2)определяют из уравнения (XVIII.45) значение 0Х;

3)вычисляют по формулам (XVIII.46)—(XVIII.49) с учетом зависимостей (XVIII.41)—(XVIII.44) безразмерные функции £, /(н), Tj', ф для ряда значе­

ний 0;

4) переходят к размерным величинам по формулам (XVI 11.37)—(XVI 11.39) и (XVIII.50).

Во ВНИИ составлена программа расчета приведенных показателей процесса разработки на языке АЛГОЛ (ГДР) для БЭСМ-6. Но можно при расчете поль­

зоваться и

готовыми графиками или таблицами функций

(6), Ух (у, 6), J2 (у,

6),

Jз (Уъ

б)» рассчитанными по формулам (XVIII.41)—(XVIII.44). На

рис.

XVII 1.5 представлена зависимость £х от (6), по которой можно определять

безразмерную координату фронта воды в момент 0 и безразмерную продолжитель­

ность безводного периода 0Х. На кривых Jj (у, 6) ( / == 1,2, 3), представленных

на рис.

XVIII.6—XVIII.9, значения

пара­

метра у

приведены вдоль кривых J : (6, 6).

Примерный

расчет зависимостей

fa (/),

1]' (/), Q (/) проведен при следующих исход­

ных данных: q =

1,9 м3/сут; Л =

10 м;

5 =

= 100 м; Lx = 100 м; s£ = 0,5;

т 2 =

0,07;

 

|Иц = 0,81 мПа-с; jlih = 2,7 мПа-с и для трех

 

вариантов коэффициента

а:

0,06267 сут 2;

 

_ 1

_ 1

 

0,02834 сут

2; 0,01618 сут

2, что соот­

 

ветствует значениям 0: 0,138; 0,675 2,074.

 

Из рис.

(XVIII.9—XVIII. 11) видно, что

 

при Р = 0,138 все рассматриваемые зависи­

Рис. XVIII.5. Зависимость функции

мости почти

совпадают

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

зависимостями для поршневого вытеснения в

I от параметра б

среде с пористостью т 2,

проницаемостью /гх

406

Рис. XVIII.8. Зависимость функции /3 от параметров б и v

Рис. XVIII.9. Зависимость доли нефти

пользования извлекаемых запасов от вре-

в добываемой жидкости от времени:

/ —поршневое вытеснение;

2 —Р =

мени.

2 —р =

0.138; 3 - Р = 0,675; 4 -

р = 2,074

1 —поршневое вытеснение;

 

 

~ 9,138; J —р —0,675; 4 -

Р = 2,0/4

408

Рис. XVI11.11. Зависимость филь­

трационного сопротивления от вре­

мени.

 

1

— поршневое вытеснение;

2

Р = 138; 3 —3 = 0,675; 4 -

0

=

2,074

 

Рис. XVIII.12. Зависимость степени

Рис. XVIII. 13. Зависимость степени

ис­

использования извлекаемых запасов

пользования извлекаемых запасов от врс-

 

 

 

 

_ _1_

 

от времени

при а = 0,06267 сут2

мени при а —0,01618 сут

2 и Р = 2,074:

и Р= 0,138:

м3/сут; 2 — q —

1 —q —1,9

м3/сут; 2 —q —3,8 м'/сут;

I —q = 1,9

= 3,8 м3/сут; 3 —q —5,7

м3/сут

3 —q = 5,7 м’/сут.

 

 

и остаточной

водонасыщенностыо sj

(показаны

пунктирными линиями). На

рис. XVIII.12

и XVIII.13

приведены графики

зависимости

степени

исполь­

зования запасов от времени для трех значений

скорости закачки q.

 

Из рис. XVIII. 13 видно,

что если Р достаточно велико

(в данном случае

Р > 2,074), то рассматриваемый показатель процесса разработки слабо зависит от скорости закачки. Кривые, соответствующие q = 3,8 м3/сут и q = 5,7 м3/сут, почти сливаются. Таким образом, если отношение характерного времени капил­ лярной пропитки к характерному времени разработки велико, то увеличение скорости закачки не приводит к заметному увеличению добычи нефти из тре­ щиновато-пористого пласта.

409

Глава x ix

 

МЕТОД МАТЕРИАЛЬНОГО

БАЛАНСА,

ЕГО ИСПОЛЬЗОВАНИЕ И

ОГРАНИЧЕНИЯ

С помощью метода материального баланса решает ряд задач разработки нефтяных, газовых и газоконденсатных месторожденийЕго применяют для опре­ деления начальных запасов углеводородов в пласте Иустановления объема вторг­ шейся в продуктивный пласт воды из водоносного бассейна.

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

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

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

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

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

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

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

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

410

щиноватых и трещиновато-пористых коллекторах осуществляется исключительно на основании решения уравнений материального баланса (см. гл. XVIII).

О б о б щ е н н ы е

у р а в н е н и я

м а т е р и а л ь н о г о

б а ­

л а н с а . Жидкости и

газы движутся в продуктивном пласте к забоям пробу­

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

Если в пласте имеется газовая шапка, то нефть первоначально была насы­ щена газом, в залежи отсутствует упругое расширение. Энергия, заключенная в растворенном газе, пополняется за счет энергии газа в шапке.

Но методу материального баланса влияние газовой шапки учитывается ее относительным объемом:

^_ начальный пластовый объем свободного газа

— начальный пластовый объем нефти

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

Для

расчетов необходимы следующие данные.

1. Начальное пластовое давление и среднее пластовое давление в последо­

вательные

моменты времени

разработки эксплуатационного объекта.

2. Добыча товарной нефти

(в стандартных условиях), накопленная к соответ­

ствующим моментам времени

разработки.

3.Общая добыча газа с вычетом его объема, закачанного обратно в пласт.

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

кернов.

5. Коэффициенты пластового объема нефти и газа, а также коэффициент растворимости газа в нефти. Их получают в функции давления на основании

лабораторных исследований образцов глубинных проб.

 

 

6. Отобранное

количество

воды

из пласта.

 

 

 

 

7. Количество воды, вторгшейся в продуктивный пласт из бассейна питания

(в том числе закачанной).

 

 

начальный

объем

нефти

И з м е н е н и е

о бъ е м а неф ти . Если

в пласте равен Q36„о (6Но — объемный коэффициент нефти при начальном

пла­

стовом давлении), то к моменту времени / разработки,

когда среднее давление

в пласте составит р, объем оставшейся в пласте нефти будет (Q3 — QH3) 6Н, где

QH3 — накопленный к моменту времени t объем отобранной из

пласта нефти;

Ьн — объемный коэффициент

нефти,

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

давлению р.

 

Уменьшение объема нефти в пласте к моменту времени t определится как

разность между начальным и остаточным объемами:

 

 

 

 

Д = Qg^Ho — (Сз — Онз) &н-

 

 

 

 

(XIX.1)

И з м е н е н и е

о б ъ е м а с в о б о д н о г о

г а з а .

Отношение на­

чального объема свободного газа в пласте (не только скопившегося в шапке, но и рассеянного в виде мелких пузырьков в объеме пласта) к начальному объему

нефти можно определить так:

GcBrо

(XIX.2)

” Qa*H0

 

Тогда начальный объем свободного газа

(XIX.3)

С?сВг0 == ГРз^но.

 

411

Здесь Gc — объем свободного газа в пласте в стандартных условиях; Вг0 — объемный коэффициент газа при начальном пластовом давлении.

Объем оставшегося в пласте газа к моменту времени t определится как разность между начальным объемом газа, свободного и растворенного в нефти, и объемом газа, добытого из пласта и оставшегося в нем в растворенном виде, т. е.

ГОбъем свободного

газа! _

ГОбъем

газа,

свободного и растворенного]

[

ко времени t

J

[

при начальном давлении р0

J

_ГОбъем добытого газа]

ГОбъем

газа, оставшегося]

 

[

ко времени t

\

[в растворе ко времени t J

 

Запишем это равенство в другом виде:

 

 

Gc = £

Ь Qe^roJ

[ФизЯр] — [(Оз

Физ) Яг].

(XIX.4)

Уменьшение объема свободного газа в пласте в результате его добычи и снижения давления теперь легко определить как

[Уменьшение объема свободного] —Дг —Q3fbmгаза в пласте

— [

qH0

^3^ г0 — QmRp — (Q3 — Q113) ЯР] а Р>

(XIX.5)

где Вг — объемный

коэффициент газа (коэффициент сжимаемости)

при теку­

щем пластовом давлении; Rr0 и Rr — объем растворенного газа в единице объ­ ема нефти соответственно при начальном и текущем пластовых давлениях.

И з м е н е н и е о б ъ е м а воды . В процессе отбора нефти и ее заме­ щения водой объем воды в пределах продуктивной части пласта будет всегда увеличиваться. Это увеличение можно выразить следующим образом:

[•Увеличение объема] = Дв = (W + Wp — bBWU3) — W = Wp —

[

воды

\ —

(XIX.6)

где W — начальный

объем воды в продуктивном пласте;

Wp — объем вторг­

шейся к моменту времени t воды из питающей области; WU3 — объем добытой к моменту времени t воды; Ьв — объемный коэффициент воды.

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

Фз^но 4- QsBbuo — Q3 [bn "Ь (Яго — Rr) Вг) 4* Qh3 [^н ~г (Яго — Яг) ^г]

+ (Rp - Яго) BrQua - Q™r— bmB? = Wp - bBW*3.

(XIX.7)

 

 

£>ro

 

Обозначим

 

Я0 =

[6н +

(Яг0- Я г)Яг].

(XIX.8)

Тогда

6Н0

= Д0о-

 

Здесь В0 — двухфазный коэффициент пластового объема нефти при давлении р (физический смысл этого коэффициента понятен из приведенного соотношения); В0о — то же, при начальном пластовом давлении.

Используя (XIX.8), преобразуем уравнение (XIX.7) и решим его относи­

тельно Q3:

 

<?из [Во + (Rp - Яго) Вг]~ (Ц7Р - bBWna)

(XIX.9)

Во— В00 + —7ТГ*- (Вг — вгд)

 

&Г0

 

412

Соседние файлы в папке книги