Влияние изоморфизма на прочностные свойства карбонатных пород
Акелян Н.С., Вайнгардт А.А., Пискунова А.А., Чурков А.В.


ООО «НК «Роснефть»-НТЦ»

(ОГ ПАО «НК «Роснефть»)

В данной работе изучено возможное влияние изоморфных форм карбонатных минералов на деформационно-прочностные и фильтрационно-емкостные свойства горных пород. Проведен комплексный анализ литологических и геомеханических характеристик карбонатных коллекторов. Выявлены значительные отличия упруго-деформационных и фильтрационно-емкостных свойств (ФЕС) изоморфных форм доломита от обычного доломита. Рекомендовано выделять ожелезненный доломит в отдельный литотип и разделить разрез по механическим фациям с учетом этого выделения, что позволит оптимизировать геолого-технические мероприятия, проводимые на скважине.
Постепенная выработка высокопродуктивных терригенных пластов приводит к необходимости разработки более труднодоступных залежей, таких как карбонатные коллектора. Доля карбонатных пород, вторых по распространенности, составляет 18–20 % от общего количества осадочных пород с содержанием в них до 60 % добываемой в настоящее время нефти [1]. В то же время они значительно менее изучены по сравнению с терригенными породами в виду сложности и неоднородности их минерального состава, а также макро- и микроструктуры.
Карбонатные минералы, слагающие карбонатную породу, определяются анионом СО3-. Химический состав и структура карбонатов зависят в большей степени от катионов (таких, например, как Ca, Mg, Mn, Fe, Zn, Sr и др.). К основным породообразующим минералам относятся кальцит (и арагонит), доломит, а также более редко встречающиеся родохрозит, магнезит и сидерит [2]. Каждый из карбонатных минералов способен в довольно широком пределе изменять свой химический состав, что связано с проявляющимся изоморфизмом двухвалентных катионов. В связи с чем появляется большое многообразие минеральных соединений, обладающих определенными физико-химическими свойствами [3]. Поэтому для получения полного представления о строении и свойствах карбонатных залежей требуется более подробное изучение состава и содержания изоморфных примесей в карбонатных минералах.
Наряду с изучением вопроса о минеральном составе карбонатных коллекторов большое внимание уделяется изучению порового пространства, а также закономерностям формирования трещин в этих породах. Обширные исследования посвящены изучению связи трещиноватости пород с их физико-механическими свойствами [4–6]. Также известно, что упругие постоянные горных пород существенно зависят от многих факторов. Обобщая, можно объединить эти факторы в две группы: влияние минерального состава и внутренней структуры. Понятие «внутренняя структура» включает: морфологическую и композиционную структуры, степень связности — цементацию, сплошность, структуру порового пространства и микротрещин.
Настоящая работа посвящена изучению строения и свойств доломитов — одному из наиболее распространенных карбонатов после кальцита, а также влиянию его изоморфизма на прочностные свойства карбонатных горных пород. Объектом исследования послужили отложения турнейского яруса каменноугольной системы (С1t), приуроченные к валу Сорокина Варандей-Адзъвинской структурной зоны севера Тимано-Печерской нефтегазоносной провинции. Залежь относится к нефтенасыщенным карбонатным коллекторам порово-трещинного, реже — порового типов.
Литологическая характеристика отложений
На первом этапе было проведено подробное послойное литологическое описание породы, уточнены структурные особенности пород, выявлены их вторичные процессы, охарактеризовано пустотное пространство.
Отложения вскрыты скважиной в интервале глубин по бурению 2 359,25–2 859,56 м. Большая часть разреза представлена известняками с большими по мощности прослоями доломита и единичным прослоем аргиллита (рис. 1).
Рис. 1. Распределение литологических типов пород
в отложениях С1t
Интенсивная доломитизация привела к формированию в известковой толще различных по мощности прослоев вторичных мелко-тонкокристаллических доломитов. Контакты прослоев постепенные. Встречаются крупные гнезда белого разнокристаллического ангидрита и вкрапления пирита.
По разрезу прослеживается повышенная трещиноватость с многочисленными разнонаправленными трещинами, как открытыми, так частично и полностью выполненными глинисто-битуминозным материалом. Вниз по разрезу наблюдается значительное уменьшение следов окремнения и сульфатизации, повышается содержание глинистого материала. Текстуры пород сменяются неясно-волнисто-горизонтальнослоистыми, образованными миллиметровыми и более крупными прослоями глинистого материала, хорошо выражены следы биотурбации. Трещиноватость пород становится ниже, исчезают признаки нефтенасыщения.
Детальное изучение доломитов в шлифах на оптическом микроскопе Olympus BX51 показало, что структура пород мелко-тонкокристаллическая, иногда перекристаллизована до среднекристаллической (рис. 2).
Рис. 2. Доломит в шлифах на оптическом микроскопе, увеличение 40 (николи II, николи +)
В основной массе пород угадываются реликты перекристаллизованного органического детрита. Среди вторичных процессов отмечается интенсивная перекристаллизация, реже — окремнение и сульфатизация основной массы. Пустотное пространство пород представлено трещинами, порами и единичными кавернами. Поры редкие, преимущественно межкристаллические, размером до 0,05 мм, приурочены к участкам перекристаллизации, их количество не превышает 2 %. Трещины как открытые, так и частично и полностью выполненные глинистыми минералами и органическим веществом, реже — ангидритом, толщиной до 0,75 мм. Кроме того, вдоль крупных пустот наблюдаются процессы окремнения.
Минералогические исследования породы
На втором этапе были проведены детальные минералогические исследования порядка 280 образцов кернового материала, предварительно прошедших экстракцию в аппаратах Сокслета.
Для изучения карбонатных минералов используются различные методы анализа. Одним из относительно быстрых и простых методов является изучение минералов в шлифах с дополнительным окрашиванием. Визуальная диагностика карбонатов устанавливает текстурно-структурные особенности, а также количественные соотношения между чистыми кальцитом, доломитом и другими карбонатами. Однако метод беспомощен при наличии минералов изоморфного ряда, и, как следствие, невозможно отличить железистые доломиты от магнезита или анкерит от магнезированного сидерита по показателям преломления.
Одним из самых точных и достоверных методов диагностики карбонатных минералов и их изоморфных смесей является рентгеноструктурный анализ, который с наибольшей полнотой отражает индивидуальные особенности их строения [7]. Каждое чистое карбонатное соединение характеризуется строго определенными величинами параметров элементарной ячейки и присущим только этому минералу набором межплоскостных расстояний и их интенсивностей. Для изоморфных разностей характерны промежуточные значения параметров элементарной ячейки, положения и интенсивности рефлексов. Определение качественного и количественного минерального состава горных пород проводилось методом рентгеновской дифрактометрии на порошковых неориентированных препаратах на рентгеновском дифрактометре Ultima IV фирмы Rigaku. Стандартные параметры съемки: анод — Cu (медь), напряжение рентгеновской трубки — 40 кВ, ток — 30 мА, мощность — 1 200 Вт, углы съемки 2Θ от 2° до 90°, скорость съемки 2°/мин.
Были получены и изучены дифрактограммы исследуемых образцов горных пород с общей карбонатностью в пределах 80–100 %. Из них 201 образец представлен известняками, 73 образца — доломитами. При исследовании дифрактограмм доломитов были замечены отличия в межплоскостных расстояниях от эталонных значений (рис. 3).
Рис. 3. Результаты рентгеноструктурного анализа:
а — доломит ожелезненный,
б — доломит

Изменение положения главной диагностической линии в доломитах было подробно изучено в работе И.Г. Демчук [8, 10].
Данный эффект объясняется присутствием в доломитах изоморфных примесей железа с образованием ожелезненного доломита (Са[Mg,Fe](СО3)2). Отмечена прямая зависимость сдвига главной диагностической линии в железистых доломитах от содержания изоморфного железа, замещающего магний в кристаллической решетке, что обусловлено изменением параметров решетки, увеличением межплоскостного расстояния. С увеличением содержания массовой доли FeO главная линия смещается левее. Предельное значение замещения ионов магния на ионы железа по данным [9] составляет порядка 50 %, соотношение Fe:Mg = 1:1,1.
В исследованных нами образцах доля ожелезненного доломита от общего числа доломитов составила порядка 30 % (28 образцов). По каждому образцу прослеживается описанная выше тенденция смещения главной диагностической линии d104 с увеличением доли оксида железа. Зависимости межплоскостных расстояний кристаллической решетки доломитов от содержания железа по данным [10] и данным, полученным нами в ходе исследований дифрактограмм на пробах отложений С1t, представлены в таблице 1.
Табл. 1. Зависимость межплоскостных расстояний кристаллической решетки доломитов от содержания железа

Дифференциальным термическим анализом (ДТА) было подтверждено наличие изоморфных карбонатных минералов в образцах горных пород. Метод основан на различной температуре разложения карбонатных минералов. Исследования проводились с использованием прибора синхронного термического анализа TGA/DSC 3+ компании Mettler Toledo на измельченных образцах горных пород, которые нагревались в инертной атмосфере азота от 50 до 1 000 °С со скоростью 10 °С/мин в корундовых тиглях. Так, на термограммах чистого доломита наблюдается термический эффект при 650–800 °С [11]. Ожелезненный доломит диагностируется на термограммах по наличию характерного дополнительного узкого эндотермического пика в области 730–760 °С, что можно связать с первоначальным распадом сложного карбонатного минерала на компоненты (CaCO3, MgCO3 и FeCO3) и одновременной диссоциацией его железистой составляющей [3]. Смещение этого пика в более низкотемпературную область, а также его интенсивность прямо пропорциональны содержанию железа в составе [10]. Дальнейшее разложение протекает так же, как и в простом доломите, но при более высокой температуре (800–900 °С) и связано с диссоциацией кальцитовой составляющей. Сравнение термограмм карбонатных минералов представлено на рисунке 4.
Рис. 4. Сравнение термограмм карбонатных минералов: а — доломит ожелезненный, б — доломит

Для получения изображения поверхности образцов карбонатных горных пород при больших увеличениях, с высоким пространственным разрешением с целью изучения морфологических особенностей и структуры порового пространства использовался растровый электронный микроскоп (РЭМ) SERON AIS1800C с комбинированной системой рентгеновского энергодисперсионного микроанализа Quantax EDS «XFlash 610 Mini» фирмы Bruker (ускоряющее напряжение 20 кВ, ток зонда 1,08 нА, t = 10 сек.), который позволяет проводить спектральный анализ и полуколичественно определять содержание химических элементов. Подготовка образцов проводилась путем откалывания кусочка необходимого размера перпендикулярно напластованию породы, после чего проводилось очищение поверхности образца от смещенных частиц и напыление золотом для снятия поверхностных зарядов.
Было показано, что породы имеют кристаллическую структуру от мелко-тонкокристаллической до иногда среднекристаллической. Визуально доломит и ожелезненный доломит дифференцировать невозможно. Результаты микрозондового анализа подтвердили ожелезнение доломита (рис. 5).
Рис. 5. РЭМ-изображения карбонатных горных пород: а — известняка, б — доломита, в — Fe-доломита

Упруго-деформационные свойства
У всех исследуемых образцов предварительно были определены основные петрофизические характеристики, в том числе коэффициент открытой пористости методом жидкостенасыщения по Преображенскому в атмосферных условиях, объемная плотность расчетным методом, коэффициент открытой пористости в пластовых условиях исходя из сжимаемости образцов на установке УИК-ПС «Гло-Бел Нефтесервис» [12].
В таблице 2 приведены средние значения петрофизических показателей для образцов с карбонатностью в диапазоне 80–100 %.
Табл. 2. Основные петрофизические характеристики породы

Было установлено, что пористость ожелезненного доломита близка по значениям к известняку и лежит в широком интервале 0,01–0,17 долей ед., а не к доломиту, пористость которого имеет более низкое значение и лежит в узком диапазоне 0,01–0,05 долей ед.
Значения открытой пористости в атмосферных условиях и пластовых условиях имеют несущественную разницу между собой для каждого из литотипов, что говорит о низкой сжимаемости их порового пространства и, в целом, о высокой крепости пород.
Такое существенное отличие результатов определения открытой пористости для доломита и доломита ожелезненного, по мнению авторов, может объясняться вторичными процессами, а также структурными нарушениями при изменении объема элементарной ячейки кристаллической решетки. В ходе сложных вторичных изменений в карбонатных коллекторах происходит замещение ионов кальция в известняке (СаСО3) ионами магния (СаMg(СО3)2),
образуется значительное пустотное пространство за счет уменьшения объема элементарной ячейки. Дальнейшее замещение ионов магния железом (Са[Mg,Fe](СО3)2) приводит к незначительному увеличению объема элементарной ячейки и дестабилизации ионных связей за счет большего радиуса ионов железа по сравнению с ионами магния, что, в свою очередь, также способствует образованию микронапряжений на границах элементов породы, представленных различными изоморфными формами доломитов, и, как следствие, большей вероятности возникновения структурных дефектов, увеличению пустотного пространства и формированию большей сети трещин. В итоге изменение пустотного пространства напрямую отражается на значениях открытой пористости и объемной плотности.
Было выдвинуто предположение о влиянии изоморфных замещений в доломите на упруго-деформационные свойства горных пород. В рамках данной работы акустическим методом измерялись времена прохождения упругих волн, рассчитывались скорости упругих волн и их интервальное время на установке Autolab-1500, NER (США). По результатам определений рассчитывались динамические упругие константы, в частности модуль Юнга и коэффициент Пуассона.
В дополнение, для установления различий механических свойств представленных пород, проведены прочностные эксперименты. Определены пределы прочности при одноосном растяжении (бразильский метод, тест TS) и при одноосном сжатии (тест UCS). Прочностные тесты проводились в условиях насыщения породы водным солевым раствором, соответствующим минерализации пластовых вод, скорость нагружения подбиралась таким образом, чтобы происходило первоначальное уплотнение пород, а также для учета всех релаксационных процессов при деформации в образце, и составляла 1–5 МПа/мин.
Построенные зависимости модуля Юнга и коэффициента Пуассона от пористости выявили отличие упругих свойств ожелезненного доломита. Так, на приведенном ниже графике зависимости модуля Юнга от пористости видно, что ожелезненный доломит имеет более низкие значения модуля Юнга, чем у доломита. В то же время график зависимости коэффициента Пуассона от пористости не имеет ярко выраженной корреляции. Тем не менее, как и модуль Юнга, коэффициент Пуассона у ожелезненного доломита в среднем имеет более низкие значения (рис. 6, табл. 3).
Таким образом, можно говорить о возможности выделения из всего разреза карбонатных пород фации с ожелезненным доломитом при интерпретации акустического каротажа.
Рис. 6. Графики зависимостей динамических упругих констант от пористости

Табл. 3. Упруго-деформационные характеристики

Также была построена зависимость интервального времени пробега продольной волны (DТ) от пористости. Справочные и расчетные значения интервального времени при нулевой пористости совпадают у доломита (DТсправоч. = 142 мкс/м, DТрасчет. = 142,31 мкс/м) и у известняка (DТсправоч. = 155 мкс/м, DТрасчет. = 155,02 мкс/м). Однако у доломита ожелезненного значение интервального времени составило 149,92 мкс/м. Можно заметить, что значения DT ожелезненного доломита находятся между значениями этого параметра для доломита и известняка (рис. 7, табл. 3). Описанные зависимости хорошо коррелируют с зависимостями интервального времени от параметров пористости и плотности: с увеличением пористости и плотности породы интервальное время уменьшается.
Прочностные свойства, такие как предел прочности при одноосном растяжении (бразильский метод) и предел прочности при одноосном сжатии, также показали существенную дифференциацию исследуемых карбонатов. По этим параметрам доломит ожелезненный значительно отличается от доломита и по своим значениям также более близок к известняку. Таким образом, доломит ожелезненный возможно отнести в отдельную механическую фацию также и по данному признаку.
Рис. 7. Зависимость интервального времени от пористости

Геомеханическое 1D-моделирование устойчивости пород в стволе скважины
Как известно, упругие свойства пород напрямую связаны с прочностью пород, такими как пределы прочности на разрыв и на сжатие [13, 14]. Данные параметры учитываются при построении цифровой модели устойчивости пород. Моделирование карбонатных коллекторов представляет особую проблему, поскольку они имеют меньшую изученность по сравнению с терригенными породами, а также большой разброс физических свойств [15, 16]. Для иллюстрации существенного влияния упругих и прочностных свойств различных по составу карбонатных пород в рамках данной работы было проведено 1D-моделирование устойчивости ствола скважины с учетом разделения механических фаций доломита и доломита ожелезненного с учетом полученных петрофизических характеристик и зависимостей (рис. 8). Моделирование производилось с помощью симулятора РН-СИГМА.
Рис. 8.
1D-моделирование устойчивости ствола скважины

В модели использованы полученные в рамках данной работы зависимости основных параметров: динамические и статические упругие константы, плотности пород, пределы прочности, результаты построения паспортов прочности пород, а также зависимости между этими параметрами с разбивкой по механическим фациям (доломит, ожелезненный доломит). Также для построения модели механических свойств и расчета градиентов гидроразрыва и обрушения пород по всему разрезу скважины использованы данные плотностного, акустического широкополосного каротажей, а в интервалах их отсутствия — синтетические кривые, восстановленные из данных каротажей сопротивления. Использованы некоторые константы, характерные для рассматриваемых месторождений и в целом для региона (коэффициент пороупругости Био, региональные напряжения и их азимут), полученные и обобщенные ранее. Выявлено, что градиент давления гидроразрыва пласта для ожелезненного доломита ниже в 2–2,5 раза, а градиент обрушения стенок ствола скважины выше; уточнен и коридор безопасного давления.
Результаты моделирования, проведенного с учетом уточненных фаций, показывают, что интервалы с преимущественным содержанием ожелезненных доломитов менее устойчивы и в целом более склонны как к потенциальному обрушению стенок скважины, так и автоГРП (спонтанный гидроразрыв пласта при превышении допустимого давления бурового раствора) при строительстве скважины. С целью минимизации подобных рисков на основании построенной с учетом выделения фации ожелезненных доломитов модели выдаются рекомендации при бурении. В дополнение стоит указать на существенное влияние учета различных изоморфных модификаций доломитов на результаты расчета параметров ГРП.
Акелян Н.С., Вайнгардт А.А., Пискунова А.А., Чурков А.В.

ООО «НК «Роснефть»-НТЦ» (ОГ ПАО «НК «Роснефть»), Краснодар, Россия

akelyanns@ntc.rosneft.ru
Комплексный подход к изучению на основе литологических, минералогических, петрофизических, а также упруго-деформационных исследований образцов карбонатных горных пород.
горные породы, карбонатные минералы, изоморфизм, упруго-деформационные свойства
Акелян Н.С., Вайнгардт А.А., Пискунова А.А., Чурков А.В. Влияние изоморфизма карбонатных минералов на деформационно-прочностные свойства горных пород // Экспозиция Нефть Газ. 2024. № 6. C. 28–34. DOI: 10.24412/2076-6785-2024-6-28-34
29.08.2024
УДК 549.02, 550
DOI: 10.24412/2076-6785-2024-6-28-34

Рекомендуемые статьи
© Экспозиция Нефть Газ. Научно-технический журнал. Входит в перечень ВАК
+7 (495) 414-34-88