Научная статья на тему 'Исследование процесса нефтеизвлечения в коллекторах различного типа пустотности с использованием многомерного статистического анализа'

Исследование процесса нефтеизвлечения в коллекторах различного типа пустотности с использованием многомерного статистического анализа Текст научной статьи по специальности «Науки о Земле и смежные экологические науки»

CC BY
281
41
i Надоели баннеры? Вы всегда можете отключить рекламу.
Ключевые слова
СЛОЖНОПОСТРОЕННАЯ ЗАЛЕЖЬ / КОЛЛЕКТОР ТРЕЩИННОГО ТИПА / КОЛЛЕКТОР ПОРОВОГО ТИПА / ЗАКОНОМЕРНОСТИ НЕФТЕИЗВЛЕЧЕНИЯ / ЭЛЕМЕНТ СИСТЕМЫ РАЗРАБОТКИ / ОЦЕНКА ВЗАИМОДЕЙСТВИЯ МЕЖДУ СКВАЖИНАМИ / КОРРЕЛЯЦИЯ / КОРРЕЛЯЦИОННЫЕ ПОЛЯ / МНОГОМЕРНЫЙ РЕГРЕССИОННЫЙ АНАЛИЗ / МАТЕМАТИЧЕСКИЕ МОДЕЛИ / ФАКТОРЫ / ФОРМИРУЮЩИЕ ДЕБИТ НЕФТИ В РАЗЛИЧНЫХ КОЛЛЕКТОРАХ / ЛИНЕЙНЫЕ ДИСКРИМИНАНТНЫЕ ФУНКЦИИ / ДОСТОВЕРНОСТЬ РАСПОЗНАВАНИЯ КОЛЛЕКТОРОВ / СЦЕНАРИИ ФОРМИРОВАНИЯ ДЕБИТА НЕФТИ / ЭФФЕКТИВНОСТЬ СИСТЕМЫ ПОДДЕРЖАНИЯ ПЛАСТОВОГО ДАВЛЕНИЯ / COMPLEX ACCUMULATION / NATURALLY FRACTURED RESERVOIRS / PORE TYPE RESERVOIR / OIL PRODUCTION LAW / COMPONENTS OF DEVELOPMENT SCHEME / WELL INTERACTION ANALYSIS / CORRELATION / CORRELATION FIELDS / MULTIDIMENSIONAL REGRESSION ANALYSIS / MATHEMATICAL MODELS / FACTORS FORMING OIL RATE IN DIFFERENT RESERVOIRS / LINEAR DISCRIMINANT FUNCTIONS / AUTHENTICITY OF RESERVOIR DETECTION / OIL RATE FORMATION SCENARIOS / EFFICIENCY OF PRESSURE MAINTENANCE SYSTEM

Аннотация научной статьи по наукам о Земле и смежным экологическим наукам, автор научной работы — Галкин Владислав Игнатьевич, Пономарева Инна Николаевна, Репина Вера Андреевна

Изучены вопросы установления закономерностей осуществления процессов нефтеизвлечения, учитывающих распространение в пределах одной залежи коллекторов с различными видами пустотности. С этой целью привлечены промысловые материалы по разработке турне-фаменской карбонатной залежи нефтяного месторождения, а также для сравнения данные о разработке залежи нефти со схожими свойствами нефти, но с терригенным гранулярным коллектором порового типа. Для рассмотрения в статье приняты элементы реализованных систем разработки, представляющих собой одну нагнетательную и соседние добывающие скважины. Тип коллектора в пределах элемента разработки принимался по материалам различных исследований, в том числе по кривой восстановления давления, обработанной по методике Уоррена-Рута. На первом этапе решения поставленной задачи рассчитаны значения коэффициентов корреляции между приемистостью нагнетательной скважины и дебитами соседних добывающих для различных моментов времени, при этом расчет производился в предположении, что коэффициент корреляции является численной мерой взаимодействия между двумя скважинами. Установлено, что изменение коэффициента корреляции во времени для поровых коллекторов значительно отличается от характера его же поведения для коллектора трещинного типа. В дальнейшем проведенные исследования позволили получить многомерные математические модели, характеризующие процесс заводнения и позволяющие определять дебиты добывающих скважин, расположенных вблизи очагов нагнетания, с учетом типа пустотности коллектора на рассматриваемом участке залежи. Анализ линейных дискриминантных функций, построенных с учетом типа пустотности коллектора, позволил установить, что процесс вытеснения нефти водой в терригенных поровых, карбонатных поровых и карбонатных трещинных коллекторах происходит по различным сценариям.

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

Похожие темы научных работ по наукам о Земле и смежным экологическим наукам , автор научной работы — Галкин Владислав Игнатьевич, Пономарева Инна Николаевна, Репина Вера Андреевна

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

Study of Oil Recovery from Reservoirs of Different Void Types with Use of Multidimensional Statistical Analysis

Oil recovery laws that take into account distribution of reservoirs with different void types within the same accumulation. Carbonate field data of development of Tournaisian-Famennian oil accumulation was used. For comparison purposes data of development of the field with similar oil properties but reservoir of clastic grain rock and pore type were used. One injector and neighbor producers were used as components of applied production scheme. The type of reservoir within one development object was determined by several studies including pressure build-up curve processed by Warren-Root method. At the first stage correlation coefficients between injection and production of neighbor well were calculated. Calculation was done for different time and with assumption, according to which correlation coefficient is a quantitative measure of interactions between two wells. It is determined that use of correlation coefficient for pore reservoirs is significantly differ to the character of its behavior for fracture reservoir type. Multidimensional mathematical models that characterize flooding and allow determining producer’s rate were obtained with considered void type. Linear discriminant functions are built with considered void type of reservoir. Analysis of those functions determined that replacement of oil by water in clastic and carbonates porous rocks and carbonate naturally fractured reservoirs follow different scenarios.

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

ВЕСТНИК ПНИПУ. ГЕОЛОГИЯ. НЕФТЕГАЗОВОЕ И ГОРНОЕ ДЕЛО BULLETIN OF PNRPU. GEOLOGY. OIL & GAS ENGINEERING & MINING

ISSN 2224-9923

Том / Volume 15 №19 2016

http://vestnik.pstu.ru/geo/

УДК 553.982.2

Статья / Article © ПНИПУ / PNRPU, 2016

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

В.И. Галкин, И.Н. Пономарева, В.А. Репина1

Пермский национальный исследовательский политехнический университет (614990, Россия, г. Пермь, Комсомольский пр., 29) 1Филиал ООО «ЛУКОЙЛ-Инжиниринг» «ПермНИПИнефть» в г. Перми (614066, Россия, г. Пермь, ул. Советской Армии, 29)

STUDY OF OIL RECOVERY FROM RESERVOIRS OF DIFFERENT VOID TYPES WITH USE OF MULTIDIMENSIONAL STATISTICAL ANALYSIS

V.I. Galkin, I.N. Ponomareva, V.A. Repina1

Perm National Research Polytechnic University (29 Komsomolskii av., Perm, 614990, Russian Federation) 1PermNIPIneft branch of LUKOIL-Engineering LLC in Perm (29 Sovetskoi Armii st., Perm, 614066, Russian Federation)

Получена / Received: 01.12.2015. Принята / Accepted: 27.04.2016. Опубликована / Published: 30.06.2016

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

Изучены вопросы установления закономерностей осуществления процессов нефтеизвлечения, учитывающих распространение в пределах одной залежи коллекторов с различными видами пустотности. С этой целью привлечены промысловые материалы по разработке турне-фаменской карбонатной залежи нефтяного месторождения, а также для сравнения - данные о разработке залежи нефти со схожими свойствами нефти, но с терригенным гранулярным коллектором порового типа. Для рассмотрения в статье приняты элементы реализованных систем разработки, представляющих собой одну нагнетательную и соседние добывающие скважины. Тип коллектора в пределах элемента разработки принимался по материалам различных исследований, в том числе по кривой восстановления давления, обработанной по методике Уоррена-Рута. На первом этапе решения поставленной задачи рассчитаны значения коэффициентов корреляции между приемистостью нагнетательной скважины и дебитами соседних добывающих для различных моментов времени, при этом расчет производился в предположении, что коэффициент корреляции является численной мерой взаимодействия между двумя скважинами. Установлено, что изменение коэффициента корреляции во времени для поровых коллекторов значительно отличается от характера его же поведения для коллектора трещинного типа. В дальнейшем проведенные исследования позволили получить многомерные математические модели, характеризующие процесс заводнения и позволяющие определять дебиты добывающих скважин, расположенных вблизи очагов нагнетания, с учетом типа пустотности коллектора на рассматриваемом участке залежи. Анализ линейных дискриминантных функций, построенных с учетом типа пустот-ности коллектора, позволил установить, что процесс вытеснения нефти водой в терригенных поровых, карбонатных поровых и карбонатных трещинных коллекторах происходит по различным сценариям.

Key words:

complex accumulation, naturally fractured reservoirs, pore type reservoir, oil production law, components of development scheme, well interaction analysis, correlation, correlation fields, multidimensional regression analysis, mathematical models, factors forming oil rate in different reservoirs, linear discriminant functions, authenticity of reservoir detection, oil rate formation scenarios, efficiency of pressure maintenance system

Oil recovery laws that take into account distribution of reservoirs with different void types within the same accumulation. Carbonate field data of development of Tournaisian-Famennian oil accumulation was used. For comparison purposes data of development of the field with similar oil properties but reservoir of clastic grain rock and pore type were used. One injector and neighbor producers were used as components of applied production scheme. The type of reservoir within one development object was determined by several studies including pressure build-up curve processed by Warren-Root method. At the first stage correlation coefficients between injection and production of neighbor well were calculated. Calculation was done for different time and with assumption, according to which correlation coefficient is a quantitative measure of interactions between two wells. It is determined that use of correlation coefficient for pore reservoirs is significantly differ to the character of its behavior for fracture reservoir type. Multidimensional mathematical models that characterize flooding and allow determining producer's rate were obtained with considered void type. Linear discriminant functions are built with considered void type of reservoir. Analysis of those functions determined that replacement of oil by water in clastic and carbonates porous rocks and carbonate naturally fractured reservoirs follow different scenarios.

Галкин Владислав Игнатьевич - доктор геолого-минералогических наук, профессор, заведующий кафедрой геологии нефти и газа (тел.: +007 219 80 17, e-mail: [email protected]).

Пономарёва Инна Николаевна - кандидат технических наук, доцент, доцент кафедры нефтегазовых технологий (моб. тел.: +007 919 483 84 30, e-mail: [email protected]). Контактное лицо для переписки.

Репина Вера Андреевна - инженер 1-й категории отдела гидродинамического моделирования (моб. тел.: +007 902 809 63 88, e-mail: [email protected]).

Vladislav I. Galkin (Author ID in Scopus: 55418067700) - Doctor of Geological and Mineralogical Sciences, Professor, Head of the Department at the Department of Oil and Gas Geology (tel.: +007 219 80 17, e-mail: [email protected]).

Inna N. Ponomareva (Author ID in Scopus: 36242312300) - PhD in Engineering, Associate Professor, Associate Professor at the Department of Oil and Gas Technologies (mob. tel.: +007 919 483 84 30, e-mail: [email protected]). The contact person for correspondence.

Vera A. Repina - Engineer of 1st category of Hydrodynamic Modeling Department (mob. tel.: +007 902 809 63 88, e-mail: [email protected]).

Введение

Характерной особенностью многих турне-фаменских залежей нефти на территории Пермского края является наличие зоны распространения коллекторов трещинно-порового типа. Данный вывод получен по результатам различных исследований [1, 2], а в работах [3, 4] показаны особенности распространения коллекторов такого типа. На основе материалов обработки кривых восстановления давления в соответствии с методикой Уоррена-Рута, актуальность и достоверность которой всесторонне исследована и подтверждена в работах [5-12], установлено, что обычно зоны трещиноватых коллекторов направлены с юго-запада на северо-восток. Наличие данных зон во многом осложняет не только геологическое строение залежей, но и процессы их разработки. Актуальным представляется вопрос оценки эффективности реализованной на залежи системы поддержания пластового давления и выявление процессов вытеснения нефти водой в зависимости от типа пустотности коллектора. Значительный накопленный объем промысловых данных позволяет решать поставленную задачу методами математической статистики [13].

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

Построение и анализ математических моделей взаимодействия добывающих и нагнетательных скважин

По всем выбранным для исследования элементам разработки привлечены ежемесячные данные по закачке воды (УН 0 ), приемистости нагнетательных (Рг) и дебитам нефти ^т) добывающих скважин.

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

скв. 1 - Рг = 23,5 + 0,0293 ГН0 при г = 0,92;

скв. 2 - Рг = 34,6 + 0,0228 ГНг0 при г = 0,87;

скв. 3 - Рг = 12,8 + 0,0301 ГН0 при г = 0,90.

Отсюда видно, что величины 0 и Рг достаточно хорошо коррелируют. Для анализа также были привлечены накопленные значения приемистости Рг". В табл. 1 приведены модели изменения величины Рг" во времени (£).

Таблица 1

Зависимости PГ от t

№ скв. п/п Уравнение регрессии Коэффициент R2

1 Рг" = 242,8 - 0,0006^ + 0,001г4 - 0,187^ + + 11,04г2 - 57,74г 0,99

2 Рг" = -916,1 + 0,005^ - 1,884Г2 + 295,4Г 0,99

3 Рг" = 16,26 + 0,402^ + 84,06г 0,99

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

Для установления влияния закачки объемов воды в пласт на среднесуточные дебиты нефти в добывающих скважинах Qт высчитывались коэффициенты корреляции г между показателями УН0 , Рг, Рг" и Qт. Предполагается, что по значениям г можно выполнить оценку влияния объема закачиваемой воды в нагнетательную скважину на дебиты нефти добывающих скважин. На начальном этапе определялись значения г за период, равный 3 месяцам с начала закачки во-

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

Для скважины 1 анализ выполняется по 5 ближайшим добывающим скважинам, для которых характерен трещинный тип коллектора. Для нагнетательной скважины 2 использовались 7 ближайших добывающих скважин с поровым типом коллектора. Для скважины 3 использовались данные по ближайшим 4 добывающим скважинам с поровым типом коллектора, породы представлены терригенными отложениями. Расстояния от нагнетательных до добывающих скважин изменяются от 334 до 960 м. Максимальные различия г в процессе изменения во времени наблюдается при г(Ргп - )т). Поэтому приведем типовые кривые изменения значений г(Ргп - )т) для 3 добывающих скважин (рис. 1).

Из анализа представленной на рис. 1 диаграммы следует, что изменения значений г между Ргп и )т - г(Ргп - )т) для изучаемых скважин во всех случаях индивидуальны. При этом необходимо отметить, что различия в виде кривых изменения значений г(Ргп - )т) во времени для поровых карбонатных и поровых терригенных коллекторов существенно меньшие, чем различия в виде данных кривых, характеризующих поровые и трещинные коллекторы. Аналогичные вычисления г выполнены между УНо и )т - г(Рг - )т) и между Гя2 о и Рг - ^ - г(Рг - &).

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

ом =-124,337г (КН2о - Рг) -1,719г(Рг - ) + 57,033г(Ргп - ) +

Рис. 1. Изменение коэффициента г между Ргп и )т во времени по добывающим скважинам

+ 182,540

при Я = 0,85, р < 0,0000, ошибка прогноза составляет 6,1 т/сут. По данной формуле были вычислены значения . Графическое сопоставление значений )т и во времени представлено на рис. 2, а.

Из анализа рис. 2, а следует, что в пределах графика выделяются 4 участка. Первый участок расположен в диапазоне от 0 до 42 месяцев, здесь наблюдается хорошее совпадение значений ) . Для второго участка при 43 < ^ < 84 мес. наблюдается превышение ) над )т. Далее значения )т и 0тм слабо «контролируют» друг друга, но в их пределах выделяются еще два подпо-ля. Третий участок располагается в диапазоне 85-112 месяцев; четвертый - в интервале 113-138 месяцев. Для учета соотношений )т и 0м в пределах выделенных временных соотношений, наблюдаемых на рис. 2, а, построим многомерные модели.

Построение и анализ временных многомерных математических моделей взаимодействия добывающих и нагнетательных скважин

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

Из данных табл. 2 видно, что формирование

^ /лмв

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

По формулам были вычислены значения О по 138 данным.

Сопоставление значений О и О во времени приведено на рис. 2, б. Анализ показывает, что в целом значения О и О совпадают достаточно

хорошо. Здесь необходимо отметить, что во временном диапазоне 85-112 месяцев значения О и слабо «контролируют» друг друга, что количественно подтверждается значением г, приведенным в табл. 2. На основании данных расчетов, выполненных по общей и временным моделям, по данной добывающей скважине построим поля корреляции между О и О , О и О* (рис. 3).

Корреляционные поля между значениям О , вычисленными по общей модели, и значениями Отмв, вычисленными по временным моделям с фактическими дебитами От, несколько отличаются по силе корреляционных связей (см. рис. 3). В первом случае г = 0,86, во втором г = 0,98.

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

: a ......

□ о :.............. о : ^ ......

ftLo

Jpg

;........i.........Ш........i..... 31 n>

О oö

-20 0 20 40 60 80 100 120 140 160 Время, месяц

а

-20 0 20 40 60 80 100 120 140 160 Время, месяц

б

о Фактические значения дебитов нефти □ Модельные значения дебитов нефти Рис. 2. Изменение О и О во времени: а - многомерная модель; б - многомерная временная модель

Таблица 2

Модели для вычисления QM

г

Интервал времени, Коэффициенты при показателях Р Свободный член R

месяц Р Р

r (VH2O - Pr) r(Pr - QT) r(Prn - QT)

0-42 -81,2298 6,5303 31,6133 144,339 0,871

0,021634 0,495413 0,000000 0,000026 0,00000

43-84 -34,4540 -0,7467 94,1672 125,0308 0,815

0,469092 0,528889 0,001472 0,069976 0,0000

85-112 282,870 -101,588 50,634 -199,380 0,206

0,491273 0,631770 0,660733 0,567275 0,78420

113-138 170,859 64,840 -27,324 -137,853 0,842

0,285989 0,529289 0,180338 0,357782 0,00000

Рис. 3. Поля корреляции

оценку влияния объема закачки воды на дебиты нефти путем вычисления разностного параметра йО по следующей формуле:

йО=ом - ом)/ оТ.

Изменение значений йО для изучаемых скважин во всех случаях индивидуально.

По данным расчетов, выполненных по обобщенным и временным моделям для нагнетательной скв. 1, выполнено сопоставление От и О^ по 5 соседним добывающим скважинам (рис. 4, а).

Корреляционные связи между значениями Отм, вычисленными по обобщенным моделям, и фактическими дебитами От (г = 0,69, р = = 0,000) значительно слабее, чем между значениям Омв, вычисленными по дифференцированным моделям, и От (г = 0,99, р = 0,000). Уравнение регрессии между От и Ом, вычисленными по обобщенным моделям, имеет следующий вид:

От = 9,932 + 0,860 Отм .

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

От = -0,001 + 0,999 Отмв.

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

Для сравнения средних значений и От по скважинам будем использовать статистику t. Результаты расчетов (на примере нагнетательной скважины 1) приведены в табл. 3.

Анализ данных табл. 3 показывает, что по средним значениям наблюдаются статистические неразличия значений по От и О™ для скважин 1д, 2д, 3д, 4д. Для скважины 5д средние значения От и статистически различны. Аналогичные расчеты выполнены и по нагнетательным скважинам 2 и 3.

m cd X

X *

cd И

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

и

0

1

(N cd К К

S

и и

0

1

vo

cd К S

I

и и

0

1

S

CD

те ч о с

Л X X о

S ^

те

ч

CD

Л Л О

о S

рц

Таблица 3

Сопоставление средних значений От и О

Соседние добывающие бт, т/сут Отмо, т/сут QM , т/сут tмо tмв

скважины (общая модель) (временные модели) Р Р

0,25076 0,05076

1д 28,3 ± 22,6 27,6 ± 18,9 28,7 ± 22,6

0,802186 0,902186

0,21526 0,211665

2д 18,3 ± 19,4 17,8 ± 18,3 17,9 ± 18,5

0,829729 0,833527

3д 2,4 ± 3,2 2,3 ± 1,1 2,3 ± 2,7 0,081919 0,066063

0,9347771 0,947376

0,295935 0,285291

4д 19,3 ± 26,5 18,4 ± 23,3 18,4 ± 25,3

0,767506 0,775638

17,12945 0,011606

5д 56,2 ± 11,4 18,4 ± 23,3 56,2 ± 10,9

0,000000 0,990748

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

Данные на рис. 4, б свидетельствуют о том, что корреляция в первом случае (г = 0,72, р = 0,000) несколько слабее, чем во втором (г = 0,81, р = 0,000).

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

Корреляция в первом варианте (г = 0,69, р = 0,000) несколько слабее, чем во втором (г = 0,78, р = 0,000). Сравнение корреляционных полей, приведенных на рис. 4, показывает,

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

Сравнение средних значений изучаемых показателей в различных зонах коллекторов по критерию t приведено в табл. 4.

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

Таблица 4

Сравнение средних значений показателей по скважинам

Показатель Скважина t Р

1 2 3 Скважина 1-2 Скважина 1-3 Скважина 2-3

ГVh2O - Pr) 0,91 ± 0,02 0,91 ± 0,02 0,91 ± 0,02 0,307 0,7586 0,353 0,7237 0,026 0,9791

r(Pr - QT) -0,26 ± 0,24 -0,09 ± 0,17 -0,29 ± 0,30 -6,171 0,0000 1,014 0,3109 6,411 0,0000

r(Prn - бт) -0,70 ± 0,38 -0,67 ± 0,40 -0,57 ± 0,41 -0,705 0,4809 -2,918 0,0037 -1,986 0,0479

QM, т/сут 28,3 ± 19,4 14,4 ± 21,3 5,4 ± 6,8 5,043 0,00000 6,140 0,00000 1,005 0,3156

Q™ , т/сут 28,1 ± 21,5 16,3 ±18,9 5,6 ± 23,1 4,576 0,00000 5,576 0,00000 0,877 0,3810

dQ -0,79 ± 1,45 -0,39 ± 0,71 -2,197 ± 0,81 -2,6828 0,0077 -4,193 0,00000 5,124 0,0000

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

На заключительном этапе изучения специфики влияния значений г(УН 0 - Рг), г(Рг - (т), г(Ргп - (), ОМ, бтмв , ¿б, определенных для различных зон коллекторов, на дебиты нефти используем метод дискриминантного анализа, основанный на построении линейных дискри-минантных функций (ЛДФ) [14, 15]. При построении ЛДФ трещинные карбонатные коллекторы будут иметь значения выше перечисленных показателей, вычисленных с помощью построенных моделей по ближайшим добывающим скважинам. В данном случае выполняется оценка взаимодействия между ними и нагнетательной скважиной 1. Поровые коллекторы известняков будут представлять модельные данные, которые получены по ближайшим добывающим скважинам. Здесь оценивается влияние на них нагнетательной скважины 2. Поровые коллекторы песчаников будут представлять модельные данные, полученные по 4 добывающим скважинам, расположенным вблизи к нагнетательной скважине 3.

По этим данным были построены следующие ЛДФ:

71 = 12,1046 г{УЩ0 - Рг) - 0,03597г(Рг - (т) -- 0,23894г(Ргп - бт) - 0,00209 (тм + + 0,01268 бтмв - 0,06982^6 - 8,55129 при Я = 0,970, х2 = 6726,37, р = 0,000000;

72 = 0,856605 г(УН20 - Рг) + 2,406087г(Рг - (т) + + 1,270443г(Ргп - бт) + 0,029147 (тм -- 0,069647 бтмв + 0,29827^6 + 0,686364

при Я = 0,718, х2 = 1366,184, р = 0,000000.

По данным функциям были вычислены значения 71 и 72, которые для различных типов коллекторов приведены на рис. 5.

Из рис. 5 видно, что значения 71 и 72 достаточно хорошо разделяются в пределах изучаемых групп коллекторов. Среднее значение 71 для трещинных коллекторов равно +3,034, для

поровых известняков +2,156, для поровых песчаников -6,338. Среднее значение 72 для трещинных коллекторов равно -1,318, для поровых известняков -1,133, для поровых песчаников -0,148. Достоверность распознавания коллекторов - известняков трещинного типа составляет 83,5 %, известняков порового типа - 96,4 % и песчаников порового типа -100,0 %. Все это показывает, что влияние нагнетательных скважин 1-3 на дебиты добывающих скважин характеризуется различными значениями избранных для анализа критериев. Таким образом, можно констатировать, что процессы нагнетания воды и вытеснения нефти в пластах для изучаемых типов коллекторов формируют дебиты нефти по различным сценариям. Таким образом, при анализе эффективности систем поддержания пластового давления в условиях различных типов коллекторов данное обстоятельство необходимо учитывать.

Рис. 5. Распределение значений 71 и 72 для различных по видам пустотности коллекторов

Заключение

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

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

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

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

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

1. Митрофанов В.П., Злобин А.А. Остаточная нефтенасыщенность и особенности порового пространства карбонатных пород / Филиал ООО «ЛУКОЙЛ-Инжиниринг» «ПермНИПИнефть» в г. Перми. - Пермь, 2003. - 240 с.

2. Денк С. О. Проблемы трещиноватых продуктивных объектов. - Пермь: Электронные издательские системы, 2004. - 334 с.

3. Черепанов С. С. Комплексное изучение трещиноватости карбонатных залежей методом Уоррена-Рута с использованием данных сейсмо-фациального анализа (на примере турнефамен-ской залежи Озерного месторождения) // Вестник Пермского национального исследовательского политехнического университета. Геология. Нефтегазовое и горное дело. - 2015. - № 14. -С. 6-12. Б01: 10.15593/2224-9923/2015.14.1.

4. Галкин В. И., Пономарева И. Н., Черепанов С. С. Разработка методики оценки возможностей выделения типов коллекторов по данным кривых восстановления давления по геолого-промысловым характеристикам пласта // Вестник Пермского национального исследовательского политехнического университета. Геология. Нефтегазовое и горное дело. - 2015. - № 17. -С. 32-40. Б01: 10.15593/2224-9923/2015.17.4.

5. Черепанов С.С., Мартюшев Д.А., Пономарева И. Н. Оценка фильтрационно-емкостных свойств трещиноватых карбонатных коллекторов месторождений Предуральского краевого прогиба // Нефтяное хозяйство. - 2013. - № 3. -С. 62-65.

6. Тиаб Дж., Доналдсон Э.Ч. Петрофизика: теория и практика изучения коллекторских свойств горных пород и движения пластовых флюидов: пер. с англ. - М.: Премиум Инжиниринг, 2009. - 868 с.

7. Houze O., Viturat D., Fjaere O.S. Dinamic data analysis. - Paris: Kappa Engineering, 2008. -694 p.

8. Warren J.E., Root P.J. The behavior of naturally fractured reservoirs // Soc. Petrol. Eng. J. - 1963. -Vol. 3, is. 3. - Р. 245-255. DOI: 10.2118/426-PA.

9. Tiab D. Modern core analysis. Vol. 1. Theory, core laboratories. - Houston, Texas, 1993. - 200 p.

10. Van Golf-Racht T.D. Fundamentals of fractured reservoir engeneering / Elsevier scientific publishing company. - Amsterdam - Oxford -New York, 1982. - 709 p.

11. Horne R.N. Modern well test analysis: A computer Aided Approach. - 2nd ed. - Palo Alto: Petroway Inc, 2006. - 257 p.

12. Определение параметров трещиноватости пород на основе комплексного анализа данных изучения керна, гидродинамических и геофизических исследований скважин / С.С. Черепанов, И.Н. Пономарева, А.А. Ерофеев, С.В. Галкин // Нефтяное хозяйство. - 2014. -№ 2. - С. 94-96.

13. Johnson N.L., Leone F.C. Statistics and experimental design. - New York - London -Sydney - Toronto, 1977. - 606 p.

14. Путилов И. С. Разработка технологий комплексного изучения геологического строения и размещения месторождений нефти и газа. -Пермь: Изд-во Перм. нац. исслед. политехн. ун-та, 2014. - 285 с.

15. Путилов И.С., Галкин В.И. Применение вероятностного статистического анализа для изучения фациальной зональности турне-фаменского карбонатного комплекса Сибирского месторождения // Нефтяное хозяйство. - 2007. -№ 9. - С. 112-114.

References

1. Mitrofanov V.P., Zlobin A.A. Ostatochnaia neftenasysh-chennost' i osobennosti porovogo prostranstva karbonatnykh porod [Residual

oil saturation and features of pore space of carbonate rocks]. Perm': PermNIPIneft branch of LUKOIL-Engineering LLC in Perm, 2003. 240 p.

2. Denk S.O. Problemy treshchinovatykh produktivnykh obektov [Challenges of naturally fractured reservoirs]. Perm': Elektronnye izdatel'skie sistemy, 2004. 334 p.

3. Cherepanov S.S. Kompleksnoe izuchenie treshchinovatosti karbonatnykh zalezhei metodom Uorrena-Ruta s ispol'zovaniem dannykh seismofat-sial'nogo analiza (na primere turne-famenskoi zalezhi Ozernogo mestorozhdeniia) [Integrated research of carbonate reservoir fracturing by Warren-Root method using seismic facies analysis (evidence from Tournaisian-Famennian deposit of Ozernoe field)]. Bulletin of PNRPU. Geology. Oil & Gas Engineering & Mining, 2015, no.14, pp.6-12. DOI: 10.15593/2224-9923/2015.14.1.

4. Galkin V.I., Ponomareva I.N., Cherepanov S.S. Razrabotka metodiki otsenki vozmozhnostei vydeleniia tipov kollektorov po dannym krivykh vosstanovleniia davleniia po geologo-promys-lovym kharakteristikam plasta [Development of the methodology for evaluation of possibilities to determine reservoir types based on pressure buildup curves, geological and reservoir properties of the formation (case study of Famennian deposits of Ozernoe field)]. Bulletin of PNRPU. Geology. Oil & Gas Engineering & Mining, 2015, no.17, pp. 32-40. DOI: 10.15593/2224-9923/2015.17.4.

5. Cherepanov S.S., Martiushev D.A., Pono-mareva I.N. Otsenka fil'tratsionno-emkostnykh svoistv treshchinovatykh karbonatnykh kollektorov mestorozhdenii Predural'skogo kraevogo progiba [Evaluation of the petrophysical and reservoir properties of carbonate reservoirs of the Upper Kama region on full-size core samples]. Neftyanoe Khozyaistvo - Oil Industry, 2013, no.3, pp.62-65.

6. Tiab Dzh., Donaldson E.Ch. Petrofizika: teoriia i praktika izucheniia kollektorskikh svoistv gornykh porod i dvizheniia plastovykh fliuidov [Petrophysics: theory and practice of study of reservoir properties of rocks and reservoir fluid movement]. Perevod s angliiskogo. Moscow: Premium Inzhiniring, 2009. 868 p.

7. Houze O., Viturat D., Fjaere O.S. Dinamic data analysis. Paris: Kappa Engineering, 2008. 694 p.

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

8. Warren J.E., Root P.J. The behavior of naturally fractured reservoirs. Soc. Petrol. Eng. J., 1963, vol.3, is.3, pp.245-255. DOI: 10.2118/426-PA.

9. Tiab D. Modern core analysis. Vol. 1. Theory, core laboratories. Houston, Texas, 1993. 200 p.

10. Van Golf-Racht T.D. Fundamentals of fractured reservoir engeneering. Amsterdam-Oxford-New York: Elsevier scientific publishing company, 1982. 709 p.

11. Horne R.N. Modern well test analysis: A computer Aided Approach. 2nd ed. Palo Alto: Petroway Inc, 2006. 257 p.

12. Cherepanov S.S., Ponomareva I.N., Ero-feev A.A., Galkin S.V. Opredelenie parametrov treshchinovatosti porod na osnove kompleksnogo analiza dannykh izucheniia kerna, gidrodinami-cheskikh i geofizicheskikh issledovanii skvazhin [Determination of fractured rock parameters based on a comprehensive analysis of the data core studies, hydrodynamic and geophysical well tests]. Neftyanoe Khozyaistvo - Oil Industry, 2014, no.2, pp.94-96.

13. Johnson N.L., Leone F.C. Statistics and experimental design. New York-London-Sydney-Toronto, 1977. 606 p.

14. Putilov I.S. Razrabotka tekhnologii komp-leksnogo izucheniia geologicheskogo stroeniia i razmeshcheniia mestorozhdenii nefti i gaza [Development of technologies for a comprehensive study of of geological structure and location of oil and gas fields] Perm': Izdatel'stvovo Permskogo natsional'nogo issledovatel'skogo politekhnicheskogo universiteta, 2014. 285 p.

15. Putilov I.S., Galkin V.I. Primenenie veroiatnostnogo statisticheskogo analiza dlia izucheniia fatsial'noi zonal'nosti turne-famenskogo karbonatnogo kompleksa Sibirskogo mestorozhdeniia [Use of probabilistic statistical analysis for the study of facies zonation of Tournaisian-Famennian carbonate complex of Siberian field]. Neftyanoe Khozyaistvo - Oil Industry, 2007, no.9, pp.112-114.

Просьба ссылаться на эту статью в русскоязычных источниках следующим образом:

Галкин В.И., Пономарева И.Н., Репина В.А. Исследование процесса нефтеизвлечения в коллекторах различного типа пустотности с использованием многомерного статистического анализа // Вестник Пермского национального исследовательского политехнического университета. Геология. Нефтегазовое и горное дело. - 2016. - Т.15, №19. - С.145-154. DOI: 10.15593/2224-9923/2016.19.5

Please cite this article in English as:

Galkin V.I., Ponomareva I.N., Repina V.A. Study of oil recovery from reservoirs of different void types with use of multidimensional

statistical analysis. Bulletin of PNRPU. Geology. Oil & Gas Engineering & Mining, 2016, vol.15, no. 19, pp.145-154. DOI: 10.15593/2224-9923/2016.19.5

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