Вытеснение метана из газогидратного пласта при закачке диоксида углерода

2016 Математика и механика № 6

УДК 532.546:536.421 DOI 10.17223/19988621/44/9

В.Ш. Шагалов, М.К. Хасанов, Г.Р. Рафикова

ВЫТЕСНЕНИЕ МЕТАНА ИЗ ГАЗОГИДРАТНОГО ПЛАСТА ПРИ ЗАКАЧКЕ ДИОКСИДА УГЛЕРОДА1

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

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

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

Основными способами разработки газогидратных пластов являются депресси-онное воздействие на пласт и нагрев гидратосодержащих пород. Моделирование разложения газогидрата в пластах данными способами рассмотрены в работах [6-8]. Одним из последних инновационных способов извлечения метана из состава гидрата является его замещение диоксидом углерода в газогидрате. Предложенный метод позволяет одновременно добыть метан и законсервировать углекислый газ в необходимых объемах. Экспериментальные исследования процесса замещения метана из состава гидрата двуокисью углерода проведены авторами [9-11], которыми было установлено, что процесс может происходить без высвобождения свободной воды при термобарических условиях, соответствующих стабильному существованию гидрата метана.

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

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

Постановка задачи. Основные уравнения

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

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

Запишем уравнения сохранения масс для диоксида углерода и метана, фильтрующихся через пласт:

д

dt dx gd

= - J-*.

dKС), d() 7

dt cx gm.

Здесь m - пористость скелета (для объема, занятого газом и гидратом); Sg - газонасыщенность; p0 и ugi (i= d, m) - парциальные плотности и скорости компонент газовой смеси; Jgd и Jgm - интенсивности перехода диоксида углерода в состав гидрата и вытеснения метана из состава гидрата; нижние индексы i= d, m относятся соответственно к углекислому газу и метану.

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

Sh = Shd + Shm ,

где Shi - гидратонасыщенность, которая относится к гидратам диоксида углерода (i=d) и метана (i=m). При этом

Sg + Sh = 1 или Sg + Shd + Shm = i.

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

d

dt ~Jgd, dt ~ Jgm. (

Здесь p0 и Ghi (i=d и m) плотность гидрата и массовое содержание газа в составе

hi ш

гидрата для диоксида углерода (i=d) и метана (i=m).

Для гидратов диоксида углерода и метана [12] phd = 1117 кг/м3, p0hm = 910 кг/м3 и Gd =0.29, Gm =0.13, в связи с чем выполняется соотношение

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

числу молекул диоксида углерода, переходящих в состав гидрата. Отсюда, для интенсивностей Jgd и Jgm выполняется следующее соотношение:

^^ = Jgm

Md Mm

где Mi (i=d и m) -молекулярные массы диоксида углерода и метана. Для газовой смеси в целом введем среднемассовую скорость

0 0 ,0 /00,0\ ,-гч

^ ^ = ^ ^ + Р gm Ugm , =pgd +pgm ) .

Сложим уравнения (1), тогда с учетом (4) и (5) будем иметь

д(mSgр0) д(mSgP0ug) (м Л

Х / ^ +-= Jgd -1| .

д/ дx ' ^ у Md

Для фильтрации и диффузионного перемешивания газовой смеси применим соответственно закон Дарси и закон Фика:

mSgUg =& ;

др0

р^ =-P0gmwgm = ,

Здесь Dg - коэффициент диффузионного перемешивания смеси метана и диоксида углерода, wgi (i=d и m) - диффузионная скорость компонентов газовой смеси.

Будем полагать, что газовая смесь является калорически совершенной и подчиняется закону Дальтона

Pgd = р^ , Pgm = р0gm ^ , Г = Pgd + Pgm .

Приведенные выше уравнения дополним уравнением энергии

дТ „ 0 дТ д ( дТ Л ( 0 дShd „, д8ы Л

рс -Т+^ & = дХ Ы +т {^«Цт)

Рс = (1 - т)р! с*к + т (Sgр0 сg + Shdр°и сЫ + ShmPPhm ст) ,

р0 ^ = р^ сgd +Р0» сgm , Х = (1 - т)Хк + т (SgЪg + ShdЪhd + Shmkhm ) .

Здесь с, 1 - удельная теплоемкость и теплопроводность фаз (р = g, h), рс , 1 -удельная объемная теплоемкость и теплопроводность системы, , 1Ш - удельная теплота разложения и образования гидрата метана и диоксида углерода соответственно, отнесенные на единицу массы.

Систему уравнений (1), (2), (6) и (10) дополним граничными и начальными условиями:

х = o, / > 0: р = ре, pgd = Pgde, Т = те, х = х0, г > 0: p = То, дpgd/дх = 0, дТ/дх = 0,

0 < х < х0 , г = 0: p = Po, Pgd = 0, Т = To, где х0 - протяженность пористого пласта.

Кинетика вытеснения метана из состава газогидрата диоксидом углерода

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

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

Л2 Са, ,Л2

$Нт = 1

\тй)

=

=1 — а

Для процесса переноса диффундирующего газа через слой гидрата диоксида углерода запишем уравнение диффузии

Ф'

Яй

оы д

дг

дг

а(тй) > Г > ая

где - коэффициент диффузии углекислого газа в гидрате С02. При граничных условиях

Р0й'

о

Ряй (я

Р яй

' = 0, г = а,

квазистационарное решение уравнения (13) для плотности подвижного углеки слого газа в гидрате С02 имеет вид

. =р0 1П

Ряй = Р ей

\тй)

> г > а„

Здесь р

ей

1п

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

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

Запишем выражение для потока массы подвижного диоксида углерода к поверхности контакта между гидратами диоксида углерода и метана, отнесенные к единице ее площади как

1 ей

= - А

ы

дР

ей

дг

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

ей

0

Ряй ( я) АЬй

а(тй) 1п(а(тй)/ае ) '

0

Чтобы связать поток с интенсивностью потребления диоксида углерода

3й, необходимо определить полную площадь поверхности контакта в единице

объема пористой среды (удельная поверхность контакта) между гидратами СО2 и СН4. В рамках принятой схемы для этой площади можем записать

т

т = 2па

па

где П(тй) - число каналов радиусом а(тй) в единице объема пористой среды.

По аналогии с законом Генри примем, что насыщенная концентрация свободного газа СО2 в составе гидрата пропорциональна величине парциального давления СО2 в газовой фазе, т.е.

Р0й,)=ГйРей ,

где Г й - параметр Генри. Тогда для интенсивности 3потребления СО2 на образование гидрата в единице объема, связанной с как 3= , используя (17) и (18), получим

Р0 А

Рей (г) АЬй

3ей = 2т

а

Отметим, что р0^ ^ и Аы являются неизвестными эмпирическими параметрами. В выражение (20), а следовательно, и в общую систему уравнений данные параметры входят в виде произведения. Поэтому для дальнейшего введем один эмпирический параметр, приведенный коэффициент диффузии для СО2, отвечающий в силу вышепринятых допущений, за кинетику образования гидрата СО2 в виде

А = .

Рей

Подставляя (19) и первое уравнение из (9) в (21), получим

А = Аш ГйЯйГ.

Если, по аналогии с жидкостями, коэффициент диффузии Аш и параметр Генри Гй считать достаточно консервативными величинами (слабо зависящими от давления Рей , например) и поскольку относительное изменение температуры по шкале Кельвина для рассматриваемых процессов слабо меняется (дт/Т < 10-2), то введенный коэффициент А , имеющий размерность коэффициента диффузии, можно принять за постоянный эмпирический модельный параметр. Наличие небольшой доли примесей других газов в составе углекислого газа и метана может повлиять на процессы гидратообразования и разложения [1, 13, 14]. Эмпирический параметр А зависит как от структуры гидрата, так и от особенностей скелета пористой среды.

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

к £ = к

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

Результаты

На основе системы уравнений (1), (2), (6) и (10) с учетом (7) - (9) и (20) с начальными и граничными условиями (11) были проведены численные расчеты. Для параметров системы приняты следующие значения:

Т0 = 273 К, Тв = 273 К, р0 = 2.8 МПа, рв = 3.2 МПа, £йт0 = 0.15, х0 = 10 м, т = 0.1, к = 10-11 м2, рс = 2.5 -106 Дж/(м3-К), Х = 2 Вт/(мК), 1кт = 4.37-105 Дж/кг, 1М = 4.8-105 Дж/(кг-К), р0 = 910 кг/м3, р0 = 1117 кг/м3,

пгп пи Ьт кй

От = 0.13, Ой = 0.29 , К8т = 520 Дж/(кг-К), К*а = 189 Дж/(кг-К), с8т = 1560 Дж/(кг-К), с ^ = 731 Дж/(кг-К), = 1.12 X10-5 Пас, = 1.5-10-5 Па с, Б*, = 10-5 м2/с, Б = 10-15 м2/с.

На рис. 1 кривыми 1, 2, 3 и 4 представлены распределения общего и парциальных давлений для диоксида углерода и метана (а, Ъ, с), температуры (й), гидрато-насыщенностей С02 и СН4 (в, /) по координате в моменты времени 5, 25, 300 и 1900 с соответственно. Точечная кривая соответствует равновесной температуре при текущем давлении системы.

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

Массовые расходы углекислого газа и метана, отнесенные на единицу площади поперечного сечения пласта,

Эй = Р^*и , = Р^я"я , Чцт = Р°ттБг"я . Здесь д яйв - расход закачиваемого углекислого газа; дгй , д - расходы диоксида углерода и метана на выходе.

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

t t t

т*йв = | д , т*й = | д , т*т = | дгтй1:.

0 0 0

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

3.1 3.02.92.8-

Pgd, ЛПа

2.8.

2.1 1.4 0.7 4

2.5

—Г~

5.0

—1—

7.5

Т, К 275

274

273

1 2

"1-1-1-1-1 1 I

л, м 0 2.5 5.0 7.5 м

2.8 2.1 1.4

0.7

Рис. 1. Эволюция полей общего давления p

Fig. 1. Evolution of the (a) total pressure fields p, (b, c) partial pressure of carbon dioxide pgd and methane pgm, respectively, (d) temperature T, and (e, f) saturation of the carbon dioxide hydrates Shd and methane Shm, respectively

d

4

3

qgde, qgd, кг/(м2-с) mgde, mgd, кг/м2

a b

Рис. 2. Динамика массовых расходов и массы закачиваемого и выходящего углекислого газа (а, b) и извлекаемого метана (c, d) Fig. 2. Dynamics of the mass flow and mass of: (a, b) injected and discharged carbon dioxide and (c, d) recoverable methane

Видно, что на начальном этапе вытесняется свободный метан, затем происходит скачок расходов, часть диоксида углерода поглощается, и часть метана выделяется вследствие начала процесса газозамещения. Следующий этап характеризуется вытеснением метана, полученного в результате процесса замещения метана из состава гидрата углекислым газом. После полного извлечения метана из пласта расходы диоксида углерода на входе и выходе совпадают. В результате процесса извлекается около 35 кг/м2 метана и консервируется около 105 кг/м2 углекислого газа.

Заключение

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

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

ЛИТЕРАТУРА

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

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

3. Shagapov V.Sh., Musakaev N.G., Khasanov M.K. Formation of gas hydrates in a porous medium during an injection of cold gas // Int. J. Heat and Mass Transfer. 2015. V. 84. P. 1030. DOI: 10.1016/j.ijheatmasstransfer.2015.01.105.

4. Шагапов В.Ш., Хасанов М.К., Гималтдинов И.К., Столповский М.В. Численное моделирование образования газогидрата в пористом пласте конечной протяженности при продувке его газом // Прикладная математика и техническая физика. 2011. Т. 52. № 4. С. 116.

5. Рафикова Г.Р. Образование газогидрата в замкнутом объеме, заполненном водонасы-щенной пористой средой // Вестник Кемеровского государственного университета. 2015. Вып. 2(62). Т. 2. С. 122.

6. Васильев В.И., Попов В.В., Цыпкин Г.Г. Численное исследование разложения газовых гидратов, сосуществующих с газом в природных пластах // Изв. РАН. Механика жидкости и газа. 2006. № 4. С. 127.

7. Шагапов В.Ш., Хасанов М.К., Гималтдинов И.К., Столповский М.В. Особенности разложения газовых гидратов в пористых средах при нагнетании теплого газа // Теплофизика и аэромеханика. 2013. Т. 20. № 3. С. 347.

8. Шагапов В.Ш., Насырова Л.А. О нагреве пористой среды, частично заполненной газогидратом, в условиях непроницаемых границ // Теплофизика высоких температур. 1999. Т. 37. № 5. С. 784.

9. McGrail B.P., Zhu Г., Hunter R.B., White M.D., Patil S.L., Kulkarni A.S. A New method for enhanced production of gas hydrates with CO2 // AAPG Hedberg Conference «Gas hydrates: energy resource potential and associated geologic hazards». Vancouver, Canada, 2004.

10. Espinoza D.N., Santamarina J.C. P-wave monitoring of hydrate-bearing sand during CH4-CO2 replacement // Int. J. Greenhouse Gas Control. 2011. V. 5. P. 1032. DOI: 10.1016/j.ijggc.2011.02.006.

11. Jung J. W., Nicolas Espinoza D. and Carlos Santamarina J. Properties and phenomena relevant to CH4-CO2 replacement in hydrate-bearing sediments // J. Geophysical Research. 2012. V. 115. B10102. DOI: 10.1029/2009JB000812.

12. Макогон Ю. Ф. Газогидраты. История изучения и перспективы освоения // Геология и полезные ископаемые Мирового океана. 2010. № 2. С. 5.

13. Rehder G., Brewer P.W., Peltzer E.T., Friederich G. Enhanced lifetime of methane bubble streams within the deep ocean // Geophysical Research Letters. 2002. No. 29. P. 21. DOI: 10.1029/2001GL013966.

14. 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? // J. Geophysical Research. 2006. V. 111. P. 382. DOI: 10.1029/2005JC003183.

Статья поступила 20.06.16 г.

Shagapov V. Sh., Khasanov M. K., Rafikova G. R. DISPLACEMENT OF METHANE FROM A GAS HYDRATE RESERVOIR IN THE PROCESS OF CARBON DIOXIDE INJECTION. Tomsk State University Journal of Mathematics and Mechanics. 6(44). pp. 104-114

DOI 10.17223/19988621/44/9

In this paper, a numerical simulation of the process of methane displacement from a gas hydrate reservoir during the injection of carbon dioxide followed by methane substitution by carbon dioxide in the composition of the hydrate was performed. The process of gas replacement was assumed to occur in a stable field of methane hydrates and carbon dioxide with CO2 in a gaseous state. The intensity of carbon dioxide hydrate formation was determined by diffusion of carbon dioxide across the hydrate layer formed between a stream of the gas mixture and the methane hydrate in the porous channels.

The general pressure of the system, the partial pressures of the carbon dioxide and methane, the saturations of hydrates of carbon dioxide and methane, and the temperature as functions of time were obtained. In the initial stage of the process, free methane was forced out, and subsequently the methane received after the gas replacement process was extracted. This resulted in a complete transition of methane hydrate into carbon dioxide hydrate and a complete extraction of methane from the porous reservoir.

Keywords: gas hydrate reservoir, substitution of methane with carbon dioxide in the hydrate composition, diffusion, formation of gas hydrate.

SHAGAPOV Vladislav Shaykhulagzamovich (Academician of Academy of Sciences of RB, Doctor of Physics and Mathematics, Birsk Branch of Bashkir State University, Birsk, Russian Federation; Researcher Institute of Mechanics and Engineering, Kazan, Russian Federation; Researcher Institute of Mechanics, Kazan, Russian Federation). E-mail: Shagapov@rambler.ru

KHASANOV Marat Kamilovich (Candidate of Physics and Mathematics, Sterlitamak Branch of Bashkir State University, Sterlitamak, Russian Federation). E-mail: hasanovmk@mail.ru

RAFIKOVA Guzal Rinatovna (Birsk Branch of Bashkir State University, Birsk, Russian Federation).

E-mail: rafikova_guzal@mail.ru

REFERENCES

1. Makogon Yu.F. (1974) Gidraty prirodnykh gazov [Hydrates of natural gases]. Moscow: Nedra.

2. Istomin V.A., Yakushev V.S. (1992) Gazovye gidraty v prirodnykh usloviyakh [Gas hydrates in the natural environment]. Moscow: Nedra.

3. Shagapov V.Sh., Musakaev N.G., Khasanov M.K. (2015) Formation of gas hydrates in a porous medium during an injection of cold gas. International Journal of Heat and Mass Transfer. 84. p. 1030. DOI 10.1016/j.ijheatmasstransfer.2015.01.105.

4. Shagapov V.Sh., Khasanov M.K., Stolpovskii M.V., Gimaltdinov I.K. (2011) Numerical modeling of formation of a gas hydrate in a finite-length porous bed purged by a gas. Journal of Applied Mechanics and Technical Physics. 52(4). p. 599. DOI 10.1134/ S0021894411040134.

5. Rafikova G.R. (2015) Obrazovanie gazogidrata v zamknutom ob'eme, zapolnennom vodonasyshchennoy poristoy sredoy [Formation of gas hydrate in a closed volume filled with porous water-saturated media]. Vestnik Kemerovskogo gosudarstvennogo universiteta -Bulletin of Kemerovo State University. 2(62). p. 122.

6. Vasil'ev V.I., Popov V.V., Tsypkin G.G. (2006) Chislennoe issledovanie razlozheniya gazovykh gidratov, sosushchestvuyushchikh s gazom v prirodnykh plastakh [Numerical study

of the decomposition of gas hydrates coexisting with a gas in natural reservoirs], Izv. Ross. Akad. Nauk. Mekhanika zhidkosti i gaza - Fluid Dynamics. 4, pp, 127,

7, Shagapov V.Sh., Khasanov M.K., Gimaltdinov I.K., Stolpovsky M.V. (2013) The features of gas hydrate dissociation in porous media at warm gas injection, Thermophysics and Aeromechanics. 20(3), p, 339, DOI 10,1134/S0869864313030104,

8, Shagapov V,Sh,, Nasyrova L,A, (1999) O nagreve poristoy sredy, chastichno zapolnennoy gazogidratom, v usloviyakh nepronitsaemykh granits [The heating of porous medium partly filled with gas hydrate in the presence of impermeable boundaries], Teplofizika vysokikh temperatur - High Temperature. 37(5), p, 784,

9, McGrail B,P,, Zhu T,, Hunter R,B,, White M,D,, Patil S,L,, Kulkarni A,S, (2004) A new method for enhanced production of gas hydrates with CO2, AAPG Hedberg Conference «Gas hydrates: energy resource potential and associated geologic hazards», Vancouver,

10, Espinoza D,N,, Santamarina J,C, (2011) P-wave monitoring of hydrate-bearing sand during CH4-CO2 replacement, International Journal of Greenhouse Gas Control. 5, p, 1032, DOI 10,1016/j,ijggc,2011,02,006,

11, Jung J,W,, Nicolas Espinoza D,, Carlos Santamarina J, (2010) Properties and phenomena relevant to CH4-CO2 replacement in hydrate-bearing sediments, Journal of Geophysical Research. 115, B10102, DOI 10,1029/2009JB000812,

12, Makogon Yu,F, (2010) Gazogidraty, Istoriya izucheniya i perspektivy osvoeniya [Gas hydrates, The history of exploration and development prospects], Geologiya i poleznye iskopaemye Mirovogo okeana, 2, p, 5,

13, Rehder G, Brewer P,W,, Peltzer E,T,, Friederich G, (2002) Ehnaced lifetime of methane bubble streams within the deep ocean, Geophysical Research Letters. 29, p, 21, DOI 10,1029/2001GL013966,

14, MgGinnis D,F,, Greinert J,, Artemov Y,, Beaubien S, E,, Wuest A, (2006) Fate of rising methane bubbles in stratified waters: How much methane reaches the atmosphere? Journal of Geophysical Research, 111, p, 382, DOI 10,1029/2005JC003183,