Аналитическая модель притока жидкости для многоствольной скважины в анизотропном пласте | Вестник Томского государственного университета. Математика и механика. 2022. № 77. DOI: 10.17223/19988621/77/8

Аналитическая модель притока жидкости для многоствольной скважины в анизотропном пласте

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

An analytical model of the fluid influx for a multilateral well in an anisotropic formation.pdf Введение Истощение и ухудшение структуры запасов на месторождениях Западной Сибири требует вовлечения в разработку и эксплуатацию залежей углеводородов, представленных коллекторами с вторичными фильтрационно-емкостными свойствами и сложным и неоднородным строением. Вовлечение в разработку подобных коллекторов требует применения новых подходов и технологий, позволяющих увеличить добычу нефти путем увеличения охвата дренируемых запасов и одновременно сократить экономические издержки. Одной из таких технологий, успешно внедряемых в добывающих компаниях, является технология с применением многоствольных скважин (скважин со сложным профилем), которые в условиях анизотропии (неоднородности) подобных коллекторов экономически эффективнее, чем вертикальные скважины [1]. Разработка залежей с использованием скважин подобной конфигурации осложнена выбором оптимальной схемы их размещения и системы разработки для конкретного месторождения. В связи с этим нашей целью являлась выработка общего подхода к формированию оптимальной системы разработки месторождений нефти, включающей подбор схемы размещения и наиболее эффективных режимов эксплуатации скважин. Для решения данной проблемы целесообразно разработать методику аналитического расчета, позволяющую проводить экспресс-оценку притока жидкости к двуствольной скважине, не прибегая к многовариантным и затратным по времени расчетам с использованием гидродинамических моделей. Объект исследования На сегодняшний день фильтрация жидкости к вертикальной скважине достаточно полно изучена как российскими, так и зарубежными специалистами. Однако, как уже упоминалось выше, разработка и эксплуатация залежей, приуроченных к коллекторам с вторичными фильтрационно-емкостными свойствами, вертикальными скважинами неэффективна. Приток к скважинам сложного профиля (с несколькими стволами) остается предметом изучения современных исследователей. Одними из первых формулу притока жидкости для многоствольной скважины получили Ю.П. Борисов, В.П. Пилатовский и В.П. Табаков [2]: 102 Орлов А.А., Русских А.С. Аналитическая модель притока жидкости для многоствольной скважины Q 2пеЛР ln Rk - ln L • sin a 4 h , h • sin a +-ln- L • n 2nrc (1) _ k • h H где 0 - дебит жидкости многоствольной скважины, к - проницаемость пласта, h - толщина пласта, ДР - депрессия на пласт, д - вязкость жидкости, е - коэффициент гидропроводности пласта, Гс - радиус скважины, Rk - радиус контура питания, a - угол наклона горизонтального ствола от вертикали, L - длина стволов, n - количество стволов. В работе [3] задача расчета фильтрации жидкости к многоствольной скважине решена без учета анизотропии пласта в вертикальном направлении. По этой причине, как отмечают авторы методики, расчет может давать заниженные результаты по сравнению с фактическими данными о работе скважин. Одной из наиболее применяемых на практике аналитических моделей притока жидкости к многоствольной скважине является модель В.Г. Григу лецкого и Б.А. Никитина [2]: Q (2) 2тЛР Р •h Р•h , 4Rk L • n 2 nr ln-- c L где P - анизотропия проницаемости пласта. Методика [3] представляет собой частный случай решения [2] той же задачи, вследствие чего обе методики имеют существенный недостаток - не учитывают влияние области дренирования стволов друг на друга. Для расчета притока жидкости к многоствольной скважине также используется полуэмпирическая формула В.П. Меркулова [4]: Q = - 2тЛР nb* + h • (ln | ^+^ | + X) +1 ln I--r 1 2c J I n(a + b )8n-1 1- R„ 2n (3) где 5* - относительное смещение горизонтального ствола от среднего положения 1 л У**" в пласте, а значения параметров a , b , с , к определяются следующим образом: a* = l + 2h; b* = 2lh + 4h2; c* = l; X = 0,426/ h + 4,45 . (4) Недостатком методики [3] является то, что она полуэмпирическая, следовательно, не обладает необходимой универсальностью. При выборе конфигурации многоствольной скважины основными факторами являются длина стволов и угол между ними. Расчет максимальных дебита и продуктивности двуствольной скважины затрудняет интерференция, связанная с взаимным расположением стволов скважины. Выбор оптимального варианта расположения стволов скважины путем многовариантных расчетов с использованием гидродинамической модели занимает значительное количество времени. В известных аналитических методиках [2-4] не рассмотрено влияние угла между стволами, а также длины стволов на интерференцию (взаимовлияние), и, следовательно, на продуктивность многоствольной скважины. 103 Механика / Mechanics Аналитическая методика расчета притока жидкости к двуствольной скважине Рассмотрим приток к двуствольной скважине, у которой оба ствола расположены в одной плоскости, на одном уровне и имеют длину L. Для оценки притока область вокруг скважины схематично разделена на две зоны: зона I - внешняя часть двуствольной скважины, зона II - ее внутренняя часть, для которой характерна интерференция стволов (рис. 1). Рис. 1. Геометрическое представление двуствольной скважины Fig. 1. Geometry of a dual-lateral well Каждая из зон фильтрации имеет соответствующие граничные условия (5). Зона I Зона II О II о II Р | = Рс = const * Р = Р Р = Р 1 c ,1 \\ Ф=0=2ж 1 c (5) О? II йн Р = Р r=rt Рк где Pc - давление в скважине, Pk - давление на контуре питания, п - радиус контура питания, ф - угол между стволами. Уравнение, описывающее приток жидкости Q к двуствольной скважине, является решением системы уравнений (5) [5, 6] (6) _k_ ЭР ц dr б = ^ I s где v - скорость фильтрации, S - площадь фильтрации. Решение системы уравнений (6) для описания притока жидкости к двуствольной скважине имеет вид: 104 Орлов А.А., Русских А.С. Аналитическая модель притока жидкости для многоствольной скважины Q = 4kLAP < ln| ^ Л-' 2 г. b l ln 1 Рһ 'j ' 1 2 Р JJ Д-ОРҺ , B = 2b 2 1 + sin a 2r 2kh3AP exp(B) al 2|л 1 2/, l Br„ ln l (7) (8) где a - большая полуось эллипса дренирования, b - малая полуось эллипса дренирования (рис. 2), rc - радиус скважины. Рис. 2. Схема области дренирования стволов скважины Fig. 2. Scheme of a drainage area of wellbores Депрессия на пласт АР вычисляется как разница между пластовым давлением и давлением в скважине [7, 8]: AP = Рпл - Pc , (10) Рпл =Р8 l# + L sin 0j • (11) где p - плотность жидкости; g - ускорение свободного падения; H - глубина залегания пласта; Ѳ - угол наклона горизонтального ствола к вертикали пласта; L -эффективная мощность продуктивного пласта. Учет взаимовлияния стволов друг на друга в модели описан коэффициентами. А и В (8), которые учитывают угол между стволами и радиус их дренирования соответственно. Полученные результаты и их анализ Для расчета притока были использованы данные, представленные в таблице 1. Верификация предложенной модели проведена с использованием фактической геолого-геофизической информации, характерной для Неокомских продуктивных 105 Механика / Mechanics отложений Северо-Островного месторождения Ханты-Мансийского автономного округа (см. табл. 1); она показала, что результаты проведенных нами расчетов коррелируют с результатами расчетов, полученных с помощью методик [2] и [3], использующих другие подходы к описанию притока жидкости к многоствольной скважине (рис. 3). Таким образом, можно сделать вывод, что наша модель адекватно описывает приток жидкости к двуствольной скважине и может быть использована для дальнейших расчетов. Таблица 1 Исходные данные для расчета Параметр Размерность Значение Глубина залегания пласта м 2200 Мощность пласта м 10 Проницаемость пласта в латеральном направлении пласта мД 14.2 Проницаемость пласта в вертикальном направлении пласта мД 7.1 Вязкость жидкости сПз 1.3 Радиус скважины м 0.1 Давление на забое скважины МПа 15 Новая методика---Методика [2] -------Методика [3] Рис. 3. Зависимость дебита от длины стволов двуствольной скважины при расчете разными методами Fig. 3. Well flow rate as a function of the length of bores for a dual-lateral well calculated by different methods 106 Орлов А.А., Русских А.С. Аналитическая модель притока жидкости для многоствольной скважины С использованием предложенной методики проведена оценка влияния изменения длины стволов скважины на ее дебит жидкости (угол между стволами принят в расчетах 45° - наиболее распространенное значение при бурении двуствольных скважин) [9]. В результате проведенных расчетов определена эффективная длина стволов скважины, которая составила 750 м (см. рис. 3). При увеличении длины ствола до 1 000 м прирост дебита жидкости составляет менее 3%, а затраты на бурение увеличиваются на 17%. Следовательно, использование стволов длиной более 750 м экономически неэффективно. После выбора длины стволов скважины исследовано влияние угла между стволами на ее продуктивность. Результаты расчетов с учетом критерия минимального прироста дебита жидкости (< 3%) представлены в табл. 2. Рекомендуемый угол между стволами скважины составил 50°. При его увеличении до 180° прирост дебита жидкости незначительный (табл. 3) - около 1%. Т аблица 2 Зависимость дебита жидкости от изменения длины стволов Длина стволов скважины, м Дебит жидкости, м3/сут Прирост дебита жидкости, % 15 197 - 25 270 37.1 50 347 28.5 75 388 11.8 100 415 7.0 150 453 9.2 200 472 4.2 300 491 4.0 400 508 3.5 500 524 3.1 600 541 3.2 700 558 3.1 750 575 3.0 800 583 1.4 900 590 1.2 1000 596 1.0 Таблица 3 Зависимость дебита жидкости от углов между стволами Угол между стволами, град Дебит жидкости, м3/сут Прирост дебита жидкости, % 0 166 - 10 184 10.8 20 202 9.8 30 225 11.4 35 246 9.3 40 300 22.0 45 431 43.7 50 510 18.3 55 523 2.5 60 528 1.0 107 Механика / Mechanics Окончание табл. 3 Угол между стволами, град Дебит жидкости, м3/сут Прирост дебита жидкости, % 80 538 1.9 100 549 2.0 120 558 1.6 140 567 1.6 160 572 0.9 180 578 1.0 Таким образом, наибольшей эффективностью обладает двуствольная скважина с длиной стволов 750 м и углом между стволами 50°. Заключение Разработана методика аналитического расчета притока жидкости к двуствольной скважине, учитывающая угол между стволами многоствольной скважины и интерференцию между ними, позволяющая проводить экспресс-оценку притока жидкости к двуствольной скважине, которая отличается от известных методик тем, что позволяет без использования многовариантных расчетов на гидродинамическом симуляторе определить с параметры скважины, которые обеспечат ее максимальную продуктивность и экономическую эффективность. На основании численных исследований с использованием разработанной нами аналитической модели определена эффективная конфигурация двуствольной скважины с максимальной эффективностью, которая учитывает интерференцию между стволами скважины: длина стволов скважины 750 м, угол между ними 50°.

Ключевые слова

фильтрация жидкости, многоствольная скважина, приток жидкости, контур питания, аналитическая модель

Авторы

ФИООрганизацияДополнительноE-mail
Орлов Алексей АлексеевичТомский политехнический университетдоктор технических наук, профессор Отделения ядернотопливного цикла Инженерной школы ядерных технологийorlovaa@tpu.ru
Русских Анастасия СергеевнаООО «Недра-Консалт»начальник отдела гидродинамического моделированияanastasiya_russkih@bk.ru
Всего: 2

Ссылки

Гусейнова Е.Л., Гусейнов Э.М. Технология бурения многоствольных скважин // Совре менные технологии в нефтегазовом деле - 2017 : сб. тр. Междунар. науч.-техн. конф. : в 2 т. Уфа : Изд-во УГНТУ, 2017. Т. 1. С. 239-242.
Борисов Ю.П., Пилатовский В.П., Табаков В.П. Разработка нефтяных месторождений горизонтальными и многозабойными скважинами. М. : Недра, 1964. 156 с.
Никитин Б.А., Григулецкий В.Г. Стационарный приток нефти к одиночной горизон тальной скважине в анизотропном пласте // Нефтяное хозяйство. 1992. № 8. С. 1012.
Меркулов В.П. О дебите наклонных и горизонтальных скважин // Нефтяное хозяйство. 1960. № 2. С. 51-56.
Басниев К.С. Подземная гидромеханика : учебник для вузов. Ижевск : Ин-т компьютер ных исслед., 2006. 488 с.
Гиматудинов Ш.К. Физика нефтяного и газового пласта : учебник для вузов. М. : Альянс, 2005. 311 с.
Баренблатт Г.И. Движение жидкостей и газов в природных пластах. М. : Недра, 1984. 208 с.
Маскет М. Течение однородных жидкостей в пористой среде. М. : Гостоптехиздат, 1949. 590 с.
Фрайя Х, Омер Э, Пулик Т., Джардон М., Кайя М, Паэс Р, Сотомаор П.Г., Умуджоро К. Новые подходы к строительству многоствольных горизонтальных скважин // Нефтегазовое обозрение. 2003. Весна. С. 44-67.
 Аналитическая модель притока жидкости для многоствольной скважины в анизотропном пласте | Вестник Томского государственного университета. Математика и механика. 2022. № 77. DOI: 10.17223/19988621/77/8

Аналитическая модель притока жидкости для многоствольной скважины в анизотропном пласте | Вестник Томского государственного университета. Математика и механика. 2022. № 77. DOI: 10.17223/19988621/77/8