Научная статья на тему 'Численное моделирование закачки углекислого газа в истощенное месторождение углеводородов'

Численное моделирование закачки углекислого газа в истощенное месторождение углеводородов Текст научной статьи по специальности «Науки о Земле и смежные экологические науки»

CC BY
499
105
i Надоели баннеры? Вы всегда можете отключить рекламу.
Ключевые слова
ПОРИСТАЯ СРЕДА / ГАЗОГИДРАТЫ / УГЛЕКИСЛЫЙ ГАЗ / ФИЛЬТРАЦИЯ / ГИДРАТООБРАЗОВАНИЕ / POROUS MEDIUM / GAS HYDRATES / CARBON DIOXIDE / FILTRATION / HYDRATE FORMATION

Аннотация научной статьи по наукам о Земле и смежным экологическим наукам, автор научной работы — Хасанов Марат Камилович

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

i Надоели баннеры? Вы всегда можете отключить рекламу.

Похожие темы научных работ по наукам о Земле и смежным экологическим наукам , автор научной работы — Хасанов Марат Камилович

iНе можете найти то, что вам нужно? Попробуйте сервис подбора литературы.
i Надоели баннеры? Вы всегда можете отключить рекламу.

Now underground utilization of the carbon dioxide worked out by industrial facilities is one of directions of solving the problem of prevention of further increase in greenhouse gas concentration in the Earth’s atmosphere. At the same time the safest and the most reliable method of prevention of greenhouse gas exit to the surface is their conversion into the gas hydrate state, which allows storing the huge amount of gas in small volume and at low pressures. The aim of the work is the theoretical research of hydrodynamic and thermophysical processes at underground utilization of greenhouse gases in the gas hydrate state. Based on the methods of mechanics of continuous medium the author has built the mathematical model of carbon dioxide injection in a natural layer, saturated with water and methane in the initial state, followed by hydrate formation. For the case of injection of warm carbon dioxide (with the temperature exceeding the initial temperature of the layer) the author constructed the self-similar solutions of the plane-rational task describing dynamics of temperature and hydrodynamic fields in the layer. As a result of numerical calculations it was determined that formation of carbon dioxide gas hydrate can happen both on the frontal surface, and in the extended zone. The author defined the limiting values of the mass flow rate of carbon dioxide injection, corresponding to occurrence of an extended gas hydrate formation zone. It is shown that the mode with formation of gas hydrate in extended area is implemented at small values of mass flow rate of injection, as well as at high values of permeability and initial water saturation of layer. Implementation of the mode with extended area of hydrate formation in layers with high values of permeability and initial water saturation, and also at small values of the mass flow rate of injection is caused by the fact that under existing conditions the pressure on phase transitions border is below the equilibrium pressure of formation of carbon dioxide gas hydrate corresponding to the temperature on this border.

Текст научной работы на тему «Численное моделирование закачки углекислого газа в истощенное месторождение углеводородов»

УДК 532.546:536.421

ЧИСЛЕННОЕ МОДЕЛИРОВАНИЕ ЗАКАЧКИ УГЛЕКИСЛОГО ГАЗА В ИСТОЩЕННОЕ МЕСТОРОЖДЕНИЕ УГЛЕВОДОРОДОВ

Хасанов Марат Камилович,

hasanovmk@mail.ru

Стерлитамакский филиал Башкирского государственного университета, Россия, 453103, г. Стерлитамак, пр. Ленина, 37.

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

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

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

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

Ключевые слова:

Пористая среда, газогидраты, углекислый газ, фильтрация, гидратообразование.

Введение

В настоящее время в рамках борьбы с глобальным потеплением изучается потенциальная возможность подземного захоронения выработанного промышленными объектами углекислого газа [1-7]. Поскольку при долгосрочном хранении углекислого газа в виде флюида существует опасность его эмиссии в атмосферу, то в настоящее время изучается возможность подземного захоронения углекислого газа в газогидратном состоянии [8, 9]. Поскольку при одинаковых условиях в единице объема газогидрата содержится значительно больше газа, чем в свободном состоянии, то газогидратное состояние позволяет хранить большое количество углекислого газа при относительно низких давлениях [10, 11]. Это позволяет уменьшить риск выхода углекислого газа на поверхность и фактически увеличить емкость подземных резервуаров для хранения диоксида углерода. Кроме того, в природных пористых пластах площадь контакта воды и газа принимает огромные значения. Это свойство пористых сред благоприятно для интенсификации процесса образования гидрата и позволяет относительно быстро и экономично осуществить перевод углекислого газа в газогидратное состояние, например, путем его закачки в отработанные нефтяные или газовые коллекторы, частично насыщенные водой. Поскольку соответствующие технологические идеи должны быть подкреплены расчетами и

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

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

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

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

Равновесные значения температуры Т и давления р образования газогидрата достаточно хорошо описываются уравнением [11]:

Т = Т0 + Т*1п( р/р,0), (1)

где р80 - соответствующее температуре Т0 равновесное давление, Па; Т* - эмпирический параметр, К.

С учетом (1) начальное термобарическое состояние пласта удовлетворяет условиям:

г=0: Т=Т0, р=р0, р0>р80, (г>0,*=0).

В данной работе будем рассматривать случай, когда процесс гидратообразования лимитируется не кинетическими механизмами, а скоростью фильтрационного массопереноса в пористой среде. Это справедливо для масштабов времени, значительно превышающих характерное время кинетики процесса. Экспериментальные значения времени кинетики гидратообразования в пористых средах сильно зависят от характерных размеров пор, но в среднем составляют порядка несколько часов. Так, в работе [16], в которой исследовалось образование и разложение гидрата углекислого газа в пористой среде, данное время составило порядка одного часа. Эти характерные времена, как правило, очень малы по сравнению со временем распространения диоксида углерода в протяженном природном пласте.

В рассматриваемой задаче будем пренебрегать диффузионным перемешиванием углекислого газа и метана, т. к. вследствие непрерывной закачки диоксида углерода в пласт интенсивность массопе-реноса, обусловленного фильтрацией в проницаемой пористой среде, значительно превышает интенсивность массопереноса, обусловленного диффузией [13]. Кроме того, фронт вытеснения метана диоксидом углерода можно считать устойчивым вследствие ламинарного характера течений в пористых средах, а также большей по сравнению с метаном вязкостью углекислого газ [14]. Тогда, учитывая, что в рассматриваемой задаче начальное состояние пласта соответствуют условиям образования газогидрата диоксида углерода, можно полагать, что образование газогидрата диоксида углерода происходит на границе вытеснения метана углекислым газом. Следовательно, в рассматриваемом случае в пласте образуются две характерные области. В первой (ближней) области (/=1) поры насыщены углекислым газом и его газогидратом, а во второй (дальней) области (/=2) -метаном и водой. Соответственно газогидрат диоксида углерода образуется из углекислого газа и воды на подвижной межфазной поверхности, разделяющей первую и вторую области.

Основные уравнения

Система основных уравнений, описывающая процессы тепломассопереноса в пласте и представляющая собой законы сохранения масс и энергии, закон Дарси и уравнение состояния, в плоскорадиальном случае в каждой из областей имеет вид [17, 18]:

^ (р^ ) +1А ( грфДц1) = 0,

д1 г дг

I (рст) -.рс^дТ-=-дГ- (г^

01: дг г дг V дг

= -

к, др

И дг Р= Р/^Т.

Здесь I - время, с; г - радиальная координата, м; ф - пористость, д. ед; р -давление, Па; Т - температура, К; нижние индексы Ь=й,т относятся соответственно к параметрам углекислого газа и метана; р - плотность, кг/м3; к - фазовая проницаемость, м2; ц - действительная средняя скорость газа, м/с; С1 - удельная теплоемкость, Дж/(К-кг); И - динамическая вязкость газа, Па с; Я - газовая постоянная, Дж/(К-кг); Б1 - насыщенность пор пласта газом, д. ед; рС и Я - эффективные значения объемной теплоемкости, Дж/(К-м3) и коэффициента теплопроводности, Вт/(м-К) для пласта, насыщенного газом, водой или газогидратом. Поскольку основной вклад в значения рС и Я вносят соответствующие параметры горной породы, то в дальнейшем будем считать их постоянными величинами.

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

Зависимость коэффициента фазовой проницаемости к для углекислого газа и метана от их насыщенности Б1 и абсолютной проницаемости пласта к0 можно задать на основе формулы Козени [19]:

к, = кД (I = ё,т).

На границе между областями выполняются соотношения, следующие из условий баланса массы и теплоты:

ка дР1 =ф

Иа дг

= ф В^д а + о : у

V

ра

- К др2 = фог,

^ т т п?

Ит дг

ф$р„ (1 - а)Уп =ФД, ру ,

Я°Т-Яд-Т = фр*^уп.

дг дг

(3)

Здесь рк и р - плотность газогидрата и воды, кг/м3; Бк - насыщенность пор пласта газогидратом, д. ед.; Ьк - удельная теплота образования газогидрата углекислого газа, Дж/кг; нижние индекс /=1,2 относятся соответственно к параметрам в первой и второй области; Б1 - насыщенность пор водой, д. ед.; О - относительная массовая концентрация углекислого газа в газогидрате, д. ед.; Уп -скорость движения границы гидратообразования, м/с. Температуру и давление на этой границе будем считать непрерывными. Здесь и далее нижний индекс п относится к параметрам на границе между областями.

Из третьего уравнения системы (3) для величины гидратонасыщенности в первой области имеем:

рД1 О

Д =

р„ (1 - а)

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

массовым расходом Q и температурой Tw инжектируется углекислый газ. В результате его инжекции образуется область, насыщенная углекислым газом и его газогидратом. В данной работе будем рассматривать достаточно большие значения времени после начала инжекции углекислого газа, при которых радиус первой зоны гп значительно превышает радиус скважины тш. Тогда можно полагать, что величина тт практически не оказывает влияния на особенности протекания процесса гидрато-образования в пласте. Тогда условия на стенках скважины с учетом закона Дарси можно записать в виде:

г = к.: - ■

к,я ( др21

Т = Т (г ^ 0, г > 0),

где Q - массовый расход инжекции углекислого газа, кг/(м-с).

На основе системы (2) можно получить уравнения теплопроводности и пьезопроводности и в виде:

Р = х(рр)1 1

дг ' г дг\

Г дг У

дТ,

др2 дТ1

1 д( дТ 1

~^=Х(Т)X ^~^+х(Т)Л—I I •

дг дг дг г дг [ дг)

(4)

где

Iр) =-

каР1

Ус1Ф (1 - ^ )

I?) =

ктР2

Я,ф(1 - $ 0)

XI =

РаСа кс

Х = РтСт кт

X

%(Т) =_.

2Х У-тР! РС

2Х У/Р1 Начальные условия имеют вид:

*=0: Т=Т0, р=р0 (г>0,^=0).

Автомодельное решение

Введем автомодельную переменную:

Для этой переменной уравнения пьезопровод-ности и теплопроводности (4) примут вид:

ср2 х(р) 1 с ( с1р21

= 2\ , (1 = 1,2), с/ ) с/)

= 2хД ^ + 2 I.

(5)

-СТП=2Х Ср1 Т.+2С(/Т

1 с/ с/ / с/ г с/

Здесь уравнение пьезопроводности является нелинейным уравнением, т. к. в коэффициент пье-зопроводности х!р) входит неизвестная функция р. Рассматривая случай небольших перепадов давления в пласте, можно применить линеаризацию Лейбензона [20], т. е. переменное давление р] в параметре ц<р) принять постоянным и равным исходному пластовому давлению р0.

После интегрирования уравнений пьезопровод-ности и теплопроводности (5) решение для распре-

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

2 2 + / 1 А = рп +-т

як.

1-1 ( /

iНе можете найти то, что вам нужно? Попробуйте сервис подбора литературы.

Т = Тп +

Т - Тп) |/ехр - X, р2)с/

|/ехр [-£ - X! р?\с%

0 </</

(6)

(р - 4гхр[-^

г 1 ( /,

11ехр I-4)

(Тп - То)}/ехр- X2р22)

р2 = р0 +

Т = Т + -

2 10 т

}/ехр(-/-X2р2) с/

(7)

£п </<Г,

где ^=Х>У/Т>0=1,2).

На основе соотношений (3) с учетом решений для давления и температуры (6) и (7) получим уравнения для нахождения координаты границы ги-дратообразования /(п) и значений параметров р(п) и

Т(п) на ней:

QУсRЛ як.

ехр [-|)= Л р&

(р" -р2)ехр(-^1 Г 1 ( / 1 ^

!гхр I-Й)

(Тп - Т,)ехр-X!рг

= 4 рп/2

(8)

(9)

1 ( / 1 |_ехр [-^ - Х. р2 ]

где

Т - Т(п))ехр- X2рп2

г 1 ( 1 _ |_ехр^-Х2р2]

= (рж_+1 - ]

кс [ Р0с )

(10)

кс [ р0с

А = (1 - в = к

2рС

Записанная система уравнений в работе решалась следующим образом. Выразив из уравнения (8) величину рп (как функцию 4) и подставив данное значение рп в уравнение (9), получаем трансцендентное уравнение с одной неизвестной £п, которое в работе решено методом половинного деления. Затем из уравнений (8) и (10) определяем значения давления рп и температуры Тп на границе £п между областями.

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

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

На рис. 1 представлены распределения температуры и давления в пласте при двух разных значениях расхода закачки углекислого газа 0=0,03 кг/(м-с) (а) и 0,01 кг/(м-с) (Ь). Для осталь-

ных параметров приняты следующие значения: ф=0,2, $10=0,1, р0=2,9 МПа, Т=278 К, Т0=276 К, к)=10-14м2, 0=0,28, Я=2 Вт/(м-К), рС=2,5-106Дж/(К-кг), рй=1100 кг/м3, С,=890 Дж/(Ккг), Ст=1560 Дж/(Ккг), Я,=189 Дж/(К-кг), Ят=520 Дж/(Ккг), и=1,3-10-5Па-с, Ит=105Па-с, 4=4-105Дж/кг, Т*=7,6 К, р,0=1,68 МПа.

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

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

Рис. 1. Распределение температуры, давления и равновесного давления (штриховая кривая) при инжекции углекислого газа с расходом 0,03 кг/(м-с) (а) и 0,01 кг/(м-с) (b)

Fig. 1. Distribution of temperature, pressure and equilibrium pressure (dashed curve) when carbon dioxide is injected at a flow rate of 0,03 kg/(m-s) (a) and 0,01 kg/(m-s) (b)

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

T . >„ = T+teLA / PC = To + .fGcV

(1 - G) pC

Минимальное значение для величины равновесной температуры образования газогидрата углекислого газа определяется выражением: Ts = T0 + T* ln( pj ps 0),

где Ts - температура образования газогидрата углекислого газа, вычисляемая по формуле (1) и соответствующая давлению p0.

Тогда необходимое условие реализации режима с протяженной областью гидратообразования имеет вид:

T*ln( Po / Pso)

(1 - G) pC

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

Необходимое и достаточное условие реализации режима с протяженной областью гидратообра-зования определяется неравенством: T < T

где Tn - температура на границе гидратообразования, определяемая из уравнения (10).

Рис. 2. Зависимость предельного значения расхода от проницаемости при $о=0,10 (кривая 1) и S0=0,11 (кривая 2)

Fig. 2. Limit mass flow rate versus the permeability for S=0,10 (curve 1) and Sl0=0,11 (curve 2)

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

На рис. 2 приведена зависимость предельного значения расхода закачиваемого углекислого газа от проницаемости пласта при разных значениях начальной водонасыщенности пласта £10=0,10 (кривая 1) и £10=0,11 (кривая 2). Значения остальных параметров соответствуют рис. 1. Из рис. 2 видно, что с повышением проницаемости пласта предельное значение расхода увеличивается, причем тем быстрее, чем выше начальная водонасы-щенность пласта. Таким образом, режим с протяженной областью гиратообразования реализуется в высокопроницаемых пластах, а также в случае высоких значений исходной водонасыщенности пласта. Это обусловлено тем, что с одной стороны увеличение проницаемости при неизменном массовом расходе требует уменьшения давления на скважине, что приводит к снижению давления в пласте. С другой стороны увеличение исходной во-донасыщенности приводит к росту тепловыделения на границе фазовых переходов и соответственно к повышению температуры пласта и равновесного давления гидратообразования. Таким образом, реализация режима с протяженной областью гидратообразования в случаях высоких значений проницаемости пласта и его исходной водонасы-щенности обусловлена тем, что при данных условиях давление на границе фазовых переходов оказывается ниже равновесного давления образования газогидрата углекислого газа, соответствующего температуре на этой границе.

Выводы

Из анализа результатов численных расчетов следует, что в случае инжекции углекислого газа в пласт, насыщенный в исходном состоянии водой и метаном и имеющий проницаемость в диапазоне &0=10~14-10~12м2, в зависимости от расхода закачки и исходной водонасыщенности пласта образование газогидрата диоксида углерода может происходить как на фронтальной поверхности, так и в протяженной области. Для величины расхода закачки углекислого газа существует предельное значение, равное примерно 1 кг/(м-с) при проницаемости £0=1013м2и водонасыщенности Я10=0,11, ниже которого образование газогидрата происходит в протяженной области. Величина данного предельного значения расхода закачки углекислого газа увеличивается примерно прямо пропорционально проницаемости пласта.

СПИСОК ЛИТЕРАТУРЫ

1. Bradshaw J., Cook P. Geological sequestration of carbon dioxide // Environmental Geosciences. - 2001. - V. 8. - № 3. - P. 149-151.

2. Karacan C.O., Mitchell G.D. Behavior and effect of different coal microlithotymes during gas transport for carbon dioxide sequestration into coal seams // Int. J. Coal Geology. - 2003. - V. 53. -P. 201-217.

3. Oldenburg C.M., Pruess R., Benson S.M. Process modeling of CO2 injection into natural gas reservoirs for carbon sequestration and enhanced gas recovery // Energy Fuel. - 2001. - V. 15. -P. 293-298.

4. Benson S.M, Cole D.R. CO2 sequestration in deep sedimentary formations // Elements. - 2008. - V. 4. - P. 325-331.

5. Akinnikawe O., Ehlig-Economides C.A. Geologic model and fluid flow simulation of Woodbine aquifer CO2 sequestration // International Journal of Greenhouse Gas Control. - 2016. - V. 49. -P. 1-13.

6. Lombardi S., Altunina L.K., Beaubien S.E. Advances in the geological storage of carbon dioxide. - Berlin: Springer Publ., 2006. - 362 p.

7. IPCC special report on carbon dioxide capture and storage / B. Metz, O. Davidson, H.C. Coninck, M. Loos, L.A. Meyer. -Cambridge: Cambridge University Press, 2005. - 442 p.

8. Чувилин Е.М., Гурьева О.М. Экспериментальное изучение образования гидратов СО2 в поровом пространстве промерзающих и мерзлых пород // Криосфера Земли. - 2009. - Т. 13. -№ 3. - С. 70-79.

9. Оценка возможности захоронения углекислого газа в криоли-тозоне Западной Сибири / А.Д. Дучков, Л.С. Соколова, Д.Е. Аюнов, М.Е. Пермяков // Криосфера Земли. - 2009. -Т. 13. - №4. - С. 62-68.

10. Sloan E.D., Koh C.A. Clathrate Hydrates of Natural Gases. - Boca Raton: CRC Press, 2007. - 752 p.

11. Бык С.Ш., Макогон Ю.Ф., Фомина В.И. Газовые гидраты. -М.: Химия, 1980. - 296 с.

12. Хасанов М.К. Исследование режимов образования газогидратов в пористой среде, частично насыщенной льдом // Теплофизика и аэромеханика. - 2015. - Т. 22. - № 2. - С. 255-266.

13. Theoretical research of the gas hydrate deposits development using the injection of carbon dioxide / V.Sh. Shagapov, M.K. Khas-anov, N.G. Musakaev, Ngoc Hai Duong // International Journal of Heat and Mass Transfer. - 2017. - V. 107. - P. 347-357.

14. Цыпкин Г.Г. Образование гидрата углекислого газа при его инжекции в истощенное месторождение углеводородов // Изв. РАН. Механика жидкости и газа. - 2014. - № 6. - С. 101-108.

15. Цыпкин Г.Г. Образование гидрата при инжекции жидкой двуокиси углерода в пласт, насыщенный метаном и водой // Изв. РАН. Механика жидкости и газа. - 2016. - № 5. -С. 99-107.

16. MRI measurements of CO2 hydrate dissociation rate in a porous medium / M. Yang, Y. Song, Y. Zhao, Y. Liu, L. Jiang, Q. Li // Magn. Reson. Imaging. - 2011. - V. 29. - P. 1007-1013.

17. Хасанов М.К., Мусакаев Н.Г., Гималтдинов И.К. Особенности разложения газогидратов с образованием льда в пористой среде // Инженерно-физический журнал. - 2015. - Т. 88. - № 5. -С. 1022-1031.

18. Гималтдинов И.К., Хасанов М.К. Математическая модель образования газогидрата при инжекции газа в пласт, частично насыщенный льдом // Прикладная математика и механика. -2016. - Т. 80. - № 1. - С. 80-90.

19. Баренблатт Г.И., Ентов В.М., Рыжик В.М. Движение жидкостей и газов в природных пластах. - М.: Недра, 1982. - 208 c.

20. Басниев К.С., Кочина И.Н., Максимов В.М. Подземная гидромеханика. - М.: Недра, 1993. - 416 c.

Поступила 15.06.2017 г.

Информация об авторах

Хасанов М.К., кандидат физико-математических наук, доцент кафедры прикладной информатики и программирования Факультета математики и информационных технологий Стерлитамакского филиала Башкирского государственного университета.

UDC 532.546:536.421

NUMERICAL SIMULATION OF CARBON DIOXIDE INJECTION INTO A DEPLETED HYDROCARBON DEPOSIT

Marat K. Khasanov,

hasanovmk@mail.ru

Sterlitamak Branch of Bashkir State University, 37, Lenin Avenue, Sterlitamak, 453103, Russia.

Now underground utilization of the carbon dioxide worked out by industrial facilities is one of directions of solving the problem of prevention of further increase in greenhouse gas concentration in the Earth's atmosphere. At the same time the safest and the most reliable method of prevention of greenhouse gas exit to the surface is their conversion into the gas hydrate state, which allows storing the huge amount of gas in the small volume and at low pressures.

The aim of the work is the theoretical research of hydrodynamic and thermophysical processes at underground utilization of greenhouse gases in the gas hydrate state.

Based on the methods of mechanics of continuous medium the author has built the mathematical model of carbon dioxide injection in a natural layer, saturated with the water and methane in the initial state, followed by hydrate formation. For the case of injection of the warm carbon dioxide (with the temperature exceeding the initial temperature of the layer) the author constructed the self-similar solutions of the plane-rational task describing dynamics of temperature and hydrodynamic fields in the layer.

As a result of numerical calculations it was determined that formation of carbon dioxide gas hydrate can happen both on the frontal surface, and in the extended zone. The author defined the limiting values of the mass flow rate of carbon dioxide injection, corresponding to occurrence of an extended gas hydrate formation zone. It is shown that the mode with formation of gas hydrate in extended area is implemented at small values of mass flow rate of injection, as well as at high values of permeability and initial water saturation of layer. Implementation of the mode with extended area of hydrate formation in layers with high values of permeability and initial water saturation, and also at small values of the mass flow rate of injection is caused by the fact that under existing conditions the pressure on phase transitions border is below the equilibrium pressure of formation of carbon dioxide gas hydrate corresponding to the temperature on this border.

Key words:

Porous medium, gas hydrates, carbon dioxide, filtration, hydrate formation.

REFERNCES

1. Bradshaw J., Cook P. Geological sequestration of carbon dioxide. Environmental Geosciences, 2001, vol. 8, no. 3, pp. 149-151.

2. Karacan C.O., Mitchell G.D. Behavior and effect of different coal microlithotymes during gas transport for carbon dioxide sequestration into coal seams. Int. J. Coal Geology, 2003, vol. 53, pp. 201-217.

iНе можете найти то, что вам нужно? Попробуйте сервис подбора литературы.

3. Oldenburg C.M., Pruess R., Benson S.M. Process modeling of CO2 injection into natural gas reservoirs for carbon sequestration and enhanced gas recovery. Energy Fuel, 2001, vol. 15, pp. 293-298.

4. Benson S.M, Cole D.R. CO2 sequestration in deep sedimentary formations. Elements, 2008, vol. 4, pp. 325-331.

5. Akinnikawe O., Ehlig-Economides C.A. Geologic model and fluid flow simulation of Woodbine aquifer CO2 sequestration. International Journal of Greenhouse Gas Control, 2016, vol. 49, pp. 1-13.

6. Lombardi S., Altunina L.K., Beaubien S.E. Advances in the geological storage of carbon dioxide. Berlin, Springer Publ., 2006. 362 p.

7. Metz B., Davidson O., Coninck H.C., Loos M., Meyer L.A. IPCC special report on carbon dioxide capture and storage. Cambridge, Cambridge University Press, 2005. 442 p.

8. Chuvilin E.M., Guryeva O.M. Experimental investigation of CO2 gas hydrate formation in porous media of frozen and freezing sediments. Earth's Cryosphere, 2009, vol. 13, no. 3, pp. 70-79. In Rus.

9. Duchkov A.D., Sokolova L.S., Ayunov D.E., Permyakov M.E. Assessment of potential of west Siberian permafrost for the carbon dioxide storage. Earth's Cryosphere, 2009, vol. 13, no. 4, pp. 62-68. In Rus.

10. Sloan E.D., Koh C.A. Clathrate Hydrates of Natural Gases. Boca Raton, CRC Press, 2007. 752 p.

11. Byk S.Sh., Makogon Yu.F., Fomina V.I. Gazovye gidraty [Gas hydrates]. Moscow, Khimiya Publ., 1980. 296 p.

12. Khasanov M.K. Investigation of regimes of gas hydrate formati-

on in a porous medium, partially saturated with ice. Thermophys-ics and Aeromechanics, 2015, vol. 22, no. 2, pp. 255-265. In Rus.

13. Shagapov V.Sh., Khasanov M.K., Musakaev N.G., Ngoc Hai Duong Theoretical research of the gas hydrate deposits development using the injection of carbon dioxide. International Journal of Heat and Mass Transfer, 2017, vol. 107, pp. 347-357.

14. Tsypkin G.G. Formation of carbon dioxide hydrate at the injection of carbon dioxide into a depleted hydrocarbon field. Fluid Dynamics, 2014, vol. 49, no. 6, pp. 789-795. In Rus.

15. Tsypkin G.G. Formation of hydrate in injection of liquid carbon dioxide into a reservoir saturated with methane and water. Fluid Dynamics, 2016, vol. 51, no. 5, pp. 672-679. In Rus.

16. Yang M., Song Y., Zhao Y., Liu Y., Jiang L., Li Q. MRI measurements of CO2 hydrate dissociation rate in a porous medium. Magn. Reson. Imaging, 2011, vol. 29, pp. 1007-1013.

17. Khasanov M.K., Musakaev N.G., Gimaltdinov I.K. Features of the Decomposition of Gas Hydrates with the Formation of Ice in a Porous Medium. Journal of Engineering Physics and Ther-mophysics, 2015, vol. 88, no. 5, pp. 1052-1061. In Rus.

18. Gimaltdinov I.K., Khasanov M.K. Mathematical model of the formation of a gas hydrate on the injection of gas into a stratum partially saturated with ice. Journal of Applied Mathematics and Mechanics, 2016, vol. 80, no. 1, pp. 57-64. In Rus.

19. Barenblatt G.I., Entov V.M., Ryzhik V.M. Dvizhenie zhidkostey i gazov v prirodnykh plastakh [Movement of liquids and gases in natural formations]. Moscow, Nedra Publ., 1982. 208 p.

20. Basniev K.S, Kochina I.N., Maksimov V.M. Podzemnaya gidro-mekhanika [Underground hydromechanics]. Moscow, Nedra Publ., 1993. 416 p.

Received: 15 June 2017.

Information about the author

Marat K. Khasanov, Cand. Sc., associate professor, Sterlitamak Branch of Bashkir State University.

i Надоели баннеры? Вы всегда можете отключить рекламу.