Comparison of various approaches in numerical modeling of well stimulation methods

Cover Page

Cite item

Abstract

This paper a numerical simulation of methods of intensification of flow to the well by using various approaches. To simulate the hydrochloric acid impact, an approach was applied based on changing the well productivity factor, as well as an approach that consisted of using a chemical reaction in a hydrodynamic model. The criteria for qualitative forecasting of technological development indicators was the actual data on one of the wells of the analogous field of the object. As a result of calculations on the example of a actual field in the conditions of extended horizontal wells, increments of additional oil production were obtained with applying various approaches of modeling the process of hydrochloric acid treatment. It was revealed that in the conditions of extended horizontal wells, the use of negative values of skin factors increases oil production by several times in comparison with the approach of compositional modeling with the occurrence of chemical reactions. Sensitivity analysis to the volume and concentration of the injected acid was carried out by using specialized software. It was concluded that as a result of the chemical reaction in compositional modeling, the effect of hydrochloric acid treatment significantly depends on the composition of the rock, the reaction rate, the concentration and volume of the injected acid. Multi-stage hydraulic fracturing was modeled by using a planar fracture system tool and discrete fracture system model. A slight discrepancy between the results of calculating the hydrodynamic model between these methods for modeling multi-stage hydraulic fracturing was identified.

Full Text

Введение

Солянокислотная обработка (далее – СКО) скважин предназначена для воздействия на призабойную зону пласта с целью повышения проницаемости. Моделирование данного процесса и прогнозирование технологической и экономической эффективности является важнейшим и необходимым этапом проектирования данной технологии. В основном для моделирования процесса СКО применяют упрощенный подход – задание отрицательных скин-факторов или изменение продуктивности скважин вдоль всего ствола скважин. Однако в условиях протяженных горизонтальных скважин (далее – ГС) данный метод недостаточно применим в связи со сложностью определения скин-фактора по длине ГС, что может привести к завышенным оценкам технологической эффективности. В данной работе опробован подход моделирования, включающий в себя протекание химических реакций. Химическая реакция закладывалась между соляной кислотой и твердым компонентом (известняк), продуктом реакции являлись соль, растворимая в воде, и углекислый газ. Также в качестве исходных данных использовались начальная концентрация твердого компонента (известняк) в составе общей горной породы, константа скорости химической реакции, концентрация и объем закачанной кислоты. В качестве зависимости проницаемости от пористости в гидродинамической модели (далее – ГДМ) использовалось уравнение Кармена-Козени (1).

Целью работы являлись оценка эффективности процесса СКО и многостадийного гидравлического разрыва пласта (далее – МГРП) в условиях протяженных ГС с помощью инструментов композиционного моделирования и оценка объема закачки кислотного раствора.

Моделирование СКО

В используемом гидродинамическом симуляторе Stars возможны 3 способа моделирования СКО в зависимости от целей исследования:

  • изменение коэффициента продуктивности скважины или задание отрицательного скин-фактора: отражение результирующего эффекта от СКО без учета физики процесса (менее точный, но быстрый способ) [1]. Моделирование проводилось путем задания скин-фактора до значения -2, как среднего значения по анализу гидродинамических исследований (далее – ГДИ) с СКО месторождений- аналогов на протяжении 6 мес;
  • использование пользовательских реакций: определение двух компонентов HCl и твердого вещества и добавление одной простой реакции между этими компонентами (компромисс между скоростью и точностью).

В зависимости от способа моделирования требуется свой набор исходных данных:

  • изменение коэффициента продуктивности скважины: промысловые данные о дебитах скважин до и после обработки;
  • использование пользовательских реакций: стандартные исследования СКО на керне (зависимость перепада давления от прокачанных объемов и объем прокачки до образования каверны в зависимости от расхода агента).

Моделирование СКО было выполнено на объекте месторождения N, характеризующегося карбонатным типом коллектора. Краткое описание геолого-физических характеристик представлено в табл. 1.

 

Таблица 1. Характеристика объектов

Параметры

Пласт P1

Пласт C2

Пласт C1

Пласт C1S

Абсолютная отметка кровли, м

-820

-948

-962

-1030

Абсолютная отметка ВНК, м

1632,2

1632,2

Тип коллектора

карбонатный, поровый, кавернозный

карбонатный, поровый, кавернозный

карбонатный, поровый, кавернозный

карбонатный, поровый, кавернозный

Средняя общая толщина, м

91

458

Средняя эффективная нефтенасыщенная толщина, м

41,0

57,1

Проницаемость по ГДИ, 10-³ мкм²

53

53

53

53

Начальная пластовая температура, °С

30,3

31

33,5

35

Начальное пластовое давление, МПа

13,9

14,4

15

15,4

Вязкость нефти в пластовых условиях, мПа∙с

14,00

14,00

14,00

14,00

Плотность нефти в пластовых условиях, кг/м³

870

870

870

870

Плотность нефти в пов. условиях, кг/м³

906

906

906

906

Объемный коэффициент нефти, д. ед.

1,090

1,090

1,090

1,090

Газосодержание, м³/т

40,0

40,0

40,0

40,0

 

Константа скорости реакции задавалась согласно данным месторождения-аналога им. Р. Требса 4,9∙10-5 1/с [2]. Для учета зависимости проницаемости от пористости использовалось уравнение Кармена-Козени:

K(ф) = K0∙[ф/ф0]2∙[(1-ф0)/(1-ф)]2 (1)

В ГДМ задавались следующие химические реакции:

2HCl + CaCO2 = CaCl2 + Н2О + СО2 (2)

4HCl + CaМg(CO3)2 = CaCl2 +

+ МgCl2 + 2Н2О + 2СО2 (3)

где хлористый кальций (CaCl₂) – соль, хорошо растворимая в воде, углекислый газ (CO₂) при соответствующем давлении (свыше 7,6 МПа) также растворяется в воде.

Начальная концентрация CaCO₃ – 0,5 м³/м³.

Моделирование выполнялось на примере проектной скважины 4ПР длиной горизонтального участка 1000 м.

Концентрация кислоты варьировалась от 5 до 15 %. Время закачки – от 1 до 5 сут.

Для оценки длительности обработки призабойной зоны (далее – ОПЗ) были рассмотрены результаты ОПЗ по 4 объектам-аналогам структурной зоны Вала Сорокина Тимано-Печорской нефтегазоносной провинции.

Для моделирования среднее время эффекта от ОПЗ принято равным 6 мес, как среднее значение по аналогам.

 

Таблица 2. Результаты расчетов

Дополнительная добыча нефти за 6 мес, т

Концентрация, %

Время закачки, сут

5

10

15

Средний объем закачки кислотного состава на 1 пог. м длины скважины, м³/пог. м

1

1134

1172,7

1235,8

1,5

3

2465,6

2677,9

3157,2

4,1

5

3590,9

3912,3

4116,6

6,4

скин-фактор -2

3252,4

 

 

Наиболее близкий результат к скин-фактору -2 (по месторождениям-аналогам) показал вариант с закачкой концентрации кислоты 15% и объемом закачки 4,1 м³/пог. м длины ГС.

Моделирование МГРП

Из рассмотренных месторождений достаточный для статистического анализа объем операций гидравлического разрыва пласта (далее – ГРП) имеется только по одному месторождению (9 скв.). По остальным месторождениям-аналогам единичные ГРП по технологическим причинам работы не были выполнены в полном объеме, что не позволяет сделать вывод об их эффективности. Прирост дебита нефти после ГРП составил от 6,1 до 51,7 т/сут, среднее значение – 25,2 т/сут. Рост обводненности по скважинам составил до 16%. Дополнительная добыча на 01.03.2020 г. составила 52,4 тыс. т. Продолжительность эффекта от ГРП составляет более 2 лет, по 8 из 9 скв. эффект продолжается. По проведённым ГРП средняя оценка полудлины трещины составляет 93 м, ширина трещины – 3,2 мм. Средняя оценка безразмерной проводимости составляет 1 ед. Для обеспечения полудлины трещины около 100 м объем закачки пропанта изменяется в диапазоне от 35 до 60 т.

В симуляторе Stars для задания трещин МГРП предлагается 2 способа: планарная система трещин и модель дискретных систем трещин (далее – DFN).

Согласно статистике МГРП на месторождениях-аналогах на секторной ГДМ по ГС 4prod заданы трещины МГРП, полудлина трещины 100 м, ширина 3 мм, высота трещины 70 м, как среднее значение по анализу ГРП по аналогам. Оценка проницаемости проппанта получена с учётом средних оценок безразмерной проводимости FCD = 1: Kf = Kматр∙Xf/w = 35∙100/0,003 = = 1166667 мД. Работа трещин ГРП моделировалась на протяжении трех лет после запуска скважин.

Модель планарных трещин моделируется путем измельчения сетки до ячейки необходимых геометрических размеров (ширина, длина). Ограничением использования планарных трещин является задание направления трещины только по направлению расчетной сетки (по направлению I, J).

Модель DFN (сеть дискретных трещин) используется, чтобы смоделировать трещины как непосредственные ячейки. Идет прямое соединение ячеек модели с трещинами. Задаются направления и формы трещин. Преимуществом модели DFN перед планарными трещинами является возможность учета произвольного направления трещины без привязки к ориентации основной сетки ГДМ.

Как видно из рис. 1, уровни добычи нефти при использовании модели планарных систем трещин и модели DFN находятся на сопоставимом уровне.

 

Рисунок 1. Динамика добычи нефти

 

Рисунок 2. Накопленная добыча нефти

 

Рисунок 3. Время расчета

 

Выводы

Выполнен обзор инструментов моделирования СКО и МГРП в симуляторе Stars на примере проектной скважины. Выполнено моделирование СКО с использованием химических реакций композиционного моделирования и задания отрицательного скин-фактора на проектной скважине 4prod. Получены уровни добычи, сопоставимые с добычей нефти на скважинах после проведения СКО месторождений-аналогов. Преимуществом композиционного моделирования является получение расчетным путем необходимого объема кислотного раствора.

Выполнено моделирование МГРП с использованием двух различных моделей трещин: планарной системы трещин и модели дискретных трещин. Получены сопоставимые уровни добычи. Время расчета также находится на сопоставимом уровне. Преимуществом модели DFN является возможность учета произвольного направления задания трещин по результатам дополнительных специализированных исследований.

Из рассмотренных методов интенсификации притока к скважине технологически наиболее эффективно применение МГРП. Прирост накопленной добычи нефти технологии МГРП по сравнению с СКО составляет 25 тыс. т за 3 года.

Рассмотренные технологии моделирования методов воздействия на призабойную зону ГС рекомендуется использовать при обосновании проектных решений.

×

About the authors

A. A. Terentiyev

“SamaNIPINeft” LLP

Email: terentevaa1@samnipi.rosneft.ru

ведущий инженер отдела разработки месторождений высоковязкой нефти

Russian Federation, Samara city

S. M. Durkin

“SamaNIPINeft” LLP

Email: durkinsm@samnipi.rosneft.ru

главный эксперт по разработке

Russian Federation, Samara city

Konstantin V. Pchela

“SamaNIPINeft” LLP

Author for correspondence.
Email: PchelaKV@samnipi.rosneft.ru

главный менеджер проекта

Russian Federation, Samara city

References

  1. Computer Modeling Group Ltd. STARS, CMG: Version 2020 User’s Guide. – Calgary, Alberta, Canada, 2020.
  2. Фоломеев А.Е. Совершенствование технологии кислотного воздействия на высокотемпературные карбонатные коллекторы. Дисс. канд. техн. наук. – Уфимский государственный нефтяной технический университет, Уфа, 2020.

Supplementary files

Supplementary Files
Action
1. Figure 1. Dynamics of oil production

Download (90KB)
2. Figure 2. Cumulative oil production

Download (47KB)
3. Figure 3. Calculation time

Download (37KB)

Copyright (c) 2022 Terentiyev A.A., Durkin S.M., Pchela K.V.

Creative Commons License
This work is licensed under a Creative Commons Attribution-NonCommercial-NoDerivatives 4.0 International License.

This website uses cookies

You consent to our cookies if you continue to use our website.

About Cookies