Научная статья на тему 'Численная реализация механизма термогазового воздействия на двумерной модели'

Численная реализация механизма термогазового воздействия на двумерной модели Текст научной статьи по специальности «Энергетика и рациональное природопользование»

CC BY
375
84
i Надоели баннеры? Вы всегда можете отключить рекламу.
Ключевые слова
БАЖЕНОВСКАЯ СВИТА / ТЕРМОГАЗОВОЕ ВОЗДЕЙСТВИЕ / ДРЕНИРУЕМЫЕ И НЕДРЕНИРУЕМЫЕ ПЛАСТЫ / КЕРОГЕН / КОКС / ВНУТРИПЛАСТОВОЕ ГОРЕНИЕ / МОДЕЛИРОВАНИЕ / 2D МОДЕЛЬ / МЕТОДЫ УВЕЛИЧЕНИЯ НЕФТЕОТДАЧИ / DEPOSITS OF THE BAZHENOV SUITE / THERMAL GAS TECHNOLOGY / DRAINABLE AND NONDRAINABLE FORMATIONS / KEROGEN / COKE / INTERBEDDING COMBUSTION / MODELING / 2D MODEL / ENHANCED OIL RECOVERY METHODS

Аннотация научной статьи по энергетике и рациональному природопользованию, автор научной работы — Шахмаев А. М., Симаков Я. О., Пятибратов П. В., Мосесян А. А.

В настоящее время на большинстве месторождений России добыча нефти из пластов с традиционными запасами снижается, и, несмотря на имеющиеся тенденции развития альтернативных источников энергии, потребление углеводородных ресурсов, согласно различным опубликованным материалам, только будет расти [1]. Одним из источников нетрадиционных запасов углеводородов являются пласты баженовской свиты. По оценке разных экспертов ее ресурсы оцениваются примерно в 100 млрд тонн. Эта оценка не включает углеводородный ресурс керогена (несформировавшейся нефти), содержащийся в баженовской свите. Данная статья представляет описание на двумерной модели технологии моделирования термогазового воздействия (далее ТГВ) на пласты баженовской свиты, с учетом особенностей последней и допущений, связанных с ограничениями симулятора. Материалы и методы На сегодняшний момент рассматриваются следующие методы: разработка баженовской свиты на естественном режиме; гидроразрыв пласта; сода-ПАВ-полимерное воздействие на пласт (ASP-заводнение); ТГВ. ТГВ является наиболее перспективным, с технологической точки зрения, способом разработки. В то же время ТГВ сложный физико-химический процесс, который требует реализации на гидродинамической модели для оценки технологических показателей. Для проведения данных работ выбран гидродинамический симулятор CMG STARS, который позволяет моделировать тепловые и газовые методы увеличения нефтеотдачи и считается одним из лучших симуляторов в этой области. Итоги На двумерной модели поэтапно проведена численная реализация ключевых процессов и явлений, протекающих при ТГВ применительно к пластам баженовской свиты. Перед численной реализацией ТГВ для условий баженовской свиты, были решены вопросы, касающиеся значений ФЕС и их изменения в процессе воздействия. На первом этапе моделирования выделены основные зоны, характеризующие процесс, и дано физическое обоснование отдельным показателям. На втором этапе в линейную модель была дополнительно введена первоначально непроницаемая толща, характеризующая матрицу. Проведены расчеты вариантов с вовлечением матрицы в процесс дренирования. Выводы При реализации ТГВ на двумерной модели в условиях баженовской свиты, с учетом лабораторных исследований керна и флюидов пластов баженовской свиты, были получены качественные результаты. Отличием от предыдущих работ, выполненных на эту тему, является то, что авторами был описан подход к моделированию ТГВ и реализован механизм динамического изменения ФЕС от изменения термобарических условий в недренируемой толще, то есть процесс ТГВ численно реализован полностью в рамках одного гидродинамического симулятора CMG STARS. Стоит отметить, что полученные результаты приводятся в единицах относительного прироста, малые абсолютные числа прироста связаны с размером моделируемого объекта линейная модель. В дальнейшем авторами планируется реализовать процесс ТГВ на полномасштабной модели, с воспроизведением реальных промысловых данных. Данная работа может послужить хорошей базой для оценки характерных зон ТГВ и получения сопоставимых результатов.

i Надоели баннеры? Вы всегда можете отключить рекламу.

Похожие темы научных работ по энергетике и рациональному природопользованию , автор научной работы — Шахмаев А. М., Симаков Я. О., Пятибратов П. В., Мосесян А. А.

iНе можете найти то, что вам нужно? Попробуйте сервис подбора литературы.
i Надоели баннеры? Вы всегда можете отключить рекламу.

Numerical implementation of the thermal gas technology mechanism in the 2D model

A decline in the oil production from the Russian nonshale reservoirs has currently been happening; regardless of the existing alternate energy source development trends,hydrocarbon resource consumption will expand, according to several released reports [1]. Deposits of the Bazhenov suite are one of the alternate hydrocarbon reserves.This article describes a 2D modeling of thermal gas technology (hereinafter TGT) for the deposits of the Bazhenov suite considering peculiarities of the latter, and allowances related to the simulator limitations. Materials and methods The methods of the deposits of the Bazhenov suite depletion, hydraulic fracturing, sodium surfactant polymer treatment of formation (ASPflooding), TGT are considered at the moment. Results The 2D model shows а step-by-step numerical implementation of the key processes and occurrences ongoing during the thermal gas technology with regard to the deposits of the Bazhenov suite.The first phase of the modeling allots main areas characterizing the process, and gives a physical rational to the certain indicators.The second phase introduces a primarily impermeable sequence into the model characterizing the matrix. The options are calculated with the matrix involved in the drainage process. Conclusions Qualitative results are obtained during the TGT at the deposits of the Bazhenov suite conditions considering the 2D model, core and formation fluids laboratory analysis. This work differs from those previously released in the fact that the authors describe a certain approach to the TGT modeling, and implement a mechanism of the dynamic change of the reservoir porosity and permeability from changing the pressure and temperature conditions in the nondrainable sequence; i.e. the TGT process is numerically implemented in full by one hydrodynamic simulator, which is CMG STARS.

Текст научной работы на тему «Численная реализация механизма термогазового воздействия на двумерной модели»

добыча

УДК 622.276

Численная реализация механизма термогазового воздействия на двумерной модели

А.м. шахмаев

аспирант1

[email protected] Я.О. Симаков

заместитель начальника по моделированию департамента разработки месторождений2 [email protected]

П.В. Пятибратов

к.т.н., доцент кафедры разработки и эксплуатации нефтяных месторождений1 [email protected]

А.А. мосесян

инженер кафедры бурения нефтяных и газовых скважин1

[email protected]

1РГУ нефти и газа (НИУ) имени И. М. Губкина, г. Москва, Россия 2ООО «Петек», г. Москва, Россия

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

Данная статья представляет описание на двумерной модели технологии моделирования термогазового воздействия (далее — ТГВ) на пласты баженовской свиты, с учетом особенностей последней и допущений, связанных с ограничениями симулятора.

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

Стоит отметить, что термогазовое воздействие является перспективным методом вовлечения в разработку углеводородного потенциала керогена, который позволяет использовать преимущества известных методов увеличения нефтеотдачи (далее - МУН). Прогнозный КИН для ТГВ достигает 40% [2, 3].

Пласты баженовской свиты обладают высоким энергетическим потенциалом. Начальные пластовые условия являются благоприятными для инициации внутрипластовых окислительных процессов [4]. Повсеместное наличие керогена в скелете породы позволяет, при дальнейшем развитии процесса, получить высокотемпературный пиролиз с образованием тяжелого остатка (кокса), который создаст благоприятные условия для развития высокотемпературного фронта горения [5, 6].

Рассмотрим основные особенности ба-женовской свиты, накладывающие отпечаток на моделирование ТГВ:

• большая доля глинистых и кероген-глинистых, практически непроницаемых толщ, формирующих недренируемую зону («матрицу») — отношение по толщине к дренируемой зоне в среднем 4:1;

• наличие нефтегенерационного материала (керогена), распределенного по всему объему дренируемой и недренируемой зон;

• наличие пустотного пространства матрицы, занятого жидкими углеводородами (далее — УВ), запертыми в непроницаемых породах (общие запасы жидких УВ могут в 2-3 раза превышать запасы дренируемой зоны).

Рис. 1 показывает основные зоны, формирующиеся при реализации ТГВ.

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

• процесс термодеструкции керогена при температурах свыше 300-350 оС, приводящий к выходу жидкой фазы (так называемая «синтетическая» нефть), попутного газа и тяжелого остатка (кокса);

• горение кокса — как основной источник энергии, генерирующийся в пласте;

• процесс жидкофазного окисления тяжелых компонентов нефти и керогена;

• изменение фильтрационно-емкостных свойств (далее — ФЕС) коллектора при разложении керогена и последующем сгорании твердого остатка (кокса), с последующим увеличением пустотного пространства, занятого флюидами;

• изменение ФЕС коллектора в процессе теплового воздействия и изменение давления (растрескивание породы при нагреве, проявление упругих сил, геомеханические эффекты, характерные для упругого и упруго-пластичного сжатия/расширения). К этим процессам, в том числе, относится вовлечение матрицы в разработку за счет охвата температурным воздействием ее приграничных областей и последующего частичного раскрытия микротрещин, образовавшихся в процессе нефтегенерации и пребывающих в настоящее время в «запечатанном» состоянии.

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

Рис. 1 — Общий вид зон при реализации ТГВ Fig. 1 — General view of zones of TGT

материалы и методы

На сегодняшний момент рассматриваются следующие методы: разработка баженовской свиты на естественном режиме; гидроразрыв пласта; сода-ПАВ-полимерное воздействие на пласт (ASP-заводнение); ТГВ. ТГВ является наиболее перспективным, с технологической точки зрения, способом разработки. В то же время ТГВ — сложный физико-химический процесс, который требует реализации на гидродинамической модели для оценки технологических показателей. Для проведения данных работ выбран гидродинамический симулятор CMG STARS, который позволяет моделировать тепловые и газовые методы увеличения нефтеотдачи и считается одним из лучших симуляторов в этой области.

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

баженовская свита, термогазовое воздействие, дренируемые и недренируемые пласты, кероген, кокс, внутрипластовое горение, моделирование, 2D модель, методы увеличения нефтеотдачи

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

Поэтапная реализация математического описания процесса ТГВ

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

• формализации подходов к определению параметров, характеризующих пустотное пространство пород баженовской свиты и насыщающих его веществ (пустотность, трещиноватость, свободная нефть, Total Organic Carbon — степень насыщения органикой нефтематеринских пород и т. д.);

• изменения ФЕС первоначально непроницаемой миктротрещиноватой матрицы при воздействии на нее процессов термической обработки с сопутствующим изменением пластового давления.

Далее процесс ТГВ на породы баженовской свиты можно условно разделить на несколько составляющих:

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

2. Вовлечение в разработку нефтекероген-содержащих зон пласта за счет тепловой обработки с последующим пиролизом и высвобождением более легких жидких и газообразных УВ.

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

В работе [7] приведено описание основных особенностей применения термогазового воздействия на баженовскую свиту и отмечено, что компьютерное моделирование подтвердило перспективность

применения ТГВ и целесообразность его развития.

В работах [8-14] достаточно подробно описаны математические подходы к численной реализации вышеперечисленных пунктов. Однако, особенную трудность у авторов вызывала задача вовлечения в процесс дренирования непроницаемых толщ матрицы. В работах [7-12] авторы не описывали достаточно подробно реализацию данного эффекта, сконцентрировав свои работы на исследовании других, не менее важных процессов ТГВ. В работе [13] был предложен коэффициент для учета добычи нефти из не-дренируемой зоны.

В работе [14] процесс ТГВ реализован с помощью гидродинамического симулятора CMG STARS и внешней среды для программирования, которая и учитывает эффект вовлечения матрицы в процесс разработки. Отличительная особенность данной работы заключается в том, что механизм ТГВ численно реализован полностью в рамках одного гидродинамического симулятора. В статье описаны возможности CMG STARS по учету изменения ФЕС от изменения термобарических условий в недренируемой зоне, и сделан наиболее оптимальный выбор с точки зрения времени счета и воспроизведения эффекта. Стоит сразу же сказать, что успешная реализация данной задачи должна обеспечить важнейшую составляющую эффекта от ТГВ в условиях баженовской свиты.

Эффект вовлечения недренируемой зоны связан с тем, что запасы свободной нефти литотипов матрицы (микротрещиноватой породы) существенно превышают аналогичные запасы дренируемой зоны (по отдельным оценкам, более чем в 2 раза). Запасы керогена, связанные с глинистыми породами матрицы, также кратно больше и достигают 40% от объема породы (при 4-10 % в карбонатах дренируемой зоны) [14].

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

Т.

хг X

Рис. 2 — Кривая изменения температуры при сухом внутрипластовом горении для

прямолинейного пласта [15] Fig. 2 — Curve of temperature change during dry in-situ combustion for straight-line formation [15]

500

0 20 40 60 SO 100 120 140

Расстояние, m

Рис. 3 — Кривая изменения температуры при внутрипластовом горении, полученная на двумерной модели

Fig. 3 — 2D interbedding in-situ combustion temperature curve

Размерность

Количество Размеры

Начальная пористость

500x1x5

ячеек ячеики матрица

2500 1x1x1 0,09

Дренируемая Нефтенасыщенность зона

0,062 0,85

Таб. 1 — Основные параметры размерности линейной модели Tab. 1 — Basic 2D dimension parameters

дренируемая зона

КЕРОГЕН

0,097 (доли ед.)

116,4 (кг/м3) 1830 (г/моль/м3)

жидкие УВ Сумма (потенциал)

0,062

0,159

матрица

КЕРОГЕН жидкие Сумма УВ (потенциал)

0,115 (доли ед.)

200 (кг/м3) 0,09

1685 (г/моль/м3)

0,205

Таб. 2 — Средние параметры пустотности, использовавшиеся в линейной модели Tab. 2 — Average 2D model voids parameters

Показатель

Коэффициент теплопроводности, Дж/-м-сут-°С Объемная теплоемкость, Дж/м3-°С

Горные породы Коллектор Окружающие Кероген Нефть породы

1,8-105 1,8-105

Пластовые жидкости Вода

1,0105 1,15104

5,35-104

2,3-106

2,3-106

Таб. 3 — Теплофизические свойства пластовых жидкостей и породы коллектора Tab. 3 — Thermal properties of the fluids and reservoir rock

Рис. 4 — Зависимость проницаемости от температуры, в недренируемой зоне Fig. 4 — Permeability versus temperature in the non-drainable zone

Формула компонента

Реакция

Энтальпия Энергия Константа реакции, активации, Аррениуса, Дж/моль Дж/моль 1/(сутки*кПа)

Низкотемпературное окисление -••> HO+ +CO 1,2107

1,2107

100 000 100 000

2-108 2-108

Термолиз (постоянная скорость реакции вне зависимости от температуры)

iНе можете найти то, что вам нужно? Попробуйте сервис подбора литературы.

Kerogen Kerogen.....•> C2-C4+ +Synt_Oil + Coke

36 000

Высокотемпературное окисление

Coke Coke + -O2.....•>......^H20+C02 9,0^ 107 21 000

Kerogen Kerogen + O2.....•>.....•> H20 +^CO2 9,0407 21 700

100

10,0 10,0

Таб. 4 — Базовые уравнения и основные параметры химических реакций Tab. 4 — Basic equations and basic parameters of chemical reactions

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

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

Температура Т* — положение фронта горения с координатой Хф. На этой кривой виден излом на расстоянии Хт от начала координат. Это связано с конвективным переносом тепла. Сечение пласта с приближенной координатой Х = Хт называется фронтом конвекции. В процессе внутрипластового горения фронты горения и конвекции с координатами Хф и Хт перемещаются по ходу движения закачиваемых в пласт веществ [15].

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

моделирование матрицы

Основной целью вовлечения матрицы является увеличение КИН за счет дополнительного притока легкой нефти, «закупоренной» в матрице. Для начала необходимо решить задачу о динамическом изменении ФЕС, которое обеспечит вовлечение матричной толщи в процесс фильтрации. Кроме того, не стоит забывать о постоянном изменении ФЕС дренируемой зоны в результате колебания давления и температуры в процессе внутрипластового горения. Учет всех этих эффектов имеет важное значение для более реальной оценки эффекта внутрипластового горения на пласты баженовской свиты.

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

Основные параметры размерности модели представлены в таб. 1, параметры пустотности для дренируемой и недренируемой зон представлены в таб. 2, теплофизические параметры пласта и флюидов представлены в таб. 3.

Возможности STARS по динамическому изменению ФЕС. Выбор механизма

Существует несколько вариантов реализации в CMG STARS механизмов изменения ФЕС коллекторов:

• встроенная модель упругого расширения/ сжатия, (Dilation/Recompaction);

• механизм разрыва глинистой перемычки;

С18+

С5С7

с5с7+ o2.....h2o+co2

Рис. 5 — Начальное распределение проницаемости в двумерной модели с матрицей Fig. 5 — Initial permeability distribution in the 2D model with a matrix

Рис. 6а — Изменение основных параметров процесса ТГВ в дренируемой зоне и в матрице — температура (50 лет закачки) Fig. 6a — Change of the basic TGTparameters in the drainable zone and in the matrix -temperature (50 years of injection)

Рис. 6б — Изменение основных параметров процесса ТГВ в дренируемой зоне и в матрице — проницаемость (50 лет закачки) Fig. 6б — Change of the basic TGT parameters in the drainable zone and in the matrix — permeability (50 years of injection)

• модель растрескивания породы Бартона-Бэндиса [16].

Перед тем как выбирать механизм динамического изменения ФЕС, необходимо четко сформулировать, какие эффекты, связанные с процессом ТГВ, приводят к подобным изменениям.

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

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

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

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

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

- сгорания топлива (твердого остатка при реакции пиролиза).

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

Предварительно, необходимо отметить общие ограничения симулятора.

Ячейка, проинициализированная нулевой проницаемостью, не будет участвовать в процессе массопереноса, но будет проводить и аккумулировать тепло, т.е. изменение её первоначально нулевой проницаемости полностью исключено. Это касается всех методов, кроме механизма разрыва глинистой перемычки. Изменение проницаемости ячеек необратимо (кроме модели Бартона-Бэндиса)

Отсутствует возможность задания табличной функции зависимости проницаемости от давления и температуры. Только один из двух механизмов может применяться для каждого литотипа.

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

механизм разрыва глинистой перемычки

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

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

Зависимость проницаемости от температуры, используемая в модели, приведена на рис. 4.

Основные химические реакции, протекающие в пласте при ТГВ представлены в таб. 4.

Реализация ТГВ на двумерной модели с дренируемой зоной и непроницаемой матрицей

Рис. 5 показывает начальное распределение проницаемости в линейной модели ТГВ с матрицей.

На рис. 6 (а, б, в, г) отражены изменения основных параметров, за счет которых

Расстояние,!*

I I I I I \ I I I I I I I I I I I I I * I I I I t I I I I I I I I t I I I I * p I I I I I I I t I I I

0 ICQ 200 300 400 ЯЮ

Рис. 6в — Изменение основных параметров процесса ТГВ в дренируемой зоне и в матрице — пористость (50 лет закачки) Fig. 6в — Change of the basic TGT parameters in the drainable zone and in the matrix — porosity (50 years of injection)

Расстояние. M

I p I I I I I I i I I I I J I i p I I I I I I I I I I I I I j I I I I I p I I I I I I p I I I I I I I

о 100 2oa 300 taa 500

Рис. 6г — Изменение основных параметров процесса ТГВ в дренируемой зоне и в матрице — содержание твердого остатка (50 лет закачки) Fig. 6г — Change of the basic TGT parameters in the drainable zone and in the matrix — solid residue content (50 years of injection)

достигается главный эффект от ТГВ — дополнительная добыча жидких УВ. Проиллюстрировано изменение температуры, проницаемости и пористости, а также содержания твердого остатка (кокса) на конечный момент времени, с вовлечением матрицы в процесс ТГВ. На рис. 7 приведен график накопленной добычи нефти, иллюстрирующий дополнительный эффект от вовлечения матрицы.

Фронт горения определяется как зона с наибольшей температурой, а также максимальным содержанием кокса, который образовывается перед зоной наибольшей температуры. [17]. Рис. 6а и 6г показывают, что за 50 лет фронт горения продвинулся на 25 метров, но стоит отметить, что его положение сильно зависит от множества факторов, таких как:

• интенсивность окисления - скорость реакции;

• объем закачиваемого окислителя, его концентрация;

• скорость его фильтрации;

• характеристика системы нефть - коллектор, а именно содержание керогена (топлива) для реакции горения.

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

На рис. 8 представлено распределение начальной нефтенасыщенности, на рис. 9 — нефтенасыщенность через 1,5 года добычи. Показан резкий рост накопленной добычи на рис. 7 за первые полтора года — отбор начальных запасов нефти (начальной нефтенасыщенности) дренируемой зоны. Последующий рост и стабилизация накопленной добычи нефти (рис. 6) связаны с приходом очередных валов нефти, за счет реализации термогазового воздействия ТГВ. На рис. 9 это показано полосой более насыщенного зеленого цвета.

Расчетный технологический эффект за 50-летний период реализации ТГВ на модели (с вовлечением матрицы в разработку) составил около 50% дополнительной добычи нефти, относительно варианта ТГВ без вовлечения матрицы.

Итоги

На двумерной модели поэтапно проведена численная реализация ключевых процессов и явлений, протекающих при ТГВ применительно к пластам баженовской свиты. Перед численной реализацией ТГВ для условий ба-женовской свиты, были решены вопросы, касающиеся значений ФЕС и их изменения в процессе воздействия. На первом этапе моделирования выделены основные зоны, ха-растеризующие процесс, и дано физическое обоснование отдельным показателям. На втором этапе в линейную модель была дополнительно введена первоначально непроницаемая толща, характеризующая матрицу. Проведены расчеты вариантов с вовлечением матрицы в процесс дренирования.

Рис. 7 — Сопоставление накопленной добычи нефти с вовлечением и без вовлечения

матрицы, за 50 лет разработки Fig. 7 — Correlation between cumulative oil production with and without the matrix during 50

years of the development

Выводы

При реализации ТГВ на двумерной модели в условиях баженовской свиты, с учетом лабораторных исследований керна и флюидов

пластов баженовской свиты, были получены качественные результаты. Отличием от предыдущих работ, выполненных на эту тему, является то, что авторами был описан подход к моделированию ТГВ и реализован механизм динамического изменения ФЕС от изменения термобарических условий в недренируемой толще, то есть процесс ТГВ численно реализован полностью в рамках одного гидродинамического симулятора - CMG STARS. Стоит отметить, что полученные результаты приводятся в единицах относительного прироста, малые абсолютные числа прироста связаны с размером моделируемого объекта - линейная модель. В дальнейшем авторами планируется реализовать процесс ТГВ на полномасштабной модели, с воспроизведением реальных промысловых данных. Данная работа может послужить хорошей базой для оценки характерных зон ТГВ и получения сопоставимых результатов.

Авторы выражают огромную благодарность А.А. Боксерману за ценные советы и консультации по методу термогазового воздействия.

Список литературы

1. World Energy Outlook 2014 // International Energy Agency. 2014. [Электронный ресурс]. Режим доступа: https://www.iea.org/publications/ freepublications/publication/WEO2014. pdf

2. Алекперов В.Ю., Грайфер В.И., Николаев Н.М., Карпов В.Б., Кокорев

B.И. и др. Новый отечественный способ разработки месторождений баженовской свиты (часть 1) // Нефтяное хозяйство. 2013. №12.

C. 100-105.

3. Алекперов В.Ю., Грайфер В.И., Николаев Н.М., Карпов В.Б., Кокорев

B.И. и др. Новый отечественный способ разработки месторождений баженовской свиты (часть 2) // Нефтяное хозяйство. 2014. №1.

C. 50-53.

4. Сонич В.П., Батурин Ю.Е., Малышев

A.Г., Зарипов О.Г., Шеметилло В.Г. Проблемы и перспективы освоения Баженовской свиты // Нефтяное хозяйство. 2001. №9. С. 63-68.

5. Палий А.П. Результаты и перспективы развития термогазового способа разработки месторождений с нетрадиционными запасами. XVI конференция молодых специалистов, работающих в организациях, осуществляющих деятельность, связанную с использованием участков недр на территории Ханты-мансийского автономного округа - Югры, тезисы докладов. Ханты-Мансийск, 2016.

6. Батурин Ю.Е., Сонич В.П., Малышев А.Г., Зарипов О.Г., Шеметилло В.Г. Оценка перспектив применения гидротермовоздействия в пласте ЮО месторождений ОАО «Сургутнефтегаз» // Интервал. 2002. №1. С. 17-36.

7. А.А. Боксерман, Кокорев В.И., Плынин

B.В., Ушакова А.С. Современное состояние и перспективы применения термогазового метода увеличения

Расстояний, и

| I I I I Т Г Г-Т-Т-у Т-Г-1 I J I I Т-Т J-r-L I I I I I I J I I I I I I I ■ I I I I I ■ I -Г Т-1-Г-Г" у

С IDC ЗАО 2CL Ш bJJ

Рис. 8 — Распределение начальной нефтенасыщенности Fig. 8 — Initial oil saturation distribution

Расстояние. hi

'I "Г Г П Т 'I I T I "I Г Г""Г"| ' f' l 1 "T I Г"Г V I I rfTTfT'rrn'f t Г I' f t 1 I' f-l-r'T^T—r-p о icc m ЭСС ш w

Рис. 9 — Распределение нефтенасыщенности через 1,5 года добычи Fig. 9 — Oil saturation distribution in 1.5 years of the development

нефтеотдачи на месторождениях баженовской свиты. II Международная конференция. Наноявления при разработке месторождений углеводородного сырья: от наноминералогии и нанохимии к нанотехнологиям, тезисы докладов. Москва, 2010.

8. Ямбаев М. Ф. Основные особенности термогазового увеличения нефтеотдачи применительно к условиям сложнопостроенных коллекторов. М., 2006. 153 с.

9. Афанаскин И. В. Повышение технологической эффективности метода направленной закачки воздуха в нефтяные пласты на основе численного моделирования и результатов гидродинамических исследований скважин. М., 2013. 273 с.

10. Диева Н. Н. Гидродинамическое моделирование термохимического воздействия на пласты трудноизвлекаемых углеводородов. М., 2015. 113 с.

11. Щеколдин К. А. Обоснование технологических режимов термогазового воздействия на залежи

баженовской свиты. М., 2016. 126 с.

iНе можете найти то, что вам нужно? Попробуйте сервис подбора литературы.

12. Шевченко А.В. Разработка и реализация численных методов моделирования многокомпонентной неизотермической фильтрации. М., 2015. 96 с.

13. Кокорев В.И. Технико-технологические основы инновационных методов разработки месторождений

с трудноизвлекаемыми и нетрадиционными запасами нефти. М., 2010. 318 с.

14. Алексеев Ю.В., Ерофеев А.А., Пачежерцев А.А., Меретин А.С., Никитин Р.Н. Перспективы применения термохимических методов воздействия при разработке залежей баженовской свиты // Нефтяное хозяйство. 2015. №9. С. 93-97.

15. Желтов, Ю.П. Разработка нефтяных месторождений. М.: Недра, 1986. 332 с.

16. Бурже Ж., Сурио П., Комбарну М. Термические методы повышения нефтеотдачи пластов. М.: Недра, 1989. 422 с.

17. Computer Modeling Group LTD. Stars user guide. Advanced processes and thermal reservoir simulator, 2016.

ENGLISH

OIL PRODUCTION

Numerical implementation of the thermal gas technology mechanism in the 2D model

Authors:

Artur M. Shakhmaev — post-graduate student1; [email protected]

Yaroslav O. Simakov — deputy director for production engineering2; [email protected]

Petr V. Pyatibratov — associate professor1; [email protected]

Ashot A. Mosesyan — engineer1; [email protected]

1Gubkin Russian State University of Oil and Gas (National Research University), Moscow, Russian Federation 2LLC "Petek", Moscow, Russian Federation

UDC 622.276

Abstract

A decline in the oil production from the Russian nonshale reservoirs has currently been happening; regardless of the existing alternate energy source development trends, hydrocarbon resource consumption will expand, according to several released reports [1]. Deposits of the Bazhenov suite are one of the alternate hydrocarbon reserves. This article describes a 2D modeling of thermal gas technology (hereinafter - TGT) for the deposits of the Bazhenov suite considering peculiarities of the latter, and allowances related to the simulator limitations.

Materials and methods

The methods of the deposits of the Bazhenov suite depletion, hydraulic fracturing, sodium surfactant polymer treatment of formation (ASP

flooding), TGT are considered at the moment. Results

The 2D model shows a step-by-step numerical implementation of the key processes and occurrences ongoing during the thermal gas technology with regard to the deposits of the Bazhenov suite.

The first phase of the modeling allots main areas characterizing the process, and gives a physical rational to the certain indicators. The second phase introduces a primarily impermeable sequence into the model characterizing the matrix. The options are calculated with the matrix involved in the drainage process.

Conclusions

Qualitative results are obtained during the TGT

at the deposits of the Bazhenov suite conditions considering the 2D model, core and formation fluids laboratory analysis. This work differs from those previously released in the fact that the authors describe a certain approach to the TGT modeling, and implement a mechanism of the dynamic change of the reservoir porosity and permeability from changing the pressure and temperature conditions in the nondrainable sequence; i.e. the TGT process is numerically implemented in full by one hydrodynamic simulator, which is CMG STARS.

Keywords

deposits of the Bazhenov suite, thermal gas technology, drainable and nondrainable formations, kerogen, coke, interbedding combustion, modeling, 2D model, enhanced oil recovery methods

References

1. World Energy Outlook 2014, International Energy Agency, 2014. Available at: https:// www.iea.org/publications/freepublications/ publication/WEO2014.pdf

2. Alekperov V.Yu., Grayfer V.I., Nikolaev N.M., Karpov V.B., Kokorev V.I. and oth. Novyy otechestvennyy sposob razrabotki mestorozhdeniy bazhenovskoy svity (chast'

1) [A new domestic method for developing deposits of the Bazhenov suite (part 1). Oil industry, 2013, Issue 12, pp. 100-105.

3. Alekperov V.Yu., Grayfer V.I., Nikolaev N.M., Karpov V.B., Kokorev V.I. and oth. Novyy otechestvennyy sposob razrabotki mestorozhdeniy bazhenovskoy svity (chast'

2) [A new domestic method for developing deposits of the Bazhenov suite (part 2). Oil industry, 2014, Issue 1, pp. 50-53.

4. Sonich V.P., Baturin Yu.E., Malyshev A.G., Zaripov O.G., Shemetillo V.G. Problemy

i perspektivy osvoeniya Bazhenovskoy svity [Problems and prospects for the development of the Bazhenov suite]. Oil industry, 2001, issue 9, pp. 63-68.

5. Paliy A.P. Rezul'taty iperspektivyrazvitiya termogazovogo sposoba razrabotki mestorozhdeniy s netraditsionnymi zapasami [Results and prospects of development of the thermogas method of development of deposits with non-traditional reserves]. XVI conference of young specialists, who working in company which develop oil field in Khanty-Mansiisk autonomous district —Yugra. Khanty-Mansiysk, 2016.

6. Baturin Yu.E., Sonich V.P., Malyshev A.G., Zaripov O.G., Shemetillo V.G.

Otsenka perspektivprimeneniya gidrotermovozdeystviya vplaste YuO mestorozhdeniy OAO «Surgutneftegaz» [Estimation of the prospects for the use

of hydroelectricity in the reservoir of the South fields of OJSC "Surgutneftegas"]. Interval, 2002, issue 1, pp. 17-36.

7. A.A. Bokserman, Kokorev V.I., Plynin V.V., Ushakova A.S. Sovremennoe sostoyanie i perspektivy primeneniya termogazovogo metoda uvelicheniya nefteotdachi na mestorozhdeniyakh bazhenovskoy svity [Current state and prospects

for the application of the thermogas method of increasing oil recovery at the deposits of the Bazhenov suite]. II International conference. Nanophenomena during development of oil field: from nanominerology and nanochemistry to nanotechnology. Moscow, 2010.

8. Yambaev M. F. Osnovnye osobennosti termogazovogo uvelicheniya nefteotdachiprimenitel'no kusloviyam slozhnopostroennykh kollektorov [The main features of the thermogas increase in oil recovery in the context of complex reservoirs]. Moscow, 2006,153 p.

9. Afanaskin I. V. Povyshenie tekhnologicheskoy effektivnosti metoda napravlennoy zakachki vozdukha v neftyanye plasty na osnove chislennogo modelirovaniya i rezul'tatov gidrodinamicheskikh issledovaniy skvazhin [Increase of technological efficiency of the method of directed injection of air into

oil reservoirs on the basis of numerical simulation and the results of hydrodynamic studies of wells]. Moscow, 2013, 273 p.

10. Dieva N. N. Gidrodinamicheskoe modelirovanie termokhimicheskogo vozdeystviya na plasty trudnoizvlekaemykh uglevodorodov [Hydrodynamic modeling of thermochemical impact on formations of hard-to-recover hydrocarbons]. Moscow, 2015, 113 p.

11. Shchekoldin K. A. Obosnovanie

tekhnologicheskikh rezhimov termogazovogo vozdeystviya na zalezhi bazhenovskoy svity [Justification of technological modes of thermogas effect on deposits of the Bazhenov suite]. Moscow, 2016, 126 p.

12. Shevchenko A.V. Razrabotka i realizatsiya chislennykh metodov modelirovaniya mnogokomponentnoy neizotermicheskoy fil'tratsii [Development and implementation of numerical methods for modeling multicomponent nonisothermal filtration]. Moscow, 2015, 96 p.

13. Kokorev V.I. Tekhniko-tekhnologicheskie osnovy innovatsionnykh metodov razrabotki mestorozhdeniy s trudnoizvlekaemymi i netraditsionnymi zapasami nefti [Technical and technological fundamentals of innovative methods for developing deposits with hard-to-recover and unconventional oil reserves]. Moscow, 2010, 318 p.

14. Alekseev Yu.V., Erofeev A.A., Pachezhertsev A.A., Meretin A.S., Nikitin R.N. Perspektivy primeneniya termokhimicheskikh metodov vozdeystviya pri razrabotke zalezhey bazhenovskoy svity [Prospects of application of thermochemical methods of influence in the development of deposits of the Bazhenov suite]. Oil industry, 2015, issue 9, pp. 94-97.

15. Zheltov, Yu.P. Razrabotka neftyanykh mestorozhdeniy [Development of oil deposits]. Moscow: Nedra, 1986, 332 p.

16. Burzhe Zh., Surio P., Kombarnu M. Termicheskie metody povysheniya nefteotdachi plastov [Thermal methods for enhanced oil recovery]. M.: Nedra, 1989, 422 p.

17. Computer Modeling Group LTD. Stars user guide. Advanced processes and thermal reservoir simulator, 2016.

i Надоели баннеры? Вы всегда можете отключить рекламу.