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

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

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

ДАВЛЕНИЕ,

Рис. 9.14. Аппроксимация непрерывного снижения пластового давления

рядом дискретных ступеней давления (упражнение 9.2)

В табл. 9.3 приведены значения накопленной добычи нефти Ыр и среднего газового фактора Кр на конец года, а также характеристики РУТ, соответствующие значениям давления, указанным во втором столбце табл. 9.2.1*

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

Порядок решения задач этого типа таков:

1)Рассчитать приток воды с применением теории ^установившего­ ся притока Херста и ван Эвердингена при геС) =10.

2)Интерпретировать уравнение материального баланса как уравне­ ние прямой линии, используя метод Гавлена и Оде (глава 3, раздел 7). Для залежи без начальной газовой шапки уравнение материального баланса запишется следующим образом:

Время, годы

Ыр, млн

Кр, ст. м3 /

Во, пл. м3 /

К, ст. м3 /

В , пл. м3 /

ст. м3

ст. м3

ст. м3

ст. м3

8

3

 

 

 

 

ст. м3

 

0

 

 

П5,7 (К)

1,404 (В„.)

115,7(11..)

0,00522 (Вс.)

 

1

 

1,26

135,28

1,374

105,38

0,00551

 

2

 

2,93

150,41

1,349

97,01

0,00601

 

3

 

4,63

163,76

1,329

90,25

0,00657

 

4

 

6,47

173,55

1,316

83,84

0,00719

 

5

 

7,97

182,45

1,303

78,68

0,00781

 

6

 

9,29

189,57

1,294

74,40

0,00843

 

7

 

10,40

194,91

1,287

70,84

0,00899

 

8

 

11,25

199,36

1,280

68,17

0,00955

 

9

 

11,85

203,81

1,276

66,04

0,00989

 

10

 

12,31

206,48

1,273

64,79

0,01022

 

 

 

 

Таблица9*3

 

 

 

N

р 1 О ' р

 

з' В*

Г(В0-В 0|) + (К,-К)В81

(сА г + д

Др ]+ \У,

 

01

 

1 - 5

[В + (К

- К ) В ] = ЫВ

 

 

 

 

 

где \Уе - суммарный приток воды. Если полное снижение давления за 10 лет равно 8,96 МПа, то член, учитывающий сжимаемость,

 

 

 

 

 

(СА с + С()

Аг

будет равен

 

 

 

1 - 5 \УС

*

 

 

 

 

 

 

(с 5

+ с ()

А

0,435 х 0,05 + 0,58

х 10'3х 10’6х 8,96 х Ю 6^ 0,06.

4 УУ \УС

Г7

Д г -ч —

______

7

7

1 - 5

Р

 

0,95

 

 

 

 

 

Из приведенной ниже табл. 9.6 следует, что максимальное значе­ ние члена, учитывающего расширение нефти с первоначально рас­ творенным в ней газом, составляет 0,277 (Ео / Во.)тах = 0,3895 / 1,404 = 0,277. Отношение этих двух значений составит всего лишь 2 %, и поэтому сжимаемостями породы и воды можно пренебречь. Тогда

уравнение материального баланса сведется к следующему:

 

N [В +(К

- К ) В ] = Ы [ ( В

- В. ) + ( К . - К ) В ] + \ У

е

р 1 О У р

З7

О ОV 4 31 з 7

или

Р = ЫЕ +\У.

 

о е

График зависимости между Р / Е и

/ Ео должен быть линейным,

отсекать на оси ординат (при / Ео = 0) отрезок Р / Ео = N и иметь

угловой коэффициент, равный единице.

3)

Если при геЕ) = 10 такой график будет нелинейным, то нужно варьи­

ровать значение этого параметра, чтобы получить прямую линию.

1)

Расчет при геС = 10

 

Поскольку были выбраны интервалы протяженностью в один год,

удобно определять коэффициент при безразмерном времени, выра­ зив время 1; в годах (уравнение (9.7)).

_ 31536000 к!

~ФН&2

=31536000 х 200 х 1015 / (0,25 х 0,55 х 10'3 х 1,015 х 10'9 х 28042) 1 = 5,751. Постоянный параметр круговой водоносной области (уравнение (9.8))

V = 2п Г<рЬсг2.

Здесь { = 140° / 360° = 0,3889 и И = 6,28 х 0,3889 х 0,25 х 30,48 х 1,015 х 10-3 х 28042 = 148516 м3 / МПа.

Теперь можно рассчитать приток воды, используя значения Др. и из табл. 9.4. Результаты приведены в табл. 9.5.

Время, годы

Безразмерное

Ар, МПа

 

^ 0

время

(г„ = 1 0 )

 

 

(г.о = 5>

 

 

 

 

0

0

0,83

 

 

1

5,67

1,55

4,95

4,88

2

11,34

1,35

8,12

7,46

3

17,01

1Д7

10,90

9,10

4

22.68

1,01

13,50

10,09

5

28,35

0,85

15,90

10,83

6

34,02

0,72

18,10

11,27

7

39,69

0,58

20,20

11,52

8

45,36

0,44

22,20

11,69

9

51,03

0,32

24,00

11,81

10

56,70

 

25,70

11,89

Т,

 

 

п -

1

\У, млн

 

 

= и 1

д Р)\ \ ( т с - 1 0.)

годы

 

 

пл. м3

1

0,14852

(0,83

х 4,95) х 106

 

0,61

2

0,14852

(0,83

х 8,12 + 1,55 х 4,95) х 106

2,14

3

0,14852

(0,83

х 10,90 + 1,55 х 8,12 + 1,35 х 4,95) х 106

4,21

4

0,14852

(0,83

х 13,50 + 1,55 х 10,90 + 1,35 х 8,12 + 1,17 х 4,95) х 106

6,66

5

0,14852

(0,83

х 15,90 + 1,55 х 13,50 + 1,35 х 10,90 + 1,17 х 8,12 + 1,01 х

9,41

4,95) х 106

 

 

 

 

 

 

6

0,14852

(0,83

х 18,10 + 1,55 х 15,90 + 1,35 х 13,50 + 1,17 х 10,90 + 1,01 х

12,34

8,12 + 0,85 х 4,95) х 106

 

 

 

 

7

0,14852

(0,83

х 20,20 + 1,55 х 18,10 + 1,35 х 15,90 + 1,17 х 13,50 + 1,01 х

15,38

10,90 + 0,85 х 8,12 + 0,72 х 4,95) х 106

 

 

8

0,14852 (0,83

х 22,20 + 1,44 х 20,20 + 1,35 х 18,10 + 1,17 х 15,90 + 1,01 х

18,14

13,50 + 0,85 х 10,90 + 0,72 х 8,12 + 0,58 х 4,95) х 106

 

 

9

0,14852 (0,83

х 24,00 + 1,44 х 22,20 + 1,35 х 20,20 + 1,17 х 18,10 + 1,01 х

21,18

15,90 + 0,85 х 13,50 + 0,72 х 10,90 + 0,58 х 8,12 + 0,44 х 4,95) х 106

 

 

 

0,14852 (0,83

х 25,70 + 1,44 х 24,00 + 1,35 х 22,20 + 1,17 х 20,20 + 1,01

24,13

10

х 18,10 + 0,85 х 15,90 + 0,72 х 13,50 + 0,58 х 10,90 + 0,44 х 8,12 + 0,32 х

 

4,95) х 106

 

 

 

Таблица 9.5

Время, годы

Р, млн пл. м3

Ео, пл. м3 /

Р / Е , млн

 

\У /Е о

ст. м3

О

(*->=«>)

(геС = 5)

 

 

ст. м3

1

12,124

0,0268

452,4

142,9

140,9

2

30,761

0,0574

535,9

234,5

223,8

3

52,826

0,0923

572,3

286,7

260,3

4

79,798

0,1411

565,5

297,2

253,5

5

105,964

0,1881

563,3

314,8

251,3

6

132,292

0,2380

555,8

362,2

243,8

7

157,080

0,2862

548,8

338,2

236,8

8

179,177

0,3299

543,1

352,5

231,2

9

196,654

0,3630

541,7

373,6

229,7

10

210,743

0,3895

541,1

396,4

229,1

2) Расчет материального баланса Р = N (В + (К - К ) В ),

Ео =((Во - Во.) + (К. - К) Вр.Как видно на рис. 9.15, после первого года зависимость между Р / Ео и / Ео отклоняется от теоретической прямой с единичным угловым коэффициентом. Это говорит о том, что правильное значение геГ) должно быть несколько меньше. Поэтому

расчет повторяется с ге0 = 5.

 

3) Расчет

при ге0 = 5

 

Значения

при геГ) = 5 приведены в табл. 9.7. Для их расчета ис­

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

из столбца 5 табл. 9.4.

Соответствующие значения / Ео, требуемые для расчета матери­ ального баланса, приведены в табл. 9.6. По этим значениям построен график, показанный на рис. 9.15. Как можно видеть, здесь все точки

т,

\У = 1 7 П1 Д р .^ 0 (Т0 - 1

)

\Уе> млн м3

годы

(млн барр.)

}-о

>

1

0,14852 х (0,83 х 4,88) х 106

 

0,602 (3,79)

2

0,14852 х (0,83 х 7,46 + 1,55 х4,88) х 106

 

2,043

 

(12,85)

 

 

 

3

0,14852 х (0,83 х 9,10 + 1,55 х7,46 + 1,35 х 4,88) х 106

3,818

(24,01)

 

 

 

4

0,14852 х (0,83 х 10,09 + 1,55 х 9,10 + 1,35 х 7,46 + 1,17 х 4,88) х 106

5,682

(35,74)

 

 

 

5

0,14852 х (0,83 х 10,83 + 1,55 х 10,09 + 1,35 х 9,10 + 1,17 х 7,46 +

7,511

1,01 х 4,88) х 106

 

(47,24)

 

 

6

0,14852 х (0,83 х 11,27 + 1,55 х 10,83 + 1,35 х 10,09 + 1,17 х 9,10 +

9,222 (58)

1,01 х 7,46 + 0,85 х 4,88) х 106

 

 

 

 

7

0,14852 х(0,83 х 11,52 + 1,55 х 11,27+ 1,35 х 10,83+ 1,17 х 10,09 +

10,768

1,01 х 9,10 + 0,85 х 7,46 + 0,72 х 4,88) х 106

 

(67,72)

 

 

8

0,14852 х (0,83 х 11,69 + 1,55 х 11,52 + 1,35 х 11,27 + 1,17 х 10,83 +

12,115

1,01 х 10,09 + 0,85 х 9,10 + 0,72 х 7,46 + 0,58 х 4,88) х 106

(76,19)

 

9

0,14852 х (0,83 х 11,81 + 1,55 х 11,69 + 1,35 х 11,52 + 1,17 х 11,27 + 1,01

13,248

х 10,83 + 0,85 х 10,09 + 0,72 х 9,10 + 0,58 х 7,46 + 0,44 х 4,88) х 106

(83,32)

 

 

0,14852 х (0,83 х 11,89 + 1,55 х 11,81 + 1,35 х 11,69 + 1,17 х 11,52 +

14,170

10

1,01 х 11,27 + 0,85 х 10,83 + 0,72 х 10,09 + 0,58 х 9,10 + 0,44 х 7,46 +

(89,12)

 

0,32 х 4,88) х 106

 

 

 

 

лежат на прямой линии с единичным угловым коэффициентом. Эта прямая отсекает на оси ординат отрезок Ир = 49,6 млн ст. м3 (312 млн ст. барр.). Таким образом, подтверждается предположение, что ге0 = 5 - правильное значение отношения радиусов водоносной области и залежи.

Из рисунка видно, что после третьего года происходит обращение тренда, при этом и / Ео, и Р / Ео начинают со временем убывать. Гавлена и Оде указали4,5, что такого обращения тренда и следует ожидать в случае активной ограниченной водоносной области, но при пра­ вильном значении ге0 все точки должны располагаться на прямой ли­ нии. При некорректном выборе модели водоносной области (ге0 =10) расчетные значения / Ео постоянно возрастают при уменьшении Р / Е , как видно на рис. 9.15.

100

200

300

400

 

\Л/е/Е омлн

ст.барр.

 

Рис* 9*15* Подгонка модели законтурной водоносной области по методу Гавлена и Оде

9 .4 . П РИ БЛ И Ж ЕН Н А Я ТЕОРИЯ Ф ЕТКО ВИ ЧА П РИ ТО К А

ВО Д Ы В ЗА Л Е Ж Ь Д Л Я СЛУЧАЯ ОГРА Н И ЧЕН Н О Й

ВО Д О Н О СН О Й ОБЛАСТИ

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

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

(9.18)

где - дебит притока, ] - коэффициент продуктивности водоносной области,

р - пластовое давление, то есть давление на ВНК или ГВК, ра - среднее давление в водоносной области.

Для определения среднего давления в водоносной области исполь­

зуется простое уравнение материального баланса

 

=

р ,- р а).

(9.19)

где р. - начальное давление в водоносной области и в залежи. Это со­ отношение можно выразить иначе:

(9.20)

Здесь \У. = с\У.р. определяется как максимально возможный объ­ ём притока воды в залежь из водоносной области и характеризует максимально возможное расширение водоносной области. Диффе­ ренцирование уравнения (9.20) по времени дает

аш

XV ар

а*

(9.21)

р. а*

Подставляя уравнение (9.21) в (9.18) и разделяя переменные, получаем

______’А ц

р. - р

Теперь проинтегрируем это уравнение, приняв в качестве началь­ ного условия, что в момент 1 = 0 (\^е - 0, ра = р .) на границе залежи создается депрессия Др = р. - р. В течение всего рассматриваемого периода давление на границе р остается постоянным, так что

1“ (Р.-Р,) = - М + с.

е!

Здесь С - произвольная постоянная интегрирования, которую можно определить с учетом начальных условий как С = 1п (р. - р). Отсюда

Ра - Р = (р> - р) е~ 1р‘*'

(9.22)

Подстановка этого выражения в уравнение притока (9.18) дает

 

<№

. -Ы /У Г .

 

___е

= Нр. -р)е

(9.23)

а

 

 

 

Наконец, интегрируя уравнение (9.23) при принятых начальных условиях получаем следующее выражение:

Ш

(р, - р) (1 - е"

1

).

(9.24)

 

Р.

 

 

 

Как можно видеть, если Xстремится к бесконечности, то

\У.

\\г = — - (р4- р) = Щ (р4- р) р*

Это максимально возможный объём притока воды в залежь при распространении зоны депрессии р. - р по водоносной области.

Вообще-то уравнение (9.24) нельзя считать очень полезным, по­ скольку оно получено при условии постоянства давления на вну­ тренней границе. Для использования этого уравнения в практиче­ ских расчетах, когда давление на границе постоянно изменяется во времени, следует также применить принцип суперпозиции. Однако Феткович показал, что можно использовать уравнение (9.24) в дру­ гой записи, благодаря чему отпадает необходимость в суммирова­ нии решений. Для притока на первом интервале уравнение (9.24) можно записать так:

ДШ

(Р4-Р1)(1 -е-1Р* Д^/^еО.

(9.25)

е 1

 

 

Здесь рх- среднее пластовое давление на границе залежи на первом интервале. Для второго интервала Д*

^

- 1р д*,/\у.ч

(9.26)

Д^е2 = ---” (р&1“ Р2) (1_е

^ 2

е‘)>

Р1

 

 

 

где ра1 - среднее давление в водоносной области в конце первого ин­ тервала, которое определяется из уравнения (9.20)

Ра, = Р,

V/е).

')

(9.27)

 

 

 

В общем, для п-го интервала

 

 

 

 

Ш = — 1а (ра

- Р„) (1 - е

Й>. А1п/ Шй),

(9.28)

р,

1,-1

 

 

 

 

п -

1

 

 

 

I

Ш

 

 

где

 

 

 

(9.29)

Значения среднего давления на границе залежи рп рассчитываются, как описано в разделе 9.3 (уравнение (9.15)):

Феткович показал, что, последовательно применяя уравнения (9.28) и (9.29), можно рассчитать приток воды при различных геоме­ триях водоносной области. Полученные результаты хорошо согласу­ ются с результатами расчета с использованием теории ^установив­ шегося притока Херста и ван Эвердингена для случая ограниченной водоносной области.

Формулы для расчета коэффициента продуктивности водонос­ ной области ], зависящего и от геометрии и от условий фильтрации, приведены в табл. 9.8. Выражения для ] для круговой геометрии при квазиустановившейся и установившейся фильтрации имеют такую же форму, как выражения для коэффициентов продуктивности при притоке жидкости в скважину, описываемом зависимостями из табл. 6.1 в главе 6. Единственное различие заключается в том, что вместо радиуса скважины г^ берется радиус залежи г . Следует отметить также, что, в то время как выражения для ] для условий квазиустано­ вившейся фильтрации (9.30) используются совместно с уравнениями Фетковича (9.28) и (9.29), выражения для условий установившейся фильтрации (9.31) используются иным образом. При работе с этими выражениями принимают, что вода, притекающая в залежь из водо-

Режим фильтрации

Водоносная область

Водоносная область

 

круговой геометрии

линейной геометрии

 

 

1 м3 / с / Па

м3 / с / Па

 

Квазиустановившаяся филь­

1кЪ

кЬ\\г

 

трация (с выражением депрес­

/

г

з \

(9.30)

3

сии в виде ра - р)

-ГГт)

цЬ

 

 

 

 

 

 

 

 

 

 

Установившаяся фильтрация (с

1кЪ

кЬ\\г

 

выражением депрессии в виде

 

,

г

(9.31)

ц

цЬ

р ,-р )

1п

го