1209
.pdfИспользуем связь между горным давлением по вертикали рг, средним нормальным напряжением сг и внутрипоровым (пластовым) давлением р, определяемую формулой (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 йтк„зэки— |
=©„мш. |
|
|||
ш |
|
|
|
|
|