Актуальные проблемы нефти и газа ■ Вып. 2(21) 2018 ■ http://oilgasjournal.ru
ВЫЧИСЛЕНИЕ КОЭФФИЦИЕНТА ПРОВОДИМОСТИ СКВАЖИННОГО СОЕДИНЕНИЯ - ПОЛУАНАЛИТИЧЕСКИЙ МЕТОД СТЕНФОРДСКОГО УНИВЕРСИТЕТА
Э.С. Закиров, С.Н. Закиров, И.М. Индрупский, Д.П. Аникеев ИПНГ РАН, e-mail: [email protected]
В данной статье обозревается полуаналитический метод Стенфордского университета для вычисления коэффициента проводимости скважинного соединения скважины произвольной траектории. Она является продолжением статей* данного номера журнала. В некотором смысле, рассматриваемая статья завершает обсуждение методов вычисления исследуемой величины для прогрессивных скважин, ставя «жирную точку» в этом вопросе. Всем гидродинамикам, осуществляющим моделирование скважин сложной траектории, рекомендуется использовать данный метод для вычисления коэффициента проводимости.
Описываемый ниже подход является плодом труда целого коллектива исследователей в Стенфордском университете (в том числе магистрантов, докторантов). Можно отметить работы [1, 2, 3, 4, 5], составившие основу рассматриваемого метода. На сегодняшний день предлагаемый подход реализован в программном продукте AdWell, используемом рядом компаний [6]. Программа разработана и поддерживается Dr. Huanquan Pan [7].
Данный метод разработан для точной оценки проводимости ячейки, вскрытой скважиной нетрадиционного типа в сложных сеточных областях. Основной недостаток -использование метода предполагает режим истощения пластовой энергии, но вполне можно считать, что данное условие не является сильно ограничивающим.
Метод основан на решении однофазной фильтрационной задачи в локальной околоскважинной области. Обычно эта локальная область представляет собой зону дренирования скважины. Из-за особенностей полуаналитического решения данная
* См. статьи:
Закиров Э.С., Закиров С.Н., Индрупский И.М., Аникеев Д.П. «О представлении скважины в 3D гидродинамической модели»;
Закиров Э.С., Закиров С.Н., Индрупский И.М., Аникеев Д.П. «Вычисление коэффициента проводимости скважинного соединения. Метод Писмена и близкие подходы»;
Закиров Э.С., Закиров С.Н., Индрупский И.М., Аникеев Д.П. «Вычисление коэффициента проводимости скважинного соединения наклонной или горизонтальной скважины»;
Закиров Э.С., Закиров С.Н., Индрупский И.М., Аникеев Д.П. «Вычисление коэффициента проводимости скважинного соединения на сетке Вороного» в данном выпуске.
область должна иметь форму прямоугольного параллелепипеда.
Отличительной особенностью данной задачи является ее существенная зависимость от условий, заданных на скважине. Граничные условия могут представлять собой или постоянный заданный потенциал, или условие непротекания. Как оказалось, вычисляемая проводимость скважинной ячейки слабо зависит от типа применяемых граничных условий. С вычислительной точки зрения наиболее удобным оказалось задание суммарного дебита скважины и условий непротекания на всех внешних границах, за исключением нижней грани модели, на которой задается постоянное давление. Поскольку в этом случае скважина выходит на стационарный режим и удается избежать нестационарных вычислений.
При решении задачи вычисления коэффициента проводимости скважинного соединения дискретизации подвергается траектория скважины, а не пласт. При этом ствол скважины разбивается на сегменты, что снижает потребности в вычислительных ресурсах.
Контрольное решение, включающее забойное давление, определяется полуаналитически, основываясь на функциях Грина. Как только соответствующее решение найдено, специальный алгоритм определения пересечений скважины с сеточными блоками идентифицирует сеточные блоки, вскрытые скважиной, и осуществляет распределение дебита скважины по сеточным интервалам.
Это необходимый этап, поскольку сегмент может не совпадать в точности с одной сеточной ячейкой. Затем используется традиционное конечно-разностное моделирование в той же самой околоскважинной области с теми же граничными условиями. При этом для неоднородных коллекторов мелкомасштабная проницаемость должна быть промасштабирована на целевую сетку. Отличие численной задачи от аналитической состоит в том, что источники и стоки определены при помощи полуаналитического решения. Соответственно, совместное использование сеточных давлений (как результатов конечно-разностных вычислений) и дебитов вместе с забойным давлением (из полуаналитического контрольного решения) позволяет вычислять точные проводимости скважинной ячейки при помощи уравнения
w
Т = д Д (1)
Рг - Рг
Здесь д* - объемный расход интервала г -й скважины, ^ - вязкость флюида, pi -давление в сеточном блоке, р* - давление в скважине.
Данный метод отличается наибольшей общностью и повышенной функциональностью. Он позволяет точно моделировать очень сложные (например, многозабойные) скважины или даже группы скважин, которые пересекают сетку произвольным образом. Может использоваться сетка любого типа, если только на ней можно проводить однофазные вычисления конечно-разностным методом.
Кроме того, метод вычисления проводимости скважинной ячейки основывается на вычислениях по определенной пользователем локальной околоскважинной области, а не по отдельным сеточным блокам, профильным разрезам или бесконечным областям. В этом он выгодно отличается от предыдущих подходов. При этом данный способ не ограничен контрольным решением. Однако, это аналитическое решение доказало свою максимальную точность для однородных полей и приемлемую точность для случаев неоднородных полей проницаемости. Этот метод учитывает сферичность притока к концам скважины, а также взаимодействие с другими скважинами и границами.
Значительно позднее в работе [8] был представлен близкий полуаналитический подход, основанный на функциях Грина и его третьем тождестве.
Итак, полуаналитический метод вычисления проводимости скважинной ячейки основан на использовании полученных полуаналитическим способом следующих переменных: забойного давления и распределения притока к скважине, а также пластовых давлений, полученных численным конечно-разностным моделированием. По определению, проводимость скважинной ячейки может быть вычислена напрямую с использованием формулы (1).
Первый шаг в данном подходе состоит в точном вычислении притоков к скважине и забойных давлений для каждого сегмента ствола скважины. При помощи полуаналитического метода, основанного на функциях Грина, для решения однофазной фильтрационной задачи:
/ \ дФ
У-(кУФ ) = т/ис1— . (2)
Уравнение (2) соответствует однофазному течению слабосжимаемого флюида. Условные обозначения общеизвестны. Данное уравнение записывается относительно потенциала Ф = р + для учета эффекта действия гравитационных сил. Дополнительно к уравнению (2) необходимо задать соответствующие начальные и граничные условия.
Для получения искомого решения только скважина дискретизируется на сегменты,
без дискретизации пласта. После того, как решение получено для каждого сегмента, контрольные дебиты каждого блока, вскрытого скважиной, определяются по специальному алгоритму детекции пересечений блоков со скважиной. По этим дебитам вычисляют пластовые давления с использованием однофазного конечно-разностного симулятора (рис. 1).
Контрольное решение (полуаналитическое)
ч!' р!
(7 - сегмент скважины)
Отображение на конечно-разностную сетку
, р7 ^ ч!, Р!
ч!, р! ^ рЬ'оск
(I - скважинный блок)
Рис. 1. Аналитическое решение (приток к скважине и давление) и его связь с конечно-разностной сеткой (давления в блоках)
Наконец, полуаналитическую проводимость скважинной ячейки получают комбинацией дебита скважины и забойного давления из полуаналитического решения и соответствующего давления в сеточном блоке из однофазного конечно-разностного моделирования при помощи формулы (1). При нахождении единственного значения Ж1 для каждого блока дебиты и давления должны быть получены для стационарных или псевдостационарных условий.
Для включения в модель околоскважинной неоднородности в работе [9] предложен подход, называемый «£ - к * моделью» [10]. В рамках данной модели неоднородность проницаемости представляется при помощи постоянной фоновой проницаемости к* и локального скин-эффекта £. Фоновую проницаемость получают степенным осреднением или любым другим методом ирБсаН^'а [11]. Например, решением линейной фильтрационной задачи на полной геологической сетке. Скин-эффект изменяется вдоль траектории скважины. Околоскважинная проницаемость измененной зоны ка . для
каждого сегмента у вычисляется через интегрирование с весами значений мелкомасштабной проницаемости вдоль сегмента по области с эффективным радиусом г .. Применение формулы Хоукинса для скин-фактора [12] позволяет вычислить значение
скин-фактора £. для каждого сегмента:
£ =
( 1* Л к -1
V ка, 7 У
1п ^ . (3)
Тогда значение скин-эффекта включается в систему линейных уравнений, описывающую приток в индивидуальные сегменты скважины. Использование модели £ - к * незначительно усложняет базовый полуаналитический подход.
Публикации [9, 13, 14] подтверждают приемлемую точность результатов применения описанной методики для ряда тестовых примеров на основании их сравнения с контрольными конечно-разностными вычислениями.
Все сказанное позволяет точно моделировать скважины с различными траекториями, пересекающими сеточную область произвольным образом. Рассматриваемый метод учитывает сферичность притока к концам скважины, а также взаимодействие с другими скважинами и границами.
Сама вычислительная часть реализована в симуляторе AdWell в рамках проекта БиРЫ-НМ в Стенфордском университете.
От словесного описания на уровне понимания идеи метода перейдем к деталям ее реализации на уровне формул.
Предполагается, что проницаемость к является диагональным тензором. Тогда уравнение (2) может быть записано в следующем виде:
, д2Ф , д2Ф , д2Ф дФ (л,
х у г а
Скважины могут иметь произвольную конфигурацию и траекторию (рис. 2). Каждая скважина из общего числа скважин п^ может иметь щ (¡м!) боковых стволов, каждый из которых состоит из щ(¡м!) сегментов. Общее число сегментов равно N :
щ )
N. = Е 5>.(К,Ь) . (5)
Г
Ъ-М 1
Рис. 2. Произвольная конфигурация скважины (из работы [15])
Перемасштабированием координатных осей отротропно-анизотропную среду можно заменить изотропной. Соответствующие формулы преобразования и пересчета радиуса скважины можно найти в работе [16]. При этом скважина из круга превращается в эллипс. В новой системе координат уравнение (4) записывается следующим образом:
д2Ф д2Ф д2Ф дФ
- +
- +
У2 /2 7'2
(6)
к дг
Затем осуществляется обезразмеривание всех переменных уравнения (6), когда все размерные переменные превращаются в безразмерные величины в диапазоне от 0 до 1. Обезразмеривание приводит к переменным гв, чв ,ФВ, хп, уа, га и уравнению:
V 2ф = Ф .
Совместная модель пласт-скважина имеет 2Ы,, неизвестных. Они суть следующие: безразмерные дебиты притока (оттока) ) и депрессия (^, / ) в средней
точке М(гм!,/,^) каждого сегмента скважины ¡п (см. рис. 2). Следовательно, для корректного описания системы необходимо такое же число уравнений. Эти уравнения
получают из условия баланса массы, соотношений для депрессии с учетом гидравлической характеристики течения на забое.
Закон сохранения массы требует, чтобы расход каждой скважины равнялся сумме вкладов (расходов) всех сегментов скважины:
X X Чо , ь,)= Чо )
1 ь 1
(7)
где ^ = 1,...,^ . В соответствии с принципом суперпозиции значение депрессии в любой
точке пространства может быть выражено в виде суммы вкладов от всех источников/стоков в системе. При рассмотрении депрессии в каждом сегменте системы
принцип суперпозиции дает:
) п*(] м Л1 )
Фо Ом ЛI, Л ) = X X X Чо О* Л1, Л Рп Ом , 7, Л )[М, Ь , )] ,
(8)
1м= 1 11 =1 =1
где Ф01,7,7 )] представляет собой безразмерную депрессию в середине
сегмента 1 бокового ствола 7 скважины 1, вызванную источником/стоком единичной
интенсивности сегмента ^ бокового ствола /, скважины .
Безразмерная величина депрессии может быть вычислена при помощи функций мгновенного источника/стока. Например, для одномерной системы с непроницаемыми границами в безразмерных координатах 0 и 1 безразмерный перепад давления в точке х в момент времени t, вызванный мгновенным источником единичной интенсивности в точке X в начальный момент времени, может быть вычислен по формуле:
ад
(х, х0, ^ ) = 2X ехр (-к ) Бт кжх соб кжх0
к=1
(9)
■Е
к=-ад
ехр
( х 2 к х0) 4 ^
2Л
+ ехр
(х - 2 к + х0) 4 ^
2 V
В [15] показывается, что оба представления формулы (9) эквивалентны. Оба выражения требуют суммирования бесконечных рядов. Однако, экспоненциальная форма сходится быстрее для малых t. Тогда как тригонометрический ряд имеет преимущество для больших t. Используя разделение переменных, решение трехмерной задачи может быть получено умножением одномерных решений типа (9). Интегрирование по времени решения трехмерной задачи дает непрерывное решение:
1
tD (
Ф.„
(jw, jl, )[м (К> ь, Ь)] =-1 £
XDeУDeZDe 0
XD XD 0 tD
(
V XDe XDe XDe У
УD УD 0 tD
V УDe УDe УDe У
(10)
zD zD0
V ZDe ZDe ZDe У
Наконец, интегрирование точечных источников/стоков по скважинному сегменту [М0, М1 ] дает выражение для значения депрессии в середине каждого сегмента:
МБ
Ф (М», , jl , )[М (^ , Ь, )]= |
Мп
ФwD (jw , М , js )[М (iw , Ь, К )]"
4жг
wD
ШГ) +
(11)
Мт 1
Г-.
4жг
Мв0 ' wD
Член добавлен и вычтен одновременно во избежание проблем со
сходимостью при интегрировании (см. работу [17]), поскольку вблизи ствола скважины (на расстоянии радиуса скважины) решение сингулярно. Действительно, около точечного источника течение сферично, а давление изменяется обратно пропорционально расстоянию.
Потенциал представляет собой решение уравнения:
4л '
,.2 Ф
Згп
(12)
Тем самым прием с вычитанием члена позволяет снизить сингулярность
подынтегрального выражения первого интеграла. Второй интеграл можно вычислить аналитически. Допустим, скважина представлена сегментом [Мдо, Мд х], а вклад
интеграла оценивается в точке М. Обозначим угол Моо-М,15 МтМ) символом р, а длину [МД0М] символом X. Тогда
MD1 1
I = Г —
* Аггг
МП„ 4^ 4лЬ
1п
1 - Х соб В + Ь
X2 2 X
Ь ь
соб р +1
X
Ь
(1 - соб В)
(13)
Если сегмент скважины представлен не прямой линией, то рассматриваемый интеграл можно вычислить разбиением кривой на части.
К данному моменту большинство вычислительных шагов прояснено с достаточной степенью подробности. За дальнейшими нюансами использования метода
заинтересованный читатель отсылается к основополагающим работам, отвечающим на вопросы: как учесть гидравлические характеристики течения в стволе скважины, как решить полученную систему уравнений, как вычислить скин-фактор, как произвести учет забойных устройств, управляющих притоком и т.д.
Казалось бы, использование столь продвинутой вычислительной техники является крайне далеким вопросом при моделировании скважин нетрадиционных типов. Однако, как показывают значительное количество примеров работ [1, 2, 3, 4, 5], достаточно импортировать рассчитанные указанным сложным способом проводимости скважинных ячеек в коммерческий симулятор, чтобы сохранить точность аналитического моделирования в обычном гидродинамическом симуляторе.
Поэтому рекомендуется использовать AdWell в качестве предпроцессора подготовки проводимостей скважинных ячеек.
Статья написана в рамках выполнения государственного задания (тема «Научное обоснование новых экологически чистых технологий разработки месторождений углеводородов в сложных горно-геологических условиях на основе SD-компьютерных экспериментов», № АААА-А16-116022510270-1).
ЛИТЕРАТУРА
1. Fukagawa H. Semianalytical Modeling of Nonconventional Well Performance in Realistic Reservoirs // Master's report. Stanford University, 2002.
2. Serve J. An Enhanced Framework for Modeling Complex Well Configurations // Master's Report. Stanford University, 2002.
3. Valvatne P.H. A Framework for Modeling Complex Well Configurations // M.S. Thesis. Stanford University, 2000.
4. Wolfsteiner C., Aziz K., Durlofsky L.J. Modeling Conventional and Non-Conventional Wells // Paper presented at 6th International Forum on Reservoir Simulation. Salzburg, Austria. 3-7 September 2001. 11 р.
5. Wolfsteiner C., Durlofsky L.J., Aziz K. Calculation of Well Index for Nonconventional Wells on Arbitrary Grids // Comp. Geosciences. 2003. Vol. 7. Р. 61-82.
6. Durlofsky L.J., Aziz K. Advanced Techniques for Reservoir Simulation and Modeling of Nonconventional Wells // Final Report. 2004. August. 213 р.
7. Pan H.Q. AdWell 2.1 User Manual. Stanford University, 2005.
8. Kurtoglu B., Medeiros F.Jr., Ozkan E., Kazemi H. Semianalytical Representation of Wells and Near-Well Flow Convergence in Numerical Reservoir Simulation // Paper SPE 116136 prepared for presentation at the SPE Annual Technical Conference and Exhibition. Denver, Colorado, USA, 21-24 September 2008. 20 р.
9. Wolfsteiner C., Durlofsky L.J., Aziz K. Approximate Model for Productivity of Nonconventional Wells in Heterogeneous Reservoirs // SPEJ. 2000. June. Р. 218-226.
10. Durlofsky L.J. An Approximate Model for Well Productivity in Heterogeneous Porous Media // Mathematical Geology. 2000. Vol. 32. Р. 421-438.
11. Закиров Э.С. Upscaling в 3D компьютерном моделировании. М.: ЗАО «Книга и Бизнес», 2007. 344 c.
12. Hawkins M.F. A Note on the Skin Effect // Trans. AIME. 1956. Vol. 207. Р. 356-357.
13. Valvatne P.H., Durlofsky L.J., Aziz K. Semi-analytical modeling of the performance of intelligent well completions // Paper SPE 66368 prepared for presentation at the SPE Reservoir Simulation Symposium. Houston, Texas, 11-14 February 2001. 10 р.
14. Wolfsteiner C., Durlofsky L.J., Aziz K. Efficient estimation of the effects of wellbore hydraulics and reservoir heterogeneity on the productivity of nonconventional wells // Paper SPE 59399 prepared for presentation at the SPE Asia Pacific Conference. Yokohama, Japan, 25-26 April 2000. 10 р.
15. Ouyang L.-B., Aziz K. A General Single-Phase Wellbore/Reservoir Coupling Model for Multilateral Wells // SPERE&E. 2001. Vol. 4. Р. 327-335.
16. Besson J. Performance of Slanted and Horizontal Wells on an Anisotropic Medium // Paper SPE 20965 prepared for presentation at the European Petroleum Conference. The Hague, Netherlands, 22-24 October 1990. 14 р.
17. Maizeret P.D. Well Indices for Nonconventional Wells // Master's report. Stanford University, 1996.