Об устойчивости фронта вытеснения нефти в условиях термогазового воздействия на пласт

DOI: 10.21122/1029-7448-2016-59-2-141-150 УДК 622.01; 622.276.654

Об устойчивости фронта вытеснения нефти в условиях термогазового воздействия на пласт

К. В. Добрего1)

1)Белорусский национальный технический университет

© Белорусский национальный технический университет, 2016 Ве1агш1ап National Technical University, 2016

Реферат. Термогазовый метод воздействия на пласт является перспективным методом повышения нефтеотдачи и в настоящее время осваивается в Беларуси. Устойчивость температурного фронта и фронта вытеснения нефти - принципиальный вопрос данной технологии. В статье методом малых возмущений решается задача об устойчивости фронта вытеснения нефти при термогазовом воздействии на пласт. Показано, что инкремент роста возмущений отличается от случая фильтрационного горения газа и определяется волновым числом волны возмущения, разностью скоростей фильтрации дутья и фронта тепловыделения, отношением плотностей, коэффициентов фильтрации, сжимаемостей среды дутья и вытесняемой нефти и другими величинами. Проанализированы основные параметры, влияющие на устойчивость фронта. Даны рекомендации по методам и процедурам, повышающим устойчивость (компенсирующим неустойчивость) фронта вытеснения. Механизмами подавления или компенсации возмущений фронта являются: снижение вязкости нефтяной смеси и повышение вязкости дутья; увеличение сжимаемости дутья; увеличение толщины теплового фронта термогазового воздействия прежде всего за счет повышения его температуры и изменения химизма (использование низкотемпературных катализаторов); снижение скорости волны и/или соответственно подачи дутья. Применение воздушного или водовоз-душного дутья, а также насыщение воды плохо растворимыми газами обеспечивают относительно большую устойчивость фронта по сравнению с использованием составов на основе воды. Результаты и рекомендации предлагается использовать в рамках общих методик и процедур управления процессом термогазового воздействия в условиях опытной и опытно-промышленной нефтедобычи. Дальнейшие исследования методов и механизмов компенсации неустойчивости фронта, в том числе его динамики на нелинейной стадии, могут быть проведены путем детального многомерного расчета системы.

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

Для цитирования: Добрего, К. В. Об устойчивости фронта вытеснения нефти в условиях термогазового воздействия на пласт / К. В. Добрего // Энергетика. Изв. высш. учеб. заведений и энерг. объединений СНГ. 2016. Т. 59, № 2. C. 141-150

On Stability of the Oil Displacement Front

under Conditions of Thermogas Impact on the Oil Bearing Layer

K. V. Dobrego1

1)Be1arusian National Technical University

Abstract. The method of thermogas impact on the oil bearing layer is a perspective one to improve oil recovery. It is now under experimental implementation in Belarus. Stability of

Адрес для переписки Address for correspondence

Добрего Кирилл Викторович Dobrego Kirill V.

Белорусский национальный технический университет Be1arusian National Technical University

просп. Независимости, 65/2, 65/2 Nezavisimosty Ave.,

220013, г. Минск, Республика Беларусь 220013, Minsk, Republic of Belarus

Тел.: +375 17 292-42-32 Tel.: +375 17 292-42-32

ef@bntu.by ef@bntu.by

the thermal front and that of the oil displacement is a point of principle for the given technology. The article offers solution based on the method of small perturbations for the problem of the oil-displacement front stability under thermogas impact on the oil layer. The perturbation growth increment is shown to be quite different from that of the gas filtration combustion. Being determined by the perturbation wave number (or wave length), by difference in speed of the blowing filtration and that of the heat development front, by the ratios of densities, filtration coefficients, compressibility of the displacement agent and the displaced oil mass, and by other parameters. The author analyses the main parameters influencing the stability of the front. Recommendations are given on possible methods and procedures improving stability of the displacement front (neutralizing perturbation growth). The mechanisms of suppression or compensation of the front perturbation growth are as follows: the crude oil viscosity reduction and increase of that of the displacement agent, the displacement agent compressibility increase, increase of the thermogas impact heat-front width in the first instance at the expense of the temperature growth and alteration in the chemism (low temperature catalysts utilization, etc.), reduction in speed of the front propagation and/or correspondingly of the displacement agent feed. Utilization of the gas or water-gas displacement agent as well as the agent impregnation with hydrophobic gases provides a relatively better stability of the front as compared to that employing water based compounds. Analytical data and derived recommendations are offered for utilization within the framework of general techniques and procedures of the thermogas impact process management in the context of developmental and experimental-industrial oil extraction. Further research into methods and mechanisms of the displacement front perturbation compensation, including the front dynamics at nonlinear stage, might be conducted via detailed 2D and 3D numerical simulation of the system.

Keywords: thermogas method, oil recovery, stability, displacement front

For citation: Dobrego K. V.

Введение

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

Экспериментальные исследования, проведенные П. Г. Саффманом, Г. И. Тейлором, Р. Л. Чоуком [1-3] и другими, показали, что развитие возмущений плоского фронта вытеснения в пористой среде при нарушении устойчивости происходит в виде неограниченно разрастающихся языков вытесняющей жидкости. Эксперименты на насыпных пористых средах [4] свидетельствуют о том, что нарушение устойчивости происходит при отношении вязкости взаимодействующих жидкостей, превышающем критическое значение 10-15. В то же время при малых скоростях вытеснения возмущения затухают даже при отношениях вязкостей, больших критического значения. Математическая модель развития языков вытесняющей жидкости предложена Г. И. Баренблаттом [5]. Влиянию капиллярных сил на процесс двухфазной фильтрации посвящены работы В. И. Пеньковского [6, 7]. Влияние вибраций на устойчивость плоского фронта вытеснения в пористой среде рассмотрено Д. В. Любимовым [8]. Обобщение модели двух-

фазной фильтрации Рапопорта - Лиса на случай анизотропных пористых сред представлено в работе М. Н. Дмитриева [9]. При этом устойчивость фронтов вытеснения, распространяющихся совместно с тепловыми и концентрационными полями, исследована значительно хуже.

Устойчивость фронта вытеснения нефти при ТГМ может быть описана и исследована традиционными методами, апробированными в теории горения и гидродинамике. В данной статье воспользуемся методом малых возмущений аналогично исследованию устойчивости фронта при фильтрационном горении газов (ФГГ) [10], учитывая физические особенности, отличающие процесс вытеснения нефти при ТГМ.

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

Рассмотрим плоский фронт вытеснения, сформированный равномерно подаваемым дутьем при квазистационарном горении внутри пласта (рис. 1).

Т 4

hi h /ii

Рис. 1. Схема плоского фронта вытеснения нефти при ТГМ: /l - фронт тепловыделения; fii - среднемассовый фронт вытеснения нефтяной смеси; /2i - передний фронт соприкосновения остаточных углеводородов и дутья

Fig. 1. Diagram of the plane front of oil displacement at TGM: /l - emission of heat front; /li - average weight displacement front of the crude oil; /2i - front contact edge of the residual hydrocarbons and the displacement agent

Ввиду того что закачиваемый газ при высоком давлении и температуре может находиться в сверхкритическом состоянии, не будем употреблять термин «газ», а используем термин «флюид», относя его как к дутью, так и к вытесняемой смеси - нефти с водой дутьем и продуктами горения. В отсутствие горения среднемассовый фронт вытеснения будет распространяться со скоростью движения флюида дутья Vi, выполняющего роль поршня, а сам фронт размываться со временем за счет как диффузии флюидов, так и механической дисперсии Ddis (задержание и перемешивание нефти в погранслоях и застойных областях пор). Скорость фронта соприкосновения /21 может быть связана со скоростью фильтрации дутья u / 2i = const • Vi.

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

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

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

РАД^ - и/1] = Р2^и/1

^ТГМ

где С0 - концентрация кислорода в дутье; ^ - массовый стехиометриче-ский коэффициент для горючего вещества; 10 - характерный линейный масштаб дисперсии и диффузии, относительно которого оценивается уменьшение содержания нефти в слое вытеснения 10/7 ТГМ.

Особенностью данной задачи является разделение фронта тепловыделения (области максимального тепловыделения или максимального градиента температурного фронта) и среднемассового фронта вытеснения. Для гидродинамической задачи границей раздела сред будем считать сред-немассовый фронт вытеснения. Примем, что на границе имеет место резкое изменение параметров среды и фильтрации - фронты, разделяющие относительно однородные флюиды, узки по сравнению со всеми характерными размерами пласта. Слева и справа от границы раздела флюиды характеризуются скоростями движения У1 и У2, плотностями р1 и р2, температурами Т1 и Т2, вязкостью V и у2, коэффициентами сжимаемости в 1 и в 2 (рис. 1). Рассматривается двумерная система, где вертикальная координата пласта не учитывается.

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

и% = -кУр;

^ + ^у

от

где к - коэффициент фильтрации, к = £0/(тц). Предполагается, что вязкость не зависит от давления.

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

Р = Ро

где р0, р0 - плотность и давление при фиксированных

условиях; в = 1/р0(5р/ф) - коэффициент сжимаемости.

Коэффициент сжимаемости нефти зависит от температуры, количества растворенного в нефти газа, фракционного состава. Для дегазированной нефти он колеблется от 4 • 10-10 до 7 • 10-10 Па-1, а для легких нефтей, в которых растворено значительное количество газа, достигает 140 • 10-10 Па-1. С повышением температуры сжимаемость растет. Для большинства месторождений нефти сжимаемость не выходит за пределы (7-30) • 10-10 Па-1 [11] и для рассматриваемых оценок принята в = 10-9 Па-1.

Подставляя (2) и (4) в (3) и принимая во внимание, что вр0 << 1, приходим к уравнению нестационарной фильтрации

дР = КД

& в

Связь между деформацией фронта горения и возмущением поля фильтрации устанавливается исходя из известной зависимости между скоростью волны ФГГ и скоростью фильтрации флюида [10]

им. = им. ).

В отличие от случая газовой фильтрации, уравнение (5) - линейное. Следовательно, суперпозиция его решений также является решением, и возмущение должно удовлетворять (5) независимо от вида основного невозмущенного решения. Для определенности считаем, что невозмущенное решение характеризуется скоростью движения фронта ык, давлением р = = С1 + хС2 (стационарное решение (5)), нормальной к фронту скоростью фильтрации ¥1 в области перед фронтом и У2 = (р1 / р2)^ - за фронтом.

Найдем функцию возмущения давления, удовлетворяющую (5) и отвечающую требованиям обращения в нуль на бесконечности и периодичности по поперечной координате:

р ~ехр

где О - временной инкремент роста возмущения О = (к / в)(к2х — к у).

Согласно (7), для квазистационарных процессов кх1 » ку1 = ку2 » кх 2,

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

р1 = А ехр^у у + kх х + ^);

V =-к1kх А ехр^у у + kx х + 00;

и1 = Aexp(iky у + kx х + Qt);

р2 = 5 ехр(/куу - кхх + О);

у2 = к 2кхВ ехр(/куу - кхх + О);

и2 = -/'к2куВ ехр(/куу - кхх + О);

/ = С ехр

где V, и - нормальная и касательная к начальному фронту составляющие возмущения скорости фильтрации; / - функция смещения (деформации) фронта вытеснения/горения.

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

«

А П Л / Г, \

к2в1 к1в

и1 +

2 У

к2в1 к1в

2 У

Для записи третьего условия используем уравнение (6) в предположении, что у1 << V1 :

/;=^,

где а = р1/р2; и № - скорость невозмущенного фронта волны термогазового

воздействия; в = ди№ /ди - отклик фронта горения на изменение скорости

фильтрации (определяемый тепломассообменом во фронте горения).

Заметим, что уравнение (10) не следует из условия сохранения потока импульса или галилеевой инвариантности явления, подобно случаю свободного пламени, и не соответствует принятому в [12]. Его можно получить исходя из (2) и непрерывности давления при фильтрации. Действи-

тельно, из непрерывности давления на границе фронта горения (непрерывность изобар на рис. 2) и из геометрических соображений следует [Ур^т а1 = |Ур2|эт а2. Тогда тангенциальная компонента скорости до

разрыва и1 = (к1/р1)[Ур^та1, а после разрыва и2 = (к2/р2)|Ур2|эта2 = = (к2Р1/к1Р2)и1, т. е. сохраняется с точностью до множителя к 2Р1/к1Р2. Это поясняет рис. 2.

Рис. 2. Пояснение к выводу (10): V1, V2, а1, а2 - скорости и углы отклонения скорости фильтрации от нормали до фронта вытеснения и после него соответственно; штриховые линии - изобары

Fig. 2. Explanation for the conclusion (10): V1, V2, ab a2 - speeds and angles of the filtration speed deviation from the normal before the displacement front and after it respectively; dashed lines - isobar

a2

Пренебрегая величинами возмущений скорости (но не скоростью движения фронта) по сравнению с V, а также учитывая, что при малых амплитудах деформации ехр(кхх) « 1, подставляем (8) в уравнения (9)-(11):

aQC + a

f \ к

Pi

kA = QC -

v^i

f \ K2

V^ 2

kB;

К2в1 V К1в2 J

[V - uw ]kC -

К2в1 VK1P2 J

f \ К1

vP1 J

kA = a[V - uw]kC -

/ л к.

kB;

VH 2 J

QC = -e

f \ к

P

kA.

V^ 1 J

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

a

Q = ek[V -uw]-

( п ^

K2P1

V К1в2 J

a +

f n > '

K2P1

V К1в2 J

Для случаев ФГГ и ТГМ значения параметров в (15) существенно различаются. Если в первом варианте а = р1/р2 = Т2/Т1 >> 1, то во вто-

ром а = р/р2~1. Величина к2/к1 ~ \]Т1/Т2 < 1 в первом случае и также

меньше единицы во втором из-за более высокой вязкости нефти. Таким образом, формула (15) приводит к различным выводам для упомянутых случаев - абсолютной неустойчивости фронта ФГГ с О ~ ки№ и условной неустойчивости фронта ТГМ при

— > 0 или Р1К1Р2 >Р2К2Р1.

Р2 К1в2

Учитывая определение коэффициента фильтрации к = к0/(тц) и кинематической вязкости у = ц / р, условие неустойчивости (16) сведется к требованию у2Р2 > у1Р1. Таким образом, предпосылкой для неустойчивости фронта вытеснения при ТГМ являются высокая кинематическая вязкость и сжимаемость нефтяного флюида по сравнению с дутьем. Отсюда следует, что вытеснение несжимаемым водным раствором способствует неустойчивости, а использование газового или водогазового дутья - стабилизации фронта вытеснения. Из данного анализа следует, что при прогреве нефтяного флюида и образовании теплового фронта устойчивость фронта вытеснения повышается.

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

К = ^,0 Л".

Коэффициент у оценивается через толщину теплового фронта Ах и безразмерную скорость его движения

у = 2 Ах

г \

и,

V и у

В случае ТГМ [3] ширина фронта может составлять несколько метров в силу низкой максимальной температуры, а безразмерная скорость может примерно соответствовать скорости тепловой волны и№/и{ ~ 1. Инкремент роста возмущения имеет вид

а

^ о >

К2Р1

V К1в2 У

О = гк[Ух - п№ ]-р^ - ущк2.

к2р,

а +

V к1в2 У

Оценивая порядок уи{к2, можно показать, что неоднородности с характерными длинами волн около нескольких толщин фронта будут подавлять-

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

ВЫВОДЫ

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

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

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

ЛИТЕРАТУРА

1. Saffman, P. G. The Penetration of a Fluid Into a Porous Medium or Hele-Shaw Cell Containing a More Viscous Liquid / P. G. Saffman, G. I. Taylor // Proc. R. Soc. London, 1958. A245, Р. 321-329.

2. Saffman, P. G. Viscous Fingering in Hele-Shaw Cells / P. G. Saffman // J. Fluid Mech. 1986. No 173. Р. 73-94.

3. Chouke, R. L. Instability of Slow, Immiscible, Viscous Liquid-Liquid Displacements in Permeable Media / R. L. Chouke, P. Meurs, C. Poel // Pet. Trans. AIME. 216. 1959. P. 188-194.

4. Рыжик, В. М. Исследование устойчивости движения раздела воды и нефти / В. М. Рыжик, Б. Е. Кисленко // ВНК. Физико-геологические факторы при разработке нефтяных и нефтегазоконденсатных месторождений. М.: Недра, 1969. С. 82-92.

5. Баренблатт, Г. И. Движение жидкостей и газов в природных пластах / Г. И. Баренблатт, В. М. Ентов, В. М. Рыжик. М.: Недра, 1984. 207 с.

6. Разработка целиков нефти при воздействии на пласт химических реагентов / В. И. Пень-ковский [и др.] // Прикладная механика и техническая физика. 2013. Т. 54, № 3. С. 87-94.

7. Остаточные нефтенасыщенные зоны продуктивных пластов и способы воздействия на них с целью вовлечения в разработку / В. И. Пеньковский [и др.] // Физико-технические проблемы разработки полезных ископаемых. 2012. № 5. С. 41-51.

8. Любимов, Д. В. Влияние вибраций на устойчивость плоского фронта вытеснения в пористой среде / Д. В. Любимов, Г. А. Седельников // Известия РАН. Механика жидкости и газа. 2006. Т. 41, № 1. С. 6-14.

9. Дмитриев, М. Н. Модель двухфазной фильтрации Рапопорта - Лиса в анизотропных средах / М. Н. Дмитриев // Известия РАН. Механика жидкости и газа. 2011. Т. 46, № 2. С. 136-144.

10. Добрего, К. В. Физика фильтрационного горения газов / К. В. Добрего, С. А. Жданок. Минск: Институт тепло- и массообмена НАН Беларуси, 2002. 204 с.

11. Щелкачев, В. Н. Основы и приложения теории неустановившейся фильтрации / В. Н. Щел-качев. М.: Нефть и газ, 1995. 493 с.

12. Минаев, С. С. О неустойчивости фронта пламени при фильтрационном горении газов / С. С. Минаев, С. И. Потытняков, В. С. Бабакин // Физика горения и взрыва. 1994. Т. 30, № 6. С. 49-54.

Поступила 11.01.2016 Подписана в печать 12.03.2016 Опубликована онлайн 30.03.2016

REFERENCES

1. Saffman P. G., Taylor G. I. (1958) The Penetration of a Fluid Into a Porous Medium or Hele-Shaw Cell Сontaining a More Viscous Liquid. Proceedings of the Royal Society of London, A245, 321-329.

2. Saffman P. G. (1986) Viscous Fingering in Hele-Shaw Cells. Journal of Fluid Mechanics, 173, 73-94.

3. Chouke R. L., Meurs P., Poel C. (1959) Instability of Slow, Immiscible, Viscous LiquidLiquid Displacements in Permeable Media. Trans. AIME, 216, 188-194.

4. Ryzhik V. М., Kislenko B. Ye. (1969) Investigation of Dynamic Stability of Oil-Water Interface. Physical-Geological Factors at Development of Oil and Oil-Gas Condensate Fields. Moscow, Nedra, 82-92 (in Russian).

5. Barenblatt G. I., Yentov V. М., Ryzhik V. М. (1984) Liquid and Gas Motion in Natural Fo-mations. Moscow, Nedra. 207 (in Russian).

6. Penkovskii V. I., Korsakova N. K., Altunina L. K., Kuvshinov V. A. (2013) Development of Unrecovered Oil Affecting the Formation by Chemical Reagents. Journal of Applied Mechanics and Technical Physics, 54 (3), 415-422 DOI: 10.1134/S0021894413030103.

7. Penkovskii V. I., Korsakovа N. K.. Simonov B. F., Savchenko A. V. (2012) Residual Oil-Saturated Zones of Productive Strata and Methods of Affecting them Towards Involvement Into Development. Fiziko-Tekhnicheskie Problemy Razrabotki Poleznykh Iskopaemykh [Phy-sicotechnical Challenges of Mining], (5), 41-51 (in Russian).

8. Lubimov D. V., Sedelnikov G. А. (2006) Effect of Vibration on the Stability of a Plane Displacement Front in a Porous Medium. Fluid Dynamics, 41 (1), 3-11. DOI: 10.1007/s10697-006-0016-0.

9. Dmitriev М. N. (2011) Rapoport - Leas Two-Phase Flow Model for Anisotropic Porous Media. Fluid Dynamics, 46 (2), 291-298. DOI: 10.1134/S0015462811020128.

10. Dobrego K. V., Zhdanok S. А. (2002) Physics of Filtrational Combustion of Gases. Minsk: Heat and Mass Transfer Institute of the NAS of Belarus. 204 (in Russian).

11. Shchelkachev V. N. (1995) Basis and Applications of the Unsteady Flow Theory. Moscow, Neft i gaz. 493 (in Russian).

12. Minaev S. S., Potytnyakov S. I., Babakin V. S. (1994) On Instability of Flame Front During Filtrational Combustion of Gases. Vizika Goreniya i Vzryva [Combustion, Explosion, and Shock Waves], 30 (6), 761-763. DOI: 10.1007/BF00755247.

Recеived: 11 January 2016 Accepted: 12 March 2016 Published online: 30 March 2016