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

книги / Основы разработки нефтяных и газовых месторождений

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

значением последнего дебита ^п. Разница между ними заключается в том, что если и I и I » Д1, то точки на графике А, соответствующие определенным значениям р^$, смещены по горизонтали на величину

г

+д*

1 + М

18 п

м

18 М

по отношению к таким же точкам на графике В. Поэтому графики на рис. 7.20 смещены один относительно другого по вертикали на вели­ чину ш 1§ (1п /1) для заданного А!. Это смещение можно интерпрети­ ровать как разность

Р,‘п - р* = Ш 1§ 1а.,

(7.57)

где р*1пи р* - полученные экстраполяцией значения р^з (ын) при Д{ = °о, соответственно, для фактической и эффективной продолжи­ тельности работы скважины. Кроме того, принято, что при обыч­ ном исследовании последний период работы скважины достаточно продолжителен, чтобы произошел переход к квазиустановившейся фильтрации. Тогда уравнение Мэтьюза, Бронса и Хейзбрека (7.44), по которому можно рассчитать р, имеет вид

 

Ро(МВН) (*сл) = 471

кН

(7-58)

 

^ ~ (Р' - Р) = 2>303 !§ (САи

 

а цВ

(7.59)

ИЛИ

р*-р = 0,186 ^

- 21§(СА1ОА) = т1§(СА11)А.

Уравнение (7.59) пригодно при использовании эффективной про­ должительности работы скважины, а для фактической продолжи­ тельности применяется зависимость

Р,'п-Р (п = т18 (СаЧ ,)-

(7-60)

Вычитая уравнение (7.59) из (7.60), получаем

(р ',п - Р*) ~ (Р,п - Р) =

К

(7.61)

- р

Сравнивая полученную зависимость с уравнением (7.57), можно видеть, что р( = р. Поэтому определение среднего пластового дав­ ления по методу Мэтьюза, Бронса и Хейзбрека дает одинаковые ре­ зультаты независимо от того, используется при анализе фактическая или эффективная продолжительность работы скважины.

Рассуждая подобным образом, несложно показать, что среднее пла­ стовое давление, определенное по результатам исследования, не зави­ сит от того, какая продолжительность работы скважины используется при анализе. Это утверждение справедливо при условии, что продол­ жительность работы скважины равна или превышает время 1$$$, требу­ емое для перехода к квазиустановившейся фильтрации в дренируемом объеме, и что в анализе используется последний дебит. Иллюстрацией этого утверждения является график С на рис. 7.20, построенный для продолжительности работы скважины, равной 18$$. В этом случае гра­ фик С смещается по горизонтали относительно графика В, построен­ ного для эффективной продолжительности работы скважины, так что зависимости, эквивалентные выражениям (7.57) и (7.61), выглядят как

 

Р'-Р^ = т18 ^

 

и

(Р " Р») - - Р8Й) = т

 

 

858

Все это показывает, что метод анализа Мэтьюза, Бронса и Хейз­ брека дает одно и то же значение р, используется ли для построения КВД время I или 1888. К такому же выводу пришли в своих работах Пинсон (Ршзоп)12 и Каземи (Кахегш)13.

Следует также отметить, что значение скин-фактора, определен­ ное при анализе, также не зависит от того, какая продолжительность работы скважины применяется. Это объясняется тем, что значение Р^8(ыы) 1ь> используемое для расчета 5 (уравнение 7.52), не зависит от того, какая продолжительность работы скважины используется при анализе, и одинаково для графиков А, В и С на рис. 7.20.

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

на момент остановки скважины существовал квазиустановившийся режим фильтрации. И даже если это условие не выполняется точно, ошибка будет небольшой. Использование такого сочетания дебита и продолжительности работы скважины может быть неприемлемым при первом исследовании, когда скважина работала относительно недолго и с непостоянным дебитом. Оде (ОйеЬ) и Зелиг (8еБ§)14 пред­ ставили метод анализа данных восстановления давления в таких условиях, который может дать более точные результаты. При даль­ нейшем описании анализа восстановления давления будет исполь­ зоваться исключительно эффективная продолжительность работы скважины, обозначаемая I, и последний дебит, обозначаемый ^ При­ мер использования 1888 в анализе восстановления давления при ис­ следовании газовых скважин будет дан в главе 8 (раздел 11).

Построив график Хорнера по измеренным значениям давления, можно применить метод МБХ для определения р, действуя следую­ щим образом:

1) Экстраполировать начальный прямолинейный участок до зна-

I + Д{

чения 1§ = 0 и определить значение р*. Используя уравнение (7.50), определить к по угловому коэффициенту прямой.

2) Разбить пласт на отдельные дренируемые объемы, так чтобы

где - дебит 1-той скважины, дренирующей объем V., а ^тот и Утот - соответственно суммарный дебит и полный объем пласта. Эта зависимость, представленная в главе 5 (раздел 5.3), приемлема для скважин, дренирующих пласт при квазиустановившемся режиме фильтрации. Однако Мэтьюз, Броне и Хейзбрек утверждают, что эту зависимость можно использовать с достаточной точностью незави­ симо от того, какой режим фильтрации существует. Она позволяет определить V.. Отсюда можно получить и площадь области, дрениру­ емой скважиной (А.) - при допущении, что средняя толщина пласта в области дренирования равна толщине, измеренной в скважине. С помощью структурной карты пласта можно приблизительно опре­ делить форму области дренирования и положение скважины отно­ сительно ее границы, и по этим данным выбрать соответствующий график МБХ (рис. 7.11-15).

3) Определить безразмерное время (уравнение 7.49)

*0А= 3600 —- (часы),

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

4) На соответствующем графике МБХ, рис. 7.11 - 15, найти кривую для геометрической формы, в наибольшей степени схожей с опреде­ ленной по геологическим материалам, и снять ординату р0(мвн) (*0А), соответствующую расчетному значению безразмерного эффектив­ ного времени (продолжительности) работы скважины 1ОА. Расчетное значение р0(мвн) (*0А) равно (уравнение 7.58).

Поскольку на первом этапе было определено р*, можно напрямую рассчитать р. Следует отметить, что координаты 1ОАи р0(мвн) графи­ ков МБХ могут быть выражены в единицах разных систем, так как оба эти параметра безразмерны.

Известен и другой метод определения р, представленный Дитцем9.

Он предусматривает расчет значения

I + Д{

, при котором следует

определять значение р по экстраполированной прямой на графике Хорнера, как показано на рис. 7.21.

Пусть А*. - продолжительность остановки скважины, при которой гипотетическое давление на экстраполированной прямой равняется среднему пластовому давлению. Тогда р^5(ш0 в уравнении (7.48) рав­ но р, и его можно переписать следующим образом:

(р* - р) = 1,1511§

+ р0 (1в) - 0,51п

5

*

Левую часть этой зависимости можно определить с помощью уравнения материального баланса (7.38). Тогда

Р*- Р)

Рис. 7.21. Метод Дитца, позволяющий определить среднее пластовое

давление р и динамическое давление в ячейке сетки р6

Если выразить р0 в этой зависимости в общем виде, с использова­ нием уравнения (7.42), то

2,3031§

1 + А*

 

= 4 п кЬ

 

Д(.

$

— Рэ(МВН)

(Р* - Р)

 

 

 

 

ЧРВ„

 

 

ч

1 + А{

(р -р)

 

или

_____$

(7.62)

 

м

Ш

Это уравнение прямой, угловой коэффициент которой ш = 0,186 ^I^В0 / кЬ, демонстрирует эквивалентность методов Дитца и МБХ. Об этом же свидетельствует рис. 7.21.

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

Ро(мвн) ~ 2,303 1§ (СА1оа),

и поэтому

1 + А*

' =Ч(Са 1оа).

(7.63)

Ч

Физическая непроницаемая граница

Граница ячейки при численном моделировании

' Рис. 7*22. Модель пласта с непроницаемой границей области, дре­

нируемой скважиной А, и наложенная квадратная сетка для численного

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

 

По этой зависимости можно определить значения 1§

для ра­

боты с графиком Хорнера.

8

Усовершенствованный метод Дитца определения р часто исполь­ зуется для сравнения измеренных значений давления со средним давлением в ячейке сетки, полученным на численной модели.

Допустим, что модель для численного моделирования содержит несколько ячеек, находящихся внутри естественной непроницаемой границы области дренирования (рис. 7.22). При моделировании в конце каждого временного шага рассчитывается среднее давление в каждой ячейке. Поэтому применяя интерполяцию во времени дав­ ления, полученного путем моделирования, можно относительно просто определять давление в отдельных ячейках, соответствующее моменту исследования скважины А независимо от того, совпадает этот момент с окончанием временного шага или нет. Есть два пути равнения измеренного давления с давлением в ячейках, полученным путем моделирования.

Первый из этих путей - расчет среднего давления внутри непро­ ницаемой границы на момент исследования методом МБХ или Дит­ ца и сравнение полученных результатов с усредненным по объему давлением во всех ячейках в пределах непроницаемой границы. Это довольно утомительная процедура. Есть более простой приближен­ ный метод, предложенный ван Поленом (уап РооИеп)15 и развитый Эрлафером16. Этот метод предусматривает использование графиков Хорнера в сочетании с методом Дитца для расчета так называемого

«динамического давления в ячейке» ра, которое представляет собой просто среднее давление в ячейке, в которой расположена скважи­ на, на момент исследования. Целью анализа является определение

щее давление, определенное‘по гипотетическому прямолинейному графику Хорнера, равнялось динамическому давлению, то есть р^5 (ыы) = Ра* Здесь также можно применить зависимость (7.63), но в этом случае нужно определять 1ОАисходя из площади ячейки сетки, а не всей площади в пределах непроницаемой границы. Тогда САприни­ мает постоянное значение 19,1. Такой выбор объясняется тем, что граница ячейки сетки не является непроницаемой границей. Поэто­ му граничные условия в большей степени соответствуют условиям установившейся фильтрации, для которых Дитц представил только один случай, когда скважина расположена в центре круговой обла­ сти дренирования. Соответствующий этому случаю коэффициент формы равен 19,1 (рис. 6.4). Таким образом, прямоугольная ячейка сетки аппроксимируется кругом, площадь которого равна площади ячейки. Поэтому на графике Хорнера берется абсцисса

(7.64)

и считывается соответствующее значение ра, как показано на рис. 7.21. Возможность использования уравнения (7.64) также зависит от того, работала ли скважина до начала исследования при квазиустановившемся или установившемся режиме фильтрации. Обычно в том и дру­ гом случае 1: » Д^. Используя это допущение, ван Полей предложил выражение для расчета продолжительности остановки скважины, при которой р^з (уы) = ра. Его можно легко получить из уравнения (7.64) в виде

 

срцсА

_

фцстгг^

Л*4

3600 х 19,1 к

 

68760 к

или

 

 

(ррсг^

МА= 4,57 х 10'5 — :----,

 

41

 

к

где ге - радиус круга, площадь которого равна площади ячейки сетки.

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

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

УПРАЖНЕНИЕ 7.6. АНАЛИЗ РЕЗУЛЬТАТОВ ИССЛЕДОВАНИЯ МЕТОДОМ ВОССТАНОВЛЕНИЯ ДАВЛЕНИЯ. БЕСКОНЕЧНЫЙ ПЛАСТ

Скважина-первооткрывательница работала до остановки для пер­ вого исследования методом восстановления давления приблизитель­ но 100 часов. Ниже приведены оценочные характеристики пласта и пластовых флюидов и дебит.

ЯсТ

=

19,56 ст. м3 / сут

ср

=

0,2

 

 

 

N

=

79,5 ст. м3

И0

=

1 мПа с

 

 

 

р

 

 

 

 

 

 

Ь

=

6,1 м

В .

=

1,22 пл.м3 / ст. м3

 

 

 

 

 

 

01

 

 

 

 

 

г

=

100,6 мм

С

=

2,9 х 10'3 / МПа = (с 5

+ с 5 + сг)

 

 

 

 

 

1

4 О О

№№

V

А* 121,4 га

Значения давления, измеренные в ходе исследования, приведены в табл. 7.5

Требуется определить:

Каково начальное пластовое давление?

Какова фазовая проницаемость (известно, что скважина вскрыла пласт на всю толщину)?

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

Давление в

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

Давление в

остановки скважи­

остановленной

остановки скважи­

остановленной

ны, А1 (часы)

скважине, МПа

ны, Д1 (часы)

скважине, МПа

0,0

31,07 (р^)

 

 

0,5

32,23

3,0

32,84

0,66

32,44

4,0

32,86

1,0

32,63

6,0

32,89

1,5

32,75

8,0

32,91

2,0

32,80

10,0

32,92

2,5

32,83

12,0

32,94

Каково значение скин-фактора?

Каков дополнительный перепад давлений из-за скин-эффекта?

Правомерно ли приравнивать р. к р*, если изначально принято, что скважина расположена в центре круговой области дрениро­ вания?

УПРАЖНЕНИЕ 7.6. РЕШЕНИЕ

1) Эффективная продолжительность работы скважины

I

х 24= (79,5 / 19,56) х 24 = 97,6 часа.

 

*4йпа]

В табл. 7.6 указаны точки, по которым был построен график Хор­ нера в координатах

1 + Д1

_____§

М

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

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

1 + ДГ

,

1 + ДС

р ,МПа

остановки скважины,

м

‘8

Д(

(фунт / дюйм2)

А1 (часы)

 

 

 

 

 

0

 

 

 

31,07(4506) (р_)

0,5

196,2

 

2,29

32,23

(4675)

0,66

148.9

 

2,17

32,44 (4705)

1,0

98,6

 

1,99

32,63

(4733)

1,5

66,1

 

1,82

32,75

(4750)

2,0

49,8

 

1,70

32,80

(4757)

2,5

40,04

 

1,60

32,83

(4761)

3,0

33,53

 

1,52

32,84

(4763)

4,0

25,40

 

1,40

32,86

(4766)

6,0

17,27

 

1,24

32,89

(4770)

8,0

13,20

 

1,12

32,91 (4773)

10,0

10,76

 

1,03

32,92

(4775)

12,0

9,13

 

0,96

32,94 (4777)

Р«в

1+Д{д*

Рис. 7.23. График Хорнера для случая бесконечного пласта

,1 + М _ Л

д|-----и дает р* = 33,10 МПа. Принимая, что при исследова­

нии скважины-первооткрывательницы можно использовать уравне­ ние (7.54), получаем р. = 33,10 МПа.

2) Угловой коэффициент прямолинейного участка КВД ш = 168,9 кПа

/единица логарифма. Поскольку скважина вскрыла пласт на всю тол­ щину, фазовая проницаемость пласта равна

 

0,186

.

0,186 х 19,56 х 1 х 103 х 1,22

 

к "

гаЬ

=

168,9 х 1<Я х 86400 x 6,1 “ 50 *

’ *“ *Л

3) Для определения скин-фактора можно использовать уравнение (7.52). Гипотетическое значение р„5(иы),.Ьг = 32,77 МПа в этом уравне­ нии получено экстраполяцией прямолинейного участка КВД до зна­ чения Д1 = 1 час (рис. 7.23).

5= 1,51

ГФ\У5 (ЫЫ) 1 — Ьг Р ууР _ ^ ^

 

 

(32,77-31,07)

^

^

•0,35

168,9 х 10'3 "

 

ш

 

ФИСГ„

 

 

 

50x10 15

-0,35

 

- 1§ 3600 - 1§ 0,2 х 1 х 10'3 х 2,9 х 103 х 106 х 0,00836

 

 

 

5 ~

6.