Численное моделирование процессов массопереноса на
призабойном участке нефтяной скважины с учетом
неравномерности притока нефтеводяной смеси с несущей
нефтяной фазой по толщине слоисто-неоднородного пласта
ВВЕДЕНИЕ 3
1 Математическая модель двухфазного потока на забое скважины с учетом
фильтрационного притока из пласта 5
1.1 Уравнение неразрывности однофазного и двухфазного потоков 5
1.2 Уравнения нестационарного нефтеводяного потока в скважине 8
2 Численное моделирование нестационарных процессов 12
2.1 Разностная схема расчета нестационарного потока в скважине 12
2.2 Численный алгоритм решения задачи методом простой итерации 13
2.3 Численный алгоритм решения задачи методом Ньютона 14
2.4 Получение критерия Куранта для нелинейной задачи переноса 18
3 Описание программы 21
4 Результаты расчетов 26
4.1 Сравнение итерационных методов решения задачи 26
4.1.1 Симметричный пласт 26
4.1.2 Несимметричный пласт сложной структуры 30
4.1.3 Симметричный пласт с изолированным участком 34
4.2 Исследование сходимости разностной схемы для метода простой итерации и
устойчивости по критерию Куранта 38
4.3 Влияние структуры потока на его параметры 42
ЗАКЛЮЧЕНИЕ 56
СПИСОК ИСПОЛЬЗОВАННЫХ ИСТОЧНИКОВ 57
ПРИЛОЖЕНИЕ. Листинг программы 58
В процессе освоения месторождений, при пуске и остановке для проведения геофизических, гидродинамических исследований или ремонта скважина работает в нестационарном режиме [1]. Наиболее распространены пласты слоистого строения, состоящие из неоднородных абсолютной проницаемости и по толщине пропластков с различными фильтрационно-емкостными характеристиками. Процесс освоения скважины осложняется также неоднородностью смеси, в состав которой входят нефтяная и водяная фазы с существенно-различными физико-химическими свойствами, которые поступают из пласта на перфорированном участке обсадной колонны. Поэтому актуальными задачами анализа процессов работы скважины являются: создание обобщенной модели
нестационарной работы скважины в слоисто-неоднородном пласте, определение степени влияния данного столба на нестационарную работу скважины, исследование внутрискважинных перетоков жидкости между пропластками с различными фильтрационными свойствами.
На рисунке 1 в качестве примера схематично показан фрагмент вертикального разреза симметричного горизонтального пласта, который имеет сложную слоистонеоднородную трещиновато-пористую структуру.
Пласт состоит из N слоев, имеющих различную толщину Н = yt- , абсолютную
проницаемость Kt > 0 и динамическую пористость тг, I = 1,2,...,N . Параметры трещин
и пористых блоков обозначены на рисунке индексами «f» и «Ь». Общая толщина
коллектора H = H + Н2 +... + Нд,. Координатная ось Or направлена вдоль его кровли, а
ось Oz - вдоль оси добывающей скважины, боковая поверхность забоя которой
расположена на левой границе пласта при r = r0. Начало координат (z = 0, r = 0)
находится на уровне подошвы пласта. Поверхности z = 0 и z = Н его кровли и подошвы
непроницаемы. Границы уг (/ = 1,2,...,JV-1) пропластков проницаемы, так что они
гидродинамически связаны между собой. Скважина может вскрывать как все слои, так только некоторые из них. Вскрытые перфорацией "проницаемые" (yL ,уя) и невскрытые
"непроницаемые" (Г2, Гл ) участки левой и правой боковых границ пласта при r = r0 и
r = Rr показаны на рисунке штриховыми и сплошными линиями соответственно. Стрелки,
расположенные в окрестности левой границы yt, иллюстрируют направление потока
водонефтяной смеси, поступающего из пласта в обсадную колонну добывающей скважины.
Целью работы является:
- создание обобщенной модели нестационарных процессов на призабойном участке скважины с учетом притока нефтеводяной смеси с несущей нефтяной фазой из слоистонеоднородного пласта;
- разработка разностной схемы и алгоритма численной реализации модели;
- разработка программы на основе численной и алгоритмической моделей массопереноса на призабойном участке нефтяной скважины с учетом неравномерности притока нефтеводяной смеси по толщине слоисто-неоднородного пласта;
- численное исследование на основе многовариантных вычислительных экспериментов эффекта внутрискважинных перетоков жидкости между пропластками с различными фильтрационными свойствами и оценки продолжительности переходных процессов при освоении добывающей скважины, вскрывающей слоисто-неоднородноый пласт.
В данной работе удалось создать математическую модель неравномерного притока в
скважину нефтеводяной смеси с несущей нефтяной фазой из нефтяного пласта сложной
структуры, оценить неоднородности потока, исследовать выход нестационарного
процесса на квазистационарный режим.
Также была проведена разработка разностной схемы и алгоритма численной
реализации модели с использованием двух методов – метода простой итерации и
итерационного метода Ньютона. Разработана программа на основе численной и
алгоритмической моделей массопереноса на призабойном участке нефтяной скважины с
учетом неравномерности притока нефтеводяной смеси с несущей нефтяной фазой по
толщине слоисто-неоднородного пласта.
На основе проведенных экспериментов можем сделать такие выводы, как:
разностная схема является устойчивой и сходится «в себе»;
применение критерия Куранта позволяет на каждом временном слое находить
значения временного шага, что обеспечивает условие устойчивости решения задачи
переноса. При этом время в десятки раз сокращается, нежели при расчетах с
фиксированными временными шагами, где соблюдается устойчивость схемы во всем
временном диапазоне решения задачи;
многовариантные расчеты можно проводить при достаточно больших
пространственных шагах сетки с использованием критерия Куранта для определения
временного шага разностной схемы.
На основе многовариантных вычислительных экспериментов был выявлен эффект
внутрискважинных перетоков жидкости между пропластками с различными
фильтрационными свойствами и оценка продолжительности переходных процессов при
освоении добывающей скважины, вскрывающей слоисто-неоднородный пласт.
Все поставленные цели были достигнуты