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

• Автор: Колдоба, Елена Валентиновна
• Шифр специальности: 05.13.18
• Научная степень: Кандидатская
• Год защиты: 2005
• Место защиты: Москва
• Количество страниц: 113 с.

Оглавление

Введение

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 Разрывные решения уравнений двухкомпонентной фильтрации, графический метод анализа
термодинамического
условия
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]. В рамках модели Баклея-Леверетта, система уравнений двухфазной фильтрации расщепляется на гиперболи-

График невыпуклой функции W(t) f dx К I
r dx
/rw=w-P2- <24>
Так как p > P2, то полная скорость фильтрации W > 0.
В третьем уравнении (12) перейдем к новой независимой переменной z, которую введем из соотношения dz — W(t)dt/m. Тогда это уравнение примет вид
1+ = ° <“> с начальными условиями 5(0,х) = 0. На левой границе: 5(2,0) = 1.
В характеристической форме уравнение (25) имеет вид

лS“°-
Значит на каждой характеристике s = const, и тогда из первого соотношения вытекает, что все характеристики уравнения (25) суть прямые линии в (ж, г)-плоскости.