автореферат диссертации по информатике, вычислительной технике и управлению, 05.13.18, диссертация на тему:Математическое моделирование изотермической многокомпонентной многофазной фильтрации с фазовыми переходами

Введение.

1. Математические модели многокомпонентной фильтрации.

1.1 Уравнения многокомпонентной фильтрации.

1.1.1 Уравнения двухкомпонентной фильтрации.

1.2 Модель Баклея-Леверетта.

1.3 Некоторые математические свойства моделей фильтрации.

1.4 Задача Баклея-Леверетта.

2.Термодинамика многокомпонентных систем и уравнения состояния.

2.1 Фазовое равновесие.

2.2 Уравнения состояния.

2.3 Соотношения, связывающие потенциал Гиббса и кривые фазового равновесия.

2.4 Соотношения, связывающие уравнения состояния и кривые фазового равновесия.

2.5 Модельные уравнения состояния и потенциал Гиббса.

3. Разрывные решения уравнений многокомпонентной фильтрации.

3.1 Баланс свободной энергии.

3.2 Термодинамическое условие на разрывах.

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

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

4.1 Вычислительные алгоритмы для уравнений фильтрации несжимаемых флюидов.

4.2 Вычислительные алгоритмы для уравнений фильтрации сжимаемых флюидов.

4.3 Тестирование вычислительных алгоритмов для уравнений фильтрации сжимаемых флюидов с фазовыми переходами.

Настоящая работа посвящена математическому моделированию изотермической многокомпонентной многофазной фильтрации с фазовыми переходами. Исследуется модель, которая широко используется для прогнозирования разработки нефте- и газосодержащих пластов [1, 9, 41, 50, 51, 58, 21, 70, 69, 76, 84, 85]. Однако, фильтрационные течения растворов, сопровождающиеся фазовыми переходами (особенно с изменением числа фаз), изучены недостаточно.

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

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

2) установить условия термодинамического согласования модели, т.е. критерии непротиворечивости уравнений состояния фаз и кривых фазового равновесия;

3) построить термодинамически согласованные модели растворов, которые были бы достаточно просты для эффективного численного моделирования сложных фильтрационных течений и в тоже время достаточно точно передавали бы фазовое поведение растворов в некотором диапазоне давлений и концентраций;

4) разработать методы отбора физически недопустимых разрывных решений;

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

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

Модели многокомпонентной многофазной фильтрации изучаются, в основном, в связи с проблемами повышения эффективности добычи природных углеводородов — нефти, газа, газоконденсата. Модель, описывающая двухфазную фильтрацию несжимаемых флюидов без фазовых переходов, была предложена в 1941 г. Баклеем и Левереттом [71] и до сих пор широко используется для моделирования задач вытеснения [19, 20, 51, 41, 70, 76, 50, 84, 85, 88]. В рамках модели Баклея-Леверетта, система уравнений двухфазной фильтрации расщепляется на гиперболическое уравнение для фазовой насыщенности и эллиптическое уравнение для давления.

Свойства системы уравнений, описывающих многокомпонентную многофазную фильтрацию, также хорошо изучены в предположениях, что все фильтрующиеся фазы несжимаемы и для всех компонентов при смешении выполнено правило Амаго (аддитивность парциальных объемов) [70]. В этих предположениях из уравнений, выражающих законы сохранения количеств каждого компонента, вытекает условие бездивергентности поля скорости фильтрации. В одномерной геометрии указанное свойство означает, что полный (объемный) поток фильтрующихся флюидов одинаков через все сечения, хотя, вообще говоря, может зависеть от времени. В связи с этим говорят, что такие модели описывают течения с постоянным полным потоком. При этом уравнения многокомпонентной фильтрации расщепляются на эллиптическое уравнение для давления и систему уравнений первого порядка для концентраций компонентов, в которую давление не входит. Если число компонентов невелико, удается сформулировать условия гиперболичности последней системы и провести анализ допустимости разрывных решений, исходя из условий эво-люционности и/или существования структуры скачка [7, 17]. Однако, правило Амаго выполняется далеко не всегда. Условие несжимаемости фаз также не всегда является удачным приближением. Например, для "газированной нефти" закон Амаго не выполняется, а выделяющийся газ нельзя считать несжимаемым [3, 8, 65].

Нефть, конденсат, природные газы — это залегающие в недрах Земли многокомпонентные растворы (смеси), состоящие, в основном, из углеводородных соединений с добавкой неуглеводородных веществ. Сложное поведение таких смесей, когда нельзя рассматривать каждый компонент в отдельности, описывает физика растворов [4, 40, 47, 63, 64, 66]. Растворами называются физически однородные (гомогенные) смеси двух или нескольких веществ. Из-за взаимодействия молекул при растворении происходят реакции, которые обычно сопровождаются выделением или поглощением тепла, что делает растворы похожими на химические соединения. От химических соединений растворы отличаются тем, что в химические соединения вещества вступают в строго определённых пропорция, тогда как относительные количества веществ в растворах могут меняться в более или менее широких пределах. В зависимости от термобарических условий и от величины полной молярной концентрации раствор находится в однофазном состоянии или расслаивается на несколько фаз. Фазы разделены поверхностями раздела. Сложное фазовое поведение отражают фазовые диаграммы, которые определяют не только границы фазовых областей, но и задают концентрации в фазах.

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

Для моделирования фильтрационных течений таких смесей разработана и применяется иерархия моделей разной степени сложности и подробности. Эти модели описывают фильтрацию многокомпонентных растворов в тех или иных предположениях о характере течения (перепаде давления, диапазоне изменения состава раствора и т.п), термодинамических свойствах фаз, влиянии капиллярных эффектов, свойствах коллектора и т.д. Ясно, что чем подробнее модель, тем сложнее она поддается теоретическому анализу и исследованию с привлечением методов численного моделирования. В настоящее время наиболее исследованными являются модели, объединенные условием сохранения полного потока фильтрующегося флюида. Такие модели описывают фильтрацию квазиидеальных растворов, т.е. для каждой фазы выполнен закон Амаго (парциальные объемы компонентов не зависят от давления). Предельным случаем таких моделей являются модели типа Баклея-Леверетта, описывающие фильтрацию несжимаемых несмешивающихся флюидов [5].

В настоящей работе рассматривается модель многокомпонентной фильтрации, в которой относительно свойств парциальных объемов компонентов никаких предположений не делается за исключением выполнения неравенств, гарантирующих термодинамическую устойчивость системы. Однако, принята гипотеза, согласно которой все компоненты в том или ином количестве присутствуют во всех фазах [47, 56]. В этом случае условия фазового равновесия выражаются равенством химпотен-циалов компонентов во всех фазах. Отказ от этой гипотезы означал бы, что в некоторой фазе не растворен какой-либо компонент смеси, как это принимается в модели "газированной нефти" [1]. В этом случае химпо-тенциал нерастворимого в данной фазе компонента не равен химпотен-циалу этого компонента в других фазах. Однако, это обстоятельство не влияет ни на справедливость правила фаз Гиббса, ни на полученные в работе результаты: баланс свободной энергии, условия допустимости разрывов уравнений многокомпонентной фильтрации. В рассматриваемой модели также полагается, что течение является изотермическим и достаточно медленным, чтобы успело установиться локальное фазовое равновесие. Кроме того, считается, что капиллярным скачком давления между фазами можно пренебречь и для скорости фильтрации каждой фазы выполнен закон Дарси.

Изучению свойств растворов, в частности углеводородных, — расслоение на фазы, уравнения состояния фаз, изменение состава фаз при изменении термобарических условий, ограниченная растворимость компонентов и т.д. — посвящено множество работ [4, 40, 47, 64, 66], тем не менее термодинамическое описание сложных многокомпонентных растворов сталкивается с большими трудностями, как теоретического, так и технического характера. Исследования показали, что в гомогенной области раствор ведёт себя как гипотетическое чистое вещество с некоторыми псевдопараметрами. Каждая из фаз многокомпонентного раствора описывается своим уравнением состояния. Разработано много методик, каждая из которых с той или иной точностью описывает раствор, в том или ином диапазоне давлений и температур [4, 43, 44, 64]. Например, широко используется подход, позволяющий по единому уравнению состояния (Ван-дер-Ваальса, Пенга-Робинсона и др.) вычислять термодинамические потенциалы и определять кривые фазового равновесия (конденсации и кипения). Однако, для реального нефтегазового раствора, состоящего из десятков компонент, находящегося под давлением сотни атмосфер, рассчитать коэффициенты единого уравнения часто проблематично. Поэтому в некоторых диапазонах изменения параметров смеси такой подход дает значительные ошибки.

Некоторые авторы полагают [53], что в настоящее время возможно самый перспективный путь — создавать для каждого месторождения свое уравнение состояния, которое бы содержало информацию не только о флюиде, но и о коллекторе (т.е. учитывало влияние пористости и материала коллектора на уравнение состояния).

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

Один из вопросов, изучаемый в настоящей работе, состоит в следующем: какую информацию о свойствах модели изотермической многокомпонентной фильтрации, которая будет описана в главе I, можно получить, опираясь на общие термодинамические соотношения, имеющее место для гетерогенных систем. Так при разработке вычислительных алгоритмов важны "гиперболические" качества системы уравнений. Надежду на результативность такой постановки вопроса дает теория С.К.Годунова, развитая для систем первого порядка, обладающих дополнительным балансным соотношением [15]. Применительно к рассматриваемой здесь модели, выражающей законы сохранения количеств каждого из компонентов, участвующих в фильтрации, таким дополнительным соотношением является уравнение баланса свободной энергии. В отличие от теории С.К.Годунова уравнения сохранения для компонентов не образуют систему первого порядка. В частности по этой причине система не является гиперболической. Тем не менее она обладает некоторыми "гиперболическими" качествами, которые проявляются в конечной скорости распространения возмущений концентрации компонентов и фазовых насыщенностей. Возмущение же давления распространяются мгновенно. Важной задачей теории является изучение разрывных решений уравнений изотермической многокомпонентной фильтрации. В рамках рассматриваемой здесь модели сильные разрывы трактуются как поверхности, на которых претерпевают скачки концентрации компонентов раствора и объемные насыщенности фаз, в то время как давление (одинаковое во всех фазах) остается непрерывным. При этом значения гидродинамических и термодинамических величин по обе стороны от разрывов связаны соотношениями Гюгонио. На самом деле, сильные разрывы концентраций компонентов раствора и фазовых насыщенностей — суть узкие зоны неравновесных фазовых переходов. Процессы, протекающие в пределах этих переходных зон и определяющие их внутреннюю структуру, здесь не рассматриваются. В качестве критерия допустимости разрыва принимается условие невозрастания на нем свободной энергии смеси.

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

Как известно из теории систем гиперболических уравнений, не все скачки, удовлетворяющие условиям Гюгонио, являются физически допустимыми, часть из них должна быть запрещена. Для этого необходимо использовать дополнительные условия. Таких условий в газовой динамике может быть сформулировано, по крайней мере,три [14, 45, 42, 59]:

1) эволюционность скачка;

2) существование структуры скачка;

3) энтропийное условие.

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

Эволюционность скачка. В моделях многокомпонентной фильтрации с постоянным полным потоком широко используется условие "априорной эволюционности скачка" или "условие Лакса" [45, 82], которое связывает количество приходящих на разрыв и уходящих с разрыва характеристик.

Существование структуры скачка. В работах [7, 17, 70] для модели двухфазной фильтрации несмешивающихся флюидов с межфазным обменом активной примесью условия существования структуры скачка используются для выделения допустимых разрывов решений. Наряду с учетом капиллярного скачка давлений между фазами, приводящего к появлению диффузионных слагаемых в уравнениях переноса компонентов, в работе [7] учитывается также неравновесная адсорбция. Энтропийное условие. В подземной гидродинамике оно не было сформулировано и поэтому ранее не использовалось. Это условие не зависит от типа уравнений, описывающих систему, и выражает термодинамический закон — неубывание энтропии: неравновесные процессы могут увеличивать (создавать) энтропию, но не могут уменьшать (уничтожать её). Как отмечено в работе [56]: "Основной проблемой, стоящей перед термодинамикой неравновесных процессов, является точное вычисление возрастания энтропии". Отсутствие таких точных выражений приводит к тому, что энтропийное условие обычно формулируется в виде неравенств, что несколько снижает эффективность отбора нефизических решений. Для изотермических задач вместо энтропийного условия рассматривается условие на свободную энергию: неравновесные процессы могут уменьшать свободную энергию, но не могут увеличивать её [47, 56].

Хотя уравнения многокомпонентной фильтрации сжимаемых флюидов не попадают под классификацию на эллиптические, параболические или гиперболические, они демонстрируют некоторые свойства этих типов уравнений. Поэтому в настоящей работе был принят подход, в рамках которого вычислительные алгоритмы для численного интегрирования этих уравнений, строились по аналогии с алгоритмами, разработанными для решения классических типов уравнений. Однако, как отмечено авторами монографии [45]: "прямое применение даже хорошо проверенных численных методов для моделирования недостаточно понятных физических проблем может привести к неверным результатам". В этом случае особое значение приобретает тестирование вычислительных алгоритмов, которое, по крайней мере в принципе, позволяет выявить их недостатки и установить область их надежного применения. Тестирование вычислительных алгоритмов естественно проводить на относительно простых задачах, решения которых тем не менее содержат особенности, характерные для решений общего случая. К таким особенностям относится присутствие в решении сильных и слабых разрывов концентраций компонентов и фазовых насьпценностей. Появление разрыва решения сопровождаются потерей локальной аппроксимации, что приводит к необходимости использования консервативных разностных схем. Другой особенностью уравнений многокомпонентной фильтрации является невыпуклость функции потока и соответственно распространение разрывов в режиме Жуге. Эта особенность рассматриваемого класса моделей содержится уже в уравнении Баклея-Леверетта. Поэтому вычислительный алгоритм, разработанный для общего случая, должен правильно описывать фильтрацию несжимаемых несмешивающихся флюидов в модели Баклея-Леверетта, когда система расщепляется на эллиптическое уравнение для давления и гиперболическое уравнение для насыщенности.

Классической задачей теории двухфазной фильтрации является задача Баклея-Леверетта о вытеснении одного несжимаемого флюида другим. Эта задача допускает решение в квадратурах, хотя оно и не является автомодельным в обычном смысле. Для построения этого решения важно, что модель Баклея-Леверетта относится к классу моделей с постоянным полным потоком и расщепляется на эллиптическое уравнение для давления и гиперболическое уравнение для насыщенности. Указанное решение для насышенности становится автомодельным после соответствующей перенормировки времени и может быть построено стандартными методами, а давление затем находится из условия постоянства полного потока. Этот прием расщепления системы на эллиптическую и гиперболическую подсистемы широко используется для построения и анализа решений уравнений многокомпонентной фильтрации в моделях с постоянным полным потоком. В рамках такого подхода исследовались процессы вытеснения несжимаемых флюидов, были построены автомодельные (в указанном выше смысле) решения уравнений многокомпонентной фильтрации о распаде разрыва [11, 8, 2, 49]. За пределами этого подхода находится работа [8], в которой построены и проанализированы автомодельные решения задачи взаимного вытеснения двух несмешива-ющихся фаз с учетом сжимаемости одной из них.

Основные принципы, положенные в основу построения вычислительного алгоритма, состояли в следующем:

1) консервативность разностной схемы для каждого компонента раствора;

2) второй порядок пространственной аппроксимации в областях гладкости решения;

3) аппроксимация доли компонент в потоке против скорости фильтрации;

4) аппроксимация градиента давления на неявном временном слое.

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

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

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

Актуальность.

Актуальность темы диссертации обусловлена тем обстоятельством, что фильтрационные течения растворов, сопровождающиеся фазовыми переходами (особенно с изменением числа фаз), изучены недостаточно. В то же время, даже небольшое выпадение новой фазы в какой-то области пласта изменяет режим течения, может резко увеличить или уменьшить гидродинамическое сопротивление, вызвать колебания на скважине [53, 55, 86]. Выделение пузырьков газа из добывамой нефти в мелких порах иногда приводит к полной закупорке коллектора [62]. Непредсказуемые изменения состава, возникновение колебаний на скважине — это типичные проблемы разработки залежей "газированной" нефти и газоконденсата. Вследствие этого в пласте безвозвратно теряются значительные запасы углеводородов. На Вуктыльском конденсатном месторождении 70% конденсата осталось в пласте. При понижении давления на скважинах Уренгойского газоконденсатного месторождения, через некоторое время после начала эксплуатации наблюдалось значительное облегчение и осветление добываемого конденсата и т.д.

Цель работы.

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

Научная новизна.

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

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

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

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

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

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

Практическая ценность

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

Апробация работы.

Основные результаты докладывались автором на IX Всероссийском совещании по проблемам построения сеток для решения задач математической физики, посвященной памяти академика А.Ф. Сидорова, 2002; на XIV Всероссийской конференции "Теоретические основы и конструирование численных алгоритмов для решения задач математической физики", посвященной памяти К.И.Бабенко, 2002; на конференции "Ломоносовские чтения", МГУ, 2003; in International conference "New Trends in Continuum mechanics", Constatza, Ruminia, 2003; на XI школе - семинаре "Современные проблемы аэрогидромеханики", Сочи, "Буревестник", 2003; на конференции "Ломоносовские чтения", МГУ, 2004; на XII школе - семинаре "Современные проблемы аэрогидромеханики", Сочи, "Буревестник", 2004; на X Всероссийском совещании по проблемам построения сеток для решения задач математической физики, посвященной памяти академика А.Ф. Сидорова, 2004; на конференции "Ломоносовские чтения", МГУ, 2005; на XIII школе - семинаре "Современные проблемы аэрогидромеханики", Сочи, "Буревестник", 2005.

Публикации по теме диссертации

Результаты диссертация изложены в 15 печатных работах автора, приведенных в списке литературы. Из них журнальных статей — 2, статей в сборниках или трудах конференций — 9, препринтов — 4.

Объем и структура работы.

Диссертация состоит из введения, четырех глав и заключения. Объем диссертации составляет 113 стр., включает 32 рисунка, 3 таблицы, список литературы из 90 наименований.

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

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

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

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

1. Список литературы

2. Азиз X., Сеттари Э. Математическое моделирование пластовых систем. — М.: Недра, 1982, 408с.

3. Бабейко А.Ю., Динариев О.Ю. Моделирование ретроградной конденсации при стационарной радиальной фильтрации. — МЖГ, 1994, N 6, с.92-97.

4. Басниев К. С., Бедриковецкий П.Г. Многофазное вытеснение смешивающихся жидкостей и газов из пористых сред. — Итоги науки и техники. Комплексные и специальные разделы механики. М.: ВИНИТИ, 1988, т.З, с. 81-162.

5. Баталии О.Ю., Брусиловский А.И., Захаров М.Ю. Фазовые равновесия в системах природных углеводородов. — М.: Недра, 1992.

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

7. Баренблатт Г.И. Подобие, автомодельность, промежуточная ассимптотика: Теория и приложения к геофизической гидродинамике.— Л.: Гидрометеоиздат, 1982, 254 с.

8. Бедриковецкий П.Г., Лурье М.Б. Устойчивость и допустимость разрывов в системах уравнений двухфазной фильтрации. — ПММ, 1983, т.17, вып.4, с.590-600.

9. Бедриковецкий П.Г., Каневская Р.Д., Лурье М.В. Автомодельные решения задач двухфазной фильтрации с учетом сжимаемости одной из фаз. — МЖГ, 1990, N 4, с.71-80.

10. Вабищевич П.Н., Павлов А.Н., Чурбанов А.Г. Методы расчета нестационарных несжимаемых течений в естественных переменных на неразнесенных сетках. — Мат.моделирование, 1996, т.8, N 7, с.81-108.

11. Волосевич П.П., Леванов Е.И. Автомодельные решения задач газовой динамики и теплопереноса. — М.: Из-во МФТИ, 1997г.

12. Восков Д.В., Ентов В.М. К задаче о вытеснении нефти смесями газов. — Изв. РАН. Сер. МЖГ, 2001, N 2, с. 112-121.

13. Вязников К.В., Тишкин В.Ф., Фаворский А.П. Построение монотонных разностных схем повышенного порядка аппроксимации для1. СПИСОК ЛИТЕРАТУРЫ 107систем уравнений гиперболического типа.— Математическое моделирование, 1989, N 1, с.95.

14. Годунов С.К., Рябенький B.C. Разностные схемы (введение в теорию). — М.: Наука, 1977.

15. Годунов С.К. Элементы механики сплошной среды. — М.: Наука, 1978.

16. Годунов С.К., Роменский Е.И. Элементы механики сплошных сред и законы сохранения. — Университетская серия 4, Научная книга, Новосибирск, 1998.

17. Городецкий Е.Е., Куликов В.Д.,Федюнина JI.B., Анисимов М.А. Изоморфное описание двухфазной области критических бинарных растворов. — ЖЭТФ, 1997, т. Ill, N 1, с. 1-7.

18. Зазовский А.Ф. Структура скачков в задачах вытеснения нефти химреагентами, влияющих на фазовое равновесие. — МЖГ, 1985, N 5, с. 116-126.

19. Зазовский А.Ф. Двухфазная трехкомпонентная фильтрация с переменным суммарным потоком. — Изв. АН СССР, МЖГ, 1985, N 3, с. 113-120.

20. Ентов В.М., Зазовский А.Ф. Гидродинамика процессов повышения нефтеотдачи. — М.: Недра, 1988, 231 с.

21. Ентов В.М. Математические модели повышения нефтеотдачи нагнетанием газов. — Наука и технология углеводородов, 2000, N 5, с. 166-172.

22. Закиров С.Н., Сомов Б.Е., Гордон В.Я. и др. — Многомерная многокомпонентная фильтрация : Справочное пособие. — М.: Недра, 1988, 335 стр.

23. Кабанов Н.И., Панфилов М.Б. Ретроградное расслоение газового кондесата на две жидкие углеводородные фазы. — Газовая промышленность, 1996, N 5, с.57-58.

24. Колдоба A.B., Повещенко Ю.А., Самарская Е.А., Тишкин В.Ф. Методы математического моделирования окружающей среды. — М.: "Наука", 2000, 254 с.

25. Колдоба Е.В. О термодинамике газоконденсатных смесей в пористой среде. — Препринт N 1 ИПМ РАН, 1999, 28с.1. СПИСОК ЛИТЕРАТУРЫ 108

26. Колдоба Е.В. Автомодельные решения двухкомпонентной фильтрации с фазовыми переходами. — Тез. докл. науч. конф. "Ломоносовские чтения", МГУ, 2005, с. 118.

27. Колдоба А.В., Колдоба Е.В. Численное моделирование двухкомпонентной фильтрации. — Препринт N84 ИПМ РАН, 1999, 22с.

28. Колдоба А.В., Колдоба Е.В. Разрывные решения многокомпонентной фильтрации.— Препринт N85 ИПМ РАН, 1999, 20с.

29. Колдоба А.В., Колдоба Е.В. Тестирование и численное моделирование многокомпонентной фильтрации с фазовыми переходами. — Препринт N39 ИПМ РАН, 2002, 22с.

30. Koldoba A.V., Koldoba E.V. The Discontinuous Solutions of the Transport Equations for Compositional Flow in Porous Media. — Transport in Porous Media, vol.52, no.2, August 2003, p. 267-277.

31. Колдоба А.В., Колдоба Е.В. Модельное уравнение состояния и потенциал Гиббса для численного расчета задач многокомпонентной фильтрации с фазовыми переходами. — Геохимия, N5, 2004, с. 573576.

32. Koldoba А. V., Koldoba E.V. The Free Energy Condition for Compositional Flow in Porous Media.— International conference New Trends in Continuum mechanics. Constatza, Ruminia, 2003, p. 116.

33. Koldoba A.V., Koldoba E.V. Procedure of "Entropy" Correction for the Transport Equitions for Compositional Flow in Porous Media.— International conference New Trends in Continuum mechanics. Constatza, Ruminia, 2003, p. 117.1. СПИСОК ЛИТЕРАТУРЫ 109

34. Колдоба A.B., Колдоба Е.В. Уравнение баланса свободной энергиидля многокомпонентной многофазной фильтрации и его применение. — Тез. докл. науч. конф. "Ломоносовские чтения", МГУ, 2003, с. 78.

35. Колдоба A.B., Колдоба Е.В. Распространение слабых разрывов концентрации при фильтрации многокомпонентной смеси. — Тез. докл. XI школы семинара " Современные проблемы аэрогидродинамики", Сочи, "Буревестник", МГУ, 2003, с.93.

36. Колдоба A.B., Колдоба Е.В. О термодинамическом согдасовании моделей многокомпонентной фильтрации с фазовыми переходами. — Тез. докл. XII школы семинара "Современные проблемы аэрогидродинамики", Сочи, "Буревестник", МГУ, 2004, с.97.

37. Колдоба A.B., Колдоба Е.В. Автомодельные решения типа бегущей волны для двухкомпонентной фильтрации с фазовыми переходами. — Тез. докл. XII школы семинара "Современные проблемы аэрогидродинамики", Сочи, "Буревестник", МГУ, 2005, с.54.

38. Колдоба Е.В. Автомодельные решения двухкомпонентной фильтрации с фазовыми переходами. — Тез. докл. науч. конф. "Ломоносовские чтения", МГУ, 2005, с. 118.

39. Колесниченко A.B., Максимов В.М. Методы неравновесной термодинамики для моделирования многофазного многокомпонентного континуума. — Сб.науч.Трудов, МИНГ, 1986, Вып. 200, с. 191.

40. Коновалов А.Н. Задачи фильтрации многофазной несжимаемой жидкости. — Новосибирск: Наука, Сиб.отд-е, 1988, 166с.

41. Курант Р. Уравнения с частными производными. — М.: Мир, 1964, 830с.

42. Куропатенко В.Ф. Обмен импульсом и энергией в неравновесных многокомпонентных средах//Прик. механика и техническая физика. 2005. Т. 46, 1. С. 7-15.

43. Куропатенко В.Ф. Уравнения состояния в математических моделях механики и физики //Экстремальные состояния вещества: Сб.науч.тр. АН СССР. ИВТАН, 1991. С. 3-38.

44. Куликовский А.Г., Погорелое Н.Г., Семенов А.Ю. Математические вопросы численного решения гиперболических систем уравнений. — М.: Физматлит, 2001, с.608.1. СПИСОК ЛИТЕРА ТУРЫ 110

45. Куликовский А.Г., Свешникова Е.И. Нелинейные волны в упругихсредах. — Московский лицей, 1998, с.412.

46. Ландау Л.Д., Лифшиц Е.М. Статическая физика, том 5. — М: Наука, 1995, с.584.

47. Ландау Л.Д., Лифшиц Е.М. Гидродинамика, том 6. — М: Наука, 1986, с.736.

48. Лурье М.В., Максимов В.М., Филинов М.В. Исследование различных случаев взаимного вытеснения несмешивающихся жидкостей в пористой среде. — Инж.-физ. жур., 1981, Т.41, N 4, с.656-662.

49. Максимов М.М., Рыбицкая Л.П. Математическое моделирование процессов разработки нефтяных месторождений. — М.: Недра, 1976, 264с.

50. Николаевский В.Н., Бондарев Э.А., Миркин М.И., Степанова Г.С., Терзи В.П. Движение углеводородных смесей в пористой среде. — М.: Недра, 1968.

51. Николаевский В.Н. Механика пористых и трещиноватых сред. — М.: Недра, 1984, с. 242.

52. Панфилов М.Б., Кайгородова М.В., Бороздлк О.И. Прогнозирование разработки газоконденсатных месторождений. — Науч. техн. обзор ВНИИЭГазпром, М., 1988, вып.7, РиЭГГМ, 36с.

53. Панфилов М.В.,Панфилова И.В. Осредненные модели фильтрационных процессов с неоднородной внутренней структурой. — М.: Наука, 1996, 383с.

54. Панфилов М.В. Фильтрация жидкостей с сильно различающимися подвижностями при наличии фазовых переходов: явления пограничного слоя, пространственные фазовые структуры, неустойчивость течений. — Изв. РАН, МЖГ, 2002, N 2, с.124-135.

55. Пригожин И., Дефэй Р. Химическая термодинамика. — Новосибирск: Наука, 1966, 509 с.(перевод с анг.: I.Prigogine and R.Defay Chemical thermodynamics. Longmans Green and Co, London-New York-Toronto, 1954)

56. Pud P., Праусниц Дж., Шервуд Т. Свойства газов и жидкостей. — JL: Химия, 1982.1. СПИСОК ЛИТЕРАТУРЫ 111

57. Розенберг М.Д., Кундин С.А. Многофазная многокомпонентнаяфильтрация при добыче нефти и газа. — М.: Недра, 1978.

58. Рождественский Б.Л., Яненко Н.Н. Системы квазилинейных уравнений и их приложение к газовой динамике. — М.: Наука, 1978.

59. Самарский А.А. Введение в численные методы. — М.: Наука, 1982, 272с.

60. Самарский А.А., Вабищевич П.Н. Разностные схемы для уравнений переноса. I — Дифференциальные уравнения, 1998, Т.34, с. 1675-1685.

61. Селлков В.И., Кадет В.В. Перколяционные модели процессов переноса в микронеоднородных средах.— М.: Недра, 1995.

62. Сивухин Д.П. Общий курс физики, т.2, Термодинамика и молекулярная физика. — М.: Наука, 1979.

63. Уэйлес С. Фазовые равновесия в химической технологии: в 2-х частях. — М.: "Мир", 1989.

64. Филинов М.В., Марон В.И., Рохлин И.М. Об учете сжимаемости при фильтрации многофазной жидкости. — Тр. Моск. ин-та нефте-хим. и газ. пром-ти, 1969, Вып. 79, с. 37-45.

65. Хейфец Л.И., Наймарк А.В. Многофазные процессы в пористых средах. — М.: Химия, 1982.

66. Христианович С.А. Неустановившееся течение жидкости и газа в пористой среде при резких изменениях давления во времени при больших градиентах пористости. ФТПРПИ N 1, с.3-18, Новосибирск, Наука, 1985, с. 46.

67. Voronov V.P., Belyakov M.Yu., Gorodetskii E.E., Kulikov V.D., Muratov A.R., Nagaev V.B. Phase behavior of methane-pentan mixture in bulk and porous media. — Transport in Porous Media, vol.52, no.2, August 2003, p.123-140.

68. Barenblatt, G.I., Entov, V.M., and Ryzhik, V.M. Theory of Fluid Flows Through Natural Rocks. — Dordrecht, Kluwer Academic Publishers, 1990.

69. Bedrikovetsky, P. ^.Mathematical Theory of Oil & Gas Recovery (with applications to ex-USSR oil & gas condensate fields). — Kluwer Academic Publishers, London-Boston-Dordrecht, 1993.1. СПИСОК ЛИТЕРАТУРЫ 112

70. Buckley S.E. and Leverett M.S. Mechanizm of Fluid Displacement insands — Trans. AIME, 1942, v.146, p.107-116 .

71. Courant R., Friedrichs K.O. Supersonic Flow and Dhock Waves. — Applied Mathematical Sciences 21, Springer, New York. Рус. пер.: Курант P., Фридрихе К. Сверхзвуковое течение и ударные волны. — Иностр. лит., Москва, 1950.]

72. David G.Schaefer and Michael Shearer The classification of 2x2 system of non-strictly hyperbolic conservation laws, with application to oil recovery. — Communications on Pure and Applied Mathimatics, V. XL, 141-178 (1987).

73. Einfeld B. On Godunov-type methods for gas dynamics. — SI AM J. Numer. Anal., 25, 294, 1988.

74. Entov V.M. Nonlinear waves in physicochemical hydrodynamics of enhanced oil recovery. Multicomponent flows. — Porous media: Physics, Models, Simulation: Proc Inter. Conf, Moscow, Russia, 1997.

75. Ewing R.E. Simulation of multi-phase flows in porous media. — Transport in Porous Media, N 6, p.49-63, 1991.

76. Harten A., Hyman J.M., Lax P.D. On finite-difference approximations and entropy conditions for shocks. — Comm. Pure Appl. Math., No. 3, 297-322, 1976.

77. Johns R.T., Dindoruk В., Orr F.M. Analytical theory of combined condensing/ vaporising gas drives. — Soc. Petr. Engrs Adv. Technol. Series, 1993, v. 1, N 2, p. 7-16.

78. Hirsch C. Numerical Computation of internal and external flows. — Computational methods for iinviscous and viscous flow, 1990, v.2, 6911. P

79. Lake L. W. Enhanced oil recovery. — Englewood Cliffs: Prentice Hall, 1989, 550p.

80. Lax P.D. Hyperbolic Conservation Laws II,— Comm. Pure Appl.Math.,1957, 10, p.537-566.

81. Lax P.D. Hyperbolic Systems of Conservation Laws and the Mathematical Theory of Shock Waves. — no 11; SIAM: Philadelphia, 1973.1. СПИСОК ЛИТЕРАТУРЫ (113

82. Leung S.S., Griffiths R.B. Thermodinamic properties near the liquicvapour critical line in mixtures of Hez and He4. — Phys.Rev.A, 1973, v.8(5), p.2673-2686.

83. Muskat M., Meres M.W. Flow of heterogeneous fluid though porous media — Physics, 1936, v.7, N 9, p. 346-363.

84. Muskat M. Physical principles of oil production— N.Y.:McGraw-Hill, Inc. 1949, 922 (Рус.перев.: Маскет M. Физические основы технологии добычи нефти. М.: Наука, 1969, 545с.)

85. Mitlin V.S. Two-phase multicomponent filtration: instabilities, autowaves and retrograde phenomena. — J. Fluid Mech., 1990, v. 220, p. 369-395.

86. Orr F.M., Johns R.T., Dindoruk B. Development of misciblity in four-component CO2 floods. — Soc. Petr. Res. Eng., 1993, v. 8, N 2, p. 135142.

87. Orr L. Theory of Gas Injection Processes. — Stanford University, California, 1999.

88. Schottky W., Ulich H., Wagner C. Thermodynamik. — Berlin, 1929.