Способ определения коэффициента вытеснения нефти башкирских карбонатных отложений Соликамской депрессии
METHOD FOR DETERMINING THE COEFFICIENT OF OIL DISPLACEMENT OF THE BASHKIR CARBONATE DEPOSITS OF THE SOLIKAMSK DEPRESSION

Изобретение относится к области нефтяной геологии и разработки нефтяных месторождений и является петрофизической основой подсчета извлекаемых запасов нефти из пород, представленных башкирскими карбонатными отложениями Соликамской депрессии.
Для понимания существа вопроса, следует указать, что Соликамская депрессия представляет собой крупную наложенную отрицательную структуру, сформировавшуюся в раннепермское время (около 300 млн лет назад). В тектоническом отношении приурочена к зоне Предуральского краевого прогиба. Расположена в северной части Пермского края. Отложения башкирского яруса представлены двумя типами разрезов: карбонатным и терригенно-карбонатным. Карбонатный тип разреза имеет наиболее широкое распространение и в нем преобладают органогенные, органогенно-обломочные и оолитовые известняки. Абсолютные отметки кровли башкирского яруса в пределах Соликамской депрессии (Предуральский краевой прогиб) находятся в интервале от - 1400 м до - 1800 м [Лядова Н.А., Яковлев Ю.А., Распопов А.В. Геология и разработка нефтяных месторождений Пермского края. - М: ОАО ВНИИОЭНГ, 2010. - 335 с.].
Известны универсальные способы определения коэффициента вытеснения нефти, применимые для различных отложений (Методические рекомендации по определению коэффициента вытеснения нефти водой расчетным способом / Михневич В.Г., Тульбович Б.И. - Пермь, 1980. - 12 с. и Методические рекомендации по определению коэффициента вытеснения нефти водой расчетным способом / Михневич В.Г., Тульбович Б.И., Хижняк Г.П. - Пермь, 1988. - 12 с.), включающие отбор образцов, определение проницаемости продуктивной породы по газу, вязкости нефти, а также определение дополнительно количества смол силикагелевых и асфальтенов в нефти. При этом значение коэффициента вытеснения нефти Квт водой, согласно указанным известным рекомендациям, рассчитывают по соответствующей математической формуле, с учетом измеренных значений проницаемости продуктивной породы по газу, вязкости нефти, количества в нефти смол силикагелевых и асфальтенов.
Недостатками указанного известного способа являются техническая сложность, потребность в дорогостоящем специальном оборудовании для моделирования пластовых условий. Кроме того, эти известные способы не отличаются высокой точностью (относительная погрешность - 5,9-6,2%; абсолютная погрешность - 3,2-3,4%).
н
н
Также известен способ определения коэффициента вытеснения нефти водой в лабораторных условиях, описанный в автореферате к диссертации «Комплексное решение проблемы оценки коэффициента вытеснения нефти в различных геолого-технологических условиях» (автор Хижняк Г.П., 2012 (http://www.dissercat.com/content/kompleksnoe-reshenie-problemy-otsenki-koeffitsienta-vytesneniya-nefti-v-razlichnykh-geologo-), согласно которому производят отбор образцов, определяют проницаемость продуктивной породы по нефти (Кпрн), вязкость (μ) нефти, а значение коэффициента вытеснения нефти для образца устанавливают с использованием отношения указанных выше параметров Кпрн/μ.
Недостатком указанного известного способа является недостаточная точность результатов по установлению коэффициента вытеснения нефти, а именно: относительная погрешность - находится в диапазоне 7,5-12,5%; абсолютная погрешность - 4,9-8,4%.
Наиболее близким к предлагаемому изобретению является способ определения коэффициента вытеснения нефти водой в лабораторных условиях (ОСТ 39-195-86 Нефть. Метод определения коэффициента вытеснения нефти водой в лабораторных условиях), который можно использовать для любых видов горных пород, в том числе и для башкирских карбонатных отложений Соликамской депрессии. Согласно известному способу производят изготовление из керна стандартных цилиндрических образцов, экстрагирование содержащейся нефти из изготовленных образцов, высушивание их до постоянной массы, и измерение размеров цилиндрических образцов, после чего производят определение коэффициента пористости (Кп) образца, его коэффициента проницаемости по газу (Кпрг), далее выполняют водонасыщение под вакуумом проэкстрагированных образцов пластовой водой, создают остаточную водонасыщенность (Ков), соответствующую начальному максимальному нефтенасыщению, производят донасыщение нефтью и выдерживают несколько часов при заданных давлении и температуре. Вытеснение нефти из образца производится с постоянной заданной скоростью, посредством прессов или насосов. Коэффициент вытеснения нефти Квт в известном способе определяют по формуле:
Квт=Vн/Vн.нач,
где Квт - коэффициент вытеснения нефти водой, д.ед.;
3
Vн.нач - объем нефти, первоначально содержащейся в образце, определяемый по разности объемов пустот и остаточной воды, см;
3
Vн - объем вытесненной из образца нефти, см.
Недостатками указанного известного способа являются техническая сложность, потребность в дорогостоящем оборудовании, многократная сборка и разборка кернодержателя после каждого испытания, что усложняет способ.
Кроме того, точность определения Квт указанным известным способом составляет величину 5% и более, что может привести к ошибочным выводам при подсчете запасов нефти в пласте.
Таким образом, опыт определения коэффициента вытеснения нефти показал, что применяемые в настоящее время известные способы определения не обладают достаточной точностью и надежностью в плане обеспечения достоверности результатов.
Технический результат, достигаемый предлагаемым изобретением, заключается в повышении точности определения коэффициента вытеснения нефти водой для башкирских карбонатных отложений Соликамской депрессии, за счет увеличения при его определении числа информационных показателей и определенной взаимосвязи между ними.
Дополнительным техническим результатом является упрощение способа, за счет исключения длительных и трудоемких лабораторных операций по определению коэффициента вытеснения нефти, и за счет того, что при его реализации используют параметры образца, определяемые в массовом порядке при стандартных лабораторных исследованиях керна.
н
Указанный технический результат достигается предлагаемым способом определения коэффициента вытеснения нефти башкирских карбонатных отложений Соликамской депрессии, согласно которому из керна реальной породы башкирских карбонатных отложений Соликамской депрессии изготовляют стандартные цилиндрические образцы, экстрагируют их от нефти, высушивают до стабилизации массы, взвешивают в сухом состоянии, определяют коэффициент пористости Кп, коэффициент абсолютной проницаемости по газу Кпрг и вязкость нефти μ,, под вакуумом производят 100%-ное насыщение образца пластовой водой или ее моделью, насыщенный образец взвешивают в воздухе, методом капилляриметрии производится вытеснение воды из образцов до величины остаточной водонасыщенности и определяют коэффициент остаточной водонасыщенности Ков, и далее устанавливают коэффициент вытеснения нефти Квт, новым при этом является то, что после 100%-ного насыщения водой образца под вакуумом дополнительно производят взвешивание его в воде, затем определяют объемную плотность ρ горной породы, из которой состоит образец, методом дискриминантного анализа рассчитывают значение дискриминантной функции Z образца по формуле:
Z=-101,442+191,381⋅Кп-29,490⋅Ков-13,620⋅Кпрг+35,118⋅ρ,
где Кп - коэффициент пористости, д.ед.;
Ков - коэффициент остаточной водонасыщенности образца, д.ед.;
2
Кпрг - коэффициент абсолютной проницаемости по газу, мкм;
3
ρ - объемная плотность горной породы, г/см;
по указанной дискриминантной функции Z определяют класс горной породы, из которого выполнен образец, исходя из следующего:
- при Z>0 - горную породу относят к первому классу;
- при Z<0 - ко второму классу;
М1
далее для образца, отнесенного к первому классу, рассчитывают коэффициент вытеснения нефти Квт по формуле:
М1
н
н
Квт=А0+А1⋅Кпрг+А2⋅μ+А3⋅(Кпрг/μ)+А4⋅ρ+А5⋅Ков,
где А0, Al, А2, A3 и А4 - коэффициенты и равны: А0=1,1483,
А1=-5,6251, А2=0,1718, А3=16,1795, А4=-0,4404, А5=-0,1534.
М1
Таким образом, для образца, отнесенного к первому классу, рассчитывают коэффициент вытеснения нефти Квт по формуле:
М2
А для образца, отнесенного ко второму классу, рассчитывают коэффициент вытеснения нефти Квт по следующей формуле:
М2
н
Квт=В0+В1⋅Кпрг+В2⋅(Кпрг/μ),
где В0, B1, В2 - коэффициенты и равны: В0=0,5712, В1=0,1914, В2=0,2823.
М2
Таким образом для образца, отнесенного ко второму классу, рассчитывают коэффициент вытеснения нефти Квт по следующей формуле:
М2
н
Квт=0,5712+0,1914⋅Кпрг+0,2823⋅(Кпрг/μ),
М1
где Квт - коэффициент вытеснения нефти образца первого класса, д.ед;
М2
Квт - коэффициент вытеснения нефти образца второго класса, д.ед;
н
μ - вязкость нефти, мПа⋅с.
н
Значение коэффициента вытеснения нефти Квт для конкретного продуктивного пласта рассчитывают по средним для всех образцов обоих классов значениям коэффициентов пористости Кп, абсолютной проницаемости по газу Кпрг и известному значению вязкости нефти μ.
Поставленный технический результат достигается за счет следующего.
Для понимания существа вопроса следует пояснить, что башкирские карбонатные отложения Соликамской депрессии, расположенной в северной части Пермского края, являются одними из основных объектов нефтедобычи в этих регионах. Коэффициент вытеснения нефти, входя в основную формулу подсчета извлекаемых запасов нефти, является петрофизической основой их подсчета в породах, представленных башкирскими карбонатными отложениями Соликамской депрессии.
н
Благодаря расширению количества используемых в предлагаемом способе информационных показателей (помимо определения стандартных характеристик: коэффициента остаточной водонасыщенности Ков, коэффициента пористости Кп, коэффициента проницаемости образца по газу Кпрг и вязкости нефти μ, дополнительно определяется объемная плотность ρ горной породы), а также благодаря разнесению пород по классам, возникает возможность дифференцированно оценить коэффициент вытеснения нефти водой различных зон продуктивных пластов и достоверно оценить извлекаемые запасы в башкирских карбонатных отложениях Соликамской депрессии.
За счет того что экспериментальным путем были получены определенные зависимости для математического расчета Квт для горных пород первого и отдельно для пород второго класса, обеспечивается дополнительная информационная связь указанных выше параметров образца между собой, что делает предлагаемый способ достоверным и точным.
М
Для доказательства такого вывода, приводим иллюстрацию в виде двух рисунков. На Рис. 1 приведена зависимость модельного, т.е. математически рассчитанного по уравнению (1) значения Квт, от экспериментального значения Квт без разделения данных на классы, полученная при реализации предлагаемого способа с использованием параметра «объемная плотность ρ горной породы»:
М1
М2
На Рис. 2 - приведена полученная при реализации предлагаемого способа аналогичная зависимость с разделением исходных данных на классы. Модельные (расчетные) значения Квт (для горной породы первого класса) и Квт (для горной породы второго класса), рассчитаны по уравнениям (2) и (3):
Корреляционное поле фактических и рассчитанных по уравнениям (2) и (3) значений коэффициента вытеснения (рис. 2) показывает, что они хорошо контролируют друг друга, при средней абсолютной ошибке 0,013 д.ед. и средней относительной
Таким образом, предлагаемый способ позволяет с высокой достоверностью определить значение Квт по данным стандартных исследований керна. Статистически обоснованы регрессионные уравнения для оценки Квт в башкирских карбонатных отложениях месторождений Соликамской депрессии.
Благодаря исключению ряда лабораторных операций, влекущих за собой необходимость использования дорогостоящего оборудования и требующих значительных временных затрат, способ является простым в реализации.
При реализации предлагаемого способа выполняли следующие операции в нижеуказанной последовательности (для наглядности, операции способа совмещены с примером конкретного осуществления).
1. Из керна продуктивных карбонатных пластов, представленных башкирскими карбонатными отложениями Логовского и Мысьинского месторождений Соликамской депрессии Пермского края изготовили стандартные цилиндрические образцы диаметром 0,30 см: 9 шт. - для Логовского месторождения, 10 шт. - для Мысьинского месторождения.
1
2. Затем указанные образцы подвергали экстрагированию от нефти и высушивали до постоянной массы, определяя вес P (г) каждого сухого проэкстрагированного образца.
2
3. У каждого образца определяли коэффициент пористости (Кп, д.ед.) (ГОСТ 26450.1-85. Породы горные. Метод определения коэффициента открытой пористости методом жидкостенасыщения) и абсолютную проницаемость по газу (Кпрг, мкм) (ГОСТ 26450.2-85. Породы горные. Метод определения коэффициента абсолютной газопроницаемости при стационарной и нестационарной фильтрации), например, с помощью автоматического прибора для определения пористости и проницаемости «АР-608» (США, «Coretest Systems»).
2
3
3
4. Далее под вакуумом на Автоматической установке насыщения керна AST-600 (США, «Core Lab Instruments») производили 100%-ное насыщение всех образцов жидкостью: пластовой водой или ее моделью, с последующим взвешиванием каждого насыщенного образца как в насыщающей жидкости Р (г), так и в воздухе Р (г), и с использованием этих параметров определяли объемную плотность ρ (г/см) горной породы, из которой выполнен образец, по формуле:
ж
3
где ρ - плотность жидкости, г/см.
При этом в качестве пластовой воды использовали ее модель - 4-х нормальный раствор хлорида натрия NaCl:
3
- Логовское месторождение - плотность - 1,177 г/см;
3
- Мысьинское месторождение - плотность - 1,159 г/см.
4
5. Затем методом капилляриметрии производили вытеснение воды из образцов, например, на Капилляриметре групповом В32-32 (Россия, НТЦ «Амплитуда»), до значения остаточной водонасыщенности, соответствующей начальной нефтенасыщенности продуктивного пласта, и определяли вес образца с остаточной водой Р (г). С учетом полученных показателей рассчитывали коэффициент остаточной водонасыщенности Ков по формуле:
6. С помощью линейного дискриминантного анализа [Девис Дж. Статистика и анализ геологических данных. М.: Мир, 1977 - 353 с.] по средним значениям коэффициента остаточной водонасыщенности Ков, плотности ρ горной породы и коэффициента пористости Кп рассчитывали дискриминатную функцию Z по формуле:
По знаку числового значения указанной дискриминантной функции Z определяли класс [Девис Дж. Статистика и анализ геологических данных. М.: Мир, 1977 - 353 с.] модели пласта исходя из следующего:
при Z>0 - горную породу относят к первому классу;
при Z<0 - горную породу относят ко второму классу.
Расчет показал, что дискриминантная функция Z образца Логовского месторождения имеет положительное значение (Z>0), т.е. образец относится к первому классу, дискриминантная функция Z образца Мысьинского месторождения имеет отрицательное значение (Z<0) - т.е. образец относится ко второму классу.
В таблице 1 приведены характеристики образцов пласта башкирских отложений. Причем следует указать, что из образцов одного месторождения, приведенных в таблице 1, компонуется модель пласта этого месторождения, и вместо лабораторного, длительного и трудоемкого определения на ней Квт по прототипу, этот Квт рассчитывается по вышеуказанной математической формуле (2) или (3) по средним значениям для образцов модели параметров Кп, Кпрг и др.
В таблице 2 приведены данные о средних характеристиках моделей пласта башкирских карбонатных отложений Логовского и Мысьинского месторождений Соликамской депрессии Пермского Прикамья и результаты дискриминантного анализа.
Значение P(Z)_1=0.998 для модели Логовского месторождения подтверждает правильность отнесения данной модели к 1-му классу, а значение P(Z)_2=1 для модели Мысьинского месторождения подтверждает правильность отнесения данной модели ко 2-му классу.
М1
М2
Для обоих классов методом регрессионного анализа обоснованы многомерные регрессионные уравнения для расчета значений коэффициента вытеснения нефти для породы первого класса Квт и для породы второго класса Квт по формулам:
- для 1 класса:
- для 2 класса:
Рассчитанные таким образом значения коэффициента вытеснения нефти равны:
- 1 класс, модель пласта Логовского месторождения - Квт=0,531;
- 2 класс, модель пласта Мысьинского месторождения - Квт=0,640.
М1
М2
абс
отн
Сравнение определенных по предлагаемому способу значений Квт и Квт с определенными экспериментально по прототипу для моделей башкирских карбонатных отложений Соликамской депрессии значениями Квт свидетельствует о высокой точности предлагаемого способа. Абсолютная ошибка расчета (ΔКвт) не превысила 0,006 ед., относительное отклонение (ΔКвт) - не более 0,85%. Данные приведены в таблице 3.
М
М1
М2
Достоверность предлагаемого способа в сравнении с прототипом представлена на рис. 1 и рис. 2. На рис. 1 приведена зависимость фактических Квт и модельных Квт значений коэффициента вытеснения башкирских карбонатных отложений Соликамской депрессии, рассчитанных без разделения моделей пласта на классы. На рис. 2 приведена зависимость фактических Квт и модельных Квт, Квт значений коэффициента вытеснения после разделения исходной выборки на классы.
М1
М2
Данные, приведенные на рис. 1, показывают, что при сопоставлении модельных и фактических значений Квт поле корреляции состоит из двух частей: при значениях Квт меньших 0,53 - модельные значения преимущественно превышают фактические, при Квт больших 0,53 - модельные значения существенно отличаются от экспериментальных, и корреляционное поле имеет значительный разброс вокруг линии равных значений. Это позволило сделать предположение о неоднородности выборки и наличии обособленных групп значений. Данные, приведенные на рис. 2, показывают, что фактические Квт и модельные (по классам) Квт, Квт значения коэффициента вытеснения хорошо контролируют друг друга, при средней абсолютной ошибке 0,013 д.ед. и средней относительной .
Погрешность относительная по прототипу составляет 5%, в то время как по предлагаемому способу погрешность выявлена в пределах 0-0,85% (табл. 3).
Благодаря тому, что предлагаемый способ позволяет проводить определение Квт с точностью, превышающей в 5,9 раз точность определения этого показателя по прототипу, обеспечивается наиболее точное представление о распределении запасов нефти. Это позволяет не только уточнить подсчет запасов, но и внести коррективы при осуществлении проекта рациональной разработки нефтяной залежи.
Примечание:
P(Z)_1 - вероятность отнесения образца пласта к 1-му классу,
P(Z)_2 - вероятность отнесения образца пласта ко 2-му классу.
абс
ΔКвт - абсолютная ошибка расчета;
отн
ΔКвт - относительное отклонение;
Квт - коэффициент вытеснения нефти, определенный по прототипу, д.ед.;
М
Квт - коэффициент вытеснения нефти, определенный по предлагаемому способу, д.ед.