Нагнетание газа в водонасыщенный пористый пласт с образованием гидрата в диффузионном режиме

Вестник Самарского университета. Естественнонаучная серия. 2016. № 1-2 51

МЕХАНИКА

УДК 532.546:536.421

Г.Р. Рафикова, В.Ш. Шагапов1

НАГНЕТАНИЕ ГАЗА В ВОДОНАСЫЩЕННЫЙ ПОРИСТЫЙ ПЛАСТ С ОБРАЗОВАНИЕМ ГИДРАТА В ДИФФУЗИОННОМ РЕЖИМЕ2

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

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

Введение

Гидратное состояние является выгодным для консервации и хранения углеводородных газов. В 1 м3 газогидрата метана при T=0oC и p=3 МПа содержится около 160 м3 метана. Для того чтобы такую массу газа хранить в свободном состоянии в одном кубическом метре, давление необходимо повышать более чем до 16 МПа. Поэтому представляется привлекательным создавать такие условия хранения газов в подземных резервуарах, чтобы значительная его часть находилась в составе газогидрата.

Исследования процессов образования и разложения газовых гидратов в пористых средах проводились такими авторами, как Г.Г. Цыпкин, Э.А. Бондарев, И.К. Гималтдинов, М.К. Хасанов, О.Р. Нурисламов и др. [1-5]. В работах получены автомодельные и численные решения, позволяющие изучать процессы образования и разложения газогидратов в пористых пластах с движущейся фрон-

!© Рафикова Г.Р., Шагапов В.Ш., 2016

Рафикова Гузаль Ринатовна (rafikova_guzal@mail.ru), кафедра высшей и прикладной математики, Бирский филиал Башкирского государственного университета, 452453, г. Бирск, Интернациональная, 10.

Шагапов Владислав Шайхулагзамович (shagapov@rambler.ru), лаборатория «Механика сплошной среды», Институт механики и машиностроения Казанского научного центра РАН, 420111, г. Казань, Лобачевского, 2/31.

2Исследование выполнено за счет гранта Российского научного фонда (проект №15-11-20022).

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

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

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

P°sk = pW = Р0 = const, Usk = Vw = vh = 0.

Здесь и в дальнейшем нижние индексы sk, wh, h и g соответствуют параметрам, относящимся к скелету, воде, газогидрату и газу.

Пусть m и Si — пористость скелета и относительная объемная насыщенность пор водой, гидратом и газом (i = w,h,g) . В рамках вышеперечисленных допущений уравнения сохранения масс для плоскоодномерного течения запишем в виде:

д

dt = Jw' dt = Jh dt + dx = Jg.

Здесь vg - скорость газа, Jw и Jg - интенсивности потребления воды и газа на образование гидрата с интенсивностью Jh в единице объема пористой среды. Газогидрат является клатратным соединением [6], для которого на одну молекулу метана приходится шесть молекул воды. Поэтому массовое содержание метана G в составе гидрата постоянное число ( G « 0.12). В этой связи для Jw, Jg и Jh должны выполняться следующие кинематические соотношения:

Jw =

С учетом этих соотношений из первого и второго уравнений (2) получим

mSwpW +

Полагая, что для исходного состояния

На основе кинематического условия

Sw + Sh + Sg = 1

с учетом (5) можем получить

с* _ (1 — Sw0 — Sg) q _ q ph f! ГЛ Q

Sh = ^-f.--Г, Sw = Sw0--

(1 — ph (1 — G)/pw ) pw

Следовательно, текущие величины гидратонасыщенности Sh и водонасыщенно-сти Sw однозначно определяются значением газонасыщенности Sg.

Полагая, что температура для скелета пористой среды и всех ее составляющих совпадает (Tsk =Tw =Th =Tg=T), и пренебрегая баротермическим эффектом по

сравнению с тепловыми эффектами за счет гидратообразования, запишем уравнение притока тепла:

дТ о д^Т д2Т

— + р с3 шБд ид — = Лдх^

Рс^7 + Р°° сд тБд = + Ы,

рс = (1 - т) р^Свк + т Р03]С3, Л = (1 - т) Лзк + т ^ Б^Л^ .

j=w,g,h j=w,g,h

Здесь ^, Лj - удельная теплоемкость и теплопроводность фаз (у = 1л,д,к) и рс и Л - удельная объемная теплоемкость и теплопроводность системы. Второе слагаемое в правой части (8) выражает тепловой эффект гидратообразования, где I - удельная теплота, отнесенная на единицу массы гидрата, которая выделяется при его образовании.

Для газовой фазы примем уравнение Клайперона — Менделеева

Р = Р0 К Т.

Для процесса фильтрации примем обобщенный закон Дарси

",Бд ид = - £ дР-

где цд — динамическая вязкость газа, т — "живая" проницаемость, зависящая от "живой" пористости тБд и характерных "живых" радиусов пор.

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

Т = То, р = р8

При Ь ^ 0 через левую границу пласта (х=0) происходит нагнетание газа при давлении ре с температурой газа Те, соответствующим термобарическим условиям образования гидрата.

Т = Те, р = ре,

Будем полагать, что потоки массы и тепла через правую границу (х = хо) отсутствуют

дТ др

— =0, дх =0

2. Кинетика образования гидрата

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

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

Б- = 1 -

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

дЬ г дг

' ^ д (др1

Далее будем использовать квазистационарное решение, полагая др!' /дЬ = 0 и применяя следующие граничные условия:

р!' = Р1м > г = а

получаем решение уравнения (15) в виде:

, 0 1п(г/а„) ( ч

р - = Р^Ш^^!-/-Г ' <г<0,д) ,

1 ;1п

здесь - средняя плотность подвижной жидкости в составе гидрата.

Отметим, что второе условие из (16) означает мгновенный переход подвижной воды в состав гидрата при достижении границы контакта с газом.

Выражение для потока массы подвижной воды к поверхности контакта между

газом и гидратом, отнесенного к единице ее площади, запишем в виде:

=др-)

\ /ад

Подставляя сюда решение (17), получим

=

ш ;

ад 1п \а

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

шБд

Яд = 2пад Пд , Пд = -^ ,

па2

где пд - число каналов радиусом ад в единице объема пористой среды. Тогда из выражений (20) для величины удельной поверхности контакта имеем

2тБ.

=

д

Для интенсивности потребления воды на образование гидрата в единице объема газа, связанной с как ЛШ = яд, с учетом (19) и (21) получим:

р1,,

, \ ш

2тБд 2 ""-^

а

а

д

Здесь отметим, что рW(s) и Dw являются неизвестными эмпирическими параметрами, которые в полученные решения для интенсивности потребления воды входят в виде произведения.

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

Dw ро( ,

D = р0""

Выражение (22) с учетом (23) примет вид:

ро D

^ = 2тБд 2 а / . .

ад 1п

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

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

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

кд = к

где к - коэффициент абсолютной проницаемости скелета.

3. Результаты исследования

Для параметров, характеризующих систему "пористая среда - газогидрат -газ - вода" здесь и в дальнейшем приняты следующие значения m = 0.1, к = = 10~9 м2, рд = 10~5 Па-с, pc = 2.5-106 Дж/(кг-К), Л = 2 Вт/(м-К), l = = 5105 Дж/кг, p0h = 900 кг/м3, G =0.12, p0w = 1000 кг/м3, Rg = 520 Дж/(кг-К), cg = 1560 Дж/(кг-К), D = 510^11 м2/с. Протяженность пористого пласта берем равной xo = 100 м. Отметим, что принятые значения для приведенного коэффициента D по порядку величин близки к опытным данным [8, 12] по интенсивности роста гидратной корки при миграции метанового пузырька в морской воде, в термобарических условиях образования гидрата.

Для исходных и граничных температур и давлений принято T0 = 273 К, Te = = 273 К и po = ps (To), pe = 10 МПа, начальная водонасыщенность полагалась равной Sw0 = 0.3.

На рис. 1 линиями 1, 2 и 3 представлены распределения давления (а), температуры (б) и гидратонасыщенности (в) в моменты времени t =0.5, 1 и 3 с соответственно. Здесь и в дальнейшем штриховая линия соответствует равновесной температуре для текущих значений давлений пласта.

Рис. 1. Распределения давления р

Видно, что при данном режиме в пористой среде наблюдается достижение однородного давления р = ре (за секунды), то есть за достаточно короткий период устанавливается гомобарический режим [13]. Текущая температура при этом не успевает проследить за равновесным значением температуры для текущего давления. Затем в остальное время после достижения давления в пористом пласте граничного значения р = ре накопление гидратной массы происходит в гомобари-ческом режиме.

На рис. 2 представлены осциллограммы давления (а), температуры (б) и гид-рато- и водонасыщенности (в) в точке с координатой х =50 м при различных значениях водонасыщенности БWо =0.3 и 0.5 (линия 1 и 2 соответственно). Точечные линии иллюстрируют эволюцию водонасыщенности.

При начальной водонасыщенности Б^ =0.3 гидратообразование проходит интенсивнее (примерно в полтора раза) по сравнению со значением большей водо-

насыщенности Б^ =0.5, что объясняется более высокой удельной площадью контакта между гидратом и газом.

Рис. 2. Динамика давления р (а), температуры Т (б) и гидратонасыщенности (в) в точке х = 50 м для различных значений начальной водонасыщенности

Sw0

На рис. 3 линиями 1, 2, 3 представлена динамика давления, температуры, гидратонасыщенности в точке с координатой х =50 м при начальном значении водонасыщенности Б^ =0.3 при следующих значениях приведенного коэффициента диффузии В = 5-10~10, 5-10~и и 5-10 12 м2/с соответственно. Правая часть рис. 3 представляет собой увеличенные фрагменты для начального этапа.

Определяющими параметрами характерных времен фильтрации 1(р) (время выравнивания давления в пласте до граничного значения) и диффузии воды (накопление гидрата) являются коэффициент пьезопроводности х(р) (р) ~ х0/х(р)) и диффузии В « а?/В). Из рис. 3 следует, что если для коэффициента диффузии принять В = 5-10~10 м2/с за время полного перехода воды в состав гидрата, давление в резервуаре не успевает достичь значения р = ре. При дальнейшем

Рис. 3. Осцилограммы давления p

уменьшении коэффициента диффузии (D ~ 510-11 м2/с) наблюдается примерно одинаковое значение характерных времен t(D) « t(p). При малых значениях коэффициента диффузии процесс установления однородного давления (p « pe) будет происходить гораздо быстрее процесса полного перехода воды в пористой среде в состав гидрата t(D) <t(p).

Заключение

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

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

Литература

[1] Нурисламов О.Р., Шагапов В.Ш. Нагнетание газа во влажную пористую среду с образованием газогидрата // Прикладная математика и механика. 2009. Т. 73. № 5. С. 809-823.

[2] Особенности разложения газовых гидратов в пористых средах при нагнетании теплого газа / В.Ш. Шагапов [и др.] // Теплофизика и аэромеханика. 2013. Т. 20. № 3. С. 347-354.

[3] Численное моделирование образования газогидрата в пористом пласте конечной протяженности при продувке его газом / В.Ш. Шагапов [и др.] // Прикладная математика и техническая физика. 2011. Т. 52. № 4. С. 116-126.

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

[5] Бондарев Э.А., Максимов А.М., Цыпкин Г.Г. К математическому моделированию диссоциации газовых гидратов // Докл. АН СССР. 1989. Т. 308. № 3. С. 575-578.

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

[7] Макогон Ю.Ф. Гидраты природных газов. М.: Недра, 1974. 208 с.

[8] Ehnaced lifetime of methane bubble streams within the deep ocean / G. Rehder [et al.] // Geophysical research letters. 2002. № 29. P. 21.

[9] Fate of rising methane bubbles in stratified waters: How much methane reaches the atmosphere? / D.F. MgGinnis [et al.] // Journal of geophysical research. 2006. V. 111. P. 382.

[10] Булейко В.М., Вовчук Г.А., Григорьев Б.А. Экспериментальное исследование термодинамических свойств гидратов углеводородов алканового ряда // Вести газовой науки. 2012. № 3(11). С. 282.

[11] Нестеров А.Н. Применение поверхностно-активных веществ для интенсификации процессов образования гидратов в технологиях транспорта и хранения газа // Современное состояние газогидратных исследований в мире и практические результаты для газовой промышленности. М: ООО ИРЦ Газпром, 2004. С. 66.

[12] Тазетдинов Б.И. Особенности образования и разложения газогидратов в водных и газовых средах: дис. ... канд. физ.-мат. наук. Уфа: БашГУ, 2014. 113 с.

[13] Нигматулин Р.И. Динамика многофазных сред: в 2 ч. М: Наука, 1987. Ч. 1. 464 с. Ч. 2. 359 с.;

References

[1] Nurislamov O.R., Shagapov V.Sh. Nagnetanie gaza vo vlazhnuiu poristuiu sredu s obrazovaniem gazogidrata [Gas injection into a moist porous medium with the formation of a gas hydrate]. Prikladnaia matematika i mekhanika [Applied Mathematics and Mechanics], 2009, V. 73, no. 5, pp. 809-823 [in Russian].

[2] Shagapov V.Sh., Hasanov M.K., Gimaltdinov I.K., Stolpovskii M.V. Osobennosti razlozheniia gazovykh gidratov v poristykh sredakh pri nagnetanii teplogo gaza [The features of gas hydrate dissociation in porous media at warm gas injection]. Teplofizika i aeromekhanika [Thermophysics and Aeromechanics], 2013, V. 20, no. 3, pp. 347-354 [in Russian].

[3] Shagapov V.Sh., Hasanov M.K., Gimaltdinov I.K., Stolpovskii M.V. Chislennoe modelirovanie obrazovaniia gazogidrata v poristom plaste konechnoi protiazhennosti pri produvke ego gazom [Numerical modeling of formation of a gas hydrate in a finite-length porous bed purged by a gas]. Prikladnaia matematika i tekhnicheskaia fizika [Journal of Applied Mechanics and Technical Physics], 2011, V. 52, no. 4, pp. 116-126 [in Russian].

[4] Chypkin G.G. Techeniia s fazovymi perekhodami v poristykh sredakh [Flows with phase transitions in porous media]. M.: Fizmatlit, 2009, 232 p. [in Russian].

[5] Bondarev E.A., Maksimov A.M., Tsypkin G.G. K matematicheskomu modelirovaniiu dissotsiatsii gazovykh gidratov [On the mathematical modeling of dissociation of gas hydrates]. Dokl. AN SSSR [Proceedings of the USSR Academy of Sciences], 1989, V. 308, no. 3, pp. 575-578 [in Russian].

[6] Istomin V.A., Yakushev V.S. Gazovye gidraty v prirodnykh usloviiakh [Gas hydrates in nature]. M.: Nedra, 1992, 236 p. [in Russian].

[7] Makogon Yu.F. Gidraty prirodnykh gazov [Hydrates of natural gases]. M.: Nedra, 1974, 208 p. [in Russian].

[8] Rehder G, Brewer P.W., Peltzer E.T., Friederich G. Ehnaced lifetime of methane bubble streams within the deep ocean. Geophysical research letters, 2002, no. 29, pp. 21 [in English].

[9] MgGinnis D.F., Greinert J., Artemov Y., Beaubien S.E., Wuest A. Fate of rising methane bubbles in stratified waters: How much methane reaches the atmosphere? Journal of geophysical research, 2006, V. 111, pp. 382 [in English].

[10] Buleiko V.M., Vovchuk G.A., Grigoriev B.A. Eksperimental'noe issledovanie termodinamicheskikh svoistv gidratov uglevodorodov alkanovogo riada [Experimental study of thermodynamic properties of hydrocarbons alkane series]. Vesti gazovoi nauki [Gas Science News], 2012, no. 3(11), pp. 282 [in Russian].

[11] Nesterov A.N. Primenenie poverkhnostno-aktivnykh veshchestv dlia intensifikatsii protsessov obrazovaniia gidratov v tekhnologiiakh transporta i khraneniia gaza [Use of surfactants for the intensification of the processes of formation of hydrates in the technologies of transportation and storage of gas] in Sovremennoe sostoianie gazogidratnykh issledovanii v mire i prakticheskie rezul'taty dlia gazovoi promyshlennosti [Modern state of gas hydrates researches in the world and practical results for gas industry]. M.: OOO IRTs Gazprom, 2004, pp. 66 [in Russian].

[12] Tazetdinov B.I. Osobennosti obrazovaniia i razlozheniia gazogidratov v vodnykh i gazovykh sredakh: dis. . . . kand. fiz.-mat. nauk [Features of formation and decomposition of gas hydrates in aqueous and gaseous environments: Candidate's of Physical and Mathematical Sciences thesis]. U.: BashGU, 2014, 113 p. [in Russian].

[13] Nigmatulin R.I. Dinamika mnogofaznykh sred. V 2 ch. [Dynamics of multiphase medium. In 2 parts]. M.: Nauka, 1987. Part 1, 464 p. Part 2, 359 p. [in Russian].

G.R. Rafikova, V.Sh. Shagapov3

INJECTING GAS INTO THE WATER-FILLED POROSITY RESERVOIRS WITH THE FORMATION OF HYDRATES IN THE DIFFUSION MODE4

The paper considers the problem of injecting gas into the porous layer initially saturated with gas and water, accompanied by formation of hydrate. We have investigated case where the intensity of hydrate is limited by diffusion of water through the hydrate layer formed between water and gas in the core of the pore of the channel. As part of this scheme kinetics of hydrate is determined one empirical parameter D having the dimension of the diffusion coefficient (m2/s). We have analyzed effect of the value of this parameter on the features of the process of hydrate formation, depending on the parameters that define the initial state of porous reservoir and reservoir characteristics.

Key words: porous media, filtration, gas hydrate, diffusion, formation of hydrate.

Статья поступила в редакцию 3/Л/2016. The article received 3/II/2016.

3 Rafikova Guzal Rinatovna (rafikova_guzal@mail.ru), Department of Higher and Applied Mathematics, Birsk branch of the Bashkir State University, 10, Internasionalnaia Street, Birsk, 452453, Russian Federation.

Shagapov Vladislav Shayhulagzamovich (shagapov@rambler.ru), laboratory "Continuum mechanics", Institute of Mechanics and Engineering, Kazan Scientific Center, Russian Academy of Sciences, 2/31, Lobachevskiy Street, Kazan, Russian Federation.

4The study was performed by a grant from the Russian Science Foundation (project No.15-11-20022).