Научная статья на тему 'ОПРЕДЕЛЕНИЕ АНИЗОТРОПИИ ПРОНИЦАЕМОСТИ ПО РЕЗУЛЬТАТАМ ГИДРОДИНАМИЧЕСКИХ ИССЛЕДОВАНИЙ СКВАЖИН'

ОПРЕДЕЛЕНИЕ АНИЗОТРОПИИ ПРОНИЦАЕМОСТИ ПО РЕЗУЛЬТАТАМ ГИДРОДИНАМИЧЕСКИХ ИССЛЕДОВАНИЙ СКВАЖИН Текст научной статьи по специальности «Науки о Земле и смежные экологические науки»

CC BY
249
32
i Надоели баннеры? Вы всегда можете отключить рекламу.
Ключевые слова
АНИЗОТРОПИЯ ПЛАСТА / ВЕРТИКАЛЬНАЯ ПРОНИЦАЕМОСТЬ / ГИДРОДИНАМИЧЕСКИЕ ИССЛЕДОВАНИЯ СКВАЖИН / ВЕРТИКАЛЬНОЕ ГИДРОПРОСЛУШИВАНИЕ / 3D ГИДРОПРОСЛУШИВАНИЕ / FORMATION ANISOTROPY / VERTICAL PERMEABILITY / DYNAMIC WELL TESTING / VERTICAL INTERFERENCE TEST / 3D INTERFERENCE TEST

Аннотация научной статьи по наукам о Земле и смежным экологическим наукам, автор научной работы — Аникеев Д.П., Цаган-Манджиев Т.Н.

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

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

Похожие темы научных работ по наукам о Земле и смежным экологическим наукам , автор научной работы — Аникеев Д.П., Цаган-Манджиев Т.Н.

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

DETERMINATION OF PERMEABILITY ANISOTROPY FROM HYDRODYNAMIC WELL TESTING STUDIES

Transition to 3D computer modeling has sharply raised the issue of saturating models with accurate primary information about reservoir properties throughout entire reservoir volume. Fundamental problem of the difference in core and reservoir scales does not allow relying only on core research results. There are a number of examples demonstrating that well data cannot be directly transferred to the entire 3D flow model. This article presents several methods for improved estimating in-situ reservoir permeability, taking into account its anisotropy, from well test data for vertical and horizontal wells. Various methods of carrying out specialized studies aimed at assessing directional strata permeability are delivered.

Текст научной работы на тему «ОПРЕДЕЛЕНИЕ АНИЗОТРОПИИ ПРОНИЦАЕМОСТИ ПО РЕЗУЛЬТАТАМ ГИДРОДИНАМИЧЕСКИХ ИССЛЕДОВАНИЙ СКВАЖИН»

Актуальные проблемы нефти и газа ■ Вып. 2(21) 2018 ■ http://oilgasjournal.ru

ОПРЕДЕЛЕНИЕ АНИЗОТРОПИИ ПРОНИЦАЕМОСТИ ПО РЕЗУЛЬТАТАМ ГИДРОДИНАМИЧЕСКИХ ИССЛЕДОВАНИЙ СКВАЖИН

Д.П. Аникеев, Т.Н. Цаган-Манджиев ИПНГ РАН, e-mail: anikeev@ogri.ru

Введение

Переход к 3D компьютерному моделированию остро поставил вопрос насыщения модели точной первичной информацией о коллекторских свойствах пласта.

Фундаментальная проблема различия в масштабах керна и пласта не позволяет опираться на результаты керновых исследований без их вдумчивого анализа. Характерный размер кернов составляет сантиметры, а типичный размер стороны сеточной ячейки в полномасштабной гидродинамической модели может лежать от 25 м до 100 м в плане. Поэтому, например, неоднородность коллектора с характерным размером в метры не будет «чувствоваться» на керне, но в то же время будет влиять на достоверность гидродинамической модели. Можно привести и другие примеры, демонстрирующие, что скважинные данные невозможно напрямую перенести на весь объем 3D гидродинамической модели.

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

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

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

• Этап I. На этом этапе скважина эксплуатируется с целью добычи или закачки флюида в течение некоторого периода времени.

• Этап II. Скважина закрывается для общепринятого снятия кривой восстановления забойного давления (КВД) или кривой падения давления (КПД) в зависимости от режима

на первом этапе. При этом не формулируется требование обязательного восстановления забойного давления до уровня начального пластового.

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

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

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

Технологии исследования скважин по определению параметров пласта

в пластовых условиях

В данном типе задач параметры пласта определяются на основе решения обратной задачи. Алгоритм их решения предусматривает последовательное выполнение следующих этапов:

• Решение прямой фильтрационной задачи.

Целью решения прямой задачи является получение динамик изменения замеряемых показателей в процессе проведения исследований при заданных текущих значениях управляющих параметров.

Для решения нелинейной прямой задачи применяется неявная двухслойная по времени конечно-разностная (конечно-объемная) численная схема. На каждом временном слое возникает следующая нелинейная система алгебраических уравнений:

Ру(ху,ху~\и) = О у = 1,2,...,ТУ, (1)

где х1 - вектор независимых переменных прямой задачи (фазовых переменных) на у'-ом временном слое, N - номер последнего шага прямой задачи. В качестве х1 при

двухфазной фильтрации выступает вектор значений давления в нефтяной фазе и значений водонасыщенности в каждом сеточном блоке на временном слое у'. Через и обозначен вектор идентифицируемых параметров (параметров пласта, включая ОФП), которые далее носят название управляющих параметров.

Граничные условия соответствуют отсутствию потока через внешнюю границу, а также специальному условию, выбираемому в зависимости от режима эксплуатации, по потоку в/из вскрытых перфорацией ячеек сеточной модели [2]. Начальные условия соответствуют решению специальной задачи гравитационно-капиллярного равновесия.

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

(хт) 8хт = -Р}(х]<у)) , (2)

где V - номер текущей итерации метода Ньютона-Рафсона, через дх'<л> обозначена величина очередной добавки к решению нелинейной системы (1), а Р^ - матрица частных производных (матрица Якоби) уравнений (1) для у'-ого временного слоя прямой задачи по фазовым переменным х1.

Таким образом, в результате решения прямой задачи определяются значения фазовых переменных х1 (давление в нефтяной фазе и водонасыщенности) для всех узлов сетки и на всех временных слоях 7 = 1,2,...,N при текущем наборе значений управляющих параметров й.

• Вычисление функционала при полученных значениях замеряемых показателей.

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

Функционал качества в векторном виде представлен следующим выражением:

N т

J{й) = YJ(yJ{xJ{й\й)-YJ) а(у](х\й\й)-У]), (3)

7=1

где7 обозначает номер замера на скважине, N - число замеров за время исследования, у' и У7 - соответственно векторы расчетных и фактически замеренных значений показателей

работы скважины на момент у'-ого замера; Т - знак транспонирования. В качестве компонентов векторов у7 и У может выступать некоторый набор из замеряемых значений, множество которых зависит от конкретной задачи. Диагональная матрица О обеспечивает нормировку разноразмерных измерений и задает удельное влияние каждого из замеров на значение критерия качества. В частности, величинам, замеряемым с меньшей погрешностью, должны соответствовать большие значения соответствующих элементов на диагонали О.

• Для определения производных функционала качества по управляющим параметрам находится решение дополнительной сопряженной краевой задачи, формулировка которой следует из методологии теории оптимального управления [2].

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

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

Рассмотренные технологии были протестированы на синтетических примерах. Процедура тестирования традиционна для данного класса обратных задач. Задаются некоторые значения искомых (управляющих) параметров, которые считают «истинными». На их основе решается прямая задача. Полученные динамики показателей работы скважины (забойного давления, дебита по нефти, водонасыщенности околоскважинной зоны) воспринимаются в качестве замеренных при промысловом исследовании (или дополнительно «зашумляются» путем добавления случайных величин с заданным математическим распределением). Они используются как исходные данные для решения обратной задачи. При этом «истинные» значения параметров «забываются». Затем задается произвольное начальное приближение и реализуется итерационный процесс. Тогда появляется возможность сопоставить полученные в результате работы

алгоритма величины искомых параметров с их «истинными» значениями. В результате такой процедуры удается протестировать реализацию алгоритма решения обратной задачи, а также отчасти преодолеть неизбежное свойство обратных задач -некорректность в смысле существования, единственности и устойчивости решения. Созданные технологии отразились в полученных патентах на новые способы исследования скважин и пластов [3, 4].

Технология вертикального гидропрослушивания

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

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

Для описанной процедуры исследования в случае возбуждения за счет добычи нефти прямая задача состоит в моделировании однофазной фильтрации в осесимметричной 2D постановке с учетом различия значений проницаемости слоев в горизонтальном (радиальном) кт и вертикальном к2 направлениях [5]:

т от I от) ог V ог) ст ог

Здесь кг и к2 - коэффициенты эффективной проницаемости (проницаемости для нефти при остаточной водонасыщенности) в горизонтальном и вертикальном направлениях, / - коэффициент вязкости флюида (нефти), 0* = 0 + т00ж - коэффициент упругоемкости пласта, т0 - коэффициент эффективной пористости пласта при начальном

давлении, вж и вс - коэффициенты сжимаемости (упругости) жидкости (нефти) и эффективной пористой среды (с учетом остаточной водонасыщенности) соответственно [2], р0 - плотность нефти при начальном давлении.

Процедура моделирования скважины описана в статье*. Граничные условия к уравнению (4) соответствуют непротеканию флюида через верхнюю и нижнюю границы пласта и через условную круговую внешнюю границу, а также заданному дебиту (расходу) на активном интервале скважины. Начальное условие предполагает заданное (например, невозмущенное гидростатическое) распределение давления в пределах моделируемого участка на начало исследования.

На основе реализованного численного алгоритма решения прямой задачи (с учетом логарифмического шага разностной сетки по оси г) проведены расчеты по определению зоны чувствительности двух вариантов исследования - вертикального гидропрослушивания по схеме, представленной на рис. 1, и известного исследования с несколькими точками замера давления в открытом стволе с использованием пластоиспытателя. В последнем случае разнесение активного и реагирующих интервалов по высоте не превышает 1,5-2 м, а объем отбираемого флюида сильно ограничен.

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

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

* См. статью Закиров Э.С., Закиров С.Н., Индрупский И.М., Аникеев Д.П. «О представлении скважины в 3Б гидродинамической модели» в данном выпуске.

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

Таблица 1

Исходные данные сравнительного примера

Упругоемкость пласта, атм-1 1,45^10-4 Эффективная пористость пласта при нач. давл., доли ед. 0,2

Радиус скважины, м 0,1 Нач. пласт.давление, атм 344,5

Гориз.проницаемость кг, мД 20 Толщина пласта к, м 21,5

Верт. проницаемость к2, мД 5 Скин-фактор, доли ед. 1

Дебит, м3/сут 1,6 Вязкость нефти, сПз 1

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

Функционал качества в задаче вертикального гидропрослушивания имеет вид:

М ■ Л 9 М -А, -9

Т(и) =у С (Р]ф -Рр )2 + У С (Р]ф -Р]'р )2 (5)

^ V / акт\ акт акт / ¿¿^реаг\ реаг реаг / ~> \ /

где М - количество замеров давления, Рj - давление в точке замера на момент времени у, верхние индексы ф и р отвечают фактическим и расчетным величинам, а нижние индексы акт и реаг - активному и реагирующему интервалу соответственно, Сакт и Среаг -весовые коэффициенты, и - вектор идентифицируемых (управляющих) параметров, представленный значениями кг и к2 нескольких выделенных зон пласта, включая скин-зону.

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

проницаемости) при общей протяженности пласта 500 м. На рис. 3 сопоставлены расчетные (для начального приближения и при найденных значениях параметров) и «фактические» динамики изменения давления на реагирующем интервале. Аналогичная степень совпадения имеет место и на активном интервале.

Таблица 2

Оцениваемые параметры и их значения

Проницаемость, мДарси Фактическое значение До идентификации После идентификации

3 10 3

кг1 10 30 10

кг2 15 35 15

кю (не идентиф.) 2 2 2

кл 5 10 5

кп (не идентиф.) 10 10 10

1, час

Рис. 3. Динамики давления на реагирующем интервале в рассматриваемом примере

Технология 3Б гидропрослушивания на основе горизонтальных скважин

Одним из следствий концепции эффективного порового пространства (ЭПП) явилась технология вертикально-латерального заводнения [6]. Потребности ее проектирования и реализации привели к обоснованию метода 3D гидропрослушивания на основе горизонтальных скважин для уточнения степени сообщаемости пласта вдоль вертикальной координаты. Схема элемента разработки и соответствующего исследования приведена на рис. 4.

Технология 3Б гидропрослушивания апробирована на одном из месторождений Западной Сибири. Опытный участок (см. рис. 4) состоял из двух горизонтальных скважин -добывающей (реагирующей) у кровли пласта и нагнетательной (возбуждающей) - у подошвы пласта [7].

вид в плане

вид в разрезе

И- ~ 300 м -Н

Рис. 4. Схема 3Б гидропрослушивания и вертикально-латерального заводнения

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

Для интерпретации результатов гидропрослушивания создана 3D секторная гидродинамическая модель участка. Из-за отсутствия надежной 3D геологической модели

500 м

низкопроницаемый ««*. прослой

т

300 м

пласта использовалась слоисто-неоднородная модель. Параметры слоев модели (включая начальные значения проницаемости) приняты по результатам ГИС в пилотном стволе добывающей скважины.

Согласно представлениям концепции ЭПП [6], все пропластки в разрезе модели являются проницаемыми. Для прослоев, отнесенных по ГИС к неколлекторам, в модели принято начальное значение горизонтальной проницаемости 0,1 мД. Расчеты на секторной гидродинамической модели выполнены в двухфазной (нефть-вода) постановке из-за высокой обводненности добываемой продукции на данном участке вследствие предшествующей разработки.

В качестве идентифицируемых параметров рассматривались два параметра. А именно - единый множитель для проницаемости всех слоев секторной модели во всех координатных направлениях, а также множитель для вертикальной проницаемости (по вариантам - для всех слоев модели или только для низкопроницаемых (бывших «неколлекторов»), что сказывается на интегральной характеристике анизотропии проницаемости модели).

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

192 190

Я /

Я /

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

аТ '

В 188 /

к ^

С /

в '

я

ч

о 186

В

182 --

1600 1800 2000 Время с начала исследования, ч

Рис. 5. Расчетная и фактическая динамики давления в реагирующей скважине

Средневзвешенные (с учетом толщин слоев) значения проницаемости составили в горизонтальном направлении 12,9 мД в пересчете на эффективную проницаемость (0,7 мД для проницаемости по воде при остаточной нефтенасыщенности), в вертикальном направлении - 1,6 мД (0,08 мД) соответственно; соотношение горизонтальной и вертикальной проницаемости составило 8,1.

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

Технология ЭБ гидропрослушивания на основе вертикальных скважин Схема размещения возбуждающей и реагирующих скважин, представленная на рис. 6, является примером соответствующей технологии 3Б гидропрослушивания [6, 8, 9]. Здесь реагирующие скважины в плане расположены «веером» вокруг нагнетательной, и их интервалы вскрытия пласта относятся к разным глубинным отметкам.

Рис. 6. Схема 3Б гидропрослушивания: скв. 7 - возбуждающая, 1-6 - реагирующие Решение прямой задачи осуществляется численно, а ее постановка формулируется для 3Б однофазного случая, с учетом заполненной матрицы коэффициентов тензора абсолютной проницаемости К:

д{ш(Р)р(Р))

М Р)

к(ур+р(Р)gVh) = 4 + 6(0, (6)

дг

где учитывается зависимость плотности и вязкости нефти и эффективной пористости р, ^, т от давления Р, И - высотная отметка данной точки пласта относительно

условного уровня, 6(г) - интенсивность источников/стоков, связанных с работой возбуждающей скважины (скважин). Начальные и граничные условия отражают заданное

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

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

/ \ м к -9

J(и)=Е -р> )2 . (7)

7=1 №=1 Р ф

Здесь у - порядковый номер соответствующего момента измерений, М - количество таких моментов, N - количество скважин, по которым осуществляются замеры давления,

с - весовые коэффициенты, и - вектор идентифицируемых параметров, р , р -

^ р ^ ф

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

Результаты тестирования алгоритма решения обратной задачи на синтетических примерах показали возможность успешной идентификации всех 6 компонент заполненной симметричной матрицы тензора проницаемости совместно с емкостными параметрами коллектора [8].

Промысловый эксперимент по 3D гидропрослушиванию проведен на Памятно-Сасовском месторождении Волгоградской области. Продуктивный пласт сложен рифогенным карбонатным коллектором и характеризуется значительными нефтенасыщенными толщинами. Особенностью исследования явилось использование законтурной нагнетательной скважины (14-ПС) в качестве возбуждающей и нескольких добывающих скважин (126, 2, 121 и 13) в роли реагирующих [6, 8, 9]. Схема участка проведения эксперимента представлена на рис. 7. На рис. 8 дается профильный разрез пласта с отмеченными интервалами вскрытия его скважинами.

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

Рис. 7. Границы участка залежи (обозначены красным цветом) и направления координатных осей, принятые при интерпретации результатов 3Б гидропрослушивания

Начальный

Рис. 8. Схема вскрытия пласта скважинами

Математические эксперименты

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

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

- в каждой из зон. При этом рассматривались и другие варианты разбиения 3Б секторной модели на зоны.

На рис. 9 и 10 в качестве примера приводятся фактические и расчетные динамики забойного давления в одной из реагирующих и возбуждающей скважинах. В результате итерационной процедуры по уточнению коллекторских свойств пласта достигнута удовлетворительная степень совпадения расчетных и фактических динамик забойных давлений. При этом конечная величина функционала качества (7) оказалась почти в 130 раз меньшей начальной. Наиболее близко от замеренных значений отстоят расчетные забойные давления по скв. 126. Наиболее удалены от истинных значений давления по скв. 121 и 14ПС, что можно объяснить неучетом фактора многофазности. Имеющиеся расхождения между фактическими и садаптированными значениями забойных давлений связаны с недостаточной степенью достоверности исходной геологической модели объекта.

223.2 ■

222.7 ■

а я

о

§ 222.2 ■

я

ч

о о

Я

«

ю 221.7 ■

я РО

221.2

♦ ♦ ♦фактические замеры - - - до идентификации - после идентификации

320

н 310

я

аГ а я

о Ч в

?! зоо

ш о

Я

«

о »о я

м 290

о

280

♦ ♦ ♦фактические замеры - - - до идентификации

после идентификации

0

Рис. 9

0.5

2

1 1.5

Время, сут

Динамики забойных давлений для реагирующей скважины 126

2.5

0

10

30

20

Время, сут

Рис. 10. Динамики забойных давлений для возбуждающей скважины 14-ПС

40

В результате решения обратной задачи получены результаты, представленные в табл.3.

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

Таблица 3

Значения управляющих параметров

Управляющие параметры До идентификации После идентификации

Зона 1 (скв. 2) 1080 140,21

Зона 2 (скв. 13) 4880 127,58

К, мД Зона 3 (скв. 121) 3590 183,88

Зона 4 (скв. 12б) 1540 397,54

Зона 5 (скв. 14-ПС) 23 1571,88

Зона 1 (скв. 2) 1080 34,11

Зона 2 (скв. 13) 4880 618,98

куу, мД Зона 3 (скв. 121) 3590 47,89

Зона 4 (скв. 12б) 1540 241,27

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

Зона 5 (скв. 14-ПС) 23 2,67

Зона 1 (скв. 2) 108 1424,04

Зона 2 (скв. 13) 488 34,83

К, мД Зона 3 (скв. 121) 359 5321,59

Зона 4 (скв. 12б) 154 1516,42

Зона 5 (скв. 14-ПС) 2,3 770,59

Зона 1 (скв. 2) 0,1 0,123

Зона 2 (скв. 13) 0,15 0,0098

то , д. ед. Зона 3 (скв. 121) 0,1 0,4

Зона 4 (скв. 12б) 0,1 0,014

Зона 5 (скв. 14-ПС) 0,15 0,01

Использование подхода на основе решения обратных задач позволило решить и некоторые другие задачи. Так в статьях [10, 11] представлены результаты исследований, проведенных в условиях аномально низкой приемистости пласта.

Статья написана в рамках выполнения государственного задания (тема «Научное обоснование новых экологически чистых технологий разработки месторождений углеводородов в сложных горно-геологических условиях на основе 3D-компьютерных экспериментов», № АААА-А16-116022510270-1).

ЛИТЕРАТУРА

1. Закиров Э.С., Тарасов А.И., Индрупский И.М. Новый подход к исследованию газовых скважин и интерпретации получаемых результатов // Газовая промышленность. 2003. № 9. С. 61-63.

2. Закиров Э.С. Трехмерные многофазные задачи прогнозирования, анализа и регулирования разработки месторождений нефти и газа // М.: Изд. дом «Грааль», 2001. 302 с.

3. Закиров С.Н., Закиров Э.С. Пат. RU 2258137 С1. Способ исследования и интерпретации результатов исследований скважины, вскрывшей нефтегазовую залежь

№ 2004119525/03; Заявл. 29.06.2004; Опубл. 10.08.2005 // Изобретения. Полезные модели. 2005. Бюл. № 22. - Режим доступа: http://www1.fips.ru

4. Закиров С.Н., Закиров Э.С., Индрупский И.М., Аникеев Д.П. Пат. RU 2213864 C2. Способ исследования и интерпретации результатов исследований скважин и пластов № 2001132910/03; Заявл. 06.12.2001; Опубл. 10.10.2003 // Изобретения. Полезные модели. 2005. Бюл. № 28. - Режим доступа: http://www1.fips.ru

5. Лейбензон Л.С. Движение природных жидкостей и газов в пористой среде. М.-Л.: Гостехиздат, 1947. 244 с.

6. Закиров С.Н., Индрупский И.М., Закиров Э.С., Закиров И.С. и др. Новые принципы и технологии разработки месторождений нефти и газа. Часть 2. М.-Ижевск: Институт компьютерных исследований, 2009. 484 с.

7. Индрупский И.М. Интерпретация результатов 3D гидропрослушивания на основе секторного моделирования // Геология, геофизика и разработка нефтяных и газовых месторождений. 2008. № 11. С. 37-41.

8. Брадулина О.В. 3D гидропрослушивание на Памятно-Сасовском месторождении // Геология, геофизика и разработка нефтяных и газовых месторождений. 2008. № 12. С. 42-46.

9. Закиров С.Н., Закиров Э.С., Индрупский И.М., Левченко В.С., Брадулина О.В., Цаган-Манджиев Т.Н. Вертикальное и 3D гидропрослушивание продуктивных пластов // Новые технологии освоения и разработки трудноизвлекаемых запасов нефти и газа и повышения нефтегазоотдачи: Сб. тр. VII Междунар. технологич. симпозиума. М., 2008. С. 49-63.

10. Закиров Э.С., Индрупский И.М., Васильев И.В., Аникеев Д.П., Цаган-Манджиев Т.Н., Родионов А.Е., Лачугин Д.С., Афанасьев В.С., Афанасьев С.В., Антонович А.А. Проведение комплексных исследований по оценке относительных фазовых проницаемостей для нефти и воды и коэффициента вытеснения в условиях аномально низкой приемистости пласта (часть 1) // Нефтяное хозяйство. 2016. № 9. С. 56-60.

11. Закиров Э.С., Индрупский И.М., Васильев И.В., Аникеев Д.П., Цаган-Манджиев Т.Н., Родионов А.Е., Лачугин Д.С., Афанасьев В.С., Афанасьев С.В., Антонович А.А. Проведение комплексных исследований по оценке относительных фазовых проницаемостей для нефти и воды и коэффициента вытеснения в условиях аномально низкой приемистости пласта (часть 2) // Нефтяное хозяйство. 2017. № 10. С. 90-93.

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