Файл: Закиров, С. Н. Проектирование и разработка газовых месторождений учебное пособие.pdf

ВУЗ: Не указан

Категория: Не указан

Дисциплина: Не указана

Добавлен: 15.10.2024

Просмотров: 151

Скачиваний: 0

ВНИМАНИЕ! Если данный файл нарушает Ваши авторские права, то обязательно сообщите нам.

Таким образом, при решении численными методами необходимо трационные сопротивления от R c

Дюпюи:

_ 2nkh '

задач подземной гидродинамики учитывать дополнительные филь­ ф до R c, например, по формуле

ро pc

/gx

р, R c . ф

R c

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

Если учесть, что при фильтрации газа одно из условий подобия записывается в виде:

CRCq . 2ратСр

то получим, что при моделировании газовых скважин величина добавочного сопротивления определяется той же формулой (5). Радиус фиктивной скважины при фильтрации газа (по закону Дарси) также равен 0,2077 шага сеточной области.

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

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

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

Уравнение притока газа к реальной скважине может быть пред­ ставлено в виде двучленной формулы

р 2 — p l = a q + b q z.

Здесь коэффициенты фильтрационных сопротивлений а и Ъ

относятся к области, заключенной между радиусом, равным Дж, и радиусом действительной скважины R c.

Уравнение (2) перепишем в виде:

Р2~ Рс

Чa-\-bq

Отсюда следует, что

i —•Cqq Cq

P2— -Pc

(9)

 

a-\-bq

 

141


Тогда, сопоставляя уравнения (1) и (9) и учитывая второе усло­

вие подобия в случае фильтрации газа

c qCR

1,

получаем

 

 

 

 

2ратСр

 

 

r

Р2— Ро.^ f

Ро — Рс _

2ратСр

p2 — pl

(10)

р

_1_ д

р -йдоб

 

Сд

a -)- bq

 

По правилу производных пропорций из соотношения (10) имеем

_______ Ро — Рс_______

_ Ро— Рс

_£й _ (а + б9) _ 1 Л

ЛД°б

2рат

4

 

Следовательно, добавочное сопротивление при нарушении закона Дарси в призабойной зоне будет

R доб :

(aJr b q ) — - j - R .

( И )

 

2рат

 

Нетрудно видеть, что из этого соотношения в частном случае, если положить b = 0, получаем формулу (5).

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

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

Величины фильтрационных сопротивлений а л Ъ, входящие в фор­

мулу (11), трудно вычислить с достаточной достоверностью. Учиты­ вая, что основные потери давления при притоке газа к скважине приходятся на призабойную зону, рекомендуется вместо коэффи­ циентов а и Ъ подставлять в уравнение (И ) величины фильтрацион­ ных сопротивлений А и В , определяемые по данным интерпретации

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

по

формуле

 

 

 

■^доб

 

и

при нелинейном законе

сопротивления — по формуле

 

* д о б —

^

^ ~~

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

Используем уравнение притока реального газа к скважине

Р г — Pi = Л (p*z)cp q + B z cvq2

(12)

142


или при учете также реальных свойств пласта приближенно запишем

р2— p l =

A ^ ^ - ^ ^ g + B zcpq2.

(13)

Зависимость ф = ф (р)

перестраиваем в зависимость ср =

ф (р2) .

Эта зависимость хорошо аппроксимируется двумя (тремя или более)

прямолинейными отрезками. Пусть зависимость ф = ф (р2)

аппрокси­

мируется следующими

уравнениями:

 

 

 

<P =

<?iP2+

£ i ,

p i «£ Р*

; р 1;

(14)

Ф =

C2p 2 +

D 2,

р\ s= р 2

;р *;

(15)

 

Ф = С3р 2,

0 s£ p 2

■ p i

 

(16)

Тогда для первого интервала изменения квадрата давления имеем

 

р 2 =

 

Ci

R l

 

 

 

 

 

 

Ci

 

 

 

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

и рс приходятся на интервал

[р2, р|].

В

этом случае,

например,

уравнение (12) принимает вид:

 

 

 

 

 

Ф — Tc = ^C'1 (ii*z)cpg +

5 C 1zcpgr2.

(17)

Отсюда

 

 

 

 

 

 

 

___ ____________________ ф — ф с ___________________

^АС у (p*z)Cp-|-.6CizCp<7

Тогда вместо соотношения (10) имеем

г

Ф-Фо _

P*tC9

ф— фс

Г

Фо— фс

, ЛО\

Ф

I n

CR

ACi(ii*z)cp^-BCyzcpq - Vrf

Лдоб •

^

 

4

 

 

 

 

 

Следовательно,

 

 

 

 

 

 

^доб =

[АСу <ji*z)cр + BCyZcpq\ -

1

R .

(19)

Так же определяется 7?доб для двух других выделенных интерва­ лов изменения квадратов давления: Су заменяется соответственно на С 2 или С 3.

Чтобы не получалось принадлежности р 2 и р 2 разным интервалам аппроксимации, можно поступать следующим образом. Как только р\

достигнет величины р 2, производится новая аппроксимация зависи­ мости ф = ф (р2). Новая аппроксимирующая линия проводится, например, так, как показано на рис. 44 (пунктирная линия а).

Данная аппроксимирующая зависимость используется в расчетах до тех пор, пока р 2 не достигает левого конца интервала. После этого расчеты продолжаются с использованием аппроксимации (15) и т. д.

При нарушении закона Дарси непрерывный процесс решения задачи возможен лишь в случае пренебрежения реальными

143


свойствами газа и если дебиты скважин не изменяются во времени. В этом случае величина добавочного сопротивления зависит от де­ бита скважины q [см. формулу (8)]. В условиях реального газа

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

 

 

 

 

 

Р и р с.

 

 

 

 

 

 

 

 

При достаточно малом шаге по

 

 

 

 

 

времени и монотонном изменении де­

 

 

 

 

 

битов

величину Д доб для интервала

 

 

 

 

 

времени

U,

t +

Д£] можно вычис­

 

 

 

 

 

лять

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

 

 

 

 

 

решения задачи на момент времени t.

 

 

 

 

 

Более точные результаты могут быть

Рис. 44.

К

пояснению

идеи

получены

при проведении расчетов

во втором приближении. Тогда при

«скользящей»

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

 

 

 

 

 

вычислении 7?доб для интервала вре­

 

 

 

 

 

мени

[t,

t +

АН

используются ре­

шение задачи на момент

времени t и результаты решения задачи в

первом приближении на момент времени

t +

A t .

При этом, напри­

мер, z

вычисляется по

значению

давления,

определенному по

формуле

 

 

 

 

 

 

 

 

 

Р = J

{ у [ р (0 +

Рс (О ] + j [ р (t +

At) - f pc (t + Д*)]j •

 

§

9. Расчеты

по

разработке

газовой

залежи

 

 

на электрических моделях с учетом

 

 

общей депрессионной воронки

 

При разработке месторождений природного газа образуется той или иной «глубины» общая депрессионная воронка. Это означает, что в каждый момент времени давление в пласте от точки к точке изменяется. Образование депрессионной воронки в процессе разра­ ботки Северо-Ставропольского месторождения показано на рис. 45.

При рассмотрении (в § 7) возможности применения электрических моделей для интегрирования дифференциальных уравнений неустановившейся фильтрации газов вводилась новая временная пере­ менная т. Ее введение (линеаризация уравнения) основывалось на допущении, что давление (в формулах 19 и 38 § 7) в каждый момент времени в любой точке пласта близко к среднему пластовому давле­ нию по залежи на соответствующий момент времени. При наличии депрессионной воронки значительной «глубины» это допущение может приводить к ощутимым погрешностям.

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

144


Н . Г. Степанова и одного из авторов книги. Эти методы основаны на линеаризации исходной задачи на отдельных временных шагах.

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

Внешний нонтур

Внешний нонтур

газоносности

газоносности

Рис. 45. Формирование общей депрессионной «воронки»:

Г — 20/XII

1957 г.; I I — 15/V 1958 г.; I l l — 10/Х 1958 г.;

IV — 25/ХП

1958

г.;

V —

20/IV

1959

г.;

VI — 10/VIII 1959 г; VII — 5/ХП 1959 г.; VIII — 20/Ш

1960

г.;

IX —

20/VI

1960

г.;

X — 15/IX

1960 г.; X I — 15/Ш

1961 г.; X I I — 20/IX 1961 г.: 1 — линия

распределения

Тоэф^эф!

тэф— эффективная

пористость;

эффективная

мощность

 

 

 

на

Северо-Ставропольском месторождении

 

 

 

дебитами. В случае идеального газа задача сводится к решению урав­ нения

 

£[■ £-£]+£[

 

 

 

<*>

при следующих условиях:

 

 

 

 

 

 

 

 

 

 

р =

р н = const; t = 0;

 

 

(2)

 

 

 

=

 

 

 

* =

1> 2> • •

п>

 

(3)

 

 

 

 

4

f

= 0, (X, у)ег.

 

 

(4)

Представим

уравнение

 

(1)

следующим

образом:

 

д

Г kh

dp2 “I

,

 

д

Г й

др2 ~|__

amh

др2

,г,

А'г

I

м

й'г !

^

АII I

и

Fill

n

A f *

' '

Ю Заказ 1013

 

 

 

 

 

 

 

 

 

145