Russian Federation
Russian Federation
Russian Federation
Russian Federation
Syktyvkar, Syktyvkar, Russian Federation
UDK 552.54 Карбонатные породы
UDK 552.14 Диагенез. Образование осадочных пород
The search for hydrocarbon deposits in the sediments of the Lower Paleozoic complex of the Timan-Pechora oil and gas province, which are associated with the main prospects for increasing reserves, is a relevant topic for research. It is necessary to understand the formation processes of carbonate reservoirs within the complex under consideration to develop a strategy for predicting the structure of the complex reservoirs. Obtaining data on the parameters of the fluid regime of the formation of secondary dolomites in the Paleozoic carbonate reservoirs of the Khasyrey field (Gamburtsev swell) was the aim of the study. We studied primary fluid inclusions in dolomite and anhydrite cements filling vuggy voids in reservoir rocks. We showed that the inclusions contained up to 10 vol. % of the gas phase, which was represented by methane and water vapor. The liquid phase was characterized by a substantially chloride-calcium composition, the salinity is 13.5–23.5 wt. % NaCl-eq. These data might indicate the formation of cavities and cements of secondary dolomites from moderate-temperature solutions of high salinity, which were formed as a result of mixing sedimentation (primary marine) waters and brines of hydrothermal vents. The products of evaporite dissolution in the presence of hydrocarbons and sulfate reduction processes were the most likely source and supplier of the fluids, both the formation of leaching caverns in the rock and the cements performing them. Fluids could enter the upper floors of the section from hydrothermal sources in the basement zones and lower horizons of the sedimentary cover through a system of faults and cracks during tectonic activation.
homogenization, cryothermometry, fluid inclusions, secondary dolomite, carbonate reservoir
Введение
Породы-коллекторы продуктивных толщ S2–D1 северо-востока Тимано-Печорской нефтегазоносной провинции представлены практически исключительно вторичными доломитами, в формировании пустотного пространства которых принимают участие поры, каверны и трещины (Майдль, 1989; Душин и др., 2016; Жемчугова, Маслова, 2020; Ульныров, Майдль, 2023). Коллекторские параметры в доломитовых породах сильно варьируются, и зависят от прохождения процессов доломитизации и/или рекристаллизации относительно исходной кристаллической структуры и образования вторичной пористости. Палеозойские доломиты часто более пористые, чем связанные с ними известняки. Последовательность процессов и относительное время доломитизации важны для понимания происхождения пористости в любом конкретном теле доломитов.
Крупные доломитовые тела — обычное явление в горных породах, и большинство исследователей считают их результатом реакции замещения с участием поровых вод, привносом Mg и удалением Ca (Land, 1985; Hardie, 1987). В последние годы возрастает понимание того, что замещающая доломитизация требует не только благоприятных геохимических условий, но и эффективного механизма циркуляции потоков жидкости. И большинство современных моделей генезиса массивных доломитов по существу являются гидрологическими (Machel, 2004).
Проблема прогноза геометрии доломитовых тел, их связности и распределения в них коллекторских свойств решается с помощью двух- и трехмерных количественных моделей реактивного транспорта (RTM). Они объединяют геохимические условия и реакции с характером флюидных потоков и основаны на исследовании процессов доломитизации и связанных с ними диагенетических реакций карбонатов. Существующие модели можно дифференцировать по трем фундаментальным драйверам, которые управляют циркуляцией флюидов, необходимых для доломитизации: высотный (топографический) градиент напора метеорной и/или морской воды; градиент плотности жидкости из-за изменений солености и/или температуры; градиент давления седиментационного и/или тектонического уплотнения. Потоки жидкости (циркуляция) обычно являются продуктом одновременного воздействия нескольких различных движущих сил. К настоящему времени разработаны модели доломитизации для четырех различных гидрогеологических систем: 1) рефлюкса рассола, 2) зоны смешения, 3) геотермальной циркуляции и 4) гидротермальной циркуляции, контролируемой разломами (Machel, 2004; Whitaker, Smart, Jones, 2004; Benjakul R. et al., 2020).
Характер доломитизации пород определяет их коллекторские свойства, поэтому моделирование активно внедряется в практику разведки и разработки нефтяных и газовых месторождений. От качества «первичных» данных о параметрах флюидов зависят и выбор концепции моделирования эпигенетических процессов, и точность прогноза распространения «улучшенных» коллекторов. Одним из главных поставщиков таких данных являются включения первичных растворов в растущих кристаллах, исследование которых позволяет судить о составе и температуре растворов, участвующих в преобразованиях породы-коллектора.
Представленная работа является первым опытом микротермометрического анализа флюидных включений в доломитовых коллекторах нашего региона. Целью работы являлось установление состава и температурного режима минералообразующих сред, сохранившихся в доломитовом и ангидритовом цементах трещиновато-кавернозных доломитов Хасырейского месторождения вала Гамбурцева.
Методы исследований
Исследования литологического состава и пустотного пространства образцов пород проводились с помощью микроскопа Leica DM-2700 М, а также на электронном микроскопе TESCAN VEGA3 LMH, оснащенном энергодисперсионным спектрометром (ЭДС) «Oxford Instruments X-Max 50 мм2» (ЦКП «Геонаука», г. Сыктывкар). Анализ изотопов углерода и кислорода проведен в ЦКП «Геонаука» Института геологии ФИЦ Коми НЦ УрО РАН на масс-спектрометре DELTA V Advantage. Значения изотопных коэффициентов определялись по стандартам PDBNBS18 и NBS19 (TS-limestone) для углерода и SMOW — для кислорода. Погрешность определения обоих коэффициентов не превышала ±0.1 ‰.
Изучение флюидных включений проведено в двух образцах из керна скв. 35-Хасырейская, отобранных из нижнедевонской и верхнесилурийской толщ.
Анализировались включения ромбоэдрического доломитового и позднего ангидритового цементов. Для выяснения генезиса вторичной минерализации в двусторонне полированных пластинах проводилось изучение флюидных включений методами гомогенизации и криометрии с использованием термо-криостолика THMSG600 фирмы Linkam. Погрешность измерений ±0.2 °С при отрицательных температурах и ±0.5 °С при положительных. Соленость растворов во включениях измерялась по температуре плавления льда (Bodnar, Vityk, 1994). Солевой состав включений определялся по температуре эвтектики водно-солевой системы (Борисенко, 1977). Газовый состав индивидуальных включений изучался частично на высокоразрешающем рамановском спектрометре LabRam HR800 (Horiba Jobin Yvon) при комнатной температуре. Для регистрации спектров применялась решетка спектрометра 600 ш/мм, размер конфокального отверстия составлял 300 и 500 мкм, щель — 100 мкм, мощность возбуждающего излучения He-Ne лазера (длина волны 632.8 нм) — 20 мВт, Ar+ лазера — 120 мВт (514,5 нм). Изучение флюидных включений проводилось в двусторонне полированных пластинах толщиной 0.3–0.5 мм. Часть спектров получена на рамановском микроскопе Ramos M520 (SOL instruments, Беларусь). Параметры регистрации спектров: твердотельный лазер (выходная мощность 80 мВт, l = 532 нм), решетка спектрометра 600 ш/мм, размер конфокального отверстия — 100 мкм, объектив ×50 (микроскоп Nexcope NE910, Китай), время накопления сигнала — 10 секунд, количество измерений на одном участке спектрального диапазона — 1. Регистрация спектров осуществлялась при комнатной температуре.
Геологическое положение и литологическая характеристика
Хасырейское месторождение приурочено к одноименному локальному поднятию вала Гамбурцева, расположенному в центральной части Варандей-Адзьвинской структурно-тектонической зоны (ВАСЗ) (рис. 1). Ее современный структурный план сложился в результате разнонаправленных субвертикальных и субгоризонтальных движений, тесно связанных с заложением и геологической эволюцией Уральской и Пай-Хойско-Новоземельской складчатых систем (Малышев, 2002). Вал Гамбурцева представляет собой серию расположенных кулисами чрезвычайно узких складок с относительно плоскими сводами, ограниченными дизъюнктивными нарушениями с запада и востока (Сотникова, 2009).
В составе нижнедевонских отложений (D1) вала Гамбурцева выделяются три пачки (снизу вверх): доломитовая (D1dol) в овинпармском горизонте, глинисто-доломитовая, ангидрито-доломитовая в сотчемкыртинском горизонте. Покрышкой залежей являются плотные аргиллиты, мергели тиманского и саргаевского горизонта верхнего девона.
На Хасырейском месторождении в составе пачки D1dol выделяется ряд интервалов, породы которых различаются морфологией пустотного пространства и литолого-петрофизическими свойствами (рис. 2).
Развитая в толще сеть трещин вторичной пористости и каверн способствует не только связности пластов пород-коллекторов овинпармского горизонта, но и (исходя из гидродинамики залежи) связывает их с силурийскими отложениями. Наиболее продуктивными являются отложения двух литологических типов. Первый тип представляют среднезернистые доломиты с низким содержанием глинистого материала и развитой межзерновой пористостью с доломитизацией (рис. 2, a), часто приуроченной к зонам развития микротрещин (рис. 2, с).
Второй тип пород-коллекторов формируют доломиты замещения отложениями нижней сублиторали. Они имеют первично известковый состав с крупными остатками колоний и раковин «рассеянных» в иловой матрице (рис. 2, b, с), относимых к отложениям нижней сублиторали. Доломитизация ила проходила относительно равномерно и привела к его замещению плотным агрегатом тонко-мелкозернистого доломита. Раковины и крупные обломки табулят, строматопорат обычно не образуют скоплений с «каркасной структурой», а погружены в иловую массу. Они преобразовывались по иным сценариям и в более поздние сроки. Их остатки либо полностью, либо частично замещены доломитом и кремнием, либо выщелочены до образования каверн, немного заполненных цементами крупнокристаллического доломита (рис. 2, d-h), иногда с кристаллами ангидрита. Аналогичные цементы заполняют также зоны дробления пород и трещины (рис. 2, i). Кристаллы доломита часто имеют зональное, подчеркнутое обилием включений, или блочное строение. Встречаются кристаллы «седловидной» формы (рис. 2, f). В ряде случаев грани кристаллов демонстрируют отчетливые следы растворения (рис. 2, e, h).
Среди вторичных доломитов выделяются несколько морфоструктурных разновидностей, сохраняющих реликты или полностью замещающих первичные структуры известняков. Микроскопическое и электронно-микроскопическое исследования позволяют выделить две-три генерации доломитов, отличающихся размером и морфологией, а часто и составом примесей. Доломиты первой и отчасти второй генерации характеризуются несколько повышенным содержанием кальция. Кристаллы второй генерации, слагающие зоны и участки рекристаллизации, демонстрируют незначительную примесь железа при недостатке магния. Содержание атомов железа составляет 0.2–0.5 ат. %. Третья генерация, представленная крупными кристаллами доломитовых цементов, часто содержит уже от 1.4 % до 9.8 ат. % железа, обогащающих внешние зоны кристаллов (рис. 3, a). В меньшем количестве в цементах заполнения присутствуют кристаллы ангидрита, обычно более поздние, чем доломитовые (рис. 3, b, c).
Цементы лишь частично заполняют пустоты, оставляя некоторый объем пространства свободным. Достаточно часто встречаются микроскопления кристаллов галита, выполняющих микротрещины и микропоры в доломитах (рис. 3, e, f). Это позволяет связывать их генезис с глубинными растворами, подпитываемыми рассолами погруженных соленосных толщ ордовика. Выделения пирита иногда имеют необычную пластинчатую форму зерен, выполняющих узкие промежутки между кристаллами доломита. Иногда сохраняются зерна с равным содержанием атомов железа и серы, возможно, реликты-свидетели изменений окислительно-восстановительных обстановок и редукции сульфатов.
Цементы, частично заполняющие каверны и пустоты расширения трещин из «рыхлых» агрегатов крупных кристаллов доломита с ровными или искривленными «седловидными» поверхностями граней, большинство исследователей связывают с геотермальными или гидротермальными растворами. Любые растворы, поступающие из нижних горизонтов осадочного чехла, имеют повышенную относительно вмещающих пород температуру и могут рассматриваться как геотермальные (Whitaker, Xiao, 2010). Подтверждением этому считают явное пространственное совпадение процессов цементации, пустот выщелачивания и трещиноватости пород (Davies, Smith, 2006; Du et al., 2018; Warren, 2019). Для суждения о генезисе доломитовых цементов, отнесении их к продуктам геотермальных или гидротермальных флюидов нужны данные, позволяющие судить о температурном режиме их образования. С этой целью нами были проведен изотопный анализ доломитовых фаз разных генераций и исследованы флюидные включения, обнаруженные в крупно-грубозернистых кристаллах доломита и ангидрита из карбонатных пород-коллекторов в керне скв. 35-Хасырейская.
Изотопные отношения d13C и d18O доломитов и доломитовых цементов
Было проведено изучение состава изотопов в разных доломитовых фазах кернов скважин 35-Хасырейская и 25-Черпаюская. Проанализировано 7 образцов доломитовых пород гребенского и овинпармского горизонтов. В каждом образце опробовались: крупнокристаллический молочно-белый доломит выполнения пустот и трещин (1), среднезернистый светло-серый доломит рекристаллизации матрицы (2) и мелкозернистый серый доломит менее измененных участков матрицы (3) (табл. 1). Все образцы демонстрируют облегчение изотопного состава по обоим параметрам от более раннего мелкозернистого доломита (3) к более позднему крупнозернистому доломиту (1). Амплитуда облегчения в отдельных образцах составляет 0.6–1.4 ‰ для значений d13С и 2.5–3.8 ‰ для величины d18О.
При этом более легкими по изотопам являются пробы из образца тектонической брекчии с признаками сдвиговой деформации. В менее нарушенных образцах наиболее легкими по изотопам являются доломитовые цементы поздней генерации. В целом пробы разных генераций демонстрируют облегчение изотопов по мере укрупнения размера кристаллов — от ранних фаз к фазам более поздним. Аналогичные результаты получены для доломитов овинпармского горизонта ручья Дэршор (Майдль, Нечаев, 2023). В них крупнозернистый доломит прожилков и цементов был обеднен (d13C от 0.0–1 до 5 ‰, а d18O от 1–3 до 6–7 ‰) по сравнению с доломитом диффузной доломитизации матрикса. Наибольшее облегчение фиксировалось для трещиноватых пород. При этом, совокупность проб крупнозернистых доломитовых цементов демонстрирует значимую ковариацию параметров d13С и d18О, в выборке проб матричного доломита линейная связь между параметрами отсутствует (Майдль, Нечаев, 2023).
Как известно, величина d18O широко используется в качестве индикаторов диагенеза и катагенеза карбонатных пород, хотя из-за влияния как температуры флюида, так и изотопного состава воды, трактовка значений этого параметра неоднозначна. Низкие значения d18О могут быть следствием: а) пресноводного диагенеза и инфильтрационных вод, обедненных d18О; б) растворов повышенных температур в зоне катагенеза (Swart P.K., 2015). Следовательно, карбонаты, сформированные при высоких температурах, могут иметь пониженные значения d18O, сопоставимые с воздействием метеорного диагенеза, и их изотопный состав не дает однозначного ответа о природе флюидов.
Результаты исследования флюидных включений
Для определения условий минералообразования нами изучены первичные включения в доломите и ангидрите, цементирующих каверну по выщелоченному органогенному обломку (рис. 4, а, обр. 3/3) и доломитовому цементу выполнения внутренних пустот замещенного остатка колонии табулят (рис. 4, b, обр. 13/3). К первичным отнесены те включения, которые встречаются поодиночке или небольшими группами.
Для доломита (обр. 3/3) характерны двухфазовые первичные включения с газовой фазой до 10 об. %. Они имеют разнообразную, иногда с элементами огранки форму, их размер не превышает 20 мкм (рис. 4, a). Температура эвтектики в жидкой фазе изменяется в пределах –62…–64 °С, что характерно для водных растворов хлорида кальция с примесями. Плавление происходит при –19…–19.3 °С, что свидетельствует о том, что соленость раствора равна 21–22.2 мас. % NaСl-экв. Гомогенизация включений происходит при 102–120 °С (табл. 2, рис. 5).
Доломит (обр. 13/3) полупрозрачный и отличается очень мелкими, не более 10 мкм, включениями неправильной формы. Двухфазовые первичные включения содержат газовую фазу до 10 об. %. Температура эвтектики жидкой фазы равна –68 °С, на основании этого можно предположить присутствие солей хлорида кальция с примесями солей лития. Температура плавления, равная –21.5 °С, свидетельствует о том, что концентрация солей близка 23.5 мас. % NaСl-экв.
Флюидные включения в ангидрите (обр. 3/3) отличаются крупным, часто более 100 мкм, размером и сильно вытянутой формой (рис. 5, b). Газовая фаза присутствует в них в небольшом количестве, менее 5 об. %. Температура эвтектики (–48…–60 °С) характерна для хлоридов кальция с примесями. Плавление последней льдинки происходит при температурах –9…–18 °С, что свидетельствует о том, что соленость раствора равна 13.5–21.0 мас. % NaСl-экв. Температура гомогенизации первичных флюидных включений равна 97–165 °С. (табл. 2, рис. 6). При визуальном изучении отмечается большое количество вскрытых включений. Вскрытие, вероятно, произошло во время изготовления пластинок. Кроме того, довольно часто наблюдаются включения с разным соотношением газовой и жидкой фазы (рис. 5, с). Для термометрии такие включения непригодны, так как образовались, скорее всего, в результате гетерогенного захвата, либо при контакте пластинок ангидрита с водой и спиртом на стадии изготовления.
Для определения газового состава флюидных включений проведено рамановское исследование. В доломитах наблюдается очень большой фон (рис. 7, а). Виден пик воды, редко — метана. Пики других газов отсутствуют или перекрыты фоном. Возможно, некоторые газы присутствуют в количествах ниже предела определения. В ангидритах фон заметно ниже, но, как и в доломитах, диагностируется вода, реже с метаном (рис. 7, b).
Обсуждение результатов
По морфологическим и геохимическим особенностям крупнозернистый ромбоэдрический или «седловидный» доломит выполнения пустот не отличается от аналогичных доломитов, описанных во многих разрезах катагенетически измененных карбонатных пород, включающих залежи углеводородов, воды древних водоносных горизонтов и рудные сульфидные месторождения типа долины Миссисипи — (MVT). Характерными признаками для них являются: послойное обогащение Fe, Mn и Ca, изменение изотопных отношений d13С в диапазоне от слегка положительных до умеренно отрицательных значений и умеренно отрицательные значения d18О (Warren, 2019).
В исследованных доломитах флюидные включения гомогенизируются при температуре от 102 до 120 °C. Соленость жидкой фазы флюидных включений равна 21.0–23.5 мас. % NaCl-экв, что значительно превышает соленость морской воды. В ангидрите, который является более поздним по времени образования, интервал температур гомогенизации включений несколько шире и составляет 97–165 °C. Соленость жидкой фазы изменяется в пределах от 13.5 до 21.0 мас. % NaCl-экв. Температуры эвтектики жидкой фазы включений в доломите и ангидрите предполагают наличие хлоридов кальция, но более низкая температура во включениях доломита может быть следствием примеси солей лития.
Полученные данные свидетельствуют, что кристаллизация цементов происходила из рассолов относительно умеренных температур с признаками прямой метаморфизации (накопления кальция) морских талассогенных вод. Как известно, метаморфизация сопровождается увеличением с глубиной общей минерализации подземных рассолов и увеличением концентраций ряда микроэлементов, достигающих максимальных величин в глубоких метаморфизованных рассолах хлоркальциевого типа (Краснов и др. 2004). В свете современных представлений о формировании седиментационных рассолов считается, что образование в них избыточных концентраций редких щелочных элементов происходит в результате взаимодействия с вмещающими породами первичных хлормагниевых рассолов, вследствие чего происходит метаморфизация их в хлоркальциевые рассолы (Крайнов и др., 2004), а вмещающие известняки постепенно замещаются доломитами. Максимальные же концентрации редких щелочных элементов приурочены к наиболее метаморфизованным рассолам, обладающим максимальными содержаниями кальция.
В рассольных водах литий накапливается в значительных количествах. Наиболее обогащены им внутрисолевые и надсолевые рассолы галогенных формаций (Крайнов и др., 2004).
Учитывая, что состав вод комплекса вала Гамбурцева и сопредельных территорий гряды Чернышева хлоридно-натриевый с высоким содержанием кальция, рассолы, законсервированные во включениях, являются более метаморфизованными, чем воды вмещающих пород. Присутствие в породах галита свидетельствует о напорном (илизионном) режиме поступления рассолов из нижних горизонтов осадочного чехла и их генетической связи с рассолами соленосных толщ ордовика. Метановый состав газовой фазы включений может быть связан с процессами генерации и термальной деградации углеводородов.
Интересным является тот факт, что температура гомогенизации включений в более позднем ангидрите выше, чем в доломите (рис. 3, 5). Вероятно, это связано с активизацией гидротермального процесса на более поздней стадии.
Важным является также то, что крупнозернистый доломит выполняет каверны и пустоты трещин, то есть его кристаллизацию предваряет процесс избирательного растворения — выщелачивания. Это накладывает требование — циркулирующие рассолы изначально должны быть агрессивны к карбонатным породам, насыщаясь минералами по мере циркуляции.
Растворение эвапоритов в присутствии углеводородов и продуктов термального восстановления сульфатов считается одним из наиболее вероятных источников рассолов, образующих седловидный доломит. Эти рассолы (как и рассолы магматического происхождения) отличаются повышенной кислотностью и могут быть химически очень агрессивны для известняков (Warren, 2019; Machel, 2001). При достаточно большом объеме растворения пород может формироваться гидротермальный карст (Gao et al., 2020).
Заключение
В изученных нами верхнесилурийско-нижнедевонских разрезах скважин вала Гамбурцева морфологические и минералогические особенности вторичных преобразований пород-коллекторов дают основание полагать значительное участие в этом процессе термальных агрессивных растворов, поступавших по трещинным зонам из подстилающих отложений. Формирование зон трещинной проницаемости и «раскрытие» глубинных водно-углеводородных систем, вероятно, было приурочено к основным для провинции рубежам тектогенеза P/T — T/J. Флюиды, принимавшие участие в образовании жильной минерализации продуктивных коллекторов вала Гамбурцева, характеризуются невысокой температурой 97–165 °C и повышенной соленостью. В них преобладают соли кальция, что соответствует характеристикам геотермальных рассолов, более измененным, чем воды вмещающих пород, включающим, вероятно, продукты растворения соленосных толщ ордовика и продуктов термальной сульфатредукции. Присутствие в породах галита свидетельствует о напорном (илизионном) режиме поступления рассолов из нижних горизонтов осадочного чехла и их генетической связи с рассолами соленосных толщ ордовика.
1. Belonin M. D., Budanov G. F., Danilevsky S. A., Prishchepa O. M., Teplov E. L. Timan-Pechora province: geological structure, oil and gas content and development prospects. St. Petersburg: Nedra, 2004, 396 p. (in Russian)
2. Borisenko A. S. Study of the salt composition of solutions of gas-liquid inclusions in minerals using cryometry. Geology and Geophysics, 1977, No. 8, pp. 16–27. (in Russian)
3. Dushin A. S., Melnikov A. V., Fedorov A. I., Rykus M. V. New fluid dynamic model of the carbonate reservoir of the R. Trebs field based on the synthesis of geological and field data. Oil and Gas Business, 2016, V. 14, No. 2, pp. 13–23. (in Russian)
4. Zhemchugova V. A., Maslova E. E. Secondary dolomitization as a factor determining the reservoir properties of Lower Devonian deposits on the eastern side of the Khoreyver depression (Timan-Pechora oil and gas basin). Bulletin of Moscow University. Series 4: Geology, 2020, No. 3, pp. 47–56. (in Russian)
5. Krainov S. R., Ryzhenko B. N., Shvets V. M. Geochemistry of groundwater: Theoretical, applied and environmental aspects. Moscow: Nauka, 2004, 677 p. (in Russian)
6. Maydl T. V. Lithology and reservoirs of productive carbonate deposits of the Lower Devonian of the Gamburtsev swell. Geology and resources of fossil fuels of the European North of the USSR (Proceedings of the Institute of Geology Komi SC UB RAS). Syktyvkar, 1989, pp. 34–46.(in Russian)
7. Maydl T. V., Nechaev M. S. Paleokarst, thermal karst of the Lower Devonian carbonate reservoirs of the Timan-Pechora province: nature of manifestation and diagnostics. Proc. of the All-Russian Scientific Conference (with foreign participation). Fundamental problems in the study of volcanic-sedimentary, terrigenous and carbonate complexes (Litol 2023). Moscow: GEOS, 2023, pp. 104–107. (in Russian)
8. Malyshev N. A. Tectonics, evolution and oil and gas content of sedimentary basins of the European north of Russia. Yekaterinburg: UB RAS, 2002, 271 p. (in Russian)
9. Sotnikova A. G. Varandey-Adzva aulacogen: zones of oil and gas accumulation in carbonate Middle Ordovician-Lower Devonian deposits and priority directions of geological exploration for hydrocarbon raw materials (onshore, Pechora Sea shelf). Geology of Oil and Gas, 2009, No. 2, pp. 10–21. (in Russian)
10. Ulnyrov I. L., Maydl T. V. Lithology, sedimentation conditions and formation of reservoir rocks of the Upper Silurian carbonate deposits of the Gamburtsev swell (Timan-Pechora oil and gas province). Oil and Gas Geology. Theory and practice, 2023, V. 18, No. 3. http://www.ngtp.ru/rub/2023/26_2023.html (in Russian)
11. Benjakul R., Hollis C., Robertson H. A., Sonnenthal E. L., Whitaker F. F. Understanding controls on hydrothermal dolomitisation: insights from 3D reactive transport modelling of geothermal convection // Solid Earth. 2020. Vol. 11. No. 6. P. 2439–2461. DOI:https://doi.org/10.5194/se-11-2439-2020
12. Bodnar R. J., Vityk M. O. Interpretation of microterhrmometric data for H2O-NaCl fluid inclusions. Fluid inclusions in minerals: methods and applications. Pontignano: Siena. 1994. P. 117—130.
13. Davies G. R., Smith L. B. Structurally controlled hydrothermal dolomite reservoir facies: An overview // Bulletin American Association Petroleum Geologists. 2006. V. 90. P. 1641—1690.
14. Du Y., Fan T., Machel H. G, Gao Z. Genesis of Upper Cambrian-Lower Ordovician dolomites in the Tahe Oilfield, Tarim Basin, NW China: Several limitations from petrology, geochemistry, and fluid inclusions // Marine and Petroleum Geology. 2018. V. 91. P. 43—70. DOI:https://doi.org/10.1016/j.marpetgeo.2017.12.023
15. Gao B., Tian F., Pan R. F., Zheng W. H., Li R., Huang T. J., Liu Y. S. Hydrothermal dolomite paleokarst reservoir development in wolonghe gasfield, Sichuan Basin, revealed by seismic characterization. // Water. 2020. V. 12. №. 2. 579 p. DOI:https://doi.org/10.3390/w12020579
16. Hardie L. A. Dolomitization: A critical view on some current views // Journal of Sedimentary Petrology. 1987. Vol. 57. P. 166–183.
17. Land L. S. The origin of massive dolomite // Journal of geological education. 1985. Vol. 33. P. 112–125.
18. Machel H. G. Bacterial and thermochemical sulfate reduction in diagenetic settings-old and new insights // Sedimentary Geology. 2001. V. 140. P. 143–175. DOI:https://doi.org/10.1016/S0037-0738(00)00176-7
19. Machel H. G. Concepts and models of dolomitization: a critical reappraisal. Geological Society Special Publications, 2004, vol. 235, pp. 7–63.
20. Swart P. K. The geochemistry of carbonate diagenesis: The past, present and future // Sedimentology. 2015. Vol. 62. P. 1233–1304. DOI:https://doi.org/10.1111/sed.12205
21. Warren J. Saline dolomites: Ancient — part 4 of 4 — saddle dolomite, calcite and anhydrite spar burial salts. Salty matters. 2019. P. 1–22. URL: http://www.saltworkconsultants.com/blog-salty-matters
22. Whitaker F. F., Smart P. L., Jones G. D. Dolomitization: from conceptual to numerical models // Special Publication of the Geological Society. 2004. Vol. 235. P. 99–139. DOI:https://doi.org/10.1144/GSL.SP.2004.235.01.05
23. Whitaker F. F., Xiao Y. Reactive transport modeling of early burial dolomitization of carbonate platforms by geothermal convection // AAPG Bulletin. 2010. Vol. 94. P. 889–917.
