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

книги из ГПНТБ / Ширковский, А. И. Добыча и подземное хранение газа учебное пособие

.pdf
Скачиваний:
68
Добавлен:
20.10.2023
Размер:
8.32 Mб
Скачать

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

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

§ 44. Процессы, происходящие в нефтяной залежи во время хранения газа

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

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

Авторы работы [33] провели тщательные экспериментальные исследования на модели. Они исследовали влияние плотности неф­ ти, давления, температуры, отношения коэффициентов динамиче­ ской вязкости газа и нефти на процессы вытеснения и испарения газа. В качестве хранящегося (сухого) газа использовался метан с небольшим количеством этана, пропана и более тяжелых компо­ нентов (меньше 1% объемн.).

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

В экспериментах максимальный объем вытесненной нефти со­ ставлял 76,5%, испарившейся— 18,9% от начального объема то­ варной нефти в образце, максимальный объем добытой нефти со­ ставлял 93,8% от того же объема.

Увеличение угла наклона модели пласта с 15 до 30° к гори­

150

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

§ 45. Определение максимальной емкости подземного хранилища

Рассмотрим схему частично выработанного нефтяного место­ рождения, изображенную на рис. 35. Общий объем газа в храни­ лище складывается из трех частей: объема свободного газа в га­ зовой шапке; объема газа, рас­ творенного в остаточной нефти:

окклюдированного (рассеянного в форме отдельных пузырьков в

массе нефти) газа.

 

нефти

в

Массу остаточной

пласте

VQ выразим

как

разность

между

начальной

массой \ ’я

и

массой

добытой нефти

Уд

 

^0 = ^ 3 - ^ ,

Уз = FhmQ(l Sa) РBg,

(203) Рис. 35. Схематически» разрез нефтянон залежи массивного типа в кон-

(204)це разработки

где F — площадь нефтенасыщенного коллектора в м2; h — средне­ взвешенная по площади F нефтенасыщенная мощность в м; /п0 — коэффициент абсолютной пористости в долях единицы; S„ — объем

•связанной воды в долях единицы; р„ — плотность нефти при стан­ дартных условиях (р = 760 мм рт. ст. и ^=20°С).

При начальном пластовом давлении р„ и пластовой темпера­

туре t„ в 1 т нефти растворяется а м3 газа.

 

Масса газа

(в кг), растворенного в 1 т нефти

 

 

GP = aA -1,205 (273 +

20)

 

 

 

(273 +

/„)

 

где А — относительная плотность газа по воздуху.

в жид­

Кажущаяся плотность газа рк (при известных А и р„)

кой фазе можно определить по графикам работы [3].

 

Объем газа

(в м3) в единице объема жидкой фазы

 

 

Уг.ж = Ог/рк£.

 

(205)

Общий объем нефти (в м3),

насыщенной газом

 

 

Ун =

1+Уг . ж.

 

(206)

Общая масса этого объема

 

 

 

'

Gr ж =

Gr -f- рн£.

(207)

151

Плотность нефти, насыщенной газом (в кг/м3)

(208)

По графикам работы [17] можно найти поправку Арр на сжи­ маемость при пластовом давлении и поправку Др« на увеличение объема при пластовой температуре.

Реальная плотность нефти, насыщенной газом, в пластовых условиях

Р н г = Р ' . г +

Дрр — Ар/.

(209)

Пластовый коэффициент

 

 

 

ь==3

^ '

(2Ю)

 

Рн.г g

 

Объем порового пространства, занимаемый оставшейся нефтью

в пласте

 

 

 

О __

( У 3

У д ) Ь

(211)

ЫШК

 

Pllg

 

 

 

При повышении давления

с

р к до конечного

давления ртах

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

Удb

f Ртал

Рк \

(212)

п.,0

\

/ '

Рн8 \ 2тах

 

Объем растворенного в оставшейся нефти газа

(в м3)

Qр ;

(УзУд).а.

 

(213)

 

Pag

 

 

Объем газа (в м3), который будет

закачан в

газовую шапку

=

гтах

 

(214)

 

 

 

Общий объем газа (в м3), который можно закачать в частично выработанное нефтяное месторождение при постоянном объеме порового пространства

Qx= Q o + Qp+ Qm-

(215)

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

Ртах—85 кгс/см2;

а = 60 м3/м3;

*=16 °С;

р„ = 0,867

т/м3;

р„ = 0,38 т/м3;

Дрр= 0,005 т/м3;

ДР1 =0,001 т/м3;

У3 = 552000

т

или

637 000

м3; z mai= 0,77;

z„ = l; р„ = 5 кгс/см2; Qr=860 000 м3; Уд= 384 600 т; Д=0,635.

 

,

0Г = 60-1,205-0,635" (273 +

^ --

= 4 6 ,5 кг.

 

 

(273 +

16)

 

 

 

152

1 + ж = 46,5:360 = 0,129 м3; V„ = 1 + 0 ,1 2 9 = 1,129 м3;

рн г = 0,81 + 0,005 — 0,001 = 0,814 т/м3; Ь =

913,5

= 1,12;

Q0 = (637 000-1,12 — 216 500)^-

 

 

Qp =

(552 000 — 384 000)

• 60 =

11,6-Ю6 м3;

0,867

 

 

<3Ш=

860 000

=

94,9-106 м3;

Qx = (52,1 + 11,6 + 94,9). 106 =

158,6- 10е м3.

§46. Оценка производительности нагнетательно­ эксплуатационных скважин при отборе газа

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

Р (Рил + Рз)/^>

(216)

получим

(217)

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

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

Глава VI

ПОДЗЕМНОЕ ХРАНЕНИЕ ГАЗА В ЛОВУШКАХ ВОДОНАСЫЩЕННЫХ КОЛЛЕКТОРОВ

§ 47. Методы определения герметичности кровли ловушки, выбранной для создания подземного хранилища

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

Поэтому в процессе разведки и опытной закачки газа необ­ ходимо:

доказательство герметичности кровли ловушки; расчет коэффициента проницаемости водонасыщенного коллек­

тора; определение остаточной водонасыщенности при вытеснении воды

газом; измерение или вычисление объемной газонасыщенности обвод­

ненной зоны при отборе газа; определение продуктивных характеристик эксплуатационных

скважин; изучение прочности свойств газопасыщениого коллектора и

разработка мероприятий по укреплению призабойных зон скважин.

О п р е д е л е н и е г е р м е т и ч н о с т и к р о в л и

л о в у ш к и

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

154

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

в единице объема жидкости и состав

7

3

15 6

_____

растворенного в воде газа различны.

Если приведенные к одной плоскости

 

 

 

 

 

отсчета давления, солевой и газовый

 

 

 

 

 

составы одинаковы,

есть

основания

 

 

 

 

 

полагать, что эти пласты сообщаются

 

 

 

 

 

между собой.

 

 

 

 

 

 

 

 

 

 

Герметичность кровли должна быть

 

 

 

 

 

установлена над предполагаемой пло­

 

 

 

 

 

щадью газоносности создаваемого хра­

 

 

 

 

 

нилища (заштрихованная площадь на

 

 

 

 

 

рис. 36).

методы

пробных

отка­

 

 

 

 

 

 

Используя

 

 

 

 

 

чек (или закачек) жидкости из объек­

 

 

 

 

 

та

II последовательно

через скв. 1, 3,

 

 

 

 

 

5

фиксируют

изменения

положения

 

 

 

 

 

уровней жидкости в скв. 2, 4, 6, 7. При

 

 

 

 

 

этом обязательно регистрируют изме­

 

 

 

 

 

нения барометрического давления ат­

 

 

 

 

 

мосферы.

 

 

 

 

Рис. 36. Схематический разрез

 

Если скв. 2, 4, 6, 7 не реагируют на

изменение давления в скв.

1, 3, 5, мож­

и структурная карта

подзем­

ного хранилища

газа,

созда­

но

предполагать, что

кровля

непро­

ваемого

в ловушке водонасы­

ницаема для жидкости. Этот метод,

 

щенного

пласта

 

однако, не дает надежных результатов,

репрессии

незначительны,

поскольку создаваемые депрессии или

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

непроницаемость покрышки по воде, а не по газу.

 

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

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

Преимущества этого метода весьма существенны: 1) опреде­ ляется герметичность кровли для газа; 2) величина депрессии и репрессии может быть существенно больше, чем при откачках и закачках жидкости; 3) вследствие большого различия в вязко­ стях и плотностях газа и воды закачиваемый газ будет насыщать небольшую мощность пласта и распространяться на значительное расстояние по площади, при этом можно уменьшить число наблю­ дательных скважин, вскрывших пласт I, для определения герме­ тичности покрышки, а также сократится время для проведения исследований.

155

При закачке и отборе воздуха в скв. 1, 3 и 5 (последова­ тельно) фиксируют изменение давления (уровня) в скв. 2, 4, 6, 7. Если скв. 2, 4, 6, 7 не реагирует на изменение давления в пла­ сте II, покрышка ловушки считается герметичной.

§ 48. Определение объемной газонасыщенности обводненной зоны при отборе газа

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

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

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

уравнение баланса газа в залежи

 

 

 

Qr + a (QoQr) -ggL =

i£ ? ~

,

(218)

zBpB

Pb

 

 

где Qr и Qo — соответственно текущий

и начальный объемы поро­

вого пространства газонасыщенной части залежи

в м3; а — коэф­

фициент объемной газонасыщенности обводненной зоны в долях от эффективной пористости; pD/zB— средневзвешенное приведенное давление в обводненной части залежи в кгс/см2; pBfzB— средне­ взвешенное приведенное давление в необводненной части залежи

в кгс/см2; Q3 и Qa— соответственно

начальные запасы газа и

объем отобранного газа из залежи,

приведенные

к давлению

1 кгс/см2 и пластовой температуре, в м3.

 

Объем вторгшейся воды в залежь

 

 

Q „^Q 0— -(^ ~ Qa)PaZB .

(219)

Рв2а

 

 

Объем добытого газа из обводненной части залежи за счет вторжения воды

156

( 220)

где pajza— начальное приведенное давление в залежи в кгс/см2. Объем оставшегося газа в обводненной зоне

Q0 = QB ^ - - A Q a = a (Q0- Q r) ^ .

(221)

Объем чисто газонасыщенной части залежи

Qr = (Q3 — Од) zb

VQb zq AQff ) 2b

(222)

P b

 

P b

 

Коэффициент объемной

газонасыщенности обводненной

зоны

( qb — -AQz ) z*

(223)

а = Д----

?°--------

J— .

 

(20 ^г) Рв

 

Коэффициент газоотдачи обводненной зоны

 

p =

i

а Рв2о

(224)

 

 

*вРо

 

Газ, остающийся в обводненной зоне, является составной частью буферного газа. Величина коэффициента объемной газонасыщен­ ности обводненной зоны и масса газа в этом объеме существенно зависят от геолого-физических параметров пласта, темпов отбора газа из хранилища, физических свойств воды и газа, литологии газонасыщенного коллектора (т. е. типа горных пород, слагаю­ щих коллектор: пески, песчаники, известняки, доломиты и др.) [28, 31].

Пример 22. Определить коэффициент остаточной объемной газонасыщенности

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

газоотдачи Аманакского

подземного

хранилища

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

исходных данных:

Ро

= 30,3 кгс/см2;

Рв

кгс/см2; Qo=9,95 106 м3;

— = 15,6

 

 

 

го

 

 

zB

 

 

Од = 208,1 • 10е м3; ~г~ = 22,9 кгс/см2;

 

=25 кгс/см2.

 

 

 

 

гв

 

 

 

гво

 

 

 

 

(9,95.106-30,3 — 208,1 -10е)

м3.

 

Ов =

9,9510е — ----------------——----------------- = 4 млн.

 

 

 

 

 

 

15,6

 

 

 

Д0д =

208,1 — 9,95.10е (30,3— 15,6) =

61,9-Ю3 м3.

 

 

0о =

4-10в.30,3 — 61,9-10е = 59,3-10е м3.

 

Л

(9,95-103-30,3 — 20 8 ,Ы 0 3)

(4-103.30,3 — 61,9-10»

 

в

, - -----------------

-------------- ------------------- ( 5 J — =-----------2.I6-10-

 

 

 

 

59,3-10е

 

 

 

 

 

 

 

 

= 0,332.

 

 

 

 

 

(9,95 — 2,15) -10е- 22,9

 

 

 

 

 

 

Р = 1 - 0 ,3 3 2 -

 

22,9

 

 

 

 

 

 

 

0,75.

 

 

 

 

 

 

 

30,3

 

 

157

§ 49. Закачка газа в купольную часть (ловушку) водонапорной системы в условиях упругого режима

П е р в а я ф а з а н е у с т а н о в и в ш е й с я ф и л ь т р а ц и и

Гидродинамическая схема для расчета приведена на рис. 37. Изменение давления в хранилище при закачке газа в пласт с

постоянным темпом можно определить по формуле [26]

Pi'

Рп

 

 

 

Q0

 

Р г ~

QrPat - h p 0Q 0

 

 

=

0.

 

(225)

 

 

 

2x2

 

 

 

2xl

 

 

 

 

 

 

 

 

 

4Y,t

л

 

/

*5

4y.t

 

 

 

 

 

 

 

 

 

In

 

 

 

In

 

 

 

 

 

 

 

 

 

R0

 

 

 

Ro

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Зная

P r = p r ( t ) ,

легко

рассчитать

Q= Q(/),

а по графику Q=

= Q(h),

построенному по структурной

карте ловушки,

 

определить

A = h(t).

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Пример 23. Определить, как

 

 

 

 

 

 

 

 

изменились в течение первой фазы

 

 

 

 

 

 

 

 

неустановившейся

 

 

фильтрации

 

 

 

 

 

 

 

 

объем

газоносной

части

 

пласта

 

 

 

 

 

 

 

 

Q= Q (2),

мощность

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

 

 

 

 

 

 

 

 

h = h (2),

давление

 

в

хранилище

 

 

 

 

 

 

 

 

P r - P r ( t ) ,

расход

жидкости

q B =

 

 

 

 

 

 

 

 

= q B

(t)

при

создании

хранилища

 

 

 

 

 

 

 

 

газа

в

ловушке

пластовой

водо­

 

 

 

 

 

 

 

 

напорной системы (при упругом

 

 

 

 

 

 

 

 

режиме и постоянном расходе за­

 

 

 

 

 

 

 

 

качиваемого

газа

рг const).

 

 

 

 

 

 

 

 

 

щие

Для расчета приняты следую­

 

 

 

 

 

 

 

 

исходные

данные:

 

т = 0,2;

 

 

 

 

 

 

 

 

£=1,5Д ;

р в = 1

сПз;

мощность

 

 

 

 

 

 

 

 

пласта

ht= 20 м;

рп = 60

 

кгс/см2;

 

 

 

 

 

 

 

 

число

нагнетательных

скважин

 

 

 

 

 

 

 

 

п0 = 20;

радиус окружности, на ко­

 

 

 

 

 

 

 

 

торой

расположены

нагнетатель­

 

 

 

 

 

 

 

 

ные

скважины,

У?б =

500

м;

р„ =

Рис.

37. Гидродинамическая

схема

водо-

= 1000

кг/м3;

высота

ловушки

Н = 5 5

м;

средний

радиус

отступ­

«асыщеыиого

пласта

 

при расчете закачки

ления

воды,

принимаемый

посто­

газа

в

процессе

создания

хранилища:

янным,

/?о= 500 м; постоянный рас­

 

 

 

г ,

До + *о

 

 

ход

газа

qB—10е м3/сут; коэффи-

 

 

 

*'0 —

 

2

 

 

циент объемного

упругого

сжатия

 

 

 

 

 

 

 

 

жидкости

Р в = 4,5 5

1 0 ~ 5

 

см2/кгс;

сжатия пород

рс = 1,55

10~5

см2/кгс;

 

коэффициент

объемного

упругого

радиус контура

пластовой

 

водонапорной

системы /?к = 19,47 км.

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Схема пластовой водонапорной системы приведены на рис. 37, а кривые за­

висимостей П= 12 (/г),

Fm = Fm(h) представлены на рис. 38.

 

 

 

 

 

 

 

 

 

Результаты вычислений при Й0=0

приведены в табл.

30.

 

 

 

 

 

 

 

158

Т а б л и ц а 30

Основные данные, характеризующие процесс создания хранилища газа

в куполообразной ловушке пластовой водонапорной системы

в условиях

упорного

 

 

 

режима

 

 

 

 

t , сут

R , Юз м

Р Г , кгс/см2

Д, 10' м '

h , м

Чв,

м '/сут

Q, 10е м»

30

7,95

81,0

0,371

3,5

|12 400

30

60

11,25

83,0

0,723

5,0

12 050

60

90

13,75

84,0

1,07

6,0

11 800

90

120

15,8

84,8

1,41

6,8

11 720

120

150

17,8

85,5

1,75

7,6

11 630

150

180

19,47

86,0

2,09

8,35

 

11 600

180

По данным табл.

30

построены кривые зависимостей p r = P r ( t) ,

h = h ( t )

и Q= Q(t),

приведенные на рис.

39.

Из данных таблицы

 

 

 

 

 

е,Аз

 

 

 

ругаем Л,м Я,10'V

 

6

8

9,1о У

 

/I (раза Е (раза

О

 

 

100 -10 5

Т

 

 

80 - 8 -0 Л

О

 

 

 

 

во - В ■3

А1

8

 

 

^ 9

 

00 - 0

Z

 

12

-

 

Гт\

 

20 -1

-1

 

16

 

 

 

 

 

 

fl,М

Рис. 38. Кривые зависимостей Q и Fin от газонасыщеиной мощ­ ности пласта h

0 00 120 180 ZOO 300 t,cупг

Рис.

39.

Кривые

зависимостей

Q=

= Q (0 >

h — h ( t )

и P r = p r {t) при

за­

качке

газа

в ловушку с постоянным

темпом в течение первой и второй фазнеустановившейся фильтрации воды

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

В т о р а я ф а з а н е у с т а н о в и в ш е й с я ф и л ь т р а ц и и ж и д к о с т и при з а д а н н о м во в р е м е н и р а с х о д е з а к а ч и в а е м о г о г а з а q v = q T( t )

Пластовая водонапорная система рассматривается как замкну­ тый подземный резервуар [26].

159

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