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

290 Вестник СамГУ — Естественнонаучная серия. 2008. №3(62).

УДК 532.546; 536.421

ИНЖЕКЦИЯ ГАЗА В ПОРИСТУЮ СРЕДУ, СОПРОВОЖДАЮЩАЯСЯ ОБРАЗОВАНИЕМ ГАЗОГИДРАТА

© 2008 М.К. Хасанов1

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

Ключевые слова: пористая среда, образование газогидрата.

Введение

Образование газовых гидратов в пористых структурах в настоящий момент имеет широкие промышленные перспективы, связанные, в первую очередь, с возможностью хранения газа в гидратном состоянии [1]. Рядом исследователей предлагается подземная газогидратная консервация парниковых газов, которая обеспечивает высокий уровень безопасности хранения и не требует больших энергетических затрат [2]. В основе гидратного способа хранения газа положено то обстоятельство, что при одинаковых условиях в единице объема в гидратном состоянии содержится значительно больше газа, чем в свободном состоянии [3].

Вопросам математического моделирования процессов разложения газовых гидратов при тепловом и депрессионном воздействии на пористую среду, насыщенную газовым гидратом и газом, посвящена работа [4]. В работе [5] исследованы некоторые закономерности процесса образования газо-гидрата в пористой среде при закачке газа в пористую среду, насыщенную газом и водой. В настоящей работе проанализированы условия существования двух качественно различных режимов образования гидрата в пористой среде при инжекции газа.

1 Хасанов Марат Камилович (hasanovmk@mail.ru), кафедра прикладной математики и механики Стерлитамакской государственной педагогической академии, 453103, Россия, Башкортостан, г. Стерлитамак, пр. Ленина, 37.

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

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

Здесь С —массовое содержание газа в гидрате; т —пористость; р — давление, Т —температура; ру и S у — истинная плотность и насыщенность пор г]-й фазы (] = Н, I, д); индексы д, I, Н относятся к параметрам газа, воды и гидрата соответственно; , к§ и ^ —скорость, проницаемость и динамическая вязкость для газовой фазы; Ьь — удельная теплота гидратообра-зования; рс и X — удельная объемная теплоемкость и коэффициент теплопроводности системы. Анализ показывает, что для случая, когда исходная водонасыщенность пор менее половины, можно пренебречь переменностью удельной объемной теплоемкости и коэффициента теплопроводности системы.

Зависимость коэффициента проницаемости для газа к§ от газонасыщен-ности можно задать на основе формулы Козени:

дт_д_/1дт

§ § дх дх \ дх

Р = PgRgT■

где ко — абсолютная проницаемость пласта.

Значения температуры и давления в области образования гидрата связаны c условием фазового равновесия [1]:

Г, = Г0 + Т, 1п

Р*о

где То — исходная температура системы, р^о — равновесное давление, соответствующее исходной температуре, Т* — эмпирический параметр, зависящий от вида газогидрата.

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

[т (5йрй(1 - О) + 5грО Х«] = 0,

[т (р^Си^ - Хм) - ркБ нО Хм)] = 0,

. дТ

дх

= \тр]гТ]г5 ■

Здесь [у] —скачок параметра у на границе Х(^ между областями; Х^) —скорость движения этой границы. Температуру и давление на этих границах будем полагать непрерывными.

Будем полагать, что пласт в начальный момент времени насыщен газом и водой, давление ро и температура То которых в исходном состоянии соответствуют термодинамическим условиям существования их в свободном состоянии (ро ^ Рsо) и изначально одинаковы во всем пласте. Эти условия могут быть записаны следующим образом:

г = 0 : Т = То, р = ро (х ^ 0) ■

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

Х = о : Т = Те, р = ре (г > 0) ■

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

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

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

Сформулированная задача имеет автомодельное решение. Введем автомодельную переменную:

V

? =

где к(Т) — коэффициент температуропроводности пласта. Тогда решение задачи может быть записано в виде:

1«

(Р2е-Рі)) I' еХР(--^— Кь

2 2 Р(1) = РМ +

1

-2

4П(1)

0 < 1 < ?(*),

I

о

!«

(Ге - Г(*)^ еХР 1

ехр | — -------| сЩ

Г(1) = Гм +

4П(1)

?2 Ре(1) 7

4 2рі

1«

/

0 < 1 < ?(*)>

еХР

1 Ре(1) 2

4 2А Р<>'

01

(РМ-

то

Р20>/

22 Р (2) = Р0 +

1

ехр|"і|;и

1

ехр | — ------| аь,

1« < 1 < то ,

то

(Г« - Го) J

еХР

4П(2)

12 Ре(2) 2

4 2Й

01

Г(2) = Г0 +

1

1«

еХР

12 Ре(2) 2

1« < 1 < то ,

4 2^'

01

к,:

(р)

где П(<-) =

©

к;:/ =

р£0С£ &(і) р0

з

%) = .

к(Г^ (і) " 'ч‘) Л ^

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

0

1

Рис. 1. Распределение температуры пласта (1), равновесной температуры (2) и давления. а — непротиворечивый случай, b — случай переохлаждения за фронтом.

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

На рис. 1, а представлены распределения температуры и давления при нагнетании метана под давлением pe = 5 МПа и с температурой Te = 278 К в пласт с начальным давлением ро = 4 МПа и исходной водонасыщенно-стью Sio = 0,2. Для остальных параметров, характеризующих систему, приняты следующие значения: m = 0,1, G = 0,12, To = 280 К, T* = 10 К, ps0 = 5,5 МПа, k0 = 10-14 м2, Lh = 5 ■ 105 Дж/кг, cg = 1560 Дж/(К-кг). Пунктирная линия 2 показывает равновесную температуру, соответствующую полученному распределению давления. Из данного рисунка видно, что температура пласта перед фронтом гидратообразования ниже равновесной температуры, а за фронтом — выше этой температуры, что является физически непротиворечивым. Следовательно, в этом случае можно говорить о существовании решения с фронтальной поверхностью образования гидрата.

На рис. 1, b приведены распределения температуры и давления при нагнетании газа под давлением pe = 6 МПа и с температурой Te = 278 К. Видно, что температура пласта за фронтом гидратообразования опускается ниже равновесной температуры (пунктирная линия 2), т.е. в этой области наблюдается переохлаждение воды. Таким образом, при данном давлении

р„МПа

ь

р., МПа

а

20-

15-

10-

Рис. 2. Зависимость критического значения давления нагнетания газа от абсолютной проницаемости пласта при разных значениях: а — начального давления в пласте ро = 3 МПа (1), ро = 5 МПа (2); Ь — начальной температуры пласта

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

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

Т0 = 280 К (1), Т0 = 285 К (2).

ратообразования.

Заключение

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

Литература

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

[2] Jadhawar, P. Subsurface carbon dioxide storage through clathrate hydrate formation / P. Jadhawar, A.H. Mohammadi, J. Yang // Advances in the Geological Storage of Carbon Dioxide. - Springer. Printed in the Netherlands. - 2006. - P. 111-126.

[3] Истомин, В.А. Газовые гидраты в природных условиях / В.А. Истомин, В.С. Якушев. - М.: Недра, 1992. - 236 c

[4] Нигматулин, Р.И. Автомодельная задача о разложении газогидра-тов в пористой среде при депрессии и нагреве / Р.И. Нигматулин,

B.Ш.Шагапов, В.Р. Сыртланов // ПМТФ. - 1998. - Т. 39. - №3. -

C. 111-118.

[5] Шагапов, В.Ш. Нагнетание газа в пористый резервуар, насыщенный газом и водой / В.Ш.Шагапов, Н.Г. Мусакаев, М.К.Хасанов // Теплофизика и аэромеханика. - 2005. - Т. 12. - №4. - С. 645-656.

Поступила в редакцию 15/XI/2007; в окончательном варианте — 15/XI/2007.

GAS INJECTION IN THE POROUS MEDIUM WITH FORMATION OF GASHYDRATE

© 2008 M.K. Khasanov2

The particularities of the formation gas hydrate at injection of the gas into porous media, in initial condition filled by gas and water are considered. The hydrate can be forming on frontal surface and in longitude area. The conditions, under which each hydrate forming mode is realized are defined.

Keywords: porous medium, formation of gashydrate.

Paper received 15/XI/2007. Paper accepted 15/XI/2007.

2Khasanov Marat Kamilovich (hasanovmk@mail.ru), Dept. of Applied Mathematics and Mechanics, Sterlitamak state teacher training academy, Sterlitamak, 453103, Russia