Моделирование процессов образования и разложения газовых гидратов в пористой среде при депрессионном воздействии

Вестник МГТУ, том 16, № 4, 2013 г.

стр.803-807

УДК 532.546

Моделирование процессов образования и разложения газовых гидратов в пористой среде при депрессионном воздействии

М.Х. Хайруллин, П.Е. Морозов, А.И. Абдуллин, М.Н. Шамсиев

Институт механики и машиностроения КазНЦ РАН

Аннотация. На основе численного моделирования двухфазной неизотермической фильтрации в пористой среде, насыщенной газом, водой и гидратом, проведён анализ депрессионного метода воздействия на газогидратную залежь.

Abstract. In this paper the analysis of depressurization in gas hydrate reservoir has been performed. It’s based on numerical simulations of nonisothermal two-phase filtration in porous medium saturated with gas, water and hydrate.

Ключевые слова: газовые гидраты, пористая среда, диссоциация, моделирование, депрессионное воздействие Key words: gas hydrates, porous medium, dissociation, simulation, depressurization

1. Введение

Г азовые гидраты представляют собой твёрдые кристаллические соединения молекул газа и воды, устойчивые при определённых давлениях и температурах. При диссоциации одного кубического метра газового гидрата выделяется до 180 м3 газа и 0,8 м3 воды (Makogon, 2010). Природные газовые гидраты, состоящие в основном из гидрата метана, рассматриваются как перспективные источники углеводородного сырья. К настоящему времени газогидраты обнаружены в районах, расположенных вдоль побережий Северной и Южной Америки, Евроазиатского континента, в Австралии, Индии, Японии, под Чёрным, Каспийским и Средиземным морями, озером Байкал и др. По различным оценкам ресурсы природных газовых гидратов значительно превышают суммарные ресурсы всех горючих ископаемых. Поэтому актуальным является вопрос о поиске экономически рентабельных технологий разработки залежей газовых гидратов (Makogon, 2010; Birchwood et al., 2010). В настоящее время рассматриваются три основных метода добычи газа из гидратного пласта: понижение давления ниже равновесного, нагрев гидратосодержащих пород выше температуры разложения гидрата и их комбинация

(рис. 1).

Рис. 1. Фазовая диаграмма системы "газ - вода - гидрат" и возможные методы добычи газа из гидратов. 1 - депрессионное воздействие,

2 - тепловое воздействие, 3 - комбинированное воздействие

803

Хайруллин М.Х. и др. Моделирование процессов образования и разложения...

2. Постановка задачи

Для математического моделирования процесса объёмной диссоциации газовых гидратов в пористой среде было предложено две модели - равновесная и кинетическая. В основе первой модели лежит предположение о том, что в каждом элементарном объёме пористой среды выполняется условие термодинамического равновесия системы "газ - вода - гидрат" (Бондарев и др., 1976; Максимов, Цыпкин, 1990; Максимов, 1992; Бондарев, Попов, 2002; Нигматуллин и др., 1998). Кинетическая модель основана на экспериментальной зависимости скорости диссоциации газогидрата от давления, температуры и размеров его частиц (Kim et al, 1987). В работе (Kowalsky, Moridis, 2007) показано, что при численном моделировании процесса диссоциации газовых гидратов в пласте равновесная и кинетическая модели дают близкие результаты. Кроме того, отмечается, что использование равновесной модели предпочтительнее, так как она требует меньше вычислительных ресурсов.

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

Основные уравнения, описывающие динамику равновесной смеси газа, воды и гидрата в пористой среде, выводятся из законов сохранения массы и энергии (Максимов, 1992; Бондарев, Попов, 2002; Нигматуллин и др., 1998; Цыпкин, 2009). В радиальных координатах система уравнений имеет следующий вид (Бондарев, Попов, 2002):

p1 f -mqp it -m(1 - Sw»-S‘)rc-his -rffrXiT

+

p c — + p c fw

r g g r w w

Mg Mw

+ k

m—f (1 - S )(1 - Sh )-0 = — it { wA h zRT) r fr

m i (Sw (1 - Sh)) = -1 i-( rXM fp) - m(1 -e)P

ft r fr { m dr ) p

dp fT = 0 .

fr fr

k (Sh )f p fp M zRT fr

Л

~mPh£-

fSlL

ft

fSh

ft

{ Mw fr) Pw

где p - давление; T - температура; m - пористость; pg, pw - вязкость газа и воды; Sw - водонасыщенность; Sh - гидратонасыщенность; fg, fw - функции относительных фазовых проницаемостей для газа и воды; pg, pw, ph - плотность газа, воды и гидрата; е - массовое содержание газа в гидрате; z - коэффициент сверхсжимаемости газа; (pC)e - объёмная теплоёмкость скелета и вмещающих флюидов; cg, cw -теплоёмкость газа и воды; q - теплота фазового перехода гидрата; п - коэффициент адиабатического расширения; KJT - коэффициент Джоуля-Томсона; X - коэффициент теплопроводности скелета и вмещающих флюидов.

Условие термодинамического равновесия смеси "газ - вода - гидрат" описывается эмпирическим уравнением (Бондарев и др., 1976):

T = alnp + b,

где a и b константы, определяемые экспериментально.

Предполагается, что в начальный момент времени пласт имеет гидратонасыщенность Sh0, водонасыщенность Sw0, равновесное давление p0 и температуру T0 = alnp0 + b. На контуре питания пласта поддерживается постоянное давление p0 и температура T0. На забое скважины задаётся постоянное

p

Из решения системы нелинейных уравнений (1-4) находятся неизвестные значения функций p, T, Sh и Sw. Массовый дебит газа и воды в пластовых условиях определяется в виде:

Q, = 2Pk(Sh )pg Mr 'fp ;

m frr=,

Qw = 2pHk(Sh)pw ^r 'fp ,

Mw frr-,

804

Вестник МГТУ, том 16, № 4, 2013 г.

стр.803-807

где rc - радиус скважины; H - толщина пласта.

Для численного решения данной задачи используется метод конечных разностей со сгущающейся в окрестности скважины сеткой узлов. Построение такой сетки осуществляется с помощью преобразования координат u = lnr (Азиз, Сеттари, 1982). Конечно-разностные аналоги уравнений (1-3) имеют вид:

(РС).

Tn - ГГ1 и „.V. ^ р" - рТ

Dt

1--Mhcg (1 - )(1 - S" )p^‘ -

Tn - 2Tn + Tn f" / —г

-Ae-2u ------i---^ + k№.)p".c ^-KITe-2u

h2 g m it h 2

Dt

n

(p" - p- )2

-k (sr jpc, m+p„cm m

= mqph

?-2u

pn - pn, Tn - T."

Г i fi-1 i i-1

h

h

rt n rt n-1 _____

hi - hi , i = 1, M ; Dt.

m

Dt

f p" p"-1 ^

(1 - S" )(1 - S" ^— (1 - S"-1 )(1 - S"-1 )^—

\ wi J\ hi / zRt T" ' wi ' hi ' zRt T "-1

--mPh£

rt n rt n-1 ________

hi - hi i = 1, M ;

Dt.

Dt.

n n-1

(s: (1 - S"M)- S"-1 (1 - Sr))+ m(1 - e)^^—hi

Pw Dtn

m

k (S„ )fw

Mw

i+12

k (S„ )fw

Mw

pn ~*n

Pi-1

-12

n

2u

2u

= e

- e

2

h

h

Система нелинейных уравнений (6-8) решается итерационно. Вычислительный алгоритм заключается в следующем:

1. Из системы уравнений (7) методом потоковой прогонки находится распределение давления р.

2. Из решения системы уравнений (8) определяется водонасыщенность Sw.

3. Если S"1 > 0, то с учётом условия термодинамического равновесия (4) из уравнений (6) определяется гидратонасыщенность на новом временном слое. В случае S"1 = 0 из уравнений (6) находится распределение температуры. 3

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

Расчёты проводились при следующих значениях параметров: р0 = 10 МПа; pc = 5 МПа; Sh0 = 0,39; Sw0 = 0,21; rc = 0,1 м; Rk = 500 м; H = 10 м; pw = 1000 кг/м3; ph = 920 кг/м3; m = 0,2; k0 = 0,01 мкм2; N = 2;

= 1 мПа/сек; fg = 0,013 мПа/сек; s = 0,147; a = 7,28; b = 169,7 K; (pC)e = 1,48 Дж/(м3-К); X = 1,71 Вт/(м-К); cg = 2093 Дж/(кгК); cw = 4200 Дж/(кгК); ch = 3210 Дж/(кгК); KIT = 0,4 К/МПа; n = 0,14 К/МПа; q = 5,1 105 Дж/кг. Значения параметров для газа соответствуют чистому метану.

Коэффициент сверхсжимаемости газа вычисляется по формуле Латонова-Гуревича (Бондарев, Попов, 2002):

Г

z( P,T ) =

0,17376ln

f t ^

T

\ кр

+ 0,73

+ 0,1-

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

(Максимов, Цыпкин, 1990):

кР

805

Хайруллин М.Х. и др. Моделирование процессов образования и разложения...

fg

(1 - Sw /0,9)3'5 (1 + 3Sw), 0 < Sw < 0,9, 0, Sw > 0,9'

fw

((Sw - 0'2)/0'8)3'5'0,2 < Sw < 1,

0,0 < S < 0,2.

w

Зависимость проницаемости от гидратонасыщенности принимается в виде (Бондарев и др., 1976; Нигматуллин и др., 1998):

k(Sh) = *с(1 - Sh)N,

где к0 - абсолютная проницаемость породы в отсутствии гидрата.

На рис. 2 представлены графики распределения гидрато- и водонасыщенности в различные моменты времени. В результате разложения гидрата с течением времени водонасыщенность увеличивается по всему пласту, принимая максимальное значение за фронтом диссоциации.

Рис. 2. Распределение гидрато- и водонасыщенности Рис. 3. Массовый расход газа

на моменты времени: 1 - 1 сут, 2 - 10 сут, 3 - 100 сут

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

Общий объём разложившегося гидрата определяется по формуле:

Vh = 2лИш t г S0 - Sh

Расчёты показали, что в результате диссоциации газового гидрата в течение ста суток из пласта дополнительно добывается порядка 106 м3 газа.

В процессе добычи газа при определённых термобарических условиях возможно вторичное образование газовых гидратов в окрестности скважины (Бондарев, Попов, 2002). На рис. 4 представлены графики распределения гидрато- и водонасыщенности при начальном пластовом давлении р0 = 20 МПа и забойном давлении pc = 15 МПа. В результате охлаждения газа вблизи скважины гидратонасыщенность начинает расти до тех пор, пока вся свободная вода не перейдёт в гидрат. В удалённой зоне пласта, как и в предыдущем примере, происходит диссоциация газогидрата. С ростом гидратонасыщенности в окрестности скважины фильтрационные свойства призабойной зоны ухудшаются. Зависимость массового расхода газа от времени представлена на рис. 5.

4. Заключение

В работе показано, что при депрессионном воздействии объём разложившегося гидрата незначителен по сравнению с общим объёмом гидрата в пласте. Поэтому для разработки газогидратной

806

Вестник МГТУ, том 16, № 4, 2013 г.

стр.803-807

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

на моменты времени: 1 - 1 сут, 2 - 10 сут, 3 - 100 сут

Работа выполнена при финансовой поддержке РФФИ (грант № 11-05-00878-а).

Литература

Birchwood R., Dai J., Shelander D. Developments in gas hydrate. Oilfield review, Spring, p. 18-33, 2010.

Kim H.C., Bishnoi P.R., Heidemann R.A., Rizvi S.S.H. Kinetics of methane hydrate decomposition. Chem. Eng. Science, v. 42, N 7, p. 1645-1653, 1987.

Kowalsky M.B., Moridis G.J. Comparison of kinetic and equilibrium reaction models in simulating gas hydrate behaviour in porous media. Energy Conversion and Management, N 48, p. 1850-1863, 2007.

Makogon Y.F. Natural gas hydrates - a promising source of energy. J. Natural Gas Science and Eng., N 2, p. 49-59, 2010.

Азиз Х., Сеттари Э. Математическое моделирование пластовых систем. М., Недра, 407 с., 1982.

Бондарев Э.А., Бабэ Г.Д., Гройсман А.Г., Каниболотский М.А. Механика образования гидратов в газовых потоках. Новосибирск, Наука, 157 с., 1976.

Бондарев Э.А., Попов В.В. Динамика образования гидратов при добыче природного газа. Выч. технологии, т. 7, № 5, с. 28-33, 2002.

Максимов А.М. Математическая модель объёмной диссоциации газовых гидратов в пористой среде: учёт подвижности водной фазы. ИФЖ, т. 62, № 1, с. 76-81, 1992.

Максимов А.М., Цыпкин Г.Г. О разложении газовых гидратов, сосуществующих с газом в природных пластах. Изв. АН СССР. МЖГ, т. 25, № 5, с. 84-87, 1990.

Нигматуллин Р.И., Шагапов В.Ш., Сыртланов В.Р. Автомодельная задача о разложении газогидратов в пористой среде при депрессии и нагреве. ПМТФ, т. 39, № 3, с. 111-118, 1998.

Цыпкин Г.Г. Течения с фазовыми переходами в пористых средах. М., Физматлит, 232 с., 2009.

807