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

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

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

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

Добавлен: 15.10.2024

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

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

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

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

Определив по уравнению (1) зависимость р = р (t), из уравнения технологического режима эксплуатации, записанного в виде рс (t) =

= р (t) — б, находим зависимость изменения во времени забойного давления в «средней» скважине рс = рс (<).

При известных зависимостях р = р (t) и рс = рс (t) уравнение (2) притока газа к скважине представляет квадратное уравнение отно­ сительно дебита «средней» скважины. Решая его, имеем

4(f) =

A (p*z)Cp ,

/

(

A (p*z)cp

у

I P2 (t) — P l(t)

(8)

2>Bzcp

у

\

2Bzcp

)

Bzcp

По найденным величинам пластового и забойного давления вы­

числяются

значения ([i*z)cp

и zcp. Вычисленные значения

(p*z)cp

и zcp, а

также

пластового

и

забойного

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

в формулу (8)

и определяются

 

величины дебитов «средней» сква­

жины q в различные моменты времени.

 

 

 

Зная изменение во времени среднего дебита газовых скважин q = = q (t), по формуле (7) определяют изменение во времени потреб­ ного числа газовых скважин п = п (t).

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

Дифференциальное уравнение истощения газовой залежи запи­ сывается в виде:

nq (t) =

a Q „

 

p (t )

(9)

Рат

dt

z q,)

 

 

Решение дифференциального уравнения (9) возможно только по шагам. Разделим переменные в уравнении (9) и проинтегрируем их

по времени от t до t -f At и

по -Дг- от

-£=-)

до (-^=-1

:

 

z(p)

\z (р) It

\z(p)/t+At

 

<+fд<n (t)q (t)d t= - A Ш---------

 

(10)

aQH J

z [p («)J

z [ p ( t + A«)]

 

Индекс или аргумент t или t + Ai в выражении (10) означает принадлежность соответствующей величины времени t или t -j- At. Считается, что в момент времени t решение задачи известно, и реше­ ние отыскивается на момент времени t + At.

i Здесь и в некоторых других параграфах методика расчетов приводится применительно к технологическому режиму эксплуатации скважин при допусти­ мой депрессии на пласт. Для иных технологических режимов порядок расчетов аналогичен рассмотренному в § 1 данной главы.

102


При достаточно малом шаге At зависимость nq = п (t) q (t) близка к прямолинейной. Тогда, заменяя площадь криволинейной трапеции площадью прямолинейной трапеции, имеем

Рат

п (t) д (t) + n p -f-Лг) д (t+A t) At :

Р 00

р (t + At)

(И)

21р Р+ д01

afiH

2

zip (HI

 

Согласно уравнению (2) и уравнению технологического режима эксплуатации скважин при допустимой депрессии на пласт х, запи­ санным для момента времени t -j- At, находим

p(t + At) =

b(t + At)

4(P*z)Cp£+A£

•q(t +

,

B z c p t + A t

q2(t + At).

 

 

26 {t

 

At)

 

A*) > 26 (t + At)

 

 

 

 

 

 

 

 

 

 

 

 

(12)

Подставляя выражение для p (t -j- At) в (И ),

получаем

- - / -g -

[n(t)q(t) + n(t + At)q (*+Д*)] =

_ 4 Э

---------- 6 £ + **> _ ■

2aQ„

1

w v v / T

V T

/

ч \

 

л

z[p(0]

2z[p(£ + A£)]

4

(P

z)Cp t + A t

■q(t + At)— ^

 

Bz,cp t+ A t

q2 (t + At).

26 (£+

 

 

 

At) z [p (< + дг)]

Д£) z [p (t -f Д£)]

 

 

 

26 (t +

 

Решим полученное квадратное уравнение относительно q (t + At):

q(t + At) = —

n(t + At)paT .

,

-4 (P*z)cp иМ

X

aQH

■At +

 

 

 

 

 

 

 

 

6 ( г + Д г )

z [p (i-j-Ai)]

 

 

 

X 6 (£-|- At) z [p (i + Дг)] .

 

 

 

 

 

 

 

2Bz,cp t+ A t

 

 

 

 

 

+

 

n (t

At) рат

AH

 

-4 (P*s)cp t + A t

 

 

 

 

 

 

 

 

6 (t +

At) z [p (г+дг)]

 

6 (t -f- Д£) z [p (£+ Д£)]

 

+

2

p (0

n (г) рат At

9 ( 0 -

X

 

2Bz,'cp t + A t

 

z [P (<)]

 

aQH

 

 

 

 

 

 

 

 

 

 

6 (t+ A i)

\

6 (£ +

Д£)г[р(£ +

Д£)1

(13)

 

 

 

z [p {t + A 0I j

 

 

Bz4pt+At

 

 

 

 

 

 

 

 

 

 

Расчеты по формуле (13) проводятся в следующей последователь­ ности. В первом приближении все величины, относящиеся к вре­ менному слою t + A t , вычисляются по результатам решения задачи на момент времени t. По этим приближенным величинам с исполь­ зованием формулы (13) вычисляется средний дебит газовых скважин

1 В общем случае приходится для периода падающей добычи газа рассмат­ ривать различные варианты по изменению во времени депрессии на пласт. Поэ­ тому здесь и в дальнейшем предполагается, что зависимость 6 = 6 (г) известна (задана).

103


к моменту t + At. Определяется отобранное количество газа за рас-

сматриваемыи интервал времени At ------> ' — 1— - A

Затем по формуле (1) вычисляется значение пластового р (t + At) давления в момент времени (t + At) (значение пластового давления можно определять также с использованием формулы 12), а по фор­ муле

рс(t -f- At) = р (t -f At) —6 (t 4- At)

определяется величина забойного давления. По найденным значе­

ниям р (t + At) и рс (t -f At) уточняются величины в формуле (13), относящиеся ко времени t -f- At. Вновь вычисляются дебит сква­

жины

q (t +

At), отобранное количество газа за время At, р (t +

+ At),

рс (t +

At) и т. д. до тех пор, пока значения дебитов сква­

жин q (t +

At)

в последнем и предпоследнем приближениях не бу­

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

ющей добычи газа: р = р (t); рс = р с (t); q = q (t); Q = Q (t). После окончания расчетов с шагом по времени At расчеты повто­ ряются при шаге At/2. Если при этом разница, например, в дебитах скважин окажется не больше е, то шаг At не завышен. В противном случае расчеты проводятся с шагом Дг/4 и т. д. в описанном выше

порядке.

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

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

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

ных свойств газа записываются

в виде

(при

режиме

допустимой

депрессии на пласт):

 

 

 

 

р ( * ) = ( - т - ~

Рат% - " —

Н

”);

(14)

\ zh

ocQH

/

 

 

p ( t ) - P c ( t ) = S ;

 

 

(15)

p 2 (0 - Pl (t) =

A q (t) + B q 2 (t)]

(16)

n (t)

^ ^

 

 

(17)

ПW

q (t)

 

 

 

Порядок решения системы уравнений (14)—(17) очевиден и не нуждается в особых пояснениях.


Для периода падающей добыли газа справедлива формула

Рат

п (t) q (t)+ n (t+ A t) q (f + A Q

_ p (t)

P(*+AQ

айн

2

z[jp(f)]

z [p (i + Ai)l

Из уравнения притока газа к скважине (16) с улетом (15) имеем

jf J_ д^ч

Aq (t +

At)

,

Bq2(t + At) .

6 (t -t-AQ

(19)

P V , t s i ) -

2 6 (t

Д г )

'

2 6 ( г + Л г ) *

2

 

Подстановка (19) в (18) дает

Рат

П (t) q (t) + n (t +

At) q (t + At)

^ _

P(t)

aQu

2

 

 

z fp(«)]

 

Aq (t + At) + Bq2(t + At)

6(<+ At)

 

26 (t- f А г ) z [ p

( г + Д г ) ]

2z [/> ( г - } - А г ) ]

Решая данное уравнение относительно дебита газовой скважины в момент времени t + At периода падающей добыли газа, полулаем

q (t + At) =

П ( t + А г ) р а т

+ /

/,и (0 Р ат9(*)д

п (t + At) Рат At-

X

 

ocQH

6 (t +

At) z [p ((4-Ai)]

X

6 (t +

At) z [p (i-j- A#)]

+

 

 

2 В

 

At

 

А

6 ( г + Д г ) z [ р ( г + А г ) ]

 

 

2В

• 6(t + A t)z [p (t + At))

^ _ 9

p (t)

6 (t+ A t)

\6 (t + A i)z[p (t + Ai)]

\

айн

z[p(<)]

z[p(i+Ai)] )

В

(20)

Раслеты по формуле (20) и последующие расчеты по определению других показателей разработки при падающей добыче газа прово­ дятся аналогично ранее рассмотренному случаю.

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

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

105


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

момент

времени

t + At.

По

найденным значениям р (t

At) и

Р а (t +

At) (i =

1, 2, . .

n)

определяются (p*z)cp ,.<+д< и zcp iUAt .

Затем с использованием уравнения притока газа к каждой скважине определяются дебиты отдельных скважин на момент времени t + At. Суммирование дебитов по отдельным скважинам дает величину от­ бора газа из месторождения в целом на рассматриваемый момент времени. По формуле

<?,об (* + д о = е доб ( о +

ы

(21)

уточняется значение добытого количества

газа на

момент вре­

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

Нетрудно видеть, что предлагаемый процесс последовательных приближений является сходящимся. Если при расчетах в первом приближении задано завышенное значение добытого количества газа <?доб (t -j- At) на момент времени t + At, то найденное с исполь­ зованием уравнения материального баланса среднее пластовое давле­ ние окажется заниженным. Соответственно будут заниженными зна­ чения дебитов скважин и величина Q (t + At) в (21). Поэтому в рас­ четах второго приближения будет задаваться величина @доб (t -|- At), вычисленная согласно (21), меньшая, чем использованное значе­ ние Qдоб (t + At) при проведении расчетов в первом приближении, и т. д.

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

 

П

 

n

 

 

 

 

2

(P*z)cpf l i

n2

BiZCp.q?

 

А =

i=l

_______

 

п

2

(22)

 

n

 

 

(P*z)cp 2 91

zcp

2

9i

 

 

 

i=1

 

. 1=1

 

 

При этом уравнение притока газа к «средней» скважине записы­

ваем в виде:

 

Р2 (0 - Р! (*)'= А (Р*2)сР? (*) + Bzcpq2 (t).

(23)

106