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

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

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

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

Добавлен: 15.10.2024

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

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

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

а уравнение материального баланса для газовой залежи соответ­ ственно запишется в виде:

V-QuPuZа

аЙ (0 р (t) гат

(9)

 

Рат/(?ДОб (О-

 

2 (Р)

 

Здесь

1

p(t)-. J pa dQ.

ай (t) £2(0

Нетрудно видеть, что дифференциальное уравнение истощения залежи при водонапорном режиме имеет вид:

<?*(0 =

2ат

А

' р (г) ай (г)'

(Ю)

 

Part

dt

. 2 (р) .

 

Уравнения (9) и (10) получены в предположении, что вода

пол­

ностью вытесняет газ из обводняемого объема пласта. На

самом

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

остается «защемленный» газ.

 

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

Так как обводненный объем пласта равен QH— Q (t), то в этом объеме при среднем значении коэффициента остаточной газонасыщенности аост находится газ в количестве

Мобв(0 = рат[QH- Q(t )]ао с

т ■- .

(И)

2 [ Л

(г)] Рат

 

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

аЙнЯн^ат__ ай (() р (t) zaT

| ^ fQ* б (t) -}-

 

+ [QH- Q (t)] aOCT (pB) - ^ 2 1 .

(12)

 

2[pB(0]

 

Здесь pB— среднее давление в

обводненном объеме

пласта;

z (ре) — коэффициент сверхсжимаемости при рв и Тпл; аост — отно­

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

При среднем значении коэффициента остаточной газонасыщен-

ности аост (рв) суммарное количество воды QB(t), поступающей в залежь к некоторому моменту времени t, распределится в объеме

31


- — ^ „ ■■■. Тогда газонасыщенный объем (внутри контура газ —

аКост (Рв)

вода) пласта ко времени t составит

aQ (t) = a Q „ -

 

Qb (t)

(13)

a

ct0CT {Pq)

 

 

Таким образом, под текущим газонасыщенным объемом в (12) понимается его выражение согласно (13).

Из уравнения материального баланса (12) вытекает следующее дифференциальное уравнение истощения залежи:

Q * ( t ) =

г ат

А

кй (t) Р (ty

2ат d_

а0сТ(Рв)[Йв-Я (01

Part

dt

2 (р) .

Part dt

(14)

Уравнения (12) и (14) (с учетом (13)) более строгие, чем соответ­ ствующие уравнения (9) и (10). Принципиальных затруднений для использования (12), (14) и (13) при определении показателей разра­ ботки газовых месторождений в условиях водонапорного режима не имеется [71]. Однако использование указанных формул приводит к усложнению методики расчетов, что объясняется необходимостью

определения величины рБ(t) и изменения значения коэффициента остаточной газонасыщенности. Кроме того, при анализе фактических

данных затрудняется определение зависимости рв = рв (t). Расчеты значительно упрощаются, если в (12) принять

pB(t)^p(t).

(15)

Условие (15) характеризует допущение о том, что газ защемляется

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

в залежи,

и изменение коэффициента остаточной газонасыщенности опреде­ ляется изменением во времени среднего пластового давления. В этом

случае из (12)

с учетом (15) получаем 1

 

 

.[Я (* )](-^ -Л * < ? д о б (о )

 

 

aQ(f)+-[Q„-Q(0]«ocT

 

или с учетом

(13)

 

 

 

 

 

<‘ Ч

Соответствующее дифференциальное уравнение истощения газо­

вой залежи записывается в виде (10), при этомосй (t) =

ocQH— QB(t).

1 Значение

коэффициента zar близко

к единице [26].

Поэтому здесь

и в дальнейшем для простоты принимаем zaT =

1. Примемздесь и вдальнейшем,

что дебиты, отборы и т. д. приведены к атмосферному давлению и пластовой температуре.

32


Важность полученных частных уравнений (17) и

(10)

состоит

в том, что для использования их благодаря допущению (15)

не тре­

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

обводненной

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

гать десятков процентов.

Проведенные в последние годы, в основном лабораторные, исследованпя

обнаружили влияние следующих

факторов

на

особенности фильтрации газа

в пласте: сорбционных процессов

[49 и др.];

деформации скелета горных пород

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

[18, 53]; наличия остаточной

нефти в пласте.

 

 

 

Уравнение материального баланса для залежи в условиях газового режима при учете десорбции газа в процессе снижения пластового давления может быть записано в виде:

а^Рп =

a ^ t ) + р ^

б (/) _

{t)

(18)

ZH

г (/;)

 

 

 

Здесь <?д (г) — суммарное количество десорбированного газа ко времени t,

приведенное к р ат и ТПл,

в м*.

 

 

 

Опубликованные в литературе сведения относительно адсорбции н десорбции газа пористой средой при существующих способах интерпретации результатов экспериментальных исследований малопригодны применительно к уравнению (18). Для наших целей наиболее подходящими являются результаты исследова­ ний А. М. Кулиева, А. И. Алиева и Э. В. Григоряна (1971) по адсорбционной способности образца породы по отношению к сухому газу. Содержание метана в газе составило 93%. Адсорбция проводилась при повышении давления до 20 кгс/см2. После этого проводилась десорбция при снижении давления до

атмосферного.

указанных экспериментов показала, что

 

Обработка результатов

 

 

“ =

0,098/4-0,41 ’

(19)

где а — количество газа в

см3/г,

адсорбированного па поверхности

образца

породы массой 1 г.

При давлении более 3 кгс/см2 формула (19) дает отклонение от эксперимен­ тальных данных не более чем на 7,9% .

Таким образом, можем записать, что

QnPn

P(t)

Qa(0108

( 20)

0,098p (t) + 0,41

Здесь Йп — геометрический объем

залежи в м3; рп — плотность породы

в г/м3; ап — количество газа, адсорбированного на поверхности образца породы

массой 1 г при начальном давлении р-л в см3/г.

залежь с площадью газонос­

Для

примера рассмотрпм гппотетическую

ности 50

км2, мощностью 100 м п коэффициентом пористости,

равным 0,3. Кол­

лектором является песчаник с рп = 2,3 г/см3;

рн = 100 кгс/см2.

Результаты расчетов по формулам (18) и (20) при отмеченных исходных

данных

показывают, что зависимость —

=

/ [<?доб (01

ПРИ десорбции

 

z(p)

 

 

 

3 Заказ 1013

33


располагается выше аналогичной зависимости без учета десорбции газа. Для принятых условий определяем начальные запасы, которые получаются на 8% больше запасов, подсчитанных без учета явления десорбции. Этот процент может быть и иным. В рассмотренном примере он оказался завышенным. При наличии связанной воды количество адсорбированного газа будет меньше полученного в результате расчетов по формуле (20). Кроме того, принято, что уравнение (19) справедливо для больших давлений, чем в проведенных экспериментах, которые следует рассматривать пока как частные.

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

материального баланса записывается в виде (при

принятии а =

1):

P(t)

ехр[ ат(ри— P(t))] =

Рн

Рат<?доб (0

(21)

2 (Р)

Zh

Пн

 

 

 

 

где ат — коэффициент сжимаемости пор в

1/(кгс/см2).

 

p/z(p), кгс/смг

Рис. 9. Зависимости p/z (р) = f [@доб (01 для деформируемого (2) и недеформируемого (1) коллекторов при одинако­ вых запасах газа в залежи

П р и м е ч а н

и е . Если согласно § 2

гл. XII произвести

экстраполяцию начально­

го участка зависимости p/z (р) = f [<2доб (01

до оси абсцисс (линия з), то оцениваем завы­ шенные начальные запасы газа в пласте

Уравнение (21) следует из принятия экспоненциального характера измене­ ния коэффициента пористости при снижении давления [53], т. е. когда

m= m0exp [ — ат (Рн— р (*))]•

Здесь т0 — коэффициент пористости при давлении рн.

Для примера рассмотрим пласт с начальными запасами газа Qзап = 100 млрд. м3 при рн = 300 кгс/см2. Содержание метана в газе 98%, пластовая температура

равняется 50я С.

Согласно [53], для гранулярных

коллекторов принято ат =

=

10-4 1/кгс/см2,

а для трещиноватых ат = Ю_3

1/(кгс/см2).

приведены

 

Результаты расчетов для условий трещиноватого коллектора

на

рис. 9. Вследствие деформации пласта зависимость — Кг— =

/ (<?доб (0)

2(р)

(линия 2) располагается выше соответствующей зависимости при отсутствии деформации (линия 1), что объясняется уменьшением во времени норового объема

залежи. При р = 0 линии 1 ж 2 сходятся в одну точку, так как независимо от того, является ли пласт деформируемым или нет, добытое количество газа к мо­

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

34


количества газа для месторождений Чирен (БНР) и Битковского аналогичны зависимости 2 на рис. 9.

Если по начальным точкам линии 2 попытаться определить начальные запасы газа в пласте (линия 3), то они окажутся завышенными (для рассма­ триваемого трещиноватого пласта) на 45%. Для гранулярного коллектора это завышение составляет около 5% .

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

Примем для упрощения, что растворимость газа в остаточной нефти под­ чиняется закону Генри1. Тогда уравнение материального баланса с учетом выделения из нефти газа и при пренебрежении усадкой нефти представится в виде (18). Суммарное количество газа, выделившегося из нефти к моменту

времени t,

приведенное

к рат и Гпл, будет

 

 

 

(?д (О = аннПну [Гн—Р(01-

(22)

Здесь

сснн — среднее

для залежи значение коэффициента нефтенасыщен-

ности; у — коэффициент

растворимости

в м3/м3 •(кгс/см2).

 

В примере расчета

было принято

о Нц = 0,16; ри = 420

кгс/см2; у =

= 0,438 м3/м® ■(кгс/см2).

по формулам (18) и (22) показывают,

что вследствие

Результаты расчетов

дегазации нефти зависимость p/z (р) = } (<?д0б (0) располагается выше соот­ ветствующей зависимостп при отсутствии остаточной нефтенасыщенности. В рассматриваемом примере дегазация увеличивает начальные запасы газа

впласте на 9,3%.

§5. Дифференциальные уравнения неустановившейся фильтрации газа в пористой среде

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

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

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

Лабораторные и промысловые исследования показывают зависимость коэф­ фициентов проницаемости и пористости от изменения горного (пластового)

1 При наличии соответствующих экспериментальных данных коэффициент растворимости в уравнении (22) принимается зависящим от давления, т. е.

у = у (р) [13].

3*

35