CC BY
31
СТРУКТУРНО-ДИНАМИЧЕСКИЙ АНАЛИЗ
ИМПУЛЬСНОИ ЯМР-СПЕКТРОСКОПИЕЙ
В ИССЛЕДОВАНИИ ОСТАТОЧНЫХ
НЕФТЯНЫХ СИСТЕМ РАЗЛИЧНОЙ ПРИРОДЫ
КЕМАЛОВ А.Ф., КЕМАЛОВ Р.А. Научно-технологический центр КГТУ Казань
МУЛЛАХМЕТОВ Н.Р., ГАЗИЗОВ Р.Р. «Природные битумы» e_maN. kemalov@mail ru
ФАИЗРАХМАНОВ А.Т.
КЛЮЧЕВЫЕ СЛОВА Тяжелые нефтяные остатки, импульсная ЯМР-спектроскопия, битумы, парафи-
KEYWORDS: ны, физико-химические свойства, глубокая переработка, нефтяные дисперсные
системы. Heavy oil residues, pulsed NMR- spectroscopy, bitumen, paraffins, Physic-ochemical characteristics, deep conversion, petroleum disperse systems
Согласно положениям физико-химической механики, разработанным акад. П.А. Ребиндером, правильное представление о структуре позволяет научно решать вопросы создания материала с заданными свойствами.
Возникает необходимость изучения структуры битумов в широком диапазоне соотношений его структурных компонентов и влияния на образование структур природы этих компонентов, в частности углеводородных, ассортимент которых среди остаточных продуктов нефтепереработки весьма разнообразен. В ранее проводимых работах отсутствуют какие-либо количественные характеристики УВ компонентов битумов, используя которые можно было бы заранее прогнозировать свойства битумов из различных нефтей. Нет также сведений о пределах, в которых можно изменять свойства битумов путем регулирования качества УВ компонента битумов при различных соотношениях асфальтенов, смол и УВ.
According to positions of physicochemical mechanics, developed by academician P.A. Rebinder, correct understanding of the structure allows to solve competently issues of the creation a material with desired properties. There is a need to research the structure of bitumen over a wide range of ratios of its structural components and the influence to formation structures of nature of these components, particularly hydrocarbons, which assortment among the range of residual products of oil refinery is very diverse. In previously conducted works, there are no quantitative characteristics of hydrocarbon components of bitumen, using of which would be possible to predict in advance the properties of bitumen from different crude oils. There is also no information about limit, in which possible to change properties of bitumen by regulating quality hydrocarbon component of bitumen under different ratios asphaltenes, resins and hydrocarbons.
Ядерная магнитная релаксация (ЯМР) относится к числу современных инструментальных методов физикохимии и является фундаментальным свойством ядерного магнетизма, характеризующим динамику системы ядерных спинов в нефтяных дисперсных системах (НДС) /1-8/.
Открытие явления ЯМР в конденсированных средах позволило импульсной ЯМР-спектроскопии стать основой мощного метода исследования строения НДС и их свойств /1-14/.
В настоящее время трудно указать такую область естественных наук, где бы в той или иной степени не использовался импульсный ЯМР Изучение движения молекул и атомов в твердых телах, жидкостях и газах, исследование реакций ком-плексообразования, химического обмена, сольватации и межмолекулярных взаимодействий в растворах, определение параметров электронного строения металлов и сплавов, изучение электрон-ядерных взаимодействий, строение и свойств молекул вот неполный перечень задач, доступных решению с использованием ЯМР-спектроскопии.
Импульсный метод ЯМР все шире находит свое применение в экспресс анализах
нефтей и остаточных фракций /7-10/, особенно при оценке группового химического состава остаточного нефтяного сырья (ОНС) и битумов /8-11 /. Этот метод, как подчеркивает в своих работах Р.З. Сафиева, позволяет оценивать тонкие эффекты взаимного превращения масел, смол и асфальтенов, не определяемые традиционной хроматогра-фической методикой ВНИИНП адсорбци-онно-жидкостного разделения остаточных нефтепродуктов с применением экстраген-тов из-за большого временного интервала между отобранными пробами.
Таким образом, метод импульсного ЯМР привлекает сегодня не только физиков и химиков, но и специалистов в области технологии, что обусловлено, в первую очередь, практической значимостью данных о ЯМР, процессы которой зависят от многих физико-химических свойств исследуемой системы и превосходно иллюстрируют динамические явления на молекулярном уровне.
Следует отметить, что изучение спектров ЯМР основано на резонирующей системе ядер, которые чрезвычайно чувствительны к магнитному окружению, где локальные магнитные поля вблизи резонирующего ядра зависят от внутри- и
Тип системы Структурно-динамические характеристики системы Степень ТДО в НДС
(I) Тм - малое значение Т1А = Тм (или несколько выше)
Рм - высокое (порядка 80% масс.) Р1А>Р2А (до 40-50% масс.)
(II) Тд - высокие значения Т1А = Т2А (или несколько выше)
Рм - высокое (порядка 80% масс.) Р1А =>Р2А (до 50-70% масс.)
Табл. 1 Структурно-динамические состояния НДС в начальный момент времени (Тэксп = Ысх)
межмолекулярных эффектов. Это и определяет ценность ЯМР- спектроскопии для исследования строения и поведения многоэлектронных (молекулярных) систем.
Исследования НДС с помощью ЯМР-релаксометрии в первую очередь позволяют изучать температурное влияние на структурно-динамические состояния НДС на молекулярном уровне с последующим изменением термодинамических характеристик спин-спиновых составляющих НДС, таким образом, можем достоверно предположить структурные особенности НДС в целом.
Известно, что работа магнитного поля и переменных электромагнитных импульсов, действующих по методу Кара-Парсе-ла-Мейбума-Гила (КПМГ), способствует возникновению в спиновой системе явления намагничивания ядер, вследствие чего они начинают прецессировать с определенной частотой амплитуд колебаний, т.е. увеличивается молекулярная подвижность парамагнитных ядер. Это явление выражается временем спин-спиновой релаксации Т2|.
В определенный момент времени, при постоянной температуре эксперимента (Тэксп) насыщение намагниченностью ядер фазы I достигает своей постоянной дискретности и происходит «обмен» ядерной намагниченностью с окружающими спиновыми системами.
Это способствует увеличению молекулярной подвижности наиболее активных звеньев исследуемой фазы. Регистрируемый «термодинамический отклик» (ТДО) соседних сложных структурных единиц (ССЕ), выражается сегментарной ►
прецессиеи или амплитудами частот колебаний наиболее активных магнитных ядер и называется временем спин-решеточной релаксации Т1Г
Отметим, что в исследованиях импульсной ЯМР-спектроскопии спин-решеточной системой фазы I выступает совокупность спин-спиновых систем £1 /10/.
Однако следует отметить тот факт, что на ТДО структурной единицы решетки должна накладываться ее собственная внутренняя прецессия магнитных ядер - характеризуемая Т2|. Это явление заключается в том, что исследуемая НДС находится под действием внешних сил - постоянного магнитного поля, переменных электромагнитных импульсов и Т . Таким образом,
1 эксп ~ '
наблюдается прямо пропорциональная зависимость между временами спин-спиновой и спин-решеточной релаксации ССЕ и НДС в целом.
В связи с этим наличие постоянного магнитного поля при изучении структурно-динамических параметров в тяжелых нефтяных остатках (ТНО) объясняется обязательным наличием равенства значений времен спин-спиновой релаксации Т2| у ССЕ, находящихся в одной из фаз. В противном случае действие переменного магнитного поля изменит ТДО за счет уменьшения Т2|, и в итоге величина Т1| будет стремиться к нулю (Т^0).
Основываясь приведенными рассуждениями, экстремумы ТДО в НДС будут
наблюдаться у наиболее подвижных фаз А
и ^ Те. Т1д > Т1В > Т1С.
Возможные структурно-динамические состояния НДС в начальный момент времени = ^ . 25-3°°С)
Знания о структурно-динамических параметрах, состояниях термодинамической стабильности смолисто-асфальтеновых веществ (САВ) в НДС позволяют ее охарактеризовать «экспресс анализом», т.е. используя данные ЯМР-релаксометрии исходного состояния нефтяной системы (времена релаксаций, населенностей протонов фаз, межпротонные расстояния и коэффициенты диффузии, а также расчетные термодинамические характеристики НДС - энергия активации, энтропия, энтальпия и т.д.) позволяют выдвигать теоретические предпосылки ее физико-химической, термодинамической активности при деструктивной термообработке.
Таким образом, основываясь значениями амплитуд частот термоколебаний магнитных ядер различных молекул, насе-ленностей окружающих их протонов, можно установить тип сегментарного строения углеродной цепи молекулы.
Таблица 1 Структурно-динамические состояния НДС в начальный момент времени (Т = t )
4 эксп исх'
Следует отметить, что трудности в определении степени замещенности, а также длины алифатической цепи фазы А
Рис. 1 - Температурные зависимости времен релаксации гудронов
Показатели Гудрон-1 Гудрон-2
Плотность Вязкость условная ВУ80 0,9686 22,96 0,9878 51,76
Содержание, % масс.
САВ 18,25 28,23
Сера 0,492 0,887
Парафины <2 до 10
Асфальтены/смолы 0,64 0,45
Табл. 2 - Физико-химические свойства гудронов
в начальный момент времени объясняются ее «замороженным» состоянием. При высоких значениях времени релаксации и населенности протонов наиболее подвижной фазы А в спин-спиновой системе Т2д и Р2А ( II ), согласуя это с населенностью протонов водорода и временем релаксации спин-решеточной системы - Р1А и Т1А, где ТДО численно равен количеству протонов в спин-спиновой системе Р1А ~ Р2А (либо незначительно ниже), строение заместителей фазы А можно охарактеризовать, с определенной долей вероятности, как короткоцепочное с высокой степенью разветвленности.
Наличие в алкильных углеродных цепях двойных связей можно определить из значений Т2А и Р2А, учитывая, что длина С=С связи (0,134 нм) меньше длины С-С (0,154 нм), а энергия связи 6С=С (605 кДж/ моль) больше 6С-С (350 кДж/моль). Это обстоятельство подтверждается тем, что при увеличении количества С=С двойных связей, значения Т2А будут резко уменьшаться /1/ и, таким образом численные характеристики времени релаксации фазы А спин-спиновой системы Т2А будут близки к значениям времени релаксации менее подвижной фазы В спин-спиновой системы Т2В, в этом случае численность протонов фазы В - Р2В будет иметь сравнительно малые значения. Вместе с тем значения времени релаксации фазы А спин-решеточной системы - Т1А могут быть кратно ниже значений времени спин-спиновой системы Т2А, так как наличие двойных связей в ОНС способствует «упорядочиванию», т.е. проявляет склонность к дальнейшему межмолекулярному структурированию.
Согласно системе (I) в табл.1 можно сделать предположение, что за отклик Р1А ответственны наиболее активные замещенные алкильные цепи (преимущественно короткоцепочные, которыми могут быть метил-, этил заместители), но, учитывая малые значения Т2А Т1А, следует, что причиной отклика могут являться стерические затруднения углеродной цепи. Таким образом, основная сегментарность звеньев фазы А выражается длиноцепочным парафиновым характером.
Рассматривая ОНС как термодинамическую систему фаз А, В и С, необходимо отметить, что исходя из молекулярной подвижности и структурно-динамических параметров фаз, можно с большой долей уверенности предполагать строение их углеродных звеньев. Ранее было установлено, что в структурно-динамическом анализе (СДА) нефтяных дисперсных систем присутствуют фазы А, В и С, которые ввиду различного их содержания и молекулярной подвижности условно отнесены к маслам, смолам и асфальтенам соответственно.
В качестве исходного сырья процесса окисления исследовались на релаксометре ЯМР 08БК/РС, изготовляемого по ТУ 25-4823764.0031 -90, структурно-динамические состояния ОНС нафтено-ароматического - Карабашского НБЗ (гудрон-1) и парафино-нафтенового - Елховского НПУ (гудрон-2) оснований при различных температурных режимах, физико-химические свойства которых представлены в табл.2. ►
Для анализа гудронов использовались следующие режимы измерений: период запуска Т = 500 мс - 2 с, интервал между 900 и 1800 импульсами N = 5 - 1000, число накоплений п = 3 - 50. Время анализа не превышало 3 мин.
Значения энергий активации (ЕА) молекулярного движения в предположении соотношения Аррениуса для времен корреляции с = 0ехр(ЕАЖТ) получены с использованием формулы (1) /1/:
где Т(1), Т(2) - температуры, при которых определяются значения времен релаксации Т2(2) и Т2(1).При этом учитывается лишь внутримолекулярный вклад в релаксацию, обусловленный диполь-дипольными взаимодействиями от вращательной диффузии, который в обычно выполняющихся для НДС условиях высокотемпературного приближения ш 0<<1, для протонов со спином I = % записывается в виде:
Огибающая спада сигналов спин-эхо, как гудронов, так и полученных на их основе окисленных битумов может быть разделена на несколько компонент | и описана формулой:
Согласно рис.1 в гудронах наблюдается прямолинейная зависимость между температурой эксперимента (Тэксп) и временами релаксации Т2|, т.е. 103/Т прямо пропорциональна Т2|, однако Т2С для гу-дрона-1 значительно отличается от Т2С для гудрона-2.
Первоначально на рис.1 можно обнаружить незначительные различия в молекулярной подвижности компонентов исследуемых гудронов, характеризуемые временами спин-спиновой релаксации Т2|.
Однако более детальный анализ представленных экспериментальных данных показывает, что Тэксп, при которой проявляется фаза С, для каждого гудрона различна. Так для гудрона-1 это 720С, а для гудрона-2 - 49,60С.
Высокая Т , при которой начинает
эксп
проявляться наименее подвижная фаза С, объясняется стерическими затруднениями молекулярной структуры при ее движении, то есть можно предположить наличие более упорядоченной структуры молекул асфальтенов гудрона-2, чем в гудроне-1. Вместе с этим из рис.1а видно, что на всем интервале Тэксп у гудрона-1 не наблюдается резких изменений времени релаксации Т2С, что показывает термодинамическую стабильность ее структуры, на наш взгляд, состоящую преимущественно из звеньев (сегментов) полициклической нафтено-ароматической структуры с невысокой степенью замещенности алкиль-ными цепями. При сравнительном анализе рис.1а видно, что при увеличении Тэксп выше 60 0С отчетливо проявляется молекулярная подвижность замещенных звеньев фазы С; однако, если фаза С в гудроне-1, который при 720С переходит в ньютоновскую жидкость, характеризуется временем Т2, которая напоминает «дрожание», то есть прецессии парамагнитных ядер с незначительной амплитудой частот колебаний, то в гудроне-2 при Тэксп выше 600С фаза С характеризуется термодинамической нестабильностью, то есть ССЕ при изменении внешних условий (увеличение Тэксп) не стремится к термодинамическому равновесию - к минимуму энергии взаимодействия. Это может быть вызвано как строением самой дисперсной фазы, так и степенью замещенности, разветвленнос-тью и длиной цепи заместителей. По нашему мнению, характер замещения фазы С выражен длинными алкильными цепями с малой степенью разветвления, что полагаем, объясняется значениями населен-ностей протонов Р2С (рис.3), где показана
невысокая концентрация в гудроне-1 фазы С (в пределах 20% мас.). Таким образом, из сравнительного анализа значений времени релаксации и населенности протонов водорода фазы С - Р2С и Т2С в гудроне-2 и гудроне-1 можно судить о наличии малой степени ароматичности и большого содержания парафино-нафтеновых структур.
На основе зависимостей Т2А и Т2В от 103/Т отмечаем, что время Т2А для гудро-на-2 составляет 2,5 мс, что заметно превышает время Т2А гудрона-10,65 мс. Вместе с этим Тэксп, при которой проявляется наиболее подвижная фаза А, для гудро-на-1 составило 28,20С и оказалось соответственно меньше, чем для гудрона-2 и составило 31,80С. Также следует отметить, что при одинаковых Тэксп 103/Т для обоих гудронов - 3,28К наблюдаются различные времена релаксации Т2А. Так для гудрона-1 это 1 мс, а для гудрона-2 составляет 2,5 мс. Следует отметить, что при 720С>Т>-55 0С Т2А для гудрона-1 больше, чем для гудрона-2 в среднем на 8% мас. Однако при Т>720С наблюдается обратная зависимость, объясняемая наличием фазовых переходов (ФП), связанных со структурно-динамической перестройкой молекулярной морфологии. Таким образом, можно заключить, что наиболее подвижная фаза А гудрона-2 характеризуется длинноцепоч-ными парафиновыми звеньями, а гудрона-1 - преимущественно короткоцепочными алифатическими звеньями, являющиеся структурными сегментами более высокомолекулярной части дисперсионной среды и соответственно наименее подвижной - смол, характеризуемых временем релаксации Т2В.
Основываясь на экспериментальных данных молекулярной подвижности фазы В можно отметить, что на всем интервале 103/Т время релаксации Т2В для гудрона-1 меньше, чем для гудрона-2. Этот факт служит характеристикой химического строения фазы В исследуемых ОНС, которая, как мы считаем, в гудроне-1 состоит ►
Рис. 2 - Зависимость энергии активации от температуры
Рис. 3 - Температурная зависимость населенностей протонов фаз в гудроне-2
преимущественно из полициклических нафтено-ароматических структур, при этом фаза В в гудроне-2 характеризуется в основном полициклическими парафи-но-нафтеновыми структурами, что также подтверждается значениями времен спин-решеточной релаксации Т^ (рис.3) и насе-ленностей протонов Р,. (рис.3).
Исследования наиболее подвижных фаз А и В - Р2А и Р2В показывают симбат-ность графических зависимостей (рис.3). В начальный момент времени измерения в интервале Т - 103/Т 3,13-3,165К
г эксп ' '
Р2А у гудрона-2 (=47,5%мас.) значительно превышает значения населенностей протонов гудрона-1 (=62%мас.), при этом наблюдаемое скачкообразное изменение Р2А гудрона-2, стабилизирующееся при температуре 550С до 50%мас., изменяется в дальнейшем в незначительной степени в пределах 40%мас.
При сравнительном анализе времен спин-спиновой релаксации Т2А, В и населенности протонов Р2А,В для гудрона-1 наблюдается явление, когда при температурах 103/Т 3,28К и 3,09К соответственно в момент времени Т2 для фаз А и В - населенности протонов обнаруживаются при Т 103/Т 3,165К и 2,96К.
г эксп ' '
Этот факт можно объяснить исходя из того, что при заданной Тэксп внутримолекулярному движению подвержены только некоторые, наиболее подвижные, сегменты молекул фаз.
Таким образом, при рассмотрении изменений населенностей протонов Р2С (рис.4) можно предположить, что фаза С гудрона-1 представлена полициклической нафтено-ароматической структурой, а в связи с тем, что на рис.4а Р2С характеризуется термодинамической нестабильностью, в отличие от Р2С гу-дрона-2 (рис.3), выражающейся экстремумами зависимости Р2С от 103/Т, то можно предположить наличие высокого химического сродства фазы С с компонентами смол.
Соизмеряя величину термодинамического «отклика», то есть отношения времен Т^ к населенностям протонов решетки Р1М со значениями времен спин-решеточной релаксации Т„, населенностей протонов Р^ у исследуемых ОНС, а, также учитывая энергии активации (ЕА), можно охарактеризовать сегментарную термодинамическую активность элементов решетки, вместе с этим с помощью величин Р^ можно дать оценку степени упорядоченности спин-решеточной системы. При этом, предварительно изучив с помощью Т2. и Р2. спин-спиновую систему магнитных ядер ТНО, исследование параметров Т^ и Р^ может выявить количественное распределение фаз в спин-решеточной системе.
В связи с этим, графические зависимости для гудрона-1 времен спин-решеточной релаксации Т1. находятся в плоскопараллельном направлении согласно зависимостям спин-спиновой релаксации Т2. (рис.1), также как и населенности протонов Р1. лежат параллельно друг другу. Наличие экстремумов (рис.4б) для Р1А и Р1С объясняется наличием ФП. Увеличение времени Т1С в интервале температур 2,6-2,75 К может являться также причиной наличия ФП, связанного со структурно-динамической перестройкой молекулярной морфологии, а также вовлечения в молекулярную подвижность фазы С некоторых структурных звеньев молекул смол, проявляющих наибольшее к ней химическое сродство.
Анализ рис.2 - зависимости значений энергии активации от температуры нагрева гудронов, подтверждает предпосылки о возможном химическом строении компонентов в ССЕ гудронов. Так из рис.2 видно, что при низких Тэксп энергии активации в гудроне-1 выше, а при более высоких значениях ниже, чем в гудроне-2. При этом максимальные значения ЕА близки к значениям, характерным для потенциальных
барьеров при стерических затруднениях молекулярного движения длинноце-почных полимеров типа полиэтилена с ЕАпол~12-14 ккал/моль /11/. Поведение Еа для гудрона-2 в зависимости от Тэксп характеризуется при низких Тэксп «замороженным» состоянием сегментов длинных алифатических цепей и высокой подвижностью концевых и боковых цепей, а при высоких Т - стерическими
эксп
затруднениями молекулярной сегментальной подвижности.
ИСПОЛЬЗОВАННАЯ ЛИТЕРАТУРА:
1. Вашман А.А., Пронин И.С. Ядерная магнитная релаксация и ее применение в химической физике. М.: Наука, 1979, 236 с.
2. Воронцов В.К. Ядерный магнитный резонанс. Соросовский образовательный журнал №10, 1996.
3. Керингтон А., Мак-Лечлан Э. Магнитный резонанс и его применение в химии. М.:, 1970. 447с.
4. Туманян Б.П. Научные и прикладные аспекты теории нефтяных дисперсных систем. М.: Техника, 2000. - 335с.
5. Сафиева Р.З. Физикохимия нефти. М.: Химия, 1998. - 448с.
6. Кемалов, А.Ф. Научно-практические основы физико- химической механики и статистического анализа дисперсных систем (учебное пособие)/ А.Ф. Кемалов, Р.А. Кемалов. - Казань: КГТУ. 2008.-472 с/264с.
7. Кемалов Р.А. Модифицированные специальные битумы и лакокрасочные материалы на их основе: Ав-тореф.дисс.канд.тех.наук.-Казань, 2003, 20с.
8. Кемалов А.Ф. Интенсификация производства окисленных битумов и модифицированные битумные материалы на их основе: Автореф.дисс. докт.тех.наук.-Казань, 2005, 42 с.
9. Кемалов А.Ф., Прокопьев В.П. // В сб. тр. междунар. конф. «Проблемы комплексного освоения трудноизвлека-емых запасов нефти и природных битумов (добыча и переработка)». Казань, ИОФХ КНЦ РАН, 4-8 октября 1994. Т5. с.1608-1610.
10. Кемалов А.Ф. и др. Связь между ЯМР параметрами и эксплуатационными характеристиками битумов // Химия и технология топлив и масел. - 1999. - №2. - с.37.
11. А.с. СССР №1583823 «Способ определения группового состава нефтяных остатков» / Кемалов А.Ф., Фах-рутдинов Р.З., Прокопьев В.П. и др. // Бюлл. изобр. - 1990. - № 29
12. Кемалов А.Ф и др. Динамика структурного упорядочения в нефтяных остатках и наполненном углеродом полиэтилена. Международная научно-практическая конференция «Нефть и газ - 99». Казань 1999, Т2, с.362-367.
13. Прокопьев В.П. и др. Исследование импульсным методом ЯМР процесса структурирования полиэтилена. Высокомолекулярные соединения.-1984, Т.ХХУ1Б., с.115. ■
Рис.4 - Зависимость населенностей протонов фаз Pi гудрона-1 от температуры
