ПРИМЕНЕНИЕ ИМПУЛЬСНЫХ И ВЫСОКОЧАСТОТНЫХ МЕТОДОВ ЭЛЕКТРОННОГО ПАРАМАГНИТНОГО РЕЗОНАНСА ДЛЯ ИССЛЕДОВАНИЯ НЕФТЯНЫХ ДИСПЕРСНЫХ СИСТЕМ Текст научной статьи по специальности «Нанотехнологии»

обзорная статья

DOI: https://doi.org/10.18599/grs.2020.4.2-14 ~ УДК 538.955: 547.912

Применение импульсных и высокочастотных методов электронного парамагнитного резонанса для исследования

нефтяных дисперсных систем

М.Р. Гафуров1*, А.А. Пономарев2, Г.В. Мамин1, А.А. Родионов1, Ф.Ф. Мурзаханов1,

Т. Араш1, С.Б. Орлинский1

1 Казанский федеральный университет, Казань, Россия 2Тюменский индустриальный университет, Тюмень, Россия

Методами стационарного и импульсного электронного парамагнитного резонанса (ЭПР) в двух частотных диапазонах (X и W, с частотами микроволнового излучения около 9 ГГц и 95 ГГц, соответственно) исследованы спектральные и релаксационные характеристики «свободных» органических радикалов (СР) и ванадил-порфи-риновых комплексов (ВПК) в различных образцах нефтяных дисперсных систем (битумы, нефти, их высокомолекулярные компоненты и растворы). Рассмотрены отличительные особенности применения импульсных методов (электронное спиновое эхо, изучение модуляции спада сигнала электронного спинового эха, времен электронной релаксации) и высокочастотного ЭПР в изучении углеводородных систем. Применение ЭПР в W-диапазоне позволяет спектрально разрешить линии от разных парамагнитных центров, точнее определить их спектральные характеристики. Показано, что электронное спиновое эхо (в том числе и в сильных магнитных полях с индукцией 3.4 Тл) можно наблюдать и при комнатных температурах, что потенциально позволяет применять разнообразный арсенал методик импульсного ЭПР в полевых условиях и удешевить проведение импульсных измерений. Анализ кривой спада поперечной намагниченности на компонентах ВПК позволяет идентифицировать электрон-ядерные взаимодействия с ядрами окружения 14N и 1H в условиях in situ, в то время как в спектрах ЭПР указанные сверхтонкие взаимодействия обнаружить не удается. Выявлено, что форма линии сигнала ЭПР для СР в W-диапазоне не описывается изотропной одиночной линией, как в X-диапазоне. Описан эффект увеличения скоростей электронной поперечной релаксации в асфальтенах в рамках модели спектральной диффузии между быстро- и медленно релаксирующими парамагнитными центрами в надмолекулярных комплексах асфальтенов.

Ключевые слова: электронный парамагнитный резонанс, нефтяные дисперсные системы, асфальтены

Для цитирования: Гафуров М.Р., Пономарев А.А., Мамин Г.В., Родионов А.А., Мурзаханов Ф.Ф., Араш Т., Орлинский С.Б. (2020). Применение импульсных и высокочастотных методов электронного парамагнитного резонанса для исследования нефтяных дисперсных систем. Георесурсы, 22(4), с. 2-14. DOI: https://doi.org/10.18599/ grs.2020.4.2-14

1. Введение

С позиции коллоидной химии нефть может рассматриваться как сложная многокомпонентная смесь, которая в зависимости от внешних условий проявляет свойства молекулярного раствора или дисперсной системы (Сюняев, 1980; Сюняев и др., 1990; Сафиева, 2004). В объеме нефтяных дисперсных систем (НДС) распределены смолы, асфальтены, склонные к самоассоциации, и другие твердые частицы. Изучение свойств и комплексного поведения такой сложноорганизованной многокомпонентной системы продолжает оставаться актуальной и трудноразрешимой задачей (Сюняев, 1980; Сюняев и др., 1990; Сафиева, 2004).

НДС могут содержать до 1020 парамагнитных центров (ПМЦ) на 1 грамм вещества (Yen, Chilingarian, 1994, 2000). Большинство ПМЦ сконцентрировано в высокомолекулярных компонентах НДС, таких как асфальтены, смолы и полициклические ароматические углеводороды. Массовое содержание таких высокомолекулярных

* Ответственный автор: Марат Ревгерович Гафуров

E-mail: marat.gafuroy@kpfu.ru

© 2020 Коллектив авторов

компонент может достигать 45% для нефтей и 75% для битумов (Yen, Chilingarian, 1994, 2000). Это означает, что в асфальто-смолистых веществах (АСВ) со среднечисловой молекулярной массой 1000 а.е.м. на каждую молекулу может приходиться до одного неспаренного электрона (Yen, Chilingarian, 1994). Очевидно, такая высокая концентрация ПМЦ способна влиять как на магнитные, так и на электрофизические свойства вещества (Dolomatov et al., 2018). Исследование как углеводородных, так и других типов ПМЦ может предоставить дополнительные сведения о наличии и концентрации углеводородного сырья в нефтесодержащих породах, структуре и свойствах НДС, что обусловлено практической необходимостью эффективной разведки, добычи, транспортировки и более глубокой переработки трудноизвлекаемых углеводородных запасов - тяжелой (высоковязкой) нефти, битумов и др. (Martyanov et al., 2017).

Существование ПМЦ в НДС обуславливается наличием d-металлов (в основном V, Ni, Fe, Mn) и стабильных органических «свободных» радикалов (СР). Как правило, СР сконцентрированы в асфальтенах, благодаря наличию делокализованных п-электронов ароматических колец и углеродных радикалов боковых алкановых цепочек.

Точное положение и структура CP в НДС до сих пор являются предметом научных дискуссий. Присутствие гетероатомов (S, N, O) существенно расширяет номенклатуру ПМЦ, усложняя их однозначную идентификацию. В этом вопросе исследователи традиционно ссылаются на работы по исследованию ПМЦ в антраците (Yen, Chilingarian, 2000).

Часто, особенно в тяжелой (высоковязкой) нефти и битумах, детектируется наличие ванадия (до десятых долей процента). Считается установленным (Piccinato et al., 2012), что ванадий присутствует в НДС в основном в виде ванадил-порфириновых комплексов (ВПК) и сконцентрирован по большей части в АСВ с преобладанием в асфальтенах (рис. 1).

Рис. 1. Схематичная структура остова молекулы ВПК. Значения величины индукции внешнего магнитного поля, параллельного плоскости порфирина (Bt) или перпендикулярного ему (B2), определяются из спектров ЭПР

Агрегирование асфальтенов и образование асфаль-теновых отложений представляют собой серьезную проблему при добыче, транспортировке и переработке НДС, обуславливая интерес к изучению их структуры и свойств. Сведения об участии ВПК в процессах агрегирования асфальтенов достаточно противоречивы (Zhao et al., 2015). Разделяют «свободные» ВПК (имеющие небольшие размеры и простейшую структуру), которые достаточно легко удаляются из асфальтенов, и «сильносвязанные» ванадилпорфирины, которые связаны с молекулами асфальтенов либо ковалентно, либо посредством множественных межмолекулярных взаимодействий (Zhao et al., 2015; Абызгильдин и др., 1977). Ванадий же, в свою очередь, считается основным элементом, отравляющим катализаторы, металлом, мешающим крекингу нефти, вызывающим износ оборудования (при сгорании ванадий-содержащих топлив адгезионно- и коррозионно-активные неорганические соединения ванадия являются одной из главных причин интенсивного золового заноса и коррозии высокотемпературных поверхностей). Дополнительный интерес к изучению структурных и динамических характеристик петропорфиринов и других ПМЦ в последнее время связан с возможностью использования нативных и примесных ПМЦ в качестве поляризующих агентов для усиления сигнала ядерного магнитного резонанса (ЯМР) и ядерного магнитного каротажа (ЯМК) протонов НДС при помощи различных механизмов динамической поляризации ядер (ДПЯ) для целей геологоразведки и фундаментального изучения строения асфальтенов (Gizatullin et al., 2018, 2019; Alexandrov et al., 2019; Sapunov et al., 2016, 2019; Davydov et al., 2018; Zhang et al., 2020).

Электронный парамагнитный резонанс (ЭПР) давно известен как один из наиболее информативных методов, позволяющих детектировать и идентифицировать парамагнитные центры в НДС, измерить их концентрацию,

установить структуру парамагнитных комплексов (Khasanova et al., 2017). Метод не требует дополнительной пробоподготовки (например, путем разбавления, как для оптических измерений), является неразрушающим (в отличие, например, от хроматографических/масс-спектрометрических методов), позволяет исследовать твердые (битумы) и жидкие (нефть) образцы. ЭПР-исследования можно проводить как на компонентах НДС и их растворах, так и в условиях in situ, нативных НДС.

Основным инструментарием ЭПР-исследований на сегодняшний день являются спектрометры Х-диапазона с рабочей микроволновой частотой (сверхвысокая частота, СВЧ) v = 9-10 ГГц, оперирующие в стационарном режиме (СВЧ-излучение с небольшой мощностью постоянно подается на образец). Коммерчески выпускаются как приборы, работающие в других частотных диапазонах, так и использующие импульсные методики измерений (Qin, Wamcke, 2015a, 2015b). Удивительно, но технические возможности современных спектрометров недостаточно широко используются в исследованиях нефти, ее компонент и производных (по крайней мере, если судить по анализу открытых источников), в то время как активно применяются для расшифровки структур сложных биомолекул (Qin, Wamcke, 2015a, 2015b). До сих пор подавляющее большинство ЭПР-экспериментов на НДС ограничивалось методиками, впервые представленными в 1950-1970 гг. (Yen, Chilingarian, 1994; Zhang et. al., 2020; Trukhan et al., 2014; Wang et al., 2016; Dolomatov et al., 2016; Гарифьянов, Козырев, 1956; Gutowsky et al., 1958; O'Reilly, 1958; Ильясов и др., 1961; Ильясов, 1962; Yen et al., 1962; Муравьев и др., 2006; Yakubov et al., 2017). Лишь в последние годы возрос интерес к изучению не-фракционированных НДС методами высокочастотного и импульсного ЭПР (Ramachandran et al., 2015; Tayeb Ben et al., 2015; Mamin et al., 2016; Gracheva et al., 2016).

Целью данной работы являлось исследование возможностей импульсного (X- и W-диапазоны) и высокочастотного (W-диапазон) ЭПР для изучения НДС. Работа проведена в стенах Института физики Казанского федерального университета (г. Казань).

2. Материалы и методы

2.1. Материалы

Из множества НДС для представленного обзора были выбраны несколько образцов, сильно различающихся по вязкости и содержанию высокомолекулярных компонентов. В таблице 1 представлены реологические параметры и компонентный состав нефти с Ашальчинского месторождения (Республика Татарстан), Салымского месторождения (Ханты-Мансийский автономный округ, баженовский горизонт), Бока де Харуко (Республика Куба) и нефти с китайского месторождения, исследованные в данной работе. Детали проведения групповых и реологических измерений, результаты других аналитических методов исследования образцов приведены в работах (Gizatullin et al., 2018, 2019; Vakhin et al., 2020; Нестеров и др., 2019; Пономарев, 2019).

Объектами исследования были также 12 различных порошков асфальтенов, выделенных из технических битумов с двух нефтеперерабатывающих заводов (образцы N1 и N2), и тяжелых нефтей Аканского месторождения,

НАУЧНО-ТЕХНИЧЕСКИЙ ЖУРНАЛ

www.geors.ru ГЕОРЕСУРСЫ

Плотность при t = 20°С, кг/м3 Вязкость при 20°С, Пас Асфальтены, % Смолы, % Насыщенные и ароматические углеводороды, %

Ашальчинская 968 2.48 3.9 31.0 65.1

Салымская 827 0.005 0.1 3.3 96.4

Кубинская 1029 (t = 36°С) 270 17.3 42.8 39.9

Китайская 920 1.18 0.6 10.4 89

Табл. 1. Реологические параметры и компонентный состав образцов Ашальчинской, Салымской, кубинской и китайской нефти, исследованные в данной работе

расположенного на территории Нижневартовского района Ханты-Мансийского автономного округа (образцы N3 и N4). Начальная фракция асфальтенов (обозначаемая как АШ4) была получена из исходных образцов нефти и битумов экстракцией в аппарате Сокслета с помощью петролейного эфира (температура кипения 40-70°С). Далее (по схеме дробного осаждения асфальтенов из растворов в толуоле) исходные асфальтены АШ4 были разделены на две фракции: А1, выпадающая из толуольного раствора при добавлении 65% петролейного эфира как осадителя; и А2, выпадающая при максимальном (90%) добавлении петролейного эфира (табл. 2). Предполагается, что фракция А1 состоит преимущественно из молекул асфальтенов «континентального» типа, а А2 - типа архипелаг (Acevedo е! а1., 2010).

Образец Содержание Количество фракций

асфальтенов, экстрагированных

масс. % асфальтенов, масс. %

A1 A2

N1 34 83 15

N2 30 66 6.5

N3 9.2 60 7

N4 11.0 60 10

Табл. 2. Содержание асфальтенов в образцах N1-N4

2.2. Методы

Все описанные в данной работе измерения проведены с использованием возможностей спектрометров ЭПР X- и W-диапазонов фирмы Bruker (ФРГ) ESP-300 и ElexSys 680, расположенных в Институте физики КФУ Для проведения температурных измерений использовались азотные и гелиевые продувки, позволяющие изменять температуру на образце в диапазонах от 77 до 800 К и от 6 до 300 К, соответственно.

Основные элементы конструкции ЭПР-спектрометров стационарного режима остаются неизменными со времен установки Е.К. Завойского (Zavoisky, 1945). Обычно образцы НДС помещаются в кварцевую или стеклянную трубку, которая затем вставляется в резонатор ЭПр, находящийся во внешнем магнитном поле B0. Простейшим микроволновым резонатором может служить замкнутый с двух сторон металлический волновод прямоугольной или цилиндрической формы, геометрические размеры которого соответствуют длинам волн СВЧ. Фактически резонатор - это аналог колебательного контура, применяемого в ЯМР. Его предназначение - в накоплении энергии и разделении электрической и магнитной компонент сВЧ-поля, поскольку эффект ЭПр обусловлен взаимодействием именно с магнитной составляющей электромагнитной волны, в то время как электрическая компонента может вызывать нежелательный нагрев образца (Eaton et al., 2010). Стандартный размер образца для

исследования НДС в X-диапазоне - диаметром 3-4 мм и в высоту 4-30 мм. Измерения в W-диапазоне требуют гораздо меньшего объема вещества (трубки диаметром 0.1-0.8 мм, длина образца 3-10 мм), что можно рассматривать как дополнительное преимущество при работе с материалами ограниченного размера, исследовании негомогенных образцов, серий образцов НДС после применения (паро)теплового и каталитического воздействия, обработки критической водой и т.п. (Al-Muntaser et al., 2020; Mehrabi-Kalajahi et al., 2018; Galukhin et al., 2018; Mukhamatdinov et al., 2020).

Для повышения чувствительности спектрометров ЭПР постоянное магнитное поле модулируется синусоидальным сигналом с частотой 10-100 кГц и с малой амплитудой (0.1-1.0 мТл). При этом сигнал, измеряемый детектором, амплитудно-модулирован на частоте модуляции, и амплитуда его основной гармоники, как правило, пропорциональна первой производной линии поглощения. Таким образом, концентрация ПМЦ может быть определена путем двойного интегрирования записанного спектра.

На вид спектров как стационарного, так и импульсного ЭПР оказывают влияние несколько типов взаимодействий неспаренного электрона со своим окружением:

(1) зеемановское взаимодействие между неспаренным электроном и внешним магнитным полем,

(2) спин-орбитальное взаимодействие,

(3) электрон-ядерное или сверхтонкое взаимодействие (СТВ),

(4) взаимодействие с другими неспаренными электронами или спин-спиновое взаимодействие.

В X-, а особенно в W-диапазоне, наиболее значителен вклад зеемановского взаимодействия, определяющего положение линий резонанса (Yen, Chilingarian, 1994):

Чсвч = gPB, (1)

где h - постоянная Планка, g - фактор спектроскопического расщепления, р - магнетон Бора.

Спин-орбитальное и спин-спиновое взаимодействия ведут к отклонению спектроскопического g-фактора от известной величины для свободного электрона (g =2.0023). Типичный диапазон изменения g-факторов ПМЦ в НДС составляет от 1.96 (для VO2+) до 2.008 для кислород- и серосодержащих СР (Yen, Chilingarian, 1994, 2000; Khasanova et al., 2017). Для соблюдения условия резонанса (1) необходимо приложение B0 величиной около 0.35 Тл для Х-диапазона и 3.5 Тл для W-диапазона. Таким образом, высокочастотные исследования НДС требуют использования сверхпроводящих магнитов, что увеличивает стоимость проведения измерений. В ряде нефтесодержащих пород встречаются комплексы ионов группы железа Fe3+ с характерным сигналом в области «половинных» магнитных полей (g~4.3), а также сигналы от стабильных серо- и

кислород-содержащих радикалов (Khasanova et al., 2017).

Взаимодействие неспаренного электрона с ядерным магнитным моментом ядра (ядер) приводит к расщеплению электронных энергетических уровней, образуя так называемую сверхтонкую структуру (СТС) спектра ЭПР. Каждый электронный уровень m расщепляется на группу из (2I+1) близко расположенных уровней, где I - ядерный спин. Разрешенными переходами, дающими вклад в СТС спектра ЭПР, являются только те переходы, для которых выполняются условия Ams=±1 и AmI=0. Расстояние между соседними линиями СТС спектра называют константой СТВ, обозначаемой как А.

Еще в работе (Dickson et al., 1972) было показано, что значения параметров g и А могут быть использованы для уточнения природы парамагнитных центров в НДС и определения природы ВПК. Так, для лигандного окружения ванадия VOS4 значения g-фактора выше, а значения A меньше по сравнению с VON4, V0N202 или VOO4. Этот факт продолжает использоваться в ряде работ для контроля процессов десульфуризации НДС, например, в (Cui et al., 2017, 2018), хотя наблюдаемые в экспериментах X-диапазона изменения параметров спектров ЭПР крайне незначительны.

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

Импульсный режим позволяет зарегистрировать спектры, полученные либо по сигналу спада свободной индукции (ССИ, Free induction decay, FID), либо методом электронного спинового эха (ЭСЭ, Electron spin echo, ESE). Сигнал ССИ в большинстве НДС, как и в большинстве твердых тел, быстро затухает из-за дефа-зировки спинов в условиях обусловленного естественными причинами неоднородного уширения линий, либо вообще не наблюдается по причине "звона резонатора" (так называемое "мертвое время" прибора). Поэтому импульсные эксперименты проводятся в режиме низкой добротности резонатора, что снижает чувствительность ЭПР-измерений по сравнению со стационарным режимом.

Форма записанных таким образом спектров (амплитудные или интегральные значения сигналов ССИ или ЭСЭ) обычно совпадает со спектрами поглощения. В некоторых случаях (например, при изучении легких нефтей) для детектирования спектров ЭСЭ с приемлемым соотношением сигнал-шум может потребоваться замедление процессов электронной релаксации путем понижения температуры образца.

В состав импульсного спектрометра ЭПР входит дополнительный блок - формирователь импульсов,

соединенный с генератором СВЧ. По аналогии с ЯМР-применениями, СВЧ-импульсы могут быть объединены в специальные последовательности, время задержки и длительность импульсов в которых являются настраиваемыми параметрами. Однако значения этих времен для электронов на 3-6 порядков короче, чем в ЯМР, что связано с соответствующим различием времен релаксации электронов и ядер. Ниже представлены импульсные последовательности (с импульсами прямоугольной формы), использованные нами в изучении НДС, и приведены численные значения их параметров для экспериментов в X-диапазоне.

1) Двухимпульсная последовательность п/2-т-п с длительностью п/2-импульса 16 нс и временем задержки т=200 нс использовалась для детектирования ЭСЭ.

2) Сигнал спада амплитуды первичного ЭСЭ и время поперечной релаксации Te измерялись по представленной выше двухимпульсной схеме при увеличении интервала т между импульсами с шагом Дт от 4 нс и фиксированных длительностях импульсов. Для регистрации модуляции огибающей электронного спинового эха (Electron Spin Echo Envelope Modulation, ESEEM) использовалась трехимпульсная последовательность п/2-(т1+Дт)-п/2-(т2+Дт)-п/2 (Dikanov, Tsvetkov, 1992; Дзюба, 2013), где т2=5300 нс.

3) Кривая восстановления продольной намагниченности и время продольной релаксации T1e измерялись методом инверсии-восстановления с помощью последовательности п-T -п/2- т -п-ЭСЭ при фиксированной

задержки it ±

длительности импульсов и времени т, но изменяющемся времени T .

задержки

3 результаты и обсуждение

3.1. ЭПр образцов нефти

На рисунке 2 представлены спектры стационарного ЭПР образца кубинской нефти, зарегистрированные при

СР

1

3.32 3.34 3.36 3.38 3.40

B(T)

9.6 ГГц

Рис. 2. Спектры стационарного ЭПР породы нефти с месторождения Бока де Харуко в двух частотных диапазонах, обусловленные сигналами СР, ВПК и комплексами двухвалентного марганца. Подробно вид и спектроскопические параметры указанных ПМЦ обсуждаются в тексте работы

Mn

СР

280

300

400

B0, мТл

НЮЧНО-ТЕХНИЧЕСЩЙ ЖУРНАЛ

www.geors.ru ГЕОРЕСУРСЫ

комнатной температуре в X- и W-диапазонах, как наиболее комплексные среди исследованных материалов. Спектры обусловлены (1) одиночной линией, обязанной сигналу от СР ^=2.003, ширина линии от пика до пика ДНрр~6.4 Гс); (2) линиями СТС ванадилового комплекса У02+ с g-фактором аксиальной симметрии (1=7/2 для ядер 51У ^,=1.963, g1=1.985, А=468 МГц, А±=150 МГц, раздел 3.2); (3) СТС от ионов Мп2+ в карбонатной породе с изотропным g-фактором (1=5/2 для ядер 55Мп, g = 2.001, А. =260 МГц; наблюдается расщепление на центральном переходе тонкой структуры с Дотх=(-1/2) - (+1/2)). Из-за небольшого различия значений g-факторов для приведенных спиновых систем сигналы от указанных ПМЦ частично спектрально разрешены в W-диапазоне. В образцах Салымского месторождения и китайской нефти нами зарегистрированы ЭПР-сигналы только от СР (одиночные линии). Ниже подробно будут рассмотрены сигналы от ВПК и СР, характерные для большинства изученных нами НДС. В таблице 3 представлены данные по концентрации указанных выше ПМЦ, извлеченные из измерений в стационарном режиме Х-диапазона.

3.2. Изучение спектральных и релаксационных характеристик сигналов от ВПК

Спектральные характеристики ВПК аксиальной симметрии в изученных образцах приведены выше (раздел 3.1) и практически идентичны для всех исследованных образцов. Достичь изотропного спектра ЭПР для ВПК в НДС возможно путем растворения образцов в органических растворителях или повышения температуры/ давления (Тги^ап е! а1., 2014). На рисунке 3 представлена температурная зависимость спектра стационарного ЭПР ВПК образца Ашальчинской нефти в Х-диапазоне. При Тнач<480 К наблюдается порошковый спектр (16 линий), который может быть описан в рамках вращательной модели с временем вращательной корреляции твр>5-10"7 с. При Ткон>550 К в образце детектируется только изотропный спектр (8 линий) с твр~10"10 с и изотропным значением Л—300 МГц ^аШгоу е! а1., 2018). При этом процесс является обратимым по температуре, а значения компонент g-фактора не изменяются в пределах погрешности измерений.

Допуская справедливость формулы Смолуховского-Стокса-Эйнштейна-Дебая для НДС (Тги^ап е! а1., 2014):

Рис. 3. Изменения спектра ЭПР ВПК в образце Ашальчинской нефти от температуры для Т=295—661 К в Х-диапазоне

tbd = 4pR3n(T) / 3kBT,

(2)

где R - характеристический размер молекул, п(Т) - вязкость НДС, зависящая от температуры, кв - постоянная Больцмана, и используя экспериментальные значения п(Т) из работы (Пут е! а1., 2016), можно с помощью уравнения (2) оценить характеристические размеры молекул, содержащих ВПК, как 0.23±0.05 нм при высоких температурах и 1.0±0.1 нм при низких. То есть характеристические размеры ВПК при данном переходе изменяются в 6-8 раз, что хорошо согласуются с принятой схемой агрегации молекул асфальтенов (МиШш е! а1., 2014). Полученные результаты свидетельствуют в пользу модели участия вхождения, регистрируемых методом ЭПР, ВПК в структуры асфальтеновых агрегатов. Мы не проводили подобных температурных измерений в W-диапазоне ввиду отсутствия высокотемпературного резонатора.

Необходимо отметить, что в спектрах ЭПР НДС не удается наблюдать сигналы от комплексов никеля (никель-порфиринов), несмотря на их распространенность в НДС (ГакиЬоуа е! а1., 2019). Несколько предположений могут объяснить этот факт. Во-первых, если никель находится в тетраэдрической конфигурации, сигнал ЭПР может не детектироваться. Во-вторых, если никель находится в октаэдрическом окружении, которое сильно искажено в

Концентрация Концентрация Параметры описания линии Параметры описания

СР, спин/г ВПК, спин/г СР в Х-диапазоне линии СР в '-диапазоне

g = 2.0036;

Ашальчинская (1.1±0.2) 1017 (1.2±0.2) 1018 С:Ь = 61:39; ДИрр = 0.64 мТл

Не обнаружено g = 2.0023; = 2.00235,

Салымская (8.0±0.9) 1018 С:Ь = 53:47; g± = 2.00145

ДИрр = 0.31 мТл g = 2.0018; ДИрр(Ь) = 0.38 мТл

Кубинская (4.3±0.8) 1016 (1.4±0.2) 1017 С:Ь = 34:66; ДИрр = 0.55 мТл

(2.8±0.4) 1017 Не обнаружено g = 2.0042; = 2.0047,

Китайская С:Ь = 19:81; g± = 2.0041

ДИрр = 0.51 мТл ДИрр(Ь) = 0.75 мТл

Табл. 3. Концентрации СР и ВПК в исследованных образцах нефти и параметры аппроксимации спектров ЭПР для СР в Х-диапазоне как свертки гауссовой (G) и лоренцевой (Ь) форм линий и центром аксиальной симметрии в Ж-диапазоне

осевом направлении (случай порфириновых комплексов), величина начального расщепления в нулевом магнитном поле (ZFS, zero-field splitting) может намного превышать значение для тепловой энергии kT, так что заселенным оказывается только основное состояние с ms=0 (Tayeb Ben et al., 2015).

На рисунке 4 приведен импульсный спектр ЭПР в X-диапазоне для нефти Ашальчинского месторождения и его симуляция с использованием программы EasySpin для MatLab (Stoll, Schweiger, 2006) с параметрами g±=1.9860; g =1.9661; A±=160 МГц; A =470 МГц. Для детектирования ЭСЭ-спектра ВПК в исходной НДС потребовалось понижение температуры до T=250 K, в то время как ЭСЭ в асфальтенах и смолах, выделенных из исходной нефти, а также сигнал от СР в исходном образце, регистрируются и при комнатных температурах.

В таблице 4 представлены значения времен релаксации для СР и ВПК, демонстрирующие, что времена продольной релаксации ВПК в X-диапазоне при температурах, близких к комнатным, на 1-2 порядка короче аналогичных величин для СР, поэтому в ряде образцов сигнал ЭСЭ от ВПК при комнатных, а тем более повышенных температурах (особенно в сильных магнитных полях) может и не наблюдаться. Полученные нами данные для большого ряда изученных образцов контрастируют с результатами работы (Raghunathan, 1991), в которой исследовались асфальтены, выделенные из битума месторождения Атабаски (Канада). В работе (Raghunathan, 1991) было показано, что ЭСЭ для ВПК в исследованных асфальтенах в X-диапазоне наблюдается только при T<20 K, что требует применения дорогостоящего оборудования и хладагентов для достижения таких низких температур. Мы полагаем, что данное утверждение связано, скорее всего, с несовершенством самодельной техники импульсного ЭПР, использованной в работе (Raghunathan, 1991).

При понижении температуры (с целью замедления процессов релаксации) на компонентах спектра ВПК (но не для СР) возможно наблюдение осцилляций, вызванных электрон-ядерными взаимодействиями (рис. 5). Фурье-преобразованный сигнал трехимпульсного варианта ESEEM (Electron Spin Echo Envelope Modulation,

ВПК(Б1)

ЭСЭ ЭПР T = 250 K

exp

280 300 320 340 360 380 400 420 B, мТл

Рис. 4. Спектр ЭПР образца Ашальчинской нефти при Т=250 К, зарегистрированный по амплитуде ЭСЭ (ехр), и симуляция спектра ВПК ^т). Показаны вклады от компонент СТС с различными значениями т1

ВПК (Bi), мкс СР (B2), мкс

Tie 1.7(2) 23.6(8)

T2e 0.31(4) 0.30(4)

(Dikanov, Tsvetkov, 1992; Дзюба, 2013; Schweiger, Jeschke, 2001) представлен на рисунке 6. В спектре ESEEM отчетливо проявляются сигналы, обязанные взаимодействию электрона с протонами (1=1/2, сигнал на ларморовской частоте 14.7 МГц) и частично разрешенные пики в диапазоне частот 0-10 МГц, приписываемые взаимодействию c ядрами 14N с I=1 (Deligiannakis et al., 2000), в то время как в спектрах ЭПР изученных нами образцов расщеплений, вызванных этими ядрами, не наблюдалось (работа (Гилинская, 2008), где обсуждается редкая возможность наблюдения расщеплений линий ЭПР в X-диапазоне, вызванных ядрами азота). Таким образом, применение техники ESEEM позволяет определить тип лигандного окружения (в данном случае наличие ядер 14N) для вана-дил-порфириновых комплексов в НДС без их экстракции.

Относительные интенсивности и точные положения пиков на рисунке 6 зависят от выбора ориентации (значения B0) и типа (типов) преобладающих ВПК. Сложный вид спектра для 14N обусловлен наличием не только сверхтонкого, но и квадрупольного взаимодействия, характеризуемого тензором Q. Расшифровке результатов ESEEM-исследований модельных (синтетических) порфиринов посвящена работа (Reijerse et al., 1998). определению компонент тензора Q методом двойного электрон-ядерного резонанса (ДЭЯР) непосредственно в образцах НДС без их обработки и экстракции посвящено исследование (Gracheva et al., 2016). В настоящей работе методы двойных резонансов для изучения структуры и динамики НДС не рассматриваются.

На рисунке 7 представлен типичный спектр ЭСЭ асфальтеновой фракции, зарегистрированный в W-диапазоне при T=100 К (известно, что T в твердых

X-диапазон T= 100 K

Ф 0,75^ I, B0=B1= 345.1 мТл

О О

с

п м

С

4 6

2т, мкс

Рис. 5. Сигнал спада амплитуды ЭСЭ для ВПК в магнитном

поле B=345.1 мТл

0 5 10 15 20 25 30 Частота, МГц

Табл. 4. Значения Т1е и Т2е из импульсных измерений при Т=250 К в Ашальчинской нефти в магнитных полях В1 и В2 (рис. 4)

Рис. 6. Амплитудное Фурье-преобразование временной зависимости трехимпульсного ESEEM для ВПК в магнитном поле

B1=345 мТл

НАУЧНО-ТВШИЧЕСКИЙ ЖУРНАЛ

www.geors.ru ГЕОРЕСУРСЫ

VO2+

i. W-диапазон. T = 100 K

СР 1

J ß \ Эксперимент

:у Симуляция VO+СР 4-ч

B, мТл

Рис. 7. Спектр ЭПР асфальтенов Ашальчинской нефти в импульсном режиме в Ж-диапазоне при Т=100 К и Т =500 мкс

у г г повторения

и симулированный спектр, полученный сложением порошкового спектра VO2+ ^=1.984, g1=1.984, А]=16.8 тТ, А±=6.0 тТ) и одиночной линии свободного радикала ^=2.0036). Маркеры СР и VO2+ указывают на значения В при которых измерялись времена электронной релаксации для СР и ВПК

телах уменьшается с увеличением значения внешнего магнитного поля (Абрагам, Блини, 1973); поэтому для регистрации ЭСЭ для ВПК в W-диапазоне иногда требуется применение более низких по сравнению с Х-диапазоном температур). Значение времени повторения импульсной последовательности Т =500 мкс выбрано с целью

повторения г

увеличения интенсивности компоненты ВПК спектра относительно СР за счет разности во временах продольной релаксации Т1 (ср. с данными табл. 4). Небольшое отличие g-факторов СР и ВПК позволяет также спектрально разрешить эти ПМЦ на высоких частотах (рис. 7). Кроме того, благодаря более высокому спектральному разрешению, значения компонент g-тензора для ВПК можно

3,35 3,40 3,45 3,50

Магнитное поле, Тл

Рис. 8. Энергетические уровни и соответствующий спектр поглощения ЭПР для ВПК, рассчитанные для микроволновой частоты п=94 ГГц, g|=1.963, gi=1.985, А^ =470 МГц, Л±=150МГц. Отдельные вклады от каждого ЭПР-перехода выделены различными цветами

определить с лучшей точностью. Рисунок 8 поясняет вид порошкового спектра ЭПР для ВПК в W-диапазоне, стрелочками показаны значения B0, соответствующие значениям g и g± для ВПК.

3.3 сигналы ср

Известно, что спектры СР в НДС в Х-диапазоне в большинстве случаев представляют собой одиночную линию с шириной AHpp=(0.4-0.7) мТл и изотропным g-фактором в диапазоне 2.003±0.001 (Yen, Chilingarian, 1994, 2000; Gizatullin et al., 2018), которая стандартно описывается комбинацией гауссовой (G) и лоренцевой (L) форм линий с одинаковыми g-факторами и примерно одинаковой шириной (Mukhamatdinov et al., 2018). Для примера, на рис. 9 представлена аппроксимация линии ЭПР в программе Origin функцией Фойгта (свертка лоренцевой и гауссовой форм линий) для нефти Салымского месторождения. Численные параметры аппроксимаций для СР приведены в табл. 3.

Наш опыт исследований различных НДС в W-диапазоне показывает (Volodin et al., 2013), что спектры СР в подавляющем большинстве случаев имеют несимметричную форму линии (левое, слабополевое крыло «затянуто»). Таким образом, типично принятое описание спектров ЭПР для СР в X-диапазоне не позволяет описать спектр ЭПР на более высоких частотах.

На асимметричность линии ЭПР «свободных» радикалов в асфальтенах с арабских и колумбийских месторождений в стационарном режиме W-диапазона обратили внимание авторы работы (Di Mauro et al., 2005). Им не удалось описать экспериментальные данные единственным ПМЦ с анизотропным g-фактором, и, как следствие, было предположено, что спектр ЭПР обязан наличию двух (или более) различных ПМЦ с близкими по значению g-факторами, линии ЭПР которых частично разрешены в сильных магнитных полях. Соглашаясь с такой интерпретацией, необходимо отметить, что для некоторых образцов НДС возможно описание сигнала СР единственным ПМЦ.

Например, при исследовании образцов Салымской и китайской нефти нам удалось описать спектры СР в двух частотных диапазонах единым набором параметров для центра аксиальной симметрии (табл. 3, рис. 10). Дополнительными аргументами в пользу описания сигнала ЭПР единственным ПМЦ являются (1) практическая независимость релаксационных характеристик (кривых Т1СР и Т2СР) от выбора точного значения магнитного поля вдоль спектра ЭПР и (2) независимость вида спектра

X-диапазон д

\ exp Г К-——' sim

W-диапазон T = 298 K

339 340 341 342 343 344 345 346 347 Магнитное поле, мТл

3354 3356 3358 3360 3362 Магнитное поле, мТл

Рис. 9. Слева - спектр ЭПР образца Салымской нефти ЭПР в Х-диапазоне (ехр) и его аппроксимация формой линии Фойгта (вт) с параметрами, описанными в тексте работы (табл. 3). Справа - вид спектра ЭПР данного образца в Ж-диапазоне

3350

3400

3450

3500

Салымская нефть ЭСЭ W-диапазон T=300K

Китайская нефть ЭСЭ W-диапазон T=300K

33500 33550 33600

Магнитное поле, мТл

3346 3348 3350 3352 Магнитное поле, мТл

Рис. 10. Эхо-детектированные спектры ЭПР образцов Са-лымской (слева) и китайской (справа) нефти в Ш-диапазоне при комнатной температуре и их описание центром аксиальной симметрии с £ц и g±, значения которых приведены в табл. 3. Экспериментальные данные обозначены как ехр, кривые аппроксимации как sm

ЭПР от изменения температуры (в диапазоне 20-300 К). Обращает на себя внимание тот факт, что для всех исследованных в данной работе образцов нефти (не говоря уже об асфальтенах) сигнал ЭСЭ в W-диапазоне можно наблюдать при комнатных температурах.

Нам не удалось удовлетворительно описать единственным ПМЦ аксиальной и даже ромбической симметрии спектры СР для Ашальчинской и кубинской нефти. Обращает на себя внимание тот факт, что в этих образцах, в отличие от Салымской и китайской нефти, имеются и другие типы ПМЦ (рис. 2, табл. 3). Известно также, что увеличение концентрации гетероатомов (в первую очередь кислорода и серы) приводит к смещению линий ЭПР в более слабые магнитные поля (Gizatullin et а1., 2018; Khasanova et а1., 2017; и цитируемая литература). Детальных исследований влияния гетероатомов и других ПМЦ на спектральные характеристики СР в данной работе не проводилось.

3.4. Измерение времен электронной релаксации ВПК и ср во фракциях асфальтенов в W-диапазоне

Как уже отмечалось выше, позиция, форма и интенсивность линий ЭПР свободного радикала и ВПК в экспериментах по стационарному ЭПР в Х-диапазоне, а также их соотношение являются своеобразными метками для классификации углеводородного сырья из различных месторождений и контроля процессов воздействия на НДС ^еп, ШИ^апап, 1994, 2000; МеЬгаЫ-Ка^аЫ et а1., 2018). Логично предположить, что и электронные релаксационные характеристики являются параметрами, чувствительными к изменению химического и координационного окружения ПМЦ.

Объектами исследования были 12 различных порошков асфальтенов (раздел 2.1), выделенных из технических битумов и тяжелой нефти. Измерения проводились при Т=300 К. Как и для случая исследований в Х-диапазоне, было обнаружено, что для СР времена продольной релаксации Т1СР находятся в интервале 20-100 мс, значительно изменяясь от образца к образцу, в то время как для ВПК значения Т1ВПК короче и лежат в пределах 1.0-1.8 мс (рис. 11).

Времена поперечной релаксации ВПК в асфальтенах также существенно короче, чем для СР. Во всех образцах, не учитывая эффектов модуляции, наблюдался моноэкспоненциальный распад поперечной намагниченности с характерными временами Т2ВПК в диапазоне (80-220) нс.

Характер распада поперечной намагниченности Т2СР

100 80 60 40 20

б)

Рис. 11. Времена продольной релаксации Т1 для СР (а) и У02+ (б) при Т=297 К в Ш-диапазоне для 12 исследованных образцов асфальтенов

« 101

СО О

о

о я

х х о я

и К

5

10

10

10

500

1000

2500

3000

3500

1500 2000 2т, нс

Рис. 12. Сигнал спада амплитуды ЭСЭ для СР в Ш-диапазоне для образца №4 при комнатной температуре. Точки представляют собой экспериментальные данные, сплошными линиями показаны результаты аппроксимации, используя выражение (2) с т=0 в растворе толуола (верхняя линия, растворение в соотношении А1:толоул=1:10 по объему) и т=3.6*10-6 нс-2 для нерастворенной фракции (нижняя линия)

отличается от одноэкспоненциального (рис. 12) и для всех исследованных образцов хорошо описывается выражением (Матш et а1., 2016): / ~ >

1есЪо=МРК-ЫР

2 т

12 т у

■ехр(-/и-т2),

(3)

где М - коэффициент, пропорциональный концентрации СР, т - параметр, учитывающий спектральную диффузию в СР. В таблице 5 приведены данные по измерению концентрации ПМЦ в Х-диапазоне.

Все полученные релаксационные данные хорошо группируются около прямых (рис. 13):

т Гюр (3.0 ± 0.4уг'зю,

Т'гср ос (0.43 ± 0М)-Г'2уО-

(4)

(5)

Из анализа времен релаксации установлено (Матш et а1., 2016), что ускорение процесса поперечной релаксации

НАУЧНО-ТЕХНИЧЕСКИЙ ЖУРНАЛ

www.geors.ru ГЕОРЕСУРСЫ

3342 3344

3354

Образец Ainit A1 A2

N1 1.9 ± 0.1 3.6 ± 0.2 1.4 ± 0.1

N2 49 ± 2 60 ± 3 28 ± 2

N3 0.73 ± 0.05 0.7 ± 0.1 0.7 ± 0.1

N4 1.15 ± 0.09 1.0 ± 0.1 0.50 ± 0.05

Табл. 5. Относительные концентрации СР и ВПК (ССР/СВПК) в исследованных образцах асфальтенов (из измерений в X-диапазоне при T=300 К)

согласно (3) вызвано наличием в непосредственной близости (1-3 нм) от CP молекул ванадил-порфирина. Данное открытие, по нашему мнению, представляет собой еще один шаг к пониманию структуры асфальтенов и роли, которую комплексы ванадия играют в их агрегации. Отметим, что при растворении образцов асфальтенов в органических растворителях (например, в толуоле в соотношении асфальтен:толуол >1:10 по объему) кривые распада намагниченности приобретают типичный моноэкспоненциальный вид (рис. 12). Полученные данные еще раз свидетельствуют в пользу необходимости изучения НДС в условиях, близких к нативным (in situ), а не в растворах или с использованием модельных образцов -получаемые динамические характеристики на модельных образцах и растворах могут сильно отличаться от таковых для реальных систем. Здесь уместно снова подчеркнуть особенности (преимущества) методик ЭПР, позволяющих исследовать нативные НДс, их компоненты и растворы.

Подобные «ускорения» кривых T2CP не были описаны для НДС до наших экспериментов (Mamin et al., 2016). Тот факт, что нам не удалось наблюдать их в Х-диапазоне, по-видимому, связан с перекрыванием спектров CP и ВПК на более низких частотах (ср. рисунки 2, 4 и 7).

заключение

В настоящей работе мы попытались продемонстрировать некоторые возможности импульсного и высокочастотного ЭПР для изучения нефтяных дисперсных систем. Среди результатов работы можно выделить следующее.

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

Анализ кривой спада поперечной намагниченности (ESEEM) на компонентах парамагнитных ванадил-порфи-риновых комплексов НДС позволяет идентифицировать электрон-ядерные взаимодействия с ядрами окружения 14N и 'H в условиях in situ, в то время, как в спектрах ЭПР нативных образцов НДС указанные сверхтонкие взаимодействия обнаружить не удается из-за неоднородного уширения линий.

Применение высокочастотного ЭПР позволяет (по крайней мере, частично) спектрально разделить вклады

40

30

о 20

10

0

4

3

2

1

Т -i

'2 VO ,

9

106c-1

Рис. 13. Зависимость произведения времени продольной релаксации СР на коэффициент спектральной диффузии от скорости поперечной релаксации ванадил-порфиринов (красные квадратики, левая ось ординат) и зависимость скорости поперечной релаксации СР от скорости поперечной релаксации ванадил-порфиринов (открытые кружки, правая ось ординат) для 12 исследованных образцов

от различных ПМЦ, точнее определить их спектральные и релаксационные характеристики.

Обнаружено, что форма линии сигнала ЭПР «свободных» органических радикалов в исследованных образцах нефти и асфальтенов в W-диапазоне не описывается изотропной одиночной линией, как в X-диапазоне. Таким образом, результаты и выводы ряда работ по анализу изменения формы линии ЭПР в НДС в X-диапазоне при внешнем воздействии должны быть подвергнуты критическому переосмыслению.

Вместе с тем, в отличие от результатов работы (Di Mauro et al., 2005), показано, что в ряде НДС линии ЭПР «свободных» органических радикалов могут быть приписаны одиночным ПМЦ аксиальной симметрии.

Используя возможности спектрального разрешения линий различных ПМЦ в НДС в сильных магнитных полях (~3.4 Тл, W-диапазон), обнаружен эффект увеличения скоростей электронной поперечной релаксации нативных для асфальтенов НДС парамагнитных центров. Полученные данные интерпретированы в рамках модели спектральной диффузии между быстро- и медленно релаксирующими ПМЦ в надмолекулярных комплексах асфальтенов.

Авторы выражают надежду, что дальнейшие исследования смогут более полно раскрыть динамические и структурные свойства изучаемых объектов, в том числе используя и другие техники импульсного ЭПР (Qin, Warncke, 2015a, 2015b).

Благодарности/Финансирование

Авторы благодарят Д.Т. Ситдикова, М.А. Володина, А.В. Вахина, Ю.М. Ганееву и Т.Н. Юсупову за помощь в проведении измерений и подготовке образцов, Т.Б. Биктагирова за ряд проведенных теоретических расчетов, В.И. Морозова за обсуждение полученных результатов, А.В. Дуглава за внимательное прочтение рукописи и ценные комментарии.

Авторы благодарят рецензентов за возможность дополнить и исправить работу.

Работа выполнена при поддержке гранта РНФ 19-12-00332.

0

Литература

Абрагам А., Блини Б. (1973). Электронный парамагнитный резонанс переходных ионов. Том 2. Перевод с англ. М.: Мир, 349 с.

Абызгильдин Ю.Н., Михайлюк Ю.Н., Яруллин К.С., Ростовская А.А. (1977). Порфирины и металлопорфириновые комплексы нефтей. М.: Наука, 88 с.

Гарифьянов Н.С., Козырев Б.М. (1956). Парамагнитный резонанс в антраците и других содержащих углерод веществах. ЖЭТФ, 30(2), с. 272-276.

Гилинская Л.Г. (2008). Спектры ЭПР комплексов V (IV) и структура нефтяных порфиринов. Журнал структурной химии, 49(2), с. 245-254.

Дзюба С.А. (2013). Изучение структуры биологических мембран с помощью ESEEM спектроскопии спиновых меток и дейтериевого замещения. Журнал структурной химии, 54(S1), с. 5-18.

Ильясов А.В. (1962). Определение содержания ванадия в нефтях и нефтепродуктах методом ЭПР. Химия и технология топлив и масел, 9, с. 63-67.

Ильясов А.В., Гарифьянов Н.С., Рыжманов Ю.С. (1961). Электронный парамагнитный резонанс в некоторых сортах природной нефти и ее тяжелых фракциях. Химия и технология топлив и масел, 1, с. 28-31.

Муравьев Ф.А., Винокуров В.М., Галеев А.А., Булка Г.Р., Низамутдинов Н.М., Хасанова Н.М. (2006). Парамагнетизм и природа рассеянного органического вещества в пермских отложениях Татарстана. Георесурсы, 2(19), с. 40-45.

Нестеров И.И., Александров В.М., Пономарев А.А., Заватский М.Д., Лободенко Е.И., Кобылинский Д.А., Кадыров М.А. (2019). Экспериментальные исследования радикальных реакций преобразования углеводородного сырья. Известия высших учебных заведений. Нефть и газ, 4, с. 57-69. https://doi.org/10.31660/0445-0108-2019-4-57-69 Пономарев А.А. (2019). Механизм крекинга углеводородов в электромагнитных полях-к вопросу об образовании баженовской нефти. Известия высших учебных заведений. Нефть и газ, 1, с. 14-18. https:// doi.org/10.31660/0445-0108-2019-1-14-18

Сафиева Р.З. (2004). Химия нефти и газа. Нефтяные дисперсные системы: состав и свойства (часть 1). М.: РГУ нефти и газа им. И.М. Губкина, 112 с.

Сюняев З.И., Сафиева Р.З., Сюняев Р.З. (1990). Нефтяные дисперсные системы. М.: Химия, 224 с.

Сюняев З.И. (1980). Концентрация сложных структурных единиц в нефтяных дисперсных системах и методы ее регулирования. Химия и технология топлив и масел, 16(7), с. 484-489.

Acevedo S., Guzman K., Ocanto O. (2010). Determination of the number average molecular mass of asphaltenes (Mn) using their soluble A2 fraction and the vapor pressure osmometry (VPO) technique. Energy & Fuels, 24(3), pp. 1809-1812. https://doi.org/10.1021/ef9012714

Al-Muntaser A.A., Varfolomeev M.A., Suwaid M.A. et al. (2020). Hydrothermal upgrading of heavy oil in the presence of water at sub-critical, near-critical and supercritical conditions, Journal of Petroleum Science and Engineering, 184, 106592. https://doi.org/10.1016/j.petrol.2019.106592

Alexandrov A.S., Ivanov A.A., Archipov R.V., Gafurov M.R., Tagirov M.S. (2019). Pulsed NMR spectrometer with dynamic nuclear polarization for weak magnetic fields. Magnetic Resonance in Solids, 21(2), pp. 19203 (1-6). https://doi.org/10.26907/mrsej-19203

Cui Q., Ma X., Nakano K., Nakabayashi K., Miyawaki J., Al-Mutairi A. et al. (2018). Hydrotreating reactivities of atmospheric residues and correlation with their composition and properties. Energy & Fuels, 32(6), pp. 6726-6736. https://doi.org/10.1021/acs.energyfuels.8b01150

Cui Q., Nakabayashi K., Ma X., Ideta K., Miyawaki J., Marafi A.M. et al. (2017). Examining the molecular entanglement between V=O complexes and their matrices in atmospheric residues by ESR. RSC advances, 7(60), pp. 37908-37914. https://doi.org/10.1039/C7RA06436E

Davydov V.V., Dudkin V.I., Myazin N.S., Rud' V.Yu. (2018). Peculiarity of the Nuclear Magnetic Resonance Method Application for the Liquid Medium Flow Parameters Control. Applied Magnetic Resonance, 49(7), pp. 665-679. https://doi.org/10.1007/s00723-018-0994-1

Deligiannakis Y., Louloudi M., Hadjiliadis N. (2000). Electron spin echo envelope modulation (ESEEM) spectroscopy as a tool to investigate the coordination environment of metal centers. Coordination Chemistry Reviews, 204(1), pp. 1-112. https://doi.org/10.1016/S0010-8545(99)00218-0 Di Mauro E., Guedes, C.L.B., Nascimento O.R. (2005). Multifrequency (X-band to W-band) CW EPR of the organic free radical in petroleum asphaltene. Applied Magnetic Resonance, 29(4), pp. 569-575. https://doi. org/10.1007/BF03166333

Dickson F.E., Kunesh, C.J., McGinnis E.L., Petrakis L. (1972). Use of electron spin resonance to characterize the vanadium (IV)-sulfur species in petroleum. Anal. Chem., 44(6), pp. 978-981. https://doi.org/10.1021/ac60314a009

Dikanov S.A., Tsvetkov Y.D. (1992). Electron Spin-Echo Envelope Modulation (ESEEM) Spectroscopy. USA: CRC Press, 432 p.

Dolomatov M., Gafurov M., Rodionov A., Mamin G., González L.M., Vakhin A., Petrov A., Bakhtizin R., Khairudinov I., Orlinskii, S. (2018). Low-temperature thermal decomposition of heavy petroleum distillates: interconnection between the electrical properties and concentration of paramagnetic centres. IOP Conf. Ser.: Earth Environ. Sci., 155, 012007. https://doi.org/10.1088/1755-1315/155/1Z012007

Dolomatov M.U., Rodionov A.A., Gafurov M.R., Petrov A.V., Biktagirov T.B., Bakhtizin R.Z., Makarchikov S.O., Khairudinov I.Z. and Orlinskii S.B. (2016). Concentration of paramagnetic centres at low-temperature thermal destruction of asphaltenes of heavy petroleum distillates. Magnetic Resonance in Solids, 18, 16101. http://mrsej.kpfu.ru/contents.html#16101

Eaton, G.R., Eaton, S.S., Barr, D.P., Weber, R.T. (2010). Quantitative EPR. Vienna: Springer-Verlag Wien. https://doi.org/10.1007/978-3-211-92948-3 Gafurov M.R., Volodin M.A., Rodionov et al. (2018). EPR study of spectra transformations of the intrinsic vanadyl-porphyrin complexes in heavy crude oils with temperature to probe the asphaltenes' aggregation. Journal of Petroleum Science and Engineering, 166, pp. 363-368. https:// doi.org/10.1016/j.petrol.2018.02.045

Galukhin A., Bolmatenkov D., Osin Y. (2018). Heavy oil oxidation in the nano-porous medium of synthetic opal, RSC Adv., 8, pp. 18110-18116. https://doi.org/10.1039/C8RA02822B

Gizatullin B., Gafurov M., Vakhin et al. (2019). Native Vanadyl Complexes in Crude Oil as Polarizing Agents for In Situ Proton Dynamic Nuclear Polarization. Energy & Fuels, 33(11), pp. 10923-10932. https://doi. org/10.1021/acs.energyfuels.9b03049

Gizatullin B., Gafurov, M., Rodionov A., Mamin, G., Mattea, C., Stapf, S., Orlinskii, S. (2018). Proton-Radical Interaction in Crude Oil - A Combined NMR and EPR Study. Energy & fuels, 32(11), pp. 11261-11268. https://doi. org/10.1021/acs.energyfuels.8b02507

Gracheva I, Gafurov M., Mamin G., Biktagirov T., Rodionov A., Galukhin V., Orlinskii S.B. (2016). ENDOR Study of Nitrogen Hyperfine and Quadropole Tensors in Vanadyl Porhyrins of Heavy Crude Oils. Magnetic Resonance in Solids, 18, 16102. http://mrsej.kpfu.ru/contents.html#16102

Gutowsky H., Roger Ray B., Rutledge R., Unterberger R. (1958). Carbonaceous Free Radicals in Crude Petroleum. J. Chem. Phys., 28, pp. 744-745. https://doi.org/10.1063/L1744250

Ilyin S.O., Arinina M.P., Polyakova M.Y. et al. (2016). Rheological comparison of light and heavy crude oils. Fuel, 186, pp. 157-167. https:// doi.org/10.1016/j.fuel.2016.08.072

Khasanova N.M., Gabdrakhmanov D.T., Kayukova G.P., Morozov V.P., Mikhaylova A.N. (2017). EPR study of hydrocarbon generation potential of organic-rich domanik rocks, Magnetic Resonance in Solids, 19(1), 17102.

Mamin G., Gafurov M., Yusupov R., Gracheva I., Ganeeva Y., Yusupova T., Orlinskii S.B. (2016). Toward the Asphaltene Structure by Electron Paramagnetic Resonance Relaxation Studies at High Fields (3.4 T). Energy & Fuels, 30(9), pp. 6942-6946. https://doi.org/10.1021/acs.energyfuels.6b00983 Martyanov O.N., Larichev Y.V., Morozov E.V., Trukhan S.N., Kazarian S.G. (2017). The stability and evolution of oil systems studied via advanced methods in situ. Russ. Chem. Rev, 86, pp. 999-1023. https://doi.org/10.1070/ RCR4742

Mehrabi-Kalajahi S.S, Varfolomeev M.A, Yuan C. et al. (2018). EPR as a complementary tool for the analysis of low-temperature oxidation reactions of crude oils. Journal of Petroleum Science and Engineering, 169, pp. 673-682. https://doi.org/10.1016/j.petrol.2018.05.049

Mukhamatdinov I., Gafurov M., Kemalov A., et al. (2018). Study of the oxidized and non-oxidized bitumen modified with additive «Adgezolin» by using electron paramagnetic resonance. IOP Conf. Ser.: Earth Environ. Sci., 155, 012004. https://doi.org/10.1088/1755-1315/155/1/012004

Mukhamatdinov I.I., Salih I.Sh.S., Rakhmatullin I.Z., Sitnov S.A., Laikov A.V., Klochkov V.V., Vakhin A.V. (2020). Influence of Co-based catalyst on subfractional composition of heavy oil asphaltenes during aquathermolysis, Journal of Petroleum Science and Engineering, 186, 106721. https://doi. org/10.1016/j.petrol.2019.106721

Mullins O., Pomerantz A.E., Zuo J., Dong C., Annu J. (2014). Downhole fluid analysis and asphaltene science for petroleum reservoir evaluation. Rev. Chem. Biomol. Eng., 5, pp. 325-345. https://doi.org/10.1146/ annurev-chembioeng-060713-035923

O'Reilly, D. (1958). Paramagnetic Resonance of Vanadyl Etioporphyrin I. J. Chem. Phys., 29(5), pp. 1188-1189. https://doi.org/10.1063/L1744684

Piccinato M., Guedes C., Di Mauro E. (2012). Petroleum Asphaltenes. Crude Oil Emulsions — Composition Stability and Characterization. Ed. M. Abdul-Raouf. Rjeka: InTech, pp. 147-168.

Qin P, Warncke K. (2015a). Electron Paramagnetic Resonance Investigations of Biological Systems by Using Spin Labels, Spin Probes,

НЮЧНО-ТЕХНИЧЕСЩЙ ЖУРНАЛ

www.geors.ru ГЕйРЕСУРСЫ

and Intrinsic Metal Ions. Part A. Methods in Enzymology, 563, pp. 2-684.

Qin P., Warncke K. (2015b). Electron Paramagnetic Resonance Investigations of Biological Systems by Using Spin Labels, Spin Probes, and Intrinsic Metal Ions. Part B. Methods in Enzymology, 564, pp. 2-613.

Raghunathan P. (1991). Evidence for fractal dimension in asphaltene polymers from electron-spin-relaxation measurements. Chem. Phys. Lett. 182, pp. 331-335. https://doi.org/10.1016/0009-2614(91)80224-L

Ramachandran V., van Tol, J., McKenna A., Rodgers R., Marshall A., Dalal N. (2015). High Field Electron Paramagnetic Resonance Characterization of Electronic and Structural Environments for Paramagnetic Metal Ions and Organic Free Radicals in Deepwater Horizon Oil Spill Tar Balls. Anal. Chem., 87(4), pp. 2306-2313. https://doi.org/10.1021/ac504080g Reijerse E.J., Tyryshkin A.M., Dikanov S.A. (1998). Complete determination of nitrogen quadrupole and hyperfine tensors in an oxovanadium complex by simultaneous fitting of multifrequency ESEEM powder spectra. Journal of Magnetic Resonance, 131(2), pp. 295-309. https:// doi.org/10.1006/jmre.1997.1339

Sapunov V.A., Denisov A.Y., Saveliev D.V. et al. (2016). New vector/ scalar Overhauser DNP magnetometers POS-4 for magnetic observatories and directional oil drilling support. Magnetic Resonance in Solids, 18(2), 16209.

Sapunov V.A., Kashin I.V., Ushakov V.A. et al. (2019). Little-known aspects of Overhauser DNP at zero and low magnetic fields stimulated by parallel electron pumping of nitroxide radicals solutions. AIP Conf. Proc., 2174, 020112. https://doi.org/10.1063/L5134263

Schweiger A., Jeschke G. (2001). Principles of Pulse Electron Paramagnetic Resonance. Oxford: OUP.

Stoll S., Schweiger A. (2006). EasySpin, a comprehensive software package for spectral simulation and analysis in EPR. ournal of Magnetic Resonance, 178, pp. 42-55. https://doi.org/10.1016/j.jmr.2005.08.013

Tayeb Ben K., Delpoux O., Barbier J., Marques J., Verstraete J., Vezin H. (2015). Applications of Pulsed Electron Paramagnetic Resonance Spectroscopy to the Identification of Vanadyl Complexes in Asphaltene Molecules. Part 1: Influence of the Origin of the Feed. Energy & Fuels, 29(7), pp. 4608-4615. https://doi.org/10.1021/acs.energyfuels.5b00733

Trukhan S., Yudanov V., Gabrienko A., Subramani V., Kazarian S., Martyanov O. (2014). In Situ Electron Spin Resonance Study of Molecular Dynamics of Asphaltenes at Elevated Temperature and Pressure. Energy & Fuels, 28(10), pp. 6315-6321. https://doi.org/10.1021/ef5015549

Vakhin A.V., Aliev F.A., Mukhamatdinov I.I., Sitnov S.A., Sharifullin A.V., Kudryashov S.I., Afanasiev I.S., Petrashov O.V., Nurgaliev D.K. (2020). Catalytic Aquathermolysis of Boca de Jaruco Heavy Oil with Nickel-Based Oil-Soluble Catalyst. Processes, 8(5), 532. https://doi.org/10.3390/pr8050532 Volodin M.A., Mamin G.V., Izotov V.V., & Orlinskii S.B. (2013). High-frequency EPR study of crude oils. J. Phys.: Conf. Ser., 478, 012003. https:// doi.org/10.1088/1742-6596/478/1/012003

Wang W., Ma Y., Li S., Shi J., Teng J. (2016). Effect of Temperature on the EPR Properties of Oil Shale Pyrolysates. Energy & Fuels, 30(2), pp. 830-834. https://doi.org/10.1021/acs.energyfuels.5b02211

Yakubov M.R., Milordov D.V., Yakubova S.G., Morozov V.I. (2017). Vanadium and paramagnetic vanadyl complexes content in asphaltenes of heavy oils of various productive sediments. Petroleum Science and Technology, 35(14), pp. 1468-1472. https://doi.org/10.1080/10916466.2017.1344708

Yakubova S.G., Abilova G.R., Tazeeva, E.G., Borisova Y.Y., Milordov D.V., Mironov N.A., Yakubov M.R. (2019). Distribution of Vanadium and Nickel in the Case of Two-Step Solvent Fractionation of Asphaltenes of Heavy Oils. Petroleum Chemistry, 59(1), pp. S30-S36. https://doi.org/10.1134/ S0965544119130140

Yen T., Chilingarian G. (1994). Asphaltenes and asphalts. 1. Developments in petroleum science. 40A. New York: Elsevier. https://doi. org/10.1016/S0376-7361(09)70248-1

Yen T., Chilingarian G. (2000). Asphaltenes and asphalts, 2. Developments in petroleum science. 40 B. New York: Elsevier.

Yen T.F., Erdman J.G., Saraceno A.J. (1962). Investigation of the Nature of Free Radicals in Petroleum Asphaltenes and Related Substances by Electron Spin Resonance. Analytical Chemistry, 34(6), pp. 694-700. https:// doi.org/10.1021/ac60186a034

Zavoisky E. (1945). Paramagnetic Relaxation of Liquid Solutions for Perpendicular Fields. Journal of Physics (Academy of Sciences of the USSR), 9(3), pp. 211-216.

Zhang Y., Siskin M., Gray M.R., Walters C.C., Rodgers R.P. (2020). Mechanisms of Asphaltene Aggregation: Puzzles and a New Hypothesis. Energy & Fuels, 34(8), pp. 9094-9107. https://doi.org/10.1021/acs. energyfuels.0c01564

Zhao X., Xu, C., Shi Q. (2015). Porphyrins in Heavy Petroleums: A Review. In: Xu C., Shi Q. (ed.) Structure and Modeling of Complex Petroleum Mixtures. Structure and Bonding, 168. Springer, Cham. https:// doi.org/10.1007/430_2015_189

сведения об авторах

Марат Ревгерович Гафуров - доктор физ.-мат. наук, ведущий научный сотрудник, заместитель директора по научной деятельности, Институт физики, Казанский федеральный университет Россия, 420008, Казань, ул. Кремлевская, д. 16а

Андрей Александрович Пономарев - аспирант, Тюменский индустриальный университет

Россия, 625000, Тюмень, ул. Володарского, д. 38

Георгий Владимирович Мамин - канд. физ.-мат. наук, доцент кафедры квантовой электроники и радиоспектроскопии, Институт физики, Казанский федеральный университет Россия, 420008, Казань, ул. Кремлевская, д. 18

Александр Александрович Родионов - электроник 1 категории кафедры квантовой электроники и радиоспектроскопии, Институт физики, Казанский федеральный университет Россия, 420008, Казань, ул. Кремлевская, д. 18

Фадис Фанилович Мурзаханов - аспирант, Институт физики, Казанский федеральный университет

Россия, 420008, Казань, ул. Кремлевская, д. 18

Таджик Араш - аспирант, Институт геологии и нефтегазовых технологий, Казанский федеральный университет Россия, 420008, Казань, ул. Кремлевская, д. 4

Сергей Борисович Орлинский - канд. физ.-мат. наук, доцент кафедры квантовой электроники и радиоспектроскопии, Институт физики, Казанский федеральный университет Россия, 420008, Казань, ул. Кремлевская, д. 18

Статья поступила в редакцию 09.06.2020;

Принята к публикации 03.11.2020; Опубликована 11.12.2020

Application of pulsed and high-frequency electron paramagnetic resonance techniques to study petroleum disperse systems

M.R. Gafurov1*, A.A. Ponomarev2, G.V. Mamin1, A.A. Rodionov1, F.F. Murzakhanov1, T. Arash1, S.B. Orlinskii1

1Kazan Federal University, Kazan, Russian Federation 2Tyumen Industrial University, Tyumen, Russian Federation Corresponding author: Marat R. Gafurov, e-mail: marat.gafurov@kpfu.ru

Abstract. The spectral and relaxation characteristics of "free" organic radicals (FR) and vanadyl-porphyrin (VP) complexes in various petroleum disperse systems (PDS) like bitumen, petroleum, their high-molecular components and solutions were studied using

stationary (conventional) and pulsed electron paramagnetic resonance (EPR) techniques in two frequency ranges (X- and W-bands, with the microwave radiation frequencies of about 9 GHz and 95 GHz, respectively). The features of the pulsed approaches (electron spin

echo, modulation of the electron spin echo signal decay, electronic relaxation times) and high-frequency EPR for PDS investigations were examined. W-band EPR allows to resolve spectrally the lines from the different paramagnetic centers and more accurately determine their spectral characteristics. It is shown that the electron spin echo can be observed at room temperatures even at high magnetic fields of 3.4 T demonstrating the potential of application of pulsed EPR techniques for the low-cost oilfield measurements. Analysis of the VP transverse magnetization decay curve permits to identify electron-nuclear interactions with the 14N and 1H nuclei in situ while in the EPR spectra these hyperfine interactions usually cannot be detected. It is found from the W-band EPR measurements that FR lineshape cannot be fitted with isotropic parameters in contrast to the established X-band results. The observed effect of increasing the rates of electronic transverse relaxation in asphaltenes is described in the framework of a model of spectral diffusion between the fast- and slow-relaxing paramagnetic centers in supramolecular complexes of asphaltenes.

Keywords: electron paramagnetic resonance, petroleum disperse systems, asphaltenes

Recommended citation: Gafurov M.R., Ponomarev A.A., Mamin G.V., Rodionov A.A., Murzakhanov F.F., Arash T., Orlinskii S.B. (2020). Application of pulsed and high-frequency electron paramagnetic resonance techniques to study petroleum disperse systems. Georesursy = Georesources, 22(4), pp. 2-14. DOI: https:// doi.org/10.18599/grs.2020.4.2-14

Acknowledgments

The authors are thankful to D.T. Sitdikov, M.A. Volodin, A.V. Vakhin, Yu.M. Ganeeva and T.N. Yusupova for their help in carrying out measurements and preparing samples, T.B. Biktagirov for a number of theoretical calculations, V.I. Morozov for discussing the results obtained, A. V. Dooglav for careful reading of the manuscript, valuable comments and for help with the text translation.

The authors thank the reviewers for the opportunity to correct and improve the manuscript.

The work was supported by the RSF grant 19-12-00332.

References

Abragam A., Bleaney B. (1970). Electron Paramagnetic Resonance of Transition Ions. Oxford: Clarendon press.

Abyzgildin Yu.N., Mikhailyuk Yu.N., Yarullin K.S., Rostovskaya A.A. (1977). Porphyrins and metal porphyrin complexes of oils. Moscow: Nauka, 88 p. (In Russ.)

Acevedo S., Guzman K., Ocanto O. (2010). Determination of the number average molecular mass of asphaltenes (Mn) using their soluble A2 fraction and the vapor pressure osmometry (VPO) technique. Energy & Fuels, 24(3), pp. 1809-1812. https://doi.org/10.1021/ef9012714

Al-Muntaser A.A., Varfolomeev M.A., Suwaid M.A. et al. (2020). Hydrothermal upgrading of heavy oil in the presence of water at sub-critical, near-critical and supercritical conditions, Journal of Petroleum Science and Engineering, 184, 106592. https://doi.org/10.1016/j.petrol.2019.106592

Alexandrov A.S., Ivanov A.A., Archipov R.V., Gafurov M.R., Tagirov M.S. (2019). Pulsed NMR spectrometer with dynamic nuclear polarization for weak magnetic fields. Magnetic Resonance in Solids, 21(2), pp. 19203 (1-6). https://doi.org/10.26907/mrsej-19203

Cui Q., Ma X., Nakano K., Nakabayashi K., Miyawaki J., Al-Mutairi A. et al. (2018). Hydrotreating reactivities of atmospheric residues and correlation with their composition and properties. Energy & Fuels, 32(6), pp. 6726-6736. https://doi.org/10.1021/acs.energyfuels.8b01150

Cui Q., Nakabayashi K., Ma X., Ideta K., Miyawaki J., Marafi A.M. et al. (2017). Examining the molecular entanglement between V=O complexes and their matrices in atmospheric residues by ESR. RSC advances, 7(60), pp. 37908-37914. https://doi.org/10.1039/C7RA06436E

Davydov V.V., Dudkin V.I., Myazin N.S., Rud' V.Yu. (2018). Peculiarity of the Nuclear Magnetic Resonance Method Application for the Liquid Medium Flow Parameters Control. Applied Magnetic Resonance, 49(7), pp. 665-679. https://doi.org/10.1007/s00723-018-0994-1

Deligiannakis Y., Louloudi M., & Hadjiliadis N. (2000). Electron spin echo envelope modulation (ESEEM) spectroscopy as a tool to investigate the coordination environment of metal centers. Coordination Chemistry

Reviews, 204(1), pp. 1-112. https://doi.org/10.1016/S0010-8545(99)00218-0 Di Mauro E., Guedes, C.L.B., Nascimento O.R. (2005). Multifrequency (X-band to W-band) CW EPR of the organic free radical in petroleum asphaltene. Applied Magnetic Resonance, 29(4), pp. 569-575. https://doi. org/10.1007/BF03166333

Dickson F.E., Kunesh, C.J., McGinnis E.L., Petrakis L. (1972). Use of electron spin resonance to characterize the vanadium (IV)-sulfur species in petroleum. Anal. Chem., 44(6), pp. 978-981. https://doi.org/10.1021/ ac60314a009

Dikanov S.A., Tsvetkov Y.D. (1992). Electron Spin-Echo Envelope Modulation (ESEEM) Spectroscopy. USA: CRC Press, 432 p.

Dolomatov M., Gafurov M., Rodionov A., Mamin G., González L.M., Vakhin A., Petrov A., Bakhtizin R., Khairudinov I., Orlinskii, S. (2018). Low-temperature thermal decomposition of heavy petroleum distillates: interconnection between the electrical properties and concentration of paramagnetic centres. IOP Conf. Ser.: Earth Environ. Sci., 155, 012007. https://doi.org/10.1088/1755-1315/155/1Z012007

Dolomatov M.U., Rodionov A.A., Gafurov M.R., Petrov A.V., Biktagirov T.B., Bakhtizin R.Z., Makarchikov S.O., Khairudinov I.Z. and Orlinskii S.B. (2016). Concentration of paramagnetic centres at low-temperature thermal destruction of asphaltenes of heavy petroleum distillates. Magnetic Resonance in Solids, 18, 16101. http://mrsej.kpfu.ru/contents.html#16101

Dzuba S.A. (2013). Structural studies of biological membranes using ESEEM spectroscopy of spin labels and deuterium substitution. Journal of Structural Chemistry, 54(S1), pp. 5-18. (In Russ.)

Eaton, G.R., Eaton, S.S., Barr, D.P., Weber, R.T. (2010). Quantitative EPR. Vienna: Springer-Verlag Wien. https://doi.org/10.1007/978-3-211-92948-3 Gafurov M.R., Volodin M.A., Rodionov et al. (2018). EPR study of spectra transformations of the intrinsic vanadyl-porphyrin complexes in heavy crude oils with temperature to probe the asphaltenes' aggregation. Journal of Petroleum Science and Engineering, 166, pp. 363-368. https:// doi.org/10.1016/j.petrol.2018.02.045

Galukhin A., Bolmatenkov D., Osin Y. (2018). Heavy oil oxidation in the nano-porous medium of synthetic opal, RSC Adv., 8, pp. 18110-18116. https://doi.org/10.1039/C8RA02822B

Garif'yanov N.S., Kozyrev B.M. (1956). Paramagnetic resonance in anthracite and other carbon-containing substances. Zh. Eksp. Teor. Fiz., 30(2), pp. 272-276. (In Russ.)

Gilinskaya L.G. (2008). EPR spectra of complexes V (IV) and the structure of petroleum porphyrins. Journal of Structural Chemistry, 49(2), pp. 245-254. (In Russ.)

Gizatullin B., Gafurov M., Vakhin et al. (2019). Native Vanadyl Complexes in Crude Oil as Polarizing Agents for In Situ Proton Dynamic Nuclear Polarization. Energy & Fuels, 33(11), pp. 10923-10932. https://doi. org/10.1021/acs.energyfuels.9b03049

Gizatullin B., Gafurov, M., Rodionov A., Mamin, G., Mattea, C., Stapf, S., Orlinskii, S. (2018). Proton-Radical Interaction in Crude Oil - A Combined NMR and EPR Study. Energy & fuels, 32(11), pp. 11261-11268. https://doi. org/10.1021/acs.energyfuels.8b02507

Gracheva I, Gafurov M., Mamin G., Biktagirov T., Rodionov A., Galukhin V., Orlinskii S.B. (2016). ENDOR Study of Nitrogen Hyperfine and Quadropole Tensors in Vanadyl Porhyrins of Heavy Crude Oils. Magnetic Resonance in Solids, 18, 16102. http://mrsej.kpfu.ru/contents.html#16102

Gutowsky H., Roger Ray B., Rutledge R., Unterberger R. (1958). Carbonaceous Free Radicals in Crude Petroleum. J. Chem. Phys., 28, pp. 744-745. https://doi.org/10.1063/L1744250

Il'yasov A.V. (1962). Determination of vanadium content in oils and petroleum products by EPR method. Chemistry and Technology of Fuels and Oils, 9, pp. 63-67. (In Russ.)

Il'yasov A.V., Garif'yanov N.S., Ryzhmanov Yu.S. (1961). Electron paramagnetic resonance in some types of natural oils and their heavy fractions. Chemistry and Technology of Fuels and Oils, 1, pp. 28-31. (In Russ.)

Ilyin S.O., Arinina M.P., Polyakova M.Y. et al. (2016). Rheological comparison of light and heavy crude oils. Fuel, 186, pp. 157-167. https:// doi.org/10.1016/j.fuel.2016.08.072

Khasanova N.M., Gabdrakhmanov D.T., Kayukova G.P., Morozov V.P., Mikhaylova A.N. (2017). EPR study of hydrocarbon generation potential of organic-rich domanik rocks, Magnetic Resonance in Solids, 19(1), 17102.

Mamin G., Gafurov M., Yusupov R., Gracheva I., Ganeeva Y., Yusupova T., Orlinskii S.B. (2016). Toward the Asphaltene Structure by Electron Paramagnetic Resonance Relaxation Studies at High Fields (3.4 T). Energy & Fuels, 30(9), pp. 6942-6946. https://doi.org/10.1021/acs.energyfuels.6b00983 Martyanov O.N., Larichev Y.V., Morozov E.V., Trukhan S.N., Kazarian S.G. (2017). The stability and evolution of oil systems studied via advanced methods in situ. Russ. Chem. Rev, 86, pp. 999-1023. https://doi.org/10.1070/ RCR4742

НАУЧНО-ТЕХНИЧЕСКИЙ ЖУРНАЛ ^HB

ÉdiliMÜili™ 1

reOPECyPCbl / GEORESOURCES

2020. T. 22. № 4. C. 2-14

Mehrabi-Kalajahi S.S, Varfolomeev M.A, Yuan C. et al. (2018). EPR as a complementary tool for the analysis of low-temperature oxidation reactions of crude oils. Journal of Petroleum Science and Engineering, 169, pp. 673-682. https://doi.org/10.1016/j.petrol.2018.05.049

Mukhamatdinov I., Gafurov M., Kemalov A., et al. (2018). Study ofthe oxidized and non-oxidized bitumen modified with additive «Adgezolin» by using electron paramagnetic resonance. IOP Conf. Ser.: Earth Environ. Sci., 155, 012004. https://doi.org/10.1088/1755-1315/155/1Z012004

Mukhamatdinov I.I., Salih I.Sh.S., Rakhmatullin I.Z., Sitnov S.A., Laikov A.V., Klochkov V.V., Vakhin A.V. (2020). Influence of Co-based catalyst on subfractional composition of heavy oil asphaltenes during aquathermolysis, Journal of Petroleum Science and Engineering, 186, 106721. https://doi. org/10.1016/j.petrol.2019.106721

Mullins O., Pomerantz A.E., Zuo J., Dong C., Annu J. (2014). Downhole fluid analysis and asphaltene science for petroleum reservoir evaluation. Rev. Chem. Biomol. Eng., 5, pp. 325-345. https://doi.org/10.1146/ annurev-chembioeng-060713-035923

Murav'ev F.A., Vinokurov V.M., Galeev A.A., Bulka G.R., Nizamutdinov N.M., Khasanova N.M. (2006). Paramagnetism and nature of dispersed organic matter in the Permian deposits of Tatarstan. Georesursy = Georesources, 2(19), pp. 40-45. (In Russ.)

Nesterov I.I., Alexandrov V.M., Ponomarev A.A., Zavatsky M.D., Lobodenko E.I., Kobylinskiy D.A., Kadyrov M.A. (2019). Experimental studies of radical reactions of hydrocarbons conversion. Oil and Gas Studies, 4, pp. 57-69. (In Russ.) https://doi.org/10.31660/0445-0108-2019-4-57-69 O'Reilly, D. (1958). Paramagnetic Resonance of Vanadyl Etioporphyrin I. J. Chem. Phys, 29(5), pp. 1188-1189. https://doi.org/10.1063/L1744684

Piccinato M., Guedes C., Di Mauro E. (2012). Petroleum Asphaltenes. Crude Oil Emulsions — Composition Stability and Characterization. Ed. M. Abdul-Raouf. Rjeka: InTech, pp. 147-168.

Ponomarev A.A. (2019). The mechanism of cracking hydrocarbons in electromagnetic fields - to the question of oil in Bazhenov formation. Oil and Gas Studies, 1, pp. 14-18. (In Russ.) https://doi. org/10.31660/0445-0108-2019-1-14-18

Qin P, Warncke K. (2015a). Electron Paramagnetic Resonance Investigations of Biological Systems by Using Spin Labels, Spin Probes, and Intrinsic Metal Ions. Part A. Methods in Enzymology, 563, pp. 2-684.

Qin P., Warncke K. (2015b). Electron Paramagnetic Resonance Investigations of Biological Systems by Using Spin Labels, Spin Probes, and Intrinsic Metal Ions. Part B. Methods in Enzymology, 564, pp. 2-613.

Raghunathan P. (1991). Evidence for fractal dimension in asphaltene polymers from electron-spin-relaxation measurements. Chem. Phys. Lett. 182, pp. 331-335. https://doi.org/10.1016/0009-2614(91)80224-L

Ramachandran V., van Tol, J., McKenna A., Rodgers R., Marshall A., Dalal N. (2015). High Field Electron Paramagnetic Resonance Characterization of Electronic and Structural Environments for Paramagnetic Metal Ions and Organic Free Radicals in Deepwater Horizon Oil Spill Tar Balls. Anal. Chem, 87(4), pp. 2306-2313. https://doi.org/10.1021/ac504080g Reijerse E.J., Tyryshkin A.M., Dikanov S.A. (1998). Complete determination of nitrogen quadrupole and hyperfine tensors in an oxovanadium complex by simultaneous fitting of multifrequency ESEEM powder spectra. Journal of Magnetic Resonance, 131(2), pp. 295-309. https:// doi.org/10.1006/jmre.1997.1339

Safieva R.Z. (2004). Chemistry of oil and gas. Petroleum Dispersed Systems: Composition and Properties. Part 1. Moscow: Gubkin Russian State University of Oil and Gas, 112 p. (In Russ.)

Sapunov V.A., Denisov A.Y., Saveliev D.V. et al. (2016). New vector/ scalar Overhauser DNP magnetometers POS-4 for magnetic observatories and directional oil drilling support. Magnetic Resonance in Solids, 18(2), 16209.

Sapunov VA., Kashin I.V., Ushakov V.A. et al. (2019). Little-known aspects of Overhauser DNP at zero and low magnetic fields stimulated by parallel electron pumping of nitroxide radicals solutions. AIP Conf. Proc., 2174, 020112. https://doi.org/10.1063/L5134263

Schweiger A., Jeschke G. (2001). Principles of Pulse Electron Paramagnetic Resonance. Oxford: OUP.

Stoll S., Schweiger A. (2006). EasySpin, a comprehensive software package for spectral simulation and analysis in EPR. ournal of Magnetic Resonance, 178, pp. 42-55. https://doi.org/10.1016/j.jmr.2005.08.013

Syunyaev Z.I. (1980). Concentration of complex structural units in petroleum dispersed systems and methods of its regulation. Chemistry and Technology of Fuels and Oils, 16(7), pp. 484-489. (In Russ.)

Syunyaev Z.I., Safieva R.Z., Syunyaev R.Z. (1990). Petroleum dispersed systems. Moscow: Chemistry, 224 p. (In Russ.)

Tayeb Ben K., Delpoux O., Barbier J., Marques J., Verstraete J., Vezin H. (2015). Applications of Pulsed Electron Paramagnetic Resonance Spectroscopy to the Identification of Vanadyl Complexes in Asphaltene Molecules. Part 1: Influence of the Origin of the Feed. Energy & Fuels, 29(7), pp. 4608-4615. https://doi.org/10.1021/acs.energyfuels.5b00733

Trukhan S., Yudanov V., Gabrienko A., Subramani V., Kazarian S.,

Martyanov O. (2014). In Situ Electron Spin Resonance Study of Molecular Dynamics of Asphaltenes at Elevated Temperature and Pressure. Energy & Fuels, 28(10), pp. 6315-6321. https://doi.org/10.1021/ef5015549

Vakhin A.V., Aliev F.A., Mukhamatdinov I.I., Sitnov S.A., Sharifullin A.V., Kudryashov S.I., Afanasiev I.S., Petrashov O.V., Nurgaliev D.K. (2020). Catalytic Aquathermolysis of Boca de Jaruco Heavy Oil with Nickel-Based Oil-Soluble Catalyst. Processes, 8(5), 532. https://doi.org/10.3390/pr8050532 Volodin M.A., Mamin G.V., Izotov V.V., & Orlinskii S.B. (2013). High-frequency EPR study of crude oils. J. Phys.: Conf Ser., 478, 012003. https:// doi.org/10.1088/1742-6596/478/1/012003

Wang W., Ma Y., Li S., Shi J., Teng J. (2016). Effect of Temperature on the EPR Properties of Oil Shale Pyrolysates. Energy & Fuels, 30(2), pp. 830-834. https://doi.org/10.1021/acs.energyfuels.5b02211

Yakubov M.R., Milordov D.V., Yakubova S.G., Morozov V.I. (2017). Vanadium and paramagnetic vanadyl complexes content in asphaltenes of heavy oils of various productive sediments. Petroleum Science and Technology, 35(14), pp. 1468-1472. https://doi.org/10.1080/10916466.2017.1344708

Yakubova S.G., Abilova G.R., Tazeeva, E.G., Borisova Y.Y., Milordov D.V., Mironov N.A., Yakubov M.R. (2019). Distribution of Vanadium and Nickel in the Case of Two-Step Solvent Fractionation of Asphaltenes of Heavy Oils. Petroleum Chemistry, 59(1), pp. S30-S36. https://doi.org/10.1134/ S0965544119130140

Yen T., Chilingarian G. (1994). Asphaltenes and asphalts. 1. Developments in petroleum science. 40A. New York: Elsevier. https://doi. org/10.1016/S0376-7361(09)70248-1

Yen T., Chilingarian G. (2000). Asphaltenes and asphalts, 2. Developments in petroleum science. 40 B. New York: Elsevier.

Yen T.F., Erdman J.G., Saraceno A.J. (1962). Investigation of the Nature of Free Radicals in Petroleum Asphaltenes and Related Substances by Electron Spin Resonance. Analytical Chemistry, 34(6), pp. 694-700. https:// doi.org/10.1021/ac60186a034

Zavoisky E. (1945). Paramagnetic Relaxation of Liquid Solutions for Perpendicular Fields. Journal of Physics (Academy of Sciences ofthe USSR), 9(3), pp. 211-216.

Zhang Y., Siskin M., Gray M.R., Walters C.C., Rodgers R.P. (2020). Mechanisms of Asphaltene Aggregation: Puzzles and a New Hypothesis. Energy & Fuels, 34(8), pp. 9094-9107. https://doi.org/10.1021/acs. energyfuels.0c01564

Zhao X., Xu, C., Shi Q. (2015). Porphyrins in Heavy Petroleums: A Review. In: Xu C., Shi Q. (ed.) Structure and Modeling of Complex Petroleum Mixtures. Structure and Bonding, 168. Springer, Cham. https:// doi.org/10.1007/430_2015_189

About the Authors

MaratR. Gafurov - Dr. Sci. (Physics and Mathematics), Leading Researcher, Deputy Director for Research, Institute of Physics, Kazan Federal University

16a Kremlevskaya st., Kazan, 420008, Russian Federation

AndreyA. Ponomarev - Post-graduate student, Tyumen Industrial University

38 Volodarsky st., Tyumen, 625000, Russian Federation

Georgiy V. Mamin - Cand. Sci. (Physics and Mathematics), Associate Professor, Department of Quantum Electronics and Radiospectroscopy, Institute of Physics, Kazan Federal University 18 Kremlevskaya st., Kazan, 420008, Russian Federation

Alexander A. Rodionov - Electronic engineer, Department of Quantum Electronics and Radiospectroscopy, Institute of Physics, Kazan Federal University

18 Kremlevskaya st., Kazan, 420008, Russian Federation

Fadis F. Murzakhanov - Post-graduate student, Institute of Physics, Kazan Federal University

18 Kremlevskaya st., Kazan, 420008, Russian Federation

Tajik Arash - Postgraduate Student, Institute of Geology and Petroleum Technologies, Kazan Federal University

4 Kremlevskaya st., Kazan, 420008, Russian Federation

Sergey B. Orlinskii - Cand. Sci. (Physics and Mathematics), Associate Professor, Department of Quantum Electronics and Radiospectroscopy, Institute of Physics, Kazan Federal University 18 Kremlevskaya st., Kazan, 420008, Russian Federation

Manuscript received 9 June 2020;

Accepted 3 November 2020; Published 11 December 2020