Применение методов ЯМР к исследованию геометрических характеристик порового пространства гранулярных силикатов и свойств поровых флюидов

Перепухов, Александр Максимович

Применение методов ЯМР к исследованию геометрических характеристик порового пространства гранулярных силикатов и свойств поровых флюидов тема автореферата и диссертации по физике, 01.04.17 ВАК РФ

Перепухов, Александр Максимович АВТОР

кандидата физико-математических наук УЧЕНАЯ СТЕПЕНЬ

Долгопрудный МЕСТО ЗАЩИТЫ

2013 ГОД ЗАЩИТЫ

01.04.17 КОД ВАК РФ

Диссертация по физике на тему «Применение методов ЯМР к исследованию геометрических характеристик порового пространства гранулярных силикатов и свойств поровых флюидов»

Автореферат диссертации на тему "Применение методов ЯМР к исследованию геометрических характеристик порового пространства гранулярных силикатов и свойств поровых флюидов"

На правах рукописи

Перепухов Александр Максимович

ПРИМЕНЕНИЕ МЕТОДОВ ЯМР К ИССЛЕДОВАНИЮ ГЕОМЕТРИЧЕСКИХ ХАРАКТЕРИСТИК ПОРОВОГО ПРОСТРАНСТВА ГРАНУЛЯРНЫХ СИЛИКАТОВ И СВОЙСТВ ПОРОВЫХ ФЛЮИДОВ

01.04.17 - химическая физика, горение и взрыв, физика экстремальных

состояний вещества

АВТОРЕФЕРАТ диссертации на соискание ученой степени кандидата физико-математических наук

2 8 ноя т

Долгопрудный - 2013

005540233

Работа выполнена в Московском физико-техническом институте (государственном университете)

Научный руководитель:

Официальные оппоненты:

Максимычев Александр

Витальевич

доктор физико-математических наук, профессор

Иванов Виктор Владимирович

доктор физико-математических наук, профессор, член-корреспондент РАН,

генеральный директор ООО «Метрологический центр

РОСНАНО»

Набиев Шавкат Шарифович

доктор физико-математических наук, старший научный сотрудник, начальник отдела физико-химического анализа ЦФХТ НЩ «Курчатовсий институт»

Ведущая организация: Институт органической химии

им. Н.Д. Зелинского Российской академии наук (ИОХ РАН)

Защита состоится «24» декабря 2013 года в 11-00 час. на заседании диссертационного совета Д 520.009.05 на базе НИЦ «Курчатовский институт» по адресу: 123182, г. Москва, пл. Академика Курчатова, д. 1.

С диссертацией можно ознакомиться в библиотеке НИЦ «Курчатовский институт».

Автореферат разослан «19 » ноября 2013 г.

Ученый секретарь диссертационного совета,

доктор химических наук СерикВ.Ф.

ОБЩАЯ ХАРАКТЕРИСТИКА РАБОТЫ Актуальность темы исследований

Метод ядерного магнитного резонанса (ЯМР) - это неразрушающий метод анализа, который позволяет исследовать пористые среды (например, нефтяные коллекторы) в их естественном состоянии посредством изучения поведения флюида, заполняющего поры. С другой стороны, методы ЯМР позволяют изучать флюиды, заполняющие поры, без необходимости извлечения их из порового пространства. За последние 20 лет методы изучения пористых сред с помощью ЯМР получили большое развитие. ЯМР анализаторы, помещаемые внутрь нефтяной скважины, активно используются в нефтеразведке.

Среди различных методик аттестации пористых сред с помощью ЯМР можно выделить релаксационные измерения, изучение ограниченной диффузии жидкости внутри порового пространства методом ЯМР с импульсным градиентом магнитного поля, метод DDIF (Dacay due to Diffusion in Internal Fields), ЯМР криопорометрию и твердотельный ЯМР. Первые три методики используются, в основном, для определения распределения пор по размерам в диапазоне от сотен нанометров до сотен микрометров. Метод ЯМР криопорометрии используется для изучения пористых сред с порами от нанометров до сотен нанометров. Методики твердотельного ЯМР позволяют анализировать химический состав флюида в поре.

Метод аттестации пористых сред, основанный на измерениях времён продольной (Г,) и поперечной (Т2) ЯМР релаксации ядер 'Н флюида, заполняющего поровое пространство, является одним из самых распространённых методов ЯМР, применяемых для изучения пористых сред. В релаксационном методе распределение пор по размерам получают из спектра времён релаксации (Т\ или Т2), полагая, что скорость релаксации пропорциональна отношению площади поверхности поры к её объёму (S/V):

15 / 1 \ '1,2 '

Коэффициенты пропорциональности рх и рг в соотношении (1) обычно рассматриваются как феноменологические параметры для данной пары твёрдая стенка/флюид. Большинство работ по релаксационным методам проводились или на реальных кернах (диапазон размеров пор в которых примерно 100 нм — 1 мм), или на модельных объектах — пористых стёклах с порами значительно меньшего размера пор (1-10 нм). При этом имеются существенные расхождения в данных, полученных для пористых стёкол и натуральных песчаников. Таким образом, для понимания механизмов релаксации ядер ]Н флюидов в

натуральных пористых средах необходимы исследования на модельных системах с характерным размером пор 0,1 мкм - 1 мм.

Значительный интерес представляет возможность различения водного и углеводородного порового флюидов по соотношению времён Т\1Т2. Однако, литературные данные по данному поводу весьма противоречивы.

Таким образом, интерпретация релаксационных данных в случае пористых сред сопряжена со сложностями, обусловленными зависимостью скоростей релаксации не только от размеров пор, но и от состава флюида, температуры, смачиваемости флюидом поверхности, наличия парамагнитных примесей, напряжённости постоянного магнитного поля и др. Для выявления основных закономерностей релаксационных процессов в пористых средах необходимо построение физической модели, описывающей зависимости скоростей продольной и поперечной релаксации ядер 'Н воды и углеводородов от указанных параметров. Для модельных исследований силикатных пород целесообразно использовать стеклянные шарики, поскольку диоксид кремния является основным компонентом, как стекла, так и натуральных песчаников, а характерный размер коллекторных пор (0,1 мкм — 100 мкм) легко получить, варьируя размер шариков.

Цель и задачи диссертационной работы

Основной целью работы является выявление основных закономерностей ЯМР релаксации протонов флюидов, заполняющих силикатные пористые среды, путём экспериментального исследования скоростей продольной и поперечной релаксации ядер 'Н флюида в среде стеклянных шариков и построение простой модели, обобщающей полученные данные и пригодной для интерпретации релаксационных данных в случае естественных силикатных пород.

Для достижения поставленной цели ставились следующие

задачи:

• Разработка и анализ применимости метода обратного преобразования Лапласа для обработки сложных релаксационных зависимостей.

• Определение для ядер 'Н зависимости времён продольной и поперечной релаксации ядерной намагниченности от характерного размера пор для широкого набора флюидов в модельной среде стеклянных шариков. Применение релаксационных данных, полученных для модельной среды, к анализу порового пространства реальных пород. Верификация релаксационных

данных о геометрии порового пространства независимым адсорбционным методом.

• Непосредственное исследование парамагнитных свойств поверхности силикатов методом ЭПР и химическими методами.

• Изучение влияния смачиваемости поверхности флюидом на скорости продольной и поперечной протонной релаксации флюида.

• Разработка модели поверхностной релаксации ядерной намагниченности, позволяющей анализировать зависимость скорости релаксации от состава порового флюида, размера пор, резонансной частоты.

Научная новизна

• Получены калибровочные значения коэффициентов релаксации р\ и р2 для воды и ряда углеводородов (гексан, гептан, нонан, декан, бензол) в модельной среде на основе диоксида кремния, позволяющие определять распределение пор по размерам в естественных силикатных породах. Впервые релаксационные данные о геометрии порового пространства верифицированы адсорбционным методом.

• Впервые показано, что поверхностная протонная релаксация флюидов в пористых силикатах не требует присутствия посторонних парамагнитных примесей, поскольку на поверхности силикатов присутствуют парамагнитные группы, содержащие ненасыщенные связи типа (ОзБь, 038ь0-, ОзБ^О-О-)

• В рамках предложенной модели показано, что отношение времён продольной и поперечной поверхностной релаксации Т^П^ флюида в поре определяется частотой Ларморовой прецессии ядер 'Н флюида и временем адсорбции молекул флюида на поверхности силикатов.

• Впервые обоснован алгоритм определения водного и углеводородного компонентов порового флюида на основании сопоставления времен релаксации Т\ и Т2.

Основные положения, выносимые на защиту

• Оценки для разрешающей способности алгоритма обратного преобразования Лапласа в зависимости от длины временного ряда, весовых соотношений кинетических констант и отношения сигнал/шум в исходных данных.

• Значения коэффициентов релаксации рх и р2 ядер 'Н воды и ряда углеводородов (гексан, гептан, нонан, декан и бензол), применимые для анализа распределения пор по размерам в натуральных силикатных породах.

• Обоснование с помощью физических и химических методов исследования наличия на поверхности силикатов собственных парамагнитных центров в структуре БЮ2, обусловливающих парамагнитный механизм релаксации ядер 'Н флюида на поверхности пор.

• Модель парамагнитной релаксации протонной намагниченности флюида на поверхности силикатов, описывающая зависимость отношения времён Г^/Т^ от времени адсорбции молекул флюида на поверхности пор и частоты Ларморовой прецессии ядер 'Н.

Практическая значимость

Метод обратного преобразования Лапласа применим для обработки широкого класса сложных релаксационных и кинетических зависимостей, удовлетворяющих выявленным в работе ограничениям.

Полученные в работе значения коэффициентов продольной и поперечной релаксации р\ и р2 для воды и ряда углеводородов (гексан, гептан, нонан, декан и бензол) могут быть использованы для определения распределения пор по размерам в естественных силикатных породах.

Зависимость отношения времён продольной и поперечной релаксации от времени адсорбции флюида на поверхности пор позволяет использовать отношение времён Т1Х/Т2,$ для выделения водной и углеводородной фаз флюида заполняющего силикатную пористую среду.

Зависимость отношения времён продольной и поперечной релаксации от частоты Ларморовой прецессии ядер 'Н может быть использована для выбора оптимального значения напряжённости постоянного магнитного поля для анализа состава флюида, заполняющего силикатную пористую среду.

Полученная в работе модель позволяет полуколичественно моделировать влияние на скорость релаксации состава флюида, размера пор, температуры и напряженности постоянного магнитного поля.

Апробация работы

Основные результаты диссертационной работы были представлены на научной конференции МФТИ (Долгопрудный, 2007, 2008, 2009, 2011, 2012, 2013); Зимней международной школе-конференции «Магнитный резонанс и его приложения» (Санкт-Петербург, 2009); Международной школе-конференции «Magnetic Resonance And Magnetic Phenomena In Chemical And Biological Physics.» (Новосибирск, 2012); международной конференции «Петрофизическое моделирование осадочных пород» (Петергоф, 2012).

Публикации

По материалам диссертации опубликованы монография, 3 статьи в рецензируемых журналах, входящих в перечень ВАК РФ, и 11 тезисов докладов на российских и международных конференциях.

Структура и объём диссертации

Работа состоит из введения, шести глав, выводов и списка литературы. Работа содержит 145 страниц, 56 рисунков и 8 таблиц. Список литературы включает 114 наименований.

СОДЕРЖАНИЕ РАБОТЫ

Во Введении обоснована актуальность темы работы и выбор объектов и методов исследования, сформулированы цели и задачи работы.

В Первой главе («Методы анализа пористых сред с помощью ЯМР») представлен литературный обзор методов изучения пористых сред с помощью ЯМР. Приведены основы методов ЯМР каротажа, обсуждаются экспериментальные и теоретические результаты работ, посвященных изучению релаксации протонной намагниченности флюидов в модельных и природных пористых средах.

Во Второй главе («Материалы и методы») приводится описание объектов исследования, методов подготовки образцов, а также экспериментальных методов анализа, используемых в работе.

В качестве объектов исследования использованы модельные пористые среды из стеклянных шариков, кварцевый песок (рис. 1), а также заполненный нефтью керн горной породы.

Рис. 1. Микрофотографии стеклянных шариков трёх размерных групп: 57 ± 18 мкм (а), 389 ± 39 мкм (б) и 436 ± 44 мкм (в) и песка (г).

В качестве порового флюида использовались химически чистые вода, бензол (С6Н6) и ряд предельных нормальных углеводородов (УВ): гексан (С6Н]4), гептан (СуН^), нонан (С9Н2о) и декан (С10Н22). Химическая чистота веществ, используемых в качестве порового флюида, подтверждалась методом 13С ЯМР (рис. 2).

1,10

И-1-1-1-1-1-1-'-'-1-1-1-1-'-'-1-1-'-1-'-1 '-'

55 30 25 20 15 ррщ

Рис. 2 Пример спектра 13С ЯМР декана, используемого в качестве порового флюида.

Для релаксационных измерений использовали прибор Bruker the Minispec. Определение времен Г, и Т2 проводили с применением импульсных последовательностей Saturation Recovery и CPMG,

соответственно. В экспериментах использованы два магнита, обеспечивающие рабочие частоты по протонам 10 МГц и 20 МГц для образцов диаметром 40 мм и 10 мм, соответственно. Спектры времён релаксации получали, применяя обратное преобразование Лапласа к кинетическим релаксационным зависимостям. Спектры ЯМР высокого разрешения на ядрах 'Н, 13С, 27Al и 29Si регистрировали на ЯМР-спектрометре Varian Unity Inova 500М WB с рабочей частотой 500 МГц для ядер 'Н.

Изображения исследуемых гранулярных сред получали с помощью растрового электронного микроскопа Quanta 200 (FEI). Элементный состав поверхности порообразующего материала определялся методом электронно-зондового микроанализа на установке Quanta 200 (FEI). Характеристическое излучение атомов на поверхности образцов возбуждалось электронным пучком с энергией 0,2-3 кэВ. Площадь анализируемой поверхности составляла ~1 мкм2.

Регистрацию спектров электронного парамагнитного резонанса (ЭПР) проводили на ЭПР-спектрометре BRUKER ESP300. Измерения велись на частоте v = 9,41 ГГц при комнатной температуре.

Для определения удельной поверхности образцов пористых сред использовали адсорбционный анализатор удельной поверхности Sorbi-MS.

Углы смачивания флюидами поверхности стекла измеряли с помощью специального оптического стола САМ 101 (KSV Instruments).

В Третьей главе («Алгоритм обратного преобразования Лапласа для обработки сложных релаксационных зависимостей»)

обсуждается алгоритм, реализованный для получения спектров времён релаксации из многокомпонентных релаксационных кривых.

При исследовании динамики различных физико-химических процессов часто возникает задача выделения из сложной временной зависимости экспоненциальных релаксационных компонент вида:

М(0 = М0е"1т (2)

где /-время, M - регистрируемая экспериментально величина, а Г-время релаксации системы. При параллельном протекании нескольких релаксационных процессов (в гетерогенных системах или в гомогенных средах при наличии нескольких каналов релаксации), кинетическая кривая представляет собой сумму нескольких экспоненциальных кинетических компонент. В частности, измерения времен релаксации в ЯМР часто приводят к многоэкспоненциальным зависимостям M(t). Кроме того, экспериментальная зависимость M(t)

всегда содержит шумовую составляющую E{jt). В случае аддитивного шума временная зависимость M{t) будет иметь вид:

М{1) = ^е"т'+Е{1) (3)

¿=1

(где g. - относительный вклад / -той компоненты.)

Стандартной процедурой выделения экспоненциальных компонент из сложной релаксационной зависимости (3) является представление экспериментальных данных в полулогарифмических координатах lnM(t) — t и последовательное выделение наиболее медленных компонент. Последовательное повторение этой процедуры позволяет выделить две-три экспоненциальные компоненты, если времена релаксации слагаемых в (3) различаются не менее, чем на порядок величины.

Альтернативным подходом к анализу кинетических зависимостей вида (3) является использование обратного преобразования Лапласа (LПрименение L'1 к экспоненциальной функции (2) преобразует ее в дельта-функцию д(т-Т), в силу линейности преобразования Лапласа применение L'1 к зависимости вида (3) дает спектр времён релаксации Г, исследуемой системы:

Ь'\±8,е-"т') = ±Е,3(т-Т,) (4)

ï=l 1=1

Достоинством данного подхода является то, что выполнение операции L'1 не требует априорного задания числа компонент п.

Для численного анализа данных соотношение (3) удобно переписать в матричном виде:

M = KF + Е (5)

где M(tk) - вектор значений кинетической кривой M(t) размерности и (и- количество дискретных значений tk), К - ядро преобразования, представляющее собой квадратную матрицу, элементы которой выражаются соотношением: K;j=exp{-tk/Tj), к=1,...п, j=l,...n, F(zj) -вектор значений искомого спектра времён релаксации F(t) для п точек тг

Уравнение (4) формально можно переписать в виде F=ICl(M-E). Однако, матрица К - вырождена и для неё не существует обратной матрицы. С применением сингулярного разложения матрицы К, в работе реализован алгоритм обратного преобразования Лапласа, пригодный для получения спектра времён релаксации из экспериментального ряда дискретных значений M(tk). Проведен анализ разрешающей способности предложенного алгоритма на модельных

системах в зависимости от длины временного ряда, от соотношения времён Т-, в спектре релаксации и от уровня шума в анализируемых данных.

На рис. 3 представлена кинетическая зависимость и спектр времён релаксации, полученный в результате применения обратного преобразования Лапласа к модельной кривой. Рис. 3 показывает, что в отсутствие шума преобразование Ь~' надёжно разрешает времена релаксации, отличающиеся на порядок, и точно воспроизводит весовые коэффициенты g,.

0.8 10.6

|0,4 0.2 0

5/15

4/15

rt 1 3/15

с: | 2/15

< f

I : V

0 20 40 60 80 100 0 0.001 0,01 0,1 « Т

10 100

Рис. 3. а) Модельная кинетическая кривая, задаваемая соотношением

М(/) = М0(1/15-е~'/ОЖ)| +2/15-е~"0'01 + 3/15-е-"0,1 + 4/15-е~'" +5/15-е-"10), (б)

спектр времён релаксации, полученный путем применения обратного преобразования Лапласа к данной зависимости.

Влияние длины временного ряда п на вид спектра F(г) исследовано на примера модельной зависимости:

Mit) = М0 ■ (0.5е~'/01 + 0.5e""1 °) (6)

Показано, что спектры с удовлетворительным разрешением могут быть получены даже в случае коротких (и~10) временных рядов. В качестве оптимальной длины временного ряда выбрана величина п~50. Оценки разрешающей способности алгоритма L"1 показали, что минимальное отношение разрешимых времен релаксации в отсутствие шума TJTk ~2.

Влияние шума на разрешающую способность алгоритма L'1 исследовано на модельных кинетических кривых вида (6) с добавлением «белого» шума различной мощности. Обнаружено, что при отношении сигнал/шум (SNR) равном 100 и более функция M(t) воспроизводится точно. При низких значениях SNR-10 смещение

максимумов и искажение соотношения весовых коэффициентов gi в спектре достигает 20%.

Преимущества алгоритма Ь~' по сравнению с методом вычитания «медленных» кинетических компонент экспериментально показаны на примере анализа температурной зависимости времён релаксации ядер 71л в гетерогенной мембранной системе.

Четвёртая глава («Релаксация намагниченности ядер ]Н воды и углеводородов в среде гранулярных силикатов») посвящена изучению влияния размеров гранул порообразующей среды на скорости продольной и поперечной релаксации ядер 'Н в поровых жидкостях. Для модельных сред из стеклянных шариков различного размера определены коэффициенты релаксации р\ и р2 для протонов воды и ряда УВ (гексан, гептан, нонан, декан и бензол). Показана применимость полученных значений рх и р2 для определения распределения пор по размерам в среде природных силикатов на примере песка и нефтесодержащего керна.

Помещение жидкости в пористую среду приводит к укорочению времен Т\ и Т2 по сравнению с соответствующими временами в свободной жидкости (рис. 4).

1 Т2, Вода ^ Т2, Декан

о1

ДМа

€

0.5

0 12 3 4 Время поперечной релаксации, с

Т1, Вода

0 1 2 3 4 5 Время поперечной релаксации, с

Т1, Декан

0 1 2.3 4 5 Время продольной релаксации, с

-е—57мкм —■—43!

0 1 2.3 4 5 Время продольной релаксации, с

-х- свободный флюид

Рис. 4. Спектры времён поперечной и продольной релаксации свободных воды и декана, а также свободных флюидов в среде стеклянных шариков.

Для интерпретации релаксационных данных обычно используется феноменологический подход, согласно которому, в поровом флюиде выделяется приповерхностный слой толщиной Ь, скорость релаксации 1/Т в котором существенно превышает скорость релаксации в объемной фазе жидкости. При быстром диффузионном обмене между поверхностной и объемной фазами (т0«Т, та -характерное время диффузионного обмена) наблюдаемая скорость релаксации представляет собой взвешенную сумму скоростей поверхностной (\/Т$) и объемной (1/7» релаксации:

т V V )ту V г5 4 "

где 5 - площадь поверхности поры, V — объём поры, (7*$ « Ту).

Введение коэффициента релаксации р = Ь(—---—) позволяет

ТУ

переписать выражение (7) в виде:

+ (8)

Т Ту V к '

Установлено, что ускорение релаксации флюида в пористой среде из стеклянных шариков не сопряжено с изменением характера кинетики: спад намагниченности происходит по экспоненциальному закону (аналогично релаксации в свободной жидкости) и спектры времен Тх и Т2 являются однокомпонентными (рис. 4). Однокомпонентный спектр времен релаксации жидкости в однородной по размеру гранул пористой среде подтверждает наличие быстрого диффузионного обмена протонов между поверхностным слоем и объемной компонентой порового флюида.

Используя понятие калибра поры с1=У/Б, для изменения скорости релаксации К = 1/Т флюида внутри порового пространства относительно скорости релаксации свободного флюида на основании выражения (8) получаем:

АЯ=1--±- = р- (9)

Т Ту й у '

Поры, образуемые плотноупакованными шарами диаметром Н, характеризуются калибром

= 0.12// (10)

'роге л

Варьируя размер шариков, образующих пористую среду, с использованием соотношения (9) получена зависимость АЯ(\/с1) для различных флюидов. Результаты этих экспериментов, приведенные на рис. 5, показывают линейную зависимость между изменением скорости релаксации А/? и величиной 1 /У, что доказывает справедливость выражения (8) для всех исследованных жидкостей. Линейный характер зависимости Д7?(1 /с{) позволил определить

коэффициенты релаксации р1 и р2, значения которых приведены в таблице 1.

? 4

оГ з

2 1

0,05 0.10 0.15

1 га, 1/мкм

0,05 0.10 0.15 1/С1, 1/МКМ

ого

Рис. 5. Зависимости от величины 1Ш скоростей поперечной и продольной релаксации воды и углеводородов в пористой среде из стеклянных шариков.

Обнаружено, что значения р2 для всех исследованных жидкостей имеют один и тот же порядок величины (/>> ~ 10"3 см/с), в то время как значения р\ примерно на два порядка меньше величины р2 для всех углеводородов {р\ ~ 10~5 см/с), и только для воды р\ ~ р1.

Таблица 1. Значения релаксационных коэффициентов р\ и р2 для различных флюидов в пористой среде, образованной стеклянными шариками.

Поровый флюид Рг, (см/с) Ри (см/с)

Н20 3.6-10° 3.1-10-3

с6н6 2.2-10"3 6.5-10"5

С6Н14 2.7-10"3 7.2-10"ь

с7н16 1.7-10"3 2.0-Ю-5

С9Н20 1.5-Ю-3 1.4-Ю-5

СюН22 1.4-10"3 6.0-Ю-5

Полученные для модельной среды значения р2 использованы для анализа релаксационных данных в случае флюида, заполняющего реальную породу. Кинетика спада поперечной намагниченности флюида, пропитывающего песок, в отличие от флюида в модельной среде, является многокомпонентной (рис. 6а).

Рис. 6. Спектр времен протонной релаксации Т2 для воды и декана в кварцевом песке (а) и распределение пор по размеру, полученное на основании спектров времён релаксации (б).

В спектрах водного и УВ флюида обнаруживаются три компонента, которые интерпретируются как свидетельство наличия трех характерных калибров пор. Использование соотношения (9) применительно к релаксационным спектрам дает оценку распределения пор по размерам для кварцевого песка (таблица 2).

Таблица 2. Характерный калибр пор й и объёмная доля соответствующих пор р в кварцевом песке, полученные на основании спектров времен поперечной релаксации протонов Т2.

Флюид Вода Декан

Т2 (мс) 8,7 2,7 1,2 17,6 5,2 1,3

с/ (мкм) 0,31 0,10 0,04 0,25 0,07 0,02

р 0,596 0,331 0,073 0,521 0,342 0,137

Функция распределения пор по размерам в кварцевом песке (рис. 66) с учётом плотности цельного материала и насыпной плотности позволяет оценить удельную поверхность этой пористой среды значением: = 2,11-Й,82 м2/г. Непосредственное измерение удельной поверхности песка по адсорбции азота с помощью установки Сорби М8 дает величину Бт = 1,97± 0,05 м2/г, которая в пределах погрешности 7% совпадает с рассчитанной на основании

релаксационных данных для воды.

Таким образом, измерение релаксации протонной намагниченности флюида, заполняющего пористую среду, является эффективным методом исследования структурных размеров порового пространства. Среда из стеклянных шариков, заполненных флюидом, адекватно моделирует релаксационные процессы в реальных

неконсолидированных силикатных породах. Коэффициенты релаксации, полученные на модельных средах, пригодны для определения распределение пор по размерам в реальном кварцевом песке. Удельная поверхность песка, рассчитанная на основании релаксационных данных для воды, совпадает с измеренной адсорбционным методом.

В Пятой главе («Механизм поверхностной релаксации протонной намагниченности воды и углеводородов в среде силикатов») приведены результаты исследования механизмов продольной и поперечной релаксации ядер 'Н флюидов на поверхности силикатов.

Основным каналом поверхностной релаксации в природных пористых средах считается релаксация, обусловленная примесными парамагнитными центрами (ПЦ). Спектры ЭПР, полученные для стеклянных шариков различного размера, приведены на рис. 7. В спектрах наблюдаются линии при 1585 Гс ^ = 4,24) и при 33633366 Гс (£= 1,998-5-1,999). Уменьшение размера шариков изменяет соотношение интенсивностей сигналов в пользу линии g~ 2.

В (кв)

Рис. 7. Спектры ЭПР (производная сигнала поглощения) стеклянных шариков различного диаметра и фоновый сигнал (внизу)

Линия с g = 30/7 = 4,28 соответствует примесным ионам Ре3+, находящимся в кристаллическом поле с искаженной тетраэдрической симметрией. Сигналы, характеризуемые я~2, связывают с ненасыщенностью химических связей в силикатных группировках типа ^Ю-, =8ЮО. Наиболее вероятным местом локализации

подобных структур со свободными связями является поверхность

стекла. Две линии в спектре ЭПР стеклянных шариков отвечают двум типам ПЦ: ионы Ре3+, расположенные в объеме, и поверхностные группировки диоксида кремния с ненасыщенными связями. Наблюдаемая зависимость соотношения компонент ЭПР-спектра от диаметра шариков согласуется с данным предположением, поскольку при этом массовая доля оксида железа, равномерно распределенного по объёму, не зависит от диаметра шариков Н, тогда как удельная поверхность /7/ и при постоянной массе образца количество поверхностных кремний-оксидных ПЦ увеличивается с уменьшением диаметра.

Двойное интегрирование линий в спектре ЭПР позволяет определить интенсивность сигналов и оценить долю а каждого из двух видов ПЦ в расчете на моль 8Ю2. Для ионов железа с учетом поправок, определяющих вклад в наблюдаемый сигнал с изотропным g-фактором, получаем «Ге ~ 10"4, что согласуется с паспортными данными стекла (0,2% Ре203). Оценка доли кремний-оксидных ПЦ дает ~ 10"7-10"6, что составляет 1-10% от доли поверхностных молекул 8102 в шариках диаметром 100 мкм. Таким образом, полученные оценки указывают на то, что один кремний-оксидный ПЦ приходится на 10-100 поверхностных молекул ЗЮ2.

Дополнительным подтверждением наличия ПЦ на поверхности диоксида кремния являются результаты исследования каталитического окисления 2,5-дитретбутилгидрохинона (ТБГХ) (рис. 8) поверхностью стеклянных шариков. Гидрохиноны вступают в реакцию окисления в присутствии свободных радикалов. Данные 'Н ЯМР спектроскопии высокого разрешения показывают, что в спиртовом растворе ТБГХ, помещённом в среду стеклянных шариков, образуется значительное количество молекул 2,5-дитретбутилбензохинона (ТББХ) (рис 8), который является продуктом радикального окисления ТБГХ. За 24 часа в среде стеклянных шариков около половины молекул ТБГХ превращаются в молекулы ТББХ, что свидетельствует о каталитической роли поверхности шариков, поскольку количество молекул ТББХ, образовавшихся в ходе данной реакции, более чем в 100 раз превышает число поверхностных ПЦ.

з тагл

чЧ /; I

г ТББХ

ррт

6 5-1321 ррт

Рис 8 Спектры 'Н ЯМР исходного раствора ТБГХ (а) и раствора ТБГХ, с примесью ТББХ взятого из среды стеклянных шариков (б).

Существенная зависимость времен релаксации Т\ и Т2 от характерного размера порового пространства <1 указывает на то, что протоны флюида теряют коллективную поляризацию на поверхности. С помощью ЭПР установлен механизм деполяризации в силикатах, обусловленный наличием на поверхности оксида кремния ПЦ, поверхностная доля которых в стеклянных шариках составляет /~0,01н-0,1.

Основные экспериментальные результаты рассмотрены в рамках простой модели, согласно которой молекулы флюида в процессе самодиффузии, характеризуемой коэффициентом О, за характерное диффузионное время т0~ с^Ю попадают в пристеночный слой поры толщиной Ъ и адсорбируются на стенке на время тА. В течение этого времени на протоны адсорбированных молекул действует создаваемое ПЦ магнитное поле В\, индукцию которого оценим как В] ~ //с/63 (//£, - магнитный момент электрона, оценочно принимаемый равным магнетону Бора). Влияние поля В\ приводит к тому, что спины протонов в адсорбированных молекулах приобретают дополнительный фазовый сдвиг по углу прецессии на случайную величину ±5(р\, где 8<р\ - т//В1 (у - гиромагнитное отношение протона). Примем, что после десорбции молекула флюида оказывается на

расстоянии ~Ь от поверхности. Вероятность достижения противоположной стенки в результате происходящего после этого самодиффузионного процесса \\>~Ыс1«\, поэтому молекула совершит еще п посадок на ту же стенку (п = сПЬ) за время гь для которого можно получить оценку:

пЪ2 йЬ Т!~ —~ —~и/то « т0

При этом число контактов с ПЦ составит ~/п. Для протонов воды и УВ 8ц>\ « 1, поэтому движение по углу <р носит диффузионный характер. За всю серию столкновений с одной стенкой прецессия спинов приобретает дополнительный фазовый сдвиг на угол

Затем происходит серия из ~п столкновений с противоположной стенкой и далее процесс повторяется. За время / происходит N = г/т0 диффузионных переносов молекулы от одной стенки до другой, за это время суммарный набег фазы составит

Д(р~[(8<р)2 ■ Щ1!2 (11)

Соотношение (11) позволяет оценить время поперечной релаксации, полагая, Аср~ 1:

ТГГа^2 (»2)

Выражение (12) правильно отражает обнаруженную экспериментально линейную зависимость между скоростью поперечной релаксации и величиной (1 /<Т) (рис. 5).

Продольная релаксация - это диффузия ядерных магнитных моментов по энергетическим уровням. Характерное время Т\ такого процесса обусловлено переходами с обязательным изменением энергии ядерных спинов. В то же время, поперечная релаксация обусловлена переходами как связанными с изменением энергии, так и не связанными с изменением энергии. Из общих представлений теории элементарных процессов можно заключить, что при большом времени адсорбции, когда выполняется условие % = й)0тА » 1 переходы с изменением энергии будут адиабатически маловероятны и время Т\ будет экспоненциально велико:

Г1~Г2е^ » Т2 (13)

Это соотношение справедливо и при £ < 1, в этом случае Тх ~ Т2.

Для воды, заключенной в поровом пространстве шариков, имеет место Г,~ полагая соо = 2л-20 МГц = 1,3-108 рад/с, используя соотношение (9) из условия 1 получаем тА < 1-10"9с. В случае углеводородного порового флюида Тх ~ 10Т2 и из соотношения (13) получаем тА ~ 2,2-10"8 с.

Для проверки соотношения (13) были измерены времена продольной и поперечной релаксации ядер 'Н воды и У В в среде стеклянных шариков диаметром 57 мкм с использованием трех различных магнитов, обеспечивающих рабочие частоты 10 МГц, 20 МГц и 500 МГц на ядрах !Н. Полученные результаты подтверждают адекватность соотношения (13) (таблица 3.)

Таблица 3 Значения отношения времён Г^/Т^ полученные для воды и декана для трёх значений рабочих частот спектрометра, а так же расчётные значения параметра еТАШ°. Для расчёта использованы значения га:10'9 с и 2,2*10~8с для воды и декана соответственно.

Вода Декан

со0 7Ъ/7м етАшО Т]$1Т2В е*Аш 0

10 МГц 1,6 1,06 4,6 3,98

20 МГц 1,6 1,13 14,2 15,85

500 МГц 18.9 23,10 Ю30

'Измеренное при данной частоте резонанса значение времени Т\ для декана в среде стеклянных шариков практически не отличается от значения Т1 у свободного декана.

Большее значение величины та молекул УВ на поверхности стекла, по сравнению с га молекул воды согласуется с данными измерения углов смачивания на поверхности стекла. Измеренные значения краевых углов для декана и воды на поверхности стекла составили 6,4±0,3 и 42,8±1,9, соответственно.

Таким образом, показано, что парамагнитный механизм релаксации протонов на поверхности гранул из диоксида кремния не требует присутствия инородных примесей в силу наличия на поверхности БЮ2 парамагнитных центров, образованных ненасыщенными связями в структуре диоксида кремния. Соотношение времён продольной и поперечной релаксации порового флюида при заданной резонансной частоте соо определяется временем адсорбции га молекул флюида на поверхности поры.

В Шестой главе («Анализ состава порового флюида в силикатной среде») проводится сравнение аналитических возможностей методов твердотельного ЯМР высокого разрешения и ЯМР-релаксации в определении состава порового флюида (определение долей водной и УВ фаз).

Определение состава порового флюида методом ЯМР высокого разрешения затруднено из-за уширения линий вследствие значительного увеличения скорости поперечной релаксации протонов.

Установлено, что в случае пористой среды из одинаковых по размеру гранул, уширение сигнала относительно невелико, что позволяет наблюдать разрешённые сигналы ЯМР для смеси воды и УВ стандартным методом (рис. 9а). В случае неоднородных гранул (песок) для получения разрешённых спектров необходимо применение методики «вращения под магическим углом» (MAS -ЯМР) (рис 9в).

У.в.

20 0 -20 ррт

Н20

У.В.

40 20

и,о

У.В.

О -20 40 20 0 -20 40 20 0 -20

ррт ррт ррт

Рис. 9 Спектры 'Н-ЯМР смеси декана с водой, помещённых в среду стеклянных шариков диаметром 57 мкм (а — без вращения, б — MAS - ЯМР 5 кГц), и песка (в - без вращения, г - MAS - ЯМР 5 кГц). Символом * отмечены сателлиты, обусловленные вращением образца.

Метод МА8-ЯМР требует дорогостоящего оборудования и не является в полной мере невозмущающим. Центробежное ускорение при частоте вращения 5 кГц достигает 7-10^, что приводит к седиментационному разделению компонент флюидов

Предложен альтернативный метод определения состава флюида, основанный на релаксационных данных. Для поверхностной продольной релаксации воды и УВ в силикатах выполняется соотношение: Т]3вода « Г15УВ. Анализ спектров времён Тх и Т2 флюидов в пористой среде позволяет определять не только распределение пор по размерам, но и доли водной и УВ компонент в жидкой фазе. Спектр времён Т2 независимо от состава флюида отражает распределение пор по размерам (рис. 4), поскольку коэффициенты р2 для воды и УВ являются величинами одного порядка. Наличие в спектре времён Т\ компоненты, для которой Т\~<Т2>, свидетельствует о присутствии в образце воды. Наличие компоненты Т\ » <Т2> свидетельствует о присутствии в образце УВ фазы (<Т2> - характерное значение из диапазона наблюдаемых величин Т2).

Иллюстрацией предложенного подхода являются данные, приведённые на рис. 10. В спектре времён Т\ смеси воды и декана в

среде песка присутствует два пика (рис. 10а). Меньшее значение Т\ сравнимо со значением <Т2> в спектре времён поперечной релаксации, а большее значение Т\ значительно превышает значение <Т2>. Отнесение пиков выполнено с помощью экспериментов с флюидом состава УВЯЗЛО (рис. 106). Замена Н20 на Б20 позволяет получить спектры времён релаксации ядер 'Н только УВ в образце, содержащем и углеводородную, и водную фазы.

Рис. 10. Спектры времён продольной и поперечной релаксации смеси воды и декана (а) и дейтерированной воды и декана (б) в среде песка.

Таким образом, интерпретация релаксационных данных для порового флюида на основании разработанной модели релаксации позволяет определять соотношение водной и УВ фаз в составе заполняющей поры жидкости.

Выводы

1. Установлено, что для надежного определения спектра времен релаксации с помощью обратного преобразования Лапласа экспериментальный временной ряд должен содержать —50 дискретных значений, соотношение сигнал/шум должно быть не ниже 100, при этих условиях разрешимы компоненты с соотношением времён релаксации Т,/Тк ~2.

2. Получены калибровочные значения коэффициентов релаксации Р\ и р2 для воды и ряда углеводородов (гексан, гептан, нонан, декан, бензол) в модельной среде на основе диоксида кремния, позволяющие определять распределение пор по размерам в

— Ь

10 100 1000 10000 Время релаксации, мс

естественных силикатных породах с характерным калибром пор в диапазоне 0,01 - 100 мкм.

3. Впервые физическими и химическими методами анализа показано, что парамагнитный механизм релаксации протонов на поверхности гранул из диоксида кремния не требует присутствия инородных примесей в силу наличия на поверхности Si02 парамагнитных центров, образованных ненасыщенными связями в структуре диоксида кремния.

4. Разработана физическая модель поверхностной релаксации протонной намагниченности, согласно которой соотношение времён релаксации Т\/Т2 порового флюида при заданной резонансной частоте со0 определяется временем адсорбции га молекул флюида на поверхности поры. Энергия адсорбции и, тем самым, та молекул воды и углеводородов на поверхности силикатов определяется в основном дисперсионным взаимодействием, зависящем от объёма адсорбируемой молекулы.

5. Предложен алгоритм качественного определения состава флюида в пористой среде (вода, углеводороды или их смесь), основанный на сопоставлении спектров времён Tt и Т2 релаксации протонной намагниченности, зарегистрированной на частотах 10-20 МГц.

Список опубликованных работ по теме диссертации

1. Гуденко С. В., Максимычев А. В., Перепухов А. М., Родин В. В. ЯМР - Спектроскопия // МФТИ. Москва. 2011. 123 с.

2. Перепухов А. М., Шестаков С. Л. Алгоритм обратного преобразования Лапласа для обработки сложных релаксационных зависимостей // Труды МФТИ. 2010. Т. 2. №. 2. С. 35^0.

3. Перепухов А. М., Кишенков О. В., Гуденко С. В., Максимычев А. В., Меньшиков Л. И., Александров Д. А., Коростылёв Е. В. Исследование порового пространства силикатов и протонной релаксации флюида магнитно-резонансными методами. // Труды МФТИ. 2013. Т. 5. № 3. С. 154-163.

4. Negrimovsky V., Komissarov A., Perepukhov A., Suponitsky К., Perevalov V., Lukyanets Е. Phthalocyanines with phosphonate moiety via C-nucleophilicsubstitution in phthalonitriles // Journal of Porphyrins and Phthalocyanines J. Porphyrins Phthalocyanines 2013; 17: 587-595

5. Перепухов A. M. Исследования химической и изотопной чистоты синтеза октановой кислоты. // Труды 50-й научной конференции

МФТИ Часть IV Молекулярная и биологическая физика 2007 г Т.2,

C.125-126

6. Шестаков C.JL, Перепухов A.M., Павлов A.A. Изучение механизма проводимости перфторированных ионообменных мембран в различных солевых формах методами ядерного магнитного резонанса высокого разрешения. // Труды 51-й научной конференции МФТИ «Современные проблемы фундаментальных и прикладных наук»: Часть IV. Молекулярная и биологическая физика. 2008. С. 143

7. Максимычев А. В., Анохин Е. М., Перепухов А. М., Шнер В. Ф., Эльман А. Р., Кузьмин С. Г. Применение методов -ЯМР и масс-спектрального анализа к исследованиям ,3С-меченых препаратов для диагностики гастроэнтерологических и онкологических заболеваний // Труды 52-й научной конференции МФТИ Часть IV Молекулярная и биологическая физика 2009 г Т.2, С.132-133

8. Перепухов A.M. Алгоритм обратного преобразования Лапласа для обработки сложных релаксационных зависимостей (тезисы) // Труды 52-й научной конференции МФТИ Часть IV Молекулярная и биологическая физика 2009 г Т.2, С. 136-137

9. Перепухов A.M. Алгоритм обратного преобразования Лапласа для обработки сложных релаксационных зависимостей. //6-я Зимняя молодежная школа-конференция Магнитный резонанс и его приложения. Материалы конференции 2009 г. С. 98-99

10. Перепухов А. М., Кишенков О.В., Эльман А. Р., Овсянникова Л. В., Максимычев А. В. Изучение механизма реакции карбонилирования замещённых одибромбензолов в эфиры фталевой кислоты методами двумерной спектроскопии ЯМР. Труды 54-й научной конференции МФТИ Молекулярная и биологическая физика 2011 г, С. 176-177.

11. Кишенков О. В., Перепухов А. М., Максимычев А. В. Исследование характеристик модельных пористых сред методами ЯМР. // Труды 54-й научной конференции МФТИ Молекулярная и биологическая физика 2011 г, С. 179-180.

12. Kishenkov O.V., Perepukhov A. M., Maximychev A.V., Alexandrov

D.A. The Use of Nuclear Magnetic Relaxation for the Determination of Internal Geometry of Porous Media. // Magnetic Resonance And Magnetic Phenomena In Chemical And Biological Physics. The School for Young Scientists. Novosibirsk, Russia. 2012, — P. 49.

13. Перепухов A.M., Кишенков O.B., Максимычев A.B., Гуденко C.B. Исследование механизма ядерной магнитной релаксации порового

флюида. // Труды 55-й научной конференции МФТИ Молекулярная и биологическая физика 2012 г, С.119-120.

14. Кишенков О. В., Перепухов А. М., Максимычев А. В., Александров Д. А. Релаксация флюидов в пористых средах в режиме быстрого диффузионного обмена. // Труды 55-й научной конференции МФТИ Молекулярная и биологическая физика 2012 г, С.121-122.

15. Перепухов А. М., Кишенков О.В., Эльман А. Р., Овсянникова Л. В., Максимычев А. В. Изучение механизма синтеза фталатов карбонилированием дибромбензолов с использованием 13СО. // Труды 55-й научной конференции МФТИ Молекулярная и биологическая физика 2012 г, С.182-183.

Перепухов Александр Максимович

Применение методов ЯМР к исследованию геометрических характеристик порового пространства гранулярных силикатов и свойств поровых флюидов

Авторефрерат

диссертации на соискание ученой степени кандидата физико-математических наук

Подписано в печать 15.11.2013. Формат 60x84 1/16 Усл. печ. л. 1,25. Тираж 100 экз. Заказ №289

Федеральное государственное автономное образовательное учреждение высшего профессионального образования «Московский физико-технический институт (государственный университет)» 141700, Московская обл., г. Долгопрудный, Институтский пер., 9

Отдел оперативной полиграфии «Физтех-полиграф» тел.: (495)408-84-30

Текст научной работы диссертации и автореферата по физике, кандидата физико-математических наук, Перепухов, Александр Максимович, Долгопрудный

Моковский физико-технический институт

На правах рукописи

04201 453154

Перепухов Александр Максимович

Применение методов ЯМР к исследованию геометрических характеристик порового пространства гранулярных силикатов и свойств поровых

флюидов

01.04.17 - химическая физика, горение и взрыв, физика экстремальных

состояний вещества

Диссертация на соискание ученой степени кандидата физико-математических наук

Научный руководитель Максимычев А. В., д. ф.-м. н., профессор

Долгопрудный - 2013

Оглавление

Введение..................................................................................................................................................................5

Глава 1. Методы анализа пористых сред с помощью ЯМР............................8

1.1. Пористые среды....................................................................................................................................8

1.2. Ядерный магнитный резонанс..................................................................................................9

1.3. Обзор методов анализа пористых сред с помощью ЯМР................................12

1.4. ЯМР релаксация....................................................................................................................................14

1.5. Механизмы ЯМР релаксации....................................................................................................16

1.6. ЯМР релаксация в пористых средах....................................................................................22

1.7. ЯМР каротаж............................................................................................................................................30

1.8. Обзор экспериментальных данных по релаксации в пористых

средах........................................................................................................................................................................36

1.9. Задачи исследования........................................................................................................................43

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

2.1. Модельные и натуральные пористые среды................................................................45

2.2. Поровые флюиды..................................................................................................................................46

2.3. ЯМР релаксометрия............................................................................................................................48

2.4. ЯМР высокого разрешения..........................................................................................................52

2.5. Электронный парамагнитный резонанс..........................................................................54

2.6. Растровая электронная микроскопия и электронно-зондовый микроанализ........................................................................................................................................................55

2.7. Адсорбционный метод определения площади поверхности пор..............56

2.8. Измерение смачиваемости жидкостей на поверхности стекла

методом анализа видеоизображения капли..............................................................................58

Глава 3. Алгоритм обратного преобразования Лапласа для

обработки сложных релаксационных зависимостей..............................................61

3.1. Многоэкспоненциальные кинетические зависимости........................................61

3.2. Алгоритм обратного преобразования Лапласа..........................................................63

3.3. Применение обратного преобразования Лапласа к модельным функциям................................................................................ 64

3.4. Обработка экспериментальных данных..................................... 69

3.5. Резюме.............................................................................. 71

Глава 4. Релаксация намагниченности ядер Н воды и

углеводородов в среде гранулярных силикатов............................. 73

4.1. Кинетика релаксации ядерной намагниченности в свободном

флюиде................................................................................... 73

4.2. Релаксация протонной намагниченности в поровом флюиде.......... 77

4.3. Влияние самодиффузии молекул флюида на скорость поперечной релаксации протонной намагниченности......................................... 79

4.4. Коэффициенты продольной и поперечной релаксации намагниченности 'Н воды и углеводородов..................................... 81

4.5. Распределение пор по размерам в натуральных силикатах............. 85

4.6. Верификация данных ЯМР адсорбционным методом.................... 88

4.7. Резюме.............................................................................. 92

Глава 5. Механизм поверхностной релаксации протонной

ч намагниченности воды и углеводородов в среде силикатов............. 93

5.1. Свойства поверхности стеклянных шаров.................................. 93

5.2. Парамагнитные свойства гранулярных силикатов........................ 96

, 5.3. Радикальное окисление гидрохинонов в среде гранулярных

силикатов................................................................................ 103

5.4. Модель поверхностной поперечной релаксации намагниченности порового флюида...................................................................... 108

ч 5.5. Продольная поверхностная релаксация намагниченности порового

флюида................................................................................... 112

5.6. Влияние смачиваемости на скорость продольной релаксации......... 113

5.7. Резюме.............................................................................. 116

Глава 6. Анализ состава порового флюида в силикатной среде................117

6.1. Применение метода ЯМР высокого разрешения для анализа состава порового флюида............................................................................................................................................117

6.2. Использование данных ЯМР релаксации для определения состава порового флюида............................................................................................................................................120

6.3. Сравнение спектров времён Т\ и Т2 воды и декана в среде песка............122

6.4. Влияние резонансной частоты на релаксацию порового флюида............125

6.5. Резюме............................................................................................................................................................128

Основные выводы....................................................................................................................................130

Список опубликованных работ по теме диссертации..........................................131

Благодарности..............................................................................................................................................134

Список литературы..................................................................................................................................135

Ведение

Метод аттестации пористых сред, основанный на измерениях времён продольной (71]) и поперечной (Г2) ЯМР релаксации ядер ]Н флюида, заполняющего поровое пространство, является одним из самых распространённых методов ЯМР, применяемых для изучения пористых сред. В релаксационном методе распределение пор по размерам получают из спектра времён релаксации (Т\ или Т2), полагая, что скорость релаксации пропорциональна отношению площади поверхности поры к её объёму (S/V):

1 5

- ~ Pl 2 ~ (1.1)

Т1Л ' V у J

Коэффициенты пропорциональности рх и рг в соотношении (1.1) обычно рассматриваются как феноменологические параметры для данной пары твёрдая стенка/флюид. Большинство работ по релаксационным методам проводились или на реальных кернах (диапазон размеров пор в которых примерно 100 нм - 1 мм), или на модельных объектах - пористых стёклах с порами значительно меньшего размера пор (1-10 нм). При этом имеются существенные расхождения в данных, полученных для пористых стёкол и натуральных песчаников. Таким образом, для понимания механизмов релаксации ядер ]Н флюидов в натуральных пористых средах необходимы исследования на модельных системах с характерным размером пор ОД мкм -1 мм.

Значительный интерес представляет возможность различения водного и углеводородного порового флюидов по соотношению времён Т\/Т2. Однако, литературные данные по данному поводу весьма противоречивы.

Таким образом, интерпретация релаксационных данных в случае пористых сред сопряжена со сложностями, обусловленными зависимостью скоростей релаксации не только от размеров пор, но и от состава флюида, температуры, смачиваемости флюидом поверхности, наличия парамагнитных примесей, напряжённости постоянного магнитного поля и др. Для выявления основных закономерностей релаксационных процессов в пористых средах необходимо построение физической модели, описывающей зависимости скоростей продольной и поперечной релаксации ядер 'Н воды и углеводородов от указанных параметров. Для модельных исследований силикатных пород целесообразно использовать стеклянные шарики, поскольку диоксид кремния является основным компонентом, как стекла, так и натуральных песчаников, а характерный размер коллекторных пор (0,1 мкм - 100 мкм) легко получить, варьируя размер шариков.

Цель и задачи диссертационной работы

Для достижения поставленной цели ставились следующие задачи:

• Определение для ядер зависимости времён продольной и поперечной релаксации ядерной намагниченности от характерного размера пор для широкого набора флюидов в модельной среде стеклянных шариков. Применение релаксационных данных, полученных для модельной среды, к анализу порового пространства реальных пород. Верификация релаксационных данных о геометрии порового пространства независимым адсорбционным методом.

Глава 1. Методы анализа пористых сред с помощью ЯМР (Литературный обзор)

1.1. Пористые среды

Пористые среды представляют собой среды, пронизанные системой сообщающихся между собой пустот (пор), имеющих различную форму и размеры. Существуют пористые среды как природного, так и синтетического происхождения. К природным пористым средам можно отнести горные породы (песчаники, известняки и др.), древесину, костную ткань и др. К искусственным пористым средам можно отнести углеродные сорбирующие материалы, полимерные и биополимерные материалы, пористые стёкла и т.д.

С геометрической точки зрения пористые среды делятся на две большие группы: гранулярные и трещиноватые. Поры в гранулярных пористых средах образованы пустотами между твёрдыми гранулами. Трещиноватые среды представляют собой систему развитых трещин в твёрдом материале, густота которых зависит от состава и свойств пород.

Большой интерес к пористым средам связан, прежде всего, с добычей нефти, которая заключена в поровом пространстве грунтовых гонных пород. Нефтяные месторождения чаще всего представляют собой пласты терригенных и карбонатных осадочных горных пород (песчаников, известняков, глин), представляющих собой скопления зерен минералов скрепленных цементирующим материалом. Поровое пространство осадочных горных пород — сложная система сообщающихся межгранулярных пустот. Размеры пор в песчаных породах составляют обычно единицы или десятки микрометров.

Важнейшими характеристиками нефтесодержащей породы, определяющими ее нефтеотдачу, являются 1) содержание углеводородов в заполняющем пласт флюиде, 2) геометрические характеристики порового

пространства и 3) трансляционная подвижность водной и углеводородной компонент флюида. Информация об этих характеристиках среды может быть получена с использованием различных экспериментальных методов, основу которых составляет явление ядерного магнитного резонанса (ЯМР). Метод ЯМР успешно используется в нефтегазовой промышленности для аттестации порового пространства горных пород [1]. Флюид, заполняющий пористую среду, содержит ядра 'Н, что позволяет изучать свойства флюида находящегося в пористой среде методами 'Н-ЯМР.

1.2. Ядерный магнитный резонанс

Спектроскопия ядерного магнитного резонанса (ЯМР) - один из основных методов физико-химического анализа, используемый для однозначной идентификации структуры молекул. Согласно правилам IUP АС, ЯМР является одним из обязательных методов исследования, подтверждающих структуру нового химического соединения. Спектроскопия ЯМР основана на резонансном поглощении подводимой электромагнитной энергии системой ядерных магнитных моментов, помещенных в постоянное внешнее магнитное поле. [2]

Ядерный магнитный момент связан с наличием у атомного ядра ненулевого магнитного момента, обусловленного тем, что протоны и нейтроны, образующие атомное ядро, обладают собственными механическими моментами - спинами. Магнитный момент ц ядра обладающего механическим моментом I выражается следующим соотношением Где h - постоянная Планка, у - гиромагнитное

отношение, g (g-фактор) - индивидуальная для каждого ядра величина, ßyv -ядерный магнетон Бора. Магнетон Бора - это постоянная величина, выражаемая (в системе СГС) следующим соотношением: ßN=eh/2сМр, где е -элементарный заряд, Мр - масса протона, а с - скорость света. Численное

значение ядерного магнетона равно = 5,0508 10 эрг/Гс (СГС) = 5,0508 10"27Дж/Тл (СИ).

Энергия Е магнитного момента, находящегося в постоянном магнитном поле В , выражается скалярным произведением:

Е = ~(\1,В0) = -|х Во собО = - g |3уу£0 тг.

где 0 - угол между направлениями векторов (I, и В0, а т2 - проекция спина на

ось г, совпадающую с направлением внешнего магнитного поля (рис. 1.1).

Для ядра обладающего спином I, возможно всего 2/+1 значение проекции

спина т2 от -/ до +/, т.е. 21+1 значение энергии магнитного момента.

Явление ядерного магнитного резонанса можно наблюдать только на ядрах

обладающих ненулевым спином ('Н, 2Т>, 19Р, 13С). Нулевым спином обладают

ядра, в которых содержится чётное число протонов и чётное число

нейтронов, например изотоп 12С, с этим связана сложность получения

углеродных спектров, поскольку -99% всего углерода, содержащегося в

природе - это изотоп С. [3]

Рис. 1.1. Ларморова прецессия ядерного спина.

При воздействии на ядро, помещённое в постоянное магнитное поле, переменным электромагнитным полем резонансной частоты, возможен переход между состояниями с различной проекцией магнитного момента т2 на ось квантования. Квантовая теория даёт правило отбора для таких

переходов: Дш2 = ±1. Следовательно, резонансная частота поглощения электромагнитного поля Уо, при которой происходят переходы, определяется ношением:

= АЕ/к = g рд'^о/Ь = уВ0/2п.

Атомные ядра, исследуемые методом ЯМР, входят в состав молекул и находятся в окружении электронных оболочек. Заполненные электронные оболочки стабильных молекул, как правило, содержат спаренные электроны, суммарный спин которых равен нулю. В этих условиях электронные оболочки проявляют диамагнитные свойства При наложении внешнего магнитного поля В0 , электроны заполненных молекулярных орбиталей индуцируют собственное магнитное поле Ве напряжённость которого пропорциональна напряжённости поля В0 , а направление противоположно направлению поля В0 : Ве =-о В0. Таким образом, каждое ядро молекулы, помещённой в постоянное магнитное поле В0, находится в локальном магнитном поле 2?/ос, которое складывается из поля В0 и поля Ве, создаваемого электронным окружением: [4]

В1ос - В0 +Ве= В0{ 1 -су), где о - константа экранирования. Характерная величина константы экранирования а < 10"5 для протонов, а < 10"3 для более тяжелых ядер.

Электронное экранирование приводит к тому, что резонансные частоты поглощения для ядер в молекуле смещаются относительно частоты, на которой наблюдался бы резонанс для свободного ядра и определяются локальными полями, в которых находятся ядра:

у!=уД(о,/27г = уй0(1-а)/27г.

Смещение частоты Ду = V) - = уЛ0а/2я называют химическим сдвигом. На практике химический сдвиг обычно выражается в безразмерных велиинах 5 как отношение сдвига частоты относительно частоты эталонного соединения, к частот�