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

книги из ГПНТБ / Закиров, С. Н. Проектирование и разработка газовых месторождений учебное пособие

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

деби т ск важ и н ы . В этот период р або та я депрессия на вто р о й пласт

все врем я

возрастала,

что объясн яет увеличение притока газа и з это го п ласта .

С момента

д ости ж ен и я

равенства

 

P i { t ) = P i (t)

к он троли рую щ и м показатели разработки стал второй п ласт.

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

Взаключение отметим, что предлагавшиеся ранее аналитические

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

§ 3. Определение показателей разработки многопластового месторождения единой сеткой скважин при водонапорном режиме

Предполагаем, что известны все необходимые для расчетов исход­ ные геолого-промысловые данные. Задан суммарный отбор газа из месторождения во времени Q = Q (£). Требуется определить основ­ ные показатели разработки месторождения.

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

вопределении изменения во времени средних пластовых давлений

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

Каждую залежь аппроксимируем равновеликой укрупненной скважиной. Расчеты будем проводить по шагам. Если решение за­

дачи на некоторый момент времени t At известно, то, согласно

§4 главы VI, для определения изменения дебита воды, поступающей

вкаждую укрупненную скважину спустя шаг по времени At, имеем следующие расчетные соотношения:

A?b i (0 = -

а1 1

/

a i .

С1 .

(1)

2b i

 

г щ

b i

 

 

 

Адь 2 (t) =

2Ъ2

-

Т /

Л

С2

 

(2)

 

У

461

ь 2

 

14*

211

Здесь

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

+

"‘ 1% Х л

 

ы)~

 

 

 

 

т я к р С°‘ -

|0*»-.) +

+

Ш

£

2

' 4?“

A

(,0‘ ~ V

. ) + м

W А1-

а

2 л г к ln w

 

;

 

 

 

Ь1 =

. | £ ^ р

 

/fOl_ f 0,

. ) +

 

1Д<

 

Лз, .

 

 

 

 

 

 

 

o- r - r -ln

/?х (0 ’

 

 

 

 

 

 

 

1

2nk1h1

 

1 \

1

 

 

n - l/ i

/квс x&iф-!

 

 

 

 

 

 

 

 

 

 

 

n-1

 

 

 

 

 

 

 

 

 

 

 

 

 

£. - f ■ А -

 

2

 

i=l

А?“ ‘ A

 

- V . ) - А

 

 

 

X

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Х

1п

 

 

L lp B

l V l (<) _ (

 

 

- Р а т < ? д о б 1 ( о ) 2 1 л (01;

 

 

 

 

Li = аД , 1 —(QB1(t —At) + qB1(t —At) At)-,

ioi=BM --

 

 

=

 

г)

А /

I .

P B 2

2 (t — At)

,

 

Д з2

J

Дв 2

P

/ f n

f n

\

i

г

. Й 2 Щ +

 

2п*в 2й2

 

1

 

Л2 (0

L

j 2nfc2/i2

2 (to2

t02n-i)

+

 

 

 

n-l

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Я3a

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

>

^

2я Л2Л2

2

l A?B2 j ^ 2 (f ° 2

f ° 2i - i )

+ PbJ/2 ^

 

2 22JkBя*в\h22n2ln д 2 (г)

 

 

 

j-x

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Дв 2 Д г

j n

^ з 2

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

2я/св 2й2

 

я 2 (г) »

 

 

 

С2 “ ^

2L 2

L2 2nfc2\ 2 2

 

А ?в 2

 

( f ° 2 ~

f ° 2i - l ) _

L 2 Цв

 

Аг)

X

 

 

 

 

 

 

 

 

i=i

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

X^n ~И~фТ)

-^гРв 2^2(0

\

РН2а2^Н 2

' PaiQдоб 2 ( 0 ^ ^ [Рг (0 ] >

 

 

 

 

д 2 (г)

"

 

 

w

 

ZH 2

 

 

 

 

L2=~a2QH2- ( Q B2( t - A t ) + q B2(t-At) At)-,

fo2- - g f - .

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

лз2

 

 

При расчетах по формулам (1) и (2) в первом приближении при­ нимаем:

2[рх (*)] ~

2[р[ (t - At)]-,

z [p2(f)]« 2[p2 (t -

Д<)];

 

ДО; V i ^ ^ y d t - A t y

(3)

Rx{ t ) ^ R ^ t - A t y

R2(t) ~ R 2(t - A ty

<?доб 1(t) *=><?доб! («— At) + n (t At) qx(tAt) At;

<?Добг(0 ^

<?доб 2 (t At) + n (t — At) q2 (tAt) At.

212

Тогда, аналогично описанному в § 4 главы VI,

определяем при­

ближенные значения AqBl (t),

AqB2 (t), QBl (t), Qb2 (t), yx (t), y2 (t),

P i(t),p 2 (t),

R2 (t).

Найденные значения

коэффициентов

сверхсжимаемости

газа, у г (t) и у2 (t), R 1 (t) и й 2

(t) используются

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

тации находим значение притока газа из

 

первого

пласта:

 

 

Y

ш + -

[гр п о -а п

(4)

 

в1

 

 

 

 

 

 

 

 

Далее вычисляем приближенное значение рс2 (t),

которое считаем

равным Рсх (t). Следовательно,

имеем

 

 

 

 

 

?*(0 = — А*

 

2

,

PUt)—Pct(t)

 

Y

т

1

 

В2

^

2В2 ' V

\ 2 В 2 ;

 

Приближенное значение дебита одной средней скважины в мо­

мент времени t равняется

 

 

 

 

 

 

 

q ( t ) =

q1 (t) + q2 (t),

 

 

 

( 6 )

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

n{t)

Q{t)

(7)

g(t)

 

 

Полученные показатели разработки позволяют уточнить добытые количества газа из первого и второго пластов на момент времени t по следующим формулам:

С?доб 1 (0 — С?доб1

^0

п (t At) qi (t At)-\-n (t) qi (t)

At\

(8)

2

Q ДОб 2 (0 — <?Д06 2 ^

I

n ( t - A t ) q 2 {t — At)-\-n(t)q2 (t)

 

(9)

ДО

 

 

Найденные уточненные величины (3) используются в расчетах второго приближения по формулам (1) и (2). Как обычно, процесс итераций продолжается до совпадения результатов последнего и предпоследнего приближений (по одному из параметров) с погреш­ ностью не более е.

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

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

213

Таблица 20

 

 

 

П о ка зател и

р азр аботк и м есторож дения

С еди н ой сеткой скваж и н в

у сл о в и я х в од о н ап ор н о го реж им а

 

 

 

 

 

 

 

 

 

 

 

Годы разработки

 

 

 

 

 

 

Показатели

1-й

2-й

3-й

4-й

5-й

6-й

7-й

8-й

9-й

10-й

11-й

12-й

13-й

14-й

 

 

 

 

 

QВ 1

(t ),

м лн . м 3

0,75

1,80

3,9

6,5

10,61

15,11

21

27,2

34,8

42,7

51,8

61,3

71,9

82,9

<?В2

(t),

МЛН. М3

0,42

1,06

2,2

3,9

6,4

9,4

13,1

17,3

22,3

27,7

33,9

40,6

47,9

55,8

Л х ,

км

 

 

 

12,744

12,376

12,720

12,699

12,668

12,633

12,587

12,538

12,479 12,416 12,344 12,268 12,183 12,094

Л 2,

км

 

 

 

12,576

12,571

12,560

12,546

12,253

12,501

12,464

12,433

12,392 12,345 12,293 12,235 12,172

12,103

P i,

КГС/СМ3

 

297,3

293,8

287,6

279,1

269,2

259,3

250,2

241,6

233.8

226,3

210,0

211,9

205,4

199,0

р 2 ,

кгс/см ®

 

297,8

295,1

289,8

282,3

273,3

264

255,3

247,0

239,3

232,0

224,8

217,9

211,4

205,2

(71 П ),

ты с.

m 3/ c v t

959

933

898

852

818

781

751

724

700

678

657

636

620

598

(7я (г),

ТЫС.

m 3/ c v t

711

705

694

678

659

639

620

602

585

568

552

536

522

507

<?доб (*).

м л р д , м3

5,7

14,7

25,8

40,1

57,3

74,5

91,7

109

126

143

160,5

178

195

212

ть

 

 

 

 

10

16

20

26

32

34

35

36

37

38

40

41

41

43

 

 

 

 

 

 

 

 

 

 

 

Продолжение

табл. 20

 

Показатели

15-й

16-Й

17-й

18-й

19-й

20-й

21-й

22-й

23-й

24-й

25-й

<?В1 ( г ) .

МЛН.

М3

94,7

107,0

120,0

133,4

147,3

161,6

176,2

191

206

222

237

Qв2 ( г ),

м лн .

м 3

64,3

73,4

83,0

93,3

104,4

116,7

130,1

144

159

174

190

Лх

(1),

км

 

11,997

11,896

11,788

11,676

11,558

11,436

11,310

11,179

11,044

10,305

10,763

Л 2 ( г ),

км

 

12,029

11,949

11,864

11,713

11,673

11,561

11,439

11,307

11,167

11,020

10,870

P l ,

к г с /с м 3

 

192,8

186,9

181,3

176,0

171,0

166,2

161,6

157,3

153,1

149,0

145,1

Р 2 ,

КГС/см 2

 

199,1

193,3

187,9

182,7

177,9

173,4

169,2

165,3

161,5

157,8

154,2

(7i

( г ), т ы с . м 3/ с у т

579

561

543

527

509

492

475

458

440

426

411

(72

( t),

т ы с .

м 3/ с у т

492

479

466

453

442

431

421

411

402

393

384

<?доб ( 0 . млрд. М3

229

245,5

261,5

277,1

292,2

306,8

321,1

334,8

348,2

361.2

373,8

п

 

 

 

44

44

44

44

44

44

44

44

44

44

44

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

Q(t) = n(t)q(t)

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

Q = Q(t).

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

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

Р и с. 6 6 . И зменение во времени показателей разработки м есторож дения С в у с ­ л о в и я х водонапорного реж им а (вариант с р и1 — рН2 = 300 к гс /с м 2)

рис. 66 через pi обозначено среднее пластовое давление при водо­ напорном режиме в месторождении А, имеющем одинаковые пара­ метры с первым пластом месторождения С. Существенные различия

в зависимостях р х = р х (t) и pi = pi (t) связаны с разными темпами

разработки первого пласта и месторождения А.

 

 

 

Дополнительные по сравнению с приведенными в предыдущем

параграфе исходные данные таковы: кг = 1 Д; к2 =

0,7

Д; кв1 =

=

0,2 Д;

fcB2 = 0,15 Д; Rsl =

12,75 км; i?32 = 12,58 км;

= 15 м;

h,

= 14

м;

т1 = тп2 = 0,18;

Рш = рнг =

300 кгс/см2.

двупластового

 

На рис.

67 приведены показатели

разработки

месторождения С, когда рн1 — 300 кгс/см2, рнг — 250 кгс/см2.

215

До достижения равенства пластовых давлений р г (t) и р 2 (t) допустимая депрессия на первый пласт определяла характер зави­ симостей изменения во времени основных показателей разработки месторождения С. Пока выполнялось неравенство

Pi (t)>Pa(t),

возрастала депрессия на второй пласт, увеличивался дебит притека­ ющего газа из второго пласта. В начальные моменты времени (до середины третьего года разработки) q2 (t) < 0 . Это привело к неко­ торому увеличению пластового давления во второй залежи и соответ­ ственно— к ее расширению (QK2 (t) -< 0). Интересно отметить, что Qn2 (t) стало больше нуля только на восьмом году разработки.

водонапорного реж има (вар и ан т с р И1 = 3 0 0 к гс /с м 2, р Н2 = 2 5 0 к гс /сы 2)

В последующие моменты

времени, при р 2 (t)

Допусти­

мая депрессия на второй

пласт

стала определять

изменение

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

месторождения С.

 

§ 4. Особенности формулирования задач разработки многопластовых месторождений при газодинамической связи между пластами

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

216

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

пластами Сообщение между пластами может возникать вследствие непо­

средственного их контакта (рис. 68). Решение задач разработки при наличии подобной связи между пластами в настоящее время затруднительно. Такое решение связано с необходимостью форму­ лирования задач фильтрации в трехмерном пространстве для каждого пласта, что требует значительной (и трудно определимой) информа­ ции об изменении параметров пластов в направлении осей х, у и z.

Р ис. 6 8 . С хем а непосредственного

к о н ­

Р и с.

6 9 . П родукти вн ы е пласты / и

так та п родук ти вн ы х п ластов I

и I I

I I ,

разделенны е слабопроницаемой

 

 

 

перемы чкой

На границе между пластами должны соблюдаться условия непрерыв­

ности потока

 

^2

др2

 

И dpi

 

[I

дп

дп

' '

н неразрывности давления

 

 

 

 

Pi {х,

У) = Р2(х, У).

(2)

Здесь к — коэффициент проницаемости в направлении нормали п; величины с индексом 1 относятся к нижнему пласту, с индексом 2 — к верхнему.

Решению задач подземной газогидродинамики при соблюдении условий типа (1) и (2) посвящено сравнительно небольшое число исследований (М. Т. Абасов, 3. Б. Адигамов, К. Н. Джалилов, В. А. Карпычев, Ж. Дюво, М. Дюпюи). Условия (1) и (2) приходится учитывать и при решении плоских задач, когда параметры пласта по площади или вдоль координаты х или у (одномерные задачи) изменяются скачкообразно. Исследование таких задач описывается в работах Г. И. Баренблатта, Г. Г. Вахитова, Г. Л. Говоровой, М. А. Гусейн-Заде, М. М. Саттарова, В. Н. Щелкачева. Вследствие близости названных задач методы решения одних можно исполь­ зовать применительно к другим.

1

Е сли два пласта

м ногопластового м есторож дения

им ею т р а зл и ч н ы е 'о т ­

м етки

гл у б и н залегания

Г В К , то он и изолированы д р у г

о т д р у га . П ри едином

Г В К

дл я д в у х пластов следует им еть в в и д у , что они п редставляю т собой единое

газодинам ическое

ц е л о е ,и р азр аб отк а и х будет соп р овож даться соответствую ­

щ ими

обменными

процессам и .

2 1 7

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

 

 

информации о коллекторских свойствах пла­

 

 

стов. Наличие газодинамической связи меж­

 

 

ду пластами

учитывается

соответствующей

 

 

записью дифференциальных уравнений не-

 

 

установившейся фильтрации газа (нефти).

 

 

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

 

 

коллекторским свойствам,

разделены слабо­

 

 

проницаемой

перемычкой.

Параметры пере­

 

 

мычки будем

отмечать индексом «па. Выде­

 

 

лим в пласте элементарный объем, как пока­

 

 

зано на рис. 70. Вывод дифференциальных

 

 

уравнений дадим по схеме Л. Д. Воронковой.

 

 

Изменение

массы газа

в нижнем пласте

 

 

за счет фильтрационных потоков вдоль осей

Р ис.

7 0 . Элементарный

х жу равно (см. § 5 главы И)

объем п р и наличии меж ­

— dfaT" d x d y d t — dpiji dxdydt,

д у

пластами слабопро-

ницаемой перемы чки

 

 

(3)

 

 

 

 

Предполагаем, что газ из верхнего пласта перетекает в нижний (в пределах рассматриваемого элементарного объема). Тогда масса притекающего газа (с плотностью р) за время dt равна

dx dy dt.

Накопление газа в элементе нижнего пласта с учетом перетока из верхнего пласта равно сумме (3) и (4):

dP!»l

h

dpl^l

An

P 2 — Pi

1 dx

1

dy

p

An

Вместе с тем изменение массы в рассматриваемом элементе пласта равно

a1m1h1 ~ ~ dx dy dt.

( 6)

Приравнивая (5) и (6), получаем уравнение неразрывности для нижнего пласта

An P2 Pi

( 7)

P An

218

Соответственно

для

верхнего пласта с

учетом перетока

газа

из него в нижний пласт получаем

 

 

- к

др2ц2

7.

9p2v2

d2pb2h/2 dt

Pi — Pi

(8)

дх

>h

ду

hn

Подстановка в (7) и (8) выражений для компонент векторов ско­ рости фильтрации и уравнения состояния для идеального газа дает

д*р\ дх2

52Р^

Лг2

|, 32pf

'0Р2

,

О ^3 tots 1

 

to

 

)

1

дуЪ

2 a i«iiP

dpi

2kn

p

л );

ki

dt

hihn (Pi

2сс2т гР

dpi

\ 2kn

p

 

к2

dt

1 k%

( P i ■- P i ) -

(9)

(10)

Понимая под р полусумму давлений в верхнем и нижнем пластах

(в точке с координатами х

\

 

р —

Pi "I- Рч

,

имеем

 

и у), т. е.

 

 

 

£2р|

 

Срр\ _ _

2 а х т х у

др1

 

 

к„

.

2

а ,

(И )

дх%

 

дуЪ

ki

dt

 

 

 

 

 

” 1' 1

 

 

 

 

 

 

 

d2p l

,

d2Pl

2а2т 2р

др2

.

 

кп

,

2

 

(12)

дх*

I"

ду 2

/с2

dt

k2h2hn

' " 2

" Х '

 

Уравнения (И )—(12) — искомые дифференциальные уравнения, описывающие фильтрационные процессы в пластах, разделенных слабопроницаемой перемычкой.

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

У ) К ( х , У ) ^ ~ \ + - щ ;[ } ^ х , У Ж { х , У ) у у ] =

 

= 2\ia1m1(х, y)hх(х, у)

 

(Pi —Р\У,

(На)

у) h2(x,

У )-^ * -]+ ^ -[* 2(®.

У)К(х, ? /)у у ] =

 

= 2\ш2тп2(х,

y)h2(x, у)

+

I) (Р«— Р*0-

(12а>

При учете реальных свойств газа и пластов соответствующая

система

уравнений принимает

вид:

 

 

 

 

'

к г (я, У, P i ) hi (я, у)

d p i 1

.

д

Г кг (х, у , P l ) fex (х,

у)

др\

- Г -

р(P l )

z (P i)

дх J

 

ду

L

Р (P i) z2 (P i)

 

ду

дх L

 

 

= 2а1(х, у)т1 (х,

y)h1(x,

y ) 4 r [ l j ^ ' ] '

кп (я, У, Р)

у) (р5- - pl);

Р (Р) z (р) hn (х,

 

 

 

 

 

 

 

 

 

(И б)

 

к 2 (я, У, Pi) h2 (х, у)

др\ ~]

,

д

Г к2 (я, у, р 2) h2 (я ,

у)

др?

дх I

Р (Рг) z (р 2)

дх

 

 

 

Р (Рг) z (р 2)

 

др ] -

■2а2 (х, у) т2 (х, у) h2 (х, у)

 

 

-

(я, у, р)

 

 

 

 

Р (р) z (р) hn (я, г/) Ы - р\).

 

 

 

 

 

 

 

 

 

(126)

219

Здесь

р(х, у, t) = T (p1(x, у, t) + p 2(x, у, t)).

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

В последнее время решению задач стационарной и нестационар­ ной фильтрации жидкости и газа в пластах, разделенных слабо­ проницаемой перемычкой, посвящается значительное число исследо­ ваний (М. Т. Абасов, А. Арсланов, М. А. Гусейн-Заде, Г. П. Гусей­ нов, К. Н. Джалилов, А. М. Кулиев, В. И. Мотяков, Н. Мухитди­ нов, П. Т. Шмыгля, В. Н. Щелкачев и др.). Достаточно обширную библиографию см. в работах [14, 15].

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

§ 5. Определение показателей разработки при наличии газодинамической связи

между пластами в условиях газового режима

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

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

220

Соседние файлы в папке книги из ГПНТБ