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

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

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

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

Добавлен: 15.10.2024

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

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

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

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

нение

 

 

 

 

_д Г к ( х , у , р ) h ( х , у )

 

d p i ~1_

д Г к ( х , у , p ) h ( x , у )

 

'

д х \ _

р. ( р) z ( p )

' д х J ' д у L

р(Р ) z (Р)

д у

J

 

 

= 2a{x,y)m{x,y)h{x,y)^y-^ -S^

 

 

(1)

при следующих

условиях:

 

 

 

 

 

 

 

t =

0; р =

р ( х ,

z/) =

pH= const; (х, у) 6G;

 

 

(2)

 

 

 

| г =

°«

( х , у ) е Г .

 

 

(3)

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

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

Фг (Н-**)с tqt -г т|^с {qi = C tp c t.

(4)

Непосредственное удовлетворение последнему условию при реше­

нии задачи (1)— (4)

не представляется возможным. При решении

задач, аналогичных

(1)— (4), на ЭВМ или электрических моделях

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

скважин.

 

 

 

 

 

 

 

 

Поэтому

при

решении задачи

(1)— (4)

воспользуемся связью

между дебитами и забойными давлениями реальной р с

и фиктивной

р с ф скважин,

устанавливаемой

уравнением

притока

газа

 

 

 

P i фг — Рс г = аг 0 * * 4 * i 4 i + V

c Pi3?-

 

 

Здесь at

и

Ъ{ — коэффициенты

фильтрационных

сопротивлений

зоны пласта между R c и R c ф.

 

 

 

 

 

Основные

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

при фильтрации газа

приходятся

на призабойную

зону пласта.

Поэтому уравнение

притока газа

к скважине с достаточной для практических расчетов точностью может быть записано в виде 1:

P i < s ,' i - P h = a i (F*z)cp. tqt + Я (* СР

(5)

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

167


Здесь A t и В t — коэффициенты фильтрационных сопротивлений,

определяемые по результатам исследования скважин при установив­ шихся режимах; i — порядковый номер скважины.

Пусть на момент времени t решение задачи (1)— (4)

известно.

Требуется найти решение задачи на момент времени t +

2 А t.

В первом приближении при отыскании решения задачи на момент времени t + 2 А t дебиты скважин принимаются равными дебитам

скважин определенным на момент времени t q-t (t). Значения (p*z)Cpii+At и zcpit+At также принимаем согласно давлениям,

найденным на момент времени t : соответственно (p,*z)cp it t и zcp it t-

Тогда приближенные значения давлений на забоях реальных скважин, согласно (5), можно вычислить по формуле

Здесь рс фг- (t + A t) — приближенное значение давления на

забое

фиктивной скважины в момент времени t + A t , вычисленное

в ре­

зультате решения задачи при отмеченном допущении.

 

По найденному приближенному значению р с ( t -f-

A t) согласно

уравнению технологического режима эксплуатации (4) уточняем средний дебит скважины в интервале времени [if, t + 2 A t ] по

формуле

Р а (t + A t )"

(7)

Здесь величины (p*z)Cj-i+A; и zcit+Ki вычислены с использованием

найденных значений давления р с t (t + A t) и р с>фг (t + At).

Уточненные дебиты скважин используются для отыскания нового приближения задачи (1)— (4) в момент времени t + 2 А t. Процесс

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

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

мости р с i

= рс i (t) и

= qt ( t). По известным значениям забойных

давлений

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

устьевые давления, давления на входе

в групповые сборные пункты, периоды бескомпрессорной и компрес­ сорной эксплуатации и т. д.

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

168


месторождения или выполнения некоторого экстремального кри­ терия.

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

Нетрудно заметить, что решение задачи (1)— (4) возможно и на электрических моделях. Порядок расчетов в этом случае на каждом временном слое аналогичен рассмотренному. При учете реальных свойств газа (и пористой среды) расчет добавочных сопротивлений проводится на основе уравнения притока газа к скважине (17) (см. § 8 данной главы). Уравнение технологического режима экс­ плуатации скважины при допустимом градиенте давления можно использовать в виде (4).

2. Т е х н о л о г и ч е с к и й р е ж и м э к с п л у а т а ц и и п р и п о д д е р ж а н и и п о с к в а ж и н а м д о п у с т и ­ м о й д е п р е с с и и н а п л а с т .

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

Рс. ф г- —Pc i = б/-

(8)

Таким образом, при поддержании по скважинам допустимых де­ прессий на пласт требуется решать задачу (1)— (3), (8).

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

ввиде (5).

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

приближенные значения давлений на забоях укрупненных скважин,

ас использованием (8) — приближенные значения давлений на за­ боях реальных скважин.

Уравнение (5) с учетом

(8)

записывается

для момента

времени

t + A t в виде:

 

 

 

 

 

 

 

 

 

бг (t -{- At) [ 2

ф

(t

At) — б (t —f—А4)] =

 

 

=

A t (|T*z)cp it+Atf, (f +

At) +

B tz cp it+Atq! {t +

At).

(9)

Решая уравнение (9)

относительно

qt ( t +

Д<):

 

 

 

q i ( t + A t ) =

 

A j

(p*z)Cp, i t+At

 

 

 

 

 

2Bizcp £(+д(

 

 

 

 

 

 

 

 

 

 

A i

(|X *z)cp i t+At

2

I б»

(t -j- At) (2дсф tPH-AQ

81 (t -|-At))

+

2 B iZcp it+At

 

 

 

 

BiZcpit+Af

 

(1 0 )

 

 

 

 

 

 

 

получаем уточненные значения дебитов скважин, соответствующие моменту времени t -(- A t (или, что то же, — рассматриваемому ин­ тервалу времени 2 A t , так как дебиты скважин скачком изменяются

169


в моменты времени t, t + 2 A t , t -j- 4A i, . . . и остаются постоянными в пределах каждого расчетного интервала времени 2Д t).

С уточненными значениями дебитов скважин вновь отыскивается решение задачи для момента времени t + 2A i и т. д. подобно преды­

дущему случаю. Сделанное ранее замечание относительно соответ­ ствующих расчетов на электрических моделях остается в силе и здесь.

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

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

емом для i-й скважины в виде:

 

Р а = Р а ' Л

(11)

интересующая нас задача формулируется уравнением (1) и усло­ виями (2), (3) и (И ).

Приближенное решение рассматриваемой задачи на временном слое t -г 2Дt отыскивается при задании по скважинам дебитов газа, соответствующих времени t. Для уточнения дебита i-й скважины

используется уравнение (5) при учете (11). Имеем

 

 

 

A j

(M-*z)cp й+Д{

 

 

 

 

 

2 S j ' z Cp tt+ b t

 

+

A i (p*z)cp

2

Р с. ф 1 (* + д<) — Р* 1Р + АО

(12)

2-B»zCp 1

(+д(

 

'

B i z Ср 1(+д(

 

 

Дальнейшие расчеты подобны предыдущим.

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

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