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

УДК 622.276.5.001.5 Дата подачи статьи: 26.04.18

Б01: 10.15827/0236-235Х.031.3.500-506 2018. Т. 31. № 3. С. 500-506

Интерпретация, межскважинных исследований карбонатных коллекторов методом двух режимов с помощью численных моделей

И.В. Афанаскин 1, к.т.н.., ведущий научный сотрудник, ivan@qfanaskin.ru С.Г. Вольпин 1, к.т.н.., зав. отделом, sergvolpin@gmail.com О.В. Ломакина 1, научный сотрудник, ovllom@mail.ru Ю.М. Штейнберг 1, зав. группой, yurash22@gm.ail.com

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

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

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

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

Ключевые слова: гидродинамические исследования скважин, метод двух режимов, гидропрослушивание, модель двойной пористости.

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

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

Метод двух режимов, используемый для исследования скважин [3, 10-14], обычно применяется и

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

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

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

Решение прямой задачи

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

скважинного пространства для любой геометрии пласта, расстановки и количества скважин, распределения свойств пласта по площади. Такая модель может быть только численной. Рассмотрим границы применимости описываемой математической модели фильтрации флюида в продуктивном пласте:

- фильтрация однофазная;

- жидкость слабосжимаемая;

- пласт упругий;

- гравитационными силами можно пренебречь;

- продуктивный пласт содержит поры и трещины;

- проницаемость трещин значительно больше проницаемости пор;

- пористость трещин значительно меньше пористости пор;

- участки пористой матрицы не обмениваются между собой флюидом;

- к скважине приток происходит только по трещинам;

- при снижении давления поры матрицы отдают флюид в трещины;

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

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

Тогда уравнения сохранения объема флюида в поверхностных условиях отдельно для трещин и

пор можно записать соответственно как

(

д_ dt

% ^ B

+ div

W_

B

v у

Я well + qmf ■

(1)

~4mf .

(2)

5/ Фт дг [ В

где нижние индексыf и m означают трещины и пористую матрицу соответственно; ф - пористость;

B - объемный коэффициент флюида; W - вектор скорости фильтрации; - плотность потока в скважину; - плотность перетока из пор матрицы в трещины [1, 2, 4, 15-17].

Скорость фильтрации определяется через закон Дарси:

]¥ = -—grad (р), (3)

где kf - проницаемость трещин; ц - вязкость флюида; Pf - давление в трещинах.

Уравнения (1)-(3) дополняются начальными условиями, граничными условиями и замыкающими соотношениями:

Ф / =Ф/0 [ 1 + (Р/ - Р )] , (4)

Фт = ФтО [ 1 + C,m (Рт " P )] ,

в = в0 [ 1 - C (Pf -р )].

(5)

(6)

где фр, фто и Bo - пористость трещин, пористость матрицы и объемный коэффициент флюида соответственно при начальном давлении P0; Crf, Crm и Cl - сжимаемость трещин, пор и флюида соответственно.

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

Я/ = Тт/ (Рт - РГ ) ,

(7)

ГДе Tmf =

к

mf

VBmf (

к^ =

mf

к к

fx

к„ + к

fx

kmykfy kmy + kfy

кк

fz

к„, + к

fz .

^ . Bm = 2 [ B ( Pf ) + B ( Pm )]

Т^ - коэффициент трещинно-поровой проводимости; с - фактор формы (шейп-фактор); Vb - элементарный объем пласта; ктх, кту, ^ и к^, kfy, kfz -проницаемость пор и трещин вдоль осей X, У, Ъ декартовой системы координат; п - параметр формы элемента (блока) пористой матрицы (п = 1 -плита, п = 2 - призма или цилиндр, п = 3 - куб или сфера); х2т - характерный размер блоков матрицы [1, 2, 4, 15-17].

Приток флюида (источниковое слагаемое) к вертикальной скважине в псевдоустановившемся приближении [18] запишется как

= Р1 (РГ - ) , (8)

где PI =

2 пк-h

eff eff

1

( Pf )

ln

R

кeff •Jkfickjy

+ Skin

heff - эффективная толщина пласта; Rc - радиус контура питания; rwell - радиус скважины; Skin - скин-фактор скважины.

Вычислим радиус контура питания для прямоугольной конечно-разностной вычислительной сетки:

' ] IV - " ]

где i и ] - индексы ячейки сетки вдоль осей X и У соответственно; Дх и Ду - размеры ячейки вдоль осей X и У соответственно.

Уравнения (1) и (2) с учетом (4)-(6) можно переписать в виде

, ,dPf Фf» (Cf + Ci) = ~Bo 4wdi -q„f + div

i

W

B

v у

(9)

дР -

Фто(С. + С)-. =-Б0дт/ . (10)

дt

Запишем конечно-разностный аналог уравнения (9) для трещиноватой среды. Фильтрация двумерная, в горизонтальной плоскости. Численная схема - «классики». Неявное уравнение для давления в трещинах:

Ч 4У ч (Ф , 0),, (с, + С )[( Р V-(р I

= Ain+1 B \PI-j

(P j -(Pwl)

-(/ (( Рт V/] + (11)

^ (р+ц - р+1 (V - р-И)+

+о,5 (р- рг1 -о,5 (V1 - V1)},

где 4f+l - переменный шаг по времени, определяемый в соответствии с рекомендациями [18] так, чтобы максимальное по всем блокам сетки изменение давления не превышало заданного значения (при расчетах отдельно задавалось предельное изменение давления в матрице и в трещинах),

.„ ()/+о,5V 4У; 4-;

л, _

i+0,5 j

( f )

M"B/4-0,5 j Axi+0,5 j

_ 2 (kX Lj (kX )j

( f )+1 j +(kx ) =={b [( Pf )n

= =( v +xij )

i+1j

B

(Pf y. ]} ■

i + 0,5 j

аналогично для ^0,5., ц+0>5 и xn_о,5 •

Явное уравнение для давления в трещинах:

: At B„ {PI.

)ij

ч 4у ч (ф/о )„ (с, + с )[(р,; -р);

[( р, )Г-( р- ) -(/ (№/]+ (12)

+Л+о,5; (р+1; - р )-^о,5; (^ - ) +

+л;+о,5 (р+х - р И-о,5 (р; - ;)}.

Название «классики» для разностной схемы связано со способом обхода расчетных точек [19]. На каждом шаге по времени п + 1 обход точек пространственной сетки в одной (в нашем случае горизонтальной) плоскости совершается дважды. На первом и последующих шагах по времени с нечетным п + 1 вычисляются значения давления в ячейках с нечетной суммой i + ]. Этот первый обход осу-

ществляется с помощью явной по давлению схемы. При втором обходе на том же шаге по времени вычисляют значения давления в ячейках с четной суммой i + ] с помощью неявной по давлению схемы, в соседних ячейках берутся значения давления, полученные на том же шаге по времени при первом обходе. На втором и последующих шагах по времени с четными п + 1 роли узлов с четной и нечетной суммой i + ] меняются. Кратко это можно резюмировать так: при i + ] + п + 1 четном берется уравнение (12), а при нечетном - уравнение (11).

Запишем конечно-разностный аналог уравнения (10) для давления в пористой матрице:

^4У;4-, (Фто ), (Ст + С )[(рт )" -(рт % ] =

~ -4^+х Во (, [(р. V,]. <13)

Из (13) можно легко выразить давление в матрице на новом временном шаге (р );+х через дав-

).+1 ij

i \п

ление в трещинах ( Pf )

= Ax.jAyAzj (Ф„0 )j К + Ci )(P ).

m V AXjAyAzj (ф„0 )j (Crm + Ci ) + Atn+1B0 (Tmf ) А^+Ч (Tmf )j ( Pf )

)n+1 ij

. (14)

Axv Ay yAz i:j (фт0) ij(Crm + С,) + Atn+1B0 (Tmf )П

Подставляя (14) в (11) и (12), исключаем из них неизвестное давление в матрице на новом шаге по

времени (р ^ .

При расчете перетоков (7) для плоского случая 2Б-геометрии будем считать = max {kfx,^ j,

kmz =am^VkmXkmr, Где amz - коэффиЦиент вертикальной анизотропии проницаемости пористой матрицы.

Решение обратной задачи

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

Пусть необходимо найти минимум функции многих переменных fX), где X = (xi, x2, ..., xn). В данном случае функция fX) - это невязка между расчетными и фактическими точками давления, замеренного на забое скважины. Эта задача эквива-

лентна задаче нахождения значений X, при которых градиент функцииfX) равен нулю:

gradfX)) = 0. (15)

Применим к (15) метод Ньютона: grad(fXj)) + H(Xj)(Xj+1 -Xj) = 0, (16)

где j = 1, 2, 3, ..., m - номер итерации; H(X) - гессиан функции fX).

Напомним, что гессиан функции - это симметричная квадратичная форма, описывающая поведение функции во втором порядке:

H ( Х ) = Ê ÊW;,

где a,j = d2f/dxidxf; функцияf(X) задана на и-мерном пространстве вещественных чисел.

В более удобном для вычислений виде формулу (16) можно представить как

1 = _ H _1 ( ) grad ( f ( )).

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

Пример интерпретации гидродинамических исследований межскважинного пространства трещиновато-пористого коллектора методом двух режимов с помощью описанной модели

Для проверки описанной модели в качестве фактических данных будем использовать значения давления, полученные путем численного моделирования на сетке Вороного с помощью коммерческой программы для интерпретации гидродинамических исследований скважин и пластов Saphir компании Kappa Engineering [5].

Геометрия рассматриваемой модели и расстановка скважин приведены на рисунке 1. Модель содержит три скважины. Две добывающие (возмущающие) скважины по очереди исследуются методом двух режимов (путем однократного изменения режима работы, то есть дебита [3, 10-14]). В третьей (реагирующей) скважине фиксируется изменение забойного давления в результате работы первых двух скважин. Все время исследования составляет 30 суток. Одна возмущающая скважина эксплуатируется 10 суток с дебитом 300 м3/сут., затем 20 суток с дебитом 100 м3/сут. Другая возмущающая скважина эксплуатируется 20 суток с дебитом 300 м3/сут. , затем 10 суток с дебитом 100 м3/сут.

Основные исходные данные для расчетов следующие:

_ радиус скважины - 0,1 м; _ толщина пласта - 10 м;

Y

M

11,3 км

Возмущающая скважина ^

у

500 м

< ^

\

Реагирующая скважина

Рис. 1. Геометрия рассматриваемой модели и расстановка скважин

Fig. 1. Considered model geometry and well array

- расстояние между возмущающими и реагирующей скважинами - 500 м;

- объемный коэффициент - 1,1 м3/м3;

- вязкость - 1,3 мПас;

- общая сжимаемость - 4,3 • 10-5 см2/кгс;

- начальное пластовое давление - 250 кгс/см2.

Тензор проницаемости трещиноватой породы

диагональный. Проницаемость трещинной части пласта зависит от направления фильтрации. Оси расчетной сетки сориентированы по осям тензора проницаемости. Сетка неравномерная, размеры ячеек по оси Х приведены на рисунке 2. По оси У размеры ячеек аналогичны. Реагирующая скважина находится в ячейке с координатами 120-120, которая имеет размеры 0,1x0,1 м. Поскольку по формуле (8) невозможно определить забойное давление при нулевом дебите, в качестве забойного давления в реагирующей скважине принято давление в ячейке, которую она занимает. Поэтому приняты столь малые размеры ячейки.

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

k„

kf

( HQ ) f

где 1 - пара-

(Иф С,) , +(м?с, I

метр, характеризующий перетоки между матрицей и трещинами; ю - отношение емкости трещин к общей емкости системы; И - толщина пласта; С -суммарная сжимаемость пласта и флюида [1, 2, 4, 15-17].

При этом из параметра перетоков X обычно определяют проницаемость матрицы кт либо фак-

100

: 10

0,1

—

9

—•

—

-Л—

-—•--— •___

9

• •

9 • • •

•к

10 20 30 40 50 60 70 80 90 100 110 120 130 140 150 160 170 180 190 200

Номер ячейки по оси Х Рис. 2. Размеры ячеек вдоль оси X Fig. 2. Cells dimensions along the X axis

0

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

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

1. Проницаемость трещин к/, пористость матрицы фт, отношение емкостей ю, параметр перетоков X, анизотропия проницаемости трещин по площади (кх/ку)/.

2. Отношение емкостей ю, параметр перетоков X.

3. Проницаемость трещин к/, пористость матрицы т, анизотропия проницаемости трещин по площади (кх/ку)/.

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

Результаты интерпретации приведены в таблице. Видно, что наилучшее качество уточнения

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

В качестве примера на рисунке 3 приведены результаты интерпретации для первого начального приближения.

Заключение

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

Рассмотрено несколько вариантов интерпретации для различных наборов уточняемых пара-

Начальное приближение и результаты интерпретации (уточненное приближение) Initial guess and interpretation results (refined value)

Параметр Факт Нач. прибл. 1 Уточн. нач. прибл. 1 Нач. прибл. 2 Уточн. нач. прибл. 2 Нач. прибл. 3 Уточн. нач. прибл. 3

kf, мД 350 150 211 350 350 150 375

ф m, д.ед. 0,1 02 0,08 0,1 0,1 02 0,09

ю, д.ед. 0,1 0,3 0,25 0,3 0,09 0,1 0,1

X, б/р 110-6 1-10-5 110-5 110-5 1,3 10-6 110-6 110-6

(kx/ky)f, д.ед. 0,1 0,5 0,46 0,1 0,1 0,5 0,16

Примечание: подчеркиванием выделены уточняемые параметры.

Ч в a

и

• Давление, факт

• Давление, начальное приближение 1

• Давление, уточнение приближения 1 Дебит возмущающей скважины 1 Дебит возмущающей скважины 2

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

Время, сутки

Рис. 3. Результаты интерпретации для первого начального приближения Fig. 3. Interpretation results for the first initial guess

650

1500

350

1200

050

50

300

25

50

0

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

Поэтому в качестве определяемых параметров рекомендуется выбирать проницаемость трещин, пористость матрицы и анизотропию проницаемости трещин по площади. Использовать в качестве определяемых параметров сразу все пять не рекомендуется.

Работа выполнена при поддержке Программы фундаментальных научных исследований государственных академий наук № 1.2.П27, заказ № 18-0111П, НИР № 0065-2018-0111.

Литература

1. Деева Т.А., Камартдинов М.Р., Кулагина Т.Е., Манга-зеев П.В. Гидродинамические исследования скважин: анализ и интерпретация данных. Томск: Изд-во ЦППС НД ТПУ, 2009. 243 с.

2. Кульпин Л.Г., Мясников Ю.А. Гидродинамические методы исследования нефтегазоводоносных пластов. М.: Недра, 1974. 200 с.

3. Эрлогер Р. мл. Гидродинамические методы исследования скважин. М.-Ижевск: Изд-во ИКИ, 2007. 512 с.

4. Bourdet D. Well Test analysis: the use of advanced interpretation models. Amsterdam: Elsevier Sc. B.V., 2002, 425 p.

5. Houze O., Viturat D., Fjaere O.S. Dynamic data analysis. V 5.12. Kappa Engineering, 2017, 743 p.

6. Vela S., McKinley R.M. How areal heterogeneities affects pulse-test results. Soc. Pet. Eng. J., 1970, no. 2, pp. 181-191.

7. Jahns Hans O. Rapid method for obtaining a two-dimen-

sional reservoir description from well pressure response data. Soc. Pet. Eng. J., 1966, no. 4, pp. 315-327.

8. McKinley R.M., Vela S., Carlton L.A. A field application of pulse-testing for detailed reservoir description. JPT, 1968, vol. 20, no. 3, pp. 313 - 321.

9. Woods E.G. Pulse-test response of a two-zone reservoir. Soc. Pet. Eng. J., 1970, vol. 10, no. 3, pp. 245-256.

10. Афанаскин И.В., Вольпин С.Г., Ломакина О.В., Штейн-берг Ю.М. Гидродинамические исследования вертикальных скважин в пластах с двумя границами методом двух режимов // Вестн. кибернетики. 2017. № 4. С. 45-54.

11. Афанаскин И.В., Вольпин С.Г., Штейнберг Ю.М., Ломакина О.В. Гидродинамические исследования горизонтальных скважин методом двух режимов // Вестн. кибернетики. 2017. № 2. С. 108-113.

12. Афанаскин И.В., Крыганов П.В., Вольпин С.Г., Штейнберг Ю.М., Вольпин И.А. Изучение свойств нефтяных пластов с помощью гидродинамических исследований скважин методом двух режимов - теория, моделирование и практика // Вестн. ки-бернетики.2015.№ 3. С. 86-107.

13. Афанаскин И.В., Крыганов П.В., Вольпин С.Г., Штейнберг Ю.М., Вольпин И.А. Оценка фильтрационных и энергетических параметров нефтяных пластов с помощью гидродинамических исследований скважин на двух режимах: теория, моделирование и практика // Тр. НИИСИ РАН. 2015. Т. 5. № 1. С. 41-55.

14. Russell D.G. Determination of formation characteristics from two rate flow test. JPT, 1963, vol. 15, pp. 1317-1355.

15. Бурде Д. Интерпретация результатов исследований скважин // Petroleum Engineering and Related Management Tran-ing Gubkin Academy: матер. лекций. М., 1994. 109 с. (рус.).

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

17. Райс Л. Основы разработки трещиноватых коллекторов. М.-Ижевск: Из-во ИКИ, 2012. 118 с.

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

19. Роуч П. Вычислительная гидродинамика. М.: Мир, 1980. 616 с.

Software & Systems

DOI: 10.1582l/0236-235X.031.3.500-506

Received 26.04.18 2018, vol. 31, no. 3, pp. 500-506

Carbonate reservoirs crosswell survey interpretation by a two-rate test using numerical models

I.V. Afanaskin ', Ph.D., Leading Researcher, ivan@afanaskin.ru S.G. Volpin ', Ph.D., Head of Department, sergvolpin@gmail.com O.V. Lomakina ', Researcher, ovllom@mail.ru Yu.M. Shteynberg ', Head of Group, yurash22@gmail.com

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

Abstract. The paper describes a two-rate well test method that is applied at transient flow. Due to recent development this method allows defining a set of parameters that is the as wide-spread pressure build up test. It is established that the method investigates few dozen meters area, rarely hundreds. The two-rate method doesn't require well shut-in, so oil output loss reduce during the research. A disadvantage of the method is that parameter calculation is not so precise.

The paper shows that at significant rate difference between two modes and sufficient second-rate longevity the two-rate method is a reasonable alternative to interference test. In this case, test area between wells increases to a hundreds of meters. The paper describes a dual-porosity mathematic model for crosswell carbonate fractured reservoir well test interpretation. This model is applicable for crosswell reservoir two-rate test data interpretation.

There is a "classic" original finite-difference calculation scheme for this model. The authors describe a reverse-problem subsurface hydrodynamic solution by Newton method. The method has been applied to synthetic downhole pressure curve with corresponding recommendations. Few interpretation variants presented for different reservoir parameters refinement. A relative fractured volume and matrix-fractures diffusivity coefficient have low influence on a tested well at late time pressure value. It is recommended to include fractures permeability, matrix porosity and area fractures anisotropy as defined parameters. Keywords: well test, two-rate test, interference-test, dual-porosity model.

Acknowledgements. The article has been supported by the Program of Fundamental Scientific Research of State Academies of Sciences no. I.2.P27, Order no. 18-0111P, Research no. 0065-2018-0111.

1. Deeva T.A., Kamartdinov M.R., Kulagina T.E., Mangazeev P.V. Well Test: Data Analysis and Interpretation. Tomsk, TsPPS ND TPU Publ., 2009, 243 p.

2. Kulpin L.G., Myasnikov Yu.A. Investigation of Oil and Gas Reservoirs by Hydrodynamic Methods. Moscow, Nedra Publ., 1974, 200 p.

3. Erloger R. jr. Advances in Well Test Analysis. Moscow-Izhevsk, IKI Publ., 2007, 512 p.

4. Bourdet D. Well Test Analysis: The Use of Advanced Interpretation Models. Amsterdam, Elsevier Science B.V. Publ., 2002, 425 p.

5. Houze O., Viturat D., Fjaere O.S. Dynamic Data Analysis. V 5.12. Kappa Engineering Publ., 2017, 743 p.

6. Vela Saul, McKinley R.M. How areal heterogeneities affects pulse-test results. Soc. Pet. Eng. J. 1970, pp. 181-191.

7. Jahns Hans O. Rapid method for obtaining a two-dimensional reservoir description from well pressure response data. Soc. Pet. Eng. J. 1966, pp. 315-327.

8. McKinley R.M., Vela Saul, Carlton L.A. A field application of pulse-testing for detailed reservoir description. J. Pet. Tech. 1968, pp. 313-321.

9. Woods E.G. Pulse-test response of a two-zone reservoir. Soc. Pet. Eng. J. 1970, pp. 245-256.

10. Afanaskin I.V., Volpin S.G., Lomakina O.V., Shteynberg Yu.M. Two-rate tests of vertical wells in formations with two boundaries. Proc. in Cybernetics. 2017, no. 4, pp. 45-54 (in Russ.).

11. Afanaskin I.V., Volpin S.G., Shteynberg Yu.M., Lomakina O.V. Two-rate tests of horizontal wells. Proc. in Cybernetics. 2017, no. 2, pp. 108-113 (in Russ.).

12. Afanaskin I.V., Kryganov P.V., Volpin S.G., Shteynberg Yu.M., Volpin I.A. Oil-bearing formations characterization by two-rate well tests results - theory, modeling and practice. Proc. in Cybernetics. 2015, no. 3, pp. 86-107 (in Russ.).

13. Afanaskin I.V., Kryganov P.V., Volpin S.G., Shteynberg Yu.M., Volpin I.A. Evaluation of oil-bearing formation flow and energy parameters by two-rate well tests results - theory, modeling and practice. SRISA RAS Proc. 2015, vol. 5, no. 1, pp. 41-55 (in Russ.).

14. Russell D.G. Determination of formation characteristics from two rate flow test. J. Pet. Tech. 963, pp. 1347-1355.

15. Burde D. Interpretation ofwell test. Proc. Petroleum Engineering and Related Management Traning Gubkin Academy. Moscow, 1994, 109 p.

16. Golf-Rakht T.D. Fundamental of Fractured Reservoir Engineering. Moscow, Nedra Publ., 1986, 608 p.

17. Rays L. Fundamentals of the Fractured Reservoirs Development. Moscow-Izhevsk, IKI Publ., 2012, 118 p.

18. Aziz Kh., Settari E. Petroleum Reservoir Simulation. Moscow-Izhevsk, IKI Publ., 2004, 416 p.

19. Rouch P. Computational Fluid Dynamics. Moscow, Mir Publ., 1980, 616 p.

References