К оглавлению журнала

 

УДК 622.276.43

© С.Д. Богданов , 1997

ВЛИЯНИЕ КОЭФФИЦИЕНТА ДИФФУЗИИ НА ВРЕМЯ И МЕСТО УСТАНОВКИ ГЕЛЕ ВО ГО ЭКРАНА ПРИ ПРИМЕНЕНИИ МЕТОДА СИЛИКАТНО-ЩЕЛОЧНОГО ЗАВОДНЕНИЯ

С.Д. Богданов (ВНИГНИ)

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

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

В отличие от полимерного раствора реагент доступен и значительно более дешев. При прокачке обычной воды гель может быть разрушен и высокопроницаемые слои могут быть снова подключены к разработке. Метод прошел апробацию в промысловых условиях. Получен положительный результат при закачке раствора в Казахстане (месторождение Мартыши) и Западной Сибири (Талинское месторождение и др.). Метод недостаточно изучен в части возможности его применения в карбонатных коллекторах. Требуется уточнение и влияния высоких температур (свыше 80 °С) на устойчивость процесса в связи с возможной деструкцией ряда присадок (добавок), используемых для технологических целей, так как они содержат полимерные компоненты.

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

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

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

где S1, S2 - насыщенность соответственно водной и нефтяной фаз;

С - концентрация активной примеси в растворе;

а - равновесная сорбция;

pi - давление в соответствующей фазе;

рk - капиллярное давление;

ri - плотность фаз;

mi - вязкость фаз;

Wi - скорость фаз;

fi - фазовая проницаемость;

m - пористость;

k - тензор абсолютной проницаемости среды;

g - ускорение свободного падения;

Qi - функции стоков-источников (скважины);

d - дельта-функция Дирака;

D - коэффициент диффузии активного компонента;

х - xi - координаты.

Введенный в закон Дарси для водной фазы параметр R как функция концентрации, сорбции и градиента давления представляет собой так называемый фактор сопротивления. Он вычисляется как отношение подвижности воды к подвижности раствора агента. При этом:

Большой объем исследований особенностей моделирования выполнен Б.И. Леви, В.М. Санкиным и В.Р. Баязитовой [2,3]. Известно несколько моделей процесса гелеобразования. Исследования, результаты которых изложены в данной работе, основаны на совершенствовании существующей математической модели за счет учета механизма гелеобразования при изменении концентрации силиката и рН среды, впервые предложенного в работе [4] и позднее развитого в работах В.И. Титова, Н.И. Акимова, Р. Кристенсена и др. Эффективность этого метода подтверждена лабораторными экспериментами и последующими промысловыми испытаниями [1]. Сам процесс закачки силикатно-щелочного раствора происходит следующим образом. В нагнетательную скважину закачивается раствор силиката, по реологическим свойствам близкий закачиваемой воде. По мере прохождения оторочки идет сорбция силиката на породе. Одновременно с этим происходит потеря (нейтрализация) щелочи. При расчете распределения поля концентрации щелочи первоначально считалось, что сорбция линейно зависит от порового объема пласта, через который проходит фильтрация. Однако позже экспериментально Н.И. Акимову, В.И. Титову и др. удалось установить зависимость текущей щелочности среды еще и от текущей концентрации силиката в оторочке. В модели также учитывается нейтрализация щелочи за счет ее растворимости в пластовой воде пропорционально текущей водонасыщенности в каждый момент времени и в каждом узле разностной сетки.

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

Численные исследования, проведенные на математической модели, показали, что основной эффект метода связан с изменением фактора сопротивления как в пространстве, так и во времени.

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

Были выполнены сравнительные расчеты процесса гелеобразования с учетом и без учета диффузионного члена. В расчетах использовалась модель пласта и нефти со следующими свойствами (таблица).

Остальные параметры пластовой нефти следующие: плотность 0,882 г/см3, вязкость – 1,253 мПа·с, объемный коэффициент – 1,150. Эти характеристики соответствуют одному из объектов, на котором применение метода дало положительный эффект, выразившийся в снижении текущей обводненности продукции. Расчеты были выполнены при небольшом временном шаге за период 60 сут.

Из сравнения полей факторов сопротивления с учетом и без учета влияния диффузии на процесс за первые 30 сут (рис. 1) видно, как учет влияния диффузии несколько изменяет прогнозное время начала процесса загеливания. Однако уже через 60 сут учет влияния коэффициента диффузии не меняет времени и места образования зоны загеливания (рис. 2).

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

ВЫВОДЫ

1. Усовершенствованы двух- и трехмерные модели прогноза эффективности применения силикатно-щелочного заводнения с учетом полученных новых экспериментальных зависимостей

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

ЛИТЕРАТУРА

  1. Боровец О.И., Чигина В.В., Титов В.И. Влияние минерализации и жесткости воды на процесс образования геля кремниевой кислоты //НТС Сер Технология и техника для повышения производительности скважин и нефтеотдачи пластов М ВНИИ, 1991 №108 С 44 48
  2. Леви Б.И., Санкин В.М., Баязитова В.Р. Исследование эффективности силикатно-щелочного заводнения водонефтяных зон месторождений пластового типа //Нефт хоз-во 1990 №7 С 45-49
  3. Леви Б.И., Санкин В.М., Баязитова В.Р. Математическое моделирование силикатно-щелочного заводнения с внутрипластовым осадкообразованием Уфа, 1988 Деп во ВНИИОЭШе, № 1509-нг88
  4. Щелочное воздействие на нефтяные пласты и его модификации /А.Т. Горбунов, Ю.В. Желтов, Л.Н. Бученков и др. М , 1985 40 с (Обзор информ /ВНИИОЭНГ Сер Нефтепромысловое дело , Вып 15(104)

ABSTRACT

The article covers the problems concerning the advanced procedures of forecasting apply of silicate-alkaline flooding method and studying the diffusivity effect on helium screen time and site. On the basis of additional experimental data some improvements in the existing mathematical model of efficiency prognosis using silicate-alkaline flooding method have been introduced. The improvements, in first turn, were aimed at considering some features of alkali concentrations fields formation which, to a great extent, control helium screen time and generation site. It was established that in forecasting of effects of this method application in scales of several years, diffusivity factor can be ignored in calculations that will simplify the model.

Номер слоя

Толщина, м

Коэффициент пористости

Проницаемость, мкм2

Коэффициент начальной нефтенасыщенности

Коэффициент ocтаточной нефтенасыщенности

1

20,0

0,330

2,000

0,400

0,260

2

26,0

0,270

0,045

0,800

0,260

3

20,0

0,280

0,325

0,800

0,260

4

15,0

0,260

0,050

0,800

0,260

Рис.1. ПОЛЕ ФАКТОРОВ СОПРОТИВЛЕНИЯ БЕЗ УЧЕТА (А) И С УЧЕТОМ (Б) ВЛИЯНИЯ ДИФФУЗИИ НА ПРОЦЕСС ГЕЛЕОБРАЗОВАНИЯ ЗА ПЕРВЫЕ 30 сут

Рис.2. ПОЛЕ ФАКТОРОВ СОПРОТИВЛЕНИЯ БЕЗ УЧЕТА (А) И С УЧЕТОМ (Б) ВЛИЯНИЯ ДИФФУЗИИ НА ПРОЦЕСС ГЕЛЕОБРАЗОВАНИЯ ЗА 60 сут