Файл: Закиров, С. Н. Проектирование и разработка газовых месторождений учебное пособие.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(0— 108 |
( 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 |