УДК 622.276.21
ВЛИЯНИЕ ФИЛЬТРАЦИОННО-ЕМКОСТНЫХ СВОЙСТВ ОБЪЕКТОВ РАЗРАБОТКИ НА ВЕЛИЧИНУ КОЭФФИЦИЕНТА ВЫТЕСНЕНИЯ НЕФТИ В РАЗЛИЧНЫХ ГЕОЛОГО-ФИЗИЧЕСКИХ УСЛОВИЯХ
Гладких Евгений Александрович1,
Хижняк Григорий Петрович1,
Галкин Владислав Игнатьевич1,
1 Пермский национальный исследовательский политехнический университет, Россия, 614990, г. Пермь, Комсомольский пр., 29.
Актуальность исследования обусловлена недостаточной изученностью мелких и мельчайших месторождений нефти, что создает определенные сложности при подсчете запасов, составлении документов на разработку и оценке ее рентабельности. Такие месторождения располагаются часто вблизи крупных разрабатываемых объектов, поэтому произвести некоторые расчеты, например, оценить коэффициент вытеснения нефти в условиях отсутствия керна, можно по аналогии с соседними месторождениями, используя аналитические зависимости, получение которых является актуальной задачей.
Цель исследования: изучить влияние фильтрационно-емкостных свойств на величину коэффициента вытеснения нефти в различных геолого-физических условиях.
Объекты: башкирские карбонатные и визейские терригенные продуктивные отложения месторождений Башкирского свода и Соликамской депрессии Пермского края.
Методы: метод пошагового регрессионного анализа данных лабораторных определений коэффициента вытеснения. Результаты. Показан способ аналитической обработки статистических данных лабораторных определений коэффициента вытеснения нефти методом пошагового регрессионного анализа. Установлено, что на коэффициент вытеснения в пределах структурного элемента влияют различные параметры пласта и свойства флюидов. На примере башкирских карбонатных и визейских терригенных отложений показано наличие значимых корреляционных связей коэффициента вытеснения с фильтрационно-ем-костными свойствами коллекторов и вязкостью нефти. Установлены диапазоны влияния этих параметров на коэффициент вытеснения. Полученные для изучаемых пластов регрессионные уравнения свидетельствуют о том, что процессы вытеснения в од-новозрастных отложениях Башкирского свода отличаются от таковых в Соликамской депрессии. Параметры, включенные в регрессионные уравнения, показывают, что на начальном этапе значения коэффициента вытеснения обеспечиваются начальной нефтенасыщенностью и пористостью коллектора. Дальнейший рост коэффициента, связанный, вероятно, с доотмывом подвижной нефти после прорыва фронта вытеснения, определяется фильтрационными характеристиками пласта и вязкостью нефти.
Ключевые слова:
Коэффициент вытеснения, керн, регрессионное уравнение, нефтенасыщенность, пористость, проницаемость, вязкость нефти, фильтрационно-емкостные свойства.
Введение
Эпоха открытий месторождений «гигантов» уже позади. Большая часть таких месторождений, в разное время обеспечивавших основную долю годовой добычи нефти в стране, сегодня находится на 3-й, 4-й стадиях разработки с падающей добычей. Удержание или даже увеличение уровня добычи нефти может быть связано с запуском в разработку значительного числа вновь открываемых мелких месторождений. Однако зачастую это сопряжено с рядом трудностей, среди которых недостаточная изученность залежей. Как правило, такие объекты расположены неподалеку от более крупных разрабатываемых и достаточно изученных месторождений. В этой ситуации геолого-физическая информация по ним проецируется на од-новозрастные отложения более мелких «соседей», что позволяет в значительной мере сократить время составления проектной документации и введения залежей в разработку.
Важным этапом на стадии проектирования разработки месторождений является обоснование
проектного значения коэффициента извлечения нефти (КИН). Для его расчета академиком А.П. Крыловым предложена формула [1]
КИН = К • К ,
вт охв >
где Квт - коэффициент вытеснения нефти; Кохв - коэффициент охвата пласта вытеснением. Эта формула позже была дополнена рядом коэффициентов [2], характеризующих технологии разработки залежи, но коэффициент вытеснения остался в неизменном виде, что характеризует степень значимости данного параметра.
По причине непостоянства коэффициента вытеснения нефти по площади залежи в неоднородных пластах [3] лабораторное определение Квт с использованием кернового материала и флюидов изучаемого объекта [4] желательно проводить для различных геолого-физических условий [5]. Тем не менее при освоении мелких и мельчайших месторождений оно обычно слабо освещено керновым материалом, либо он полностью отсутствует. Тогда коэффициент вытеснения может быть оценен по
аналогии с соседними разрабатываемыми месторождениями с использованием аналитических зависимостей, получение которых является актуальной задачей.
Пермский край является относительно старым районом нефтедобычи. По состоянию на 01.01.2017 на территории края открыто 225 нефтяных месторождений, 176 из которых разрабатываются. В лабораториях крупнейшего в регионе научно-исследовательского и проектного института Перм-НИПИнефть с 1970 г. выполнено более 1500 лабораторных определений Квт. Накоплен значительный объем статистических данных по разрабатываемым месторождениям нефти в крае. Авторами выполнены систематизация и анализ этих данных. Некоторые результаты этой работы изложены в данной статье.
Сведения об объектах исследования
В тектоническом отношении территория края довольно разнообразна и объединяет 19 крупных тектонических структур, которые подразделяются на структуры платформенного типа и приуроченные к зоне прогибов.
Наиболее ярким представителем платформенного типа является Башкирский свод (БС), расположенный на юге края. В его границах расположено 68 месторождений, что составляет 29 % месторождений нефти и газа Пермского края. На территории Соликамской депрессии (СД) - представителя второй группы структур, открыто 38 месторождений углеводородов (УВ) или 16 % разведанных месторождений рассматриваемой территории. Ввиду интенсивной разработки месторождений на выделенных структурных элементах - 54 (79 %) на БС и 28 (74 %) на СД, а также их перспективности к открытию более мелких месторождений УВ, данные структуры были выбраны для анализа.
Более половины разведанных запасов нефти края содержится в визейских песчаниках тульского и бобриковского горизонтов. Залежи визейской нефти открыты в 134 месторождениях, в том числе 46 месторождений на территории БС, 20 на СД.
Второе место по запасам принадлежит башкирским карбонатам, залежи нефти в которых присутствуют в 77 месторождениях - 22 на БС, 16 на СД.
Из-за особенностей геологического развития и осадконакопления одновозрастные отложения в зоне прогибов залегают на значительно больших глубинах, по сравнению со структурами платформенного типа. Так, абсолютные отметки кровли башкирского яруса на Соликамской депрессии изменяются от -1400 до -1800 м против -800...-1000 м на Башкирском своде [6]. Этим определяется большая упакованность пород на Соликамской депрессии, ее большая плотность и меньшие фильтра-ционно-емкостные свойства по сравнению с Башкирским сводом. В то же время нефти месторождений Соликамской депрессии легкие и маловязкие в отличие от БС. В работе исследовалось влияние фильтрационно-емкостных свойств (ФЕС) на фор-
мирование величины Квт в карбонатных и терри-генных породах рассмотренных тектонических структур.
Исходные данные для анализа
В работе [5] показано наличие тесных корреляционных связей Квт с пористостью и проницаемостью коллекторов, которые также связаны друг с другом [7, 8]. В данном исследовании количество изучаемых параметров было увеличено. В работе использованы результаты лабораторных определений коэффициента вытеснения нефти и соответствующие ФЕС моделей пласта. В рассматриваемых отложениях (табл. 1) выполнено 156 определений Квт. В исходные выборки были включены основные параметры, характеризующие ФЕС моделей, - пористость (Кп), газопроницаемость (Кпрг), остаточная водонасыщенность (Ков), объемная плотность горной породы (р), свойства флюида -вязкость нефти (д), комплексный параметр (Кп/р) и коэффициент подвижности (Кпрг/дн), а также соответствующие им экспериментально определенные значения Квт (табл. 1).
Таблица 1. Характеристики исходных выборок
Table 1. Characteristics of initial samples
Структурный элемент Structural element Продуктивные отложения Productive deposits Параметры выборки Sample parameters Объем выборки(N) Sample size (N)
Башкирский свод Bashkir arch башкирские карбонатные (Бш) Bashkirian carbonate deposits Кп, КПрг, Ков, р, Рн, Кп/р, Кпрг/^н, Квт ф, kan, Srw, р, ^ ф/р, k3lr/mj0, Ed 19 определений 19 cases
визейские терригенные (C1V) Visean terrigenous deposits 71 определение 71 cases
Соликамская депрессия Solikamsk depression башкирские карбонатные (Бш) Bashkirian carbonate deposits 20 определений 20 cases
визейские терригенные (C1V) Visean terrigenous deposits 46 определений 46 cases
Примечание: Квт - коэффициент вытеснения нефти; Кп - пористость; К„рг - проницаемость; Ков - остаточная водонасыщенность; р - объемная плотность породыI; ¡1„ - вязкость нефти.
Note: ED - displacement efficiency; ф - porosity; kair - air permeability; S„ - residual oil saturation; p - bulk density of rock; д - oil viscosity.
Обработка и анализ данных
Обработка данных выполнена, руководствуясь [9-15], с применением метода пошагового регрессионного анализа и изложена в [16-18]. Суть его заключалась в том, что исходная выборка ранжировалась по возрастанию значений Квт. Затем по значениям выборки строились регрессионные уравнения сначала по первым трем определениям, где зависимой переменной выступал Квт, а независи-
мыми - остальные параметры выборки. Далее добавлялось еще одно определение Квт, и уравнение регрессии перестраивалось уже по четырем определениям, и так далее. При получении уравнения на первом шаге в него включается параметр, оказывающий наибольшее влияние на зависимую переменную. На следующем шаге добавляется параметр с меньшей степенью влияния и т. д. Каждое уравнение строилось исходя из условий его наибольшей статистической значимости. Такой подход при обработке данных позволил в динамике проследить влияние параметров пласта и флюида на формирование коэффициента вытеснения нефти для пласта в целом.
Результаты и обсуждение
В процессе реализации пошагового регрессионного анализа для каждой из исходных выборок было получено N-2 линейных регрессионных уравнений, где N - объем выборки данных. Общий вид уравнений был следующим:
КМ = А+А • Кп + А • Кпрг + Аз • Ков +
+а р+а м + а •—+а •—.
Р Мн
где КМ - модельное значение коэффициента вытеснения; А0-7 - регрессионные коэффициенты, подбираемые по методу наименьших квадратов.
Порядок включения параметров в уравнение (слева направо) определяется степенью их участия в формировании значений Квт, т. е. первым (после свободного члена А0) указывается наиболее значимый параметр, затем второй по значимости и т. д.
Каждое уравнение характеризуется коэффициентом детерминации Я2.
Башкирский свод
Рассмотрим зависимость коэффициента детерминации Я2 от Квт для башкирских карбонатных и визейских терригенных отложений месторождений Башкирского свода (рис. 1).
В корреляционном поле можно выделить следующие диапазоны:
• при Квт<0,56 точность аппроксимации уменьшается для обоих типов пород;
• при 0,56<Квт<0,6 значения Я2возрастают;
• при 0,6<Квт<0,64 (0,62 для Бш) коэффициент детерминации уменьшается;
• при Квт>0,64 (только для С1т) коэффициент детерминации возрастает.
Средние значения параметров выборок Башкирского свода в выделенных диапазонах Квт представлены табл. 2.
Для выделенных диапазонов Квт были получены регрессионные уравнения:
• при Квт<0,56:
Бш:КМ = 0,6218 -0,4282К ;
С, :КМ = 0,5587 -0,1267К .
IV вт ' ' ов
при 0,56<Квт<0,6:
Бш:КМ = 0,6359 -0,3597К ;
вт " " п7
К
С, :КМ = 0,5390-0,4746—!Ч
IV вт 7 7 7
Р
(1) (2)
(3)
(4)
■ при 0,6<Квт<0,64 (0,62 для Бш): Бш: КМ = 0,6103 + 0,0011мн -0,0521 -Ков; (5)
С, : КМ = 0,7411 -0,4706К -0,3720К ; (6)
"IV вт ' ' ов ' п7 4 /
Таблица 2. Средние значения параметров выборки Table 2. Average values of sample parameters
Значения Квт Values Ed Пласт Layer Кп, д.ед. ф, u.f. Кпрг, мкм2 karr, mpm2 Ков, д.ед. Srw, u.f. p, г/см3 (g/cm3) p, мПа-с Po, mPa-s Кп/p, см3/г ф/p, cm3/g Кпрг/рн, мкм2/ мПа-с karr/p, pm2/ mPa-s Квт, д.ед. ED, u.f.
<0,56 Бш 0,154 0,1120 0,226 2,29 12,75 0,0676 0,01041 0,525
C,v 0,185 0,1404 0,164 2,15 11,95 0,0865 0,01302 0,538
0,56-0,60 Бш 0,153 0,1410 0,170 2,28 12,70 0,0672 0,01199 0,581
C,v 0,205 0,5245 0,099 2,10 13,64 0,0979 0,05498 0,585
0,60-0,64 Бш 0,184 0,2297 0,166 2,20 11,01 0,0840 0,02088 0,614
C,v 0,211 0,5702 0,084 2,09 11,95 0,1009 0,07245 0,623
0,64-0,72 C,v 0,219 1,0761 0,064 2,06 9,68 0,1068 0,17718 0,664
• при Квт>0,64:
Clv : КМ = 0,6437 + 0,0283Кпрг -0,0010 • p (7)
Анализ уравнений (1)-(7) показывает, что на коэффициент вытеснения нефти в выделенных диапазонах влияют различные параметры. В диапазоне средних значений Квт (уравнения (1) и (2)) в карбонатных и терригенных породах с относительно низкими ФЕС (табл. 2) его значение формируется под влиянием остаточной водонасыщенности Ков, т. е. начальной нефтенасыщенности (Кнн), где Кнн=1-Ков. Чуть большие значения Квт (уравнения (3) и (4)) формируются уже при участии пористости коллектора, а в терригенных породах еще и плотности. При дальнейшем увеличении Квт (уравнения (5) и (6)) его величина в башкирских карбонатных отложениях определяется главным образом вязкостью нефти и в меньшей степени остаточной водо-насыщенностью, в то время как в терригенных отложениях визейского комплекса на Квт влияют остаточная водонасыщенность и пористость. Наибольшие значения Квт визейских отложений формируются под влиянием фильтрационных характеристик коллектора и вязкости нефти (7).
Из вышесказанного можно сделать вывод, что Квт башкирских карбонатных и визейских терригенных отложений Башкирского свода до значения 0,64 обеспечивается, в первую очередь, емкостными свойствами, а для коллекторов с высокими фильтрационно-емкостными свойствами значения Квт>0,64 формируются под влиянием, преимущественно, фильтрационной составляющей - проницаемости.
Соликамская депрессия
Выборки данных по Соликамской депрессии были проанализированы аналогично Башкирскому своду. На рис. 2 приведены зависимости Я2 от Квт для башкирских карбонатных и визейских тер-ригенных продуктивных отложений месторождений Соликамской депрессии.
Корреляционное поле для башкирских карбонатных отложений визуально более однородно и описывается одним уравнением (8) с высоким коэффициентом детерминации Я2=0,84:
К11 =-0,1782 + 32,2373К -
вт ' ' п
-0,6420К - 60,2686 Кп. (8)
Р
А \ л 4 г* А - А
* * А Ж \ f ч у
2 А
0,45 0,5 0,55 0,6 0,65 0,7 0,75 0,8
Km, Л еЛ
Рис. 2. Корреляционное поле значений R2и Квтдля башкирских карбонатных и визейских терригенных отложениях месторождений Соликамской депрессии Пермского края
Fig. 2. Correlation field of the values of Rand EDin the Bashkirian carbonate and Visean terrigenous deposits of the Solikamsk depression of the Perm region
Поле значений визейских терригенных отложений состоит из двух областей, в которых Квт аппроксимируется уравнениями (9) и (10): • при Квт<0,62
КМ =-2,6703 - 0,3045К +1,3917 р +
вт ов
+51,2591^-21,4910К .
п
Р
при Квт>0,62
КМ = 0,9294 - 0,5229К -
(9)
Я2 = 0,98; КМ2 = 2,2325 - 0,0079 Мн +
+0,2005К - 3,4625^ - 0,5951р,
Р
Я2 = 0,94.
(11)
(12)
КМ1 = 0,5883 - 0,3541К +
вт ов
К
+0,0291К + 0,2358-
Мн
Я2 = 0,82;
(13)
КМ2 =-0,5486 - 0,7384К +
вт ов
+1,9497Кп + 0,3784 р,
Я2 = 0,85; (14)
К
К = 0,5501 - 0,1500К +1,7990-
вт ов
Я2 = 0,62;
Мн
(15)
К
КМ3-2 = 0,5214 + 2,0426——-0,6285К ,
вт ов
вт Мн ов
Я2 = 0,77.
Соликамская депрессия Бш:
КМ1 = 1,1483 - 5,6251К + 0,1718 м +
вт прг н
К
+16,1795——-0,4404р-0,1534К ,
Мн ов
(16)
-0,0138мн - 0,0745р. (10)
Выполненный анализ показал для отложений Соликамской депрессии, как и для Башкирского свода, наличие групп значений, в которых влияющие на Квт параметры различны. Однако полученные уравнения не могут быть использованы для оценки Квт новых мелких объектов разработки, поскольку значение Квт по ним неизвестно и выбрать то или иное уравнение не представляется возможным.
Решение этой проблемы описано в работе [16], где при помощи дискриминантного анализа [19, 20] тех же самых выборок башкирских карбонатных и визейских терригенных отложений Башкирского свода и Соликамской депрессии выделены классы, для которых получены следующие многомерные регрессионные уравнения (11)-(20) оценки Квт: • Башкирский свод
Бш:
КМ1 = 6,2282 - 0,5105К +
вт ов
+0,0122мн -2,2167р-4,3314Кп,
Я2 = 0,99;
К1"2 = 0,5712 + 0,1914К + 0,2823вт прг
К
Мн
Я2 = 0,40.
(17)
(18)
С^:
(19)
КМ1 =-0,5483 - 0,2514К + 0,4811р+
вт ов
К
+2,6959——- 1,1552К + 0,0159м ,
Мн прг н
Я2=0,79;
К
КМ2 = 0,7885 - 0,075^—прг - 0,7088К -
вт ов
вт Мн ов
-0,0203Мн + 0,0671Кпрг,
Я2=0,81. (20)
В приведенных уравнениях КМ1, КвМт2, КМ3 - модельные значения коэффициента вытеснения соответственно для 1, 2 и 3 классов; КМ3-1, КвМт3-2 - модельные значения коэффициента вытеснения соответственно для 1 и 2 подклассов 3 класса.
Стоит отметить, что уравнение (18), полученное для башкирских отложений Соликамской депрессии, имеет низкий коэффициент детерминации. Используемые в данной работе параметры не позволили получить более тесную корреляционную связь, поэтому для Квт в этих отложениях может быть рекомендована формула (8).
По уравнениям рассчитаны модельные значения коэффициента вытеснения, и их значения сопоставлены с экспериментальными (рис. 3).
Выполненное сопоставление показало высокую сходимость модельных и экспериментальных значений Квт, средние абсолютные (Др,) и относительные (Дсортн) погрешности не превысили 0,017 д.ед. и 2,77 % соответственно. Полученные уравнения могут применяться для оценки Квт мелких месторождений, открываемых в пределах рассматриваемых территорий, а также при контроле потоковых исследований.
Выводы
В рассмотренных продуктивных отложениях влияние фильтрационно-емкостных свойств на коэффициент вытеснения нефти не одинаково. При этом различие имеет место не только для разных
С
к
Рис. 3. Сопоставление модельных и экспериментальных значений коэффициента вытеснения Fig. 3. Comparison of the model and experimental values of the displacement efficiency
0,65 i
0,60 -
5
0,55 -
0,50 -
0,45
Acpa6c=0,005 д.ед.
pgr
9'o
О класс 1 • класс 2
0,45 0,50 0,55
Km, Д-еД-Бш БС
0,60 0,65
0,70 -1
0,65 -
~ 0.60 -о
я
й* 0.55 -
S«
^ 0,50 -0,45 -0,40 -
А%5с=0,0 13 д.ед. ДсРОТН=2,11%
0
О класс 1 • класс 2
0,75 -| 0,70 -0,65 -
А °-60 "
£ 0,55 -0,50 0,45
Дсрабс= 0,017 д.ед. Дсрош=2,77%
О'
* >4
А
м
Г рф
А
Окласс 1 • класс 2 Д класс 3
0,45 0,50 0,55 0,60 0,65 0,70 0,75 Кет, д.ед.
Civ БС
0,75 -0,70 -ri0,65 -
SE0.60 -т
0,55 -0,50 -0,45
/ ✓ ✓ ✓
Л^.т =0 014 л ел
Дер J отн = 2,23 %
4L
оо < ✓
п
f О ✓ ✓ >' о э класс 1 * класс 2
✓ ✓ ✓
0,40 0,45 0,50 0,55 0,60 0,65 0,70 Д.ед.
Бш СД
0,45 0,50 0,55 0,60 0,65 0,70 0,75 0,80 Квт: Д.ед. Civ СД
структур, но и для отложений одного возраста в пределах одной структуры.
Получены многомерные модели оценки коэффициента вытеснения нефти в карбонатных башкирских и визейских терригенных отложениях месторождений Башкирского свода и Соликамской депрессии Пермского края.
Установлено, что для продуктивных отложений Башкирского свода, относящегося к структу-
СПИСОК ЛИТЕРАТУРЫ
1. Крылов А.П. Состояние теоретических работ по проектированию разработки нефтяных месторождений и задачи по улучшению этих работ // Опыт разработки нефтяных месторождений и задачи по улучшению этих работ. - М.: Гостоптехиздат, 1957. - С. 116-139.
2. Закиров И.С., Корпусов В.И. Коррекция структуры формулы для расчета коэффициента извлечения нефти // Нефтяное хозяйство. - 2006. - № 1. - С. 66-68.
3. Оценка распространения коэффициента вытеснения нефти по площади в неоднородных пластах / Р.Г. Сарваретдинов, А.А. Махмутов, Р.З. Нургалиев, В.А. Ленский, М.А. Токарев //
рам платформенного типа, величина коэффициента вытеснения определяется преимущественно емкостным параметром - начальной нефтенасыщен-ностью, а также вязкостью нефти.
В продуктивных отложениях Соликамской депрессии Предуральского краевого прогиба величина нефтевытеснения зависит в большей степени от фильтрационных параметров - проницаемости и подвижности.
Геология, геофизика и разработка нефтяных и газовых месторождений. - 2016. - № 12. - С. 59-61.
4. ОСТ 39-195-86. Нефть. Метод определения коэффициента вытеснения нефти водой в лабораторных условиях. - М.: Мин-нефтепром, 1986. - 19 с.
5. Закиров С.Н., Муртазалиев А.Ш. К определению коэффициентов вытеснения для терригенных и карбонатных коллекторов // Геология, геофизика и разработка нефтяных и газовых месторождений. - 2009. - № 9. - С. 45-47.
6. Лядова Н.А., Яковлев Ю.А., Распопов А.В. Геология и разработка нефтяных месторождений Пермского края. - М.: ОАО «ВНИИОЭНГ», 2010. - 335 с.
7. Permeability dependency on stiff and compliant porosities: a model and some experimental examples / S.A. Shapiro, G.P. Khizhniak, V.V. Plotnikov, R. Niemann, P.Yu. Ilyushin, S.V. Galkin // Journal of Geophysics and Engineering. - 2015. -V. 12. - P. 376-385. DOI: 10.1088/1742-2132/12/3/376.
8. Jing-qiang L., Chao-mo Z., Zhansong Z. Combine the capillary pressure curve data with the porosity to improve the prediction precision of permeability of sandstone reservoir // Journal of Petroleum and Engineering. - 2016. - V. 139. - P. 43-48.
9. Davis J.C. Statistics and data analysis in geology. 3rd ed. - New York; Chichester; Brisbane; Toronto; Singapore: John Wiley & Sons, 2002. - 656 p.
10. Афифи А.А., Эйзен С.П. Статистический анализ. Подход с использованием ЭВМ. - М.: Мир, 1982. - 488 с.
11. Johnson N.L., Leone F.C. Statistics and experimental design. -New York; London; Sydney; Toronto: John Wiley & Sons, 1977. -606 p.
12. Montgomery D.C., Peck E.A. Introduction to linear regression analysis. - New York: John Wiley & Sons, 1982. - 504 p.
13. Montgomery D.C., Peck E.A. Introduction to Linear Regression Analysis. - New York: John Wiley & Sons, 1982. - 504 p.
14. Statistics for Petroleum Engineers and Geoscientists / J.L. Jensen, L.W. Lake, P.W.M. Corbett, D.J. Goggin / Ed. by L.W. Lake. - New Jersey, USA: Prentice Hall PTR, 1997. - 390 p.
15. Afifi A., May S., Clark V.A. Practical multivariate analysis. 5th ed. - Hoboken: Chapman and Hall/CRC, 2011. - 537 p.
16. Способ оценки коэффициента вытеснения нефти на основе стандартных исследований керна / Е.А. Гладких, Г.П. Хижняк, В.И. Галкин, Н.А. Попов // Вестник Пермского национального исследовательского политехнического университета. Геология. Нефтегазовое и горное дело. - 2017. - Т. 16. - № 3. -С. 225-237. DOI: 10.15593/2224-9923/2017.3.3.
17. Гладких Е.А., Хижняк Г.П., Галкин В.И. Оценка коэффициента вытеснения нефти на основе стандартных исследований керна месторождений Пермского края // Нефтяное хозяйство. - 2017. -№ 8. - С. 90-93. DOI: 10.24887/0028-2448-2017-8-90-93.
18. Гладких Е.А., Хижняк Г.П. Разработка модели для оценки коэффициента вытеснения нефти (на примере башкирских карбонатных отложений) // Нефтепромысловое дело. - 2017. -№5. - С. 9-14.
19. Девис Дж.С. Статистический анализ данных в геологии: в 2 кн. - М.: Недра, 1990. - Кн. 1. - 319 с.
20. Девис Дж.С. Статистический анализ данных в геологии: в 2 кн. - М.: Недра, 1990. - Кн. 2. - 426 с.
Поступила 06.04.2018 г.
Информация об авторах
Гладких Е.А., аспирант кафедры нефтегазовых технологий Пермского национального исследовательского политехнического университета.
Хижняк Г.П., доктор технических наук, заведующий кафедрой нефтегазовых технологий Пермского национального исследовательского политехнического университета.
Галкин В.И., доктор геолого-минералогических наук, заведующий кафедрой геологии нефти и газа Пермского национального исследовательского политехнического университета.
UDC 622.276.21
INFLUENCE OF FILTRATION-CAPACITIVE PROPERTIES OF DEVELOPMENT OBJECTS ON THE VALUE OF OIL DISPLACEMENT EFFICIENCY IN VARIOUS GEOLOGICAL AND PHYSICAL CONDITIONS
Evgeniy A. Gladkikh1,
Grigoriy P. Khizhnyak1,
Vladislav I. Galkin1,
1 Perm National Research Polytechnic University, 29, Komsomolskiy Avenue, Perm, 614990, Russia.
The relevance of the research is caused by the insufficient study of small and tiny oil fields, which creates certain difficulties in calculating reserves, drawing up documents for development and evaluation of its profitability. Such deposits are often located near large developed objects, therefore it is possible to make some calculations, for example, to estimate the oil displacement efficiency in the absence of core, by analogy with neighboring deposits, using analytical dependencies, the receipt of which is an actual task. The main aim of the research is to study oil displacement and to determine the characteristics of the formation influencing the value of the displacement efficiency.
Objects of the researche are Bashkirian carbonate and Visean terrigenous productive deposits of the Bashkir arch and Solikamsk depression of the Perm region.
The authors have applied the method of stepwise regression analysis of laboratory determinations of the displacement efficiency. Results. The paper shows a method of analytical processing of statistical data of laboratory determinations of the oil displacement efficiency by the method of stepwise regression analysis. It is established that the displacement efficiency for formation within structural element is influenced by various parameters of formation and the properties of the fluids. The presence of statistically significant correlations between the displacement efficiency and the reservoir properties and the viscosity of oil is shown using the example of the Bashkir carbonate and Visean terrigenous deposits. The ranges of influence of these parameters on the displacement efficiency are established. The regression equations obtained for the strata studied by the layer indicate that the processes of displacement in the same-age sediments of the Bashkir arch differ from those in the Solikamsk depression. The parameters included in the regression equations indicate that at the initial stage the average values of the displacement efficiency are ensured by the initial oil saturation and porosity. A further increase in the displacement efficiency, probably associated with pre-washing of the mobile oil after the break-through of the displacement front, is determined by the filtration characteristics of formation and oil viscosity
Key words:
Displacement efficiency, core, regression equation, initial oil saturation, porosity, permeability, oil viscosity, filtration-capacitive properties of the rock.
REFERENCES
1. Krylov A.P. Sostoyanie teoreticheskikh rabot po proektirovaniyu razrabotki neftyanykh mestorozhdeny i zadachi po uluchsheniyu ehtikh rabot [Theoretical works in development of oil fields and tasks in improving these works]. Opyt razrabotki neftyanykh. mes-torozhdeniy i zadachi po uluchsheniyu etikh rabot [Experience in development of oil fields and tasks to improve these works]. Moscow, Gostoptekhizdat Publ., 1957. pp. 116-139.
2. Zakirov I.S., Korpusov V.I. Correction of structure of the formula for calculation of oil-recovery ratio. Oil industry, 2006, no. 1, pp. 66-68. In Rus.
3. Sarvaretdinov R.G., Makhmutov A.A., Nurgaliev R.Z., Len-sky V.A., Tokarev M.A. Evaluation of the spread of oil displacement efficiency by area in heterogeneous beds. Geology, geophysics and development of oil and gas fields, 2016, no. 12, pp. 59-61. In Rus.
4. OST 39-195-86. Neft. Metod opredeleniya koeffitsienta vytesne-niya nefti vodoy v laboratornykh usloviyakh [Oil. Method for determining the coefficient of oil displacement by water in laboratory conditions]. Moscow, Minnefteprom Publ., 1986. 19 p.
5. Zakirov S.N., Murtazaliev A.S. Definition of replacement coefficients for terrigenous and carbonateous collectors. Geology, geophysics and development of oil and gas fields, 2009, no. 9, pp. 45-47. In Rus.
6. Lyadova N.A., Yakovlev Yu.A., Raspopov A.V. Geologiya i razra-botka neftyanykh. mestorozhdeniy Permskogo kraya [Geology and development of oil deposits in Perm Krai]. Moscow, OAO VNIIO-ENG Publ., 2010. 335 p.
7. Shapiro S.A., Khizhniak G.P., Plotnikov V.V., Niemann R., Ily-ushin P.Yu., Galkin S.V. Permeability dependency on stiff and compliant porosities: a model and some experimental examples. Journal of Geophysics and Engineering, 2015, vol. 12, pp. 376-385. DOI: 10.1088/1742-2132/12/3/376.
8. Jing-qiang L., Chao-mo Z., Zhansong Z. Combine the capillary pressure curve data with the porosity to improve the prediction precision of permeability of sandstone reservoir. Journal of Petroleum and Engineering, 2016, vol. 139, pp. 43-48.
9. Davis J.C. Statistics and data analysis in geology. 3rd ed. New York; Chichester; Brisbane; Toronto; Singapore, John Wiley & Sons, 2002. 656 p.
10. Afifi A.A., Eyzen S.P. Statisticheskiy analiz. Podkhod s ispolzo-vaniyem EVM [Statistical analysis. Approach to computer use]. Moscow, Mir Publ., 1982. 488 p.
11. Johnson N.L., Leone F.C. Statistics and experimental design. New York; London; Sydney; Toronto, John Wiley & Sons, 1977. 606 p.
12. Montgomery D.C., Peck E.A. Introduction to linear regression analysis. New York, John Wiley & Sons, 1982. 504 p.
13. Montgomery D.C., Peck E.A. Introduction to Linear Regression Analysis. New York, John Wiley & Sons, 1982. 504 p.
14. Jensen J.L., Lake L.W., Corbett P.W.M., Goggin D.J. Statistics for Petroleum Engineers and Geoscientists. New Jersey, Prentice Hall PTR, 1997. 390 p.
15. Afifi A., May S., Clark V.A. Practical multivariate analysis. 5th ed. Hoboken, Chapman and Hall/CRC, 2011. 537 p.
16. Gladkikh E.A., Khizhnyak G.P., Galkin V.I., Popov N.A. Technique for evaluation of oil displacement coefficient based on the conventional core study. Perm Journal of Petroleum and Mi-ningEngineering, 2017, vol. 16, no. 3, pp. 225-237. In Rus. DOI: 10.15593/2224-9923/2017.3.3.
17. Gladkikh E.A., Galkin V.I., Khizhnyak G.P. The method for estimating the oil displacement coefficient based on standard core analysis. Oil industry, 2017, no. 5, pp. 90-93. In Rus.
18. Gladkikh E.A., Khizhnyak G.P. Development of a model applied to estimate the oil displacement factor (on the example of Bashkir carbonate deposits). Oilfield engineering, 2017, no. 5, pp. 9-14. In Rus.
19. Davis J.C. Statisticheskiy analiz dannykh v geologii [Statistics and data analysis in geology]. Moscow, Nedra Publ., 1990. B. 1, 319 p.
20. Davis J.C. Statisticheskiy analiz dannykh v geologii [Statistics and data analysis in geology]. Moscow, Nedra Publ., 1990. B. 2, 426 p.
Received: 6 April 2018.
Information about the authors
Evgeniy A. Gladkikh, postgraduate student, head of the Laboratory, Perm National Research Polytechnic University.
Grigoriy P. Khizhnyak, Dr. Sc., head of the Department, Perm National Research Polytechnic University. Vladislav I. Galkin, Dr. Sc., head of the Department, Perm National Research Polytechnic University.