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

2018 Математика и механика № 55

УДК 519.6:532.546

DOI 10.17223/19988621/55/10

М.К. Хасанов

МАТЕМАТИЧЕСКАЯ МОДЕЛЬ ИНЖЕКЦИИ СЕРОВОДОРОДА В ПРИРОДНЫЙ ПЛАСТ, СОПРОВОЖДАЮЩЕЙСЯ ОБРАЗОВАНИЕМ ЕГО ГАЗОГИДРАТА1

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

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

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

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

1 Работа выполнена при финансовой поддержке РФФИ и Республики Башкортостан (проект 17-48020123 р_а).

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

Термобарические условия существования газогидрата И28 показаны на фазовой диаграмме (рис. 1) [5]. На данной диаграмме кривая 1 определяет трехфазное равновесие между газообразным сероводородом, водой и газогидратом И28, кривая 2 - равновесие между жидким сероводородом, водой и газогидратом И28, а кривая 3 - двухфазное равновесие между газообразным и жидким сероводородом. Соответственно газогидрат И28 существует левее кривых 1 и 3. В квадрупольной точке Q (Тд = 302.6 К и рд = 2.24 МПа) все четыре указанные фазы находятся в равновесии.

Т, К

Рис. 1. Фазовая диаграмма системы «сероводород - вода» Fig. 1. Phase diagram of the system "hydrogen sulfide - water"

Пусть полубесконечный горизонтальный пористый пласт (занимающий полупространство х > 0) в начальный момент времени насыщен водой и нефтью. Будем полагать, что начальная температура пласта T0 выше температуры Tq , соответствующей квадрупольной точке. Следовательно, в рассматриваемой задаче начальное состояние пласта не соответствует условиям образования газогидрата сероводорода. Положим, что через границу (х=0) закачивается жидкий сероводород, давление pe и температура Te которого соответствуют условиям существования газогидрата сероводорода.

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

сти: между первой и второй областями, где вода полностью переходит в газогид-ратное состояние (фронт гидратообразования), и между второй и третьей областями, где происходит вытеснение нефти сероводородом (фронт вытеснения).

Pe, Te /

— Гидрат H2S, H2S Вода, H2S Вода, нефть

-> К

x

x

Рис. 2. Схема задачи Fig. 2. Schematic diagram of the problem

x

0

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

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

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

д д

—{р^)+—(рт£г и-) = ^

дт дт д(.дг \

рс—+рстЬи1 — = —|л—|,

д1 дх дх V дх)

о кг др

тЬ1 иг = —---,

цг- дх

Рг = Ро- (1 + Рг (Р - Ро )) . Здесь t - время; х - координата; т - пористость; Т— температура; р — давление; нижние индексы г = I, s относятся соответственно к параметрам нефти и сероводорода; рг - плотность; иг - действительная средняя скорость; кг - фазовая проницаемость; цг- - динамическая вязкость; сг - удельная теплоемкость; Бг - насыщенность; Рг - коэффициент сжимаемости; рс и X - эффективные значения объемной

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

Зависимость коэффициента фазовой проницаемости кг от насыщенности Бг и абсолютной проницаемости пласта к0 зададим следующим образом:

к, = ко(г=s, 1).

На поверхности х = х(п), разделяющей первую и вторую области, происходит полный переход воды в гидратное состояние. Поэтому из условий баланса массы

и тепла на этой границе следует

^(1) дР(1) + /ст дР21 = т5 +Рй(1-£) - ^

ц дх Ц дх Н 1р0* ри 1 (И)'

т5йРй

дТ

= т^Р^»).

Здесь V» - скорость движения фронта образования газогидрата И28, О - массовая концентрация сероводорода в гидрате, Ьн - удельная теплота образования газогидрата И28, - начальная водонасыщенность пористой среды, - плотность воды, к,(1) - фазовая проницаемость для сероводорода в первой области, к^) - фазовая проницаемость для сероводорода во второй области. Нижний индекс п относится к параметрам на границе, разделяющей первую и вторую области, а нижние индексы 1 и 2 - к параметрам первой и второй областей соответственно. Температуру на этой границе будем считать непрерывной и равной температуре квад-рупольной точки Т(^.

На поверхности х = Х(ё), разделяющей вторую и третью области, происходит

вытеснение нефти сероводородом. Поэтому с учетом условий баланса массы нефти и сероводорода, а также баланса тепла на этой границе имеем

кя

^ дх

= т

--^^Л = т

Ц1 дх

дТ

X—^-X—^ = 0,

дх дх

где V ё) - скорость движения фронта вытеснения нефти сероводородом, к1 - фазовая проницаемость для нефти в третьей области. Нижний индекс ё относится к параметрам на границе, разделяющей вторую и третью области, а нижний индекс 3 - к параметрам третьей области.

На обеих границах давление и температуру будем считать непрерывными величинами.

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

=

Ри>5и>0

Рн

На основе уравнений системы (1), уравнения пьезопроводности и температуропроводности запишутся в виде

Т^!^), '=3;

Цг^Ш^'Ш ' = 1 2. 3.

(Р) ^0) (Р) к<2) (р) к1

где х =_—_ х =_—_ х =___

Х(1) ц,т(1 -)р, ' Х(2) ц,т(1 -^о)Р, ' Х(3) Цт(1 -^) '

(Т) = ^ = 2р0scsks (1) = 2роЛкя(2) = 2ро[С1к[

Х = , "(!) = Т О , "(2) = л о , а(3) = ^ о ,

нижние индексы '=1,2,3 относятся к параметрам соответственно первой, второй и третьей областей.

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

Начальные и граничные условия имеют вид

t = о: ^ = Т = То, р = ро (х>о).

х = о: Т = Тв, р = ре

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

Введем автомодельную переменную: р = х/-у/Х^. Для этой переменной уравнения пьезопроводности и температуропроводности (8), (9) в каждой области примут вид

р ар = 2П(') ар

(л

йр

а р

V ^ )

= а

ар

где П(') =х(р)У х(Т).

После интегрирования (11), (12) с учетом начальных и граничных условий (1о) решения для распределения давления и температуры в каждой из областей могут быть записаны в виде

р р +

р

р = р +

р

Н р

р

1

где введены обозначения интегралов:

Ь / %2 ) Ь

I(а,Ь;с) = | ехр--ё% ; 1(^ (а,Ь) = |ехр

2

4с

"а

ё % ;

Р(п) и Т(п) - давление и температура на границе между первой и второй областями; Р(О) и Т(ё) - давление и температура на границе между второй и третьей областями.

На основе соотношений (2), (4) с учетом решений (13), (14) получим уравнения для определения координаты фронта гидратообразования %(п) и значений параметров на нем:

(() - Р(п) )ехр

( %2 А

Чп)

4п

(Р(п) - Ре )еХР

( %2 А

Чп)

Ъ(2) "

1

--к,

4п

*

(Т(п) -Те)еХР

=(п) -

4 а

() -Т(п))еХР

1

з

(п) :

_тРй^й 2рс

5 %

йЬ (п) '

*0)

Т(п) = Тд,

где * = тц,Х(Т) + -1).

V Р*0 Pw )

Аналогично на основе соотношений (5) - (7) с учетом решений (14) и (15) получим систему уравнений для определения координаты фронта вытеснения нефти сероводородом и значений параметров на нем:

(Р(п) - Р(ё) )ехр

2

%

ч 4п,

'х(2) I (%(ё); П(2)) (() - Р0 )ехР

= тц,х(Т) (1 -^0)%,

(ё)'

2

%

4п

= тц х(Т) (1 - ^ )%<*);

1

(Т(й) -Т(п) )еХР

Г % 2 %

" а

(Т0 -Т(ё) )ехР

=

"а

*

'(3) (%(й ),

= 0.

Система граничных уравнений (16) - (21) в работе решена следующим образом. Вначале задается нулевое приближение искомых величин на фронте гидратообразования. Далее, решая уравнение (19), находим величину р^ (как функцию %(ё)), подставляя которую в уравнение (20) получим трансцендентное уравнение

для нахождения %(ё). Решив данное уравнение (методом половинного деления),

определим величину %(ё) (и соответственно Р(^), а затем из (21) определяем Т^.

Далее, подставляя (18) в уравнение (17), получим трансцендентное уравнение для нахождения %(п). Решая данное уравнение (методом половинного деления), определяем новое приближение величины %(п). Затем, решая уравнение (16), находим

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

На рис. 3 приведены распределения температуры и давления в зависимости от автомодельных переменных. Здесь и далее, если не оговорено иное, для параметров, характеризующих систему, приняты следующие значения: т = 0,1,

= 0.2, pe = 8.6 МПа, Te = 290 К, p0 = 8 МПа, T0 = 305 К, k0 = 10-11 м2, G = 0.24, р7 = 110-9 Па-1, р, = 310-9 Па-1, X = 2 Вт/(м-К), рс = 2106 Дж/(К-кг), V = 210-3 Па-с, ji, = 210-4 Па-с, ри = 1003 кг/м3, pw = 1000 кг/м3, р0/ = 900 кг/м3, Р0, = 890 кг/м3, a = 1900 Дж/(К-кг), Cs = 1800 Дж/(К-кг), Lh = 4.1105 Дж/кг.

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

Т, К и

300-

295-

290

0

2:з 5

р, МПа

8.59-

8.58

0 1 5(») 2 5(d) 3 5

Рис. 3. Распределения температуры и давления Fig. 3. Temperature and pressure distributions

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

коэффициента пьезопроводности х(р) =—~— значительно превышает величину

коэффициента температуропроводности х(Т) = —. Кроме того, согласно рис. 3,

рс

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

На рис. 4 приведены зависимости координат фронтов образования газогидрата H2S (кривая 1) и вытеснения нефти жидким сероводородом (кривая 2) от проницаемости пласта (а) и его начального давления (Ь). Из рис. 4 видно, что скорость фронта образования газогидрата H2S практически не зависит от начального давления пласта и его проницаемости. Это обусловлено тем, что скорость данного фронта лимитируется, прежде всего, отводом тепла, выделяющегося при образовании газогидрата H2S. Также, согласно рис. 4, скорость фронта вытеснения нефти сероводородом снижается с увеличением начального давления пласта и уменьшением его проницаемости. Это обусловлено тем, что скорость фронта вытеснения нефти сероводородом лимитируется скоростью течения в пласте, которая согласно закону Дарси пропорциональна перепаду давления в пласте и его проницаемости.

k0, м2 Ро, МПа

Рис. 4. Зависимость координат границ образования газогидрата H2S

Fig. 4. Coordinates of the boundaries of H2S gas hydrate formation (curve 1) and displacement of oil by hydrogen sulfide (curve 2) as functions of (a) permeability of the reservoir and (b) initial pressure

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

Рассматриваемый в данной работе режим реализуется в том случае, когда скорость фронта вытеснения нефти сероводородом (лимитируемая скоростью течения в пласте) превышает скорость фронта гидратообразования. Поэтому существует некоторое критическое значение давления инжекции сероводорода, ниже которого происходит слияние границ. На рис. 5 приведены зависимости величины данного давления от проницаемости при разных значениях исходного давления пластаР0 = 7.9 МПа (кривая 1) иР0 = 8 МПа (кривая 2).

Рис. 5. Зависимость предельного давления от проницаемости пласта при p0 = 7.9 МПа (кривая 1) иp0 = 8 МПа (кривая 2) Fig. 5. Limiting pressure as a function of permeability of the reservoir at po = (1) 7.9 and (2) 8 MPa

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

Заключение

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

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

ЛИТЕРАТУРА

1. Machel H.G. Geological and hydrogeological evaluation of the Nisku Q-Pool in Alberta, Canada, for H2S and/or CO2 storage // Oil Gas Sci. Technol. 2005. V. 60. P. 51-65. DOI: 10.2516/ogst:2005005.

2. Xu T., Apps J. A., Pruess K. and Yamamoto H. Numerical modeling of injection and mineral trapping of CO2 with H2S and SO2 in a sandstone formation // Chem. Geology. 2007. V. 24. No 3-4. P. 319-346. DOI: 10.1016/j.chemgeo.2007.03.022.

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

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

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

6. Шагапов В.Ш., Рафикова Г.Р., Хасанов М.К. К теории образования газогидрата в частично водонасыщенной пористой среде при нагнетании метана // Теплофизика высоких температур. 2016. Т. 54. № 6. С. 911-920. DOI: 10.1134/S0018151X16060171.

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

8. Хасанов М.К. Инжекция жидкого сероводорода в пласт, насыщенный нефтью и водой // Вестник Тюменского государственного университета. Физико-математическое моделирование. Нефть, газ, энергетика. 2017. Т. 3. № 2. С. 72-84. DOI: 10.21684/2411-79782017-3-2-72-84.

9. Хасанов М.К. Математическая модель инжекции сероводорода в пористую среду, насыщенную нефтью и водой // Изв. ТулГУ. Технические науки. 2017. Вып. 9. Ч. 1. С. 337-345.

10. 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 J. Heat and Mass Transfer. 2017. V. 107. P. 347-357. DOI: 10.1016/j.ijheatmasstransfer. 2016.11.034.

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

Khasanov M.K. (2018) MATHEMATICAL MODEL OF HYDROGEN SULFIDE INJECTION INTO A NATURAL RESERVOIR ACCOMPONIED BY ITS HYDRATE FORMATION Vestnik Tomskogo gosudarstvennogo universiteta. Matematika i mekhanika [Tomsk State University Journal of Mathematics and Mechanics]. 55. pp. 113-123

DOI 10.17223/19988621/55/10

Keywords: mathematical model, self-similar solution, porous medium, filtration, gas hydrates, hydrogen sulfide.

Based on the equations of mechanics of multiphase media, a mathematical model of liquid hydrogen sulfide injection into a porous formation saturated with oil and water is developed. To describe the processes of heat and mass transfer in the porous medium, the combined equations including mass and energy conservation laws, Darci's law, and the equation of state are used.

Self-similar solutions to the problem, in which three characteristic areas are formed in the reservoir, are obtained. In the first area (the nearest region), the porous are saturated with liquid hydrogen sulfide and its hydrate; in the second area (intermediate region), there is water and liquid hydrogen sulfide; and in the third area (distant region), the porous are saturated with oil and water.

Self-similar coordinates of the interphase borders are studied as functions of permeability of the reservoir and its initial pressure. It is established that increase in the initial pressure of the reservoir and decrease in its permeability cause a decrease in the velocity of oil displacement by hydrogen sulfide. It is shown that in conditions of low permeability of the reservoir and high initial pressure, the merge of interphase boundaries could happen. The dependences of critical values of injection pressure, corresponding to a merge of interphase boundaries, on the permeability of the reservoir and its initial pressure are found.

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

REFERENCES

1. Machel H.G. (2005) Geological and hydrogeological evaluation of the Nisku Q-Pool in Alberta, Canada, for H2S and/or CO2 storage. Oil Gas Sci. Technol. 60. pp. 51-65. DOI: 10.2516/ogst:2005005.

2. Xu T., Apps J.A., Pruess K., Yamamoto H. (2007) Numerical modeling of injection and mineral trapping of CO2 with H2S and SO2 in a sandstone formation. Chemical Geology. 24(3-4). pp. 319-346. DOI: 10.1016/j.chemgeo.2007.03.022.

3. Duchkov A.D., Sokolova L.S., Ayunov D.E., Permyakov M.E. (2009) Otsenka vozmozhnosti zakhoroneniya uglekislogo gaza v kriolitozone Zapadnoy Sibiri [Assesment of potential of West Siberian permafrost for the carbon dioxide storage]. Kriosfera Zemli - Earth's Cryos-phere. 13(4). pp. 62-68.

4. Chuvilin E.M., Guryeva O.M. (2009) Eksperimental'noe izuchenie obrazovaniya gidratov CO2 v porovom prostranstve promerzayushchikh i merzlykh porod [Experimental investigation of CO2 gas hydrate formation in porous media of frozen and freezing sediments. Krios-fera Zemli - Earth's Cryosphere. 13(3). pp. 70-79.

5. Byk S.Sh., Makogon Yu.F., Fomina V.I. (1980) Gazovye gidraty [Gas hydrates]. Moscow: Khimiya.

6. Shagapov V.Sh., Rafikova G.R., Khasanov M.K. (2016) On the theory of formation of gas hydrate in partially water-saturated porous medium when injecting methane. High Temperature. 54(6). pp. 858-866. DOI: 10.1134/S0018151X16060171.

7. Gimaltdinov I.K., Khasanov M.K. (2016) 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. 80(1). pp. 57-64. DOI: 10.1016/j.jappmathmech.2016.05.009.

8. Khasanov M.K. (2017) Inzhektsiya zhidkogo serovodoroda v plast, nasyshchennyy neft'yu i vodoy [Injection of liquid hydrogen sulfide in a layer saturated with oil and water]. Vestnik Tyumenskogo gosudarstvennogo universiteta. Fiziko-matematicheskoe modelirovanie. Neft', gaz, energetika - Tyumen State University Herald. Physical and Mathematical Modeling. Oil, Gas, Energy. 3(2). pp. 72-84. DOI: 10.21684/2411-7978-2017-3-2-72-84.

9. Khasanov M.K. (2017) Matematicheskaya model' inzhektsii serovodoroda v poristuyu sredu, nasyshchennuyu neft'yu i vodoy [Mathematical model of hydrogen sulfide injection into a porous medium saturated with oil and water]. Izvestiya Tul'skogo gosudarstvennogo universiteta - Izvestiya Tula State University. 9. pp. 337-345.

10. Shagapov V.Sh., Khasanov M.K., Musakaev N.G., Ngoc Hai Duong (2017) Theoretical research of the gas hydrate deposits development using the injection of carbon dioxide. International Journal of Heat and Mass Transfer. 107. pp. 347-357. DOI: 10.1016/ j.ijheatmasstransfer.2016.11.034.