ISSN 0236-235X (P)
ISSN 2311-2735 (E)

Journal influence

Higher Attestation Commission (VAK) - К1 quartile
Russian Science Citation Index (RSCI)

Bookmark

Next issue

4
Publication date:
09 December 2024

Dual porosity model for fractured porous reservoirs development analysis based on the superelement concept

Date of submission article: 18.06.2019
UDC: 622.276.1/4
The article was published in issue no. № 3, 2019 [ pp. 478-485 ]
Abstract:The main oil field development strategy in Russia is waterflooding (water injection into an oil reservoir for oil displacement and pressure maintenance). Nowadays, most Russian oil fields are at the 3rd and 4th development stages, which means high water cut of produced liquid (90% and more). Main objec-tive of reservoir engineers is water production reduction (if possible) and oil production increase. These conditions require significant control and oil reservoir development regulation. To implement such activities, specialists need a solution for fast simulation of significant reservoirs and for fast evaluation of multiple development scenarios for testing hypothesis for geological struc-ture, history matching and production optimization. This approach is relevant for fractured porous res-ervoirs that have significant heterogeneity of filtration-conductivity properties. This fact causes early watercut growth in producing wells and leads to limitation of field production project targets. The paper proposes a methodology of numerical simulation of fractured porous oil reservoirs de-velopment based on the superelement conception. The model simulates two-phase filtration in a dual-porosity reservoir. A numerical scheme is fully explicit. The set of conservation equations is approxi-mated on a super-elements grid. This fact increases calculation speed and simplifies model generation (as cell dimension is consistent with well spacing). Calculation accuracy check requires production history matching. The proposed calculation methodology is tested on a real field example and checked by simulation in Rubis Kappa Engineering. Good matching results have been achieved at model training stages and forecast simulation.
Аннотация:Основной технологией разработки нефтяных месторождений в России является заводнение (закачка в нефтяной пласт воды для вытеснения нефти и поддержания пластового давления). При этом большая часть нефтяных месторождений нашей страны находится на 3-й или 4-й стадии раз-работки, что означает высокую обводненность добываемой жидкости (90 % и более). Основная задача специалистов по разработке нефтяных месторождений – уменьшить добычу воды и (по возможности) увеличить добычу нефти. В этих условиях много внимания уделяется контролю и регулированию разработки нефтяных месторождений. Для выполнения этих работ специалистам необходим инструмент, позволяющий быстро строить модели значительных по размерам месторождений и оперативно рассчитывать большое количество сценариев для проверки гипотез о геологическом строении, адаптации модели и решения задач оптимизации разработки. Наиболее актуально это для трещиновато-пористых коллекторов, поскольку они характеризуются высокой неоднородностью фильтрационно-емкостных свойств. Это провоцирует опережающее обводнение добывающих скважин, что препятствует достижению проектных показателей разработки месторождений. Предлагается методика численного математического моделирования разработки нефтяных месторождений в карбонатных трещиновато-поровых коллекторах на базе концепции суперэлементов. Фильтрация двухфазная, применена концепция двойной пористости. Численная схема полностью явная. Система уравнений сохранения аппроксимирована по пространству на суперэлементной сетке. Это позволяет существенно увеличить скорость вычислений и упростить построение моделей (так как размер ячеек сопоставим с расстоянием между скважинами). Для корректных расчетов необходима адаптация на историю разработки. Предлагаемая методика тестируется на модели реального месторождения, результаты расчетов сравниваются с расчетами на коммерческом симуляторе Rubis Kappa Engineering. Получены хорошее совпадение на этапе обучения модели и удовлетворительное совпадение результатов прогнозных расчетов.
Authors: I.V. Afanaskin (ivan@afanaskin.ru) - Federal State Institution "Scientific Research Institute for System Analysis of the Russian Academy of Sciences" (SRISA RAS) (Head of Group), Moscow, Russia, Ph.D, S.G. Volpin (sergvolpin@gmail.com) - Federal State Institution "Scientific Research Institute for System Analysis of the Russian Academy of Sciences" (SRISA RAS) (Head of Department), Moscow, Russia, Ph.D, A.V. Roditelev (avrod_94@mail.ru) - Federal State Institution "Scientific Research Institute for System Analysis of the Russian Academy of Sciences" (SRISA RAS) (Leading Programmer), Moscow, Russia, A.A. Kolevatov (akolevatov@niisi.ras.ru) - Federal State Institution "Scientific Research Institute for System Analysis of the Russian Academy of Sciences" (SRISA RAS) (Deputy Head of Department), Moscow, Russia, Ph.D
Keywords: oil field development simulation, two-phase filtration, carbonate reservoir, fractured porous reservoir, dual-porosity model, super-elements
Page views: 8396
PDF version article

Font size:       Font:

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

При задании фильтрационно-емкостных свойств и визуализации результатов моделирования трещиновато-пористых коллекторов обычно считается, что пласт содержит в два раза больше ячеек (2N). Тогда первая половина сетки (ячейки с 1-й по N-ю) описывает матрицу, а вторая (ячейки с N+1-й по 2N-ю) – трещины. Каждой ячейке матрицы должна соответствовать своя активная ячейка трещины и наоборот. Это двухкратное увеличение количества расчетных ячеек сильно замедляет расчеты.

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

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

Введем физико-математическую модель двойной пористости при следующих условиях:

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

-     жидкости слабосжимаемые, пласт упругий;

-     пласт состоит из двух вложенных друг в друга сред – пористых блоков (матрицы) и трещин; пористые блоки разбиты трещинами регулярным образом на параллелепипеды;

-     проницаемость трещин намного больше проницаемости матрицы; пористость (пустотность) трещин намного меньше пористости матрицы;

-     в скважине работают только трещины;

-     трещины обмениваются жидкостью с пористой матрицей, переток жидкости псевдоустановившийся;

-     фильтрацией жидкости внутри матрицы пренебрегаем.

Тогда можно записать две системы уравнений сохранения – отдельно для трещин и для матрицы. Эти системы будут связаны между собой источниковыми слагаемыми, отвечающими за переток фаз из одной среды в дру- гую [1].

Уравнения сохранения объемов нефти и воды в стандартных условиях в трещинах запи- шем в виде:

,          (1)

;        (2)

уравнения сохранения количества движения для трещин запишем в виде обобщенного закона Дарси:

,    (3)

,   (4)

где ff и kf – пористость и абсолютная проницаемость трещин;  и  – плотность перетока нефти и воды из матрицы в трещины; So,f  и Sw,f, Bo,f  и Bw,f,  и ,  и , kro,f и krw,f, Po,f  и Pw,f, ro,f  и rw,f  – насыщенность, объемный коэффициент, вектор скорости фильтрации, плотность источника (стока), моделирующего работу скважины, относительная фазовая проницаемость (ОФП), давление и плотность в пластовых условиях нефти и воды в трещинах; mo и mw – вязкость нефти и воды; g – ускорение свободного падения; Df – глубина залегания трещины по вертикали [1–3].

Система уравнений (1)–(4) дополняется замыкающими соотношениями следующим образом:

So,f  + Sw,f  = 1,                                        (5)

,                    (6)

, a = o, w,          (7)

, a = o, w,                         (8)

,                             (9)

где ff,0 – пористость трещин при начальном пластовом давлении Pf,0; Ba,0, a = o, w – объемный коэффициент фазы a при начальном пластовом давлении; Cr,f, Co,f  и Cw,f  – сжимаемость трещин, нефти и воды; ra, a = o, w – плотность фазы a; ra,STC, a = o, w – плотность фазы a в стандартных условиях; Pc,f  – капиллярное давление в трещинах в системе нефть–вода.

Систему уравнений (1)–(4) с учетом соотношений (5)–(9), принимая Sf  º Sw,f , Pf  º Po,f и используя подход, изложенный в [4], можно преобразовать в следующую систему уравнений:

     (10)

         (11)

,    (12)

.(13)

Система уравнений (10)–(13) дополняется начальными условиями следующим образом:

Pf = Pf (x, y, z, t = 0), Sf = Sf (x, y, z, t = 0)

и граничными условиями непротекания на внешних границах. Внутренние граничные условия на скважинах учитываются с помощью слагаемых  и .

Уравнения сохранения объемов нефти и воды в стандартных условиях в порах запишем в виде:

,                              (14)

,                              (15)

где обозначения аналогичны обозначениям, принятым для системы (1)–(2), с точностью до нижнего индекса m, обозначающего матричные блоки; уравнения сохранения количества движения для пор отсутствуют, так как фильтрацией жидкости в пористой матрице пренебрегаем [1–3].

Система уравнений (14)–(15) дополняется замыкающими соотношениями:

So,m + Sw,m  = 1,                                      (16)

,                  (17)

,                           (18)

где fm,0 – пористость матрицы при начальном пластовом давлении Pm,0; Cr,m – сжимаемость пор; Pc,m – капиллярное давление в порах в си- стеме нефть–вода.

Систему уравнений (14)–(15) с учетом соотношений (16)–(18), принимая Sm º Sw,m, Pm º Po,m и используя подход, изложенный в [4], можно преобразовать к следующей системе уравнений:

         (19)

. (20)

Система уравнений (19)–(20) дополняется начальными условиями следующим образом:

Pm = Pm(x, y, z, t = 0), Sm = Sm(x, y, z, t = 0)

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

Перетоки нефти и воды между трещинами и пористой матрицей учитываются с помощью слагаемых  и . Эти слагаемые связывают системы уравнений (10)–(13) и (19)–(20).

Численная схема

Аппроксимируем уравнения (10)–(13) и (19)–(20) на суперэлементной расчетной сет- ке [5, 6]. Размеры ячеек сопоставимы с расстоянием между скважинами (300–500 м). Сами ячейки строятся вокруг скважин. Тогда количество ячеек в одном горизонтальном слое примерно равно количеству скважин (иногда в больших неразбуренных областях вводятся фиктивные скважины для сгущения сетки). Форма ячеек соответствует сетке Вороного. Благодаря малому числу ячеек и их большому размеру такая вычислительная сетка позволяет использовать полностью явную расчетную схему и проводить расчеты с большой скоростью.

Суперэлементная сетка в горизонтальной плоскости близка к разбиению залежи на зоны дренирования по методу А.П. Крылова [7]. Выделяют четыре этапа построения суперэлементной сетки в плане:

-     отображение на карте границ залежи и скважин;

-     объединение внешних скважин и построение по ним выпуклой оболочки;

-     проведение триангуляции Делоне по скважинам;

-     отображение элементов Вороного, опуская перпендикуляры из середин ребер построенной триангуляции.

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

Рассмотрим двухмерную (для упрощения записи) полностью явную (по давлению и насыщенности) численную схему в горизонтальной плоскости для уравнений (10)–(13) и (19)–(20) и выразим неизвестные:

 (21)

 (22)

      (23)

      (24)

где нижний индекс i – номер ячейки; верхний индекс n – номер шага по времени; Dtn+1 – переменный шаг по времени (определяется исходя из ограничения на максимальное изменение давления и насыщенности за один временной шаг [2, 8]); Aij – длина ребра между соседними ячейками i и j; hi – толщина ячейки i; hij – средняя толщина соседних ячеек i и j (определяется как среднее арифметическое); Ni – количество соседних ячеек для ячейки i;  – средний объемный коэффициент фазы a = o, w между соседними ячейками i и j (определяется как среднее арифметическое);  – скорость фильтрации в трещинах фазы a = o, w между соседними ячейками i и j:

, a = o, w,

,

,

,

где Lij – расстояние между центрами двух соседних ячеек.

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

,

где ; ; , a = o, w;  – подвижность фазы a = o, w в трещинах;  – коэффициент продуктивности скважины; (Rc)i – радиус контура питания; (rw)i – радиус скважины;  – скин-фактор скважины.

После определения забойного давления  дебиты по фазам определяются как

, a = o, w.

Переток фаз между матрицей и трещинами , a = o, w, осуществляется в псевдоустановившемся режиме и определяется формулой (по аналогии с приводимыми в [1])

, a = o, w,

где , ,

;

,

,

 – матрично-трещинная проводимость; CD – коэффициент Дарси; si – фактор формы (шейп-фактор); lx, ly и lz – размеры блока матрицы;  – средний объемный коэффициент фазы a = o, w между матрицей и трещинами (определяется как среднее арифметическое).

В такой модели вода появится в продукции добывающей скважины сразу же после того, как текущая водонасыщенность ячейки превысит значение насыщенности связанной водой. В реальности этого не происходит. Воде не- обходимо время для того, чтобы пройти расстояние от нагнетательной скважины до до- бывающей. Для воссоздания этого эффекта используется ступенчатое задание функций относительной фазовой проницаемости: если водонасыщенность ячейки меньше насыщенности связанной водой, то ОФП по воде равна нулю, а ОФП по нефти равна своему максимальному значению; если водонасыщенность ячейки больше насыщенности связанной водой, то ОФП по нефти и воде рассчитываются как функции водонасыщенности [5, 6, 9].

Соотношения (21)–(24) позволяют рассчитывать давление и насыщенность в порах и трещинах. Они могут быть легко обобщены на трехмерный случай.

Моделирование реального месторождения

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

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

·     Средняя проницаемость трещин – 1 600 Д, пористость трещин – 0,04 д.ед.

·     Средняя проницаемость матрицы – 1,3 мД, пористость матрицы – 0,22 д.ед.

·     Эффективная нефтенасыщенная толщина – 8 м.

·     Начальное пластовое давление – 339 бар, начальная пластовая температура – 147 °С.

·     Сжимаемость трещин – 1,5∙10-4 1/бар, сжимаемость пористых блоков – 4,5∙10-5 1/бар.

·     PVT-свойства воды: объемный коэффициент – 1,02 м3/м3, вязкость – 0,36 мПа∙с, сжимаемость – 4,7∙10-5 1/бар, плотность в стандартных условиях – 1056 кг/м3.

·     PVT-свойства нефти: объемный коэффициент – 1,55 м3/м3, вязкость – 0,397 мПа∙с, сжимаемость – 1,01∙10-4 1/бар, плотность в стандартных условиях – 815 кг/м3, давление насыщения нефти газом – 188 бар, растворимость газа в нефти – 203 м3/м3.

ОФП и капиллярные давления в системе нефть–вода приведены на графиках (рис. 1).

Всего в разработке перебывало 127 скважин, из них 29 нагнетательных и 98 добывающих.

Для проверки применимости предлагаемой модели двойной пористости были созданы две гидродинамические модели: одна с помощью предлагаемой модели, другая – коммерческого симулятора Rubis компании Kappa Engineering [10]. Построенные вычислительные сетки приведены на рисунках (см. http://www.swsys. ru/uploaded/image/2019-3/2019-3-dop/1.jpg, http://www.swsys.ru/uploaded/image/2019-3/ 2019-3-dop/2.jpg). Для сгущения сетки при построении суперэлементной модели были использованы 30 фиктивных скважин (на ри- сунке они названы пустыми). История разра- ботки месторождения – 33 года. Из них 26 лет занял этап обучения суперэлементной модели, а последующие 7 лет – прогноз. Расчеты проводились при задаваемых на скважинах фактически замеренных дебитах добываемой жидкости и расходах закачиваемой воды. В процессе обучения модели уточнялись относительные фазовые проницаемости, проницаемость матрицы и пористость трещин. Суточная и накопленная добыча нефти и воды показаны на рисунках 2 и 3. Получено хорошее совмещение кривых, рассчитанных с помощью суперэлементной модели и Rubis, на периоде обучения суперэлементной модели и удовлетворительное совмещение на периоде прогноза.

Заключение

В работе описана модель двойной пористости для изучения разработки трещиновато-пористых коллекторов на базе концепции суперэлементов. Фильтрация двухфазная при забойном давлении в добывающих скважинах выше давления насыщения. Модель предполагает отсутствие фильтрации флюида по матрице. Приток жидкости в ствол скважины только из трещин, переток между матрицей и трещинами в псевдоустановившемся режиме. Размер ячеек суперэлементной сетки сопоставим с расстоянием между скважинами. Благодаря этому использована полностью явная численная схема (по давлению и насыщенности) с переменным шагом по времени. Большеразмерные ячейки позволяют использовать упрощенные геологические представления и быстро рассчитывать даже модели с очень большим количеством скважин. Это дает возможность применения предлагаемой модели для многовариантных расчетов и решения оптимизационных задач. Из-за больших размеров сетки обязательным этапом является адаптация модели на историю разработки или обучение модели с использованием результатов расчетов на более мелкой сетке. Предлагаемая математическая модель апробирована на реальном месторождении, результаты расчетов сравниваются с вычислениями с помощью коммерческого симулятора Rubis Kappa Engineering. Получены хорошее совпадение показателей разработки на этапе обучения модели и удовлетворительное совпадение на этапе прогноза.

Работа выполнена при поддержке РФФИ, грант № 16–29–15135 офи_м.

Литература

1.    Голф-Рахт Т.Д. Основы нефтепромысловой геологии и разработки трещиноватых коллекторов. М.: Недра, 1986. 608 с.

2.    Азиз Х., Сеттари Э. Математическое моделирование пластовых систем. М.–Ижевск: Изд-во ИКИ, 2004. 416 с.

3.    Каневская Р.Д. Математическое моделирование гидродинамических процессов разработки месторождений углеводородов. М.–Ижевск: Изд-во ИКИ, 2002. 140 с.

4.    Кац Р.М., Волгин Е.Р., Афанаскин И.В. Численное моделирование двухфазной фильтрации нефти и воды // Тр. НИИСИ РАН. 2014. Т. 4. № 2. С. 141–148.

5.    Булыгин Д.В., Мазо А.Б., Поташев К.А., Калинин Е.И. Геолого-технические аспекты суперэлементной фильтрационной модели нефтяных месторождений // Георесурсы. 2013. № 3. С. 31–35.

6.    Мазо А.Б., Булыгин Д.В. Суперэлементы. Новый подход к моделированию разработки нефтяных месторождений // Нефть. Газ. Новации. 2011. № 11. С. 6–8.

7.    Хисамутдинов Н.И., Хасанов М.М., Телин А.Г., Ибрагимов Г.З., Латыпов А.З., Потапов А.М. Разработка нефтяных месторождений: в 4 т.; Т. 1. Разработка нефтяных месторождений на поздней стадии. М.: Изд-во ВНИИОЭНГ, 1994. 240 с.

8.    Эртекин Т., Абу-Кассем Дж., Кинг Г. Основы прикладного моделирования. Ижевск: Изд-во ИКИ, 2012. 1060 с.

9.    Афанаскин И.В., Ялов П.В., Гиацинтов А.М., Родителев А.В. Решение задач оптимизации при суперэлементном моделировании разработки нефтяных месторождений // Программные продукты и системы. 2017. Т. 30. № 3. С. 384–391. DOI: 10.15827/0236-235X.119.384-391.

10.  Houze O., Viturat D., Fjaere O.S. Dynamic Data Analysis. V 5.20.01. Kappa Engineering Publ., 2018, 757 p.

References

  1. Golf-Raht T.D. Fundamentals of Oilfield Geology and Fractured Reservoir Development. Moscow, Nedra Publ., 1986, 608 p.
  2. Aziz Kh., Settari A. Petroleum Reservoir Simulation. London, Applied Science Publ., 1979, 476 p. (Russ. ed.: Moscow–Izhevsk, ICS Publ., 2004, 416 p.).
  3. Kanevskaya R.D. Mathematical Modeling of Hydrodynamic Processes of Hydrocarbon Deposit Development. Moscow–Izhevsk, ICS Publ., 2002, 140 p.
  4. Kats R.M., Volgin E.R., Afanaskin I.V. Numerical simulation of two-phase oil and water filtration. SRISA RAS Proc. 2014, vol. 4, iss. 2, pp. 141–148 (in Russ.).
  5. Bulygin D.V., Mazo A.B., Potashev K.A., Kalinin E.I. Geological and technical aspects of a superelement filtration model of oil fields. Georesources. 2013, iss. 3, pp. 31–35 (in Russ.).
  6. Mazo A.B., Bulygin D.V. Superelements. A new approach to modeling the development of oil fields. Oil. Gas. Novations. 2011, iss. 11, pp. 6–8 (in Russ.).
  7. Khisamutdinov N.I., Khasanov M.M., Telin A.G., Ibragimov G.Z., Latypov A.Z., Potapov A.M. Development of Oil Fields. Vol. 1. Development of Oil Fields at the Late Stage. Moscow, VNIIOENG Publ., 1994, 240 p.
  8. Ertekin T., Abu-Kassem J., King G. Basic Applied Reservoir Simulation. Izhevsk, ICS Publ., 2012,
    1060 p.
  9. Afanaskin I.V., Yalov P.V., Giatsintov A.M., Roditelev A.V. Optimization problems solution based on super-element modeling of oil-field development. Software & Systems. 2017, vol. 30, no. 3, pp. 384–391. DOI: 10.15827/0236-235X.119.384-391 (in Russ.).
  10. Houze O., Viturat D., Fjaere O.S. Dynamic Data Analysis. V 5.20.01. Kappa Engineering Publ., 2018, 757 p.

Permanent link:
http://swsys.ru/index.php?page=article&id=4625&lang=en
Print version
The article was published in issue no. № 3, 2019 [ pp. 478-485 ]

Perhaps, you might be interested in the following articles of similar topics: