Модель двойной пористости для изучения разработки трещиновато-пористых коллекторов на базе концепции суперэлементов Текст научной статьи по специальности «Науки о Земле и смежные экологические науки»

УДК 622.276.1/4 Дата подачи статьи: 18.06.19

DOI: 10.15827/0236-235X.127.478-485 2019. Т. 32. № 3. С. 478-485

Модель двойной пористости для изучения разработки трещиновато-пористъх коллекторов на базе концепции суперэлементов

И.В. Афанаскин 1, к.т.н, ведущий научный сотрудник, ivan@pfanaskin.ru С.Г. Вольпин 1, к.тн, зав. отделом, sergvolpin@gmail.com А.В. Родителев 1, ведущий программист, rav@niisi.ras.ru А.А. Колеватов 1, к.тн, зам. зав. отделом, akolevatov@niisi.rcis.ru

1 Федеральный научный центр Научно-исследовательский институт системных исследований Российской академии наук, г. Москва, 117218, Россия

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

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

Предлагается методика численного математического моделирования разработки нефтяных месторождений в карбонатных трещиновато-поровых коллекторах на базе концепции суперэлементов. Фильтрация двухфазная, применена концепция двойной пористости. Численная схема полностью явная. Система уравнений сохранения аппроксимирована по пространству на суперэлементной сетке. Это позволяет существенно увеличить скорость вычислений и упростить построение моделей (так как размер ячеек сопоставим с расстоянием между скважинами). Для корректных расчетов необходима адаптация на историю разработки.

Предлагаемая методика тестируется на модели реального месторождения, результаты расчетов сравниваются с расчетами на коммерческом симуляторе Rubis Kappa Engineering. Получены хорошее совпадение на этапе обучения модели и удовлетворительное совпадение результатов прогнозных расчетов.

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

д_

dt

д_ dt

ГФ А, fл

Bo

Ф А, /

B

o, f

il. О л

,f

J

i ^ л

Wo, f

\ f

V

f ^ л

Ww, f B ,

lo,mf ■.

mf 5

(1)

(2)

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

Г =- — №•/ В; )] (3)

цо

кгк

Ww,f =-

f rw, f

M w

[v(Pw.f )-Pw,fgV(Df )] (4)

где ф/ и к/ - пористость и абсолютная проницаемость трещин; дот/ и - плотность перетока нефти и воды из матрицы в трещины; Бо,/ ^ ^ _ _

и Бщ/, Во,/ и В/ Ж»,/ и Ж.,/ , ц0 и , кго,/ и

к™,/, Ро,/ и Р/ Ро/ и рщ/ - насыщенность, объемный коэффициент, вектор скорости фильтрации, плотность источника (стока), моделирующего работу скважины, относительная фазовая проницаемость (ОФП), давление и плотность в пластовых условиях нефти и воды в трещинах; цо и ц» - вязкость нефти и воды; g - ускорение свободного падения; Б/ - глубина залегания трещины по вертикали [1-3].

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

Б0,/ + Бщ/ = 1, (5)

Ф/ =Ф/,о [1 + С,/ (Р/ -/)] , (6)

Ва = Ва,о [ 1 - Са(Р - Ро)], а = о, (7)

Ра=Ра^тс/Ва , а = 0, (8)

Ро./ - К/ = Р./ (Я..,), (9)

где ф/,0 - пористость трещин при начальном пластовом давлении Р/0; Ва,0, а = о, » - объемный коэффициент фазы а при начальном пластовом давлении; Сг,/, Со,/ и Сщ/ - сжимаемость трещин, нефти и воды; ра, а = о, » - плотность фазы а; ра,хтс, а = о, » - плотность фазы а в стандартных условиях; Рс,/ - капиллярное давление в трещинах в системе нефть-вода.

Систему уравнений (1)-(4) с учетом соотношений (5)-(9), принимая Б/ = Б»,/, Р/ = Ро,/ и используя подход, изложенный в [4], можно

преобразовать в следующую систему уравнений:

-Р

Ф / ,0 [ + (С, - С) + св +

+во 0 V

i — Л

Wo,f

в f

V °,f У

+ B „V

w ,0

-

Ww,f

в f

V У

(10)

= - в

Oo [?o - ? o,mf I w ,0

Kf ,0

dSf , x dPf

—f Crf + CK ) Sf—^

дt y ,f ' f dt

+B„ ,0V

i — Л Ww, f

B.

(11)

Wo,f =-

f У

fro,f ^ o

= -B„,0 [?

[V(Pf )-Po,f8ViDf )] (12)

Ww.f = -

[v( pf )-v( ^ )-pw,fgv( df )] .(b)

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

Pf = Pf (x, у, t = 0), = Sf (x, у, t = 0) и граничными условиями непротекания на внешних границах. Внутренние граничные условия на скважинах учитываются с помощью слагаемых ао/ и .

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

д_

dt

d dt

ф S

im o,.

л

= - 1o.

mf

ч s ^

т m w,m

',mf 5

(14)

(15)

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

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

Ь'о,т + Sw,m = 1, (16)

Фт =Фт,0 [ 1 + Cr m (Pm — Pm,0 )] ,

P — P = P ( S ),

o,m w,m c,m y w,m ) '

(17)

(18)

где фот,0 — пористость матрицы при начальном пластовом давлении Рт, 0; Сг,т — сжимаемость

пор; Р с, т — капиллярное давление в порах в системе нефть—вода.

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

дР

ф „[с +(С - С ) 5 + С 1—^ =

"и,0 г ,т \ ^ о / т о J ^^

(19)

= — в ей , — в

o0 1 o,mf

Ф

т m

',mf w,0 У w,mf 5

f+(C- +C- ) s. f

-— B„ f„r. (20)

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

Рт = Рт(х, у, г, t = 0), Sm = Sm(x, у, г, t = 0) и граничными условиями непротекания на внешних границах.

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

о, т/ ™,т/

вают системы уравнений (10)—(13) и (19)—(20).

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

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

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

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

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

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

- отображение элементов Вороного, опус-

—

кая перпендикуляры из середин ребер построенной триангуляции.

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

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

( /=( /

Bx,» l( 1а

4a,mf ,

-At"

N

'+1

1=1

A h {W: n ц

'i ч b v а j 1 ] и

(21)

W„,

fh (ф f,» )i |c„f+( c, - co )( sf )'+ co

i \"+1 i ( f =( Sf)r

D i N

-Аг17ШЧ i( * )"-( q-mf )"+1

Fh (Фf,» Ц ' j=i

-(C,f + C )(Sf )" [(Pf )"+1 -(Pf )" ],

( Pm )i ( Pm )j

-At"+1

(22)

Bo,» ( Qo,mf )" + К» (<w )" (23)

Fh M ;[CCr,m +(Cw - Co) ( Sm)" + Co ]'

( Sm Г =( Sm )" -At'

+ 1 Bw,»

Fh (Ф.» ),(^)Г (24)

(Cr,m + Cw )(Sm )' [(Pm )' (Pm )'

где нижний индекс I - номер ячейки; верхний индекс п - номер шага по времени; Д?п+1 - переменный шаг по времени (определяется исходя из ограничения на максимальное изменение давления и насыщенности за один временной шаг [2, 8]); Лгу - длина ребра между соседними ячейками г иу; Ъг - толщина ячейки г; Ъц - средняя толщина соседних ячеек г и у (определяется как среднее арифметическое); N - количество соседних ячеек для ячейки г;

(5а )] - средний объемный коэффициент фазы

а = о, » между соседними ячейками г и у (определяется как среднее арифметическое); (Жа)]

- скорость фильтрации в трещинах фазы а = о, » между соседними ячейками г и у:

ку (О] (Фа,/£-(Фа,/)]

(Wx)" =

, а = o, w,

Ма

2 (kf )' ( kf )

k = ' ' 1 (kf ), +(kf )j ' \ n

( КХ =

(к™,/ ); . (Фа,/ ); ^(Фа,/ )] . (к™,/ )] . (Фа,/ )] <(Фа,/ )] .

(фо,/ / =(Р )] -(Ро,/ / (В/ ), ,

(Ф/ )] =(Р )] -(Р/ )] -(Р.,/ )> (), ,

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

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

( Pw*u )Ч=( Pf )"-(q )"/\Wi: |_(4f )"+(* w^ )" 2я( kf )ihi

где WI. =

ln

m (^ ),

■; (r ) = »,14^2f ;

- Skin"

( k V

i \n V rx,f I. , \n

(Xx,f ), = -, X = 0, W; (^а,f )- подвиж-

m)",

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

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

( о;=win / )П [( pf )П - ( pell )i П, а=^^ ^^

Переток фаз между матрицей и трещинами (Qamf ). , а = o, w, осуществляется в псевдоуста-

новившемся режиме и определяется формулой (по аналогии с приводимыми в [1])

( <W ). [(Фа m ).-(Фа / ). Fh ,

а = o, w,

L

'}

a=o

где

( ^ ) =

(Tmf ) = Cd (Kf ) • 2 ( kf ) ( К ) ,

ст. = 4

1

l2 l2 l2

( kf )+( km ) '

(k™,mf ). = '

(Ка.,т )] , (Фа,т )" ЦФа,/ )" , (кга,/ )] , (фа,т )" <(фа,/ )] , (фо,т)! =( -(Ро,, (»т\ ,

(Ф м!,т I =( Рт )] -(Рс,т )] ^т )] 8 (^ ), ,

(Тт{) — матрично-трещинная проводимость;

Сэ — коэффициент Дарси; а, — фактор формы (шейп-фактор); ¡х, 1у и ¡г — размеры блока матрицы; ()" — средний объемный коэффициент фазы а = о, ^ между матрицей и трещинами (определяется как среднее арифметическое).

В такой модели вода появится в продукции добывающей скважины сразу же после того, как текущая водонасыщенность ячейки превысит значение насыщенности связанной водой. В реальности этого не происходит. Воде необходимо время для того, чтобы пройти расстояние от нагнетательной скважины до добывающей. Для воссоздания этого эффекта используется ступенчатое задание функций относительной фазовой проницаемости: если во-донасыщенность ячейки меньше насыщенности связанной водой, то ОФП по воде равна нулю, а ОФП по нефти равна своему максимальному значению; если водонасыщенность ячейки больше насыщенности связанной водой, то ОФП по нефти и воде рассчитываются как функции водонасыщенности [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/бар.

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

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

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

1,0

0,9 0,8 0,7 0,6

■i 0,5 ф

=i 0,4

1= 0,3 ©

О 0,2 0,1 0,0

4

1,2 1,0 0,8 0,6 0,4 0,2 0,0

0,1 0,2 0,3 0,4 0,5 0,6 0,7 0,8 0,9 1,0 Водонасыщенность, д.ед.

а)

krw

- •-• - - -

8,0 7,2 6,4 5,6 4,8 4,0 3,2 2,4 1,6 0,8 0,0

0,0 0,1 0,2 0,3 0,4 0,5 0,6 0,7 0,8 0,9 1,0 =

J

Водонасыщенность, д.ед. б;

Рис. 1. ОФП для нефти (kro) и воды (krw), а также капиллярное давление (Pc) в системе нефть-вода: а) для трещин, б) для пористой матрицы

Fig. 1. Relative permeability for oil (kro) and water (krw), as well as capillary pressure (Pc) in the oil-water system: а) for fractures, б) for a porous matrix

Всего в разработке перебывало 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, на периоде обучения суперэлементной модели и удовлетворительное совмещение на периоде прогноза.

2000

1800

1600

£ 1400

1200

,ы

д

о в 1000

и

и 800

тф

е н 600

а

ч ы 400

б

оД 200

0

_ X II

fl ' t

ЛЛ

IV \

/

i И V

4.» / \ v

Г ! i3 kp и H M

100

90

80

70

. n

60 u ■

,ь

50 ст о

н

40 не н

д

30 о в

б

О

20

10

0

0 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33

Время, года

-Добыча нефти - Суперэлементы Добыча воды - Rubis

Добыча нефти - Rubis -Обводненность - Суперэлементы

Добыча воды - Суперэлементы Обводненность - Rubis

Рис. 2. Суточная добыча нефти и воды Fig. 2. Daily production of oil and water

9000 8000 7000 6000 5000 4000 3000 2000 1000 0

—- Накопленная добыча нефти - Суперэлементы А Накопленная добыча нефти - Rubis Накопленная добыча воды - Суперэлементы • Накопленная добыча воды - Rubis —

ij

IfT

0 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33

Время, года

Рис. 3. Накопленная добыча нефти и воды Fig. 3. Cumulative oil and water production

Заключение

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

пользовать упрощенные геологические представления и быстро рассчитывать даже модели с очень большим количеством скважин. Это дает возможность применения предлагаемой модели для многовариантных расчетов и решения оптимизационных задач. Из-за больших размеров сетки обязательным этапом является адаптация модели на историю разработки или обучение модели с использованием результатов расчетов на более мелкой сетке. Предлагаемая математическая модель апробирована на реальном месторождении, результаты расчетов сравниваются с вычислениями с помощью коммерческого симу-лятора 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.

Software & Systems Received 18.06.19

DOI: 10.15827/0236-235X.127.478-485 2019, vol. 32, no. 3, pp. 478-485

Dual porosity model for fractured porous reservoirs development analysis based

on the superelement concept

I.V. Afanaskin 1, Ph.D. (Engineering), Leading Researcher, ivan@afanaskin.ru S.G. Volpin 1, Ph.D. (Engineering), Head of Department, sergvolpin@gmail.com A.V. Roditelev 1, Leading Programmer, rav@niisi.ras.ru

A.A. Kolevatov 1, Ph.D. (Engineering), Deputy Head of Department, akolevatov@niisi.ras.ru

1 Federal State Institution "Scientific Research Institute for System Analysis of the Russian Academy of Sciences" (SRISA RAS), Moscow, 117218, Russian Federation

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 objective 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 structure, history matching and production optimization. This approach is relevant for fractured porous reservoirs 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 development 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 approximated 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.

Keywords: oil field development simulation, two-phase filtration, carbonate reservoir, fractured porous reservoir, dual-porosity filtration mechanism, super-elements.

Acknowledgements. The paper has been supported by the program of Russian Foundation for Basic Research no. 16-29-15135 ofim.

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.

Для цитирования

Афанаскин И.В., Вольпин С.Г., Родителев А.В., Колеватов А.А. Модель двойной пористости для изучения разработки трещиновато-пористых коллекторов на базе концепции суперэлементов // Программные продукты и системы. 2019. Т. 32. № 3. С. 478-485. DOI: 10.15827/0236-235X. 127.478-485.

For citation

Afanaskin I.V., Volpin S.G., Roditelev A.V., Kolevatov A.A. Dual porosity model for fractured porous reservoirs development analysis based on the superelement concept. Software & Systems. 2019, vol. 32, no. 3, pp. 478-485 (in Russ.). DOI: 10.15827/0236-235X.127.478-485.