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

1209

.pdf
Скачиваний:
2
Добавлен:
15.11.2022
Размер:
11.79 Mб
Скачать

Используем связь между горным давлением по вертикали рг, средним нормальным напряжением сг и внутрипоровым (пластовым) давлением р, определяемую формулой (11.64).

Из формулы (11.64) следует, что при pr=const

да

др

~дГ~ ~

dt

Учитывая (III.2) и (Ш.З), получим

дт дт да

Q да

Q др

~дГ ~ ~дсГ ~дГ ~ —Рс~дГ ~

~дГ'

(Ш.З)

(Ш.4)

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

Р = Ро[1+М />—A))].

(HI.5)

где рж — сжимаемость жидкости;

р0 — плотность жидкости при

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

 

Из (III.5) имеем

 

д£___др_ _др__п О _др_

(III.6)

dt ~~ др dt — PoP« dt

 

Используя закон Дарси и считая проницаемость k и вяз­ кость жидкости р, не зависящими от координаты, имеем

£

 

 

(HI.7)

divpu= -------div р grad р.

 

 

Подставим (III.4), (Ш.6) и (III.7) в (III.1). В результате

получим следующее выражение:

 

 

рРс -^г+"ФоРж

div Р

Р

(HI-8)

Учитывая незначительную сжимаемость жидкости, в фор­ муле (III.8) можно положить р«ро. Тогда окончательно по­ лучим дифференциальное уравнение упругого режима в сле­ дующем виде:

-§Г = х div grad р;

* = -ДГ:

Р= Рс+™Рж-

(III-9)

 

Здесь к и р — соответственно пьезопроводность и упругоемкость пласта (по предложению В. Н. Щелкачева).

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

101

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

- ^ а- = р д р;

р = Ы - т р ж,

(Ш Л О )

где АУП— изменение порового объема, т. е. непосредственно упругий запас пласта объемом V; АУП и Ар — абсолютные ве­ личины.

П р и м е р III.1. Рассмотрим пример определения упругого запаса нефтяно­ го месторождения. Пусть имеем месторождение, продуктивный пласт которого имеет объем V=109 м3=1 км3. Это — довольно большое месторождение, на­ пример длиной 20 км и шириной 5 км. Толщина пласта составляет 10 м. Пред­ положим, что начальное пластовое давление равно 16 МПа, а давление насы­ щения — 6 МПа. Продуктивный пласт непосредственно за контуром нефте­ носности выклинивается, замещаясь практически непроницаемыми породами. Поэтому в нефтеносную часть месторождения из его законтурной области по­ ступает незначительное количество воды. Можно рассчитывать только на упругий запас нефтеносной части месторождения. Вопрос ставится следующим образом: сколько можно добыть жидкости из пласта при снижении пластового давления Ар на 10 МПа от начального пластового до давления насыщения?

Для ответа на этот вопрос необходимо определить упругий запас нефтя­

ного месторождения

с

учетом указанных его

особенностей. Пусть (3 =

= 10-4

1/МПа.

(ШЛО)

 

Тогда, согласно

 

ДУП=

1ф Др= 10в-10-М 0 = 10« м3.

 

Таким

образом, можно

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

10 МПа упругий запас

месторождения составляет

1 млн. м3.

§ 2. ПРОГНОЗИРОВАНИЕ ИЗМЕНЕНИЯ ДАВЛЕНИЯ НА КОНТУРЕ НЕФТЯНОГО МЕСТОРОЖДЕНИЯ ПРИ УПРУГОМ РЕЖИМЕ В ЗАКОНТУРНОЙ ОБЛАСТИ ПЛАСТА

Для разработки месторождения важно знать изменение давления во времени на условном контуре нефтеносности мес­ торождения Ркон= Ркон(0 или средневзвешенного по площади

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

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

102

/ г J 4

1 Z

Рис. 59. Схема разбиения площади

Рис. 60. Схема нефтяного месторож­

нефтяного месторождения и его за­

дения круговой формы в плане:

 

контурной водоносной

области

на

1 — условный

контур

нефтеносности;

2

ячейки:

 

 

 

аппроксимация

контура

нефтеносности

ок­

I — контур выклинивания водоносной обла­

ружностью радиусом R

 

 

сти месторождения;

2 — ячейка площадью

 

 

 

 

ДхДу\ 3 — условный

контур

нефтеносности;

 

 

 

 

4 — аппроксимация контура

нефтеносности

 

 

 

 

месторождениях с высоким

содержанием

парафина в нефти

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

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

Расчет изменения во времени средневзвешенного пластового или контурного давления при геометрически сложной конфигу­ рации контура нефтеносности с учетом реального расположе­ ния скважин на месторождении возможен только с использо­ ванием численных методов и ЭВМ или аналоговых устройств.

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

103

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

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

тяным залежам в законтурных областях во

многих

случаях

действительно близок к радиальному, происходящему

как бы

в залежи круговой формы в плане.

разрабатывается

Итак, пусть месторождение (см. рис. 60)

на естественном режиме и, вследствие сравнительно незначи­ тельного упругого запаса энергии в нефтяной залежи, будем считать количество отбираемой жидкости из месторождения qm(t) равным количеству поступающей воды к нефтяной за­ лежи из законтурной областщпласта q ^(t), т. е. <7ж(0 ~<7зв(0-

При разработке нефтяных месторождений добыча жидкости

qm(t) изменяется обычно так, как это показано

на рис.

57.

Для расчета

pl<0K(t) будем считать законтурную

область

не­

ограниченной

(-ft^rs^oo). Радиальная фильтрация

воды в этой

области описывается дифференциальным уравнением упругого режима (III.9), которое в рассматриваемом случае принимает следующий вид:

(III.11)

где p (r,t) — давление в точке А с координатой г в законтурной области пласта (см. рис. 60).

Рассмотрим вначале несколько упрощенную задачу упруго­ го режима, для которой начальное и граничное условия запи­ сываются следующим образом: р —роо при t=0, R ^ r ^ .o o \

(III. 12)

Решение этой задачи получают с использованием преобра­ зования давления p(r, t) по Лапласу

оо

р (г, s) = j р (г, t) e~stdt,

о

где p(r,s) — преобразованное давление; s — параметр разования.

(III. 13)

преоб­

104

Рис.

61.

Зависимость

/(1, т)

от

lg(l+T):

 

 

 

 

/ — точное

значение

функции /(1, т)

по

Ван

Эвердингену и

Херсту;

2 — аппрокси­

мация функции формулой (III.15)

Рис. 62. Зависимость q3B{h) от К

В общем виде это решение по имеет следующий вид:

P» - P (P. T) = S - / ( P. T);

 

 

 

 

 

 

(III.14)

 

(1 _ е - “ат) [Jj (и) Y0 (u p ) - Yx (и) J0 (up)] du

 

 

 

«‘ IVW + V M l

 

 

 

p= rJR,

т = KtlR2.

 

 

 

 

 

Здесь J0(«p), J i(u), Yo(up), Yi(a) — функции

Бесселя.

и Хер­

Функция /(р,т) была рассчитана

Ван Эвердингеном

стом.

расчета изменения во

времени давления

 

необ­

Для

p Ko n ( t )

ходимо

использовать значения

этой

функции

при p= r/R= \

(рис. 61).

f(l,t)

от lg (l+ t)

можно с не­

Оказалось, что зависимость

обходимой точностью аппроксимировать следующей достаточ­ но простой формулой:

/ (1, т) =

0,5 [1 —е -8-77 о+*>]-}-1,12 lg (1+ т )

 

или

 

 

/ (1, т) =

0,5 [1 —(1 —х)-3»81] —(—0,487 In (1 —j—т).

(III. 15)

Таким образом, для //ж = const давление pKon(t) можно рас­ считать по формуле, вытекающей из выражений (III.14) и (III.15):

< Ш 1 б >

105

Однако добыча жидкости в процессе разработки месторож­

дения, естественно, не остается неизменной во времени.

 

Рассчитать изменение рКон(7) при переменном

во времени

<73B~<73B(0 можно с помощью интеграла Дюамеля.

 

q3в=

Для получения этого интеграла будем рассматривать

= q3п(т) и считать, что q33 изменяется со временем

не

непре­

рывно, а ступенчато, причем каждая ступенька А<7ЗВ1начинает­ ся в момент времени V Используем два времени: т. исчисляе­

мое с начала

разработки

месторождения.}

и

Л, с~Ьтдельными

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

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

ступеньками

Aq33l=

= const.

образом,

дебит жидкости q3B будет

зависеть

теперь

Таким

уже не от т, а от Л,,- или просто от % (рис. 62).

 

сообра­

В соответствии

с формулой (III. 16), изложенными

жениями

и рис. 62 можно

написать следующее выражение:

 

 

 

ff3B

 

 

 

 

Ркон 00 = Р°°

2nkh

^

(^»1')~Ь

 

 

 

 

 

 

о

 

 

 

 

~bA<73Di/(l, т —^О+А^зв2/ 01т —^2) + • • •] =

 

 

 

 

^зв

 

 

 

 

 

- Р - — 5 Ш

- 2 Л?га‘/(1’ Т -Х |)-

 

 

(Ш-17)

 

 

О

 

 

 

 

 

Разделим и умножим выражение, стоящее в правой части

(III. 17) под знаком суммы,

на АЛ,. В результате получим

A - W -

Л

.

-

/ О. Т - УДЛ .

 

 

(Ш.18)

 

 

 

О

 

 

 

 

Перейдем в (III.18) к пределу, полагая АЛ,— »-0. Тогда для любого Л, (индекс i можно опустить) имеем

Интеград-ДНЛЭ) и есть интеграл Дюамеля.

При разработке нефтяпых—месторождений отбор жидкости из пласта изменяется во времени обычно таким образом, что вначале он нарастает в связи с разбуриванием месторождения

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

скважин,

а затем стабилизируется на

значительное

время

и лишь

в

конце срока разработки снижается.

 

 

из

Однако если учитывать, что

приток воды происходит

законтурной области пласта,

то

снижение

поступающего

ее

106

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

Учитывая сказанное, схематизируем изменение текущего отбора воды из законтурной области пласта во времени в об­ щем случае следующим образом:

1) <7ЗВ= аХ при 0 < Я < ^

2) <7зв = ?зв1 = const при

=

 

 

(III.20)

 

 

 

 

3) Язв= Язв1~а^ ПРИ К < х < К

 

 

 

 

4) '7зв= <7з„2 = const при

 

 

 

 

 

При этом время Я*=т* соответствует началу закалки в за­

контурную область ТюдьГ^В^йошшт времени

' к \ = Х \

месторож­

дение оказывается полностью

разбуренным и

отбор

воды

из

законтурной области стабилизируется,, В момент

т=т* начи­

нают вводить в эксплуатацию

нагнетательные

скважины

в

законтурной области и приток из нее воды, затрачиваемой

на

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

закачка

воды в законтурную область пласта может быть такова,

что она

не только компенсирует добычу жидкости из нефтяного

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

Рассмотрим вначале изменения контурного давления рКон =

=Ркоъ(Я, г)

в первом из указанных случаев, т. е. при О ^^^Л .1.

Из (III.20)

имеем

= а = const. Тогда

т

0

+ 0,487 In [1 + (т—Я.)]j dX =

+ 0,487

I„ [1 + (т —X)] <ftj

-----J (x);

J (t)-0,5t-0,178[ 1 - (I+ ')3,„ ]+

+0,487 [ ( 1 + T) In(1+t)—T ].

(Ill.21)

Чтобы

получить

изменение ркон=Ркон(т) при Я^Я.Ь необ­

ходимо из формулы

(111*21)

вычесть рК н=Ркон(т) при Т>Ть

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

В результате получим при T> TI

fa) “

 

W “ т а г J (т~ ч ) ] -

в Р»—

U fa)—^fa—Tl)b

(1П.22)

В третьем случае, т. е. при т>т*. из выражения для ршш= =*Ртъ(ъ) по формуле (111*22) необходимо вычесть решение,

соответствующее изменению в третьем случае (Ш-20). Имеем

 

^

(Ш-23)

В четвертом случае при &>т** получим

Ршаfa) *

^ fa) ^ fa

 

— З г

 

(Ш Д4)

Рассматриваемая задача может ставиться и итапыгм образом.

Задают

давлепие ртш(Ц и определяют

f 3№= f 3Jffl(TJ.

Прнмепепие современных математических методов ш шычюсднтельиых средств позволяет учесть изменение параметров в законтурной о&ласщ, ее ограниченность и другие осшшдаяговдие факторы*

Вместе с тем пе всегда можпо использовать слшшые м го - матичгешпе методы и вычислтельшуво технику. В сишуадряж» трефуютх шшучешйа быстрого ответа* прюшишлг простые, ко вдесвшййо менее точные расчетные схемы. Так* длит шршблнжешото вдршаознрованш изменения дашлешшя jp^=iPead[tt}

ш т ш считать, что месторождеше

вводился в разработку в

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

с некоторым

ШОСЛОШШШИШ ДебоПШШТ

1П1усть шзшсть нефти близка к вязкости воды, шроашваешншь и толщина власти в ело нефлешасыпоешшшй ч аст ш за шредаивд условного, среднего контура нефтешшост ((рис. ©8)) <вдшшви>-

Ш

Рис. 63. Схема нефтяного месторождения с тремя точечными стоками:

/ — условный контур нефтеносности

вы. За контурное давление р Ко и (0 будем условно принимать давление в точке А, расположенной на расстоянии В от оси х.

Для приближенного расчета изменения во времени давления рКн(Х) применим следующий прием: будем считать, что отбор жидкости из всех скважин нефтяного месторождения qm заме­

няется отбором из трех, пяти или другого

числа п точечных

стоков с дебитом qi, так что

 

?

1

(1П-25)

 

 

 

Пусть, например, согласно рис. 63

 

<7ж= ^0+^1+^2*

(III.26)

Точечный сток q0 расположен в начале координат, а стоки cf\

и q2— слева и справа от него на расстояниях соответственно — а н а. Тогда, используя соответствующую формулу гл. II, полу­

чаем выражение

для приближенного

определения

изменения

давления во времени в любой

точке

пласта на

расстоянии

г=Ух2~гу2 от начала координат:

 

 

 

(0 = Р® Р (О =

Ei (

1йг)

 

 

Н

 

 

[ _ f e ± g ± * |.

(Ш.27)

Из (III.27)

имеем формулу

для

определения

изменения

давления в точке А (см. рис. 63).

 

 

АРкое(0 — Р»

Ршш(О —

(

ba )

 

4HML

у

4эй Jr

 

 

(Ш.28)

Рассмотрим примеры расчета контурного давления.

 

 

 

П р и м е р

III.2.

Глубокозалегающее

 

 

небольшое по размерам нефтяное месторож­

 

 

дение, контур нефтеносности которого име­

 

 

ет

форму,

близкую

к форме

круга

(см.

 

 

рис. 60),

окружено

обширной

водоносной

 

 

областью,

которую можно

считать прости­

 

 

рающейся до бесконечности. Начальное пла­

 

 

стовое давление в нефтяной залежи, как и

 

 

на его

контуре, при r = R = 3-103 м состав­

 

 

ляло р«о=20 МПа. Проницаемость пласта

 

 

в

законтурной

водоносной

области

k =

 

 

= 0 ,1

мкма, вязкость

воды

ц = 1 0 -3 Па-с^

 

 

упругоемкость пласта 0 = 10-9 1/Па, толщи­

 

 

на водоносного пласта h— 10 м.

 

 

 

 

 

Количество воды, поступающей из за­

1

г

контурной части месторождения в его неф­

тенасыщенную часть,

определим

по форму­

Рис. 64, Зависимость давления

ле

(III.20). При этом Ti=2 года, т .= 4

го­

да, a i= a = 0 ,1 3 6 8 м3/сут.

 

 

 

иа контуре месторождения рко«

 

 

 

от времени t

 

 

Найдем изменение контурного давлении

 

в течение первых пяти лет разработки ме­

 

 

 

 

сторождения.

водоносного пласта. Имеем

Определим прежде всего пьезопроводность х

кП0-1»

10~»* 10~* = 10-1 м^с*

По формуле (Ш.14)

ий

10-1i

О, ш М 0-1 *= о„96- 10-s I

 

^ - =

~ =

 

где i — ® еутг.

 

 

 

 

Выгдашам ие

а &Ркш(ъ)=р<ж—Рйши(т). При £=2 года=730 сут имеем

х—офв*ш-»*»=о,7оа.

 

 

По формуле

 

 

 

 

 

 

 

=2,182-1(1»;

/ М - 0 <5 - 0 Л « » - 0 , т [ | ] +

+© „487RS +

 

lm 1„7Ш8— O.TOOfq= 0,311;

Д рю я=2ъ118Й'№'©Л11 = 6J 8

МПа.

 

Прт

 

года

следует ш г а ш т »

то формуле (Ш~22)_ Имеем

‘С »Ф ,Д О № -М 095= 1,051;

х1 =0,7ШШ;

 

^ = © „ 5 * 1 „ ш — |®i ,—m

 

) +

+ ®Ъ4 ® 7 ^ Ш

Д1ш2чШ — 1,0)51)

 

 

 

 

 

 

ш ,т 0 —-

) +

+ <М®7(1jam йтк„зэки—

=©„мш.

 

ш

 

 

 

 

 

Соседние файлы в предмете [НЕСОРТИРОВАННОЕ]