Каталитический риформинг бензинов. Теория и практика Крылов Валерий
Рис. 14. Зависимость коэффициента
эффективности от модуля Тиле
Основные реакции платформинга являются эндотермическими.
Эндотермичность приводит к снижению температуры в грануле катализатора и уменьшению скорости реакции по сравнению с таковой на поверхности гранулы. Это приводит к дополнительному снижению коэффициента эффективности, что демонстрируется на рис. 15.
Рис. 15. Изменение коэффициента эффективности в зависимость от :
= (–Н) DeCs/Ts; T = Ts, где Н – энтальпия химической
реакции, De – коэффициент диффузии, Сs – концентрация реагента на поверхности гранулы, – коэффициент теплопроводности материала гранулы, Тs – температура на поверхности гранулы,
T – градиент температуры по сечению гранулы, < 0 для эндотермических реакций и > 0 для экзотермических реакций
С учетом небольшого размера молекул, участвующих в реакциях риформинга, можно было ожидать, что диффузионные затруднения будут отсутствовать, однако использование катализаторов с высокой удельной поверхностью и, следовательно, с узкими порами, с одной стороны, и высокие скорости реакций, связанные с применением платины и повышенных температур, с другой стороны, объективно создают условия для протекания реакций риформинга во внутридиффузионном режиме.
По данным [17], дегидрирование циклогексана уже при температуре 430 С на платиновом катализаторе с диаметром гранулы 3 мм (размер гранул монометаллических катализаторов риформинга первого поколения) протекает в условиях сильного диффузионного контроля: коэффициент эффективности составил 0,38. Это означает, что более 60 % удельной поверхности катализатора используется неэффективно из-за сильного падения концентрации циклогексана в центральной части гранулы.
Снижение размера гранулы с 3,0 до 1, 6 мм, применяемых в современных катализаторах, позволяет уменьшить модуль Тиле в 1,9 раза и за счет этого увеличить эффективность использования внутренней поверхности катализатора.
В связи с тем что зависимость константы скорости химической реакции от температуры является экспоненциальной, а температурная зависимость коэффициентов диффузии представляет степенные функции (Т3/2 для молекулярной диффузии и Т1/2 для кнудсеновской диффузии), при рабочих температурах платформинга коэффициент эффективности будет еще ниже.
Оценка c учетом изменения размера гранулы и температуры дает значение для реакции дегидрирования циклогексана в условиях риформинга, равное 0,31.
Как было показано ранее, реакция замыкания кольца является самой медленной стадией дегидроциклизации и контролирует скорость ароматизации парафиновых углеводородов. Скорость реакциидегидроциклизации примерно в 30 раз меньше, чем скорость дегидрирования циклогексана.
Поскольку модуль Тиле изменяется пропорционально квадратному корню из величины скорости реакции, его значение для реакции циклизации будет в первом приближении в 5,5 раза меньше – около 0,06. Это означает, что циклизация может проходить в кинетическом режиме, и ее ускорение возможно за счет уменьшения размера пор и увеличения удельной поверхности носителя.
Снижению диффузионных затруднений способствует также экзотермичность реакции циклизации.
Вместе с тем необходимо учитывать, что быстрые реакции расширения кольца и последующего дегидрирования, являющиеся промежуточными стадиями в ароматизации, будут испытывать все большее диффузионное торможение при уменьшении размера пор, и при определенном размере пор будут снижать общее превращение.
Оптимальная поровая структура должна обеспечивать баланс между реакциями, контролируемыми кинетикой, и реакциями, находящимися под диффузионным контролем.
В работе [46] приведены результаты исследования реакции дегидроциклизации н-гептана на платинорениевых катализаторах в зависимости от величины удельной поверхности и распределения пор носителя. В работе использовался катализатор с размером частиц 0,3–0,75 мм. Целью исследования
было определение поровой структуры катализатора риформинга, обеспечивающей максимальную конверсию н-гептана
в толуол.
Некоторые выводы по материалам исследования представлены ниже. В частности, было показано, что активность катализаторов, имеющих одинаковое содержание металлов, но различную поровую структуру, может отличаться в несколько раз.
1. Отсутствует прямая корреляция между активностью
и величиной удельной поверхности катализатора; при оценке активности должна учитываться поровая структура катализатора.
2. При увеличении температуры процесса различие в активности между катализаторами возрастает, этот эффект противоположен тому, что наблюдается при протекании реакции в кинетическом режиме.
3. Кроме того меняются лидеры: так, катализатор, показавший максимальную активность при 470 С, не является таковым при температуре 510 С, и, наоборот, катализатор-аутсайдер становится лидером при более высокой температуре.
4. Увеличение температуры приводит к росту диффузионных затруднений, и реакция перемещается в поры большего размера. Так, при температуре 470 С реакция протекает в порах с радиусом 3 нм и выше, а при 510 С доступными для реакции становятся поры с радиусом 5 нм и более.
5. При увеличении радиуса пор диффузионные затруднения уменьшаются, процесс сдвигается в сторону кинетического режима, и эффективная константа скорости химического превращения приближается к величине действительной константы скорости, достигаемой в отсутствие диффузионных ограничений; это приводит к более эффективному использованию поверхности поры: удельная активность (моль/с), отнесенная к поверхности поры, увеличивается с ростом радиуса поры.
6. Необходимо иметь в виду, что в исследовании использовался катализатор с размером гранул 0,3–0,75 мм. Очевидно, что на промышленном катализаторе, имеющем средний радиус пор 5 нм, но больший размер частиц (1,6 мм), реакция дегидроциклизации н-гептана протекает во внутридиффузионном режиме.
Глава 8. Механизмы реакций риформинга
Схема Миллса для реакций на бифункциональном катализаторе риформинга.
С5– и С6-дегидроциклизация алканов.
Скелетная изомеризация. Карбениевые и карбониевые ионы
Обнаружение роли кислотного промотора привело к разработке концепции бифункционального катализа, предложенной впервые G. Mills, H. Heineman, T. Milliken, A. Oblad из лаборатории Houdry Process Corporation. Концепция оказалась плодотворной и полезной также для понимания механизмов превращений и разработки эффективных катализаторов для процессов изомеризации н-бутана и н-пентана, гидроизомеризации дизельных фракций, гидрокрекинга различных нефтяных фракций и др.
Основная идея этой концепции – ускорение химических реакций за счет синергизма, обусловленного одновременным присутствием кислотных и металлических центров.
Схема реакций платформинга с участием металлических (М) и кислотных (А) активных центров представлена ниже (рис. 16).
Современная схема дополнена реакциями изомеризации и С6-циклизации на металлических центрах.
В соответствии с предложенной схемой ряд реакций может протекать без участия кислотных центров, и для них необходимы только металлические центры, образованные платиной. К таковым относятся реакции дегидрирования и гидрирования, гидрогенолиза и С6-циклизации.
Изомеризация олефинов, циклоолефинов, гидрокрекинг олефинов и циклоолефинов и С5-циклизация, наоборот, не требуют наличия металлических центров и осуществляются на кислотных центрах катализатора.
Рис. 16. Схема реакций платформинга по Миллсу
Более сложные реакции проходят с участием активных центров того и другого типа. Например, С5-дегидроциклизация
н-парафинов включает следующие стадии: дегидрирование и образование н-олефина (М), изомеризацию с образованием изоолефина (А), циклизацию изоолефина с образованием алкилциклопентана (А), дегидрирование алкилциклопентана до циклического олефина (М), расширение цикла с образованием 6-членного циклопарафина (А), дегидрирование циклопарафина с получением ароматического углеводорода (М).
Медленными стадиями превращения являются реакции циклизации олефина с образованием 5-членного кольца и последующая реакция расширения цикла. Обе реакции проходят с участием кислотных центров.
Тот факт, что при риформинге н-гексана введение н-бутиламина приводит к снижению концентрации циклогексана при одновременном накоплении метилциклопентана, может свидетельствовать о том, что для осуществления превращений требуются разные кислотные центры [47], в частности, для реакции циклизации кислотные центры Льюиса, представленные ионами Al + 3.
Снижение содержания хлора и, следовательно, бренстедовской кислотности приводит к накоплению молекул циклического олефина (метилциклопентена в случае ароматизации гексана) на поверхности платины и более глубокому дегидрированию с образованием производных циклопентадиена, легко полимеризующегося с образованием предшественников кокса.
Альтернативный маршрут С6-дегидроциклизации, протекающий на катализаторе Pt/C и Pt/KL-цеолит (катализатор RZ-100), в условиях платформинга на кислотном носителе является второстепенным механизмом ароматизации алканов.
Изомеризация н-парафиновых углеводородов включает в себя дегидрирование на металлических центрах с образованием олефинов нормального строения, изомеризацию полученных олефинов на кислотных центрах с образованием олефинов изостроения и гидрирование на металлических центрах до изопарафина. На примере этого превращения становятся понятными преимущества бифункционального катализатора.
Промежуточным соединением в реакции скелетной изомеризации является карбениевый ион, имеющий структуру с трехкоординированным атомом углерода.
В случае олефина образование карбениевого иона происходит достаточно легко при атаке протона, предоставляемого кислотным центром Бренстеда, на -связь. При этом -связь разрывается и образуется новая С–Н-связь.
Если же исходным углеводородом является алкан, то
образование карбениевого иона включает атаку протона на
С–Н-связь. При этом вначале образуется промежуточный карбониевый ион, представляющий катион с пятикоординированным атомом углерода. Например, при протонировании 2-метилпентана образуется карбониевый ион (рис. 17, б). Распад этого катиона приводит к образованию карбениевого иона и молекулы Н2 (см. рис. 17, а).
Рис. 17. Схема образования:
а – карбениевого иона; б – карбониевого иона
В этом случае разрывается более прочная С–Н-связь.
Значительная разница в энергиях связи (394–239 = 155 кДж/моль) делает возможным образование карбениевого иона из алкана лишь при использовании твердых суперкислот, например, Н-морденитов с силикатным модулем около 20, в условиях же платформинга подобный механизм не реализуется.
Глава 9. Катализ на
d
-металлах
Природа каталитических свойств d-металлов. Связь с координационной ненасыщенностью поверхностных атомов металлов.
Донорно-акцепторная и дативная связи в металлорганических комплексах. Модель Дьюра – Чата – Дункансона.
Механизм образования связей на примере молекулы СО.
Ослабление связей в молекуле СО как результат образования донорно-акцепторной и дативной связи.
Как образуются d-зоны в металлических катализаторах.
Почему происходит химическая адсорбция.
Активация молекул как результат хемосорбции.
Принцип Сабатье и вулканообразные кривые Баландина.
Почему платина является базовым элементом катализаторов риформинга.
Почему реакции дегидрирования являются быстрыми реакциями, а реакции гидрогенолиза медленными: каталитический эффект платины
Каталитическая активность d-металлов обусловлена координационной ненасыщенностью атомов, образующих поверхностные грани металлических частиц.
Координационное число (КЧ) платины, формирующей гранецентрированную кубическую решетку (ГЦК), равно 12.
КЧ для поверхностных атомов зависят от типа поверхности и составляют 9 для граней (111), 8 для (100) и 7 для (110).
Атомы на ступенях и изломах, ребрах и углах частиц имеют еще меньшие КЧ – от 7 до 5 [54].
Вывод атома из объема металла на поверхность является сильно эндотермическим процессом, связанным с разрывом связей с соседними атомами.
Энергия, необходимая для образования поверхностного атома, прямо пропорциональна энергии когезии металла и координационной ненасыщенности поверхностного атома.
Уменьшение размера частицы также приводит к увеличению поверхностной энергии за счет увеличения доли поверхностных атомов.
Переход системы в более устойчивое состояние с меньшей энергией Гиббса достигается путем коалесценции частиц при повышенных температурах или за счет адсорбции молекул окружающей среды. Из двух видов адсорбции, физической и химической, последняя имеет ключевое значение для гетерогенного катализа, так как связана с активацией молекулы, обусловленной изменениями ее электронной структуры при адсорбции на поверхности твердого тела.
В основе современного понимания механизма химической адсорбции и катализа на d-металлах лежат идеи о координационно-донорной и дативной связях и d-зоне, которые были заимствованы из металлорганической химии и физики твердого тела соответственно.
Дьюар в 1951 году предложил модель образования соли Цейзе и ее палладиевого аналога, комплекса Караша, представляющих собой комплексы этилена и металла (рис. 18) [37].
Рис. 18. Структура комплексов Pt(Pd) c этиленом
В соответствии с этой моделью, в доработанном виде носящей название модели Дьюара – Чата – Дункансона, в образовании комплекса принимают участие два типа связей: донорно-акцепторная связь, образуемая за счет передачи электронной плотности -связи молекулы этилена на вакантную d-орбиталь атома платины, и дативная связь, которая возникает за счет перекрытия заполненной d-орбитали атома металла с разрыхляющей орбиталью молекулы этилена. Заметим, что обе связи являются примером донорно-акцепторного взаимодействия, так что выделение дативной связи сделано для удобства, это указание на то, что донором в этом случае является металл.
Атомными орбиталями, удовлетворяющими этому требованию, являются dz2– и dxz-орбитали переходного металла.
Ниже представлены схемы образования донорно-акцепторной и дативной связей d-металла и молекулы этилена (рис. 19). Донорно-акцепторная связь образуется при перекрывании -МО этилена с dz2-AO металла. В образовании дативной связи участвуют разрыхляющая *-МО этилена и dxy-АО металла.
Рис. 19. Схема образования – и -связей:
стрелками показаны направления смещения электронной плотности
Образование донорно-акцепторной связи осуществляется по -типу, а дативной связи – по -типу.
Для образования дативной связи возможны два варианта перекрывания орбиталей.
Из-за небольших стерических затруднений, возникающих при боковом перекрывании, в рассмотренных комплексах реализуется схема с расположением иона металла над или под плоскостью, в которой находятся sp2-орбитали молекулы этилена.
При отсутствии таких ограничений может реализовываться схема с боковым перекрыванием, например, при образовании связи с молекулой СО, где такому перекрыванию способствует также несимметричное распределение электронной плотности в лепестках разрыхляющей орбитали, связанное с поляризацией связи.
Прочность донорно-акцепторной и дативной связи увеличивается с уменьшением различия в энергии донорной и акцепторной орбиталей в соответствии с величиной энергии стабилизации 106; 107:
Естаб ~ S2/,
где S – интеграл перекрывания; – разница в энергии исходных орбиталей.
При переносе электронной плотности с -орбитали этилена происходит накопление положительного заряда в молекуле, что ограничивает перенос электронов, в то же время обратный перенос с занятых орбиталей металла нейтрализует этот заряд, и позволяет продолжить формирование более прочной донорно-акцепторной связи.
В свою очередь передача электронной плотности с молекулы на металл увеличивает донорные свойства металла.
В итоге имеет место синергизм, который приводит к образованию более прочной связи металла и молекулы и более значительному ослаблению связи в молекуле.
Для комплексов Цейзе и Караша связывание является слабым из-за пониженного дативного потенциала положительно заряженных ионов платины и палладия.
Результатом рассмотренных взаимодействий является уменьшение порядка и прочности углерод-углеродной связи в молекуле, что коррелирует с увеличением длины связи и со смещением пиков инфракрасного поглощения в длинноволновую область спектра.
Так, длина связи С–С увеличивается со 133,7 пм в свободной молекуле этилена до 137,0 пм в комплексе платины и этилена и до 148,0 пм в комплексе этилена и никеля Ni(CO)4 [37].
Квантово-механический анализ взаимодействия на примере связывания молекулы СО и переходных металлов первой серии (3d-металлы) представлен в [106].
Электронная структура молекулы СО представлена на рис. 20.
Рис. 20. Электронная структура молекулы СО:
HOMO – high occupated molecular orbital;
LUMO – low unoccupated molecular orbital
В образовании связей в молекуле принимают участие
1-орбиталь и две 1-орбитали.
Молекулярные орбитали 2 и 3 не вносят вклада в связывание и являются фактически не поделенными электронными парами кислорода и углерода соответственно.
На самом деле орбиталь 3 является слегка разрыхляющей орбиталью для молекулы. Эта орбиталь участвует в донорно-акцепторном связывании с d-металлом в качестве донора электронов.
Особенностью электронного строения молекулы СО является несимметричное распределение электронной плотности между лепестками - и *-МО. Несимметричность обусловлена различиями в электроотрицательности и в уровнях энергии атомов углерода и кислорода, образующих молекулу.
-МО является связующей в молекуле, и по энергии она ближе к энергии p-AO кислорода, в связи с чем электронная плотность смещена в лепесток у атома кислорода.
*-орбиталь как разрыхляющая орбиталь молекулы ближе по энергии к p-орбитали атома углерода, что обусловливает концентрацию электронной плотности орбитали в лепестке у атома углерода (рис. 21).
Рис. 21. Электронные схемы *-, -орбитали
В связи с тем что -связывание с металлом происходит через атом углерода, подобное несимметричное распределение -электронной плотности более благоприятно для связывания dxz-орбитали металла с *-орбиталью (рис. 22).
Расчетное отношение интегралов перекрывания для - и *-орбиталей составляет 1,78 в пользу *-орбитали.
Поскольку энергия стабилизации при образовании связи пропорциональна квадрату интеграла перекрывания, то выигрыш в энергии при образовании связи с *-орбиталью больше в 3,2 раза. Такой выигрыш благоприятен для проявления -акцепторного характера *-МО. Вместе с тем необходимо учитывать различие в энергиях орбиталей .
Значения представлены в табл. 5.
Рис. 22. Схема dxz-орбиталей атома металла
с *-МО и -МО молекулы СО
Таблица 5
Энергетические уровни d-металлов
и орбиталей молекулы СО [106]
Энергия
Sc
Ti
V
Cr
Mn
Fе
Co
Ni
Сu
d
–8,5
–10,8
–11,0
–11,2
–11,7
–12,6
–13,2
–13,5
–14,0
7,2
4,9
4,7
4,5
4,0
3,1
2,5
2,2
1,7
*
–0,6
1,7
1,9
2,1
2,6
3,5
4,1
4,4
4,9
П р и м е ч а н и е : d – энергия d-электронов в атоме металла.
Напомним, что нулевому энергетическому уровню соответствует состояние электрона, находящегося на таком удалении от ядра, когда можно пренебречь электростатическим взаимодействием электрона и ядра атома.
При приближении электрона к ядру его потенциальная энергия падает, поэтому чем больше по величине отрицательное значение d, тем ниже энергия электрона.
Для металлов в левой части периода d–*-взаимодействие оказывается сильнее d–-взаимодействия, и для этих металлов молекула СО выступает как акцептор электронов.
Для металлов в правой части периода разница в энергии благоприятствует проявлению донорного характера -связи, однако больший интеграл перекрывания с *-связью приводит все же к тому, что молекула остается -акцептором электронов.
При переходе к металлам второго и третьего периодов происходит подъем энергетического уровня d-электронов, что приводит к сближению *-МО и d-электронов и увеличению
-акцепторного характера молекулы СО (табл. 6).
Таблица 6
Энергетические уровни d-орбиталей первого, второго
и третьего переходного периодов [106]
Первый период
Sc
Ti
V
Cr
Mn
Fе
Co
Ni
Сu
Zn
3d
–7,92
–9,22
–10,11
–10,74
–11,14
–11,65
–12,12
–12,92
–13,46
–17,29
4d
–6,60
–7,11
–7,32
–7,45
–7,83
–7,90
–8,09