автореферат диссертации по , 05.00.00, диссертация на тему:Моделирование методом линий тока с капиллярным эффектом применительно к проблемам нефтедобычи

Эта работа была проделана в Центре Исследования Фазовых Превращений и Процессов Разделения (IVC-SEP), Факультет Химической Технологии, Датский Технический Университет с Августа 2001-го года по Август 2004-го года под руководством Профессора Эрлинга Стенбю (Erling Н. Stenby), доцента Микаэля Микельсена (Michael L. Michelsen) и доцента Александра Шапиро (Alexander A. Shapiro). Проект спонсировался Датским Энергетическим Управлением (Danish Energy Authority), а так же компаниями, Maersk Oil and Gas, ChevronTexaco и Eni-Agip. Я хотел бы выразить свою признательность Профессору Стинбю за три года в великолепной атмосфере IVC-SEP. Я очень благодарен своим научным руководителям за те знания, которыми они поделились со мной за этот период.

Я очень признателен Профессору Франклину Oppy (Franklin М. Orr, Jr.) за те чудесные 6 месяцев, которые я провел в группе SUPRI-C, Нефтепромысловый Факультет, Университет Стэнфорда. Время, проведенное в Стэнфорде, было очень важным для многих аспектов моей жизни. Для меня было большим удовольствием рабоать с Кристианом Ессеном (Kristian Jessen), Йилдираем Чинаром (Yildiray Cinar), Бредли Мэлиссоном (Bradly Malisson) и Марго Геринтсен (Margo Gerintsen). Для меня является счастьем и честью, что наша дружба и сотрудничество продолжаются до сих пор.

Я хотел бы поблагодарить моих друзей в Дании и США за дружбу и поддержку: Олега Медведева, Томаса Линдвига (Thomas Lindvig), Петра Железного, Мохаммада Карими-Фарда (Mohammad Karimi-Fard) и многих других. Я очень признателен моим друзьям в Москве за то, что они всегда были рядом независимо от расстояния разделявшего нас.

И, наконец, я бесконечно благодарен моим родителям за то, что они родили и вырастили меня, любят и поддерживают во всех моих начинаниях.

Краткое содержание

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

Метод моделирования линий тока относительно недавно появился на рынке коммерческих симуляторов с потенциалом стать одной из основных технологий моделирования коллекторов. Первый симулятор на базе линий тока появился около 10 лет назад. Преимущества метода линий тока над конечно-разностными симуляторами являются высокая скорость расчетов и результаты, менее подверженные численной дисперсии [7, 17, 19, 67]. Трассировка линий тока с учетом скорости течения Дарси [31] позволяет учесть нелинейности решения, связанные с подвижностью фаз и капиллярным давлением. Линии тока позволяют разделить сложное решение 3-х мерного уравнения насыщенности на набор одномерных решений вдоль линий тока в координатах времени-полета (time-of-flight -TOF) [17, 32]. Однако на данный момент симуляторы линий тока представляют значительно меньше возможностей по моделированию различных физических аспектов процессов вытеснения по сравнению с конечно-разностными симуляторами. Одним из основных недостатков двухфазного симулятора линий тока основанного на модели нелетучей нефти является невозможность учета капиллярных сил. В гетерогенных коллекторах с разной смачиваемостью породы капиллярные эффекты могут стать очень значимыми, если не доминирующими. История разработки симулятора на основе метода линий тока расмотренна в первой главе.

В этой работе представлен метод введения капиллярных эффектов в симулятор на базе линий тока. Первая

глава представляет собой введение в течение жидкостей в поровой среде и в теорию капиллярных давлений. Вторая

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

Четвертая

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

Этот текст завершается выводами и предложениями дальнейших направлений исследований.

Оглавление

Оглавление

Предисловие i

Краткое содержание н

Оглавление iv

1. Введение 1

1.1. Течение жидкости в пористой среде.1

1.2. Физика капиллярных эффектов.1

1.3. Введение в симуляторы на базе линий тока.1

1.3.1. Обзор симулятора 3DSL 0.25.1

1.3.2. Введение методологии линий тока и ранний этап развития.1

1.3.3. Развитие метода линий тока.1

1.3.4. Последние разработки в моделировании линий тока.1

2. Математическая модель 2

2.1. Модификация основной системы уравнений.2

2.1.1. Уравнение давления с учетом капиллярных сил.2

2.1.2. Капиллярно-вязкостной потенциал.2

2.1.3. Уравнение насыщенности с капиллярными и гравитационными эффектами.2

2.1.4. Конечный вид систем уравнений.2

2.2. Решение основных систем уравнений с учетом капиллярных и гравитационных сил.2

2.2.1. Схема модифицированного симулятора линий тока.2

2.2.2. Численное представление уравнения давления.2

2.2.3. Уравнение для скважин, учитывающее капиллярные силы.2

2.2.4. Расчет скорости потока.2

2.2.5. Решение уравнения насыщенности .2

2.2.6. Возможности для возникновения ошибки материального баланса.2

2.2.7. Выбор шага по времени.2

2.3. Иллюстрация капиллярных эффектов.2

3. Сравнение различных методов учета капиллярных сил 3-1 3.1. Данные для тестовых примеров.3

3.1.1. Первый тестовый пример.3

3.1.2. Второй тестовый пример.3

3.1.3. Третий тестовый пример.3

3.2. Определение режима вытеснения.3

3.3. Тестирование различных вариантов модификации уравнения давления .3

3.4. Сравнение различных вариантов модификации уравнения насыщенности.3

3.5. Тестирование методов автоматического выбора шага.3

3.6. Моделирование двухмерного гетерогенного примера.3

3.7. Расчет трехмерного гетерогенного примера.3

3.8. Резюме.3

4. Примеры расчетов 4

4.1. Капиллярные эффекты в поровом пространстве.4

4.1.1. Капиллярные эффекты в гидрофильной породе.4

4.1.2. Капиллярные эффекты в гидрофобной и смешанной средах.4

4.2. Определение области применения симулятора на базе линий тока с учетом капиллярных эффектов.4

4.3. Моделирование лабораторных экспериментов.4

Выводы.5

Список обозначений.6

Главной целью любого проекта по гидродинамическому моделированию пласта является прогнозирование характеристик процесса вытеснения.Симуляторы на базе линий тока позволяют бьютрее получать более точные решения задач вытеснения по сравнению с методами конечных разностей. Однако доступные в настоящее время коммерческие симуляторы на базе линий тока недостаточно хорошо описывают капиллярные эффекты.В данной главе представлен обзор базовых концепций течения жидкости в пористой среде. Рассматриваются скорость Дарси и закон сохранения вещества. Представлен метод IMPES (неявного решения уравнения давления; явного решения уравнения насыщенности), используемый в симуляторах на базе линий тока.Представлены физические основы капиллярных эффектов.Рассматривается значимость капиллярных сил для гидродинамического моделирования пласта. Приводятся возможные затруднения при вводе капиллярных эффектов в симулятор на базе линий тока.Данная глава включает в себя обзор методологии линий тока. История развития методов линий тока (трубок тока) прослеживается с ввода понятия функции линий тока [54] до развития симулятора на базе линий тока 3DSL 0.25 [8]. Симулятор на базе линий тока 3DSL 0.25 предоставлен группой SUPRI-C, нефтепромысловым факультетом Стэнфордского университета в качестве базового кода для ввода капиллярных эффектов.Представлены некоторые последние достижения в области развития линий тока, включая трёхфазный композиционный симулятор, симулятор на базе линий тока с двойной проницаемостью, использование симулятора на базе линий тока для настройки модели в соответствии с историей

5. Выводы.

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

• Двухфазный симулятор несмешивающихся, несжимаемых флюидов на базе метода линий тока (3DSL0.25 [7]) был модифицирован для учета капиллярных эффектов, как в уравнении давления, так и в уравнении насыщенности. Представлены и обсуждены несколько методов учета капиллярных эффектов.

• Созданный в этой работе симулятор CapSL был протестирован на широком спектре разнообразных расчетных задач. Представлены результаты расчетов нескольких 2-х и 3-х мерных искусственных примеров. Производительность и результаты расчетов симуляторов Eclipse и CapSL так же сравнивались на базе лабораторных экспериментов.

• Как метод прямой модификации уравнения давления (ПМУ), так и метод капиллярно-вязкостного потенциала (КВП) показывают результаты расчета сходные с симулятором Eclipse. Однако результаты, полученные по методу КВП, имеют меньшую ошибку материального баланса по сравнению с методом ПМУ. Метод КВП также быстрее. Автоматическая процедура выбора шага обычно предлагает больший размер шага при использовании метода КВП, что говорит о большей устойчивости этого метода.

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

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

• При расчетах 3-х мерной полноразмерной сетки результаты полученные на CapSL, отвечали оценкам режимов вытеснения, полученным по безразмерным параметрам Бедриковетского и Жу и др. Накопленная нефтедобыча практически не изменилась с введением капиллярных сил. Однако локально капиллярные силы оказывали значительное воздействие на профили вытеснения. В ситуации с Eclipse накопленная нефтедобыча при расчетах с и без капиллярных сил была разная, что не соответствует режиму вытеснения, оцененному по безразмерным параметрам. Такое несоответствие является результатом высокой численной дисперсии решения.

• CapSL менее подвержен влиянию численной дисперсии как на полно-размерных сетках, так и при расчетах лабораторных экспериментов. Как следствие этого, тонкие детали фронта вытеснения более точно описываются CapSL, а не Eclipse.

• CapSL полностью учитывает капиллярные эффекты, возникающие в коллекторах, содержащих породы с разной смачиваемостью.

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

0.25) это некоммерческое программное обеспечение, созданное в нескольких университетах. Код этого продукта не подвергался серьезной оптимизации с целью ускорения расчетов. Eclipse 100, с другой стороны - это высоко-оптимизированный коммерческий продукт. Однако расчеты проводились в режиме IMPES - наиболее медленном и наименее подверженном численной дисперсии алгоритме решения уравнений. Многими авторами было показано, что метод AIM незначительно увеличивает численную дисперсию решения и значительно быстрее метода IMPES.

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

• По результатам настоящей работы можно порекомендовать следующие направления развития симуляторов на базе линий:

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

2) Применение усовершенствованной процедуры трассировки линий тока и переноса насыщенностей между конечно-разностной сеткой и линиями тока [48] должно привести к ускорению расчета и уменьшению численной погрешности расчета;

3) Работа с более сложными сетками (по узловым точкам, нерегулярные сетки);

4) Введение более точных уравнений скважин;

5) Введение капиллярных эффектов в двух- и трехфазные капиллярные симуляторы.

6-1

Условные обозначения

Символы с Коэффициент С сек м2кг"1 Коэффициент капиллярной передачи для метода КВП

Сп Число Куранта

D м Глубина

Е Ошибка f Функция течения фазы в многофазном потоке

F Функция g м/с2 Ускорение свободного падения

Н, h м Высота

J Функция Леверетта к м2 Проницаемость

L, I м Длина

М Соотношение подвижностей фаз m, п Коэффициенты модели Ван-Генухтена п м"2 Количество капилляров на единицу поверхности (глава 1)

N сек м2кг"1 Коэффициент капиллярной передачи для метода ПМУ Безразмерный параметр (глава 3)

Р Па Давление

Q, q м3/с Дебит г м Радиус

R Форм-фактор t с Время

Т сек Передаточная функция (глава 1) м2кг1 Коэффициент передачи (глава 2)

U м/с Скорость

V м3 Объем х, у, z Оси

Неизвестные (секция1.3.2)

Символы греческого алфавита

6 м Расстояние от центра ячейки до поверхности контакта

Безразмерный параметр ф Пористость

Па Член в уравнении капиллярных сил, зависящий от

Ф свойств породы

Соотношение гравитационной и суммарной у подвижностей фаз

Т Па Капиллярно-вязкостной потенциал

X Па с Подвижность фазы

Л Среднее ц Вязкость

0 Угол смачивемости р кг/м плотность а Н/м Межфазное натяжение т с Время полета

Q Множитель капиллярно-вязкостного потенциала

Индексы с 9 i.j. к о ог г rowi rwor s t V w wi капиллярные, критические (в формуле1.22) гравитационные

Одна из фаз или направлений течения нефть

Остаточная нефть Относительная

Относительная по нефти при водонасыщенности

Относительная по воде при нефтенасыщенности

Поверхность

Накопленная

Вязкостной гидрофильная

Начальная вода начальной остаточной

1. Al-Huthali, A.H., and Datta-Gupta, A. Streamline Simulation of Water Injection in Naturally Fractured Reservoirs. SPE 89443 presented at the 2004 SPE/DOE Symposium on Improved Oil Recovery, Tulsa, USA. 2004.

2. Amyx, J.W., Bass, D.M., and Whiting, R.L. Petroleum Reservoir Engineering. McGraw-Hill, New York, USA. 1960.

3. Aziz, K., and Settari, A. Petroleum reservoir simulation. Elsevier applied science publishers, New York, USA. 1979.

4. Басниев, КС., Кочина И.Н., и Максимов, B.M,. Подземная гидромеханика. Москва, Россия. 1993.

5. Barenblatt, G.I., Entov, V.M., Ryzhik, V.M. Theory of fluid flows through natural rocks. Kluwer Academic publisher, Netherlands. 1990.

6. Batycky, R.P. A Three-Dimensional Two-Phase Field Scale Streamline Simulator. Ph.D. Thesis, Department of Petroleum Engineering, School of Earth Science, Stanford University, Stanford, California, USA. 1997.

7. Bedrikovetsky, P.G. Mathematical Theory of Oil and Gas Recovery With Application to ex-USSR Oil and Gas Fields. Kluwer Academic Publishers, Netherlands. 1993.

8. Bedrikovetsky, P.G., Magarshak, Т.О., and Shapiro A.A. 3D Analytical Model for Displacement of Oil Using Horizontal Wells. SPE 26996 presented at III Latin American / Caribbean Confernce held in Buenos Areas, Argentina. 1994.

9. Bedrikovetsky, P.G. New Mathematical Model for EOR Displacements Honoring Oil Ganglia. SPE 38892 presented at the SPE Annual Technical Conference and Exhibition, San Antonio, Texas, USA. 1997.

10. Bedrikovetsky, P.G. WAG Displacemetns of Oil-Condensates Accounting for Hydrocarbon Ganglia. Transport in Porous Media, Vol.52, pp. 229266. 2003.

11. Bedrikovetsky, P.G. Analytical Model for the Waterflood Honoring Capillary Pressure. SPE 36130, presented at the Fourth Latin and Caribbean Petroleum Engineering Conference held in Port-of-Spain, Trinidad and Tobago. 1996.

12. Berenblyum, R.A., Shapiro, A.A., Jessen, K., Stenby, E.H., and Orr, F.M. Jr. Black Oil Streamline Simulator With Capillary Effects. SPE 84037, presented at the SPE Annual Technical Conference and Exhibition held in Denver, Colorado, USA. 2003.

13. Berenblyum, R.A., Shapiro, A.A., Stenby E.H. A Method For Simulating A Multiphase Flow. Patent Application PA 2003 01844, Denmark. 2003.

14. Berre, I., Dahle, H.K., Karlsen, K.H., and Nordhaug, H.F. A Streamline Front Tracking Method for Two- and Three-Phase Flow Including Capillary Forces. Mathematics Subject Classification, pp 1-13. 2000.

15. Blair, P.M., and Weinlaug, C.F. Solution of Two-Phase Flow Problems Using Implicit Difference Equations. SPE 2185, Presented at 43rd Annual Fall Meeting, Houston, USA. 1969.

16. Blunt, M.J., Lui, K., Thiele, M.R. A Generalized Streamline Method to Predict Reservoir Flow. Petroleum Geoscience, vol. 2, pp. 259-269. 1996.

17. Bratvedt, F., Gimse, Т., and Tegnander, C. Streamline Computations for Porous Media Flow Including Gravity. Transport in Porous Media, vol. 25, pp. 63-78. 1996.

18. Buckley, S.E., and Leverett, M.C. Mechanisms of Fluid Displacement in Sands. Petroleum Technology, May, pp107-116. 1941.

19. Christie, M.A., and Blunt, M.J. Tenth SPE Solution Project: A comparison of Upscaling Techniques. SPE72469 in SPE Reservoir Evaluation and Engineering, August, pp. 308-317. 2001.

20. Coats, K.H. Implicit Compositional Simulation of Single-Porosity and Dual-Porosity Reservoirs. SPE 18427, presented at SPE Symposium On Reservoir Simulation, Huston, USA. 1989.

21. Coats, K.H. A Note on IMPES and Some IMPES-Based Simulation Models. SPE 65092 in SPE Journal, vol. 5, pp. 245-251. 2000.

22. Coats, K.H. IMPES Stability: The CFL Limit. SPE 66345, presented at SPE Reservoir Simulation Symposium, Houston, USA. 2001

23. Coats, K.H. IMPES Stability: The Stable Step. SPE 69225, presented at SPE Reservoir Simulation Symposium, Huston, USA. 2001.

24. Coats, K.H. IMPES Stability: Selection of Stable Timesteps. SPE 84924 in SPEJ, June, pp. 187-187. 2003.

25. Crane, M., Bratvedt, F., Bratvedt, K., Childs, P., and Olufsen, R. A Fully Compositional Streamline Simulator. SPE 63156 presented at SPE Annual Technical Conference and Exhibition, Dallas, Texas, USA. 2000.

26. Datta-Gupta, A., and King, M.J. A Semianalytic Approach to Tracer Flow Modeling in Heterogeneous Permeable Media. Advances in Water Resources, vol. 18, pp. 9-24. 1995.

27. Di Donato, G., Huang, W., and Blunt, M. Streamline-Based Dual Porosity Simulation of Fractured Reservoirs. SPE 84036 presented at SPE Annual Technical Conference and Exhibition, Denver, Colorado, USA. 2003.

28. Dullien, F.A.L. Porous media. Fluid Transport and Pore Structure. Academic Press, New York, USA. 1979.

29. Emanuel, S., and Milliken, W.J. History Matching Finite Difference Models With 3D Streamlines. SPE 49000 presented at SPE Annual Technical Conference and Exhibition, New Orleans, Louisiana, USA. 1998.

30. Farkas, E. Adaptice Implicit Volume Balance Techniques. SPE 25273 presented at 12th SPE Symposium on Reservoir Evaluation, New Orleans, Louisiana, USA. 1993.

31. Higgins, R.V., and Leighton, A.J. Performance of Five-spot Water Floods in Stratified Reservoirs Using Streamlines, SPE 57 in Trans. AIME, pp.112. 1961.

32. Higgins, R.V., and Leighton, A.J. A Computer Method to Calculate Two-Phase Flow in Amy Irregular Bounded Porous Medium. SPE 243, presented at SPE Annual California Regional Meeting, USA. 1962.

33. Higgins, R.V., and Leighton, A.J. Computer Prediction of Water Drive of Oil and Gas Mixture Through Irregularly Bounded Porous Media Three-Phase Flow. SPE 283, presented at SPE production reservoir Symposium, Tulsa, Oklahoma, USA1962.

34. Ichiro, O., Datta-Gupta, A., and King, M.J. Time step selection during streamline simulation via transverse flux correction. SPE 79688, presented at 2003 SPE reservoir simulation symposium, Huston, USA. 2003.

35. Jessen, К., Wang, Y., Ermakov, P., Zhu, J., and Orr, F.M. Fast, Approxmate Soilutions for 1D Mulitcomponent Gas-Injection Problems. SPE 74700 in SPEJ, December, pp. 442-451. 2001.

36. Jessen, K., and Orr, F.M. Compositional Streamline Simulation. SPE 77379 presented at SPE Annual Technical Conference and Exhibition, San Antonio, Texas, USA. 2002.

37. Jessen, K., and Orr, F.M. Gas Cycling and the Development of Miscibilty in Condensate Reservoirs. SPE 84070 presented at SPE Annual Technical Conference and Exhibition, Denver, Colorado, USA. 2003.

38. Jessen, K., and Orr, F.M. Gravity segregation and Compositional Streamline Simulation. SPE 89448, presented at SPE/DOE symposium on Improved Oil Recovery, Tulsa, Oklahoma, USA. 2004.

39. Kiilough, J.E. and Kossack, C.A. Fifth Comparative Solution project: Evaluation of Miscible Flood Simulators. SPE 16000 presented at Ninth SPE Symposium on Reservoir Simulation, San Antonio, Texas, USA. 1987.

40. King, M.J., and Datta-Gupta, A. Streamline Simulation: A current perspective. In Situ, vol. 22(1), pp. 91-140. 1998.

41. Leverett, M.C. Capillary behavior in Porous Solids. Trans. AIME. 1941.

42. Mallison, B.T., Gerritsen, M.G, and Matringe S.F. Improved Mappings for Streamline-Based Simulation. SPE 89352 presented at the 2004 SPE/DOE Symposium on Improved Oil Recovery, Tulsa, Oklahoma, USA. 2004.

43. Marcell-De Silva, J., and Dawe, R.A. Effects of Permeability and Wettability Heterogeneities on Flow in Porous Media. SPE 81164 presented at SPE Latin American and Caribbean Petroleum Engineering Conference, Port-of-Spain, Trinidad, West Indies. 2003.

44. Martin, J.C., and Wegner, R.E. Numerical Solution of Multiphase, Two-Dimensional Incompressible Flow Using Stream-Tube Relationships. SPE 7140 in Trans. AIME, October, pp. 313-323. 1979

45. Matanga, G.B. Stream Functions in Three-Dimensional Groundwater Flow. Water resources Research, vol. 29(9), pp 3125-3133. 1993.

46. Mattax, C.C., and Dalton, R.L. Reservoir Simulation. SPE Monogram Series. SPE, Richardson, Texas, USA. 1990.

47. Miiliken, W.J., Emanuel, A.S., and Chakravarty, A. Application of 3D Streamline Simulation To Assist History Matching. SPE 74712, presented at SPEREE December, pp. 502-508. 2001.

48. Muskat, M. Flow of Homogeneous Fluids. McGraw-Hill, Boston, USA. 1937.

49. Peaceman, D.W. Fundamentals of numerical reservoir simulation. Elsevier, New York, USA. 1977

50. Peaceman, D.W. Interpretation of Well-Block Properties in Numerical Reservoir Simulation with Nonsquare Grid Blocks and Anisotropic Permeability. SPE 10528 in SPEJ, June, pp. 531-569. 1983

51. Pollock, D.W. Semianalytical Computation of Path Lines for Finite-Difference Models. Ground Water, vol.26(6), pp. 743-750. 1988.

52. Rapoport, L.A., and Leas, W.J. Properties of Linear Waterfloods. Petroleum Transactions AIME, vol.198, pp. 139-148. 1953.

53. Renard, G. A 2D Reservoir Streamtube EOR Model With Periodical Automatic Regeneration of Streamlines. In Situ, vol. 14(2), pp. 175-200. 1990.

54. Rodriguez, P.G., Segura, M.K., and Moreno, F.J.M. Streamline methodology using an efficient operator splitting for accurate modeling of capillary and gravity effects. SPE 79693, presented at SPE Reservoir Simulation Symposium, Houston, USA. 2003.

55. Russell, T.F. Stability Analysis and Switching Criteria for Adaptive Implicit Methods Based on the CFL condition. SPE 18416 presented at SPE Symposium on Reservoir Simulation, Houston, Texas, USA. 1989.

56. Schlumberger GeoQuest. Eclipse 100 Technical Description. 2003.

57. Schlumberger GeoQuest. Eclipse 300 Reference Manual. 2003.

58. Schlumberger GeoQuest. FrontSim Technical Description. 2003.

59. Seto, C.J., Jessen, K., and Orr, F.M. Compositional Streamline Simulation of Field Scale Condensate Vaporization by Gas Injection. SPE 79690 presented at SPE Reservoir Simulation Symposium, Houston, Texas, USA. 2003.

60. Shapiro, A.A, and Stenby, E.H. Factorization of Transport Coefficients in Macroporous Media. Transport in Porous Media, vol. 41, pp. 305-323. 2000.

61. Spivak, A., and Coats, K.H. Numerical Simulation of Coning Using Implicit Production Terms. SPE 2595 in SPEJ, September, pp. 257-267. 1970.

62. Thiele, M.R. Modeling Multiphase Flow in Heterogeneous Media Using Streamtubes. PhD dissertation, Stanford University, Dept. of Petroleum Engineering, Stanford, California, USA. 1994.

63. Thiele, M.R;, Batycky, R.P., Blunt, M.J., and Orr, F.M. Jr. Simulating Flow in Heterogeneous Media Using Streamtubes and Streamlines. SPE 27834 in SPE Reservoir Engineering, vol. 10(1), pp. 5-12. 1996.

64. Thiele, M. R., Batycky, R.P., and Blunt, M.J. A Streamline-Based 3D Field-Scale Compositional Reservoir Simulator. SPE 38889 presented at SPE Annual Technical Conference and Exhibition, San Antonio, Texas, USA. 1997.

65. Thiele, M.R. Streamline Simulation. 6th international forum on Reservoir Simulation, Austria. 2001.

66. Thiele, M.R., and Batycky, R.P. Water Injection Optimization Using a Streamline-Based Workflow. SPE 84080 presented at SPE Annual Technical Conference and Exhibition, Denver, Colorado, USA. 2003.

67. Thomas, G.W., and Thumau, D.H. Reservoir Simulation Using an Adaptive Implicit Method. SPE 10120 presented at the 56th Annual Fall Technical Conference and Exhibition of the SPE of AIME, San Antonio, Texas, USA. 1981.

68. Thomas, G.W., and Thurnau, D.H. The Mathematical Basis of the Adaptive Implicit Method. SPE 10495 presented at the 6th SPE Symposium in Reservoir Simulation, New Orleans, Louisiana, USA. 1982.

69. Toth, J. Adsorption: Theory, Modeling, and Analysis. Marcel Dekker, New York, USA. 2002.

70. Yan, W., Michelsen, M.L., Stenby, E.H., Berenblyum, R.A., and Shapiro, A.A. Three-phase Compositional Streamline Simulation and Its Application to WAG. SPE 89440 presented at the 2004 SPE/DOE Symposium on Improved Oil Recovery, Tulsa, USA. 2004.

71. Young, L.C., and Russell, T.F. Implementation of an Adaptive Implicit Method. SPE 25245 presented at 12th SPE Symposium on Reservoir Simulation, New Orleans, Louisiana, USA. 1993.

72. Zhou, D., Fyaers, F.J., and Orr, F.M. Scaling multiphase flow in simple heterogeneous porous media. Multiphase Transport in Porous Media, vol. 265, pp. 25-35. 1993.